Bayesian Approach for Computing Free Energy on Perturbation Graphs with Cycles.

A common approach for computing free energy differences among multiple states is to build a perturbation graph connecting the states and compute free energy differences on all edges of the graph. Such perturbation graphs are often designed to have cycles. Because free energy is a function of states, the free energy around any cycle is zero, which we refer to as the cycle consistency condition.

CHARMM at 45: Enhancements in Accessibility, Functionality, and Speed.

Since its inception nearly a half century ago, CHARMM has been playing a central role in computational biochemistry and biophysics. Commensurate with the developments in experimental research and advances in computer hardware, the range of methods and applicability of CHARMM have also grown. This review summarizes major developments that occurred after 2009 when the last review of CHARMM was published.

Optimizing Force Fields with Experimental Data Using Ensemble Reweighting and Potential Contrasting.

Despite force field improvements over the past decades, we still encounter situations where simulation results disagree with experiments due to force field inaccuracies. Such situations provide opportunities to improve force fields. In this study, we introduce a novel framework for optimizing force fields using experimental data.

Bayesian Multistate Bennett Acceptance Ratio Methods.

The multistate Bennett acceptance ratio (MBAR) method is a prevalent approach for computing the free energies of thermodynamic states. In this work, we introduce BayesMBAR, a Bayesian generalization of the MBAR method. By integration of configurations sampled from thermodynamic states with a prior distribution, BayesMBAR computes a posterior distribution of free energies.

Transferable Implicit Solvation via Contrastive Learning of Graph Neural Networks.

Implicit solvent models are essential for molecular dynamics simulations of biomolecules, striking a balance between computational efficiency and biological realism. Efforts are underway to develop accurate and transferable implicit solvent models and coarse-grained (CG) force fields in general, guided by a bottom-up approach that matches the CG energy function with the potential of mean force (PMF) defined by the finer system. However, practical challenges arise due to the lack of analytical expressions for the PMF and algorithmic limitations in parameterizing CG force fields.

Fast free energy estimates from λ-dynamics with bias-updated Gibbs sampling.

Relative binding free energy calculations have become an integral computational tool for lead optimization in structure-based drug design. Classical alchemical methods, including free energy perturbation or thermodynamic integration, compute relative free energy differences by transforming one molecule into another. However, these methods have high operational costs due to the need to perform many pairwise perturbations independently.

Transferable Coarse Graining via Contrastive Learning of Graph Neural Networks.

Coarse-grained (CG) force fields are essential for molecular dynamics simulations of biomolecules, striking a balance between computational efficiency and biological realism. These simulations employ simplified models grouping atoms into interaction sites, enabling the study of complex biomolecular systems over biologically relevant timescales. Efforts are underway to develop accurate and transferable CG force fields, guided by a bottom-up approach that matches the CG energy function with the potential of mean force (PMF) defined by the finer system.

OpenABC enables flexible, simplified, and efficient GPU accelerated simulations of biomolecular condensates.

Biomolecular condensates are important structures in various cellular processes but are challenging to study using traditional experimental techniques. In silico simulations with residue-level coarse-grained models strike a balance between computational efficiency and chemical accuracy. They could offer valuable insights by connecting the emergent properties of these complex systems with molecular sequences.

Efficient Hi-C inversion facilitates chromatin folding mechanism discovery and structure prediction.

Genome-wide chromosome conformation capture (Hi-C) experiments have revealed many structural features of chromatin across multiple length scales. Further understanding genome organization requires relating these discoveries to the mechanisms that establish chromatin structures and reconstructing these structures in three dimensions, but both objectives are difficult to achieve with existing algorithms that are often computationally expensive. To alleviate this challenge, we present an algorithm that efficiently converts Hi-C data into contact energies, which measure the interaction strength between genomic loci brought into proximity.

OpenABC Enables Flexible, Simplified, and Efficient GPU Accelerated Simulations of Biomolecular Condensates.

Biomolecular condensates are important structures in various cellular processes but are challenging to study using traditional experimental techniques. In silico simulations with residue-level coarse-grained models strike a balance between computational efficiency and chemical accuracy. They could offer valuable insights by connecting the emergent properties of these complex systems with molecular sequences.

Efficient Hi-C inversion facilitates chromatin folding mechanism discovery and structure prediction.

Genome-wide chromosome conformation capture (Hi-C) experiments have revealed many structural features of chromatin across multiple length scales. Further understanding genome organization requires relating these discoveries to the mechanisms that establish chromatin structures and reconstructing these structures in three dimensions, but both objectives are difficult to achieve with existing algorithms that are often computationally expensive. To alleviate this challenge, we present an algorithm that efficiently converts Hi-C data into contact energies, which measure the interaction strength between genomic loci brought into proximity.

Contrastive Learning of Coarse-Grained Force Fields.

Coarse-grained models have proven helpful for simulating complex systems over long time scales to provide molecular insights into various processes. Methodologies for systematic parametrization of the underlying energy function or force field that describes the interactions among different components of the system are of great interest for ensuring simulation accuracy. We present a new method, potential contrasting, to enable efficient learning of force fields that can accurately reproduce the conformational distribution produced with all-atom simulations.

Generalizing the Discrete Gibbs Sampler-Based λ-Dynamics Approach for Multisite Sampling of Many Ligands.

In this work, the discrete λ variant of the Gibbs sampler-based λ-dynamics (-GSλD) method is developed to enable multiple functional group perturbations to be investigated at one or more sites of substitution off a common ligand core. The theoretical framework and special considerations for constructing discrete λ states for multisite -GSλD are presented. The precision and accuracy of the -GSλD method is evaluated with three test cases of increasing complexity.

DeepBAR: A Fast and Exact Method for Binding Free Energy Computation.

The fast and accurate calculation of standard binding free energy has many important applications. Existing methodologies struggle at balancing accuracy and efficiency. We introduce a new method to compute binding free energy using deep generative models and the Bennett acceptance ratio method (DeepBAR).

Stability and folding pathways of tetra-nucleosome from six-dimensional free energy surface.

The three-dimensional organization of chromatin is expected to play critical roles in regulating genome functions. High-resolution characterization of its structure and dynamics could improve our understanding of gene regulation mechanisms but has remained challenging. Using a near-atomistic model that preserves the chemical specificity of protein-DNA interactions at residue and base-pair resolution, we studied the stability and folding pathways of a tetra-nucleosome.

Computing Absolute Free Energy with Deep Generative Models.

Fast and accurate evaluation of free energy has broad applications from drug design to material engineering. Computing the absolute free energy is of particular interest since it allows the assessment of the relative stability between states without intermediates. Here, we introduce a general framework for calculating the absolute free energy of a state.

MiR-506-3p regulates autophagy and proliferation in post-burn skin fibroblasts through post-transcriptionally suppressing Beclin-1 expression.

MicroRNAs (miRNAs) is involved in diverse biological processes of cells including dermal fibroblasts that contributed to wound healing and resulted in keloid scarring. MiR-506-3p has been identified as a tumor suppressor or oncogene in fibroblasts of various cancers, while the role of miR-506-3p in regulating functions of post-burn dermal fibroblasts is poorly known. In this study, miR-506-3p was confirmed to be significantly downregulated in burned tissues and heat-stimulated dermal fibroblasts.

Accelerated CDOCKER with GPUs, Parallel Simulated Annealing, and Fast Fourier Transforms.

Fast Fourier transform (FFT)-based protein ligand docking together with parallel simulated annealing for both rigid and flexible receptor docking are implemented on graphical processing unit (GPU) accelerated platforms to significantly enhance the throughput of the CDOCKER and flexible CDOCKER - the docking algorithms in the CHARMM program for biomolecule modeling. The FFT-based approach for docking, first applied in protein-protein docking to efficiently search for the binding position and orientation of proteins, is adapted here to search ligand translational and rotational spaces given a ligand conformation in protein-ligand docking. Running on GPUs, our implementation of FFT docking in CDOCKER achieves a 15 000 fold speedup in the ligand translational and rotational space search in protein-ligand docking problems.

Deciphering protein evolution and fitness landscapes with latent space models.

Protein sequences contain rich information about protein evolution, fitness landscapes, and stability. Here we investigate how latent space models trained using variational auto-encoders can infer these properties from sequences. Using both simulated and real sequences, we show that the low dimensional latent space representation of sequences, calculated using the encoder model, captures both evolutionary and ancestral relationships between sequences.

Fast Solver for Large Scale Multistate Bennett Acceptance Ratio Equations.

The multistate Bennett acceptance ratio method (MBAR) and unbinned weighted histogram analysis method (UWHAM) are widely employed approaches to calculate relative free energies of multiple thermodynamic states that gain statistical precision by employing free energy contributions from configurations sampled at each of the simulated λ states. With the increasing availability of high throughput computing resources, a large number of configurations can be sampled from hundreds or even thousands of states. Combining sampled configurations from all states to calculate relative free energies requires the iterative solution of large scale MBAR/UWHAM equations.

CDOCKER and λ-dynamics for prospective prediction in D₃R Grand Challenge 2.

The opportunity to prospectively predict ligand bound poses and free energies of binding to the Farnesoid X Receptor in the D3R Grand Challenge 2 provided a useful exercise to evaluate CHARMM based docking (CDOCKER) and [Formula: see text]-dynamics methodologies for use in "real-world" applications in computer aided drug design. In addition to measuring their current performance, several recent methodological developments have been analyzed retrospectively to highlight best procedural practices in future applications. For pose prediction with CDOCKER, when the protein structure used for rigid receptor docking was close to the crystallographic holo structure, reliable poses were obtained.

Gibbs Sampler-Based λ-Dynamics and Rao-Blackwell Estimator for Alchemical Free Energy Calculation.

λ-dynamics is a generalized ensemble method for alchemical free energy calculations. In traditional λ-dynamics, the alchemical switch variable λ is treated as a continuous variable ranging from 0 to 1 and an empirical estimator is utilized to approximate the free energy. In the present article, we describe an alternative formulation of λ-dynamics that utilizes the Gibbs sampler framework, which we call Gibbs sampler-based λ-dynamics (GSLD).

Mechanism of Vps4 hexamer function revealed by cryo-EM.

Vps4 is a member of AAA ATPase (adenosine triphosphatase associated with diverse cellular activities) that operates as an oligomer to disassemble ESCRT-III (endosomal sorting complex required for transport III) filaments, thereby catalyzing the final step in multiple ESCRT-dependent membrane remodeling events. We used electron cryo-microscopy to visualize oligomers of a hydrolysis-deficient Vps4 (vacuolar protein sorting-associated protein 4) mutant in the presence of adenosine 5'-triphosphate (ATP). We show that Vps4 subunits assemble into an asymmetric hexameric ring following an approximate helical path that sequentially stacks substrate-binding loops along the central pore.

Improved prediction of RNA secondary structure by integrating the free energy model with restraints derived from experimental probing data.

Recently, several experimental techniques have emerged for probing RNA structures based on high-throughput sequencing. However, most secondary structure prediction tools that incorporate probing data are designed and optimized for particular types of experiments. For example, RNAstructure-Fold is optimized for SHAPE data, while SeqFold is optimized for PARS data.