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-

Extracting binding energies and binding modes from biomolecular simulations of fragment binding to endothiapepsin

Fragment Databases from Screened Ligands for Drug Discovery (FDSL-DD).

Fragment-based Drug Discovery and Molecular Docking in Drug Design.

Predicting Fragment Binding Modes Using Customized Lennard-Jones Potentials in Short Molecular Dynamics Simulations

Streamlining Computational Fragment-Based Drug Discovery Through Evolutionary Optimization Informed by Ligand-Based Virtual Prescreening

Binding Energy Distribution Analysis Method: Hamiltonian Replica Exchange with Torsional Flattening for Binding Mode Prediction and Binding Free Energy Estimation

Quantum Fragment Based ab Initio Molecular Dynamics for Proteins.

Binding Analysis Using Accelerated Molecular Dynamics Simulations and Future Perspectives

A Force Balanced Fragmentation Method for Ab Initio Molecular Dynamic Simulation of Protein

Biomolecular conformational changes and ligand binding: from kinetics to thermodynamics.

Fast and accurate predictions of binding free energies using MM-PBSA and MM-GBSA

FRAGSITE: A Fragment-Based Approach for Virtual Ligand Screening

Fragmenstein: predicting protein-ligand structures of compounds derived from known crystallographic fragment hits using a strict conserved-binding–based methodology

Binding free energy of protein/ligand complexes calculated using dissociation Parallel Cascade Selection Molecular Dynamics and Markov state model

Protein-peptide association kinetics beyond the seconds timescale from atomistic simulations

Hit Identification and Binding Mode Predictions by Rigorous Free Energy Simulations

Comparison of Methodologies for Absolute Binding Free Energy Calculations of Ligands to Intrinsically Disordered Proteins

Efficient large-scale exploration of fragment hit progression by exploiting binding-site purification of actives (B-SPA) through combining multi-step array synthesis and HT crystallography.

Predictions Of Binding Of A Diverse Set Of Ligands To Gelatinase-A By A Combination Of Molecular Dynamics And Continuum Solvent Models

MDFit: automated molecular simulations workflow enables high throughput assessment of ligands-protein dynamics