Gaussian Accelerated Molecular Dynamics in OpenMM

Multiscale simulation of fluids: coupling molecular and continuum

Computer simulation is an important tool for scientific progress, especially when lab experiments are either extremely costly and difficult or lack the required resolution. However, all of the simulation methods come with limitations. In molecular dynamics (MD) simulation, the length and time scales that can be captured are limited, while computational fluid dynamics (CFD) methods are built on a range of assumptions, from the continuum hypothesis itself, to a variety of closure assumptions. To address these issues, the coupling of different methodologies provides a way to retain the best of both methods. Here, we provide a perspective on multiscale simulation based on the coupling of MD and CFD with each a distinct part of the same simulation domain. This style of coupling allows molecular detail to be present only where it is needed, so CFD can model larger scales than possible with MD alone. We present a unified perspective of the literature, showing the links between the two main types of coupling, state and flux, and discuss the varying assumptions in their use. A unique challenge in such coupled simulation is obtaining averages and constraining local parts of a molecular simulation. We highlight that incorrect localisation has resulted in an error in the literature. We then finish with some applications, focused on the simulation of fluids. Thus, we hope to motivate further research in this exciting area with applications across the spectrum of scientific disciplines.

When MELD Meets GaMD: Accelerating Biomolecular Landscape Exploration

We introduce Gaussian accelerated MELD (GaMELD) as a new method for exploring the energy landscape of biomolecules. GaMELD combines the strengths of Gaussian accelerated molecular dynamics (GaMD) and modeling employing limited data (MELD) to navigate complex energy landscapes. MELD uses replica-exchange molecular simulations to integrate limited and uncertain data into simulations via Bayesian inference. MELD has been successfully applied to problems of structure prediction like protein folding and complex structure prediction. However, the computational cost for MELD simulations has limited its broader applicability. The synergy of GaMD and MELD surmounts this limitation efficiently sampling the energy landscape at a lower computational cost (reducing the computational cost by a factor of 2 to six). Effectively, GaMD is used to shift energy distributions along replicas to increase the overlap in energy distributions across replicas, facilitating a random walk in replica space. We tested GaMELD on a benchmark set of 12 small proteins that have been previously studied through MELD and conventional MD. GaMELD consistently achieves accurate predictions with fewer replicas. By increasing the efficacy of replica exchange, GaMELD effectively accelerates convergence in the conformational space, enabling improved sampling across a diverse set of systems.

Structural basis of antibody inhibition and chemokine activation of the human CC chemokine receptor 8

The C-C motif chemokine receptor 8 (CCR8) is a class A G-protein coupled receptor that has emerged as a promising therapeutic target in cancer. Targeting CCR8 with an antibody has appeared to be an attractive therapeutic approach, but the molecular basis for chemokine-mediated activation and antibody-mediated inhibition of CCR8 are not fully elucidated. Here, we obtain an antagonist antibody against human CCR8 and determine structures of CCR8 in complex with either the antibody or the endogenous agonist ligand CCL1. Our studies reveal characteristic antibody features allowing recognition of the CCR8 extracellular loops and CCL1-CCR8 interaction modes that are distinct from other chemokine receptor - ligand pairs. Informed by these structural insights, we demonstrate that CCL1 follows a two-step, two-site binding sequence to CCR8 and that antibody-mediated inhibition of CCL1 signaling can occur by preventing the second binding event. Together, our results provide a detailed structural and mechanistic framework of CCR8 activation and inhibition that expands our molecular understanding of chemokine - receptor interactions and offers insight into the development of therapeutic antibodies targeting chemokine GPCRs.

Flexible Gaussian Accelerated Molecular Dynamics to Enhance Biological Sampling

Molecular dynamics simulations often struggle to obtain sufficient sampling to study complex molecular events due to high energy barriers separating the minima of interest. Multiple enhanced sampling techniques have been developed and improved over the years to tackle this issue. Gaussian accelerated molecular dynamics (GaMD) is a recently developed enhanced sampling technique that works by adding a biasing potential, lifting the energy landscape up, and decreasing the height of its barriers. GaMD allows one to increase the sampling of events of interest without the need of a priori knowledge of the system or the relevant coordinates. All required acceleration parameters can be obtained from a previous search run. Upon its development, several improvements for the methodology have been proposed, among them selective GaMD in which the boosting potential is selectively applied to the region of interest. There are currently four selective GaMD methods that have shown promising results. However, all of these methods are constrained on the number, location, and scenarios in which this selective boosting potential can be applied to ligands, peptides, or protein-protein interactions. In this work, we showcase a GROMOS implementation of the GaMD methodology with a fully flexible selective GaMD approach that allows the user to define, in a straightforward way, multiple boosting potentials for as many regions as desired. We show and analyze the advantages of this flexible selective approach on two previously used test systems, the alanine dipeptide and the chignolin peptide, and extend these examples to study its applicability and potential to study conformational changes of glycans and glycosylated proteins.

Deep Boosted Molecular Dynamics (DBMD): Accelerating molecular simulations with Gaussian boost potentials generated using probabilistic Bayesian deep neural network

We have developed a new Deep Boosted Molecular Dynamics (DBMD) method. Probabilistic Bayesian neural network models were implemented to construct boost potentials that exhibit Gaussian distribution with minimized anharmonicity, thereby allowing for accurate energetic reweighting and enhanced sampling of molecular simulations. DBMD was demonstrated on model systems of alanine dipeptide and the fast-folding protein and RNA structures. For alanine dipeptide, 30ns DBMD simulations captured up to 83-125 times more backbone dihedral transitions than 1μs conventional molecular dynamics (cMD) simulations and were able to accurately reproduce the original free energy profiles. Moreover, DBMD sampled multiple folding and unfolding events within 300ns simulations of the chignolin model protein and identified low-energy conformational states comparable to previous simulation findings. Finally, DBMD captured a general folding pathway of three hairpin RNAs with the GCAA, GAAA, and UUCG tetraloops. Based on Deep Learning neural network, DBMD provides a powerful and generally applicable approach to boosting biomolecular simulations. DBMD is available with open source in OpenMM at https://github.com/MiaoLab20/DBMD/.

Predicting Biomolecular Binding Kinetics: A Review

Biomolecular binding kinetics including the association (k_on) and dissociation (k_off) rates are critical parameters for therapeutic design of small-molecule drugs, peptides, and antibodies. Notably, the drug molecule residence time or dissociation rate has been shown to correlate with their efficacies better than binding affinities. A wide range of modeling approaches including quantitative structure-kinetic relationship models, Molecular Dynamics simulations, enhanced sampling, and Machine Learning has been developed to explore biomolecular binding and dissociation mechanisms and predict binding kinetic rates. Here, we review recent advances in computational modeling of biomolecular binding kinetics, with an outlook for future improvements.