Binding Analysis Using Accelerated Molecular Dynamics Simulations and Future Perspectives
Shristi Pawnikar,Apurba Bhattarai,Jinan Wang,Yinglong Miao
DOI: https://doi.org/10.2147/AABC.S247950
2022-01-06
Advances and Applications in Bioinformatics and Chemistry
Abstract:Shristi Pawnikar, Apurba Bhattarai, Jinan Wang, Yinglong Miao Center for Computational Biology and Department of Molecular Biosciences, University of Kansas, Lawrence, KS, 66047, USA Correspondence: Yinglong Miao Center for Computational Biology, University of Kansas, 2030 Becker Drive, Room 200D MRB, Lawrence, KS, 66047, USA Email Biomolecular recognition such as binding of small molecules, nucleic acids, peptides and proteins to their target receptors plays key roles in cellular function and has been targeted for therapeutic drug design. Molecular dynamics (MD) is a computational approach to analyze these binding processes at an atomistic level, which provides valuable understandings of the mechanisms of biomolecular recognition. However, the rather slow biomolecular binding events often present challenges for conventional MD (cMD), due to limited simulation timescales (typically over hundreds of nanoseconds to tens of microseconds). In this regard, enhanced sampling methods, particularly accelerated MD (aMD), have proven useful to bridge the gap and enable all-atom simulations of biomolecular binding events. Here, we will review the recent method developments of Gaussian aMD (GaMD), ligand GaMD (LiGaMD) and peptide GaMD (Pep-GaMD), which have greatly expanded our capabilities to simulate biomolecular binding processes. Spontaneous binding of various biomolecules to their receptors has been successfully simulated by GaMD. Microsecond LiGaMD and Pep-GaMD simulations have captured repetitive binding and dissociation of small-molecule ligands and highly flexible peptides, and thus enabled ligand/peptide binding thermodynamics and kinetics calculations. We will also present relevant application studies in simulations of important drug targets and future perspectives for rational computer-aided drug design. Keywords: biomolecular recognition, accelerated molecular dynamics, thermodynamics, kinetics, drug design Cellular functions such as gene translation, protein folding, ligand binding and signaling often require conformational transitions of biological macromolecules such as proteins, nucleic acids and lipids into various low-energy states. 1–4 Free energy landscapes of biomolecules essentially determine their structure, dynamics and functions. 5,6 It is important to comprehend the dynamic nature of biomolecules. Experimental techniques including nuclear magnetic resonance (NMR), X-ray crystallography and cryo-electron microscopy (cryo-EM) have greatly advanced structural characterizations of the biomolecular binding processes, but these efforts are time-consuming and very expensive. Moreover, these techniques often provide a rather static picture of the biomolecular interactions and conformations, whereas biomolecules are dynamic in vivo to perform the necessary cellular functions. It remains challenging and critical for us to understand the mechanisms, thermodynamics and kinetics of biomolecular recognition. Molecular dynamics (MD) is a powerful "computational microscope" that helps us to visualize atomistic time evolution of biomolecules. 7 Remarkable advances have been made in computing hardware (eg, Anton supercomputers and GPUs) and algorithm developments over the past several decades for running longer and cheaper MD simulations. 8–12 With microsecond-timescale conventional MD (cMD) simulations, spontaneous ligand/drug binding was successfully captured to the Src protein kinase, 13 the β 1 - and β 2 -adrenergic receptors (β 1 AR and β 1 AR) 14 and the M 2 and M 3 muscarinic receptors 15 by the DE Shaw Research Group using the Anton supercomputer, for which the simulation predicted ligand binding poses closely matched those determined in the X-ray crystal structures. Moreover, hundreds-of-nanosecond cMD simulations were combined in Markov State Models to characterize the benzamidine ligand binding to the trypsin model protein by the De Fabritiis group. 16 In addition, microsecond-timescale cMD simulations have been successful to capture peptide binding in several peptide-protein complexes, including proline-rich motif (PRM)-SH3, 17 BAD-BclX L and p53-MDM2. 18 The "folding-upon-binding" mechanism of peptides was investigated using a measles virus nucleoprotein fragment and X domain of its phosphoprotein complex. 19 Anton supercomputer was used to run ~200 μs MD simulations at 400 K elevated temperature. The simulations could capture 70 binding and unbinding events of the peptide fragment to the protein domain. Therefore, cMD has been successful to capture binding of ligands and peptides on the microsecond timescales. However, cMD with limited simulation timescales still lacks the power to study slower biological p -Abstract Truncated-
English Else