1
|
Parwana D, Gu J, Chen S, Bethel CR, Marshall E, Hujer AM, Bonomo RA, Haider S. The Structural Role of N170 in Substrate-Assisted Deacylation in KPC-2 β-Lactamase. Angew Chem Int Ed Engl 2024; 63:e202317315. [PMID: 38227422 DOI: 10.1002/anie.202317315] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/22/2023] [Revised: 01/16/2024] [Accepted: 01/16/2024] [Indexed: 01/17/2024]
Abstract
The amino acid substitutions in Klebsiella pneumoniae carbapenemase 2 (KPC-2) that have arisen in the clinic are observed to lead to the development of resistance to ceftazidime-avibactam, a preferred treatment for KPC bearing Gram-negative bacteria. Specific substitutions in the omega loop (R164-D179) result in changes in the structure and function of the enzyme, leading to alterations in substrate specificity, decreased stability, and more recently observed, increased resistance to ceftazidime/avibactam. Using accelerated rare-event sampling well-tempered metadynamics simulations, we explored in detail the structural role of R164 and D179 variants that are described to confer ceftazidime/avibactam resistance. The buried conformation of D179 substitutions produce a pronounced structural disorder in the omega loop - more than R164 mutants, where the crystallographic omega loop structure remains mostly intact. Our findings also reveal that the conformation of N170 plays an underappreciated role impacting drug binding and restricting deacylation. The results further support the hypothesis that KPC-2 D179 variants employ substrate-assisted catalysis for ceftazidime hydrolysis, involving the ring amine of the aminothiazole group to promote deacylation and catalytic turnover. Moreover, the shift in the WT conformation of N170 contributes to reduced deacylation and an altered spectrum of enzymatic activity.
Collapse
Affiliation(s)
| | - Jing Gu
- UCL School of Pharmacy, London, UK
| | | | - Christopher R Bethel
- Research Service, Louis Stokes Cleveland Department of Veterans Affairs Medical Center, Cleveland, OH, USA
| | - Emma Marshall
- Research Service, Louis Stokes Cleveland Department of Veterans Affairs Medical Center, Cleveland, OH, USA
| | - Andrea M Hujer
- Research Service, Louis Stokes Cleveland Department of Veterans Affairs Medical Center, Cleveland, OH, USA
- Department of Medicine, Case Western Reserve University School of Medicine, Cleveland, OH, USA
| | - Robert A Bonomo
- Research Service, Louis Stokes Cleveland Department of Veterans Affairs Medical Center, Cleveland, OH, USA
- Department of Medicine, Case Western Reserve University School of Medicine, Cleveland, OH, USA
- Clinician Scientist Investigator, Louis Stokes Cleveland Department of Veterans Affairs Medical Center, Cleveland, OH, USA
- Department of Molecular Biology and Microbiology, Pharmacology, Biochemistry, and Proteomics and Bioinformatics, Case Western Reserve University School of Medicine, Cleveland, OH, USA
- CWRU-Cleveland VAMC Center for Antimicrobial Resistance and Epidemiology (Case VA CARES), Cleveland, OH, USA
| | - Shozeb Haider
- UCL School of Pharmacy, London, UK
- UCL Centre for Advanced Research Computing, London, UK
| |
Collapse
|
2
|
Vijay A, Sreyas Adury VS, Mukherjee A. Targeting RdRp of SARS-CoV-2 with De Novo Molecule Generation. ACS Appl Bio Mater 2024; 7:609-616. [PMID: 37566736 DOI: 10.1021/acsabm.3c00339] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 08/13/2023]
Abstract
Viruses are known for their extremely high mutation rates, allowing them to evade both the human immune system and many forms of standard medicine. Despite this, the RNA dependent RNA polymerase (RdRp) of the RNA viruses has been largely conserved, and any significant mutation of this protein is unlikely. The recent COVID-19 pandemic presents a need for therapeutics. We have designed a de novo drug design algorithm that generates strong binding ligands from scratch, based on only the structure of the target protein's receptor. In this paper, we applied our method to target SARS-CoV-2 RdRp and generated several de novo molecules. We then chose some drug molecules based on the structural similarity to some of our strongest binding de novo molecules. Subsequently, we showed, using rigorous all-atom explicit-water free energy calculations in near-microsecond time scales using state-of-the-art well-tempered metadynamics simulations, that some of our de novo generated ligands bind more strongly to RdRp than the recent FDA approved drug remdesivir in its active form, remdesivir triphosphate (RTP). We elucidated the binding mechanism for some of the top binders and compared it with RTP. We believe that this work will be useful both by presenting lead structures for RdRp inhibition and by delivering key insights into the residues of the protein potentially involved in the binding/unbinding of these small molecule drugs, leading to more targeted studies in the future.
Collapse
Affiliation(s)
- Amal Vijay
- Department of Chemistry, Indian Institute of Science Education and Research, Pune 411008, India
| | | | - Arnab Mukherjee
- Department of Chemistry, Indian Institute of Science Education and Research, Pune 411008, India
| |
Collapse
|
3
|
Liu Z. Accelerating Kinetics with Time-Reversal Path Sampling. Molecules 2023; 28:8147. [PMID: 38138635 PMCID: PMC10745403 DOI: 10.3390/molecules28248147] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/18/2023] [Revised: 12/07/2023] [Accepted: 12/13/2023] [Indexed: 12/24/2023] Open
Abstract
In comparison to numerous enhanced sampling methods for equilibrium thermodynamics, accelerating simulations for kinetics and nonequilibrium statistics are relatively rare and less effective. Here, we derive a time-reversal path sampling (tRPS) method based on time reversibility to accelerate simulations for determining the transition rates between free-energy basins. It converts the difficult uphill path sampling into an easy downhill problem. This method is easy to implement, i.e., forward and backward shooting simulations with opposite initial velocities are conducted from random initial conformations within a transition-state region until they reach the basin minima, which are then assembled to give the distribution of transition paths efficiently. The effects of tRPS are demonstrated using a comparison with direct simulations of protein folding and unfolding, where tRPS is shown to give results consistent with direct simulations and increase the efficiency by up to five orders of magnitude. This approach is generally applicable to stochastic processes with microscopic reversibility, regardless of whether the variables are continuous or discrete.
Collapse
Affiliation(s)
- Zhirong Liu
- Beijing National Laboratory for Molecular Sciences (BNLMS), College of Chemistry and Molecular Engineering, Peking University, Beijing 100871, China
| |
Collapse
|
4
|
Bonati L, Polino D, Pizzolitto C, Biasi P, Eckert R, Reitmeier S, Schlögl R, Parrinello M. The role of dynamics in heterogeneous catalysis: Surface diffusivity and N 2 decomposition on Fe(111). Proc Natl Acad Sci U S A 2023; 120:e2313023120. [PMID: 38060558 PMCID: PMC10723053 DOI: 10.1073/pnas.2313023120] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/30/2023] [Accepted: 10/18/2023] [Indexed: 12/17/2023] Open
Abstract
Dynamics has long been recognized to play an important role in heterogeneous catalytic processes. However, until recently, it has been impossible to study their dynamical behavior at industry-relevant temperatures. Using a combination of machine learning potentials and advanced simulation techniques, we investigate the cleavage of the N[Formula: see text] triple bond on the Fe(111) surface. We find that at low temperatures our results agree with the well-established picture. However, if we increase the temperature to reach operando conditions, the surface undergoes a global dynamical change and the step structure of the Fe(111) surface is destabilized. The catalytic sites, traditionally associated with this surface, appear and disappear continuously. Our simulations illuminate the danger of extrapolating low-temperature results to operando conditions and indicate that the catalytic activity can only be inferred from calculations that take dynamics fully into account. More than that, they show that it is the transition to this highly fluctuating interfacial environment that drives the catalytic process.
Collapse
Affiliation(s)
- Luigi Bonati
- Atomistic Simulations, Italian Institute of Technology, Genova16152, Italy
| | - Daniela Polino
- Department of Innovative Technologies, University of Applied Sciences and Arts of Southern Switzerland, Lugano6962, Switzerland
| | - Cristina Pizzolitto
- Basic Research, Research and Development Division, Casale SA, Lugano6900, Switzerland
| | - Pierdomenico Biasi
- Basic Research, Research and Development Division, Casale SA, Lugano6900, Switzerland
| | - Rene Eckert
- BU Catalysts, R&D Syngas Applications, Clariant Produkte (Deutschland) GmbH, Munich83052, Germany
| | - Stephan Reitmeier
- BU Catalysts, R&D Syngas Applications, Clariant Produkte (Deutschland) GmbH, Munich83052, Germany
| | - Robert Schlögl
- Department of Inorganic Chemistry, Fritz-Haber Institute of the Max-Planck-Society, Berlin14195, Germany
| | - Michele Parrinello
- Atomistic Simulations, Italian Institute of Technology, Genova16152, Italy
| |
Collapse
|
5
|
Zhu J, Li Z, Tong H, Lu Z, Zhang N, Wei T, Chen HF. Phanto-IDP: compact model for precise intrinsically disordered protein backbone generation and enhanced sampling. Brief Bioinform 2023; 25:bbad429. [PMID: 38018910 PMCID: PMC10783862 DOI: 10.1093/bib/bbad429] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/14/2023] [Revised: 09/21/2023] [Accepted: 11/05/2023] [Indexed: 11/30/2023] Open
Abstract
The biological function of proteins is determined not only by their static structures but also by the dynamic properties of their conformational ensembles. Numerous high-accuracy static structure prediction tools have been recently developed based on deep learning; however, there remains a lack of efficient and accurate methods for exploring protein dynamic conformations. Traditionally, studies concerning protein dynamics have relied on molecular dynamics (MD) simulations, which incur significant computational costs for all-atom precision and struggle to adequately sample conformational spaces with high energy barriers. To overcome these limitations, various enhanced sampling techniques have been developed to accelerate sampling in MD. Traditional enhanced sampling approaches like replica exchange molecular dynamics (REMD) and frontier expansion sampling (FEXS) often follow the MD simulation approach and still cost a lot of computational resources and time. Variational autoencoders (VAEs), as a classic deep generative model, are not restricted by potential energy landscapes and can explore conformational spaces more efficiently than traditional methods. However, VAEs often face challenges in generating reasonable conformations for complex proteins, especially intrinsically disordered proteins (IDPs), which limits their application as an enhanced sampling method. In this study, we presented a novel deep learning model (named Phanto-IDP) that utilizes a graph-based encoder to extract protein features and a transformer-based decoder combined with variational sampling to generate highly accurate protein backbones. Ten IDPs and four structured proteins were used to evaluate the sampling ability of Phanto-IDP. The results demonstrate that Phanto-IDP has high fidelity and diversity in the generated conformation ensembles, making it a suitable tool for enhancing the efficiency of MD simulation, generating broader protein conformational space and a continuous protein transition path.
Collapse
Affiliation(s)
- Junjie Zhu
- State Key Laboratory of Microbial Metabolism, Joint International Research Laboratory of Metabolic & Developmental Sciences, Department of Bioinformatics and Biostatistics, National Experimental Teaching Center for Life Sciences and Biotechnology, School of Life Sciences and Biotechnology, Shanghai Jiao Tong University, Shanghai, 200240, China
| | - Zhengxin Li
- State Key Laboratory of Microbial Metabolism, Joint International Research Laboratory of Metabolic & Developmental Sciences, Department of Bioinformatics and Biostatistics, National Experimental Teaching Center for Life Sciences and Biotechnology, School of Life Sciences and Biotechnology, Shanghai Jiao Tong University, Shanghai, 200240, China
| | - Haowei Tong
- State Key Laboratory of Microbial Metabolism, Joint International Research Laboratory of Metabolic & Developmental Sciences, Department of Bioinformatics and Biostatistics, National Experimental Teaching Center for Life Sciences and Biotechnology, School of Life Sciences and Biotechnology, Shanghai Jiao Tong University, Shanghai, 200240, China
| | - Zhouyu Lu
- State Key Laboratory of Microbial Metabolism, Joint International Research Laboratory of Metabolic & Developmental Sciences, Department of Bioinformatics and Biostatistics, National Experimental Teaching Center for Life Sciences and Biotechnology, School of Life Sciences and Biotechnology, Shanghai Jiao Tong University, Shanghai, 200240, China
| | - Ningjie Zhang
- State Key Laboratory of Microbial Metabolism, Joint International Research Laboratory of Metabolic & Developmental Sciences, Department of Bioinformatics and Biostatistics, National Experimental Teaching Center for Life Sciences and Biotechnology, School of Life Sciences and Biotechnology, Shanghai Jiao Tong University, Shanghai, 200240, China
| | - Ting Wei
- State Key Laboratory of Microbial Metabolism, Joint International Research Laboratory of Metabolic & Developmental Sciences, Department of Bioinformatics and Biostatistics, National Experimental Teaching Center for Life Sciences and Biotechnology, School of Life Sciences and Biotechnology, Shanghai Jiao Tong University, Shanghai, 200240, China
| | - Hai-Feng Chen
- State Key Laboratory of Microbial Metabolism, Joint International Research Laboratory of Metabolic & Developmental Sciences, Department of Bioinformatics and Biostatistics, National Experimental Teaching Center for Life Sciences and Biotechnology, School of Life Sciences and Biotechnology, Shanghai Jiao Tong University, Shanghai, 200240, China
| |
Collapse
|
6
|
Abstract
This Perspective provides a contextual explanation of the current state-of-the-art alchemical free energy methods and their role in drug discovery as well as highlights select emerging technologies. The narrative attempts to answer basic questions about what goes on "under the hood" in free energy simulations and provide general guidelines for how to run simulations and analyze the results. It is the hope that this work will provide a valuable introduction to students and scientists in the field.
Collapse
Affiliation(s)
- Darrin M. York
- Laboratory for Biomolecular
Simulation Research, Institute for Quantitative Biomedicine, and Department
of Chemistry and Chemical Biology, Rutgers
University, Piscataway, New Jersey 08854, United States
| |
Collapse
|
7
|
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] [What about the content of this article? (0)] [Affiliation(s)] [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
|
8
|
Calegari Andrade M, Car R, Selloni A. Probing the self-ionization of liquid water with ab initio deep potential molecular dynamics. Proc Natl Acad Sci U S A 2023; 120:e2302468120. [PMID: 37931100 PMCID: PMC10655216 DOI: 10.1073/pnas.2302468120] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/12/2023] [Accepted: 09/29/2023] [Indexed: 11/08/2023] Open
Abstract
The chemical equilibrium between self-ionized and molecular water dictates the acid-base chemistry in aqueous solutions, yet understanding the microscopic mechanisms of water self-ionization remains experimentally and computationally challenging. Herein, Density Functional Theory (DFT)-based deep neural network (DNN) potentials are combined with enhanced sampling techniques and a global acid-base collective variable to perform extensive atomistic simulations of water self-ionization for model systems of increasing size. The explicit inclusion of long-range electrostatic interactions in the DNN potential is found to be crucial to accurately reproduce the DFT free energy profile of solvated water ion pairs in small (64 and 128 H2O) cells. The reversible work to separate the hydroxide and hydronium to a distance [Formula: see text] is found to converge for simulation cells containing more than 500 H2O, and a distance of [Formula: see text] 8 Å is the threshold beyond which the work to further separate the two ions becomes approximately zero. The slow convergence of the potential of mean force with system size is related to a restructuring of water and an increase of the local order around the water ions. Calculation of the dissociation equilibrium constant illustrates the key role of long-range electrostatics and entropic effects in the water autoionization process.
Collapse
Affiliation(s)
- Marcos Calegari Andrade
- Chemistry Department, Princeton University, Princeton, NJ08544
- Quantum Simulations Group, Materials Science Division, Lawrence Livermore National Laboratory, Livermore, CA94550
| | - Roberto Car
- Chemistry Department, Princeton University, Princeton, NJ08544
| | | |
Collapse
|
9
|
Kienlein M, Zacharias M, Reif MM. Efficient and accurate calculation of proline cis/trans isomerization free energies from Hamiltonian replica exchange molecular dynamics simulations. Structure 2023; 31:1473-1484.e6. [PMID: 37657438 DOI: 10.1016/j.str.2023.08.008] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/18/2023] [Revised: 04/17/2023] [Accepted: 08/07/2023] [Indexed: 09/03/2023]
Abstract
Proline cis/trans isomerization plays an important role in many biological processes but occurs on time scales not accessible to brute-force molecular dynamics (MD) simulations. We have designed a new Hamiltonian replica exchange scheme, ω-bias potential replica exchange molecular dynamics (ωBP-REMD), to efficiently and accurately calculate proline cis/trans isomerization free energies. ωBP-REMD is applied to various proline-containing tripeptides and a biologically important proline residue in the N2-domain of the gene-3-protein of phage fd in the wildtype and mutant variants of the protein. Excellent cis/trans transition rates are obtained. Reweighting of the sampled probability distribution along the peptide bond dihedral angle allows construction of the corresponding free-energy profile and calculation of the cis/trans isomerization free energy with high statistical precision. Very good agreement with experimental data is obtained. ωBP-REMD outperforms standard umbrella sampling in terms of convergence and agreement with experiment and strongly reduces perturbation of the local structure near the proline residue.
Collapse
Affiliation(s)
- Maximilian Kienlein
- Center for Functional Protein Assemblies (CPA), Physics Department, Chair of Theoretical Biophysics (T38), Technical University of Munich, Ernst-Otto-Fischer-Str. 8, 85748 Garching, Germany
| | - Martin Zacharias
- Center for Functional Protein Assemblies (CPA), Physics Department, Chair of Theoretical Biophysics (T38), Technical University of Munich, Ernst-Otto-Fischer-Str. 8, 85748 Garching, Germany
| | - Maria M Reif
- Center for Functional Protein Assemblies (CPA), Physics Department, Chair of Theoretical Biophysics (T38), Technical University of Munich, Ernst-Otto-Fischer-Str. 8, 85748 Garching, Germany.
| |
Collapse
|
10
|
Fortino M, Schifino G, Pietropaolo A. Simulation workflows to predict the circular dichroism and circularly polarized luminescence of chiral materials. Chirality 2023; 35:673-680. [PMID: 36896846 DOI: 10.1002/chir.23546] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/18/2022] [Revised: 02/02/2023] [Accepted: 02/06/2023] [Indexed: 03/11/2023]
Abstract
Chiral materials are attracting considerable interest in various fields in view of their unique properties and optical activity. Indeed, the peculiar features of chiral materials to absorb and emit circularly polarized light enable their use in an extensive range of applications. Motivated by the interest in boosting the development of chiral materials characterized by enhanced chiroptical properties such as circular dichroism (CD) and circular polarized luminescence (CPL), we herein illustrate in this tutorial how theoretical simulations can be used for the predictions and interpretations of chiroptical data and for the identification of chiral geometries. We are focusing on computational frameworks that can be used to investigate the theoretical aspects of chiral materials' photophysical and conformational characteristics. We will then illustrate ab initio methods based on density functional theory (DFT) and its time-dependent extension (TD-DFT) to simulate CD and CPL signals, and we will exemplify a variety of enhanced sampling techniques useful for an adequate sampling of the configurational space for chiral systems.
Collapse
Affiliation(s)
- Mariagrazia Fortino
- Dipartimento di Scienze della Salute, Università di Catanzaro, Catanzaro, Italy
| | - Gioacchino Schifino
- Dipartimento di Scienze della Salute, Università di Catanzaro, Catanzaro, Italy
| | - Adriana Pietropaolo
- Dipartimento di Scienze della Salute, Università di Catanzaro, Catanzaro, Italy
| |
Collapse
|
11
|
Math BA, Waibl F, Lamp LM, Fernández‐Quintero ML, Liedl KR. Cross-linking disulfide bonds govern solution structures of diabodies. Proteins 2023; 91:1316-1328. [PMID: 37376973 PMCID: PMC10952579 DOI: 10.1002/prot.26509] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/12/2022] [Accepted: 04/19/2023] [Indexed: 06/29/2023]
Abstract
In the last years, antibodies have emerged as a promising new class of therapeutics, due to their combination of high specificity with long serum half-life and low risk of side-effects. Diabodies are a popular novel antibody format, consisting of two Fv domains connected with short linkers. Like IgG antibodies, they simultaneously bind two target proteins. However, they offer altered properties, given their smaller size and higher rigidity. In this study, we conducted the-to our knowledge-first molecular dynamics (MD) simulations of diabodies and find a surprisingly high conformational flexibility in the relative orientation of the two Fv domains. We observe rigidifying effects through the introduction of disulfide bonds in the Fv -Fv interface and characterize the effect of different disulfide bond locations on the conformation. Additionally, we compare VH -VL orientations and paratope dynamics between diabodies and an antigen binding fragment (Fab) of the same sequence. We find mostly consistent structures and dynamics, indicating similar antigen binding properties. The most significant differences can be found within the CDR-H2 loop dynamics. Of all CDR loops, the CDR-H2 is located closest to the artificial Fv -Fv interface. All examined diabodies show similar VH -VL orientations, Fv -Fv packing and CDR loop conformations. However, the variant with a P14C-K64C disulfide bond differs most from the Fab in our measures, including the CDR-H3 loop conformational ensemble. This suggests altered antigen binding properties and underlines the need for careful validation of the disulfide bond locations in diabodies.
Collapse
Affiliation(s)
- Barbara A. Math
- Institute of General, Inorganic and Theoretical Chemistry, and Center for Molecular Biosciences Innsbruck (CMBI)University of InnsbruckInnsbruckAustria
| | - Franz Waibl
- Institute of General, Inorganic and Theoretical Chemistry, and Center for Molecular Biosciences Innsbruck (CMBI)University of InnsbruckInnsbruckAustria
| | - Leonida M. Lamp
- Institute of General, Inorganic and Theoretical Chemistry, and Center for Molecular Biosciences Innsbruck (CMBI)University of InnsbruckInnsbruckAustria
| | - Monica L. Fernández‐Quintero
- Institute of General, Inorganic and Theoretical Chemistry, and Center for Molecular Biosciences Innsbruck (CMBI)University of InnsbruckInnsbruckAustria
| | - Klaus R. Liedl
- Institute of General, Inorganic and Theoretical Chemistry, and Center for Molecular Biosciences Innsbruck (CMBI)University of InnsbruckInnsbruckAustria
| |
Collapse
|
12
|
Van Speybroeck V, Bocus M, Cnudde P, Vanduyfhuys L. Operando Modeling of Zeolite-Catalyzed Reactions Using First-Principles Molecular Dynamics Simulations. ACS Catal 2023; 13:11455-11493. [PMID: 37671178 PMCID: PMC10476167 DOI: 10.1021/acscatal.3c01945] [Citation(s) in RCA: 5] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/30/2023] [Revised: 07/27/2023] [Indexed: 09/07/2023]
Abstract
Within this Perspective, we critically reflect on the role of first-principles molecular dynamics (MD) simulations in unraveling the catalytic function within zeolites under operating conditions. First-principles MD simulations refer to methods where the dynamics of the nuclei is followed in time by integrating the Newtonian equations of motion on a potential energy surface that is determined by solving the quantum-mechanical many-body problem for the electrons. Catalytic solids used in industrial applications show an intriguing high degree of complexity, with phenomena taking place at a broad range of length and time scales. Additionally, the state and function of a catalyst critically depend on the operating conditions, such as temperature, moisture, presence of water, etc. Herein we show by means of a series of exemplary cases how first-principles MD simulations are instrumental to unravel the catalyst complexity at the molecular scale. Examples show how the nature of reactive species at higher catalytic temperatures may drastically change compared to species at lower temperatures and how the nature of active sites may dynamically change upon exposure to water. To simulate rare events, first-principles MD simulations need to be used in combination with enhanced sampling techniques to efficiently sample low-probability regions of phase space. Using these techniques, it is shown how competitive pathways at operating conditions can be discovered and how broad transition state regions can be explored. Interestingly, such simulations can also be used to study hindered diffusion under operating conditions. The cases shown clearly illustrate how first-principles MD simulations reveal insights into the catalytic function at operating conditions, which could not be discovered using static or local approaches where only a few points are considered on the potential energy surface (PES). Despite these advantages, some major hurdles still exist to fully integrate first-principles MD methods in a standard computational catalytic workflow or to use the output of MD simulations as input for multiple length/time scale methods that aim to bridge to the reactor scale. First of all, methods are needed that allow us to evaluate the interatomic forces with quantum-mechanical accuracy, albeit at a much lower computational cost compared to currently used density functional theory (DFT) methods. The use of DFT limits the currently attainable length/time scales to hundreds of picoseconds and a few nanometers, which are much smaller than realistic catalyst particle dimensions and time scales encountered in the catalysis process. One solution could be to construct machine learning potentials (MLPs), where a numerical potential is derived from underlying quantum-mechanical data, which could be used in subsequent MD simulations. As such, much longer length and time scales could be reached; however, quite some research is still necessary to construct MLPs for the complex systems encountered in industrially used catalysts. Second, most currently used enhanced sampling techniques in catalysis make use of collective variables (CVs), which are mostly determined based on chemical intuition. To explore complex reactive networks with MD simulations, methods are needed that allow the automatic discovery of CVs or methods that do not rely on a priori definition of CVs. Recently, various data-driven methods have been proposed, which could be explored for complex catalytic systems. Lastly, first-principles MD methods are currently mostly used to investigate local reactive events. We hope that with the rise of data-driven methods and more efficient methods to describe the PES, first-principles MD methods will in the future also be able to describe longer length/time scale processes in catalysis. This might lead to a consistent dynamic description of all steps-diffusion, adsorption, and reaction-as they take place at the catalyst particle level.
Collapse
Affiliation(s)
| | - Massimo Bocus
- Center for Molecular Modeling, Ghent University, Technologiepark 46, 9052 Zwijnaarde, Belgium
| | - Pieter Cnudde
- Center for Molecular Modeling, Ghent University, Technologiepark 46, 9052 Zwijnaarde, Belgium
| | - Louis Vanduyfhuys
- Center for Molecular Modeling, Ghent University, Technologiepark 46, 9052 Zwijnaarde, Belgium
| |
Collapse
|
13
|
Wu M, Liao J, Shu Z, Chen C. Enhanced sampling in explicit solvent by deep learning module in FSATOOL. J Comput Chem 2023. [PMID: 37191088 DOI: 10.1002/jcc.27132] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/06/2023] [Revised: 04/21/2023] [Accepted: 04/27/2023] [Indexed: 05/17/2023]
Abstract
FSATOOL is an integrated molecular simulation and data analysis program. Its old molecular dynamics engine only supports simulations in vacuum or implicit solvent. In this work, we implement the well-known smooth particle mesh Ewald method for simulations in explicit solvent. The new developed engine is runnable on both CPU and GPU. All the existed analysis modules in the program are compatible with the new engine. Moreover, we also build a complete deep learning module in FSATOOL. Based on the module, we further implement two useful trajectory analysis methods: state-free reversible VAMPnets and time-lagged autoencoder. They are good at searching the collective variables related to the conformational transitions of biomolecules. In FSATOOL, these collective variables can be further used to construct a bias potential for the enhanced sampling purpose. We introduce the implementation details of the methods and present their actual performances in FSATOOL by a few enhanced sampling simulations.
Collapse
Affiliation(s)
- Mincong Wu
- Biomolecular Physics and Modeling Group, School of Physics, Huazhong University of Science and Technology, Wuhan, Hubei, China
| | - Jun Liao
- Biomolecular Physics and Modeling Group, School of Physics, Huazhong University of Science and Technology, Wuhan, Hubei, China
| | - Zirui Shu
- Biomolecular Physics and Modeling Group, School of Physics, Huazhong University of Science and Technology, Wuhan, Hubei, China
| | - Changjun Chen
- Biomolecular Physics and Modeling Group, School of Physics, Huazhong University of Science and Technology, Wuhan, Hubei, China
| |
Collapse
|
14
|
Zheng LE, Barethiya S, Nordquist E, Chen J. Machine Learning Generation of Dynamic Protein Conformational Ensembles. Molecules 2023; 28:4047. [PMID: 37241789 PMCID: PMC10220786 DOI: 10.3390/molecules28104047] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/10/2023] [Revised: 05/04/2023] [Accepted: 05/09/2023] [Indexed: 05/28/2023] Open
Abstract
Machine learning has achieved remarkable success across a broad range of scientific and engineering disciplines, particularly its use for predicting native protein structures from sequence information alone. However, biomolecules are inherently dynamic, and there is a pressing need for accurate predictions of dynamic structural ensembles across multiple functional levels. These problems range from the relatively well-defined task of predicting conformational dynamics around the native state of a protein, which traditional molecular dynamics (MD) simulations are particularly adept at handling, to generating large-scale conformational transitions connecting distinct functional states of structured proteins or numerous marginally stable states within the dynamic ensembles of intrinsically disordered proteins. Machine learning has been increasingly applied to learn low-dimensional representations of protein conformational spaces, which can then be used to drive additional MD sampling or directly generate novel conformations. These methods promise to greatly reduce the computational cost of generating dynamic protein ensembles, compared to traditional MD simulations. In this review, we examine recent progress in machine learning approaches towards generative modeling of dynamic protein ensembles and emphasize the crucial importance of integrating advances in machine learning, structural data, and physical principles to achieve these ambitious goals.
Collapse
Affiliation(s)
- Li-E Zheng
- Department of Gynecology, The First Affiliated Hospital of Fujian Medical University, Fuzhou 350005, China;
| | - Shrishti Barethiya
- Department of Chemistry, University of Massachusetts Amherst, Amherst, MA 01003, USA; (S.B.); (E.N.)
| | - Erik Nordquist
- Department of Chemistry, University of Massachusetts Amherst, Amherst, MA 01003, USA; (S.B.); (E.N.)
| | - Jianhan Chen
- Department of Chemistry, University of Massachusetts Amherst, Amherst, MA 01003, USA; (S.B.); (E.N.)
| |
Collapse
|
15
|
Tu G, Xu B, Luo D, Liu J, Liu Z, Chen G, Xue W. Multi-state Model-Based Identification of Cryptic Allosteric Sites on Human Serotonin Transporter. ACS Chem Neurosci 2023; 14:1686-1694. [PMID: 37067527 DOI: 10.1021/acschemneuro.3c00155] [Citation(s) in RCA: 5] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/18/2023] Open
Abstract
Serotonin transporter (SERT) plays a fundamental role in taking the synaptic cleft serotonin back to the presynaptic neuron. The discovery of allosteric SERT modulators represents the next-generation medication for psychiatric disorders such as depression. Here, based on the cryo-EM structures of ibogaine in complex with SERT in distinct conformations, the multiple functional structures of the transporter bound to serotonin, including outward-open (OOholo), outward-occluded (OCholo), and inward-open (IOholo and IOholo'), were carefully characterized by induced-fit docking Gaussian-accelerated molecular dynamics (IFD-GaMD) simulation and the free-energy landscape analysis. Further MM/GBSA binding free energy, per-residue contribution, and molecular interaction fingerprint calculations revealed the interaction variations of serotonin with SERT in functional structures, which confirmed the allostery of SERT during serotonin reuptake. Moreover, five unique cryptic allosteric sites, which are druggable and capable of targeting by small molecules, were identified on the characterized multistate structures. These results provide structural and energetic information for the molecular mechanism of serotonin reuptake and will provide opportunities for the development of novel therapeutics based on the identified new allosteric sites on SERT.
Collapse
Affiliation(s)
- Gao Tu
- Chongqing Key Laboratory of Natural Product Synthesis and Drug Research, School of Pharmaceutical Sciences, Chongqing University, Chongqing 401331, China
| | - Binbin Xu
- Chongqing Key Laboratory of Natural Product Synthesis and Drug Research, School of Pharmaceutical Sciences, Chongqing University, Chongqing 401331, China
| | - Ding Luo
- Chongqing Key Laboratory of Natural Product Synthesis and Drug Research, School of Pharmaceutical Sciences, Chongqing University, Chongqing 401331, China
| | - Jin Liu
- Chongqing Key Laboratory of Natural Product Synthesis and Drug Research, School of Pharmaceutical Sciences, Chongqing University, Chongqing 401331, China
| | - Zerong Liu
- Central Nervous System Drug Key Laboratory of Sichuan Province, Sichuan Credit Pharmaceutical CO., Ltd., Luzhou 646000, China
| | - Gang Chen
- Central Nervous System Drug Key Laboratory of Sichuan Province, Sichuan Credit Pharmaceutical CO., Ltd., Luzhou 646000, China
| | - Weiwei Xue
- Chongqing Key Laboratory of Natural Product Synthesis and Drug Research, School of Pharmaceutical Sciences, Chongqing University, Chongqing 401331, China
| |
Collapse
|
16
|
Staritzbichler R, Ristic N, Stapke T, Hildebrand P. SmoothT - a server constructing low energy pathways from conformational ensembles for interactive visualization and enhanced sampling. Bioinformatics 2023; 39:7108772. [PMID: 37018142 PMCID: PMC10121331 DOI: 10.1093/bioinformatics/btad176] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Subscribe] [Scholar Register] [Received: 10/20/2022] [Revised: 01/21/2023] [Accepted: 03/10/2023] [Indexed: 04/06/2023]
Abstract
MOTIVATION The SmoothT software and webservice offers the construction of pathways from an ensemble of conformations. The user provides an archive of molecule conformations in PDB format, from which a starting and a final conformation need to be selected. The individual PDB files need to contain an energy value or score, estimating the quality of the respective conformation. Additionally, the user has to provide a RMSD cutoff, below which conformations are considered neighboring. From this SmoothT constructs a graph that connects similar conformations. RESULTS SmoothT returns the energetically most favorable pathway within in this graph. This pathway is directly displayed as interactive animation using the NGL viewer. Simultaneously, the energy along the pathway is plotted, highlighting the conformation that is currently displayed in the 3D window. AVAILABILITY AND IMPLEMENTATION SmoothT is available as webservice at: http://proteinformatics.org/smoothT. Examples, a tutorial and FAQs can be found there. Ensembles up to 2 GB (compressed) can be uploaded. Results will be stored for 5 days. The server is completely free and requires no registration. The C ++ source code is available at: https://github.com/starbeachlab/smoothT.
Collapse
Affiliation(s)
- René Staritzbichler
- Institute of Medical Physics and Biophysics, University of Leipzig, 04107, Leipzig, Germany
| | - Nikola Ristic
- Institute of Medical Physics and Biophysics, University of Leipzig, 04107, Leipzig, Germany
| | - Tülin Stapke
- Institute of Medical Physics and Biophysics, University of Leipzig, 04107, Leipzig, Germany
| | - Peter Hildebrand
- Institute of Medical Physics and Biophysics, University of Leipzig, 04107, Leipzig, Germany
- Institute of Medical Physics and Biophysics, Charité Universitätsmedizin Berlin, corporate member of Freie Universität Berlin und Humboldt- Universität zu Berlin, Berlin, Germany
- Berlin Institute of Health, 10178, Berlin, Germany
| |
Collapse
|
17
|
Gu S, Shen C, Yu J, Zhao H, Liu H, Liu L, Sheng R, Xu L, Wang Z, Hou T, Kang Y. Can molecular dynamics simulations improve predictions of protein-ligand binding affinity with machine learning? Brief Bioinform 2023; 24:6995375. [PMID: 36681903 DOI: 10.1093/bib/bbad008] [Citation(s) in RCA: 5] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/08/2022] [Revised: 12/04/2022] [Accepted: 12/30/2023] [Indexed: 01/23/2023] Open
Abstract
Binding affinity prediction largely determines the discovery efficiency of lead compounds in drug discovery. Recently, machine learning (ML)-based approaches have attracted much attention in hopes of enhancing the predictive performance of traditional physics-based approaches. In this study, we evaluated the impact of structural dynamic information on the binding affinity prediction by comparing the models trained on different dimensional descriptors, using three targets (i.e. JAK1, TAF1-BD2 and DDR1) and their corresponding ligands as the examples. Here, 2D descriptors are traditional ECFP4 fingerprints, 3D descriptors are the energy terms of the Smina and NNscore scoring functions and 4D descriptors contain the structural dynamic information derived from the trajectories based on molecular dynamics (MD) simulations. We systematically investigate the MD-refined binding affinity prediction performance of three classical ML algorithms (i.e. RF, SVR and XGB) as well as two common virtual screening methods, namely Glide docking and MM/PBSA. The outcomes of the ML models built using various dimensional descriptors and their combinations reveal that the MD refinement with the optimized protocol can improve the predictive performance on the TAF1-BD2 target with considerable structural flexibility, but not for the less flexible JAK1 and DDR1 targets, when taking docking poses as the initial structure instead of the crystal structures. The results highlight the importance of the initial structures to the final performance of the model through conformational analysis on the three targets with different flexibility.
Collapse
Affiliation(s)
- Shukai Gu
- Innovation Institute for Artificial Intelligence in Medicine of Zhejiang University, College of Pharmaceutical Sciences, Zhejiang University, Hangzhou 310058, Zhejiang, China
| | - Chao Shen
- Innovation Institute for Artificial Intelligence in Medicine of Zhejiang University, College of Pharmaceutical Sciences, Zhejiang University, Hangzhou 310058, Zhejiang, China
| | - Jiahui Yu
- Innovation Institute for Artificial Intelligence in Medicine of Zhejiang University, College of Pharmaceutical Sciences, Zhejiang University, Hangzhou 310058, Zhejiang, China
| | - Hong Zhao
- Innovation Institute for Artificial Intelligence in Medicine of Zhejiang University, College of Pharmaceutical Sciences, Zhejiang University, Hangzhou 310058, Zhejiang, China
| | - Huanxiang Liu
- Faculty of Applied Science, Macao Polytechnic University, Macao, SAR, China
| | - Liwei Liu
- Advanced Computing and Storage Laboratory, Central Research Institute, 2012 Laboratories, Huawei Technologies Co., Ltd., Shenzhen 518129, Guangdong, China
| | - Rong Sheng
- Health Technology Development Dept, Huawei Device Co., Ltd., Dongguan 523808, Guangdong, China
| | - Lei Xu
- Institute of Bioinformatics and Medical Engineering, School of Electrical and Information Engineering, Jiangsu University of Technology, Changzhou 213001, China
| | - Zhe Wang
- Innovation Institute for Artificial Intelligence in Medicine of Zhejiang University, College of Pharmaceutical Sciences, Zhejiang University, Hangzhou 310058, Zhejiang, China
| | - Tingjun Hou
- Innovation Institute for Artificial Intelligence in Medicine of Zhejiang University, College of Pharmaceutical Sciences, Zhejiang University, Hangzhou 310058, Zhejiang, China
| | - Yu Kang
- Innovation Institute for Artificial Intelligence in Medicine of Zhejiang University, College of Pharmaceutical Sciences, Zhejiang University, Hangzhou 310058, Zhejiang, China
| |
Collapse
|
18
|
Seidler CA, Kokot J, Fernández-Quintero ML, Liedl KR. Structural Characterization of Nanobodies during Germline Maturation. Biomolecules 2023; 13:380. [PMID: 36830754 PMCID: PMC9953242 DOI: 10.3390/biom13020380] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/15/2022] [Revised: 01/27/2023] [Accepted: 02/14/2023] [Indexed: 02/22/2023] Open
Abstract
Camelid heavy-chain antibody variable domains (VHH), nanobodies, are the smallest-known functional antibody fragments with high therapeutic potential. In this study, we investigate a VHH binding to hen egg-white lysozyme (HEL). We structurally and dynamically characterized the conformational diversity of four VHH variants to elucidate the antigen-binding process. For two of these antibodies, not only are the dissociation constants known, but also the experimentally determined crystal structures of the VHH in complex with HEL are available. We performed well-tempered metadynamics simulations in combination with molecular dynamics simulations to capture a broad conformational space and to reconstruct the thermodynamics and kinetics of conformational transitions in the antigen-binding site, the paratope. By kinetically characterizing the loop movements of the paratope, we found that, with an increase in affinity, the state populations shift towards the binding competent conformation. The contacts contributing to antigen binding, and those who contribute to the overall stability, show a clear trend towards less variable but more intense contacts. Additionally, these investigated nanobodies clearly follow the conformational selection paradigm, as the binding competent conformation pre-exists within the structural ensembles without the presence of the antigen.
Collapse
Affiliation(s)
| | | | - Monica L. Fernández-Quintero
- Department of General, Inorganic and Theoretical Chemistry, Center for Molecular Biosciences Innsbruck (CMBI), University of Innsbruck, 6020 Innsbruck, Austria
| | - Klaus R. Liedl
- Department of General, Inorganic and Theoretical Chemistry, Center for Molecular Biosciences Innsbruck (CMBI), University of Innsbruck, 6020 Innsbruck, Austria
| |
Collapse
|
19
|
Zou Z, Beyerle ER, Tsai ST, Tiwary P. Driving and characterizing nucleation of urea and glycine polymorphs in water. Proc Natl Acad Sci U S A 2023; 120:e2216099120. [PMID: 36757888 PMCID: PMC9963467 DOI: 10.1073/pnas.2216099120] [Citation(s) in RCA: 8] [Impact Index Per Article: 8.0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/20/2022] [Accepted: 01/17/2023] [Indexed: 02/10/2023] Open
Abstract
Crystal nucleation is relevant across the domains of fundamental and applied sciences. However, in many cases, its mechanism remains unclear due to a lack of temporal or spatial resolution. To gain insights into the molecular details of nucleation, some form of molecular dynamics simulations is typically performed; these simulations, in turn, are limited by their ability to run long enough to sample the nucleation event thoroughly. To overcome the timescale limits in typical molecular dynamics simulations in a manner free of prior human bias, here, we employ the machine learning-augmented molecular dynamics framework "reweighted autoencoded variational Bayes for enhanced sampling (RAVE)." We study two molecular systems-urea and glycine-in explicit all-atom water, due to their enrichment in polymorphic structures and common utility in commercial applications. From our simulations, we observe multiple back-and-forth nucleation events of different polymorphs from homogeneous solution; from these trajectories, we calculate the relative ranking of finite-sized polymorph crystals embedded in solution, in terms of the free-energy difference between the finite-sized crystal polymorph and the original solution state. We further observe that the obtained reaction coordinates and transitions are highly nonclassical.
Collapse
Affiliation(s)
- Ziyue Zou
- Department of Chemistry and Biochemistry, University of Maryland, College Park, MD20742
| | - Eric R. Beyerle
- Institute for Physical Science and Technology, University of Maryland, College Park, MD20742
| | - Sun-Ting Tsai
- Department of Physics, University of Maryland, College Park, MD20742
| | - Pratyush Tiwary
- Department of Chemistry and Biochemistry, University of Maryland, College Park, MD20742
- Institute for Physical Science and Technology, University of Maryland, College Park, MD20742
| |
Collapse
|
20
|
Sohraby F, Nunes-Alves A. Advances in computational methods for ligand binding kinetics. Trends Biochem Sci 2022; 48:437-449. [PMID: 36566088 DOI: 10.1016/j.tibs.2022.11.003] [Citation(s) in RCA: 6] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/27/2022] [Revised: 11/16/2022] [Accepted: 11/29/2022] [Indexed: 12/24/2022]
Abstract
Binding kinetic parameters can be correlated with drug efficacy, which in recent years led to the development of various computational methods for predicting binding kinetic rates and gaining insight into protein-drug binding paths and mechanisms. In this review, we introduce and compare computational methods recently developed and applied to two systems, trypsin-benzamidine and kinase-inhibitor complexes. Methods involving enhanced sampling in molecular dynamics simulations or machine learning can be used not only to predict kinetic rates, but also to reveal factors modulating the duration of residence times, selectivity, and drug resistance to mutations. Methods which require less computational time to make predictions are highlighted, and suggestions to reduce the error of computed kinetic rates are presented.
Collapse
Affiliation(s)
- Farzin Sohraby
- Institute of Chemistry, Technische Universität Berlin, 10623 Berlin, Germany
| | - Ariane Nunes-Alves
- Institute of Chemistry, Technische Universität Berlin, 10623 Berlin, Germany.
| |
Collapse
|
21
|
Kumawat A, Namsani S, Pramanik D, Roy S, Singh JK. Integrated docking and enhanced sampling-based selection of repurposing drugs for SARS-CoV-2 by targeting host dependent factors. J Biomol Struct Dyn 2022; 40:9897-9908. [PMID: 34155961 PMCID: PMC8220434 DOI: 10.1080/07391102.2021.1937319] [Citation(s) in RCA: 5] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/16/2022]
Abstract
Since the onset of global pandemic, the most focused research currently in progress is the development of potential drug candidates and clinical trials of existing FDA approved drugs for other relevant diseases, in order to repurpose them for the COVID-19. At the same time, several high throughput screenings of drugs have been reported to inhibit the viral components during the early course of infection but with little proven efficacies. Here, we investigate the drug repurposing strategies to counteract the coronavirus infection which involves several potential targetable host proteins involved in viral replication and disease progression. We report the high throughput analysis of literature-derived repurposing drug candidates that can be used to target the genetic regulators known to interact with viral proteins based on experimental and interactome studies. In this work we have performed integrated molecular docking followed by molecular dynamics (MD) simulations and free energy calculations through an expedite in silico process where the number of screened candidates reduces sequentially at every step based on physicochemical interactions. We elucidate that in addition to the pre-clinical and FDA approved drugs that targets specific regulatory proteins, a range of chemical compounds (Nafamostat, Chloramphenicol, Ponatinib) binds to the other gene transcription and translation regulatory proteins with higher affinity and may harbour potential for therapeutic uses. There is a rapid growing interest in the development of combination therapy for COVID-19 to target multiple enzymes/pathways. Our in silico approach would be useful in generating leads for experimental screening for rapid drug repurposing against SARS-CoV-2 interacting host proteins.Communicated by Ramaswamy H. Sarma.
Collapse
Affiliation(s)
- Amit Kumawat
- Prescience Insilico Private Limited, Bangalore, India,Department of Chemical Engineering, Indian Institute of Technology, Kanpur, India
| | | | - Debabrata Pramanik
- Department of Chemical Engineering, Indian Institute of Technology, Kanpur, India
| | - Sudip Roy
- Prescience Insilico Private Limited, Bangalore, India,CONTACT Sudip Roy ;
| | - Jayant K. Singh
- Prescience Insilico Private Limited, Bangalore, India,Department of Chemical Engineering, Indian Institute of Technology, Kanpur, India,Jayant K. Singh
| |
Collapse
|
22
|
Janoš P, Magistrato A. Role of Monovalent Ions in the NKCC1 Inhibition Mechanism Revealed through Molecular Simulations. Int J Mol Sci 2022; 23:ijms232315439. [PMID: 36499764 PMCID: PMC9741434 DOI: 10.3390/ijms232315439] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/28/2022] [Revised: 11/30/2022] [Accepted: 12/02/2022] [Indexed: 12/12/2022] Open
Abstract
The secondary active Na-K-Cl cotransporter 1 (NKCC1) promotes electroneutral uptake of two chloride ions, one sodium ion and one potassium ion. NKCC1 regulates Cl- homeostasis, thus being implicated in transepithelial water transport and in neuronal excitability. Aberrant NKCC1 transport is linked to a variety of human diseases. The loop diuretic drugs bumetanide, furosemide, azosemide and ethacrynic acid target NKCC1, but are characterized by poor selectivity leading to severe side effects. Despite its therapeutic importance, the molecular details of the NKCC1 inhibition mechanism remain unclear. Using all-atom simulations, we predict a putative binding mode of these drugs to the zebrafish (z) and human (h) NKCC1 orthologs. Although differing in their specific interactions with NKCC1 and/or monovalent ions, all drugs can fit within the same cavity and engage in hydrophobic interactions with M304/M382 in z/hNKCC1, a proposed ion gating residue demonstrated to be key for bumetanide binding. Consistent with experimental evidence, all drugs take advantage of the K+/Na+ ions, which plastically respond to their binding. This study not only provides atomic-level insights useful for drug discovery campaigns of more selective/potent NKCC1 inhibitors aimed to tackle diseases related to deregulated Cl- homeostasis, but it also supplies a paradigmatic example of the key importance of dynamical effects when drug binding is mediated by monovalent ions.
Collapse
|
23
|
Suh D, Feng S, Lee H, Zhang H, Park S, Kim S, Lee J, Choi S, Im W. CHARMM-GUI Enhanced Sampler for various collective variables and enhanced sampling methods. Protein Sci 2022; 31:e4446. [PMID: 36124940 PMCID: PMC9601830 DOI: 10.1002/pro.4446] [Citation(s) in RCA: 7] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/09/2022] [Revised: 09/06/2022] [Accepted: 09/07/2022] [Indexed: 11/08/2022]
Abstract
Enhanced sampling methodologies modifying underlying Hamiltonians can be used for the systems with a rugged potential energy surface that makes it hard to observe convergence using conventional unbiased molecular dynamics (MD) simulations. We present CHARMM-GUI Enhanced Sampler, a web-based tool to prepare various enhanced sampling simulations inputs with user-selected collective variables (CVs). Enhanced Sampler provides inputs for the following nine methods: accelerated MD, Gaussian accelerated MD, conformational flooding, metadynamics, adaptive biasing force, steered MD, temperature replica exchange MD, replica exchange solute tempering 2, and replica exchange umbrella sampling for the method-implemented MD packages including AMBER, CHARMM, GENESIS, GROMACS, NAMD, and OpenMM. Users only need to select a group of atoms via intuitive web-implementation in order to define commonly used nine CVs of interest: center of mass based distance, angle, dihedral, root-mean-square-distance, radius of gyration, distance projected on axis, two types of angles projected on axis, and coordination numbers. The enhanced sampling methods are tested with several biological systems to illustrate their efficiency over conventional MD. Enhanced Sampler with carefully optimized system-dependent parameters will help users to get meaningful results from their enhanced sampling simulations.
Collapse
Affiliation(s)
- Donghyuk Suh
- Department of Biological Sciences, Chemistry, Bioengineering, and Computer Science and EngineeringLehigh UniversityBethlehemPennsylvaniaUSA
- Research Institute for Pharmaceutical Sciences, College of Pharmacy and Graduate School of Pharmaceutical SciencesEwha Womans UniversitySeoulRepublic of Korea
| | - Shasha Feng
- Department of Biological Sciences, Chemistry, Bioengineering, and Computer Science and EngineeringLehigh UniversityBethlehemPennsylvaniaUSA
| | - Hwayoung Lee
- Department of Biological Sciences, Chemistry, Bioengineering, and Computer Science and EngineeringLehigh UniversityBethlehemPennsylvaniaUSA
| | - Han Zhang
- Department of Biological Sciences, Chemistry, Bioengineering, and Computer Science and EngineeringLehigh UniversityBethlehemPennsylvaniaUSA
| | - Sang‐Jun Park
- Department of Biological Sciences, Chemistry, Bioengineering, and Computer Science and EngineeringLehigh UniversityBethlehemPennsylvaniaUSA
| | - Seonghan Kim
- Department of Biological Sciences, Chemistry, Bioengineering, and Computer Science and EngineeringLehigh UniversityBethlehemPennsylvaniaUSA
| | - Jumin Lee
- Department of Biological Sciences, Chemistry, Bioengineering, and Computer Science and EngineeringLehigh UniversityBethlehemPennsylvaniaUSA
| | - Sun Choi
- Research Institute for Pharmaceutical Sciences, College of Pharmacy and Graduate School of Pharmaceutical SciencesEwha Womans UniversitySeoulRepublic of Korea
| | - Wonpil Im
- Department of Biological Sciences, Chemistry, Bioengineering, and Computer Science and EngineeringLehigh UniversityBethlehemPennsylvaniaUSA
| |
Collapse
|
24
|
Abstract
Gaussian accelerated molecular dynamics (GaMD) is a computational technique that provides both unconstrained enhanced sampling and free energy calculations of biomolecules. Here, we present the implementation of GaMD in the OpenMM simulation package and validate it on model systems of alanine dipeptide and RNA folding. For alanine dipeptide, 30 ns GaMD production simulations reproduced free energy profiles of 1000 ns conventional molecular dynamics (cMD) simulations. In addition, GaMD simulations captured the folding pathways of three hyperstable RNA tetraloops (UUCG, GCAA, and CUUG) and binding of the rbt203 ligand to the HIV-1 Tar RNA, both of which involved critical electrostatic interactions such as hydrogen bonding and base stacking. Together with previous implementations, GaMD in OpenMM will allow for wider applications in simulations of proteins, RNA, and other biomolecules.
Collapse
Affiliation(s)
- Matthew Copeland
- Center for Computational Biology and Department of Molecular Biosciences, University of Kansas, Lawrence, KS 66047
| | - Hung N. Do
- Center for Computational Biology and Department of Molecular Biosciences, University of Kansas, Lawrence, KS 66047
| | - Lane Votapka
- Department of Chemistry and Biochemistry, University of California at San Diego, La Jolla, CA 92093
| | - Keya Joshi
- Center for Computational Biology and Department of Molecular Biosciences, University of Kansas, Lawrence, KS 66047
| | - Jinan Wang
- Center for Computational Biology and Department of Molecular Biosciences, University of Kansas, Lawrence, KS 66047
| | - Rommie E. Amaro
- Department of Chemistry and Biochemistry, University of California at San Diego, La Jolla, CA 92093
| | - Yinglong Miao
- Center for Computational Biology and Department of Molecular Biosciences, University of Kansas, Lawrence, KS 66047,To whom correspondence should be addressed:
| |
Collapse
|
25
|
Wang Y, Herron L, Tiwary P. From data to noise to data for mixing physics across temperatures with generative artificial intelligence. Proc Natl Acad Sci U S A 2022; 119:e2203656119. [PMID: 35925885 DOI: 10.1073/pnas.2203656119] [Citation(s) in RCA: 8] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [What about the content of this article? (0)] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/06/2023] Open
Abstract
Using simulations or experiments performed at some set of temperatures to learn about the physics or chemistry at some other arbitrary temperature is a problem of immense practical and theoretical relevance. Here we develop a framework based on statistical mechanics and generative artificial intelligence that allows solving this problem. Specifically, we work with denoising diffusion probabilistic models and show how these models in combination with replica exchange molecular dynamics achieve superior sampling of the biomolecular energy landscape at temperatures that were never simulated without assuming any particular slow degrees of freedom. The key idea is to treat the temperature as a fluctuating random variable and not a control parameter as is usually done. This allows us to directly sample from the joint probability distribution in configuration and temperature space. The results here are demonstrated for a chirally symmetric peptide and single-strand RNA undergoing conformational transitions in all-atom water. We demonstrate how we can discover transition states and metastable states that were previously unseen at the temperature of interest and even bypass the need to perform further simulations for a wide range of temperatures. At the same time, any unphysical states are easily identifiable through very low Boltzmann weights. The procedure while shown here for a class of molecular simulations should be more generally applicable to mixing information across simulations and experiments with varying control parameters.
Collapse
|
26
|
Ahmad K, Rizzi A, Capelli R, Mandelli D, Lyu W, Carloni P. Enhanced-Sampling Simulations for the Estimation of Ligand Binding Kinetics: Current Status and Perspective. Front Mol Biosci 2022; 9:899805. [PMID: 35755817 PMCID: PMC9216551 DOI: 10.3389/fmolb.2022.899805] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/19/2022] [Accepted: 05/09/2022] [Indexed: 12/12/2022] Open
Abstract
The dissociation rate (k off) associated with ligand unbinding events from proteins is a parameter of fundamental importance in drug design. Here we review recent major advancements in molecular simulation methodologies for the prediction of k off. Next, we discuss the impact of the potential energy function models on the accuracy of calculated k off values. Finally, we provide a perspective from high-performance computing and machine learning which might help improve such predictions.
Collapse
Affiliation(s)
- Katya Ahmad
- Computational Biomedicine (IAS-5/INM-9), Forschungszentrum Jülich, Jülich, Germany
| | - Andrea Rizzi
- Computational Biomedicine (IAS-5/INM-9), Forschungszentrum Jülich, Jülich, Germany
- Atomistic Simulations, Istituto Italiano di Tecnologia, Genova, Italy
| | - Riccardo Capelli
- Department of Applied Science and Technology (DISAT), Politecnico di Torino, Torino, Italy
| | - Davide Mandelli
- Computational Biomedicine (IAS-5/INM-9), Forschungszentrum Jülich, Jülich, Germany
| | - Wenping Lyu
- Warshel Institute for Computational Biology, School of Life and Health Sciences, The Chinese University of Hong Kong (Shenzhen), Shenzhen, China
- School of Chemistry and Materials Science, University of Science and Technology of China, Hefei, China
| | - Paolo Carloni
- Computational Biomedicine (IAS-5/INM-9), Forschungszentrum Jülich, Jülich, Germany
- Molecular Neuroscience and Neuroimaging (INM-11), Forschungszentrum Jülich, Jülich, Germany
| |
Collapse
|
27
|
Wu X, Xu LY, Li EM, Dong G. Application of molecular dynamics simulation in biomedicine. Chem Biol Drug Des 2022; 99:789-800. [PMID: 35293126 DOI: 10.1111/cbdd.14038] [Citation(s) in RCA: 15] [Impact Index Per Article: 7.5] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/29/2021] [Revised: 02/25/2022] [Accepted: 03/05/2022] [Indexed: 02/05/2023]
Abstract
Molecular dynamics (MD) simulation has been widely used in the field of biomedicine to study the conformational transition of proteins caused by mutation or ligand binding/unbinding. It provides some perspectives those are difficult to find in traditional biochemical or pathological experiments, for example, detailed effects of mutations on protein structure and protein-protein/ligand interaction at the atomic level. In this review, a broad overview on conformation changes and drug discovery by MD simulation is given. We first discuss the preparation of protein structure for MD simulation, which is a key step that determines the accuracy of the simulation. Then, we summarize the applications of commonly used force fields and MD simulations in scientific research. Finally, enhanced sampling methods and common applications of these methods are introduced. In brief, MD simulation is a powerful tool and it can be used to guide experimental study. The combination of MD simulation and experimental techniques is an a priori means to solve the biomedical problems and give a deep understanding on the relationship between protein structure and function.
Collapse
Affiliation(s)
- Xiaodong Wu
- Department of Biochemistry and Molecular Biology, Shantou University Medical College, Shantou, China
| | - Li-Yan Xu
- Key Laboratory of Molecular Biology in High Cancer Incidence Coastal Area of Guangdong Higher Education Institutes, Shantou University Medical College, Shantou, China
- Cancer Research Center, Shantou University Medical College, Shantou, China
| | - En-Min Li
- Department of Biochemistry and Molecular Biology, Shantou University Medical College, Shantou, China
- Key Laboratory of Molecular Biology in High Cancer Incidence Coastal Area of Guangdong Higher Education Institutes, Shantou University Medical College, Shantou, China
| | - Geng Dong
- Department of Biochemistry and Molecular Biology, Shantou University Medical College, Shantou, China
- Medical Informatics Research Center, Shantou University Medical College, Shantou, China
| |
Collapse
|
28
|
Romero-Rivera A, Corbella M, Parracino A, Patrick WM, Kamerlin SCL. Complex Loop Dynamics Underpin Activity, Specificity, and Evolvability in the (βα) 8 Barrel Enzymes of Histidine and Tryptophan Biosynthesis. JACS Au 2022; 2:943-960. [PMID: 35557756 PMCID: PMC9088769 DOI: 10.1021/jacsau.2c00063] [Citation(s) in RCA: 5] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 01/29/2022] [Revised: 03/16/2022] [Accepted: 03/18/2022] [Indexed: 05/16/2023]
Abstract
Enzymes are conformationally dynamic, and their dynamical properties play an important role in regulating their specificity and evolvability. In this context, substantial attention has been paid to the role of ligand-gated conformational changes in enzyme catalysis; however, such studies have focused on tremendously proficient enzymes such as triosephosphate isomerase and orotidine 5'-monophosphate decarboxylase, where the rapid (μs timescale) motion of a single loop dominates the transition between catalytically inactive and active conformations. In contrast, the (βα)8-barrels of tryptophan and histidine biosynthesis, such as the specialist isomerase enzymes HisA and TrpF, and the bifunctional isomerase PriA, are decorated by multiple long loops that undergo conformational transitions on the ms (or slower) timescale. Studying the interdependent motions of multiple slow loops, and their role in catalysis, poses a significant computational challenge. This work combines conventional and enhanced molecular dynamics simulations with empirical valence bond simulations to provide rich details of the conformational behavior of the catalytic loops in HisA, PriA, and TrpF, and the role of their plasticity in facilitating bifunctionality in PriA and evolved HisA variants. In addition, we demonstrate that, similar to other enzymes activated by ligand-gated conformational changes, loops 3 and 4 of HisA and PriA act as gripper loops, facilitating the isomerization of the large bulky substrate ProFAR, albeit now on much slower timescales. This hints at convergent evolution on these different (βα)8-barrel scaffolds. Finally, our work reemphasizes the potential of engineering loop dynamics as a tool to artificially manipulate the catalytic repertoire of TIM-barrel proteins.
Collapse
Affiliation(s)
- Adrian Romero-Rivera
- Department
of Chemistry—BMC, Uppsala University, BMC Box 576, S-751 23 Uppsala, Sweden
| | - Marina Corbella
- Department
of Chemistry—BMC, Uppsala University, BMC Box 576, S-751 23 Uppsala, Sweden
| | - Antonietta Parracino
- Department
of Chemistry—BMC, Uppsala University, BMC Box 576, S-751 23 Uppsala, Sweden
| | - Wayne M. Patrick
- Centre
for Biodiscovery, School of Biological Sciences, Victoria University of Wellington, 6012 Wellington, New Zealand
| | | |
Collapse
|
29
|
Dokainish HM, Re S, Mori T, Kobayashi C, Jung J, Sugita Y. The inherent flexibility of receptor binding domains in SARS-CoV-2 spike protein. eLife 2022; 11:e75720. [PMID: 35323112 PMCID: PMC8963885 DOI: 10.7554/elife.75720] [Citation(s) in RCA: 31] [Impact Index Per Article: 15.5] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/20/2021] [Accepted: 03/15/2022] [Indexed: 12/17/2022] Open
Abstract
Spike (S) protein is the primary antigenic target for neutralization and vaccine development for the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). It decorates the virus surface and undergoes large motions of its receptor binding domains (RBDs) to enter the host cell. Here, we observe Down, one-Up, one-Open, and two-Up-like structures in enhanced molecular dynamics simulations, and characterize the transition pathways via inter-domain interactions. Transient salt-bridges between RBDA and RBDC and the interaction with glycan at N343B support RBDA motions from Down to one-Up. Reduced interactions between RBDA and RBDB in one-Up induce RBDB motions toward two-Up. The simulations overall agree with cryo-electron microscopy structure distributions and FRET experiments and provide hidden functional structures, namely, intermediates along Down-to-one-Up transition with druggable cryptic pockets as well as one-Open with a maximum exposed RBD. The inherent flexibility of S-protein thus provides essential information for antiviral drug rational design or vaccine development.
Collapse
Grants
- FLAGSHIP 2020 project Ministry of Education, Culture, Sports, Science and Technology
- 19K06532 Ministry of Education, Culture, Sports, Science and Technology
- Dynamic Structural Biology/Glycolipidologue Initiative/Biology of Intracellular Environments RIKEN
- Priority Issue on Post-K computer Ministry of Education, Culture, Sports, Science and Technology
- Program for Promoting Researches on the Supercomputer Fugaku Ministry of Education, Culture, Sports, Science and Technology
- JPMXP1020200101 Ministry of Education, Culture, Sports, Science and Technology
- JPMXP1020200201 Ministry of Education, Culture, Sports, Science and Technology
- 19H05645 Ministry of Education, Culture, Sports, Science and Technology
- 21H05249 Ministry of Education, Culture, Sports, Science and Technology
- 20K15737 Ministry of Education, Culture, Sports, Science and Technology
- 19K12229 Ministry of Education, Culture, Sports, Science and Technology
- 21H05157 Ministry of Education, Culture, Sports, Science and Technology
- hp200135 HPCI System Research project
- hp200153 HPCI System Research project
- hp200028 HPCI System Research project
- hp210107 HPCI System Research project
- hp210177 HPCI System Research project
Collapse
Affiliation(s)
- Hisham M Dokainish
- Theoretical Molecular Science Laboratory, RIKEN Cluster for Pioneering ResearchWakoJapan
| | - Suyong Re
- Artificial Intelligence Center for Health and Biomedical Research, National Institutes of Biomedical Innovation, Health and NutritionOsakaJapan
- Laboratory for Biomolecular Function Simulation, RIKEN Center for Biosystems Dynamics ResearchKobeJapan
| | - Takaharu Mori
- Theoretical Molecular Science Laboratory, RIKEN Cluster for Pioneering ResearchWakoJapan
| | - Chigusa Kobayashi
- Computational Biophysics Research Team, RIKEN Center for Computational ScienceKobeJapan
| | - Jaewoon Jung
- Theoretical Molecular Science Laboratory, RIKEN Cluster for Pioneering ResearchWakoJapan
- Computational Biophysics Research Team, RIKEN Center for Computational ScienceKobeJapan
| | - Yuji Sugita
- Theoretical Molecular Science Laboratory, RIKEN Cluster for Pioneering ResearchWakoJapan
- Laboratory for Biomolecular Function Simulation, RIKEN Center for Biosystems Dynamics ResearchKobeJapan
- Computational Biophysics Research Team, RIKEN Center for Computational ScienceKobeJapan
| |
Collapse
|
30
|
Abstract
Chemokine receptors are key G-protein-coupled receptors (GPCRs) that control cell migration in immune system responses, development of cardiovascular and central nervous systems, and numerous diseases. In particular, the CXCR4 chemokine receptor promotes metastasis, tumor growth and angiogenesis in cancers. CXCR4 is also used as one of the two co-receptors for T-tropic HIV-1 entry into host cells. Therefore, CXCR4 serves as an important therapeutic target for treating cancers and HIV infection. Apart from the CXCL12 endogenous peptide agonist, previous studies suggested that the first 17 amino acids of CXCL12 are sufficient to activate CXCR4. Two 17-residue peptides with positions 1-4 mutated to RSVM and ASLW functioned as super and partial agonists of CXCR4, respectively. However, the mechanism of peptide agonist binding in CXCR4 remains unclear. Here, we have investigated this mechanism through all-atom simulations using a novel Peptide Gaussian accelerated molecular dynamics (Pep-GaMD) method. The Pep-GaMD simulations have allowed us to explore representative binding conformations of each peptide and identify critical low-energy states of CXCR4 activated by the super versus partial peptide agonists. Our simulations have provided important mechanistic insights into peptide agonist binding in CXCR4, which are expected to facilitate rational design of new peptide modulators of CXCR4 and other chemokine receptors.
Collapse
|
31
|
Bocus M, Vanduyfhuys L, De Proft F, Weckhuysen BM, Van Speybroeck V. Mechanistic Characterization of Zeolite-Catalyzed Aromatic Electrophilic Substitution at Realistic Operating Conditions. JACS Au 2022; 2:502-514. [PMID: 35252999 PMCID: PMC8889610 DOI: 10.1021/jacsau.1c00544] [Citation(s) in RCA: 8] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 12/03/2021] [Indexed: 05/11/2023]
Abstract
Zeolite-catalyzed benzene ethylation is an important industrial reaction, as it is the first step in the production of styrene for polymer manufacturing. Furthermore, it is a prototypical example of aromatic electrophilic substitution, a key reaction in the synthesis of many bulk and fine chemicals. Despite extensive research, the reaction mechanism and the nature of elusive intermediates at realistic operating conditions is not properly understood. More in detail, the existence of the elusive arenium ion (better known as Wheland complex) formed upon electrophilic attack on the aromatic ring is still a matter of debate. Temperature effects and the presence of protic guest molecules such as water are expected to impact the reaction mechanism and lifetime of the reaction intermediates. Herein, we used enhanced sampling ab initio molecular dynamics simulations to investigate the complete mechanism of benzene ethylation with ethene and ethanol in the H-ZSM-5 zeolite. We show that both the stepwise and concerted mechanisms are active at reaction conditions and that the Wheland intermediate spontaneously appears as a shallow minimum in the free energy surface after the electrophilic attack on the benzene ring. Addition of water enhances the protonation kinetics by about 1 order of magnitude at coverages of one water molecule per Brønsted acidic site. In the fully solvated regime, an overstabilization of the BAS as hydronium ion occurs and the rate enhancement disappears. The obtained results give critical atomistic insights in the role of water to selectively tune the kinetics of protonation reactions in zeolites.
Collapse
Affiliation(s)
- Massimo Bocus
- Center
for Molecular Modeling, Ghent University, Technologiepark 46, 9052 Zwijnaarde, Belgium
| | - Louis Vanduyfhuys
- Center
for Molecular Modeling, Ghent University, Technologiepark 46, 9052 Zwijnaarde, Belgium
| | - Frank De Proft
- Eenheid
Algemene Chemie (ALGC), Vrije Universiteit
Brussel, Pleinlaan 2, 1050 Brussels, Belgium
| | - Bert M. Weckhuysen
- Inorganic
Chemistry and Catalysis Group, Debye Institute for Nanomaterials Science, Utrecht University, Universiteitsweg 99, 3584 CG Utrecht, The Netherlands
| | | |
Collapse
|
32
|
Sun J, Raymundo MAV, Chang CA. Ritonavir and xk263 Binding-Unbinding with HIV-1 Protease: Pathways, Energy and Comparison. Life (Basel) 2022; 12:116. [PMID: 35054509 DOI: 10.3390/life12010116] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [What about the content of this article? (0)] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 12/10/2021] [Revised: 01/04/2022] [Accepted: 01/10/2022] [Indexed: 01/22/2023]
Abstract
Understanding non-covalent biomolecular recognition, which includes drug-protein bound states and their binding/unbinding processes, is of fundamental importance in chemistry, biology, and medicine. Fully revealing the factors that govern the binding/unbinding processes can further assist in designing drugs with desired binding kinetics. HIV protease (HIVp) plays an integral role in the HIV life cycle, so it is a prime target for drug therapy. HIVp has flexible flaps, and the binding pocket can be accessible by a ligand via various pathways. Comparing ligand association and dissociation pathways can help elucidate the ligand-protein interactions such as key residues directly involved in the interaction or specific protein conformations that determine the binding of a ligand under certain pathway(s). Here, we investigated the ligand unbinding process for a slow binder, ritonavir, and a fast binder, xk263, by using unbiased all-atom accelerated molecular dynamics (aMD) simulation with a re-seeding approach and an explicit solvent model. Using ritonavir-HIVp and xk263-HIVp ligand-protein systems as cases, we sampled multiple unbinding pathways for each ligand and observed that the two ligands preferred the same unbinding route. However, ritonavir required a greater HIVp motion to dissociate as compared with xk263, which can leave the binding pocket with little conformational change of HIVp. We also observed that ritonavir unbinding pathways involved residues which are associated with drug resistance and are distal from catalytic site. Analyzing HIVp conformations sampled during both ligand-protein binding and unbinding processes revealed significantly more overlapping HIVp conformations for ritonavir-HIVp rather than xk263-HIVp. However, many HIVp conformations are unique in xk263-HIVp unbinding processes. The findings are consistent with previous findings that xk263 prefers an induced-fit model for binding and unbinding, whereas ritonavir favors a conformation selection model. This study deepens our understanding of the dynamic process of ligand unbinding and provides insights into ligand-protein recognition mechanisms and drug discovery.
Collapse
|
33
|
Salas-Estrada L, Fiorillo B, Filizola M. Metadynamics simulations leveraged by statistical analyses and artificial intelligence-based tools to inform the discovery of G protein-coupled receptor ligands. Front Endocrinol (Lausanne) 2022; 13:1099715. [PMID: 36619585 PMCID: PMC9816996 DOI: 10.3389/fendo.2022.1099715] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Figures] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 11/16/2022] [Accepted: 12/12/2022] [Indexed: 12/25/2022] Open
Abstract
G Protein-Coupled Receptors (GPCRs) are a large family of membrane proteins with pluridimensional signaling profiles. They undergo ligand-specific conformational changes, which in turn lead to the differential activation of intracellular signaling proteins and the consequent triggering of a variety of biological responses. This conformational plasticity directly impacts our understanding of GPCR signaling and therapeutic implications, as do ligand-specific kinetic differences in GPCR-induced transducer activation/coupling or GPCR-transducer complex stability. High-resolution experimental structures of ligand-bound GPCRs in the presence or absence of interacting transducers provide important, yet limited, insights into the highly dynamic process of ligand-induced activation or inhibition of these receptors. We and others have complemented these studies with computational strategies aimed at characterizing increasingly accurate metastable conformations of GPCRs using a combination of metadynamics simulations, state-of-the-art algorithms for statistical analyses of simulation data, and artificial intelligence-based tools. This minireview provides an overview of these approaches as well as lessons learned from them towards the identification of conformational states that may be difficult or even impossible to characterize experimentally and yet important to discover new GPCR ligands.
Collapse
|
34
|
Abstract
INTRODUCTION Hidden allosteric sites are not visible in apo-crystal structures, but they may be visible in holo-structures when a certain ligand binds and maintains the ligand intended conformation. Several computational and experimental techniques have been used to investigate these hidden sites but identifying them remains a challenge. AREAS COVERED This review provides a summary of the many theoretical approaches for predicting hidden allosteric sites in disease-related proteins. Furthermore, promising cases have been thoroughly examined to reveal the hidden allosteric site and its modulator. EXPERT OPINION In the recent past, with the development in scientific techniques and bioinformatics tools, the number of drug targets for complex human diseases has significantly increased but unfortunately most of these targets are undruggable due to several reasons. Alternative strategies such as finding cryptic (hidden) allosteric sites are an attractive approach for exploitation of the discovery of new targets. These hidden sites are difficult to recognize compared to allosteric sites, mainly due to a lack of visibility in the crystal structure. In our opinion, after many years of development, MD simulations are finally becoming successful for obtaining a detailed molecular description of drug-target interaction.
Collapse
Affiliation(s)
- Ashfaq Ur Rehman
- Key Laboratory of Cell Differentiation and Apoptosis of Chinese Ministry of Education, Clinical and Fundamental Research Center, Renji Hospital, Shanghai Jiao-Tong University School of Medicine, Shanghai, China
| | - Shaoyong Lu
- Key Laboratory of Cell Differentiation and Apoptosis of Chinese Ministry of Education, Clinical and Fundamental Research Center, Renji Hospital, Shanghai Jiao-Tong University School of Medicine, Shanghai, China
| | - Abdul Aziz Khan
- Bio-X Institutes, Key Laboratory for the Genetics of Development and Neuropsychiatric Disorders (Ministry of Education), Shanghai Key Laboratory of Psychotic Disorders, and Brain Science and Technology Research Center, Institute of Psychology and Behavioral Sciences, Shanghai Jiao Tong University, Shanghai, China
| | - Beenish Khurshid
- Department of Biochemistry, Abdul Wali Khan University Mardan, Mardan, Pakistan
| | - Salman Rasheed
- National Center for Bioinformatics, Quaid-e-Azam University, Islamabad, Pakistan
| | - Abdul Wadood
- Department of Biochemistry, Abdul Wali Khan University Mardan, Mardan, Pakistan
| | - Jian Zhang
- Key Laboratory of Cell Differentiation and Apoptosis of Chinese Ministry of Education, Clinical and Fundamental Research Center, Renji Hospital, Shanghai Jiao-Tong University School of Medicine, Shanghai, China.,School of Pharmaceutical Sciences, Zhengzhou University, Zhengzhou, China
| |
Collapse
|
35
|
Chu X, Wang Y, Tian P, Li W, Mercadante D. Editorial: Advanced Sampling and Modeling in Molecular Simulations for Slow and Large-Scale Biomolecular Dynamics. Front Mol Biosci 2021; 8:795991. [PMID: 34869608 PMCID: PMC8633950 DOI: 10.3389/fmolb.2021.795991] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Key Words] [Track Full Text] [Download PDF] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/15/2021] [Accepted: 10/23/2021] [Indexed: 11/13/2022] Open
Affiliation(s)
- Xiakun Chu
- Department of Chemistry, State University of New York, Stony Brook, NY, United States
| | - Yong Wang
- College of Life Sciences, Shanghai Institute for Advanced Study, Institute of Quantitative Biology, Zhejiang University, Hangzhou, China
| | | | - Wenfei Li
- National Laboratory of Solid State Microstructure, Department of Physics, Nanjing University, Nanjing, China
- Wenzhou Institute, University of Chinese Academy of Sciences, Wenzhou, China
| | - Davide Mercadante
- School of Chemical Sciences, The University of Auckland, Auckland, New Zealand
| |
Collapse
|
36
|
Shu Z, Wu M, Liao J, Chen C. FSATOOL 2.0: An integrated molecular dynamics simulation and trajectory data analysis program. J Comput Chem 2021; 43:215-224. [PMID: 34751974 DOI: 10.1002/jcc.26772] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/09/2021] [Revised: 09/30/2021] [Accepted: 10/04/2021] [Indexed: 11/08/2022]
Abstract
Molecular dynamics simulation is important in the computational study of the biomolecules. In this paper, we upgrade our previous FSATOOL to version 2.0. It is no longer a plugin as before. Besides the existed enhanced sampling and Markov state model analysis module, FSATOOL 2.0 has three new features now. First, it contains a molecular dynamics simulation engine on both CPU and GPU device. The engine works with an embedded enhanced sampling module. Second, it can do the free energy calculation by various practical methods, including the weighted histogram analysis method and Gaussian mixture model. Third, it has many subroutines to process the trajectory data, such as principal component analysis, time-structure based independent component analysis, contact analysis, and Φ-value analysis. Most importantly, all these calculations are integrated into one package. The trajectory data format is compatible with all the modules. With a proper input parameter file, users can do the molecular dynamics simulation and data analysis work by only a few simplified commands. The capabilities and theoretical backgrounds of FSATOOL 2.0 are introduced in the paper.
Collapse
Affiliation(s)
- Zirui Shu
- Biomolecular Physics and Modeling Group, School of Physics, Huazhong University of Science and Technology, Wuhan, Hubei, China
| | - Mincong Wu
- Biomolecular Physics and Modeling Group, School of Physics, Huazhong University of Science and Technology, Wuhan, Hubei, China
| | - Jun Liao
- Biomolecular Physics and Modeling Group, School of Physics, Huazhong University of Science and Technology, Wuhan, Hubei, China
| | - Changjun Chen
- Biomolecular Physics and Modeling Group, School of Physics, Huazhong University of Science and Technology, Wuhan, Hubei, China
| |
Collapse
|
37
|
Nin-Hill A, Mueller NPF, Molteni C, Rovira C, Alfonso-Prieto M. Photopharmacology of Ion Channels through the Light of the Computational Microscope. Int J Mol Sci 2021; 22:12072. [PMID: 34769504 PMCID: PMC8584574 DOI: 10.3390/ijms222112072] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/10/2021] [Revised: 10/31/2021] [Accepted: 11/02/2021] [Indexed: 12/13/2022] Open
Abstract
The optical control and investigation of neuronal activity can be achieved and carried out with photoswitchable ligands. Such compounds are designed in a modular fashion, combining a known ligand of the target protein and a photochromic group, as well as an additional electrophilic group for tethered ligands. Such a design strategy can be optimized by including structural data. In addition to experimental structures, computational methods (such as homology modeling, molecular docking, molecular dynamics and enhanced sampling techniques) can provide structural insights to guide photoswitch design and to understand the observed light-regulated effects. This review discusses the application of such structure-based computational methods to photoswitchable ligands targeting voltage- and ligand-gated ion channels. Structural mapping may help identify residues near the ligand binding pocket amenable for mutagenesis and covalent attachment. Modeling of the target protein in a complex with the photoswitchable ligand can shed light on the different activities of the two photoswitch isomers and the effect of site-directed mutations on photoswitch binding, as well as ion channel subtype selectivity. The examples presented here show how the integration of computational modeling with experimental data can greatly facilitate photoswitchable ligand design and optimization. Recent advances in structural biology, both experimental and computational, are expected to further strengthen this rational photopharmacology approach.
Collapse
Affiliation(s)
- Alba Nin-Hill
- Departament de Química Inorgànica i Orgànica (Secció de Química Orgànica) and Institut de Química Teòrica i Computacional (IQTCUB), Universitat de Barcelona, 08028 Barcelona, Spain; (A.N.-H.); (C.R.)
| | - Nicolas Pierre Friedrich Mueller
- Institute for Advanced Simulations IAS-5 and Institute of Neuroscience and Medicine INM-9, Computational Biomedicine, Forschungszentrum Jülich, 52425 Jülich, Germany;
- Faculty of Mathematics and Natural Sciences, Heinrich-Heine-University Düsseldorf, Universitätsstr. 1, 40225 Düsseldorf, Germany
| | - Carla Molteni
- Physics Department, King’s College London, London WC2R 2LS, UK;
| | - Carme Rovira
- Departament de Química Inorgànica i Orgànica (Secció de Química Orgànica) and Institut de Química Teòrica i Computacional (IQTCUB), Universitat de Barcelona, 08028 Barcelona, Spain; (A.N.-H.); (C.R.)
- Institució Catalana de Recerca i Estudis Avançats (ICREA), 08020 Barcelona, Spain
| | - Mercedes Alfonso-Prieto
- Institute for Advanced Simulations IAS-5 and Institute of Neuroscience and Medicine INM-9, Computational Biomedicine, Forschungszentrum Jülich, 52425 Jülich, Germany;
- Cécile and Oskar Vogt Institute for Brain Research, University Hospital Düsseldorf, Medical Faculty, Heinrich Heine University Düsseldorf, 40225 Düsseldorf, Germany
| |
Collapse
|
38
|
Abstract
The development of enhanced sampling methods has greatly extended the scope of atomistic simulations, allowing long-time phenomena to be studied with accessible computational resources. Many such methods rely on the identification of an appropriate set of collective variables. These are meant to describe the system's modes that most slowly approach equilibrium under the action of the sampling algorithm. Once identified, the equilibration of these modes is accelerated by the enhanced sampling method of choice. An attractive way of determining the collective variables is to relate them to the eigenfunctions and eigenvalues of the transfer operator. Unfortunately, this requires knowing the long-term dynamics of the system beforehand, which is generally not available. However, we have recently shown that it is indeed possible to determine efficient collective variables starting from biased simulations. In this paper, we bring the power of machine learning and the efficiency of the recently developed on the fly probability-enhanced sampling method to bear on this approach. The result is a powerful and robust algorithm that, given an initial enhanced sampling simulation performed with trial collective variables or generalized ensembles, extracts transfer operator eigenfunctions using a neural network ansatz and then accelerates them to promote sampling of rare events. To illustrate the generality of this approach, we apply it to several systems, ranging from the conformational transition of a small molecule to the folding of a miniprotein and the study of materials crystallization.
Collapse
Affiliation(s)
- Luigi Bonati
- Department of Physics, Eidgenössische Technische Hochschule (ETH) Zürich, 8092 Zürich, Switzerland;
- Atomistic Simulations, Italian Institute of Technology, 16163 Genova, Italy
| | | | - Michele Parrinello
- Atomistic Simulations, Italian Institute of Technology, 16163 Genova, Italy;
| |
Collapse
|
39
|
González-Fernández C, Basauri A, Fallanza M, Bringas E, Oostenbrink C, Ortiz I. Fighting Against Bacterial Lipopolysaccharide-Caused Infections through Molecular Dynamics Simulations: A Review. J Chem Inf Model 2021; 61:4839-4851. [PMID: 34559524 PMCID: PMC8549069 DOI: 10.1021/acs.jcim.1c00613] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/07/2023]
Abstract
![]()
Lipopolysaccharide
(LPS) is the primary component of the outer
leaflet of Gram-negative bacterial outer membranes. LPS elicits an
overwhelming immune response during infection, which can lead to life-threatening
sepsis or septic shock for which no suitable treatment is available
so far. As a result of the worldwide expanding multidrug-resistant
bacteria, the occurrence and frequency of sepsis are expected to increase;
thus, there is an urge to develop novel strategies for treating bacterial
infections. In this regard, gaining an in-depth understanding about
the ability of LPS to both stimulate the host immune system and interact
with several molecules is crucial for fighting against LPS-caused
infections and allowing for the rational design of novel antisepsis
drugs, vaccines and LPS sequestration and detection methods. Molecular
dynamics (MD) simulations, which are understood as being a computational
microscope, have proven to be of significant value to understand LPS-related
phenomena, driving and optimizing experimental research studies. In
this work, a comprehensive review on the methods that can be combined
with MD simulations, recently applied in LPS research, is provided.
We focus especially on both enhanced sampling methods, which enable
the exploration of more complex systems and access to larger time
scales, and free energy calculation approaches. Thereby, apart from
outlining several strategies for surmounting LPS-caused infections,
this work reports the current state-of-the-art of the methods applied
with MD simulations for moving a step forward in the development of
such strategies.
Collapse
Affiliation(s)
- Cristina González-Fernández
- Department of Chemical and Biomolecular Engineering, ETSIIT, University of Cantabria, Avda. Los Castros s/n, 39005 Santander, Spain
| | - Arantza Basauri
- Department of Chemical and Biomolecular Engineering, ETSIIT, University of Cantabria, Avda. Los Castros s/n, 39005 Santander, Spain
| | - Marcos Fallanza
- Department of Chemical and Biomolecular Engineering, ETSIIT, University of Cantabria, Avda. Los Castros s/n, 39005 Santander, Spain
| | - Eugenio Bringas
- Department of Chemical and Biomolecular Engineering, ETSIIT, University of Cantabria, Avda. Los Castros s/n, 39005 Santander, Spain
| | - Chris Oostenbrink
- Institute for Molecular Modeling and Simulation, BOKU - University of Natural Resources and Life Sciences, Muthgasse 18, 1190 Vienna, Austria
| | - Inmaculada Ortiz
- Department of Chemical and Biomolecular Engineering, ETSIIT, University of Cantabria, Avda. Los Castros s/n, 39005 Santander, Spain
| |
Collapse
|
40
|
Abstract
Multiscale enhanced sampling (MSES) allows for an enhanced sampling of all-atom protein structures by coupling with the accelerated dynamics of the associated coarse-grained (CG) model. In this paper, we propose an MSES extension to replace the CG model with the dynamics on the reduced subspace generated by a machine learning approach, the variational autoencoder (VAE). The molecular dynamic (MD) trajectories of the ribose-binding protein (RBP) in both the closed and open forms were used as the input by extracting the inter-residue distances as the structural features in order to train the VAE model, allowing the encoded latent layer to characterize the difference in the structural dynamics of the closed and open forms. The interpolated data characterizing the RBP structural change in between the closed and open forms were thus efficiently generated in the low-dimensional latent space of the VAE, which was then decoded into the time-series data of the inter-residue distances and was useful for driving the structural sampling at an atomistic resolution via the MSES scheme. The free energy surfaces on the latent space demonstrated the refinement of the generated data that had a single basin into the simulated data containing two closed and open basins, thus illustrating the usefulness of the MD simulation together with the molecular mechanics force field in recovering the correct structural ensemble.
Collapse
Affiliation(s)
- Kei Moritsugu
- Graduate School of Medical Life Science, Yokohama City University, Yokohama 230-0045, Japan
| |
Collapse
|
41
|
Gong X, Zhang Y, Chen J. Advanced Sampling Methods for Multiscale Simulation of Disordered Proteins and Dynamic Interactions. Biomolecules 2021; 11:1416. [PMID: 34680048 PMCID: PMC8533332 DOI: 10.3390/biom11101416] [Citation(s) in RCA: 12] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/31/2021] [Revised: 09/22/2021] [Accepted: 09/24/2021] [Indexed: 11/16/2022] Open
Abstract
Intrinsically disordered proteins (IDPs) are highly prevalent and play important roles in biology and human diseases. It is now also recognized that many IDPs remain dynamic even in specific complexes and functional assemblies. Computer simulations are essential for deriving a molecular description of the disordered protein ensembles and dynamic interactions for a mechanistic understanding of IDPs in biology, diseases, and therapeutics. Here, we provide an in-depth review of recent advances in the multi-scale simulation of disordered protein states, with a particular emphasis on the development and application of advanced sampling techniques for studying IDPs. These techniques are critical for adequate sampling of the manifold functionally relevant conformational spaces of IDPs. Together with dramatically improved protein force fields, these advanced simulation approaches have achieved substantial success and demonstrated significant promise towards the quantitative and predictive modeling of IDPs and their dynamic interactions. We will also discuss important challenges remaining in the atomistic simulation of larger systems and how various coarse-grained approaches may help to bridge the remaining gaps in the accessible time- and length-scales of IDP simulations.
Collapse
Affiliation(s)
- Xiping Gong
- Department of Chemistry, University of Massachusetts Amherst, Amherst, MA 01003, USA; (X.G.); (Y.Z.)
| | - Yumeng Zhang
- Department of Chemistry, University of Massachusetts Amherst, Amherst, MA 01003, USA; (X.G.); (Y.Z.)
| | - Jianhan Chen
- Department of Chemistry, University of Massachusetts Amherst, Amherst, MA 01003, USA; (X.G.); (Y.Z.)
- Department of Biochemistry and Molecular Biology, University of Massachusetts Amherst, Amherst, MA 01003, USA
| |
Collapse
|
42
|
Hayes RL, Brooks CL. A strategy for proline and glycine mutations to proteins with alchemical free energy calculations. J Comput Chem 2021; 42:1088-1094. [PMID: 33844328 DOI: 10.1002/jcc.26525] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/27/2020] [Revised: 03/03/2021] [Accepted: 03/05/2021] [Indexed: 11/07/2022]
Abstract
Computation of the thermodynamic consequences of protein mutations holds great promise in protein biophysics and design. Alchemical free energy methods can give improved estimates of mutational free energies, and are already widely used in calculations of relative and absolute binding free energies in small molecule design problems. In principle, alchemical methods can address any amino acid mutation with an appropriate alchemical pathway, but identifying a strategy that produces such a path for proline and glycine mutations is an ongoing challenge. Most current strategies perturb only side chain atoms, while proline and glycine mutations also alter the backbone parameters and backbone ring topology. Some strategies also perturb backbone parameters and enable glycine mutations. This work presents a strategy that enables both proline and glycine mutations and comprises two key elements: a dual backbone with restraints and scaling of bonded terms, facilitating backbone parameter changes, and a soft bond in the proline ring, enabling ring topology changes in proline mutations. These elements also have utility for core hopping and macrocycle studies in computer-aided drug design. This new strategy shows slight improvements over an alternative side chain perturbation strategy for a set T4 lysozyme mutations lacking proline and glycine, and yields good agreement with experiment for a set of T4 lysozyme proline and glycine mutations not previously studied. To our knowledge this is the first report comparing alchemical predictions of proline mutations with experiment. With this strategy in hand, alchemical methods now have access to the full palette of amino acid mutations.
Collapse
Affiliation(s)
- Ryan L Hayes
- Department of Chemistry, University of Michigan, Ann Arbor, Michigan, USA
| | - Charles L Brooks
- Department of Chemistry, University of Michigan, Ann Arbor, Michigan, USA.,Biophysics Program, University of Michigan, Ann Arbor, Michigan, USA
| |
Collapse
|
43
|
Allison JR. Computational methods for exploring protein conformations. Biochem Soc Trans 2020; 48:1707-24. [PMID: 32756904 DOI: 10.1042/BST20200193] [Citation(s) in RCA: 30] [Impact Index Per Article: 10.0] [Reference Citation Analysis] [What about the content of this article? (0)] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/02/2020] [Revised: 07/07/2020] [Accepted: 07/09/2020] [Indexed: 12/13/2022]
Abstract
Proteins are dynamic molecules that can transition between a potentially wide range of structures comprising their conformational ensemble. The nature of these conformations and their relative probabilities are described by a high-dimensional free energy landscape. While computer simulation techniques such as molecular dynamics simulations allow characterisation of the metastable conformational states and the transitions between them, and thus free energy landscapes, to be characterised, the barriers between states can be high, precluding efficient sampling without substantial computational resources. Over the past decades, a dizzying array of methods have emerged for enhancing conformational sampling, and for projecting the free energy landscape onto a reduced set of dimensions that allow conformational states to be distinguished, known as collective variables (CVs), along which sampling may be directed. Here, a brief description of what biomolecular simulation entails is followed by a more detailed exposition of the nature of CVs and methods for determining these, and, lastly, an overview of the myriad different approaches for enhancing conformational sampling, most of which rely upon CVs, including new advances in both CV determination and conformational sampling due to machine learning.
Collapse
|
44
|
Higo J, Takashima H, Fukunishi Y, Yoshimori A. Generalized-ensemble method study: A helix-mimetic compound inhibits protein-protein interaction by long-range and short-range intermolecular interactions. J Comput Chem 2021; 42:956-969. [PMID: 33755222 DOI: 10.1002/jcc.26516] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/18/2020] [Revised: 02/25/2021] [Accepted: 03/01/2021] [Indexed: 12/30/2022]
Abstract
A heterocyclic compound mS-11 is a helix-mimetic designed to inhibit binding of an intrinsic disordered protein neural restrictive silence factor/repressor element 1 silencing factor (NRSF/REST) to a receptor protein mSin3B. We apply a generalized ensemble method, multi-dimensional virtual-system coupled molecular dynamics developed by ourselves recently, to a system consisting of mS-11 and mSin3B, and obtain a thermally equilibrated distribution, which is comprised of the bound and unbound states extensively. The lowest free-energy position of mS-11 coincides with the NRSF/REST position in the experimentally-determined NRSF/REST-mSin3B complex. Importantly, the molecular orientation of mS-11 is ordering in a wide region around mSin3B. The resultant binding scenario is: When mS-11 is distant from the binding site of mSin3B, mS-11 descends the free-energy slope toward the binding site maintaining the molecular orientation to be advantageous for binding. Then, finally a long and flexible hydrophobic sidechain of mS-11 fits into the binding site, which is the lowest-free-energy complex structure inhibiting NRSF/REST binding to mSin3B.
Collapse
Affiliation(s)
- Junichi Higo
- Graduate School of Simulation Studies, University of Hyogo, Kobe, Japan
| | - Hajime Takashima
- Department of Research and Development, PRISM BioLab Co., Ltd., Fujisawa, Japan
| | - Yoshifumi Fukunishi
- Cellular and Molecular Biotechnology Research Institute, National Institute of Advanced Industrial Science and Technology (AIST), Tokyo, Japan
| | - Atsushi Yoshimori
- Chemoinformatics & AI Research Group, Institute for Theoretical Medicine, Inc., Fujisawa, Japan
| |
Collapse
|
45
|
Nierzwicki Ł, Palermo G. Corrigendum: Molecular Dynamics to Predict Cryo-EM: Capturing Transitions and Short-Lived Conformational States of Biomolecules. Front Mol Biosci 2021; 8:698735. [PMID: 34124167 PMCID: PMC8194274 DOI: 10.3389/fmolb.2021.698735] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Key Words] [Download PDF] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/22/2021] [Accepted: 04/27/2021] [Indexed: 11/13/2022] Open
Affiliation(s)
- Łukasz Nierzwicki
- Department of Bioengineering, University of California, Riverside, CA, United States
| | - Giulia Palermo
- Department of Bioengineering, University of California, Riverside, CA, United States.,Department of Chemistry, University of California, Riverside, CA, United States
| |
Collapse
|
46
|
Nierzwicki Ł, Palermo G. Molecular Dynamics to Predict Cryo-EM: Capturing Transitions and Short-Lived Conformational States of Biomolecules. Front Mol Biosci 2021; 8:641208. [PMID: 33884260 PMCID: PMC8053777 DOI: 10.3389/fmolb.2021.641208] [Citation(s) in RCA: 10] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/13/2020] [Accepted: 02/15/2021] [Indexed: 12/21/2022] Open
Abstract
Single-particle cryogenic electron microscopy (cryo-EM) has revolutionized the field of the structural biology, providing an access to the atomic resolution structures of large biomolecular complexes in their near-native environment. Today's cryo-EM maps can frequently reach the atomic-level resolution, while often containing a range of resolutions, with conformationally variable regions obtained at 6 Å or worse. Low resolution density maps obtained for protein flexible domains, as well as the ensemble of coexisting conformational states arising from cryo-EM, poses new challenges and opportunities for Molecular Dynamics (MD) simulations. With the ability to describe the biomolecular dynamics at the atomic level, MD can extend the capabilities of cryo-EM, capturing the conformational variability and predicting biologically relevant short-lived conformational states. Here, we report about the state-of-the-art MD procedures that are currently used to refine, reconstruct and interpret cryo-EM maps. We show the capability of MD to predict short-lived conformational states, finding remarkable confirmation by cryo-EM structures subsequently solved. This has been the case of the CRISPR-Cas9 genome editing machinery, whose catalytically active structure has been predicted through both long-time scale MD and enhanced sampling techniques 2 years earlier than cryo-EM. In summary, this contribution remarks the ability of MD to complement cryo-EM, describing conformational landscapes and relating structural transitions to function, ultimately discerning relevant short-lived conformational states and providing mechanistic knowledge of biological function.
Collapse
Affiliation(s)
- Łukasz Nierzwicki
- Department of Bioengineering, University of California, Riverside, CA, United States
| | - Giulia Palermo
- Department of Bioengineering, University of California, Riverside, CA, United States
- Department of Chemistry, University of California, Riverside, CA, United States
| |
Collapse
|
47
|
Namsani S, Pramanik D, Khan MA, Roy S, Singh JK. Metadynamics-based enhanced sampling protocol for virtual screening: case study for 3CLpro protein for SARS-CoV-2. J Biomol Struct Dyn 2021; 40:7002-7017. [PMID: 33663346 DOI: 10.1080/07391102.2021.1892530] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/22/2022]
Abstract
In recent times, computational methods played an important role in the down selection of chemical compounds, which could be a potential drug candidate with a high affinity to target proteins. However, the screening methodologies, including docking, often fails to identify the most effective compound, which could be a ligand for the target protein. To solve that, here we have integrated meta-dynamics, an enhanced sampling molecular simulation method, with all-atom molecular dynamics to determine a specific compound that could target the main protease of novel severe acute respiratory syndrome coronavirus 2 (SARS-COV-2). This combined computational approach uses the enhanced sampling to explore the free energy surface associated with the protein's binding site (including the ligand) in an explicit solvent. We have implemented this method to find new chemical entities that exhibit high specificity of binding to the 3-chymotrypsin-like cysteine protease (3CLpro) present in the SARS-CoV-2 and segregated to the most strongly bound ligands based on free energy and scoring functions (defined and implemented) from a set of 17 ligands which were prescreened for synthesizability and druggability. Additionally, we have compared these 17 ligands' affinities against controls, N3 and 13b α-ketoamide inhibitors, for which experimental crystal structures are available. Based on our results and analysis from the combined molecular simulation approach, we could identify the best compound which could be further taken as a potential candidate for experimental validation.Communicated by Ramaswamy H. Sarma.
Collapse
Affiliation(s)
| | - Debabrata Pramanik
- Department of Chemical Engineering, Indian Institute of Technology, Kanpur, India
| | - Mohd Aamir Khan
- Prescience Insilico Private Limited, Bangalore, India.,Department of Chemical Engineering, Indian Institute of Technology, Kanpur, India
| | - Sudip Roy
- Prescience Insilico Private Limited, Bangalore, India
| | - Jayant Kumar Singh
- Prescience Insilico Private Limited, Bangalore, India.,Department of Chemical Engineering, Indian Institute of Technology, Kanpur, India
| |
Collapse
|
48
|
Fleetwood O, Carlsson J, Delemotte L. Identification of ligand-specific G protein-coupled receptor states and prediction of downstream efficacy via data-driven modeling. eLife 2021; 10:60715. [PMID: 33506760 PMCID: PMC7886328 DOI: 10.7554/elife.60715] [Citation(s) in RCA: 14] [Impact Index Per Article: 4.7] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/03/2020] [Accepted: 01/27/2021] [Indexed: 12/22/2022] Open
Abstract
Ligand binding stabilizes different G protein-coupled receptor states via a complex allosteric process that is not completely understood. Here, we have derived free energy landscapes describing activation of the β2 adrenergic receptor bound to ligands with different efficacy profiles using enhanced sampling molecular dynamics simulations. These reveal shifts toward active-like states at the Gprotein-binding site for receptors bound to partial and full agonists, and that the ligands modulate the conformational ensemble of the receptor by tuning protein microswitches. We indeed find an excellent correlation between the conformation of the microswitches close to the ligand binding site and in the transmembrane region and experimentally reported cyclic adenosine monophosphate signaling responses. Dimensionality reduction further reveals the similarity between the unique conformational states induced by different ligands, and examining the output of classifiers highlights two distant hotspots governing agonism on transmembrane helices 5 and 7.
Collapse
Affiliation(s)
- Oliver Fleetwood
- Science for Life Laboratory, Department of Applied Physics, KTH Royal Institute of Technology, Stockholm, Sweden
| | - Jens Carlsson
- Science for Life Laboratory, Department of Cell and Molecular Biology, Uppsala University, Uppsala, Sweden
| | - Lucie Delemotte
- Science for Life Laboratory, Department of Applied Physics, KTH Royal Institute of Technology, Stockholm, Sweden
| |
Collapse
|
49
|
Abstract
Recent advances in theory and algorithms for atomically detailed simulations open the way to the study of the kinetics of a wide range of molecular processes in biophysics. The theories propose a shift from the traditionally very long molecular dynamic trajectories, which are exact but may not be efficient in the study of kinetics, to the use of a large number of short trajectories. The short trajectories exploit a mapping to a mesh in coarse space and allow for efficient calculations of kinetics and thermodynamics. In this review, I focus on one theory: Milestoning is a theory and an algorithm that offers a hierarchical calculation of properties of interest, such as the free energy profile and the mean first passage time. Approximations to the true long-time dynamics can be computed efficiently and assessed at different steps of the investigation. The theory is discussed and illustrated using two biophysical examples: ion permeation through a phospholipid membrane and protein translocation through a channel.
Collapse
Affiliation(s)
- Ron Elber
- Oden Institute for Computational Engineering and Sciences, Department of Chemistry, The University of Texas at Austin, Austin, Texas 78712, USA;
| |
Collapse
|
50
|
Abstract
Fentanyl is a potent opioid analgesic, which for decades has been used routinely in surgical and therapeutic applications. In addition to its analgesic properties, fentanyl also possesses anesthetic properties, which are not well understood. Fentanyl is used in the general anesthesia process to induce and maintain anesthesia in combination with the general anesthetic propofol, which fentanyl is known to potentiate. As the atomic-level mechanism behind the potentiation of propofol is unclear, we have used classical molecular dynamics simulations to study the interactions of these drugs with the Gloeobacter violaceus ion channel (GLIC). This ion channel has been identified as a target for many anesthetic drugs. We identified multiple binding sites using flooding style and Gaussian accelerated molecular dynamics (GaMD) simulations, showing fentanyl acting as a stabiliser that holds propofol within binding sites. Our extensive GaMD simulations were also able to show the pathway by which propofol blocks the channel pore, which has previously been suggested as a mechanism for ion channel modulation. General anesthesia is a multi-drug process and this study provides the first insight into the interactions between two different drugs in the anesthesia process in a relevant biological environment.Communicated by Ramaswamy H. Sarma.
Collapse
Affiliation(s)
| | - Nora H de Leeuw
- School of Chemistry, Cardiff University, Cardiff, UK.,School of Chemistry, University of Leeds, Leeds, UK
| |
Collapse
|