Gaussian accelerated molecular dynamics (GaMD) is a robust computational method for simultaneous unconstrained enhanced sampling and free energy calculations of biomolecules. It works by adding a harmonic boost potential to smooth biomolecular potential energy surface and reduce energy barriers. GaMD greatly accelerates biomolecular simulations by orders of magnitude. Without the need to set predefined reaction coordinates or collective variables, GaMD provides unconstrained enhanced sampling and is advantageous for simulating complex biological processes. The GaMD boost potential exhibits a Gaussian distribution, thereby allowing for energetic reweighting via cumulant expansion to the second order (i.e., “Gaussian approximation”). This leads to accurate reconstruction of free energy landscapes of biomolecules. Hybrid schemes with other enhanced sampling methods, such as the replica‐exchange GaMD (rex‐GaMD) and replica‐exchange umbrella sampling GaMD (GaREUS), have also been introduced, further improving sampling and free energy calculations. Recently, new “selective GaMD” algorithms including the Ligand GaMD (LiGaMD) and Peptide GaMD (Pep‐GaMD) enabled microsecond simulations to capture repetitive dissociation and binding of small‐molecule ligands and highly flexible peptides. The simulations then allowed highly efficient quantitative characterization of the ligand/peptide binding thermodynamics and kinetics. Taken together, GaMD and its innovative variants are applicable to simulate a wide variety of biomolecular dynamics, including protein folding, conformational changes and allostery, ligand binding, peptide binding, protein–protein/nucleic acid/carbohydrate interactions, and carbohydrate/nucleic acid interactions. In this review, we present principles of the GaMD algorithms and recent applications in biomolecular simulations and drug design. This article is categorized under: Structure and Mechanism > Computational Biochemistry and Biophysics Molecular and Statistical Mechanics > Molecular Dynamics and Monte‐Carlo Methods Molecular and Statistical Mechanics > Free Energy Methods

Background: Chemokine GPCRs play key roles in biology and medicine. Particularly, CXCR4 promotes cancer metastasis and facilitate HIV entry into host cells. Plerixafor (PLX) is a CXCR4 drug, but the pathway and binding site of PLX in CXCR4 remain unknown. Results & methodology: We have performed molecular docking and all-atom simulations using Gaussian accelerated molecular dynamics (GaMD), which are consistent with previous mutation experiments, suggesting that PLX binds to the orthosteric site of CXCR4 as an antagonist. The GaMD simulations further revealed an intermediate allosteric binding site at the extracellular mouth of CXCR4. Conclusion: The newly identified allosteric site can be targeted for novel drug design targeting CXCR4 and other chemokine receptors.

No abstract

Significance Mutations of polycystin-1 (PC1) are the major cause (85% of cases) of autosomal dominant polycystic kidney disease (ADPKD), which is the fourth leading cause of kidney failure. PC1 is thought to function as an atypical G protein-coupled receptor, yet the mechanism by which PC1 regulates G-protein signaling remains poorly understood. A significant portion of ADPKD mutations of PC1 encode a protein with defects in maturation or reduced function that may be amenable to functional rescue. In this work, we have combined complementary biochemical and cellular assay experiments and accelerated molecular simulations, which revealed an allosteric transduction pathway in activation of the PC1 C-terminal fragment. Our findings will facilitate future rational drug design efforts targeting the PC1 signaling function.

No abstract