Pharmacological hallmarks of allostery at the M4 muscarinic receptor elucidated through structure and dynamics.

Allosteric modulation of G protein-coupled receptors (GPCRs) is a major paradigm in drug discovery. Despite decades of research, a molecular-level understanding of the general principles that govern the myriad pharmacological effects exerted by GPCR allosteric modulators remains limited. The M4 muscarinic acetylcholine receptor (M4 mAChR) is a validated and clinically relevant allosteric drug target for several major psychiatric and cognitive disorders. In this study, we rigorously quantified the affinity, efficacy, and magnitude of modulation of two different positive allosteric modulators, LY2033298 (LY298) and VU0467154 (VU154), combined with the endogenous agonist acetylcholine (ACh) or the high-affinity agonist iperoxo (Ipx), at the human M4 mAChR. By determining the cryo-electron microscopy structures of the M4 mAChR, bound to a cognate Gi1 protein and in complex with ACh, Ipx, LY298-Ipx, and VU154-Ipx, and applying molecular dynamics simulations, we determine key molecular mechanisms underlying allosteric pharmacology. In addition to delineating the contribution of spatially distinct binding sites on observed pharmacology, our findings also revealed a vital role for orthosteric and allosteric ligand-receptor-transducer complex stability, mediated by conformational dynamics between these sites, in the ultimate determination of affinity, efficacy, cooperativity, probe dependence, and species variability. There results provide a holistic framework for further GPCR mechanistic studies and can aid in the discovery and design of future allosteric drugs.

Critical Non-Covalent Binding Intermediate for an Allosteric Covalent Inhibitor of SUMO E1.

Post-translational modifications by small ubiquitin-like modifiers (SUMOs) are dysregulated in many types of cancers. The SUMO E1 enzyme has recently been suggested as a new immuno-oncology target. COH000 was recently identified as a highly specific allosteric covalent inhibitor of SUMO E1. However, a marked discrepancy was found between the X-ray structure of the covalent COH000-bound SUMO E1 complex and the available structure-activity relationship (SAR) data of inhibitor analogues due to unresolved noncovalent protein-ligand interactions. Here, we have investigated noncovalent interactions between COH000 and SUMO E1 during inhibitor dissociation through novel Ligand Gaussian accelerated molecular dynamics (LiGaMD) simulations. Our simulations have identified a critical low-energy non-covalent binding intermediate conformation of COH000 that agreed excellently with published and new SAR data of the COH000 analogues, which were otherwise inconsistent with the X-ray structure. Altogether, our biochemical experiments and LiGaMD simulations have uncovered a critical non-covalent binding intermediate during allosteric inhibition of the SUMO E1 complex.

Effects of presenilin-1 familial Alzheimer's disease mutations on γ-secretase activation for cleavage of amyloid precursor protein.

Presenilin-1 (PS1) is the catalytic subunit of γ-secretase which cleaves within the transmembrane domain of over 150 peptide substrates. Dominant missense mutations in PS1 cause early-onset familial Alzheimer's disease (FAD); however, the exact pathogenic mechanism remains unknown. Here we combined Gaussian accelerated molecular dynamics (GaMD) simulations and biochemical experiments to determine the effects of six representative PS1 FAD mutations (P117L, I143T, L166P, G384A, L435F, and L286V) on the enzyme-substrate interactions between γ-secretase and amyloid precursor protein (APP). Biochemical experiments showed that all six PS1 FAD mutations rendered γ-secretase less active for the endoproteolytic (ε) cleavage of APP. Distinct low-energy conformational states were identified from the free energy profiles of wildtype and PS1 FAD-mutant γ-secretase. The P117L and L286V FAD mutants could still sample the "Active" state for substrate cleavage, but with noticeably reduced conformational space compared with the wildtype. The other mutants hardly visited the "Active" state. The PS1 FAD mutants were found to reduce γ-secretase proteolytic activity by hindering APP residue L49 from proper orientation in the active site and/or disrupting the distance between the catalytic aspartates. Therefore, our findings provide mechanistic insights into how PS1 FAD mutations affect structural dynamics and enzyme-substrate interactions of γ-secretase and APP.

Mechanism of Tripeptide Trimming of Amyloid ß-Peptide 49 by γ-Secretase.

The membrane-embedded γ-secretase complex processively cleaves within the transmembrane domain of amyloid precursor protein (APP) to produce 37-to-43-residue amyloid ß-peptides (Aß) of Alzheimer's disease (AD). Despite its importance in pathogenesis, the mechanism of processive proteolysis by γ-secretase remains poorly understood. Here, mass spectrometry and Western blotting were used to quantify the efficiency of tripeptide trimming of wild-type (WT) and familial AD (FAD) mutant Aß49. In comparison to WT Aß49, the efficiency of tripeptide trimming was similar for the I45F, A42T, and V46F Aß49 FAD mutants but substantially diminished for the I45T and T48P mutants. In parallel with biochemical experiments, all-atom simulations using a novel peptide Gaussian accelerated molecular dynamics (Pep-GaMD) method were applied to investigate the tripeptide trimming of Aß49 by γ-secretase. The starting structure was the active γ-secretase bound to Aß49 and APP intracellular domain (AICD), as generated from our previous study that captured the activation of γ-secretase for the initial endoproteolytic cleavage of APP (Bhattarai, A., ACS Cent. Sci. 2020, 6, 969-983). Pep-GaMD simulations captured remarkable structural rearrangements of both the enzyme and substrate, in which hydrogen-bonded catalytic aspartates and water became poised for tripeptide trimming of Aß49 to Aß46. These structural changes required a positively charged N-terminus of endoproteolytic coproduct AICD, which could dissociate during conformational rearrangements of the protease and Aß49. The simulation findings were highly consistent with biochemical experimental data. Taken together, our complementary biochemical experiments and Pep-GaMD simulations have enabled elucidation of the mechanism of tripeptide trimming of Aß49 by γ-secretase.

Molecular Simulations and Drug Discovery of Adenosine Receptors.

G protein-coupled receptors (GPCRs) represent the largest family of human membrane proteins. Four subtypes of adenosine receptors (ARs), the A1AR, A2AAR, A2BAR and A3AR, each with a unique pharmacological profile and distribution within the tissues in the human body, mediate many physiological functions and serve as critical drug targets for treating numerous human diseases including cancer, neuropathic pain, cardiac ischemia, stroke and diabetes. The A1AR and A3AR preferentially couple to the Gi/o proteins, while the A2AAR and A2BAR prefer coupling to the Gs proteins. Adenosine receptors were the first subclass of GPCRs that had experimental structures determined in complex with distinct G proteins. Here, we will review recent studies in molecular simulations and computer-aided drug discovery of the adenosine receptors and also highlight their future research opportunities.

GLOW: A Workflow Integrating Gaussian-Accelerated Molecular Dynamics and Deep Learning for Free Energy Profiling.

We introduce a Gaussian-accelerated molecular dynamics (GaMD), deep learning (DL), and free energy profiling workflow (GLOW) to predict molecular determinants and map free energy landscapes of biomolecules. All-atom GaMD-enhanced sampling simulations are first performed on biomolecules of interest. Structural contact maps are then calculated from GaMD simulation frames and transformed into images for building DL models using a convolutional neural network. Important structural contacts are further determined from DL models of attention maps of the structural contact gradients, which allow us to identify the system reaction coordinates. Finally, free energy profiles are calculated for the selected reaction coordinates through energetic reweighting of the GaMD simulations. We have also successfully demonstrated GLOW for the characterization of activation and allosteric modulation of a G protein-coupled receptor, using the adenosine A1 receptor (A1AR) as a model system. GLOW findings are highly consistent with previous experimental and computational studies of the A1AR, while also providing further mechanistic insights into the receptor function. In summary, GLOW provides a systematic approach to mapping free energy landscapes of biomolecules. The GLOW workflow and its user manual can be downloaded at http://miaolab.org/GLOW.

Binding Analysis Using Accelerated Molecular Dynamics Simulations and Future Perspectives.

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.

Challenges and frontiers of computational modelling of biomolecular recognition.

Biomolecular recognition including binding of small molecules, peptides and proteins to their target receptors plays a key role in cellular function and has been targeted for therapeutic drug design. However, the high flexibility of biomolecules and slow binding and dissociation processes have presented challenges for computational modeling. Here, we review the challenges and computational approaches developed to characterize biomolecular binding, including molecular docking, Molecular Dynamics (MD) simulations (especially enhanced sampling) and Machine Learning. Further improvements are still needed in order to accurately and efficiently characterize binding structures, mechanisms, thermodynamics and kinetics of biomolecules in the future.

Gaussian accelerated molecular dynamics (GaMD): principles and applications.

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.

Positive allosteric mechanisms of adenosine A1 receptor-mediated analgesia.

The adenosine A1 receptor (A1R) is a promising therapeutic target for non-opioid analgesic agents to treat neuropathic pain1,2. However, development of analgesic orthosteric A1R agonists has failed because of a lack of sufficient on-target selectivity as well as off-tissue adverse effects3. Here we show that [2-amino-4-(3,5-bis(trifluoromethyl)phenyl)thiophen-3-yl)(4-chlorophenyl)methanone] (MIPS521), a positive allosteric modulator of the A1R, exhibits analgesic efficacy in rats in vivo through modulation of the increased levels of endogenous adenosine that occur in the spinal cord of rats with neuropathic pain. We also report the structure of the A1R co-bound to adenosine, MIPS521 and a Gi2 heterotrimer, revealing an extrahelical lipid-detergent-facing allosteric binding pocket that involves transmembrane helixes 1, 6 and 7. Molecular dynamics simulations and ligand kinetic binding experiments support a mechanism whereby MIPS521 stabilizes the adenosine-receptor-G protein complex. This study provides proof of concept for structure-based allosteric drug design of non-opioid analgesic agents that are specific to disease contexts.

Mechanism of Ligand Recognition by Human ACE2 Receptor.

Angiotensin converting enzyme 2 (ACE2) plays a key role in renin-angiotensin system regulation and amino acid homeostasis. Human ACE2 acts as the receptor for severe acute respiratory syndrome coronaviruses SARS-CoV and SARS-CoV-2. ACE2 is also widely expressed in epithelial cells of the lungs, heart, kidney, and pancreas. It is considered an important drug target for treating SARS-CoV-2 as well as pulmonary diseases, heart failure, hypertension, renal diseases, and diabetes. Despite the critical importance, the mechanism of ligand binding to the human ACE2 receptor remains unknown. Here, we have addressed this challenge through all-atom simulations using a novel ligand Gaussian accelerated molecular dynamics (LiGaMD) method. Microsecond time scale LiGaMD simulations have unprecedentedly captured multiple times of spontaneous binding and unbinding of a potent inhibitor MLN-4760 in the ACE2 receptor. With ligand far away in the unbound state, the ACE2 receptor samples distinct Open, Partially Open, Closed, and Fully Closed conformations. Upon ligand binding to the active site, conformational ensemble of the ACE2 receptor is biased toward the Closed state as observed in the X-ray experimental structure. The LiGaMD simulations thus suggest a conformational selection mechanism for ligand recognition by the highly flexible ACE2 receptor, which is expected to facilitate rational drug design targeting human ACE2 against coronaviruses and other related human diseases.

Mechanism of ligand recognition by human ACE2 receptor.

Angiotensin converting enzyme 2 (ACE2) plays a key role in renin-angiotensin system regulation and amino acid homeostasis. Human ACE2 acts as the receptor for severe acute respiratory syndrome coronaviruses SARS-CoV and SARS-CoV-2. ACE2 is also widely expressed in epithelial cells of lungs, heart, kidney and pancreas. It is considered an important drug target for treating SARS-CoV-2, as well as pulmonary diseases, heart failure, hypertension, renal diseases and diabetes. Despite the critical importance, the mechanism of ligand binding to the human ACE2 receptor remains unknown. Here, we address this challenge through all-atom simulations using a novel ligand Gaussian accelerated molecular dynamics (LiGaMD) method. Microsecond LiGaMD simulations have successfully captured both binding and unbinding of the MLN-4760 inhibitor in the ACE2 receptor. In the ligand unbound state, the ACE2 receptor samples distinct Open, Partially Open and Closed conformations. Ligand binding biases the receptor conformational ensemble towards the Closed state. The LiGaMD simulations thus suggest a conformational selection mechanism for ligand recognition by the ACE2 receptor. Our simulation findings are expected to facilitate rational drug design of ACE2 against coronaviruses and other related human diseases.

Ligand Gaussian Accelerated Molecular Dynamics (LiGaMD): Characterization of Ligand Binding Thermodynamics and Kinetics.

Calculations of ligand binding free energies and kinetic rates are important for drug design. However, such tasks have proven challenging in computational chemistry and biophysics. To address this challenge, we have developed a new computational method, ligand Gaussian accelerated molecular dynamics (LiGaMD), which selectively boosts the ligand nonbonded interaction potential energy based on the Gaussian accelerated molecular dynamics (GaMD) enhanced sampling technique. Another boost potential could be applied to the remaining potential energy of the entire system in a dual-boost algorithm (LiGaMD_Dual) to facilitate ligand binding. LiGaMD has been demonstrated on host-guest and protein-ligand binding model systems. Repetitive guest binding and unbinding in the ß-cyclodextrin host were observed in hundreds-of-nanosecond LiGaMD_Dual simulations. The calculated guest binding free energies agreed excellently with experimental data with <1.0 kcal/mol errors. Compared with converged microsecond-time scale conventional molecular dynamics simulations, the sampling errors of LiGaMD_Dual simulations were also <1.0 kcal/mol. Accelerations of ligand kinetic rate constants in LiGaMD simulations were properly estimated using Kramers' rate theory. Furthermore, LiGaMD allowed us to capture repetitive dissociation and binding of the benzamidine inhibitor in trypsin within 1 µs simulations. The calculated ligand binding free energy and kinetic rate constants compared well with the experimental data. In summary, LiGaMD provides a powerful enhanced sampling approach for characterizing ligand binding thermodynamics and kinetics simultaneously, which is expected to facilitate computer-aided drug design.

Mechanisms of γ-Secretase Activation and Substrate Processing.

Amyloid ß-peptide, the principal component of characteristic cerebral plaques of Alzheimer's disease (AD), is produced through intramembrane proteolysis of the amyloid precursor protein (APP) by γ-secretase. Despite the importance in the pathogenesis of AD, the mechanisms of intramembrane proteolysis and substrate processing by γ-secretase remain poorly understood. Here, complementary all-atom simulations using a robust Gaussian accelerated molecular dynamics (GaMD) method and biochemical experiments were combined to investigate substrate processing of wildtype and mutant APP by γ-secretase. The GaMD simulations captured spontaneous activation of γ-secretase, with hydrogen bonded catalytic aspartates and water poised for proteolysis of APP at the ε cleavage site. Furthermore, GaMD simulations revealed that familial AD mutations I45F and T48P enhanced the initial ε cleavage between residues Leu49-Val50, while M51F mutation shifted the ε cleavage site to the amide bond between Thr48-Leu49. Detailed analysis of the GaMD simulations allowed us to identify distinct low-energy conformational states of γ-secretase, different secondary structures of the wildtype and mutant APP substrate, and important active-site subpockets for catalytic function of the enzyme. The simulation findings were highly consistent with experimental analyses of APP proteolytic products using mass spectrometry and Western blotting. Taken together, the GaMD simulations and biochemical experiments have enabled us to elucidate the mechanisms of γ-secretase activation and substrate processing, which should facilitate rational computer-aided drug design targeting this functionally important enzyme.

Retrospective ensemble docking of allosteric modulators in an adenosine G-protein-coupled receptor.

BACKGROUND: Ensemble docking has proven useful in drug discovery and development. It increases the hit rate by incorporating receptor flexibility into molecular docking as demonstrated on important drug targets including G-protein-coupled receptors (GPCRs). Adenosine A1 receptor (A1AR) is a key GPCR that has been targeted for treating cardiac ischemia-reperfusion injuries, neuropathic pain and renal diseases. Development of allosteric modulators, compounds binding to distinct and less conserved GPCR target sites compared with agonists and antagonists, has attracted increasing interest for designing selective drugs of the A1AR. Despite significant advances, more effective approaches are needed to discover potent and selective allosteric modulators of the A1AR. METHODS: Ensemble docking that integrates Gaussian accelerated molecular dynamic (GaMD) simulations and molecular docking using Autodock has been implemented for retrospective docking of known positive allosteric modulators (PAMs) in the A1AR. RESULTS: Ensemble docking outperforms docking of the receptor cryo-EM structure. The calculated docking enrichment factors (EFs) and the area under the receiver operating characteristic curves (AUC) are significantly increased. CONCLUSIONS: Receptor ensembles generated from GaMD simulations are able to increase the success rate of discovering PAMs of A1AR. It is important to account for receptor flexibility through GaMD simulations and flexible docking. GENERAL SIGNIFICANCE: Ensemble docking is a promising approach for drug discovery targeting flexible receptors.

G-Protein-Coupled Receptor-Membrane Interactions Depend on the Receptor Activation State.

G-protein-coupled receptors (GPCRs) are the largest family of human membrane proteins and serve as primary targets of approximately one-third of currently marketed drugs. In particular, adenosine A1 receptor (A1 AR) is an important therapeutic target for treating cardiac ischemia-reperfusion injuries, neuropathic pain, and renal diseases. As a prototypical GPCR, the A1 AR is located within a phospholipid membrane bilayer and transmits cellular signals by changing between different conformational states. It is important to elucidate the lipid-protein interactions in order to understand the functional mechanism of GPCRs. Here, all-atom simulations using a robust Gaussian accelerated molecular dynamics (GaMD) method were performed on both the inactive (antagonist bound) and active (agonist and G-protein bound) A1 AR, which was embedded in a 1-palmitoyl-2-oleoyl-glycero-3-phosphocholine (POPC) lipid bilayer. In the GaMD simulations, the membrane lipids played a key role in stabilizing different conformational states of the A1 AR. Our simulations further identified important regions of the receptor that interacted distinctly with the lipids in highly correlated manner. Activation of the A1 AR led to differential dynamics in the upper and lower leaflets of the lipid bilayer. In summary, GaMD enhanced simulations have revealed strongly coupled dynamics of the GPCR and lipids that depend on the receptor activation state. © 2019 Wiley Periodicals, Inc.

Structural Basis for Binding of Allosteric Drug Leads in the Adenosine A1 Receptor.

Despite intense interest in designing positive allosteric modulators (PAMs) as selective drugs of the adenosine A1 receptor (A1AR), structural binding modes of the receptor PAMs remain unknown. Using the first X-ray structure of the A1AR, we have performed all-atom simulations using a robust Gaussian accelerated molecular dynamics (GaMD) technique to determine binding modes of the A1AR allosteric drug leads. Two prototypical PAMs, PD81723 and VCP171, were selected. Each PAM was initially placed at least 20 Å away from the receptor. Extensive GaMD simulations using the AMBER and NAMD simulation packages at different acceleration levels captured spontaneous binding of PAMs to the A1AR. The simulations allowed us to identify low-energy binding modes of the PAMs at an allosteric site formed by the receptor extracellular loop 2 (ECL2), which are highly consistent with mutagenesis experimental data. Furthermore, the PAMs stabilized agonist binding in the receptor. In the absence of PAMs at the ECL2 allosteric site, the agonist sampled a significantly larger conformational space and even dissociated from the A1AR alone. In summary, the GaMD simulations elucidated structural binding modes of the PAMs and provided important insights into allostery in the A1AR, which will greatly facilitate the receptor structure-based drug design.

Gaussian accelerated molecular dynamics for elucidation of drug pathways.

INTRODUCTION: Understanding pathways and mechanisms of drug binding to receptors is important for rational drug design. Remarkable advances in supercomputing and methodological developments have opened a new era for application of computer simulations in predicting drug-receptor interactions at an atomistic level. Gaussian accelerated molecular dynamics (GaMD) is a computational enhanced sampling technique that works by adding a harmonic boost potential to reduce energy barriers. GaMD enables free energy calculations without the requirement of predefined collective variables. GaMD has proven useful in biomolecular simulations, in particular, the prediction of drug-receptor interactions. Areas covered: Herein, the authors review recent GaMD simulation studies that elucidated pathways of drug binding to proteins including the G-protein-coupled receptors and HIV protease. Expert opinion: GaMD is advantageous for enhanced simulations of, amongst many biological processes, drug binding to target receptors. Compared with conventional molecular dynamics, GaMD speeds up biomolecular simulations by orders of magnitude. GaMD enables routine drug binding simulations using personal computers with GPUs or common computing clusters. GaMD and, more broadly, enhanced sampling simulations are expected to dramatically increase our capabilities to determine the mechanisms of drug binding to a wide range of receptors in the near future. This will greatly facilitate computer-aided drug design.