1
|
Wang J, Miao Y. Ligand Gaussian Accelerated Molecular Dynamics 3 (LiGaMD3): Improved Calculations of Binding Thermodynamics and Kinetics of Both Small Molecules and Flexible Peptides. J Chem Theory Comput 2024. [PMID: 39002136 DOI: 10.1021/acs.jctc.4c00502] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 07/15/2024]
Abstract
Binding thermodynamics and kinetics play critical roles in drug design. However, it has proven challenging to efficiently predict ligand binding thermodynamics and kinetics of small molecules and flexible peptides using conventional molecular dynamics (cMD), due to limited simulation time scales. Based on our previously developed ligand Gaussian accelerated molecular dynamics (LiGaMD) method, we present a new approach, termed "LiGaMD3″, in which we introduce triple boosts into three individual energy terms that play important roles in small-molecule/peptide dissociation, rebinding, and system conformational changes to improve the sampling efficiency of small-molecule/peptide interactions with target proteins. To validate the performance of LiGaMD3, MDM2 bound by a small molecule (Nutlin 3) and two highly flexible peptides (PMI and P53) were chosen as the model systems. LiGaMD3 could efficiently capture repetitive small-molecule/peptide dissociation and binding events within 2 μs simulations. The predicted binding kinetic constant rates and free energies from LiGaMD3 were in agreement with the available experimental values and previous simulation results. Therefore, LiGaMD3 provides a more general and efficient approach to capture dissociation and binding of both small-molecule ligands and flexible peptides, allowing for accurate prediction of their binding thermodynamics and kinetics.
Collapse
Affiliation(s)
- Jinan Wang
- Computational Medicine Program and Department of Pharmacology, University of North Carolina-Chapel Hill, Chapel Hill, North Carolina 27599, United States
| | - Yinglong Miao
- Computational Medicine Program and Department of Pharmacology, University of North Carolina-Chapel Hill, Chapel Hill, North Carolina 27599, United States
| |
Collapse
|
2
|
Sarkar DK, Surpeta B, Brezovsky J. Incorporating Prior Knowledge in the Seeds of Adaptive Sampling Molecular Dynamics Simulations of Ligand Transport in Enzymes with Buried Active Sites. J Chem Theory Comput 2024. [PMID: 38978395 DOI: 10.1021/acs.jctc.4c00452] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 07/10/2024]
Abstract
Because most proteins have buried active sites, protein tunnels or channels play a crucial role in the transport of small molecules into buried cavities for enzymatic catalysis. Tunnels can critically modulate the biological process of protein-ligand recognition. Various molecular dynamics methods have been developed for exploring and exploiting the protein-ligand conformational space to extract high-resolution details of the binding processes, a recent example being energetically unbiased high-throughput adaptive sampling simulations. The current study systematically contrasted the role of integrating prior knowledge while generating useful initial protein-ligand configurations, called seeds, for these simulations. Using a nontrivial system of a haloalkane dehalogenase mutant with multiple transport tunnels leading to a deeply buried active site, simulations were employed to derive kinetic models describing the process of association and dissociation of the substrate molecule. The most knowledge-based seed generation enabled high-throughput simulations that could more consistently capture the entire transport process, explore the complex network of transport tunnels, and predict equilibrium dissociation constants, koff/kon, on the same order of magnitude as experimental measurements. Overall, the infusion of more knowledge into the initial seeds of adaptive sampling simulations could render analyses of transport mechanisms in enzymes more consistent even for very complex biomolecular systems, thereby promoting drug development efforts and the rational design of enzymes with buried active sites.
Collapse
Affiliation(s)
- Dheeraj Kumar Sarkar
- Laboratory of Biomolecular Interactions and Transport, Department of Gene Expression, Institute of Molecular Biology and Biotechnology, Faculty of Biology, Adam Mickiewicz University, Uniwersytetu Poznanskiego 6, Poznan 61-614, Poland
- International Institute of Molecular and Cell Biology in Warsaw, Ks Trojdena 4, Warsaw 02-109, Poland
| | - Bartlomiej Surpeta
- Laboratory of Biomolecular Interactions and Transport, Department of Gene Expression, Institute of Molecular Biology and Biotechnology, Faculty of Biology, Adam Mickiewicz University, Uniwersytetu Poznanskiego 6, Poznan 61-614, Poland
- International Institute of Molecular and Cell Biology in Warsaw, Ks Trojdena 4, Warsaw 02-109, Poland
| | - Jan Brezovsky
- Laboratory of Biomolecular Interactions and Transport, Department of Gene Expression, Institute of Molecular Biology and Biotechnology, Faculty of Biology, Adam Mickiewicz University, Uniwersytetu Poznanskiego 6, Poznan 61-614, Poland
- International Institute of Molecular and Cell Biology in Warsaw, Ks Trojdena 4, Warsaw 02-109, Poland
| |
Collapse
|
3
|
Mehrani R, Mondal J, Ghazanfari D, Goetz DJ, McCall KD, Bergmeier SC, Sharma S. Capturing the Effects of Single Atom Substitutions on the Inhibition Efficiency of Glycogen Synthase Kinase-3β Inhibitors via Markov State Modeling and Experiments. J Chem Theory Comput 2024. [PMID: 38975986 DOI: 10.1021/acs.jctc.4c00311] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 07/09/2024]
Abstract
Small modifications in the chemical structure of ligands are known to dramatically change their ability to inhibit the activity of a protein. Unraveling the mechanisms that govern these dramatic changes requires scrutinizing the dynamics of protein-ligand binding and unbinding at the atomic level. As an exemplary case, we have studied Glycogen Synthase Kinase-3β (GSK-3β), a multifunctional kinase that has been implicated in a host of pathological processes. As such, there is a keen interest in identifying ligands that inhibit GSK-3β activity. One family of compounds that are highly selective and potent inhibitors of GSK-3β is exemplified by a molecule termed COB-187. COB-187 consists of a five-member heterocyclic ring with a thione at C2, a pyridine substituted methyl at N3, and a hydroxyl and phenyl at C4. We have studied the inhibition of GSK-3β by COB-187-related ligands that differ in a single heavy atom from each other (either in the location of nitrogen in their pyridine ring, or with the pyridine ring replaced by a phenyl ring), or in the length of the alkyl group joining the pyridine and the N3. The inhibition experiments show a large range of half-maximal inhibitory concentration (IC50) values from 10 nM to 10 μM, implying that these ligands exhibit vastly different propensities to inhibit GSK-3β. To explain these differences, we perform Markov State Modeling (MSM) using fully atomistic simulations. Our MSM results are in excellent agreement with the experiments in that they accurately capture differences in the binding propensities of the ligands. The simulations show that the binding propensities are related to the ligands' ability to attain a compact conformation where their two aromatic rings are spatially close. We rationalize this result by sampling numerous binding and unbinding events via funnel metadynamics simulations, which show that indeed while approaching the bound state, the ligands prefer to be in their compact conformation. We find that the presence of nitrogen in the aromatic ring increases the probability of attaining the compact conformation. Protein-ligand binding is understood to be dictated by the energetics of interactions and entropic factors, like the release of bound water from the binding pockets. This work shows that changes in the conformational distribution of ligands due to atom-level modifications in the structure play an important role in protein-ligand binding.
Collapse
Affiliation(s)
- Ramin Mehrani
- Department of Mechanical Engineering, Ohio University, Athens, Ohio 45701, United States
| | - Jagannath Mondal
- Center for Interdisciplinary Sciences, Tata Institute of Fundamental Research, Hyderabad 500046, India
| | - Davoud Ghazanfari
- Department of Chemical and Biomolecular Engineering, Ohio University, Athens, Ohio 45701, United States
| | - Douglas J Goetz
- Department of Chemical and Biomolecular Engineering, Ohio University, Athens, Ohio 45701, United States
- Biomedical Engineering Program, Ohio University, Athens, Ohio 45701, United States
| | - Kelly D McCall
- Biomedical Engineering Program, Ohio University, Athens, Ohio 45701, United States
- Department of Specialty Medicine, Ohio University, Athens, Ohio 45701, United States
- The Diabetes Institute, Ohio University, Athens, Ohio 45701, United States
- Molecular and Cellular Biology Program, Ohio University, Athens, Ohio 45701, United States
- Translational Biomedical Sciences Program, Ohio University, Athens, Ohio 45701, United States
| | - Stephen C Bergmeier
- Biomedical Engineering Program, Ohio University, Athens, Ohio 45701, United States
- Department of Chemistry and Biochemistry, Ohio University, Athens, Ohio 45701, United States
| | - Sumit Sharma
- Department of Chemical and Biomolecular Engineering, Ohio University, Athens, Ohio 45701, United States
| |
Collapse
|
4
|
Frasnetti E, Magni A, Castelli M, Serapian SA, Moroni E, Colombo G. Structures, dynamics, complexes, and functions: From classic computation to artificial intelligence. Curr Opin Struct Biol 2024; 87:102835. [PMID: 38744148 DOI: 10.1016/j.sbi.2024.102835] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/23/2024] [Revised: 04/14/2024] [Accepted: 04/22/2024] [Indexed: 05/16/2024]
Abstract
Computational approaches can provide highly detailed insight into the molecular recognition processes that underlie drug binding, the assembly of protein complexes, and the regulation of biological functional processes. Classical simulation methods can bridge a wide range of length- and time-scales typically involved in such processes. Lately, automated learning and artificial intelligence methods have shown the potential to expand the reach of physics-based approaches, ushering in the possibility to model and even design complex protein architectures. The synergy between atomistic simulations and AI methods is an emerging frontier with a huge potential for advances in structural biology. Herein, we explore various examples and frameworks for these approaches, providing select instances and applications that illustrate their impact on fundamental biomolecular problems.
Collapse
Affiliation(s)
- Elena Frasnetti
- Department of Chemistry, University of Pavia, via Taramelli 12, 27100 Pavia, Italy
| | - Andrea Magni
- Department of Chemistry, University of Pavia, via Taramelli 12, 27100 Pavia, Italy
| | - Matteo Castelli
- Department of Chemistry, University of Pavia, via Taramelli 12, 27100 Pavia, Italy
| | - Stefano A Serapian
- Department of Chemistry, University of Pavia, via Taramelli 12, 27100 Pavia, Italy
| | | | - Giorgio Colombo
- Department of Chemistry, University of Pavia, via Taramelli 12, 27100 Pavia, Italy.
| |
Collapse
|
5
|
Wang J, Miao Y. Ligand Gaussian accelerated Molecular Dynamics 3 (LiGaMD3): Improved Calculations of Binding Thermodynamics and Kinetics of Both Small Molecules and Flexible Peptides. BIORXIV : THE PREPRINT SERVER FOR BIOLOGY 2024:2024.05.06.592668. [PMID: 38766067 PMCID: PMC11100592 DOI: 10.1101/2024.05.06.592668] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/22/2024]
Abstract
Binding thermodynamics and kinetics play critical roles in drug design. However, it has proven challenging to efficiently predict ligand binding thermodynamics and kinetics of small molecules and flexible peptides using conventional Molecular Dynamics (cMD), due to limited simulation timescales. Based on our previously developed Ligand Gaussian accelerated Molecular Dynamics (LiGaMD) method, we present a new approach, termed "LiGaMD3", in which we introduce triple boosts into three individual energy terms that play important roles in small-molecule/peptide dissociation, rebinding and system conformational changes to improve the sampling efficiency of small-molecule/peptide interactions with target proteins. To validate the performance of LiGaMD3, MDM2 bound by a small molecule (Nutlin 3) and two highly flexible peptides (PMI and P53) were chosen as model systems. LiGaMD3 could efficiently capture repetitive small-molecule/peptide dissociation and binding events within 2 microsecond simulations. The predicted binding kinetic constant rates and free energies from LiGaMD3 agreed with available experimental values and previous simulation results. Therefore, LiGaMD3 provides a more general and efficient approach to capture dissociation and binding of both small-molecule ligand and flexible peptides, allowing for accurate prediction of their binding thermodynamics and kinetics.
Collapse
Affiliation(s)
- Jinan Wang
- Computational Medicine Program and Department of Pharmacology, University of North Carolina – Chapel Hill, Chapel Hill, North Carolina, USA 27599
| | - Yinglong Miao
- Computational Medicine Program and Department of Pharmacology, University of North Carolina – Chapel Hill, Chapel Hill, North Carolina, USA 27599
| |
Collapse
|
6
|
Schmitz B, Frieg B, Homeyer N, Jessen G, Gohlke H. Extracting binding energies and binding modes from biomolecular simulations of fragment binding to endothiapepsin. Arch Pharm (Weinheim) 2024; 357:e2300612. [PMID: 38319801 DOI: 10.1002/ardp.202300612] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/20/2023] [Revised: 12/18/2023] [Accepted: 01/10/2024] [Indexed: 02/08/2024]
Abstract
Fragment-based drug discovery (FBDD) aims to discover a set of small binding fragments that may be subsequently linked together. Therefore, in-depth knowledge of the individual fragments' structural and energetic binding properties is essential. In addition to experimental techniques, the direct simulation of fragment binding by molecular dynamics (MD) simulations became popular to characterize fragment binding. However, former studies showed that long simulation times and high computational demands per fragment are needed, which limits applicability in FBDD. Here, we performed short, unbiased MD simulations of direct fragment binding to endothiapepsin, a well-characterized model system of pepsin-like aspartic proteases. To evaluate the strengths and limitations of short MD simulations for the structural and energetic characterization of fragment binding, we predicted the fragments' absolute free energies and binding poses based on the direct simulations of fragment binding and compared the predictions to experimental data. The predicted absolute free energies are in fair agreement with the experiment. Combining the MD data with binding mode predictions from molecular docking approaches helped to correctly identify the most promising fragments for further chemical optimization. Importantly, all computations and predictions were done within 5 days, suggesting that MD simulations may become a viable tool in FBDD projects.
Collapse
Affiliation(s)
- Birte Schmitz
- Institute for Pharmaceutical and Medicinal Chemistry, Heinrich Heine University Düsseldorf, Düsseldorf, Germany
| | - Benedikt Frieg
- Institute for Pharmaceutical and Medicinal Chemistry, Heinrich Heine University Düsseldorf, Düsseldorf, Germany
- John von Neumann Institute for Computing (NIC), Jülich Supercomputing Centre (JSC), and Institute of Biological Information Processing (IBI-7: Structural Biochemistry), Forschungszentrum Jülich, Jülich, Germany
| | - Nadine Homeyer
- Institute for Pharmaceutical and Medicinal Chemistry, Heinrich Heine University Düsseldorf, Düsseldorf, Germany
| | - Gisela Jessen
- Institute for Pharmaceutical and Medicinal Chemistry, Heinrich Heine University Düsseldorf, Düsseldorf, Germany
| | - Holger Gohlke
- Institute for Pharmaceutical and Medicinal Chemistry, Heinrich Heine University Düsseldorf, Düsseldorf, Germany
- John von Neumann Institute for Computing (NIC), Jülich Supercomputing Centre (JSC), and Institute of Biological Information Processing (IBI-7: Structural Biochemistry), Forschungszentrum Jülich, Jülich, Germany
- Institute of Bio- and Geosciences (IBG-4: Bioinformatics), Forschungszentrum Jülich, Jülich, Germany
| |
Collapse
|
7
|
Xu C, Zhang X, Zhao L, Verkhivker GM, Bai F. Accurate Characterization of Binding Kinetics and Allosteric Mechanisms for the HSP90 Chaperone Inhibitors Using AI-Augmented Integrative Biophysical Studies. JACS AU 2024; 4:1632-1645. [PMID: 38665669 PMCID: PMC11040708 DOI: 10.1021/jacsau.4c00123] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 02/09/2024] [Revised: 03/15/2024] [Accepted: 03/18/2024] [Indexed: 04/28/2024]
Abstract
The binding kinetics of drugs to their targets are gradually being recognized as a crucial indicator of the efficacy of drugs in vivo, leading to the development of various computational methods for predicting the binding kinetics in recent years. However, compared with the prediction of binding affinity, the underlying structure and dynamic determinants of binding kinetics are more complicated. Efficient and accurate methods for predicting binding kinetics are still lacking. In this study, quantitative structure-kinetics relationship (QSKR) models were developed using 132 inhibitors targeting the ATP binding domain of heat shock protein 90α (HSP90α) to predict the dissociation rate constant (koff), enabling a direct assessment of the drug-target residence time. These models demonstrated good predictive performance, where hydrophobic and hydrogen bond interactions significantly influence the koff prediction. In subsequent applications, our models were used to assist in the discovery of new inhibitors for the N-terminal domain of HSP90α (N-HSP90α), demonstrating predictive capabilities on an experimental validation set with a new scaffold. In X-ray crystallography experiments, the loop-middle conformation of apo N-HSP90α was observed for the first time (previously, the loop-middle conformation had only been observed in holo-N-HSP90α structures). Interestingly, we observed different conformations of apo N-HSP90α simultaneously in an asymmetric unit, which was also observed in a holo-N-HSP90α structure, suggesting an equilibrium of conformations between different states in solution, which could be one of the determinants affecting the binding kinetics of the ligand. Different ligands can undergo conformational selection or alter the equilibrium of conformations, inducing conformational rearrangements and resulting in different effects on binding kinetics. We then used molecular dynamics simulations to describe conformational changes of apo N-HSP90α in different conformational states. In summary, the study of the binding kinetics and molecular mechanisms of N-HSP90α provides valuable information for the development of more targeted therapeutic approaches.
Collapse
Affiliation(s)
- Chao Xu
- Shanghai
Institute for Advanced Immunochemical Studies and School of Life Science
and Technology, ShanghaiTech University, 393 Middle Huaxia Road, Shanghai 201210, China
| | - Xianglei Zhang
- Shanghai
Institute for Advanced Immunochemical Studies and School of Life Science
and Technology, ShanghaiTech University, 393 Middle Huaxia Road, Shanghai 201210, China
| | - Lianghao Zhao
- Shanghai
Institute for Advanced Immunochemical Studies and School of Life Science
and Technology, ShanghaiTech University, 393 Middle Huaxia Road, Shanghai 201210, China
| | - Gennady M. Verkhivker
- Keck
Center for Science and Engineering, Graduate Program in Computational
and Data Sciences, Schmid College of Science and Technology, Chapman University, Orange, California 92866, United States
- Department
of Biomedical and Pharmaceutical Sciences, Chapman University School of Pharmacy, Irvine, California 92618, United States
| | - Fang Bai
- Shanghai
Institute for Advanced Immunochemical Studies and School of Life Science
and Technology, ShanghaiTech University, 393 Middle Huaxia Road, Shanghai 201210, China
- School
of Information Science and Technology, ShanghaiTech
University, 393 Middle Huaxia Road, Shanghai 201210, China
- Shanghai
Clinical Research and Trial Center, Shanghai 201210, China
| |
Collapse
|
8
|
Daoud S, Taha M. Protein characteristics substantially influence the propensity of activity cliffs among kinase inhibitors. Sci Rep 2024; 14:9058. [PMID: 38643174 PMCID: PMC11032345 DOI: 10.1038/s41598-024-59501-w] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/09/2023] [Accepted: 04/11/2024] [Indexed: 04/22/2024] Open
Abstract
Activity cliffs (ACs) are pairs of structurally similar molecules with significantly different affinities for a biotarget, posing a challenge in computer-assisted drug discovery. This study focuses on protein kinases, significant therapeutic targets, with some exhibiting ACs while others do not despite numerous inhibitors. The hypothesis that the presence of ACs is dependent on the target protein and its complete structural context is explored. Machine learning models were developed to link protein properties to ACs, revealing specific tripeptide sequences and overall protein properties as critical factors in ACs occurrence. The study highlights the importance of considering the entire protein matrix rather than just the binding site in understanding ACs. This research provides valuable insights for drug discovery and design, paving the way for addressing ACs-related challenges in modern computational approaches.
Collapse
Affiliation(s)
- Safa Daoud
- Department of Pharmaceutical Chemistry and Pharmacognosy, Faculty of Pharmacy, Applied Sciences Private University, Amman, Jordan.
| | - Mutasem Taha
- Department of Pharmaceutical Sciences, Faculty of Pharmacy, University of Jordan, Amman, Jordan.
| |
Collapse
|
9
|
Yuan Y, Mao X, Pan X, Zhang R, Su W. Kinetic Ensemble of Tau Protein through the Markov State Model and Deep Learning Analysis. J Chem Theory Comput 2024; 20:2947-2958. [PMID: 38501645 DOI: 10.1021/acs.jctc.3c01211] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/20/2024]
Abstract
The ordered assembly of Tau protein into filaments characterizes Alzheimer's and other neurodegenerative diseases, and thus, stabilization of Tau protein is a promising avenue for tauopathies therapy. To dissect the underlying aggregation mechanisms on Tau, we employ a set of molecular simulations and the Markov state model to determine the kinetics of ensemble of K18. K18 is the microtubule-binding domain of Tau protein and plays a vital role in the microtubule assembly, recycling processes, and amyloid fibril formation. Here, we efficiently explore the conformation of K18 with about 150 μs lifetimes in silico. Our results observe that all four repeat regions (R1-R4) are very dynamic, featuring frequent conformational conversion and lacking stable conformations, and the R2 region is more flexible than the R1, R3, and R4 regions. Additionally, it is worth noting that residues 300-310 in R2-R3 and residues 319-336 in R3 tend to form sheet structures, indicating that K18 has a broader functional role than individual repeat monomers. Finally, the simulations combined with Markov state models and deep learning reveal 5 key conformational states along the transition pathway and provide the information on the microsecond time scale interstate transition rates. Overall, this study offers significant insights into the molecular mechanism of Tau pathological aggregation and develops novel strategies for both securing tauopathies and advancing drug discovery.
Collapse
Affiliation(s)
- Yongna Yuan
- School of Information Science & Engineering, Lanzhou University, South Tianshui Road, Lanzhou 730000, Gansu, China
| | - Xuqi Mao
- School of Information Science & Engineering, Lanzhou University, South Tianshui Road, Lanzhou 730000, Gansu, China
| | - Xiaohang Pan
- School of Information Science & Engineering, Lanzhou University, South Tianshui Road, Lanzhou 730000, Gansu, China
| | - Ruisheng Zhang
- School of Information Science & Engineering, Lanzhou University, South Tianshui Road, Lanzhou 730000, Gansu, China
| | - Wei Su
- School of Information Science & Engineering, Lanzhou University, South Tianshui Road, Lanzhou 730000, Gansu, China
| |
Collapse
|
10
|
Wu Y, Cao S, Qiu Y, Huang X. Tutorial on how to build non-Markovian dynamic models from molecular dynamics simulations for studying protein conformational changes. J Chem Phys 2024; 160:121501. [PMID: 38516972 DOI: 10.1063/5.0189429] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/28/2023] [Accepted: 02/20/2024] [Indexed: 03/23/2024] Open
Abstract
Protein conformational changes play crucial roles in their biological functions. In recent years, the Markov State Model (MSM) constructed from extensive Molecular Dynamics (MD) simulations has emerged as a powerful tool for modeling complex protein conformational changes. In MSMs, dynamics are modeled as a sequence of Markovian transitions among metastable conformational states at discrete time intervals (called lag time). A major challenge for MSMs is that the lag time must be long enough to allow transitions among states to become memoryless (or Markovian). However, this lag time is constrained by the length of individual MD simulations available to track these transitions. To address this challenge, we have recently developed Generalized Master Equation (GME)-based approaches, encoding non-Markovian dynamics using a time-dependent memory kernel. In this Tutorial, we introduce the theory behind two recently developed GME-based non-Markovian dynamic models: the quasi-Markov State Model (qMSM) and the Integrative Generalized Master Equation (IGME). We subsequently outline the procedures for constructing these models and provide a step-by-step tutorial on applying qMSM and IGME to study two peptide systems: alanine dipeptide and villin headpiece. This Tutorial is available at https://github.com/xuhuihuang/GME_tutorials. The protocols detailed in this Tutorial aim to be accessible for non-experts interested in studying the biomolecular dynamics using these non-Markovian dynamic models.
Collapse
Affiliation(s)
- Yue Wu
- Department of Chemistry, Theoretical Chemistry Institute, University of Wisconsin-Madison, Madison, Wisconsin 53706, USA
| | - Siqin Cao
- Department of Chemistry, Theoretical Chemistry Institute, University of Wisconsin-Madison, Madison, Wisconsin 53706, USA
| | - Yunrui Qiu
- Department of Chemistry, Theoretical Chemistry Institute, University of Wisconsin-Madison, Madison, Wisconsin 53706, USA
| | - Xuhui Huang
- Department of Chemistry, Theoretical Chemistry Institute, University of Wisconsin-Madison, Madison, Wisconsin 53706, USA
- Data Science Institute, University of Wisconsin-Madison, Madison, Wisconsin 53706, USA
| |
Collapse
|
11
|
Smith L, Novak B, Osato M, Mobley DL, Bowman GR. PopShift: A Thermodynamically Sound Approach to Estimate Binding Free Energies by Accounting for Ligand-Induced Population Shifts from a Ligand-Free Markov State Model. J Chem Theory Comput 2024; 20:1036-1050. [PMID: 38291966 PMCID: PMC10867841 DOI: 10.1021/acs.jctc.3c00870] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/08/2023] [Revised: 11/28/2023] [Accepted: 11/29/2023] [Indexed: 02/01/2024]
Abstract
Obtaining accurate binding free energies from in silico screens has been a long-standing goal for the computational chemistry community. However, accuracy and computational cost are at odds with one another, limiting the utility of methods that perform this type of calculation. Many methods achieve massive scale by explicitly or implicitly assuming that the target protein adopts a single structure, or undergoes limited fluctuations around that structure, to minimize computational cost. Others simulate each protein-ligand complex of interest, accepting lower throughput in exchange for better predictions of binding affinities. Here, we present the PopShift framework for accounting for the ensemble of structures a protein adopts and their relative probabilities. Protein degrees of freedom are enumerated once, and then arbitrarily many molecules can be screened against this ensemble. Specifically, we use Markov state models (MSMs) as a compressed representation of a protein's thermodynamic ensemble. We start with a ligand-free MSM and then calculate how addition of a ligand shifts the populations of each protein conformational state based on the strength of the interaction between that protein conformation and the ligand. In this work we use docking to estimate the affinity between a given protein structure and ligand, but any estimator of binding affinities could be used in the PopShift framework. We test PopShift on the classic benchmark pocket T4 Lysozyme L99A. We find that PopShift is more accurate than common strategies, such as docking to a single structure and traditional ensemble docking─producing results that compare favorably with alchemical binding free energy calculations in terms of RMSE but not correlation─and may have a more favorable computational cost profile in some applications. In addition to predicting binding free energies and ligand poses, PopShift also provides insight into how the probability of different protein structures is shifted upon addition of various concentrations of ligand, providing a platform for predicting affinities and allosteric effects of ligand binding. Therefore, we expect PopShift will be valuable for hit finding and for providing insight into phenomena like allostery.
Collapse
Affiliation(s)
- Louis
G. Smith
- Departments
of Biochemistry & Biophysics and Bioengineering, University of Pennsylvania, Philadelphia, Pennsylvania 19104, United States
| | - Borna Novak
- Department
of Biochemistry and Molecular Biophysics, Washington University in St. Louis, St. Louis, Missouri 63130, United States
- Medical
Scientist Training Program, Washington University
in St. Louis, St. Louis, Missouri 63130, United
States
| | - Meghan Osato
- School
of Pharmacy and Pharmaceutical Sciences, University of California, Irvine, Irvine, California 92697, United States
| | - David L. Mobley
- School
of Pharmacy and Pharmaceutical Sciences, University of California, Irvine, Irvine, California 92697, United States
| | - Gregory R. Bowman
- Departments
of Biochemistry & Biophysics and Bioengineering, University of Pennsylvania, Philadelphia, Pennsylvania 19104, United States
| |
Collapse
|
12
|
Sisk TR, Robustelli P. Folding-upon-binding pathways of an intrinsically disordered protein from a deep Markov state model. Proc Natl Acad Sci U S A 2024; 121:e2313360121. [PMID: 38294935 PMCID: PMC10861926 DOI: 10.1073/pnas.2313360121] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/10/2023] [Accepted: 11/22/2023] [Indexed: 02/02/2024] Open
Abstract
A central challenge in the study of intrinsically disordered proteins is the characterization of the mechanisms by which they bind their physiological interaction partners. Here, we utilize a deep learning-based Markov state modeling approach to characterize the folding-upon-binding pathways observed in a long timescale molecular dynamics simulation of a disordered region of the measles virus nucleoprotein NTAIL reversibly binding the X domain of the measles virus phosphoprotein complex. We find that folding-upon-binding predominantly occurs via two distinct encounter complexes that are differentiated by the binding orientation, helical content, and conformational heterogeneity of NTAIL. We observe that folding-upon-binding predominantly proceeds through a multi-step induced fit mechanism with several intermediates and do not find evidence for the existence of canonical conformational selection pathways. We observe four kinetically separated native-like bound states that interconvert on timescales of eighty to five hundred nanoseconds. These bound states share a core set of native intermolecular contacts and stable NTAIL helices and are differentiated by a sequential formation of native and non-native contacts and additional helical turns. Our analyses provide an atomic resolution structural description of intermediate states in a folding-upon-binding pathway and elucidate the nature of the kinetic barriers between metastable states in a dynamic and heterogenous, or "fuzzy", protein complex.
Collapse
Affiliation(s)
- Thomas R. Sisk
- Department of Chemistry, Dartmouth College, Hanover, NH03755
| | - Paul Robustelli
- Department of Chemistry, Dartmouth College, Hanover, NH03755
| |
Collapse
|
13
|
Madsen JJ, Yu W. Dynamic Nature of Staphylococcus aureus Type I Signal Peptidases. BIORXIV : THE PREPRINT SERVER FOR BIOLOGY 2024:2024.01.23.576923. [PMID: 38328037 PMCID: PMC10849702 DOI: 10.1101/2024.01.23.576923] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 02/09/2024]
Abstract
Molecular dynamics simulations are used to interrogate the dynamic nature of Staphylococcus aureus Type I signal peptidases, SpsA and SpsB, including the impact of the P29S mutation of SpsB. Fluctuations and plasticity- rigidity characteristics vary among the proteins, particularly in the extracellular domain. Intriguingly, the P29S mutation, which influences susceptibility to arylomycin antibiotics, affect the mechanically coupled motions in SpsB. The integrity of the active site is crucial for catalytic competency, and variations in sampled structural conformations among the proteins are consistent with diverse peptidase capabilities. We also explored the intricate interactions between the proteins and the model S. aureus membrane. It was observed that certain membrane-inserted residues in the loop around residue 50 (50s) and C-terminal loops, beyond the transmembrane domain, give rise to direct interactions with lipids in the bilayer membrane. Our findings are discussed in the context of functional knowledge about these signal peptidases, offering additional understanding of dynamic aspects relevant to some cellular processes with potential implications for drug targeting strategies.
Collapse
Affiliation(s)
- Jesper J. Madsen
- Department of Molecular Medicine, Morsani College of Medicine, University of South Florida, Tampa, Florida 33612, United States of America
- Center for Global Health and Infectious Diseases Research, Global and Planetary Health, College of Public Health, University of South Florida, Tampa, Florida 33612, United States of America
| | - Wenqi Yu
- Department of Molecular Biosciences, College of Arts and Sciences, University of South Florida, Tampa, Florida 33612, United States of America
| |
Collapse
|
14
|
Singh K, Reddy G. Excited States of apo-Guanidine-III Riboswitch Contribute to Guanidinium Binding through Both Conformational and Induced-Fit Mechanisms. J Chem Theory Comput 2024; 20:421-435. [PMID: 38134376 DOI: 10.1021/acs.jctc.3c00999] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/24/2023]
Abstract
Riboswitches are mRNA segments that regulate gene expression through conformational changes driven by their cognate ligand binding. The ykkC motif forms a riboswitch class that selectively senses a guanidinium ion (Gdm+) and regulates the downstream expression of proteins which aid in the efflux of excess Gdm+ from the cells. The aptamer domain (AD) of the guanidine-III riboswitch forms an H-type pseudoknot with a triple helical domain that binds a Gdm+. We studied the binding of Gdm+ to the AD of the guanidine (ykkC)-III riboswitch using computer simulations to probe the specificity of the riboswitch to Gdm+ binding. We show that Gdm+ binding is a fast process occurring on the nanosecond time scale, with minimal conformational changes to the AD. Using machine learning and Markov-state models, we identified the excited conformational states of the AD, which have a high Gdm+ binding propensity, making the Gdm+ binding landscape complex exhibiting both conformational selection and induced-fit mechanisms. The proposed apo-AD excited states and their role in the ligand-sensing mechanism are amenable to experimental verification. Further, targeting these excited-state conformations in discovering new antibiotics can be explored.
Collapse
Affiliation(s)
- Kushal Singh
- Solid State and Structural Chemistry Unit, Indian Institute of Science, Bengaluru 560012 Karnataka, India
| | - Govardhan Reddy
- Solid State and Structural Chemistry Unit, Indian Institute of Science, Bengaluru 560012 Karnataka, India
| |
Collapse
|
15
|
Paul SK, Guendouzi A, Banerjee A, Guendouzi A, Haldar R. Identification of approved drugs with ALDH1A1 inhibitory potential aimed at enhancing chemotherapy sensitivity in cancer cells: an in-silico drug repurposing approach. J Biomol Struct Dyn 2024:1-15. [PMID: 38189344 DOI: 10.1080/07391102.2023.2300127] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/01/2023] [Accepted: 12/21/2023] [Indexed: 01/09/2024]
Abstract
The aldehyde dehydrogenase 1A1 (ALDH1A1) also known as retinal dehydrogenase, is an enzyme normally involved in the cellular metabolism, development and detoxification processes in healthy cells. However, it's also considered a cancer stem cell marker and its high levels of expression in several cancers, including breast, lung, ovarian, and colon cancer have been associated with poor prognosis and resistance to chemotherapy. Given its crucial role in chemotherapy resistance by detoxification of chemotherapeutic drugs, ALDH1A1 has attracted significant research interest as a potential therapeutic target for cancer. Though a few synthetic inhibitors of ALDH1A1 have been synthesized and their efficacy has been proved in-vitro and in-vivo studies, none of them have passed clinical trials so far. In this scenario, we have performed an in-silico study to verify whether any of the already approved drugs used for various purposes has the ability to inhibit catalytic activity of ALDH1A1, so that they can be repurposed for cancer therapy. Keeping in mind the feasibility of repurposing in a larger population we have selected the approved drugs from five widely used drug categories such as antibiotic, antiviral, antifungal, anti diabetic and antihypertensive for screening. Computational techniques like molecular docking, molecular dynamics simulations and MM-PBSA binding energy calculation have been used in this study to screen the approved drugs. Based on the logical analysis of results, we propose that three drugs - telmisartan, irbesartan and maraviroc can inhibit the catalytic activity of ALDH1A1 and thus can be repurposed to increase chemotherapy sensitivity in cancer cells.Communicated by Ramaswamy H. Sarma.
Collapse
Affiliation(s)
- Sanjay Kumar Paul
- Department of Physiology, University of Calcutta, Kolkata, India
- Department of Zoology, Rammohan College, Kolkata, West Bengal, India
| | - Abdelmadjid Guendouzi
- Center for Research in Pharmaceutical Sciences (CRSP), Constantine, Algeria
- Ecole Normale Supérieure ENS Constantine, Constantine, Algeria
| | | | | | - Rajen Haldar
- Department of Physiology, University of Calcutta, Kolkata, India
| |
Collapse
|
16
|
Tian J, Dong X, Wu T, Wen P, Liu X, Zhang M, An X, Shi D. Revealing the conformational dynamics of UDP-GlcNAc recognition by O-GlcNAc transferase via Markov state model. Int J Biol Macromol 2024; 256:128405. [PMID: 38016609 DOI: 10.1016/j.ijbiomac.2023.128405] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/12/2023] [Revised: 11/20/2023] [Accepted: 11/22/2023] [Indexed: 11/30/2023]
Abstract
The O-linked N-acetylglucosamine (O-GlcNAc) glycosylation is a critical post-translational modification and closely linked to various physiological and pathological conditions. The O-GlcNAc transferase (OGT) functions as the only glycosyltransferase of O-GlcNAc glycosylation by transferring GlcNAc from UDP-GlcNAc to serine or threonine residues on protein substrates. The interaction mode of UDP-GlcNAc against OGT has been preliminarily revealed by the crystal structures, yet an atomic-level comprehension for the conformational dynamics of the recognition process remains elusive. Here, we construct the Markov state model based on extensive all-atom molecular dynamics (MD) simulations with an aggregated simulation time of ∼9 μs, and reveal that the UDP-GlcNAc recognition process by OGT encompasses four key metastable states, occurring within an estimated timescale of ∼10 μs. During UDP-GlcNAc recognition process, we find the pyrophosphate moiety (P2O52-) initially anchors to the active pocket via salt bridge and hydrogen bonds, facilitating subsequent binding of the uridine and GlcNAc moieties. Furthermore, the functional roles of K842 involved in the salt bridge with P2O52- were evaluated through extra mutant MD simulations. Overall, our study provides valuable insights into the UDP-GlcNAc recognition mechanism by OGT, which could further aid in mechanistic studies of O-GlcNAc glycosylation and drug development targeting on OGT.
Collapse
Affiliation(s)
- Jiaqi Tian
- School of Medical Informatics and Engineering, Xuzhou Medical University, Xuzhou, Jiangsu Province, China
| | - Xin Dong
- School of Medical Informatics and Engineering, Xuzhou Medical University, Xuzhou, Jiangsu Province, China
| | - Tianshuo Wu
- School of Medical Informatics and Engineering, Xuzhou Medical University, Xuzhou, Jiangsu Province, China
| | - Pengbo Wen
- School of Medical Informatics and Engineering, Xuzhou Medical University, Xuzhou, Jiangsu Province, China
| | - Xin Liu
- School of Medical Informatics and Engineering, Xuzhou Medical University, Xuzhou, Jiangsu Province, China
| | - Mengying Zhang
- School of Medical Informatics and Engineering, Xuzhou Medical University, Xuzhou, Jiangsu Province, China
| | - Xiaoli An
- School of Chemical Engineering, Institute of Pharmaceutical Engineering Technology and Application, Sichuan University of Science & Engineering, Xueyuan Street 180, Huixing Road, Zigong 643000, Sichuan, China.
| | - Danfeng Shi
- Warshel Institute for Computational Biology, School of Life and Health Sciences, School of Medicine, The Chinese University of Hong Kong, Shenzhen 518172, Guangdong, China.
| |
Collapse
|
17
|
Raddi RM, Voelz VA. Markov State Model of Solvent Features Reveals Water Dynamics in Protein-Peptide Binding. J Phys Chem B 2023; 127:10682-10690. [PMID: 38078851 DOI: 10.1021/acs.jpcb.3c04775] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/22/2023]
Abstract
In this work, we investigate the role of solvent in the binding reaction of the p53 transactivation domain (TAD) peptide to its receptor MDM2. Previously, our group generated 831 μs of explicit-solvent aggregate molecular simulation trajectory data for the MDM2-p53 peptide binding reaction using large-scale distributed computing and subsequently built a Markov State Model (MSM) of the binding reaction (Zhou et al. 2017). Here, we perform a tICA analysis and construct an MSM with similar hyperparameters while using only solvent-based structural features. We find a remarkably similar landscape but accelerated implied timescales for the slowest motions. The solvent shells contributing most to the first tICA eigenvector are those centered on Lys24 and Thr18 of the p53 TAD peptide in the range of 3-6 Å. Important solvent shells were visualized to reveal solvation and desolvation transitions along the peptide-protein binding trajectories. Our results provide a solvent-centric view of the hydrophobic effect in action for a realistic peptide-protein binding scenario.
Collapse
Affiliation(s)
- Robert M Raddi
- Department of Chemistry, Temple University, Philadelphia, Pennsylvania 19122, United States
| | - Vincent A Voelz
- Department of Chemistry, Temple University, Philadelphia, Pennsylvania 19122, United States
| |
Collapse
|
18
|
Patil K, Wang Y, Chen Z, Suresh K, Radhakrishnan R. Activating mutations drive human MEK1 kinase using a gear-shifting mechanism. Biochem J 2023; 480:1733-1751. [PMID: 37869794 PMCID: PMC10872882 DOI: 10.1042/bcj20230281] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/11/2023] [Revised: 09/30/2023] [Accepted: 10/20/2023] [Indexed: 10/24/2023]
Abstract
There is an unmet need to classify cancer-promoting kinase mutations in a mechanistically cognizant way. The challenge is to understand how mutations stabilize different kinase configurations to alter function, and how this influences pathogenic potential of the kinase and its responses to therapeutic inhibitors. This goal is made more challenging by the complexity of the mutational landscape of diseases, and is further compounded by the conformational plasticity of each variant where multiple conformations coexist. We focus here on the human MEK1 kinase, a vital component of the RAS/MAPK pathway in which mutations cause cancers and developmental disorders called RASopathies. We sought to explore how these mutations alter the human MEK1 kinase at atomic resolution by utilizing enhanced sampling simulations and free energy calculations. We computationally mapped the different conformational stabilities of individual mutated systems by delineating the free energy landscapes, and showed how this relates directly to experimentally quantified developmental transformation potentials of the mutations. We conclude that mutations leverage variations in the hydrogen bonding network associated with the conformational plasticity to progressively stabilize the active-like conformational state of the kinase while destabilizing the inactive-like state. The mutations alter residue-level internal molecular correlations by differentially prioritizing different conformational states, delineating the various modes of MEK1 activation reminiscent of a gear-shifting mechanism. We define the molecular basis of conversion of this kinase from its inactive to its active state, connecting structure, dynamics, and function by delineating the energy landscape and conformational plasticity, thus augmenting our understanding of MEK1 regulation.
Collapse
Affiliation(s)
- Keshav Patil
- Department of Chemical and Biomolecular Engineering, University of Pennsylvania, Philadelphia, PA, U.S.A
| | - Yiming Wang
- Department of Chemical and Biomolecular Engineering, University of Pennsylvania, Philadelphia, PA, U.S.A
| | - Zhangtao Chen
- Department of Bioengineering, University of Pennsylvania, Philadelphia, PA, U.S.A
| | - Krishna Suresh
- Department of Bioengineering, University of Pennsylvania, Philadelphia, PA, U.S.A
| | - Ravi Radhakrishnan
- Department of Chemical and Biomolecular Engineering, University of Pennsylvania, Philadelphia, PA, U.S.A
- Department of Bioengineering, University of Pennsylvania, Philadelphia, PA, U.S.A
| |
Collapse
|
19
|
Pasarkar AP, Bencomo GM, Olsson S, Dieng AB. Vendi sampling for molecular simulations: Diversity as a force for faster convergence and better exploration. J Chem Phys 2023; 159:144108. [PMID: 37823459 DOI: 10.1063/5.0166172] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/03/2023] [Accepted: 09/25/2023] [Indexed: 10/13/2023] Open
Abstract
Molecular dynamics (MD) is the method of choice for understanding the structure, function, and interactions of molecules. However, MD simulations are limited by the strong metastability of many molecules, which traps them in a single conformation basin for an extended amount of time. Enhanced sampling techniques, such as metadynamics and replica exchange, have been developed to overcome this limitation and accelerate the exploration of complex free energy landscapes. In this paper, we propose Vendi Sampling, a replica-based algorithm for increasing the efficiency and efficacy of the exploration of molecular conformation spaces. In Vendi sampling, replicas are simulated in parallel and coupled via a global statistical measure, the Vendi Score, to enhance diversity. Vendi sampling allows for the recovery of unbiased sampling statistics and dramatically improves sampling efficiency. We demonstrate the effectiveness of Vendi sampling in improving molecular dynamics simulations by showing significant improvements in coverage and mixing between metastable states and convergence of free energy estimates for four common benchmarks, including Alanine Dipeptide and Chignolin.
Collapse
Affiliation(s)
- Amey P Pasarkar
- Vertaix, Department of Computer Science, Princeton University, 35 Olden Street, Princeton, New Jersey 08544, USA
| | - Gianluca M Bencomo
- Department of Computer Science, Princeton University, 35 Olden Street, Princeton, New Jersey 08544, USA
| | - Simon Olsson
- Department of Computer Science and Engineering, Chalmers University of Technology, Rännvägen 6, 41258 Gothenburg, Sweden
| | - Adji Bousso Dieng
- Vertaix, Department of Computer Science, Princeton University, 35 Olden Street, Princeton, New Jersey 08544, USA
| |
Collapse
|
20
|
Cao S, Qiu Y, Kalin ML, Huang X. Integrative generalized master equation: A method to study long-timescale biomolecular dynamics via the integrals of memory kernels. J Chem Phys 2023; 159:134106. [PMID: 37787134 PMCID: PMC11005468 DOI: 10.1063/5.0167287] [Citation(s) in RCA: 5] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/11/2023] [Accepted: 09/18/2023] [Indexed: 10/04/2023] Open
Abstract
The generalized master equation (GME) provides a powerful approach to study biomolecular dynamics via non-Markovian dynamic models built from molecular dynamics (MD) simulations. Previously, we have implemented the GME, namely the quasi Markov State Model (qMSM), where we explicitly calculate the memory kernel and propagate dynamics using a discretized GME. qMSM can be constructed with much shorter MD trajectories than the MSM. However, since qMSM needs to explicitly compute the time-dependent memory kernels, it is heavily affected by the numerical fluctuations of simulation data when applied to study biomolecular conformational changes. This can lead to numerical instability of predicted long-time dynamics, greatly limiting the applicability of qMSM in complicated biomolecules. We present a new method, the Integrative GME (IGME), in which we analytically solve the GME under the condition when the memory kernels have decayed to zero. Our IGME overcomes the challenges of the qMSM by using the time integrations of memory kernels, thereby avoiding the numerical instability caused by explicit computation of time-dependent memory kernels. Using our solutions of the GME, we have developed a new approach to compute long-time dynamics based on MD simulations in a numerically stable, accurate and efficient way. To demonstrate its effectiveness, we have applied the IGME in three biomolecules: the alanine dipeptide, FIP35 WW-domain, and Taq RNA polymerase. In each system, the IGME achieves significantly smaller fluctuations for both memory kernels and long-time dynamics compared to the qMSM. We anticipate that the IGME can be widely applied to investigate biomolecular conformational changes.
Collapse
Affiliation(s)
- Siqin Cao
- Department of Chemistry, Theoretical Chemistry Institute, University of Wisconsin-Madison, Madison, Wisconsin 53706, USA
| | - Yunrui Qiu
- Department of Chemistry, Theoretical Chemistry Institute, University of Wisconsin-Madison, Madison, Wisconsin 53706, USA
| | - Michael L. Kalin
- Biophysics Graduate Program, University of Wisconsin-Madison, Madison, Wisconsin 53706, USA
| | - Xuhui Huang
- Department of Chemistry, Theoretical Chemistry Institute, University of Wisconsin-Madison, Madison, Wisconsin 53706, USA
| |
Collapse
|
21
|
Deeks HM, Zinovjev K, Barnoud J, Mulholland AJ, van der Kamp MW, Glowacki DR. Free energy along drug-protein binding pathways interactively sampled in virtual reality. Sci Rep 2023; 13:16665. [PMID: 37794083 PMCID: PMC10551034 DOI: 10.1038/s41598-023-43523-x] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/30/2023] [Accepted: 09/25/2023] [Indexed: 10/06/2023] Open
Abstract
We describe a two-step approach for combining interactive molecular dynamics in virtual reality (iMD-VR) with free energy (FE) calculation to explore the dynamics of biological processes at the molecular level. We refer to this combined approach as iMD-VR-FE. Stage one involves using a state-of-the-art 'human-in-the-loop' iMD-VR framework to generate a diverse range of protein-ligand unbinding pathways, benefitting from the sophistication of human spatial and chemical intuition. Stage two involves using the iMD-VR-sampled pathways as initial guesses for defining a path-based reaction coordinate from which we can obtain a corresponding free energy profile using FE methods. To investigate the performance of the method, we apply iMD-VR-FE to investigate the unbinding of a benzamidine ligand from a trypsin protein. The binding free energy calculated using iMD-VR-FE is similar for each pathway, indicating internal consistency. Moreover, the resulting free energy profiles can distinguish energetic differences between pathways corresponding to various protein-ligand conformations (e.g., helping to identify pathways that are more favourable) and enable identification of metastable states along the pathways. The two-step iMD-VR-FE approach offers an intuitive way for researchers to test hypotheses for candidate pathways in biomolecular systems, quickly obtaining both qualitative and quantitative insight.
Collapse
Affiliation(s)
- Helen M Deeks
- Center for Computational Chemistry, School of Chemistry, University of Bristol, Bristol, BS8 1TS, UK
| | - Kirill Zinovjev
- Departamento de Química Física, Universidad de Valencia, 46100, Burjassot, Spain
- School of Biochemistry, University of Bristol, Bristol, BS8 1TD, UK
| | - Jonathan Barnoud
- Center for Computational Chemistry, School of Chemistry, University of Bristol, Bristol, BS8 1TS, UK
- CiTIUS | Centro Singular de Investigación en Tecnoloxías Intelixentes da USC, Rúa de Jenaro de la Fuente, s/n, 15705, Santiago de Compostela, A Coruña, Spain
| | - Adrian J Mulholland
- Center for Computational Chemistry, School of Chemistry, University of Bristol, Bristol, BS8 1TS, UK
| | - Marc W van der Kamp
- Center for Computational Chemistry, School of Chemistry, University of Bristol, Bristol, BS8 1TS, UK.
- School of Biochemistry, University of Bristol, Bristol, BS8 1TD, UK.
| | - David R Glowacki
- CiTIUS | Centro Singular de Investigación en Tecnoloxías Intelixentes da USC, Rúa de Jenaro de la Fuente, s/n, 15705, Santiago de Compostela, A Coruña, Spain.
| |
Collapse
|
22
|
Conflitti P, Raniolo S, Limongelli V. Perspectives on Ligand/Protein Binding Kinetics Simulations: Force Fields, Machine Learning, Sampling, and User-Friendliness. J Chem Theory Comput 2023; 19:6047-6061. [PMID: 37656199 PMCID: PMC10536999 DOI: 10.1021/acs.jctc.3c00641] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/14/2023] [Indexed: 09/02/2023]
Abstract
Computational techniques applied to drug discovery have gained considerable popularity for their ability to filter potentially active drugs from inactive ones, reducing the time scale and costs of preclinical investigations. The main focus of these studies has historically been the search for compounds endowed with high affinity for a specific molecular target to ensure the formation of stable and long-lasting complexes. Recent evidence has also correlated the in vivo drug efficacy with its binding kinetics, thus opening new fascinating scenarios for ligand/protein binding kinetic simulations in drug discovery. The present article examines the state of the art in the field, providing a brief summary of the most popular and advanced ligand/protein binding kinetics techniques and evaluating their current limitations and the potential solutions to reach more accurate kinetic models. Particular emphasis is put on the need for a paradigm change in the present methodologies toward ligand and protein parametrization, the force field problem, characterization of the transition states, the sampling issue, and algorithms' performance, user-friendliness, and data openness.
Collapse
Affiliation(s)
- Paolo Conflitti
- Faculty
of Biomedical Sciences, Euler Institute, Universitá della Svizzera italiana (USI), 6900 Lugano, Switzerland
| | - Stefano Raniolo
- Faculty
of Biomedical Sciences, Euler Institute, Universitá della Svizzera italiana (USI), 6900 Lugano, Switzerland
| | - Vittorio Limongelli
- Faculty
of Biomedical Sciences, Euler Institute, Universitá della Svizzera italiana (USI), 6900 Lugano, Switzerland
- Department
of Pharmacy, University of Naples “Federico
II”, 80131 Naples, Italy
| |
Collapse
|
23
|
Rizzi V, Aureli S, Ansari N, Gervasio FL. OneOPES, a Combined Enhanced Sampling Method to Rule Them All. J Chem Theory Comput 2023; 19:5731-5742. [PMID: 37603295 PMCID: PMC10500989 DOI: 10.1021/acs.jctc.3c00254] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/06/2023] [Indexed: 08/22/2023]
Abstract
Enhanced sampling techniques have revolutionized molecular dynamics (MD) simulations, enabling the study of rare events and the calculation of free energy differences in complex systems. One of the main families of enhanced sampling techniques uses physical degrees of freedom called collective variables (CVs) to accelerate a system's dynamics and recover the original system's statistics. However, encoding all the relevant degrees of freedom in a limited number of CVs is challenging, particularly in large biophysical systems. Another category of techniques, such as parallel tempering, simulates multiple replicas of the system in parallel, without requiring CVs. However, these methods may explore less relevant high-energy portions of the phase space and become computationally expensive for large systems. To overcome the limitations of both approaches, we propose a replica exchange method called OneOPES that combines the power of multireplica simulations and CV-based enhanced sampling. This method efficiently accelerates the phase space sampling without the need for ideal CVs, extensive parameters fine tuning nor the use of a large number of replicas, as demonstrated by its successful applications to protein-ligand binding and protein folding benchmark systems. Our approach shows promise as a new direction in the development of enhanced sampling techniques for molecular dynamics simulations, providing an efficient and robust framework for the study of complex and unexplored problems.
Collapse
Affiliation(s)
- Valerio Rizzi
- School
of Pharmaceutical Sciences, University of
Geneva, Rue Michel Servet 1, 1206 Genève, Switzerland
- Institute
of Pharmaceutical Sciences of Western Switzerland (ISPSO), University of Geneva, 1206 Genève, Switzerland
- Swiss
Institute of Bioinformatics, University
of Geneva, 1206 Genève, Switzerland
| | - Simone Aureli
- School
of Pharmaceutical Sciences, University of
Geneva, Rue Michel Servet 1, 1206 Genève, Switzerland
- Institute
of Pharmaceutical Sciences of Western Switzerland (ISPSO), University of Geneva, 1206 Genève, Switzerland
- Swiss
Institute of Bioinformatics, University
of Geneva, 1206 Genève, Switzerland
| | - Narjes Ansari
- Atomistic
Simulations, Italian Institute of Technology, Via Enrico Melen 83, 16152 Genova, Italy
| | - Francesco Luigi Gervasio
- School
of Pharmaceutical Sciences, University of
Geneva, Rue Michel Servet 1, 1206 Genève, Switzerland
- Institute
of Pharmaceutical Sciences of Western Switzerland (ISPSO), University of Geneva, 1206 Genève, Switzerland
- Swiss
Institute of Bioinformatics, University
of Geneva, 1206 Genève, Switzerland
- Department
of Chemistry, University College London, WC1E 6BT London, U.K.
| |
Collapse
|
24
|
Dandekar BR, Majumdar BB, Mondal J. Nonmonotonic Modulation of the Protein-Ligand Recognition Event by Inert Crowders. J Phys Chem B 2023; 127:7449-7461. [PMID: 37590118 DOI: 10.1021/acs.jpcb.3c03946] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 08/19/2023]
Abstract
The ubiquitous event of a protein recognizing small molecules or ligands at its native binding site is crucial for initiating major biological processes. However, how a crowded environment, as is typically represented by a cellular interior, would modulate the protein-ligand search process is largely debated. Excluded volume-based theory suggests that the presence of an inert crowder would reinforce a steady stabilization and enhancement of the protein-ligand recognition process. Here, we counter this long-held perspective via the molecular dynamics simulation and Markov state model of the protein-ligand recognition event in the presence of inert crowders. Specifically, we demonstrate that, depending on concentration, even purely inert crowders can exert a nonmonotonic effect via either stabilizing or destabilizing the protein-ligand binding event. Analysis of the kinetic network of binding pathways reveals that the crowders would either modulate precedent non-native on-pathway intermediates or would devise additional ones in a multistate recognition event across a wide range of concentrations. As an important insight, crowders gradually shift the relative transitional preference of these intermediates toward a native-bound state, with ligand residence time at the binding pocket dictating the trend of nonmonotonic concentration dependence by simple inert crowders.
Collapse
|
25
|
Uppuladinne MVN, Achalere A, Sonavane U, Joshi R. Probing the structure of human tRNA 3Lys in the presence of ligands using docking, MD simulations and MSM analysis. RSC Adv 2023; 13:25778-25796. [PMID: 37655355 PMCID: PMC10467029 DOI: 10.1039/d3ra03694d] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/02/2023] [Accepted: 08/14/2023] [Indexed: 09/02/2023] Open
Abstract
The tRNA3Lys, which acts as a primer for human immunodeficiency virus type 1 (HIV-1) reverse transcription, undergoes structural changes required for the formation of a primer-template complex. Small molecules have been targeted against tRNA3Lys to inhibit the primer-template complex formation. The present study aims to understand the kinetics of the conformational landscape spanned by tRNA3Lys in apo form using molecular dynamics simulations and Markov state modeling. The study is taken further to investigate the effect of small molecules like 1,4T and 1,5T on structural conformations and kinetics of tRNA3Lys, and comparative analysis is presented. Markov state modeling of tRNA3Lys apo resulted in three metastable states where the conformations have shown the non-canonical structures of the anticodon loop. Based on analyses of ligand-tRNA3Lys interactions, crucial ion and water mediated H-bonds and free energy calculations, it was observed that the 1,4-triazole more strongly binds to the tRNA3Lys compared to 1,5-triazole. However, the MSM analysis suggest that the 1,5-triazole binding to tRNA3Lys has brought rigidity not only in the binding pocket (TΨC arm, D-TΨC loop) but also in the whole structure of tRNA3Lys. This may affect the easy opening of primer tRNA3Lys required for HIV-1 reverse transcription.
Collapse
Affiliation(s)
- Mallikarjunachari V N Uppuladinne
- High Performance Computing - Medical and Bioinformatics Applications, Centre for Development of Advanced Computing (C-DAC) Panchavati, Pashan Pune India
| | - Archana Achalere
- High Performance Computing - Medical and Bioinformatics Applications, Centre for Development of Advanced Computing (C-DAC) Panchavati, Pashan Pune India
| | - Uddhavesh Sonavane
- High Performance Computing - Medical and Bioinformatics Applications, Centre for Development of Advanced Computing (C-DAC) Panchavati, Pashan Pune India
| | - Rajendra Joshi
- High Performance Computing - Medical and Bioinformatics Applications, Centre for Development of Advanced Computing (C-DAC) Panchavati, Pashan Pune India
| |
Collapse
|
26
|
Tripathi S, Nair NN. Temperature Accelerated Sliced Sampling to Probe Ligand Dissociation from Protein. J Chem Inf Model 2023; 63:5182-5191. [PMID: 37540828 DOI: 10.1021/acs.jcim.3c00376] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 08/06/2023]
Abstract
Modeling ligand unbinding in proteins to estimate the free energy of binding and probing the mechanism presents several challenges. They primarily pertain to the entropic bottlenecks resulting from protein and solvent conformations. While exploring the unbinding processes using enhanced sampling techniques, very long simulations are required to sample all of the conformational states as the system gets trapped in local free energy minima along transverse coordinates. Here, we demonstrate that temperature accelerated sliced sampling (TASS) is an ideal approach to overcome some of the difficulties faced by conventional sampling methods in studying ligand unbinding. Using TASS, we study the unbinding of avibactam inhibitor molecules from the Class C β-lactamase (CBL) active site. Extracting CBL-avibactam unbinding free energetics, unbinding pathways, and identifying critical interactions from the TASS simulations are demonstrated.
Collapse
Affiliation(s)
- Shubhandra Tripathi
- Department of Chemistry, Indian Institute of Technology Kanpur, Kanpur 208016, India
| | - Nisanth N Nair
- Department of Chemistry, Indian Institute of Technology Kanpur, Kanpur 208016, India
| |
Collapse
|
27
|
Kozlowski N, Grubmüller H. Uncertainties in Markov State Models of Small Proteins. J Chem Theory Comput 2023; 19:5516-5524. [PMID: 37540193 PMCID: PMC10448719 DOI: 10.1021/acs.jctc.3c00372] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/03/2023] [Indexed: 08/05/2023]
Abstract
Markov state models are widely used to describe and analyze protein dynamics based on molecular dynamics simulations, specifically to extract functionally relevant characteristic time scales and motions. Particularly for larger biomolecules such as proteins, however, insufficient sampling is a notorious concern and often the source of large uncertainties that are difficult to quantify. Furthermore, there are several other sources of uncertainty, such as choice of the number of Markov states and lag time, choice and parameters of dimension reduction preprocessing step, and uncertainty due to the limited number of observed transitions; the latter is often estimated via a Bayesian approach. Here, we quantified and ranked all of these uncertainties for four small globular test proteins. We found that the largest uncertainty is due to insufficient sampling and initially increases with the total trajectory length T up to a critical tipping point, after which it decreases as 1 / T , thus providing guidelines for how much sampling is required for given accuracy. We also found that single long trajectories yielded better sampling accuracy than many shorter trajectories starting from the same structure. In comparison, the remaining sources of the above uncertainties are generally smaller by a factor of about 5, rendering them less of a concern but certainly not negligible. Importantly, the Bayes uncertainty, commonly used as the only uncertainty estimate, captures only a relatively small part of the true uncertainty, which is thus often drastically underestimated.
Collapse
Affiliation(s)
- Nicolai Kozlowski
- Department of Theoretical and Computational
Biophysics, Max-Planck-Institute for Multidisciplinary
Sciences, Göttingen 37077, Germany
| | - Helmut Grubmüller
- Department of Theoretical and Computational
Biophysics, Max-Planck-Institute for Multidisciplinary
Sciences, Göttingen 37077, Germany
| |
Collapse
|
28
|
Smith LG, Novak B, Osato M, Mobley DL, Bowman GR. PopShift: A thermodynamically sound approach to estimate binding free energies by accounting for ligand-induced population shifts from a ligand-free MSM. BIORXIV : THE PREPRINT SERVER FOR BIOLOGY 2023:2023.07.14.549110. [PMID: 37503302 PMCID: PMC10370083 DOI: 10.1101/2023.07.14.549110] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 07/29/2023]
Abstract
Obtaining accurate binding free energies from in silico screens has been a longstanding goal for the computational chemistry community. However, accuracy and computational cost are at odds with one another, limiting the utility of methods that perform this type of calculation. Many methods achieve massive scale by explicitly or implicitly assuming that the target protein adopts a single structure, or undergoes limited fluctuations around that structure, to minimize computational cost. Others simulate each protein-ligand complex of interest, accepting lower throughput in exchange for better predictions of binding affinities. Here, we present the PopShift framework for accounting for the ensemble of structures a protein adopts and their relative probabilities. Protein degrees of freedom are enumerated once, and then arbitrarily many molecules can be screened against this ensemble. Specifically, we use Markov state models (MSMs) as a compressed representation of a protein's thermodynamic ensemble. We start with a ligand-free MSM and then calculate how addition of a ligand shifts the populations of each protein conformational state based on the strength of the interaction between that protein conformation and the ligand. In this work we use docking to estimate the affinity between a given protein structure and ligand, but any estimator of binding affinities could be used in the PopShift framework. We test PopShift on the classic benchmark pocket T4 Lysozyme L99A. We find that PopShift is more accurate than common strategies, such as docking to a single structure and traditional ensemble docking-producing results that compare favorably with alchemical binding free energy calculations in terms of RMSE but not correlation - and may have a more favorable computational cost profile in some applications. In addition to predicting binding free energies and ligand poses, PopShift also provides insight into how the probability of different protein structures is shifted upon addition of various concentrations of ligand, providing a platform for predicting affinities and allosteric effects of ligand binding. Therefore, we expect PopShift will be valuable for hit finding and for providing insight into phenomena like allostery.
Collapse
Affiliation(s)
- Louis G Smith
- University of Pennsylvania, Depts. of Biochemistry & Biophysics and Bioengineering
| | - Borna Novak
- Washington University in St. Louis, Department of Biochemistry and Molecular Biophysics
- Medical Scientist Training Program, Washington University in St. Louis
| | - Meghan Osato
- University of California Irvine, School of Pharmacy and Pharmaceutical Sciences
| | - David L Mobley
- University of California Irvine, School of Pharmacy and Pharmaceutical Sciences
| | - Gregory R Bowman
- University of Pennsylvania, Depts. of Biochemistry & Biophysics and Bioengineering
| |
Collapse
|
29
|
Sisk T, Robustelli P. Folding-upon-binding pathways of an intrinsically disordered protein from a deep Markov state model. BIORXIV : THE PREPRINT SERVER FOR BIOLOGY 2023:2023.07.21.550103. [PMID: 37546728 PMCID: PMC10401938 DOI: 10.1101/2023.07.21.550103] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 08/08/2023]
Abstract
A central challenge in the study of intrinsically disordered proteins is the characterization of the mechanisms by which they bind their physiological interaction partners. Here, we utilize a deep learning based Markov state modeling approach to characterize the folding-upon-binding pathways observed in a long-time scale molecular dynamics simulation of a disordered region of the measles virus nucleoprotein NTAIL reversibly binding the X domain of the measles virus phosphoprotein complex. We find that folding-upon-binding predominantly occurs via two distinct encounter complexes that are differentiated by the binding orientation, helical content, and conformational heterogeneity of NTAIL. We do not, however, find evidence for the existence of canonical conformational selection or induced fit binding pathways. We observe four kinetically separated native-like bound states that interconvert on time scales of eighty to five hundred nanoseconds. These bound states share a core set of native intermolecular contacts and stable NTAIL helices and are differentiated by a sequential formation of native and non-native contacts and additional helical turns. Our analyses provide an atomic resolution structural description of intermediate states in a folding-upon-binding pathway and elucidate the nature of the kinetic barriers between metastable states in a dynamic and heterogenous, or "fuzzy", protein complex.
Collapse
Affiliation(s)
- Thomas Sisk
- Dartmouth College, Department of Chemistry, Hanover, NH, 03755
| | - Paul Robustelli
- Dartmouth College, Department of Chemistry, Hanover, NH, 03755
| |
Collapse
|
30
|
Qiu Y, O’Connor MS, Xue M, Liu B, Huang X. An Efficient Path Classification Algorithm Based on Variational Autoencoder to Identify Metastable Path Channels for Complex Conformational Changes. J Chem Theory Comput 2023; 19:4728-4742. [PMID: 37382437 PMCID: PMC11042546 DOI: 10.1021/acs.jctc.3c00318] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/30/2023]
Abstract
Conformational changes (i.e., dynamic transitions between pairs of conformational states) play important roles in many chemical and biological processes. Constructing the Markov state model (MSM) from extensive molecular dynamics (MD) simulations is an effective approach to dissect the mechanism of conformational changes. When combined with transition path theory (TPT), MSM can be applied to elucidate the ensemble of kinetic pathways connecting pairs of conformational states. However, the application of TPT to analyze complex conformational changes often results in a vast number of kinetic pathways with comparable fluxes. This obstacle is particularly pronounced in heterogeneous self-assembly and aggregation processes. The large number of kinetic pathways makes it challenging to comprehend the molecular mechanisms underlying conformational changes of interest. To address this challenge, we have developed a path classification algorithm named latent-space path clustering (LPC) that efficiently lumps parallel kinetic pathways into distinct metastable path channels, making them easier to comprehend. In our algorithm, MD conformations are first projected onto a low-dimensional space containing a small set of collective variables (CVs) by time-structure-based independent component analysis (tICA) with kinetic mapping. Then, MSM and TPT are constructed to obtain the ensemble of pathways, and a deep learning architecture named the variational autoencoder (VAE) is used to learn the spatial distributions of kinetic pathways in the continuous CV space. Based on the trained VAE model, the TPT-generated ensemble of kinetic pathways can be embedded into a latent space, where the classification becomes clear. We show that LPC can efficiently and accurately identify the metastable path channels in three systems: a 2D potential, the aggregation of two hydrophobic particles in water, and the folding of the Fip35 WW domain. Using the 2D potential, we further demonstrate that our LPC algorithm outperforms the previous path-lumping algorithms by making substantially fewer incorrect assignments of individual pathways to four path channels. We expect that LPC can be widely applied to identify the dominant kinetic pathways underlying complex conformational changes.
Collapse
Affiliation(s)
- Yunrui Qiu
- Department of Chemistry, Theoretical Chemistry Institute, University of Wisconsin-Madison, Madison, WI, 53706, USA
| | - Michael S. O’Connor
- Biophysics Graduate Program, University of Wisconsin-Madison, Madison, WI, 53706, USA
| | - Mingyi Xue
- Department of Chemistry, Theoretical Chemistry Institute, University of Wisconsin-Madison, Madison, WI, 53706, USA
| | - Bojun Liu
- Department of Chemistry, Theoretical Chemistry Institute, University of Wisconsin-Madison, Madison, WI, 53706, USA
| | - Xuhui Huang
- Department of Chemistry, Theoretical Chemistry Institute, University of Wisconsin-Madison, Madison, WI, 53706, USA
- Biophysics Graduate Program, University of Wisconsin-Madison, Madison, WI, 53706, USA
| |
Collapse
|
31
|
Wong CF. 15 Years of molecular simulation of drug-binding kinetics. Expert Opin Drug Discov 2023; 18:1333-1348. [PMID: 37789731 PMCID: PMC10926948 DOI: 10.1080/17460441.2023.2264770] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/20/2023] [Accepted: 09/26/2023] [Indexed: 10/05/2023]
Abstract
INTRODUCTION Drug-binding kinetics has been increasingly recognized as an important factor to be considered in drug discovery. Long residence time could prolong the action of some drugs while produce toxicity on others. Early evaluation of the binding kinetics of drug candidates could reduce attrition rate late in the drug discovery process. Computational prediction of drug-binding kinetics is useful as compounds can be evaluated even before they are made. However, simulation of drug-binding kinetics is a challenging problem because of the long-time scale involved. Nevertheless, significant progress has been made. AREAS COVERED This review illustrates the rapid evolution of qualitative to quantitative molecular dynamics-based methods that have been developed over the last 15 years. EXPERT OPINION The development of new methods based on molecular dynamics simulations now enables computation of absolute association/dissociation rate constants. Cheaper methods capable of identifying candidates with fast or slow binding kinetics, or rank-ordering rate constants are also available. Together, these methods have generated useful insights into the molecular mechanisms of drug-binding kinetics, and the design of drug candidates with therapeutically favorable kinetics. Although predicting absolute rate constants is still expensive and challenging, rapid improvement is expected in the coming years with the continuing refinement of current technologies, development of new methodologies, and the utilization of machine learning.
Collapse
Affiliation(s)
- Chung F Wong
- Department of Chemistry and Biochemistry, University of Missouri-St. Louis, St. Louis, MO, USA
| |
Collapse
|
32
|
Dandekar B, Ahalawat N, Sinha S, Mondal J. Markov State Models Reconcile Conformational Plasticity of GTPase with Its Substrate Binding Event. JACS AU 2023; 3:1728-1741. [PMID: 37388689 PMCID: PMC10302740 DOI: 10.1021/jacsau.3c00151] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 03/27/2023] [Revised: 05/01/2023] [Accepted: 05/02/2023] [Indexed: 07/01/2023]
Abstract
Ras GTPase is an enzyme that catalyzes the hydrolysis of guanosine triphosphate (GTP) and plays an important role in controlling crucial cellular signaling pathways. However, this enzyme has always been believed to be undruggable due to its strong binding affinity with its native substrate GTP. To understand the potential origin of high GTPase/GTP recognition, here we reconstruct the complete process of GTP binding to Ras GTPase via building Markov state models (MSMs) using a 0.1 ms long all-atom molecular dynamics (MD) simulation. The kinetic network model, derived from the MSM, identifies multiple pathways of GTP en route to its binding pocket. While the substrate stalls onto a set of non-native metastable GTPase/GTP encounter complexes, the MSM accurately discovers the native pose of GTP at its designated catalytic site in crystallographic precision. However, the series of events exhibit signatures of conformational plasticity in which the protein remains trapped in multiple non-native conformations even when GTP has already located itself in its native binding site. The investigation demonstrates mechanistic relays pertaining to simultaneous fluctuations of switch 1 and switch 2 residues which remain most instrumental in maneuvering the GTP-binding process. Scanning of the crystallographic database reveals close resemblance between observed non-native GTP binding poses and precedent crystal structures of substrate-bound GTPase, suggesting potential roles of these binding-competent intermediates in allosteric regulation of the recognition process.
Collapse
Affiliation(s)
| | - Navjeet Ahalawat
- Department
of Bioinformatics and Computational Biology, College of Biotechnology, CCS Haryana Agricultural University, Hisar, 125004 Haryana, India
| | - Suman Sinha
- Institute
of Pharmaceutical Research, GLA University, Mathura, 281406 Uttar Pradesh, India
| | - Jagannath Mondal
- Tata
Institute of Fundamental Research, Hyderabad, Telangana 500046, India
| |
Collapse
|
33
|
Hardie A, Cossins BP, Lovera S, Michel J. Deconstructing allostery by computational assessment of the binding determinants of allosteric PTP1B modulators. Commun Chem 2023; 6:125. [PMID: 37322137 DOI: 10.1038/s42004-023-00926-1] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/13/2023] [Accepted: 06/08/2023] [Indexed: 06/17/2023] Open
Abstract
Fragment-based drug discovery is an established methodology for finding hit molecules that can be elaborated into lead compounds. However it is currently challenging to predict whether fragment hits that do not bind to an orthosteric site could be elaborated into allosteric modulators, as in these cases binding does not necessarily translate into a functional effect. We propose a workflow using Markov State Models (MSMs) with steered molecular dynamics (sMD) to assess the allosteric potential of known binders. sMD simulations are employed to sample protein conformational space inaccessible to routine equilibrium MD timescales. Protein conformations sampled by sMD provide starting points for seeded MD simulations, which are combined into MSMs. The methodology is demonstrated on a dataset of protein tyrosine phosphatase 1B ligands. Experimentally confirmed allosteric inhibitors are correctly classified as inhibitors, whereas the deconstructed analogues show reduced inhibitory activity. Analysis of the MSMs provide insights into preferred protein-ligand arrangements that correlate with functional outcomes. The present methodology may find applications for progressing fragments towards lead molecules in FBDD campaigns.
Collapse
Affiliation(s)
- Adele Hardie
- EaStChem School of Chemistry, Joseph Black Building, University of Edinburgh, Edinburgh, EH9 3FJ, UK
| | - Benjamin P Cossins
- UCB Pharma, 216 Bath Road, Slough, UK
- Exscientia, The Schrödinger Building, Oxford Science Park, Oxford, UK
| | - Silvia Lovera
- UCB Pharma, Chemin du Foriest 1, 1420, Braine-l'Alleud, Belgium
| | - Julien Michel
- EaStChem School of Chemistry, Joseph Black Building, University of Edinburgh, Edinburgh, EH9 3FJ, UK.
| |
Collapse
|
34
|
Dominic AJ, Cao S, Montoya-Castillo A, Huang X. Memory Unlocks the Future of Biomolecular Dynamics: Transformative Tools to Uncover Physical Insights Accurately and Efficiently. J Am Chem Soc 2023; 145:9916-9927. [PMID: 37104720 DOI: 10.1021/jacs.3c01095] [Citation(s) in RCA: 4] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/29/2023]
Abstract
Conformational changes underpin function and encode complex biomolecular mechanisms. Gaining atomic-level detail of how such changes occur has the potential to reveal these mechanisms and is of critical importance in identifying drug targets, facilitating rational drug design, and enabling bioengineering applications. While the past two decades have brought Markov state model techniques to the point where practitioners can regularly use them to glimpse the long-time dynamics of slow conformations in complex systems, many systems are still beyond their reach. In this Perspective, we discuss how including memory (i.e., non-Markovian effects) can reduce the computational cost to predict the long-time dynamics in these complex systems by orders of magnitude and with greater accuracy and resolution than state-of-the-art Markov state models. We illustrate how memory lies at the heart of successful and promising techniques, ranging from the Fokker-Planck and generalized Langevin equations to deep-learning recurrent neural networks and generalized master equations. We delineate how these techniques work, identify insights that they can offer in biomolecular systems, and discuss their advantages and disadvantages in practical settings. We show how generalized master equations can enable the investigation of, for example, the gate-opening process in RNA polymerase II and demonstrate how our recent advances tame the deleterious influence of statistical underconvergence of the molecular dynamics simulations used to parameterize these techniques. This represents a significant leap forward that will enable our memory-based techniques to interrogate systems that are currently beyond the reach of even the best Markov state models. We conclude by discussing some current challenges and future prospects for how exploiting memory will open the door to many exciting opportunities.
Collapse
Affiliation(s)
- Anthony J Dominic
- Department of Chemistry, University of Colorado Boulder, Boulder, Colorado 80309, USA
| | - Siqin Cao
- Department of Chemistry, Theoretical Chemistry Institute, University of Wisconsin-Madison, Madison, Wisconsin 53706, USA
| | | | - Xuhui Huang
- Department of Chemistry, Theoretical Chemistry Institute, University of Wisconsin-Madison, Madison, Wisconsin 53706, USA
| |
Collapse
|
35
|
Gupta MN, Uversky VN. Moonlighting enzymes: when cellular context defines specificity. Cell Mol Life Sci 2023; 80:130. [PMID: 37093283 PMCID: PMC11073002 DOI: 10.1007/s00018-023-04781-0] [Citation(s) in RCA: 10] [Impact Index Per Article: 10.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/30/2023] [Revised: 03/13/2023] [Accepted: 04/15/2023] [Indexed: 04/25/2023]
Abstract
It is not often realized that the absolute protein specificity is an exception rather than a rule. Two major kinds of protein multi-specificities are promiscuity and moonlighting. This review discusses the idea of enzyme specificity and then focusses on moonlighting. Some important examples of protein moonlighting, such as crystallins, ceruloplasmin, metallothioniens, macrophage migration inhibitory factor, and enzymes of carbohydrate metabolism are discussed. How protein plasticity and intrinsic disorder enable the removing the distinction between enzymes and other biologically active proteins are outlined. Finally, information on important roles of moonlighting in human diseases is updated.
Collapse
Affiliation(s)
- Munishwar Nath Gupta
- Department of Biochemical Engineering and Biotechnology, Indian Institute of Technology, Hauz Khas, New Delhi, 110016, India
| | - Vladimir N Uversky
- Department of Molecular Medicine and USF Health Byrd Alzheimer's Research Institute, Morsani College of Medicine, University of South Florida, 12901 Bruce B. Downs Blvd., MDC07, Tampa, FL, 33612-4799, USA.
| |
Collapse
|
36
|
Ahalawat N, Sahil M, Mondal J. Resolving Protein Conformational Plasticity and Substrate Binding via Machine Learning. J Chem Theory Comput 2023; 19:2644-2657. [PMID: 37068044 DOI: 10.1021/acs.jctc.2c00932] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/18/2023]
Abstract
A long-standing target in elucidating the biomolecular recognition process is the identification of binding-competent conformations of the receptor protein. However, protein conformational plasticity and the stochastic nature of the recognition processes often preclude the assignment of a specific protein conformation to an individual ligand-bound pose. Here, we demonstrate that a computational framework coined as RF-TICA-MD, which integrates an ensemble decision-tree-based Random Forest (RF) machine learning (ML) technique with an unsupervised dimension reduction approach time-structured independent component analysis (TICA), provides an efficient and unambiguous solution toward resolving protein conformational plasticity and the substrate binding process. In particular, we consider multimicrosecond-long molecular dynamics (MD) simulation trajectories of a ligand recognition process in solvent-inaccessible cavities of archetypal proteins T4 lysozyme and cytochrome P450cam. We show that in a scenario in which clear correspondence between protein conformation and binding-competent macrostates could not be obtained via an unsupervised dimension reduction approach, an a priori decision-tree-based supervised classification of the simulated recognition trajectories via RF would help characterize key amino acid residue pairs of the protein that are deemed sensitive for ligand binding. A subsequent unsupervised dimensional reduction of the selected residue pairs via TICA would then delineate a conformational landscape of protein which is able to demarcate ligand-bound poses from unbound ones. The proposed RF-TICA-MD approach is shown to be data agnostic and found to be robust when using other ML-based classification methods such as XGBoost. As a promising spinoff of the protocol, the framework is found to be capable of identifying distal protein locations which would be allosterically important for ligand binding and would characterize their roles in recognition pathways. A Python implementation of a proposed ML workflow is available in GitHub https://github.com/navjeet0211/rf-tica-md.
Collapse
Affiliation(s)
- Navjeet Ahalawat
- Department of Bioinformatics and Computational Biology, College of Biotechnology, CCS Haryana Agricultural University, Hisar 125 004, Haryana, India
| | - Mohammad Sahil
- Center for Interdisciplinary Sciences, Tata Institute of Fundamental Research, Hyderabad 500046, India
| | - Jagannath Mondal
- Center for Interdisciplinary Sciences, Tata Institute of Fundamental Research, Hyderabad 500046, India
| |
Collapse
|
37
|
Ojha AA, Srivastava A, Votapka LW, Amaro RE. Selectivity and Ranking of Tight-Binding JAK-STAT Inhibitors Using Markovian Milestoning with Voronoi Tessellations. J Chem Inf Model 2023; 63:2469-2482. [PMID: 37023323 PMCID: PMC10131228 DOI: 10.1021/acs.jcim.2c01589] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/08/2023]
Abstract
Janus kinases (JAK), a group of proteins in the nonreceptor tyrosine kinase (NRTKs) family, play a crucial role in growth, survival, and angiogenesis. They are activated by cytokines through the Janus kinase-signal transducer and activator of a transcription (JAK-STAT) signaling pathway. JAK-STAT signaling pathways have significant roles in the regulation of cell division, apoptosis, and immunity. Identification of the V617F mutation in the Janus homology 2 (JH2) domain of JAK2 leading to myeloproliferative disorders has stimulated great interest in the drug discovery community to develop JAK2-specific inhibitors. However, such inhibitors should be selective toward JAK2 over other JAKs and display an extended residence time. Recently, novel JAK2/STAT5 axis inhibitors (N-(1H-pyrazol-3-yl)pyrimidin-2-amino derivatives) have displayed extended residence times (hours or longer) on target and adequate selectivity excluding JAK3. To facilitate a deeper understanding of the kinase-inhibitor interactions and advance the development of such inhibitors, we utilize a multiscale Markovian milestoning with Voronoi tessellations (MMVT) approach within the Simulation-Enabled Estimation of Kinetic Rates v.2 (SEEKR2) program to rank order these inhibitors based on their kinetic properties and further explain the selectivity of JAK2 inhibitors over JAK3. Our approach investigates the kinetic and thermodynamic properties of JAK-inhibitor complexes in a user-friendly, fast, efficient, and accurate manner compared to other brute force and hybrid-enhanced sampling approaches.
Collapse
Affiliation(s)
- Anupam Anand Ojha
- Department of Chemistry and Biochemistry, University of California San Diego, La Jolla, California 92093, United States
| | - Ambuj Srivastava
- Department of Chemistry and Biochemistry, University of California San Diego, La Jolla, California 92093, United States
| | - Lane William Votapka
- Department of Chemistry and Biochemistry, University of California San Diego, La Jolla, California 92093, United States
| | - Rommie E Amaro
- Department of Chemistry and Biochemistry, University of California San Diego, La Jolla, California 92093, United States
| |
Collapse
|
38
|
Wang J, Do HN, Koirala K, Miao Y. Predicting Biomolecular Binding Kinetics: A Review. J Chem Theory Comput 2023; 19:2135-2148. [PMID: 36989090 DOI: 10.1021/acs.jctc.2c01085] [Citation(s) in RCA: 8] [Impact Index Per Article: 8.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/30/2023]
Abstract
Biomolecular binding kinetics including the association (kon) and dissociation (koff) 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.
Collapse
Affiliation(s)
- Jinan Wang
- Center for Computational Biology and Department of Molecular Biosciences, University of Kansas, Lawrence, Kansas 66047, United States
| | - Hung N Do
- Center for Computational Biology and Department of Molecular Biosciences, University of Kansas, Lawrence, Kansas 66047, United States
| | - Kushal Koirala
- Center for Computational Biology and Department of Molecular Biosciences, University of Kansas, Lawrence, Kansas 66047, United States
| | - Yinglong Miao
- Center for Computational Biology and Department of Molecular Biosciences, University of Kansas, Lawrence, Kansas 66047, United States
| |
Collapse
|
39
|
Wolf S, Post M, Stock G. Path separation of dissipation-corrected targeted molecular dynamics simulations of protein-ligand unbinding. J Chem Phys 2023; 158:124106. [PMID: 37003731 DOI: 10.1063/5.0138761] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/03/2023] Open
Abstract
Protein-ligand (un)binding simulations are a recent focus of biased molecular dynamics simulations. Such binding and unbinding can occur via different pathways in and out of a binding site. Here, we present a theoretical framework on how to compute kinetics along separate paths and on how to combine the path-specific rates into global binding and unbinding rates for comparison with experimental results. Using dissipation-corrected targeted molecular dynamics in combination with temperature-boosted Langevin equation simulations [S. Wolf et al., Nat. Commun. 11, 2918 (2020)] applied to a two-dimensional model and the trypsin-benzamidine complex as test systems, we assess the robustness of the procedure and discuss the aspects of its practical applicability to predict multisecond kinetics of complex biomolecular systems.
Collapse
Affiliation(s)
- Steffen Wolf
- Biomolecular Dynamics, Institute of Physics, Albert Ludwigs University, 79104 Freiburg, Germany
| | - Matthias Post
- Biomolecular Dynamics, Institute of Physics, Albert Ludwigs University, 79104 Freiburg, Germany
| | - Gerhard Stock
- Biomolecular Dynamics, Institute of Physics, Albert Ludwigs University, 79104 Freiburg, Germany
| |
Collapse
|
40
|
Cui M, Nguyen D, Gaillez MP, Heiden S, Lin W, Thompson M, Reddavide FV, Chen Q, Zhang Y. Trio-pharmacophore DNA-encoded chemical library for simultaneous selection of fragments and linkers. Nat Commun 2023; 14:1481. [PMID: 36932079 PMCID: PMC10023787 DOI: 10.1038/s41467-023-37071-1] [Citation(s) in RCA: 7] [Impact Index Per Article: 7.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/15/2022] [Accepted: 02/28/2023] [Indexed: 03/19/2023] Open
Abstract
The split-and-pool method has been widely used to synthesize chemical libraries of a large size for early drug discovery, albeit without the possibility of meaningful quality control. In contrast, a self-assembled DNA-encoded chemical library (DEL) allows us to construct an m x n-member library by mixing an m-member and an n-member pre-purified sub-library. Herein, we report a trio-pharmacophore DEL (T-DEL) of m x l x n members through assembling three pre-purified and validated sub-libraries. The middle sub-library is synthesized using DNA-templated synthesis with different reaction mechanisms and designed as a linkage connecting the fragments displayed on the flanking two sub-libraries. Despite assembling three fragments, the resulting compounds do not exceed the up-to-date standard of molecular weight regarding drug-likeness. We demonstrate the utility of T-DEL in linker optimization for known binding fragments against trypsin and carbonic anhydrase II and by de novo selections against matrix metalloprotease-2 and -9.
Collapse
Affiliation(s)
- Meiying Cui
- B CUBE, Center for Molecular Bioengineering, Technische Universität Dresden, Dresden, Germany
| | | | - Michelle Patino Gaillez
- B CUBE, Center for Molecular Bioengineering, Technische Universität Dresden, Dresden, Germany
| | | | - Weilin Lin
- B CUBE, Center for Molecular Bioengineering, Technische Universität Dresden, Dresden, Germany
| | | | | | - Qinchang Chen
- Research Institute of Intelligent Computing, Zhejiang Lab, Hangzhou, China.
- School of Life Sciences and Technology, Tongji University, Shanghai, China.
| | - Yixin Zhang
- B CUBE, Center for Molecular Bioengineering, Technische Universität Dresden, Dresden, Germany.
| |
Collapse
|
41
|
Heckmeier PJ, Ruf J, Janković BG, Hamm P. MCL-1 promiscuity and the structural resilience of its binding partners. J Chem Phys 2023; 158:095101. [PMID: 36889945 DOI: 10.1063/5.0137239] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/11/2023] Open
Abstract
The allosteric protein MCL-1 and its natural inhibitors, the BH3-only proteins PUMA, BIM, and NOXA regulate apoptosis by interacting promiscuously within an entangled binding network. Little is known about the transient processes and dynamic conformational fluctuations that are the basis for the formation and stability of the MCL-1/BH3-only complex. In this study, we designed photoswitchable versions of MCL-1/PUMA and MCL-1/NOXA, and investigated the protein response after an ultrafast photo-perturbation with transient infrared spectroscopy. We observed partial α-helical unfolding in all cases, albeit on strongly varying timescales (1.6 ns for PUMA, 9.7 ns for the previously studied BIM, and 85 ns for NOXA). These differences are interpreted as a BH3-only-specific "structural resilience" to defy the perturbation while remaining in MCL-1's binding pocket. Thus, the presented insights could help to better understand the differences between PUMA, BIM, and NOXA, the promiscuity of MCL-1, in general, and the role of the proteins in the apoptotic network.
Collapse
Affiliation(s)
| | - Jeannette Ruf
- Department of Chemistry, University of Zurich, Zurich, Switzerland
| | | | - Peter Hamm
- Department of Chemistry, University of Zurich, Zurich, Switzerland
| |
Collapse
|
42
|
Sahil M, Sarkar S, Mondal J. Long-time-step molecular dynamics can retard simulation of protein-ligand recognition process. Biophys J 2023; 122:802-816. [PMID: 36726313 PMCID: PMC10027446 DOI: 10.1016/j.bpj.2023.01.036] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/08/2022] [Revised: 10/31/2022] [Accepted: 01/25/2023] [Indexed: 02/03/2023] Open
Abstract
Molecular dynamics (MD) simulation of biologically relevant processes at realistic time scale and atomistic precision is generally limited by prohibitively large computational cost, due to its restriction of using an ultrashort integration time step (1-2 fs). A popular numerical recipe to reduce the associated computational burden is adopting schemes that would allow relatively longer-time-step for MD propagation. Here, we explore the perceived potential of one of the most frequently used long-time-step protocols, namely the hydrogen mass repartitioning (HMR) approach, in alleviating the computational overhead associated with simulation of the kinetic process of protein-ligand recognition events. By repartitioning the mass of heavier atoms to their linked hydrogen atoms, HMR leverages around twofold longer time step than regular simulation, holding promise of significant performance boost. However, our probe into direct simulation of the protein-ligand recognition event, one of the computationally most challenging processes, shows that long-time-step HMR MD simulations do not necessarily translate to a computationally affordable solution. Our investigations spanning cumulative 176 μs in three independent proteins (T4 lysozyme, sensor domain of MopR, and galectin-3) show that long-time-step HMR-based MD simulations can catch the ligand in its act of recognizing the native cavity. But, as a major caveat, the ligand is found to require significantly longer time to identify buried native protein cavity in an HMR MD simulation than regular simulation, thereby defeating the purpose of its usage for performance upgrade. A molecular analysis shows that the longer time required by a ligand to recognize the protein in HMR is rooted in faster diffusion of the ligand, which reduces the survival probability of decisive on-pathway metastable intermediates, thereby slowing down the eventual recognition process at the native cavity. Together, the investigation stresses careful assessment of pitfalls of long-time-step algorithms before attempting to utilize them for higher performance for biomolecular recognition simulations.
Collapse
Affiliation(s)
- Mohammad Sahil
- Tata Institute of Fundamental Research, Hyderabad 500046, India
| | - Susmita Sarkar
- Tata Institute of Fundamental Research, Hyderabad 500046, India
| | - Jagannath Mondal
- Tata Institute of Fundamental Research, Hyderabad 500046, India.
| |
Collapse
|
43
|
Ligand Gaussian Accelerated Molecular Dynamics 2 (LiGaMD2): Improved Calculations of Ligand Binding Thermodynamics and Kinetics with Closed Protein Pocket. J Chem Theory Comput 2023; 19:733-745. [PMID: 36706316 DOI: 10.1021/acs.jctc.2c01194] [Citation(s) in RCA: 12] [Impact Index Per Article: 12.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/28/2023]
Abstract
Ligand binding thermodynamics and kinetics are critical parameters for drug design. However, it has proven challenging to efficiently predict ligand binding thermodynamics and kinetics from molecular simulations due to limited simulation timescales. Protein dynamics, especially in the ligand binding pocket, often plays an important role in ligand binding. Based on our previously developed Ligand Gaussian accelerated molecular dynamics (LiGaMD), here we present LiGaMD2 in which a selective boost potential was applied to both the ligand and protein residues in the binding pocket to improve sampling of ligand binding and dissociation. To validate the performance of LiGaMD2, the T4 lysozyme (T4L) mutants with open and closed pockets bound by different ligands were chosen as model systems. LiGaMD2 could efficiently capture repetitive ligand dissociation and binding within microsecond simulations of all T4L systems. The obtained ligand binding kinetic rates and free energies agreed well with available experimental values and previous modeling results. Therefore, LiGaMD2 provides an improved approach to sample opening of closed protein pockets for ligand dissociation and binding, thereby allowing for efficient calculations of ligand binding thermodynamics and kinetics.
Collapse
|
44
|
Xu D, Li C, Li W, Lin B, Lv R. Recent advances in lanthanide-doped up-conversion probes for theranostics. Front Chem 2023; 11:1036715. [PMID: 36846851 PMCID: PMC9949555 DOI: 10.3389/fchem.2023.1036715] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/05/2022] [Accepted: 01/31/2023] [Indexed: 02/11/2023] Open
Abstract
Up-conversion (or anti-Stokes) luminescence refers to the phenomenon whereby materials emit high energy, short-wavelength light upon excitation at longer wavelengths. Lanthanide-doped up-conversion nanoparticles (Ln-UCNPs) are widely used in biomedicine due to their excellent physical and chemical properties such as high penetration depth, low damage threshold and light conversion ability. Here, the latest developments in the synthesis and application of Ln-UCNPs are reviewed. First, methods used to synthesize Ln-UCNPs are introduced, and four strategies for enhancing up-conversion luminescence are analyzed, followed by an overview of the applications in phototherapy, bioimaging and biosensing. Finally, the challenges and future prospects of Ln-UCNPs are summarized.
Collapse
Affiliation(s)
| | | | | | - Bi Lin
- Engineering Research Center of Molecular and Neuro Imaging, Ministry of Education, School of Life Science and Technology, Xidian University, Xi’an, Shaanxi, China
| | | |
Collapse
|
45
|
Sobecks BL, Chen J, Shukla D. Mechanistic Basis for Enhanced Strigolactone Sensitivity in KAI2 Triple Mutant. BIORXIV : THE PREPRINT SERVER FOR BIOLOGY 2023:2023.01.18.524622. [PMID: 36712135 PMCID: PMC9882355 DOI: 10.1101/2023.01.18.524622] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Indexed: 01/22/2023]
Abstract
Striga hermonthica is a parasitic weed that destroys billions of dollars' worth of staple crops every year. Its rapid proliferation stems from an enhanced ability to metabolize strigolactones (SLs), plant hormones that direct root branching and shoot growth. Striga's SL receptor, ShHTL7, bears more similarity to the staple crop karrikin receptor KAI2 than to SL receptor D14, though KAI2 variants in plants like Arabidopsis thaliana show minimal SL sensitivity. Recently, studies have indicated that a small number of point mutations to HTL7 residues can confer SL sensitivity to AtKAI2. Here, we analyze both wild-type AtKAI2 and SL-sensitive mutant Var64 through all-atom, long-timescale molecular dynamics simulations to determine the effects of these mutations on receptor function at a molecular level. We demonstrate that the mutations stabilize SL binding by about 2 kcal/mol. They also result in a doubling of the average pocket volume, and eliminate the dependence of binding on certain pocket conformational arrangements. While the probability of certain non-binding SL-receptor interactions increases in the mutant compared with the wild-type, the rate of binding also increases by a factor of ten. All these changes account for the increased SL sensitivity in mutant KAI2, and suggest mechanisms for increasing functionality of host crop SL receptors.
Collapse
Affiliation(s)
- Briana L Sobecks
- Department of Chemical and Biomolecular Engineering, University of Illinois at Urbana-Champaign, Urbana, IL, 61801
| | - Jiming Chen
- Department of Chemical and Biomolecular Engineering, University of Illinois at Urbana-Champaign, Urbana, IL, 61801
| | - Diwakar Shukla
- Department of Chemical and Biomolecular Engineering, University of Illinois at Urbana-Champaign, Urbana, IL, 61801
- Center for Biophysics and Quantitative Biology, University of Illinois at Urbana-Champaign, Urbana, IL, 61801
- Department of Bioengineering, University of Illinois at Urbana-Champaign, Urbana, IL, 61801
- Department of Plant Biology, University of Illinois at Urbana-Champaign, Urbana, IL, 61801
| |
Collapse
|
46
|
Lu C, Peng X, Lu D. Molecular Dynamics Simulation of Protein Cages. Methods Mol Biol 2023; 2671:273-305. [PMID: 37308651 DOI: 10.1007/978-1-0716-3222-2_16] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/14/2023]
Abstract
Molecular dynamics (MD) simulations enable the description of the physical movement of the system over time based on classical mechanics at various scales depending on the models. Protein cages are a particular group of different-size proteins with hollow, spherical structures and are widely found in nature, which have vast applications in numerous fields. The MD simulation of cage proteins is particularly important as a powerful tool to unveil their structures and dynamics for various properties, assembly behavior, and molecular transport mechanisms. Here, we describe how to conduct MD simulations for cage proteins, especially technical details, and analyze some of the properties of interest using GROMACS/NAMD packages.
Collapse
Affiliation(s)
- Chenlin Lu
- Department of Chemical Engineering, Tsinghua University, Beijing, China
| | - Xue Peng
- Department of Chemical Engineering, Tsinghua University, Beijing, China
| | - Diannan Lu
- Department of Chemical Engineering, Tsinghua University, Beijing, China.
| |
Collapse
|
47
|
Aliakbar Tehrani Z, Rulíšek L, Černý J. Molecular dynamics simulations provide structural insight into binding of cyclic dinucleotides to human STING protein. J Biomol Struct Dyn 2022; 40:10250-10264. [PMID: 34187319 DOI: 10.1080/07391102.2021.1942213] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/16/2022]
Abstract
Human stimulator of interferon genes (hSTING) is a signaling adaptor protein that triggers innate immune system by response to cytosolic DNA and second messenger cyclic dinucleotides (CDNs). Natural CDNs contain purine nucleobase with different phosphodiester linkage types (3'-3', 2'-2' or mixed 2'-3'-linkages) and exhibit different binding affinity towards hSTING, ranging from micromolar to nanomolar. High-affinity CDNs are considered as suitable candidates for treatment of chronic hepatitis B and cancer. We have used molecular dynamics simulations to investigate dynamical aspects of binding of natural CDNs (specifically, 2'-2'-cGAMP, 2'-3'-cGAMP, 3'-3'-cGAMP, 3'-3'-c-di-AMP, and 3'-3'-c-di-GMP) with hSTINGwt protein. Our results revealed that CDN/hSTINGwt interactions are controlled by the balance between fluctuations (conformational changes) in the CDN ligand and the protein dynamics. Binding of different CDNs induces different degrees of conformational/dynamics changes in hSTINGwt ligand binding cavity, especially in α1-helices, the so-called lid region and α2-tails. The ligand residence time in hSTINGwt protein pocket depends on different contribution of R232 and R238 residues interacting with oxygen atoms of phosphodiester groups in ligand, water distribution around interacting charged centers (in protein residues and ligand) and structural stability of closed conformation state of hSTINGwt protein. These findings may perhaps guide design of new compounds modulating hSTING activity.Communicated by Ramaswamy H. Sarma.
Collapse
Affiliation(s)
- Zahra Aliakbar Tehrani
- Laboratory of Structural Bioinformatics of Proteins, Institute of Biotechnology of the Czech Academy of Sciences, BIOCEV, Vestec, Czech Republic
| | - Lubomír Rulíšek
- Theoretical Bioinorganic Chemistry, Institute of Organic Chemistry and Biochemistry of the Czech Academy of Sciences, Prague, Czech Republic
| | - Jiří Černý
- Laboratory of Structural Bioinformatics of Proteins, Institute of Biotechnology of the Czech Academy of Sciences, BIOCEV, Vestec, Czech Republic
| |
Collapse
|
48
|
Wehrhan L, Leppkes J, Dimos N, Loll B, Koksch B, Keller BG. Water Network in the Binding Pocket of Fluorinated BPTI-Trypsin Complexes─Insights from Simulation and Experiment. J Phys Chem B 2022; 126:9985-9999. [PMID: 36409613 DOI: 10.1021/acs.jpcb.2c05496] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/22/2022]
Abstract
Structural waters in the S1 binding pocket of β-trypsin are critical for the stabilization of the complex of β-trypsin with its inhibitor bovine pancreatic trypsin inhibitor (BPTI). The inhibitor strength of BPTI can be modulated by replacing the critical lysine residue at the P1 position by non-natural amino acids. We study BPTI variants in which the critical Lys15 in BPTI has been replaced by α-aminobutyric acid (Abu) and its fluorinated derivatives monofluoroethylglycine (MfeGly), difluoroethylglycine (DfeGly), and trifluoroethylglycine (TfeGly). We investigate the hypothesis that additional water molecules in the binding pocket can form specific noncovalent interactions with the fluorinated side chains and thereby act as an extension of the inhibitors. We report potentials of mean force (PMF) of the unbinding process for all four complexes and enzyme activity inhibition assays. Additionally, we report the protein crystal structure of the Lys15MfeGly-BPTI-β-trypsin complex (pdb: 7PH1). Both experimental and computational data show a stepwise increase in inhibitor strength with increasing fluorination of the Abu side chain. The PMF additionally shows a minimum for the encounter complex and an intermediate state just before the bound state. In the bound state, the computational analysis of the structure and dynamics of the water molecules in the S1 pocket shows a highly dynamic network of water molecules that does not indicate a rigidification or stabilizing trend in regard to energetic properties that could explain the increase in inhibitor strength. The analysis of the energy and the entropy of the water molecules in the S1 binding pocket using grid inhomogeneous solvation theory confirms this result. Overall, fluorination systematically changes the binding affinity, but the effect cannot be explained by a persistent water network in the binding pocket. Other effects, such as the hydrophobicity of fluorinated amino acids and the stability of the encounter complex as well as the additional minimum in the potential of mean force in the bound state, likely influence the affinity more directly.
Collapse
Affiliation(s)
- Leon Wehrhan
- Department of Biology, Chemistry, and Pharmacy, Freie Universität Berlin, Institute of Chemistry and Biochemistry, Arnimallee 22, Berlin14195, Germany
| | - Jakob Leppkes
- Department of Biology, Chemistry, and Pharmacy, Freie Universität Berlin, Institute of Chemistry and Biochemistry, Arnimallee 20, Berlin14195, Germany
| | - Nicole Dimos
- Department of Biology, Chemistry, and Pharmacy, Freie Universität Berlin, Institute of Chemistry and Biochemistry, Takustr. 6, Berlin14195, Germany
| | - Bernhard Loll
- Department of Biology, Chemistry, and Pharmacy, Freie Universität Berlin, Institute of Chemistry and Biochemistry, Takustr. 6, Berlin14195, Germany
| | - Beate Koksch
- Department of Biology, Chemistry, and Pharmacy, Freie Universität Berlin, Institute of Chemistry and Biochemistry, Arnimallee 20, Berlin14195, Germany
| | - Bettina G Keller
- Department of Biology, Chemistry, and Pharmacy, Freie Universität Berlin, Institute of Chemistry and Biochemistry, Arnimallee 22, Berlin14195, Germany
| |
Collapse
|
49
|
Markov field models: Scaling molecular kinetics approaches to large molecular machines. Curr Opin Struct Biol 2022; 77:102458. [PMID: 36162297 DOI: 10.1016/j.sbi.2022.102458] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/25/2022] [Accepted: 08/05/2022] [Indexed: 12/14/2022]
Abstract
With recent advances in structural biology, including experimental techniques and deep learning-enabled high-precision structure predictions, molecular dynamics methods that scale up to large biomolecular systems are required. Current state-of-the-art approaches in molecular dynamics modeling focus on encoding global configurations of molecular systems as distinct states. This paradigm commands us to map out all possible structures and sample transitions between them, a task that becomes impossible for large-scale systems such as biomolecular complexes. To arrive at scalable molecular models, we suggest moving away from global state descriptions to a set of coupled models that each describe the dynamics of local domains or sites of the molecular system. We describe limitations in the current state-of-the-art global-state Markovian modeling approaches and then introduce Markov field models as an umbrella term that includes models from various scientific communities, including Independent Markov decomposition, Ising and Potts models, and (dynamic) graphical models, and evaluate their use for computational molecular biology. Finally, we give a few examples of early adoptions of these ideas for modeling molecular kinetics and thermodynamics.
Collapse
|
50
|
Ziada S, Diharce J, Raimbaud E, Aci-Sèche S, Ducrot P, Bonnet P. Estimation of Drug-Target Residence Time by Targeted Molecular Dynamics Simulations. J Chem Inf Model 2022; 62:5536-5549. [PMID: 36350238 DOI: 10.1021/acs.jcim.2c00852] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/11/2022]
Abstract
Drug-target residence time has emerged as a key selection factor in drug discovery since the binding duration of a drug molecule to its protein target can significantly impact its in vivo efficacy. The challenge in studying the residence time, in early drug discovery stages, lies in how to cost-effectively determine the residence time for the systematic assessment of compounds. Currently, there is still a lack of computational protocols to quickly estimate such a measure, particularly for large and flexible protein targets and drugs. Here, we report an efficient computational protocol, based on targeted molecular dynamics, to rank drug candidates by their residence time and to obtain insights into ligand-target dissociation mechanisms. The method was assessed on a dataset of 10 arylpyrazole inhibitors of CDK8, a large, flexible, and clinically important target, for which the experimental residence time of the inhibitors ranges from minutes to hours. The compounds were correctly ranked according to their estimated residence time scores compared to their experimental values. The analysis of protein-ligand interactions along the dissociation trajectories highlighted the favorable contribution of hydrophobic contacts to residence time and revealed key residues that strongly affect compound residence time.
Collapse
Affiliation(s)
- Sonia Ziada
- Institut de Chimie Organique et Analytique (ICOA), UMR CNRS-Université d'Orléans 7311, Université d'Orléans BP 6759, Orléans Cedex 245067, France
| | - Julien Diharce
- Institut de Chimie Organique et Analytique (ICOA), UMR CNRS-Université d'Orléans 7311, Université d'Orléans BP 6759, Orléans Cedex 245067, France
| | - Eric Raimbaud
- Institut de Recherches Servier, 125 Chemin de Ronde, Croissy-sur-Seine78290, France
| | - Samia Aci-Sèche
- Institut de Chimie Organique et Analytique (ICOA), UMR CNRS-Université d'Orléans 7311, Université d'Orléans BP 6759, Orléans Cedex 245067, France
| | - Pierre Ducrot
- Institut de Recherches Servier, 125 Chemin de Ronde, Croissy-sur-Seine78290, France
| | - Pascal Bonnet
- Institut de Chimie Organique et Analytique (ICOA), UMR CNRS-Université d'Orléans 7311, Université d'Orléans BP 6759, Orléans Cedex 245067, France
| |
Collapse
|