1
|
Zia SR, Coricello A, Bottegoni G. Increased throughput in methods for simulating protein ligand binding and unbinding. Curr Opin Struct Biol 2024; 87:102871. [PMID: 38924980 DOI: 10.1016/j.sbi.2024.102871] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/03/2024] [Revised: 06/03/2024] [Accepted: 06/04/2024] [Indexed: 06/28/2024]
Abstract
By incorporating full flexibility and enabling the quantification of crucial parameters such as binding free energies and residence times, methods for investigating protein-ligand binding and unbinding via molecular dynamics provide details on the involved mechanisms at the molecular level. While these advancements hold promise for impacting drug discovery, a notable drawback persists: their relatively time-consuming nature limits throughput. Herein, we survey recent implementations which, employing a blend of enhanced sampling techniques, a clever choice of collective variables, and often machine learning, strive to enhance the efficiency of new and previously reported methods without compromising accuracy. Particularly noteworthy is the validation of these methods that was often performed on systems mirroring real-world drug discovery scenarios.
Collapse
Affiliation(s)
- Syeda Rehana Zia
- Department of Paediatrics and Child Health, Faculty of Health Sciences, Medical College, The Aga Khan University, Karachi, 74800, Pakistan
| | - Adriana Coricello
- Department of Biomolecular Sciences, University of Urbino Carlo Bo, Urbino, 61029, Italy.
| | - Giovanni Bottegoni
- Department of Biomolecular Sciences, University of Urbino Carlo Bo, Urbino, 61029, Italy; Institute of Clinical Sciences, College of Medical and Dental Sciences, University of Birmingham, B15 2TT, United Kingdom.
| |
Collapse
|
2
|
Akhter S, Tang Z, Wang J, Haboro M, Holmstrom ED, Wang J, Miao Y. Mechanism of Ligand Binding to Theophylline RNA Aptamer. J Chem Inf Model 2024; 64:1017-1029. [PMID: 38226603 PMCID: PMC11058067 DOI: 10.1021/acs.jcim.3c01454] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/17/2024]
Abstract
Studying RNA-ligand interactions and quantifying their binding thermodynamics and kinetics are of particular relevance in the field of drug discovery. Here, we combined biochemical binding assays and accelerated molecular simulations to investigate ligand binding and dissociation in RNA using the theophylline-binding RNA as a model system. All-atom simulations using a Ligand Gaussian accelerated Molecular Dynamics method (LiGaMD) have captured repetitive binding and dissociation of theophylline and caffeine to RNA. Theophylline's binding free energy and kinetic rate constants align with our experimental data, while caffeine's binding affinity is over 10,000 times weaker, and its kinetics could not be determined. LiGaMD simulations allowed us to identify distinct low-energy conformations and multiple ligand binding pathways to RNA. Simulations revealed a "conformational selection" mechanism for ligand binding to the flexible RNA aptamer, which provides important mechanistic insights into ligand binding to the theophylline-binding model. Our findings suggest that compound docking using a structural ensemble of representative RNA conformations would be necessary for structure-based drug design of flexible RNA.
Collapse
Affiliation(s)
- Sana Akhter
- Computational Biology Program and Department of Molecular Biosciences, University of Kansas, Lawrence, Kansas 66047, United States
| | - Zhichao Tang
- Department of Medicinal Chemistry, University of Kansas, Lawrence, Kansas 66047, United States
| | - Jinan Wang
- Computational Biology Program and Department of Molecular Biosciences, University of Kansas, Lawrence, Kansas 66047, United States
| | - Mercy Haboro
- Department of Medicinal Chemistry, University of Kansas, Lawrence, Kansas 66047, United States
| | - Erik D Holmstrom
- Department of Molecular Biosciences and Department of Chemistry, University of Kansas, Lawrence, Kansas 66045, United States
| | - Jingxin Wang
- Department of Medicinal Chemistry, University of Kansas, Lawrence, Kansas 66047, United States
| | - Yinglong Miao
- Computational Biology Program and Department of Molecular Biosciences, University of Kansas, Lawrence, Kansas 66047, United States
| |
Collapse
|
3
|
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] [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
|
4
|
Do HN, Miao Y. Deep Boosted Molecular Dynamics (DBMD): Accelerating molecular simulations with Gaussian boost potentials generated using probabilistic Bayesian deep neural network. BIORXIV : THE PREPRINT SERVER FOR BIOLOGY 2023:2023.03.25.534210. [PMID: 37034713 PMCID: PMC10081221 DOI: 10.1101/2023.03.25.534210] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Indexed: 04/22/2023]
Abstract
We have developed a new Deep Boosted Molecular Dynamics (DBMD) method. Probabilistic Bayesian neural network models were implemented to construct boost potentials that exhibit Gaussian distribution with minimized anharmonicity, thereby allowing for accurate energetic reweighting and enhanced sampling of molecular simulations. DBMD was demonstrated on model systems of alanine dipeptide and the fast-folding protein and RNA structures. For alanine dipeptide, 30ns DBMD simulations captured up to 83-125 times more backbone dihedral transitions than 1μs conventional molecular dynamics (cMD) simulations and were able to accurately reproduce the original free energy profiles. Moreover, DBMD sampled multiple folding and unfolding events within 300ns simulations of the chignolin model protein and identified low-energy conformational states comparable to previous simulation findings. Finally, DBMD captured a general folding pathway of three hairpin RNAs with the GCAA, GAAA, and UUCG tetraloops. Based on Deep Learning neural network, DBMD provides a powerful and generally applicable approach to boosting biomolecular simulations. DBMD is available with open source in OpenMM at https://github.com/MiaoLab20/DBMD/.
Collapse
Affiliation(s)
- Hung N. Do
- Center for Computational Biology and Department of Molecular Biosciences, University of Kansas, Lawrence, Kansas 66047
| | - Yinglong Miao
- Center for Computational Biology and Department of Molecular Biosciences, University of Kansas, Lawrence, Kansas 66047
- To whom correspondence should be addressed:
| |
Collapse
|
5
|
Csizi K, Reiher M. Universal
QM
/
MM
approaches for general nanoscale applications. WIRES COMPUTATIONAL MOLECULAR SCIENCE 2023. [DOI: 10.1002/wcms.1656] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 02/04/2023]
Affiliation(s)
| | - Markus Reiher
- Laboratorium für Physikalische Chemie ETH Zürich Zürich Switzerland
| |
Collapse
|
6
|
Tvaroška I, Kozmon S, Kóňa J. Molecular Modeling Insights into the Structure and Behavior of Integrins: A Review. Cells 2023; 12:cells12020324. [PMID: 36672259 PMCID: PMC9856412 DOI: 10.3390/cells12020324] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/07/2022] [Revised: 01/09/2023] [Accepted: 01/11/2023] [Indexed: 01/18/2023] Open
Abstract
Integrins are heterodimeric glycoproteins crucial to the physiology and pathology of many biological functions. As adhesion molecules, they mediate immune cell trafficking, migration, and immunological synapse formation during inflammation and cancer. The recognition of the vital roles of integrins in various diseases revealed their therapeutic potential. Despite the great effort in the last thirty years, up to now, only seven integrin-based drugs have entered the market. Recent progress in deciphering integrin functions, signaling, and interactions with ligands, along with advancement in rational drug design strategies, provide an opportunity to exploit their therapeutic potential and discover novel agents. This review will discuss the molecular modeling methods used in determining integrins' dynamic properties and in providing information toward understanding their properties and function at the atomic level. Then, we will survey the relevant contributions and the current understanding of integrin structure, activation, the binding of essential ligands, and the role of molecular modeling methods in the rational design of antagonists. We will emphasize the role played by molecular modeling methods in progress in these areas and the designing of integrin antagonists.
Collapse
Affiliation(s)
- Igor Tvaroška
- Institute of Chemistry, Slovak Academy of Sciences, Dúbravska cesta 9, 845 38 Bratislava, Slovakia
- Correspondence:
| | - Stanislav Kozmon
- Institute of Chemistry, Slovak Academy of Sciences, Dúbravska cesta 9, 845 38 Bratislava, Slovakia
- Medical Vision o. z., Záhradnícka 4837/55, 821 08 Bratislava, Slovakia
| | - Juraj Kóňa
- Institute of Chemistry, Slovak Academy of Sciences, Dúbravska cesta 9, 845 38 Bratislava, Slovakia
- Medical Vision o. z., Záhradnícka 4837/55, 821 08 Bratislava, Slovakia
| |
Collapse
|
7
|
Applications of molecular dynamics simulation in nanomedicine. Nanomedicine (Lond) 2023. [DOI: 10.1016/b978-0-12-818627-5.00007-5] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/19/2023] Open
|
8
|
Nakkala S, Modak C, Bathula R, Lanka G, Somadi G, Sreekanth S, Jain A, Potlapally SR. Identification of new anti-cancer agents against CENTERIN: Structure-based virtual screening, AutoDock and binding free energy studies. J Mol Struct 2022. [DOI: 10.1016/j.molstruc.2022.133952] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/16/2022]
|
9
|
Kaur H, Siwal SS, Saini RV, Singh N, Thakur VK. Significance of an Electrochemical Sensor and Nanocomposites: Toward the Electrocatalytic Detection of Neurotransmitters and Their Importance within the Physiological System. ACS NANOSCIENCE AU 2022; 3:1-27. [PMID: 37101467 PMCID: PMC10125382 DOI: 10.1021/acsnanoscienceau.2c00039] [Citation(s) in RCA: 8] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 08/15/2022] [Revised: 10/17/2022] [Accepted: 10/17/2022] [Indexed: 11/07/2022]
Abstract
A prominent neurotransmitter (NT), dopamine (DA), is a chemical messenger that transmits signals between one neuron to the next to pass on a signal to and from the central nervous system (CNS). The imbalanced concentration of DA may cause numerous neurological sicknesses and syndromes, for example, Parkinson's disease (PD) and schizophrenia. There are many types of NTs in the brain, including epinephrine, norepinephrine (NE), serotonin, and glutamate. Electrochemical sensors have offered a creative direction to biomedical analysis and testing. Researches are in progress to improve the performance of sensors and develop new protocols for sensor design. This review article focuses on the area of sensor growth to discover the applicability of polymers and metallic particles and composite materials as tools in electrochemical sensor surface incorporation. Electrochemical sensors have attracted the attention of researchers as they possess high sensitivity, quick reaction rate, good controllability, and instantaneous detection. Efficient complex materials provide considerable benefits for biological detection as they have exclusive chemical and physical properties. Due to distinctive electrocatalytic characteristics, metallic nanoparticles add fascinating traits to materials that depend on the material's morphology and size. Herein, we have collected much information on NTs and their importance within the physiological system. Furthermore, the electrochemical sensors and corresponding techniques (such as voltammetric, amperometry, impedance, and chronoamperometry) and the different types of electrodes' roles in the analysis of NTs are discussed. Furthermore, other methods for detecting NTs include optical and microdialysis methods. Finally, we show the advantages and disadvantages of different techniques and conclude remarks with future perspectives.
Collapse
Affiliation(s)
- Harjot Kaur
- Department of Chemistry, M.M. Engineering College, Maharishi Markandeshwar (Deemed to be University), Mullana-Ambala, Haryana 133207, India
| | - Samarjeet Singh Siwal
- Department of Chemistry, M.M. Engineering College, Maharishi Markandeshwar (Deemed to be University), Mullana-Ambala, Haryana 133207, India
| | - Reena V. Saini
- Department of Biotechnology, Maharishi Markandeshwar (Deemed to be University), Mullana-Ambala, Haryana 133207, India
| | - Nirankar Singh
- Department of Chemistry, M.M. Engineering College, Maharishi Markandeshwar (Deemed to be University), Mullana-Ambala, Haryana 133207, India
| | - Vijay Kumar Thakur
- Biorefining and Advanced Materials Research Center, Scotland’s Rural College (SRUC), Kings Buildings, Edinburgh EH9 3JG, United Kingdom
- School of Engineering, University of Petroleum & Energy Studies (UPES), Dehradun, Uttarakhand 248007, India
- Centre for Research & Development, Chandigarh University, Mohali, Punjab 140413, India
| |
Collapse
|
10
|
Lenard AJ, Mulder FAA, Madl T. Solvent paramagnetic relaxation enhancement as a versatile method for studying structure and dynamics of biomolecular systems. PROGRESS IN NUCLEAR MAGNETIC RESONANCE SPECTROSCOPY 2022; 132-133:113-139. [PMID: 36496256 DOI: 10.1016/j.pnmrs.2022.09.001] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 02/12/2022] [Revised: 09/14/2022] [Accepted: 09/15/2022] [Indexed: 06/17/2023]
Abstract
Solvent paramagnetic relaxation enhancement (sPRE) is a versatile nuclear magnetic resonance (NMR)-based method that allows characterization of the structure and dynamics of biomolecular systems through providing quantitative experimental information on solvent accessibility of NMR-active nuclei. Addition of soluble paramagnetic probes to the solution of a biomolecule leads to paramagnetic relaxation enhancement in a concentration-dependent manner. Here we review recent progress in the sPRE-based characterization of structural and dynamic properties of biomolecules and their complexes, and aim to deliver a comprehensive illustration of a growing number of applications of the method to various biological systems. We discuss the physical principles of sPRE measurements and provide an overview of available co-solute paramagnetic probes. We then explore how sPRE, in combination with complementary biophysical techniques, can further advance biomolecular structure determination, identification of interaction surfaces within protein complexes, and probing of conformational changes and low-population transient states, as well as deliver insights into weak, nonspecific, and transient interactions between proteins and co-solutes. In addition, we present examples of how the incorporation of solvent paramagnetic probes can improve the sensitivity of NMR experiments and discuss the prospects of applying sPRE to NMR metabolomics, drug discovery, and the study of intrinsically disordered proteins.
Collapse
Affiliation(s)
- Aneta J Lenard
- Gottfried Schatz Research Center for Cell Signaling, Metabolism and Ageing, Molecular Biology and Biochemistry, Research Unit Integrative Structural Biology, Medical University of Graz, 8010 Graz, Austria.
| | - Frans A A Mulder
- Interdisciplinary Nanoscience Center and Department of Chemistry, University of Aarhus, DK-8000 Aarhus, Denmark; Institute of Biochemistry, Johannes Kepler Universität Linz, 4040 Linz, Austria.
| | - Tobias Madl
- Gottfried Schatz Research Center for Cell Signaling, Metabolism and Ageing, Molecular Biology and Biochemistry, Research Unit Integrative Structural Biology, Medical University of Graz, 8010 Graz, Austria; BioTechMed-Graz, 8010 Graz, Austria.
| |
Collapse
|
11
|
Dietschreit JCB, Diestler DJ, Hulm A, Ochsenfeld C, Gómez-Bombarelli R. From free-energy profiles to activation free energies. J Chem Phys 2022; 157:084113. [DOI: 10.1063/5.0102075] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Given a chemical reaction going from reactant (R) to the product (P) on a potential energy surface (PES) and a collective variable (CV) discriminating between R and P, we define the free-energy profile (FEP) as the logarithm of the marginal Boltzmann distribution of the CV. This FEP is not a true free energy. Nevertheless, it is common to treat the FEP as the “free-energy” analog of the minimum potential energy path and to take the activation free energy, [Formula: see text], as the difference between the maximum at the transition state and the minimum at R. We show that this approximation can result in large errors. The FEP depends on the CV and is, therefore, not unique. For the same reaction, different discriminating CVs can yield different [Formula: see text]. We derive an exact expression for the activation free energy that avoids this ambiguity. We find [Formula: see text] to be a combination of the probability of the system being in the reactant state, the probability density on the dividing surface, and the thermal de Broglie wavelength associated with the transition. We apply our formalism to simple analytic models and realistic chemical systems and show that the FEP-based approximation applies only at low temperatures for CVs with a small effective mass. Most chemical reactions occur on complex, high-dimensional PES that cannot be treated analytically and pose the added challenge of choosing a good CV. We study the influence of that choice and find that, while the reaction free energy is largely unaffected, [Formula: see text] is quite sensitive.
Collapse
Affiliation(s)
- Johannes C. B. Dietschreit
- Department of Materials Science and Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, USA
| | | | - Andreas Hulm
- Chair of Theoretical Chemistry, Department of Chemistry, University of Munich (LMU), Butenandtstr. 7, D-81377 München, Germany
| | - Christian Ochsenfeld
- Chair of Theoretical Chemistry, Department of Chemistry, University of Munich (LMU), Butenandtstr. 7, D-81377 München, Germany
- Max Planck Institute for Solid State Research, Heisenbergstr. 1, D-70569 Stuttgart, Germany
| | - Rafael Gómez-Bombarelli
- Department of Materials Science and Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, USA
| |
Collapse
|
12
|
Copeland M, Do HN, Votapka L, Joshi K, Wang J, Amaro RE, Miao Y. Gaussian Accelerated Molecular Dynamics in OpenMM. J Phys Chem B 2022; 126:5810-5820. [PMID: 35895977 PMCID: PMC9773147 DOI: 10.1021/acs.jpcb.2c03765] [Citation(s) in RCA: 6] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/24/2022]
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
|
13
|
Jakhmola S, Sk MF, Chatterjee A, Jain K, Kar P, Jha HC. A plausible contributor to multiple sclerosis; presentation of antigenic myelin protein epitopes by major histocompatibility complexes. Comput Biol Med 2022; 148:105856. [PMID: 35863244 DOI: 10.1016/j.compbiomed.2022.105856] [Citation(s) in RCA: 7] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/05/2022] [Revised: 06/17/2022] [Accepted: 06/17/2022] [Indexed: 11/17/2022]
Abstract
BACKGROUND Multiple sclerosis (MS) can be induced upon successful presentation of myelin antigens by MHC I/II. Antigenic similarity between the myelin and viral proteins may worsen the immunological responses. METHODOLOGY Antigenic regions within myelin proteins; PLP1, MBP, MOG, and MAG were analyzed using SVMTrip and EMBOSS. Homology search identified sequence similarity between the predicted host epitopes and viral proteins. NetMHCpan predicted MHC I/II binding followed by peptide-protein docking through the HPEPDOCK server. Thereafter we analyzed conformational flexibility and stability of 15 protein-peptide complexes based on high docking scores. The binding free energy was calculated using conventional (MD) and Gaussian accelerated molecular dynamics simulation. RESULTS PLP1, MBP, MAG and MOG contained numerous antigenic epitopes. MBP and MOG epitopes had sequence similarity to HHV-6 BALF5; EBNA1 and CMV glycoprotein M (gM), and EBV LMP2B, gp350/220; HHV-8 ORFs respectively. Many herpes virus proteins like tegument, envelope glycoproteins, and ORFs of EBV, CMV, HHV-6, and HHV-8 demonstrated sequence similarity with MAG and PLP1. Some antigenic peptides were also linear B-cell epitopes and influenced cytokine production by T-cell. MHC I allele HLA-B*57:01 bound to PLP1 peptide and HLA-A*68:02 bound to a MAG peptide strongly. MHC II alleles HLA-DRB1*04:05 and HLA-DR1*01:01 associated with MAG- and MOG-derived peptides, respectively, demonstrating high HPEPDOCK scores. MD simulations established stable binding of certain peptides with the MHC namely HLA-B*51:01-MBP(DYKSAHKGFKGVDAQGTLSKIFKL), HLA-B*57:01-PLP1(PDKFVGITYALTVVWLLVFACSAVPVYIYF), HLA-DR1*01:01-MOG(VEDPFYWVSPGVLVLLAVLPVLLLQITVGLVFLCLQYR) and HLA-DRB1*04:05-MAG(TWVQVSLLHFVPTREA). CONCLUSIONS Cross-reactivity between self-antigens and pathogen derived immunodominant epitopes may induce MS. Our study supported the role of specific MHC alleles as a contributing MS risk factor.
Collapse
Affiliation(s)
- Shweta Jakhmola
- Infection Bioengineering Group, Department of Biosciences and Biomedical Engineering, Indian Institute of Technology Indore, India.
| | - Md Fulbabu Sk
- Department of Biosciences and Biomedical Engineering, Indian Institute of Technology Indore, India
| | - Akash Chatterjee
- Infection Bioengineering Group, Department of Biosciences and Biomedical Engineering, Indian Institute of Technology Indore, India
| | - Khushboo Jain
- Infection Bioengineering Group, Department of Biosciences and Biomedical Engineering, Indian Institute of Technology Indore, India
| | - Parimal Kar
- Department of Biosciences and Biomedical Engineering, Indian Institute of Technology Indore, India.
| | - Hem Chandra Jha
- Infection Bioengineering Group, Department of Biosciences and Biomedical Engineering, Indian Institute of Technology Indore, India.
| |
Collapse
|
14
|
Hulm A, Dietschreit JCB, Ochsenfeld C. Statistically optimal analysis of the extended-system adaptive biasing force (eABF) method. J Chem Phys 2022; 157:024110. [PMID: 35840392 DOI: 10.1063/5.0095554] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/13/2023] Open
Abstract
The extended-system adaptive biasing force (eABF) method and its newer variants offer rapid exploration of the configuration space of chemical systems. Instead of directly applying the ABF bias to collective variables, they are harmonically coupled to fictitious particles, which separates the problem of enhanced sampling from that of free energy estimation. The prevalent analysis method to obtain the potential of mean force (PMF) from eABF is thermodynamic integration. However, besides the PMF, most information is lost as the unbiased probability of visited configurations is never recovered. In this contribution, we show how statistical weights of individual frames can be computed using the Multistate Bennett's Acceptance Ratio (MBAR), putting the post-processing of eABF on one level with other frequently used sampling methods. In addition, we apply this formalism to the prediction of nuclear magnetic resonance shieldings, which are very sensitive to molecular geometries and often require extensive sampling. The results show that the combination of enhanced sampling by means of extended-system dynamics with the MBAR estimator is a highly useful tool for the calculation of ensemble properties. Furthermore, the extension of the presented scheme to the recently published Gaussian-accelerated molecular dynamics eABF hybrid is straightforward and approximation free.
Collapse
Affiliation(s)
- Andreas Hulm
- Chair of Theoretical Chemistry, Department of Chemistry, University of Munich (LMU), Butenandtstr. 7, D-81377 München, Germany
| | - Johannes C B Dietschreit
- Chair of Theoretical Chemistry, Department of Chemistry, University of Munich (LMU), Butenandtstr. 7, D-81377 München, Germany
| | - Christian Ochsenfeld
- Chair of Theoretical Chemistry, Department of Chemistry, University of Munich (LMU), Butenandtstr. 7, D-81377 München, Germany
| |
Collapse
|
15
|
Monroe JI, Shen VK. Learning Efficient, Collective Monte Carlo Moves with Variational Autoencoders. J Chem Theory Comput 2022; 18:3622-3636. [PMID: 35613327 PMCID: PMC11210279 DOI: 10.1021/acs.jctc.2c00110] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Discovering meaningful collective variables for enhancing sampling, via applied biasing potentials or tailored MC move sets, remains a major challenge within molecular simulation. While recent studies identifying collective variables with variational autoencoders (VAEs) have focused on the encoding and latent space discovered by a VAE, the impact of the decoding and its ability to act as a generative model remains unexplored. We demonstrate how VAEs may be used to learn (on-the-fly and with minimal human intervention) highly efficient, collective Monte Carlo moves that accelerate sampling along the learned collective variable. In contrast to many machine learning-based efforts to bias sampling and generate novel configurations, our methods result in exact sampling in the ensemble of interest and do not require reweighting. In fact, we show that the acceptance rates of our moves approach unity for a perfect VAE model. While this is never observed in practice, VAE-based Monte Carlo moves still enhance sampling of new configurations. We demonstrate, however, that the form of the encoding and decoding distributions, in particular the extent to which the decoder reflects the underlying physics, greatly impacts the performance of the trained VAE.
Collapse
Affiliation(s)
- Jacob I Monroe
- Chemical Sciences Division, National Institute of Standards and Technology, Gaithersburg, Maryland 20899-8320, United States
| | - Vincent K Shen
- Chemical Sciences Division, National Institute of Standards and Technology, Gaithersburg, Maryland 20899-8320, United States
| |
Collapse
|
16
|
Spiwok V, Kurečka M, Křenek A. Collective Variable for Metadynamics Derived From AlphaFold Output. Front Mol Biosci 2022; 9:878133. [PMID: 35769910 PMCID: PMC9234394 DOI: 10.3389/fmolb.2022.878133] [Citation(s) in RCA: 4] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/17/2022] [Accepted: 05/05/2022] [Indexed: 11/13/2022] Open
Abstract
AlphaFold is a neural network–based tool for the prediction of 3D structures of proteins. In CASP14, a blind structure prediction challenge, it performed significantly better than other competitors, making it the best available structure prediction tool. One of the outputs of AlphaFold is the probability profile of residue–residue distances. This makes it possible to score any conformation of the studied protein to express its compliance with the AlphaFold model. Here, we show how this score can be used to drive protein folding simulation by metadynamics and parallel tempering metadynamics. Using parallel tempering metadynamics, we simulated the folding of a mini-protein Trp-cage and β hairpin and predicted their folding equilibria. We observe the potential of the AlphaFold-based collective variable in applications beyond structure prediction, such as in structure refinement or prediction of the outcome of a mutation.
Collapse
Affiliation(s)
- Vojtěch Spiwok
- Department of Biochemistry and Microbiology, Faculty of Food and Biochemical Technology, University of Chemistry and Technology, Prague, Czechia
- *Correspondence: Vojtěch Spiwok,
| | - Martin Kurečka
- Institute of Computer Science, Masaryk University, Brno, Czechia
| | - Aleš Křenek
- Institute of Computer Science, Masaryk University, Brno, Czechia
| |
Collapse
|
17
|
Wang J, Bhattarai A, Do HN, Akhter S, Miao Y. Molecular Simulations and Drug Discovery of Adenosine Receptors. Molecules 2022; 27:2054. [PMID: 35408454 PMCID: PMC9000248 DOI: 10.3390/molecules27072054] [Citation(s) in RCA: 4] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/31/2022] [Revised: 03/18/2022] [Accepted: 03/20/2022] [Indexed: 02/02/2023] Open
Abstract
G protein-coupled receptors (GPCRs) represent the largest family of human membrane proteins. Four subtypes of adenosine receptors (ARs), the A1AR, A2AAR, A2BAR and A3AR, each with a unique pharmacological profile and distribution within the tissues in the human body, mediate many physiological functions and serve as critical drug targets for treating numerous human diseases including cancer, neuropathic pain, cardiac ischemia, stroke and diabetes. The A1AR and A3AR preferentially couple to the Gi/o proteins, while the A2AAR and A2BAR prefer coupling to the Gs proteins. Adenosine receptors were the first subclass of GPCRs that had experimental structures determined in complex with distinct G proteins. Here, we will review recent studies in molecular simulations and computer-aided drug discovery of the adenosine receptors and also highlight their future research opportunities.
Collapse
Affiliation(s)
| | | | | | | | - Yinglong Miao
- Center for Computational Biology and Department of Molecular Biosciences, University of Kansas, Lawrence, KS 66047, USA; (J.W.); (A.B.); (H.N.D.); (S.A.)
| |
Collapse
|
18
|
Pawnikar S, Miao Y. Mechanism of Peptide Agonist Binding in CXCR4 Chemokine Receptor. Front Mol Biosci 2022; 9:821055. [PMID: 35359589 PMCID: PMC8963245 DOI: 10.3389/fmolb.2022.821055] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/23/2021] [Accepted: 02/14/2022] [Indexed: 01/07/2023] Open
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
|
19
|
Wang J, Miao Y. Protein-Protein Interaction-Gaussian Accelerated Molecular Dynamics (PPI-GaMD): Characterization of Protein Binding Thermodynamics and Kinetics. J Chem Theory Comput 2022; 18:1275-1285. [PMID: 35099970 DOI: 10.1021/acs.jctc.1c00974] [Citation(s) in RCA: 18] [Impact Index Per Article: 9.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/11/2023]
Abstract
Protein-protein interactions (PPIs) play key roles in many fundamental biological processes such as cellular signaling and immune responses. However, it has proven challenging to simulate repetitive protein association and dissociation in order to calculate binding free energies and kinetics of PPIs due to long biological timescales and complex protein dynamics. To address this challenge, we have developed a new computational approach to all-atom simulations of PPIs based on a robust Gaussian accelerated molecular dynamics (GaMD) technique. The method, termed "PPI-GaMD", selectively boosts interaction potential energy between protein partners to facilitate their slow dissociation. Meanwhile, another boost potential is applied to the remaining potential energy of the entire system to effectively model the protein's flexibility and rebinding. PPI-GaMD has been demonstrated on a model system of the ribonuclease barnase interactions with its inhibitor barstar. Six independent 2 μs PPI-GaMD simulations have captured repetitive barstar dissociation and rebinding events, which enable calculations of the protein binding thermodynamics and kinetics simultaneously. The calculated binding free energies and kinetic rate constants agree well with the experimental data. Furthermore, PPI-GaMD simulations have provided mechanistic insights into barstar binding to barnase, which involves long-range electrostatic interactions and multiple binding pathways, being consistent with previous experimental and computational findings of this model system. In summary, PPI-GaMD provides a highly efficient and easy-to-use approach for binding free energy and kinetics calculations of PPIs.
Collapse
Affiliation(s)
- Jinan Wang
- Center for Computational Biology and Department of Molecular Biosciences, University of Kansas, Lawrence, Kansas 66047, United States
| | - Yinglong Miao
- Center for Computational Biology and Department of Molecular Biosciences, University of Kansas, Lawrence, Kansas 66047, United States
| |
Collapse
|
20
|
Do HN, Wang J, Bhattarai A, Miao Y. GLOW: A Workflow Integrating Gaussian-Accelerated Molecular Dynamics and Deep Learning for Free Energy Profiling. J Chem Theory Comput 2022; 18:1423-1436. [PMID: 35200019 PMCID: PMC9773012 DOI: 10.1021/acs.jctc.1c01055] [Citation(s) in RCA: 22] [Impact Index Per Article: 11.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/24/2022]
Abstract
We introduce a Gaussian-accelerated molecular dynamics (GaMD), deep learning (DL), and free energy profiling workflow (GLOW) to predict molecular determinants and map free energy landscapes of biomolecules. All-atom GaMD-enhanced sampling simulations are first performed on biomolecules of interest. Structural contact maps are then calculated from GaMD simulation frames and transformed into images for building DL models using a convolutional neural network. Important structural contacts are further determined from DL models of attention maps of the structural contact gradients, which allow us to identify the system reaction coordinates. Finally, free energy profiles are calculated for the selected reaction coordinates through energetic reweighting of the GaMD simulations. We have also successfully demonstrated GLOW for the characterization of activation and allosteric modulation of a G protein-coupled receptor, using the adenosine A1 receptor (A1AR) as a model system. GLOW findings are highly consistent with previous experimental and computational studies of the A1AR, while also providing further mechanistic insights into the receptor function. In summary, GLOW provides a systematic approach to mapping free energy landscapes of biomolecules. The GLOW workflow and its user manual can be downloaded at http://miaolab.org/GLOW.
Collapse
Affiliation(s)
- Hung N. Do
- The Center for Computational Biology and Department of Molecular Biosciences, The University of Kansas, Lawrence, Kansas 66047
| | - Jinan Wang
- The Center for Computational Biology and Department of Molecular Biosciences, The University of Kansas, Lawrence, Kansas 66047
| | - Apurba Bhattarai
- The Center for Computational Biology and Department of Molecular Biosciences, The University of Kansas, Lawrence, Kansas 66047
| | - Yinglong Miao
- The Center for Computational Biology and Department of Molecular Biosciences, The University of Kansas, Lawrence, Kansas 66047,Corresponding author:
| |
Collapse
|
21
|
Tze-Yang Ng J, Tan YS. Accelerated Ligand-Mapping Molecular Dynamics Simulations for the Detection of Recalcitrant Cryptic Pockets and Occluded Binding Sites. J Chem Theory Comput 2022; 18:1969-1981. [PMID: 35175753 DOI: 10.1021/acs.jctc.1c01177] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
The identification and characterization of binding sites is a critical component of structure-based drug design (SBDD). Probe-based/cosolvent molecular dynamics (MD) methods that allow for protein flexibility have been developed to predict ligand binding sites. However, cryptic pockets that appear only upon ligand binding and occluded binding sites with no access to the solvent pose significant challenges to these methods. Here, we report the development of accelerated ligand-mapping MD (aLMMD), which combines accelerated MD with LMMD, for the detection of these challenging binding sites. The method was validated on five proteins with what we term "recalcitrant" cryptic pockets, which are deeply buried pockets that require extensive movement of the protein backbone to expose, and three proteins with occluded binding sites. In all the cases, aLMMD was able to detect and sample the binding sites. Our results suggest that aLMMD could be used as a general approach for the detection of such elusive binding sites in protein targets, thus providing valuable information for SBDD.
Collapse
Affiliation(s)
- Justin Tze-Yang Ng
- Bioinformatics Institute, Agency for Science, Technology and Research (A*STAR), 30 Biopolis Street, #07-01 Matrix, Singapore 138671, Singapore
| | - Yaw Sing Tan
- Bioinformatics Institute, Agency for Science, Technology and Research (A*STAR), 30 Biopolis Street, #07-01 Matrix, Singapore 138671, Singapore
| |
Collapse
|
22
|
Tvaroška I. Glycosyltransferases as targets for therapeutic intervention in cancer and inflammation: molecular modeling insights. CHEMICAL PAPERS 2022. [DOI: 10.1007/s11696-021-02026-7] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
|
23
|
Wang J, Lan L, Wu X, Xu L, Miao Y. Mechanism of RNA recognition by a Musashi RNA-binding protein. Curr Res Struct Biol 2021; 4:10-20. [PMID: 34988468 PMCID: PMC8695263 DOI: 10.1016/j.crstbi.2021.12.002] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/08/2021] [Revised: 10/31/2021] [Accepted: 12/07/2021] [Indexed: 12/21/2022] Open
Abstract
The Musashi RNA-binding proteins (RBPs) regulate translation of target mRNAs and maintenance of cell stemness and tumorigenesis. Musashi-1 (MSI1), long considered as an intestinal and neural stem cell marker, has been more recently found to be over expressed in many cancers. It has served as an important drug target for treating acute myeloid leukemia and solid tumors such as ovarian, colorectal and bladder cancer. One of the reported binding targets of MSI1 is Numb, a negative regulator of the Notch signaling. However, the dynamic mechanism of Numb RNA binding to MSI1 remains unknown, largely hindering effective drug design targeting this critical interaction. Here, we have performed extensive all-atom microsecond-timescale simulations using a robust Gaussian accelerated molecular dynamics (GaMD) method, which successfully captured multiple times of spontaneous and highly accurate binding of the Numb RNA from bulk solvent to the MSI1 protein target site. GaMD simulations revealed that Numb RNA binding to MSI1 involved largely induced fit in both the RNA and protein. The simulations also identified important low-energy intermediate conformational states during RNA binding, in which Numb interacted mainly with the β2-β3 loop and C terminus of MSI1. The mechanistic understanding of RNA binding obtained from our GaMD simulations is expected to facilitate rational structure-based drug design targeting MSI1 and other RBPs.
Collapse
Affiliation(s)
- Jinan Wang
- Center for Computational Biology, University of Kansas, Lawrence, KS, 66047, USA
- Department of Molecular Biosciences, University of Kansas, Lawrence, KS, 66047, USA
| | - Lan Lan
- Department of Molecular Biosciences, University of Kansas, Lawrence, KS, 66047, USA
| | - Xiaoqing Wu
- Department of Molecular Biosciences, University of Kansas, Lawrence, KS, 66047, USA
| | - Liang Xu
- Department of Molecular Biosciences, University of Kansas, Lawrence, KS, 66047, USA
- Department of Radiation Oncology, The University of Kansas Cancer Center, Kansas City, KS, 66160, USA
| | - Yinglong Miao
- Center for Computational Biology, University of Kansas, Lawrence, KS, 66047, USA
- Department of Molecular Biosciences, University of Kansas, Lawrence, KS, 66047, USA
| |
Collapse
|
24
|
Wang J, Arantes PR, Bhattarai A, Hsu RV, Pawnikar S, Huang YMM, Palermo G, Miao Y. Gaussian accelerated molecular dynamics (GaMD): principles and applications. WILEY INTERDISCIPLINARY REVIEWS. COMPUTATIONAL MOLECULAR SCIENCE 2021; 11:e1521. [PMID: 34899998 PMCID: PMC8658739 DOI: 10.1002/wcms.1521] [Citation(s) in RCA: 95] [Impact Index Per Article: 31.7] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 12/09/2020] [Accepted: 01/28/2021] [Indexed: 12/20/2022]
Abstract
Gaussian accelerated molecular dynamics (GaMD) is a robust computational method for simultaneous unconstrained enhanced sampling and free energy calculations of biomolecules. It works by adding a harmonic boost potential to smooth biomolecular potential energy surface and reduce energy barriers. GaMD greatly accelerates biomolecular simulations by orders of magnitude. Without the need to set predefined reaction coordinates or collective variables, GaMD provides unconstrained enhanced sampling and is advantageous for simulating complex biological processes. The GaMD boost potential exhibits a Gaussian distribution, thereby allowing for energetic reweighting via cumulant expansion to the second order (i.e., "Gaussian approximation"). This leads to accurate reconstruction of free energy landscapes of biomolecules. Hybrid schemes with other enhanced sampling methods, such as the replica exchange GaMD (rex-GaMD) and replica exchange umbrella sampling GaMD (GaREUS), have also been introduced, further improving sampling and free energy calculations. Recently, new "selective GaMD" algorithms including the ligand GaMD (LiGaMD) and peptide GaMD (Pep-GaMD) enabled microsecond simulations to capture repetitive dissociation and binding of small-molecule ligands and highly flexible peptides. The simulations then allowed highly efficient quantitative characterization of the ligand/peptide binding thermodynamics and kinetics. Taken together, GaMD and its innovative variants are applicable to simulate a wide variety of biomolecular dynamics, including protein folding, conformational changes and allostery, ligand binding, peptide binding, protein-protein/nucleic acid/carbohydrate interactions, and carbohydrate/nucleic acid interactions. In this review, we present principles of the GaMD algorithms and recent applications in biomolecular simulations and drug design.
Collapse
Affiliation(s)
- Jinan Wang
- Center for Computational Biology and Department of Molecular Biosciences, University of Kansas, 2030 Becker Dr., Lawrence, KS, 66047, United States
| | - Pablo R Arantes
- Department of Bioengineering, University of California Riverside, 900 University Avenue, Riverside, CA 92512, United States
| | - Apurba Bhattarai
- Center for Computational Biology and Department of Molecular Biosciences, University of Kansas, 2030 Becker Dr, Lawrence, KS, 66047, United States
| | - Rohaine V Hsu
- Department of Bioengineering, University of California Riverside, 900 University Avenue, Riverside, CA 92512, United States
| | - Shristi Pawnikar
- Center for Computational Biology and Department of Molecular Biosciences, University of Kansas, 2030 Becker Dr., Lawrence, KS, 66047, United States
| | - Yu-Ming M Huang
- Department of Physics & Astronomy, Wayne State University, 666 W Hancock St, Detroit, MI 48207, USA
| | - Giulia Palermo
- Department of Bioengineering and Department of Chemistry, University of California Riverside, 900 University Avenue, Riverside, CA 92512, United States
| | - Yinglong Miao
- Center for Computational Biology and Department of Molecular Biosciences, University of Kansas, 2030 Becker Dr., Lawrence, Kansas 66047, United States
| |
Collapse
|
25
|
Chen H, Fu H, Chipot C, Shao X, Cai W. Overcoming Free-Energy Barriers with a Seamless Combination of a Biasing Force and a Collective Variable-Independent Boost Potential. J Chem Theory Comput 2021; 17:3886-3894. [PMID: 34106706 DOI: 10.1021/acs.jctc.1c00103] [Citation(s) in RCA: 14] [Impact Index Per Article: 4.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
Amid collective-variable (CV)-based importance-sampling algorithms, a hybrid of the extended adaptive biasing force and the well-tempered metadynamics algorithms (WTM-eABF) has proven particularly cost-effective for exploring the rugged free-energy landscapes that underlie biological processes. However, as an inherently CV-based algorithm, this hybrid scheme does not explicitly accelerate sampling in the space orthogonal to the chosen CVs, thereby limiting its efficiency and accuracy, most notably in those cases where the slow degrees of freedom of the process at hand are not accounted for in the model transition coordinate. Here, inspired by Gaussian-accelerated molecular dynamics (GaMD), we introduce the same CV-independent harmonic boost potential into WTM-eABF, yielding a hybrid algorithm coined GaWTM-eABF. This algorithm leans on WTM-eABF to explore the transition coordinate with a GaMD-mollified potential and recovers the unbiased free-energy landscape through thermodynamic integration followed by proper reweighting. As illustrated in our numerical tests, GaWTM-eABF effectively overcomes the free-energy barriers in orthogonal space and correctly recovers the unbiased potential of mean force (PMF). Furthermore, applying both GaWTM-eABF and WTM-eABF to two biologically relevant processes, namely, the reversible folding of (i) deca-alanine and (ii) chignolin, our results indicate that GaWTM-eABF reduces the uncertainty in the PMF calculation and converges appreciably faster than WTM-eABF. Obviating the need of multiple-copy strategies, GaWTM-eABF is a robust, computationally efficient algorithm to surmount the free-energy barriers in orthogonal space and maps with utmost fidelity the free-energy landscape along selections of CVs. Moreover, our strategy that combines WTM-eABF with GaMD can be easily extended to other biasing-force algorithms.
Collapse
Affiliation(s)
- Haochuan Chen
- Research Center for Analytical Sciences, College of Chemistry, Tianjin Key Laboratory of Biosensing and Molecular Recognition, State Key Laboratory of Medicinal Chemical Biology, Nankai University, Tianjin 300071, China
| | - Haohao Fu
- Research Center for Analytical Sciences, College of Chemistry, Tianjin Key Laboratory of Biosensing and Molecular Recognition, State Key Laboratory of Medicinal Chemical Biology, Nankai University, Tianjin 300071, China
| | - Christophe Chipot
- Laboratoire International Associé Centre National de la Recherche Scientifique et University of Illinois at Urbana-Champaign, Unité Mixte de Recherche n 7019, Université de Lorraine, BP 70239, 54506 Vandœuvre-lès-Nancy cedex, France.,Department of Physics, University of Illinois at Urbana-Champaign, Urbana, Illinois 61801, United States
| | - Xueguang Shao
- Research Center for Analytical Sciences, College of Chemistry, Tianjin Key Laboratory of Biosensing and Molecular Recognition, State Key Laboratory of Medicinal Chemical Biology, Nankai University, Tianjin 300071, China
| | - Wensheng Cai
- Research Center for Analytical Sciences, College of Chemistry, Tianjin Key Laboratory of Biosensing and Molecular Recognition, State Key Laboratory of Medicinal Chemical Biology, Nankai University, Tianjin 300071, China
| |
Collapse
|
26
|
Lambrughi M, Maiani E, Aykac Fas B, Shaw GS, Kragelund BB, Lindorff-Larsen K, Teilum K, Invernizzi G, Papaleo E. Ubiquitin Interacting Motifs: Duality Between Structured and Disordered Motifs. Front Mol Biosci 2021; 8:676235. [PMID: 34262938 PMCID: PMC8273247 DOI: 10.3389/fmolb.2021.676235] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/04/2021] [Accepted: 05/14/2021] [Indexed: 01/11/2023] Open
Abstract
Ubiquitin is a small protein at the heart of many cellular processes, and several different protein domains are known to recognize and bind ubiquitin. A common motif for interaction with ubiquitin is the Ubiquitin Interacting Motif (UIM), characterized by a conserved sequence signature and often found in multi-domain proteins. Multi-domain proteins with intrinsically disordered regions mediate interactions with multiple partners, orchestrating diverse pathways. Short linear motifs for binding are often embedded in these disordered regions and play crucial roles in modulating protein function. In this work, we investigated the structural propensities of UIMs using molecular dynamics simulations and NMR chemical shifts. Despite the structural portrait depicted by X-crystallography of stable helical structures, we show that UIMs feature both helical and intrinsically disordered conformations. Our results shed light on a new class of disordered UIMs. This group is here exemplified by the C-terminal domain of one isoform of ataxin-3 and a group of ubiquitin-specific proteases. Intriguingly, UIMs not only bind ubiquitin. They can be a recruitment point for other interactors, such as parkin and the heat shock protein Hsc70-4. Disordered UIMs can provide versatility and new functions to the client proteins, opening new directions for research on their interactome.
Collapse
Affiliation(s)
- Matteo Lambrughi
- Computational Biology Laboratory, Danish Cancer Society Research Center, Copenhagen, Denmark.,Department of Biotechnology and Bioscience, University of Milano-Bicocca, Milano, Italy
| | - Emiliano Maiani
- Computational Biology Laboratory, Danish Cancer Society Research Center, Copenhagen, Denmark
| | - Burcu Aykac Fas
- Computational Biology Laboratory, Danish Cancer Society Research Center, Copenhagen, Denmark
| | - Gary S Shaw
- Department of Biochemistry, Schulich School of Medicine and Dentistry, The University of Western Ontario, London, ON, Canada
| | - Birthe B Kragelund
- Structural Biology and NMR Laboratory and The Linderstrøm-Lang Centre for Protein Science, Department of Biology, University of Copenhagen, Copenhagen, Denmark
| | - Kresten Lindorff-Larsen
- Structural Biology and NMR Laboratory and The Linderstrøm-Lang Centre for Protein Science, Department of Biology, University of Copenhagen, Copenhagen, Denmark
| | - Kaare Teilum
- Structural Biology and NMR Laboratory and The Linderstrøm-Lang Centre for Protein Science, Department of Biology, University of Copenhagen, Copenhagen, Denmark
| | - Gaetano Invernizzi
- Structural Biology and NMR Laboratory and The Linderstrøm-Lang Centre for Protein Science, Department of Biology, University of Copenhagen, Copenhagen, Denmark
| | - Elena Papaleo
- Computational Biology Laboratory, Danish Cancer Society Research Center, Copenhagen, Denmark.,Cancer Systems Biology, Section for Bioinformatics, Department of Health and Technology, Technical University of Denmark, Lyngby, Denmark
| |
Collapse
|
27
|
Kingdon ADH, Alderwick LJ. Structure-based in silico approaches for drug discovery against Mycobacterium tuberculosis. Comput Struct Biotechnol J 2021; 19:3708-3719. [PMID: 34285773 PMCID: PMC8258792 DOI: 10.1016/j.csbj.2021.06.034] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/23/2021] [Revised: 06/22/2021] [Accepted: 06/22/2021] [Indexed: 12/12/2022] Open
Abstract
Mycobacterium tuberculosis is the causative agent of TB and was estimated to cause 1.4 million death in 2019, alongside 10 million new infections. Drug resistance is a growing issue, with multi-drug resistant infections representing 3.3% of all new infections, hence novel antimycobacterial drugs are urgently required to combat this growing health emergency. Alongside this, increased knowledge of gene essentiality in the pathogenic organism and larger compound databases can aid in the discovery of new drug compounds. The number of protein structures, X-ray based and modelled, is increasing and now accounts for greater than > 80% of all predicted M. tuberculosis proteins; allowing novel targets to be investigated. This review will focus on structure-based in silico approaches for drug discovery, covering a range of complexities and computational demands, with associated antimycobacterial examples. This includes molecular docking, molecular dynamic simulations, ensemble docking and free energy calculations. Applications of machine learning onto each of these approaches will be discussed. The need for experimental validation of computational hits is an essential component, which is unfortunately missing from many current studies. The future outlooks of these approaches will also be discussed.
Collapse
Key Words
- CV, collective variable
- Docking
- Drug discovery
- In silico
- LIE, Linear Interaction Energy
- MD, Molecular Dynamic
- MDR, multi-drug resistant
- MMPB(GB)SA, Molecular Mechanics with Poisson Boltzmann (or generalised Born) and Surface Area solvation
- Machine learning
- Mt, Mycobacterium tuberculosis
- Mycobacterium tuberculosis
- PTC, peptidyl transferase centre
- RMSD, root-mean square-deviation
- Tuberculosis, TB
- cMD, Classical Molecular Dynamic
- cryo-EM, cryogenic electron microscopy
- ns, nanosecond
Collapse
Affiliation(s)
- Alexander D H Kingdon
- Institute of Microbiology and Infection, School of Biosciences, University of Birmingham, Edgbaston, Birmingham B15 2TT, United Kingdom
| | - Luke J Alderwick
- Institute of Microbiology and Infection, School of Biosciences, University of Birmingham, Edgbaston, Birmingham B15 2TT, United Kingdom
| |
Collapse
|
28
|
Paissoni C, Camilloni C. How to Determine Accurate Conformational Ensembles by Metadynamics Metainference: A Chignolin Study Case. Front Mol Biosci 2021; 8:694130. [PMID: 34124166 PMCID: PMC8187852 DOI: 10.3389/fmolb.2021.694130] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/12/2021] [Accepted: 05/14/2021] [Indexed: 11/13/2022] Open
Abstract
The reliability and usefulness of molecular dynamics simulations of equilibrium processes rests on their statistical precision and their capability to generate conformational ensembles in agreement with available experimental knowledge. Metadynamics Metainference (M&M), coupling molecular dynamics with the enhanced sampling ability of Metadynamics and with the ability to integrate experimental information of Metainference, can in principle achieve both goals. Here we show that three different Metadynamics setups provide converged estimate of the populations of the three-states populated by a model peptide. Errors are estimated correctly by block averaging, but higher precision is obtained by performing independent replicates. One effect of Metadynamics is that of dramatically decreasing the number of effective frames resulting from the simulations and this is relevant for M&M where the number of replicas should be large enough to capture the conformational heterogeneity behind the experimental data. Our simulations allow also us to propose that monitoring the relative error associated with conformational averaging can help to determine the minimum number of replicas to be simulated in the context of M&M simulations. Altogether our data provides useful indication on how to generate sound conformational ensemble in agreement with experimental data.
Collapse
Affiliation(s)
- Cristina Paissoni
- Dipartimento di Bioscienze, Università degli Studi di Milano, Milan, Italy
| | - Carlo Camilloni
- Dipartimento di Bioscienze, Università degli Studi di Milano, Milan, Italy
| |
Collapse
|
29
|
Wang J, Miao Y. Peptide Gaussian accelerated molecular dynamics (Pep-GaMD): Enhanced sampling and free energy and kinetics calculations of peptide binding. J Chem Phys 2021; 153:154109. [PMID: 33092378 DOI: 10.1063/5.0021399] [Citation(s) in RCA: 66] [Impact Index Per Article: 22.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/17/2022] Open
Abstract
Peptides mediate up to 40% of known protein-protein interactions in higher eukaryotes and play an important role in cellular signaling. However, it is challenging to simulate both binding and unbinding of peptides and calculate peptide binding free energies through conventional molecular dynamics, due to long biological timescales and extremely high flexibility of the peptides. Based on the Gaussian accelerated molecular dynamics (GaMD) enhanced sampling technique, we have developed a new computational method "Pep-GaMD," which selectively boosts essential potential energy of the peptide in order to effectively model its high flexibility. In addition, another boost potential is applied to the remaining potential energy of the entire system in a dual-boost algorithm. Pep-GaMD has been demonstrated on binding of three model peptides to the SH3 domains. Independent 1 µs dual-boost Pep-GaMD simulations have captured repetitive peptide dissociation and binding events, which enable us to calculate peptide binding thermodynamics and kinetics. The calculated binding free energies and kinetic rate constants agreed very well with available experimental data. Furthermore, the all-atom Pep-GaMD simulations have provided important insights into the mechanism of peptide binding to proteins that involves long-range electrostatic interactions and mainly conformational selection. In summary, Pep-GaMD provides a highly efficient, easy-to-use approach for unconstrained enhanced sampling and calculations of peptide binding free energies and kinetics.
Collapse
Affiliation(s)
- Jinan Wang
- Center for Computational Biology and Department of Molecular Biosciences, University of Kansas, Lawrence, Kansas 66047, USA
| | - Yinglong Miao
- Center for Computational Biology and Department of Molecular Biosciences, University of Kansas, Lawrence, Kansas 66047, USA
| |
Collapse
|
30
|
Abstract
Human cytochrome P450 1B1 (CYP1B1) is an extrahepatic heme-containing monooxygenase. CYP1B1 contributes to the oxidative metabolism of xenobiotics, drugs, and endogenous substrates like melatonin, fatty acids, steroid hormones, and retinoids, which are involved in diverse critical cellular functions. CYP1B1 plays an important role in the pathogenesis of cardiovascular diseases, hormone-related cancers and is responsible for anti-cancer drug resistance. Inhibition of CYP1B1 activity is considered as an approach in cancer chemoprevention and cancer chemotherapy. CYP1B1 can activate anti-cancer prodrugs in tumor cells which display overexpression of CYP1B1 in comparison to normal cells. CYP1B1 involvement in carcinogenesis and cancer progression encourages investigation of CYP1B1 interactions with its ligands: substrates and inhibitors. Computational methods, with a simulation of molecular dynamics (MD), allow the observation of molecular interactions at the binding site of CYP1B1, which are essential in relation to the enzyme’s functions.
Collapse
|
31
|
Planas-Iglesias J, Marques SM, Pinto GP, Musil M, Stourac J, Damborsky J, Bednar D. Computational design of enzymes for biotechnological applications. Biotechnol Adv 2021; 47:107696. [PMID: 33513434 DOI: 10.1016/j.biotechadv.2021.107696] [Citation(s) in RCA: 36] [Impact Index Per Article: 12.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/07/2020] [Revised: 01/12/2021] [Accepted: 01/13/2021] [Indexed: 12/14/2022]
Abstract
Enzymes are the natural catalysts that execute biochemical reactions upholding life. Their natural effectiveness has been fine-tuned as a result of millions of years of natural evolution. Such catalytic effectiveness has prompted the use of biocatalysts from multiple sources on different applications, including the industrial production of goods (food and beverages, detergents, textile, and pharmaceutics), environmental protection, and biomedical applications. Natural enzymes often need to be improved by protein engineering to optimize their function in non-native environments. Recent technological advances have greatly facilitated this process by providing the experimental approaches of directed evolution or by enabling computer-assisted applications. Directed evolution mimics the natural selection process in a highly accelerated fashion at the expense of arduous laboratory work and economic resources. Theoretical methods provide predictions and represent an attractive complement to such experiments by waiving their inherent costs. Computational techniques can be used to engineer enzymatic reactivity, substrate specificity and ligand binding, access pathways and ligand transport, and global properties like protein stability, solubility, and flexibility. Theoretical approaches can also identify hotspots on the protein sequence for mutagenesis and predict suitable alternatives for selected positions with expected outcomes. This review covers the latest advances in computational methods for enzyme engineering and presents many successful case studies.
Collapse
Affiliation(s)
- Joan Planas-Iglesias
- Loschmidt Laboratories, Department of Experimental Biology and RECETOX, Faculty of Science, Masaryk University, Kamenice 5/A13, 625 00 Brno, Czech Republic; International Clinical Research Center, St. Anne's University Hospital Brno, Pekarska 53, 656 91 Brno, Czech Republic
| | - Sérgio M Marques
- Loschmidt Laboratories, Department of Experimental Biology and RECETOX, Faculty of Science, Masaryk University, Kamenice 5/A13, 625 00 Brno, Czech Republic; International Clinical Research Center, St. Anne's University Hospital Brno, Pekarska 53, 656 91 Brno, Czech Republic
| | - Gaspar P Pinto
- Loschmidt Laboratories, Department of Experimental Biology and RECETOX, Faculty of Science, Masaryk University, Kamenice 5/A13, 625 00 Brno, Czech Republic; International Clinical Research Center, St. Anne's University Hospital Brno, Pekarska 53, 656 91 Brno, Czech Republic
| | - Milos Musil
- Loschmidt Laboratories, Department of Experimental Biology and RECETOX, Faculty of Science, Masaryk University, Kamenice 5/A13, 625 00 Brno, Czech Republic; International Clinical Research Center, St. Anne's University Hospital Brno, Pekarska 53, 656 91 Brno, Czech Republic; IT4Innovations Centre of Excellence, Faculty of Information Technology, Brno University of Technology, 61266 Brno, Czech Republic
| | - Jan Stourac
- Loschmidt Laboratories, Department of Experimental Biology and RECETOX, Faculty of Science, Masaryk University, Kamenice 5/A13, 625 00 Brno, Czech Republic; International Clinical Research Center, St. Anne's University Hospital Brno, Pekarska 53, 656 91 Brno, Czech Republic
| | - Jiri Damborsky
- Loschmidt Laboratories, Department of Experimental Biology and RECETOX, Faculty of Science, Masaryk University, Kamenice 5/A13, 625 00 Brno, Czech Republic; International Clinical Research Center, St. Anne's University Hospital Brno, Pekarska 53, 656 91 Brno, Czech Republic.
| | - David Bednar
- Loschmidt Laboratories, Department of Experimental Biology and RECETOX, Faculty of Science, Masaryk University, Kamenice 5/A13, 625 00 Brno, Czech Republic.
| |
Collapse
|
32
|
Fadda E. Understanding the Structure and Function of Viral Glycosylation by Molecular Simulations: State-of-the-Art and Recent Case Studies. COMPREHENSIVE GLYCOSCIENCE 2021. [PMCID: PMC7834635 DOI: 10.1016/b978-0-12-819475-1.00056-0] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Indexed: 02/06/2023]
|
33
|
Gonzalez TR, Martin KP, Barnes JE, Patel JS, Ytreberg FM. Assessment of software methods for estimating protein-protein relative binding affinities. PLoS One 2020; 15:e0240573. [PMID: 33347442 PMCID: PMC7751979 DOI: 10.1371/journal.pone.0240573] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/26/2020] [Accepted: 12/07/2020] [Indexed: 11/19/2022] Open
Abstract
A growing number of computational tools have been developed to accurately and rapidly predict the impact of amino acid mutations on protein-protein relative binding affinities. Such tools have many applications, for example, designing new drugs and studying evolutionary mechanisms. In the search for accuracy, many of these methods employ expensive yet rigorous molecular dynamics simulations. By contrast, non-rigorous methods use less exhaustive statistical mechanics, allowing for more efficient calculations. However, it is unclear if such methods retain enough accuracy to replace rigorous methods in binding affinity calculations. This trade-off between accuracy and computational expense makes it difficult to determine the best method for a particular system or study. Here, eight non-rigorous computational methods were assessed using eight antibody-antigen and eight non-antibody-antigen complexes for their ability to accurately predict relative binding affinities (ΔΔG) for 654 single mutations. In addition to assessing accuracy, we analyzed the CPU cost and performance for each method using a variety of physico-chemical structural features. This allowed us to posit scenarios in which each method may be best utilized. Most methods performed worse when applied to antibody-antigen complexes compared to non-antibody-antigen complexes. Rosetta-based JayZ and EasyE methods classified mutations as destabilizing (ΔΔG < -0.5 kcal/mol) with high (83-98%) accuracy and a relatively low computational cost for non-antibody-antigen complexes. Some of the most accurate results for antibody-antigen systems came from combining molecular dynamics with FoldX with a correlation coefficient (r) of 0.46, but this was also the most computationally expensive method. Overall, our results suggest these methods can be used to quickly and accurately predict stabilizing versus destabilizing mutations but are less accurate at predicting actual binding affinities. This study highlights the need for continued development of reliable, accessible, and reproducible methods for predicting binding affinities in antibody-antigen proteins and provides a recipe for using current methods.
Collapse
Affiliation(s)
- Tawny R. Gonzalez
- Institute for Modeling Collaboration and Innovation, University of Idaho, Moscow, Idaho, United States of America
| | - Kyle P. Martin
- Institute for Modeling Collaboration and Innovation, University of Idaho, Moscow, Idaho, United States of America
- Department of Physics, University of Idaho, Moscow, Idaho, United States of America
| | - Jonathan E. Barnes
- Institute for Modeling Collaboration and Innovation, University of Idaho, Moscow, Idaho, United States of America
- Department of Physics, University of Idaho, Moscow, Idaho, United States of America
| | - Jagdish Suresh Patel
- Institute for Modeling Collaboration and Innovation, University of Idaho, Moscow, Idaho, United States of America
- Department of Biological Sciences, University of Idaho, Moscow, Idaho, United States of America
| | - F. Marty Ytreberg
- Institute for Modeling Collaboration and Innovation, University of Idaho, Moscow, Idaho, United States of America
- Department of Physics, University of Idaho, Moscow, Idaho, United States of America
| |
Collapse
|
34
|
Medeiros Selegato D, Bracco C, Giannelli C, Parigi G, Luchinat C, Sgheri L, Ravera E. Comparison of Different Reweighting Approaches for the Calculation of Conformational Variability of Macromolecules from Molecular Simulations. Chemphyschem 2020; 22:127-138. [DOI: 10.1002/cphc.202000714] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/17/2020] [Revised: 09/14/2020] [Indexed: 11/07/2022]
Affiliation(s)
- Denise Medeiros Selegato
- Magnetic Resonance Center (CERM) and Interuniversity Consortium for Magnetic Resonance of Metallo Proteins (CIRMMP) Via L. Sacconi 6 50019 Sesto Fiorentino Italy
- Dipartimento di Chimica “Ugo Schiff” Università degli Studi di Firenze Via della Lastruccia 3 50019 Sesto Fiorentino Italy
- Present address: Fundación MEDINA, Centro de Excelentia en Investigación de Medicamentos Innovadores and Andalucía MSD España Granada Spain
| | - Cesare Bracco
- Dipartimento di Matematica e Informatica “U. Dini” Università degli Studi di Firenze Viale Morgagni 67/a 50134 Florence Italy
| | - Carlotta Giannelli
- Dipartimento di Matematica e Informatica “U. Dini” Università degli Studi di Firenze Viale Morgagni 67/a 50134 Florence Italy
| | - Giacomo Parigi
- Magnetic Resonance Center (CERM) and Interuniversity Consortium for Magnetic Resonance of Metallo Proteins (CIRMMP) Via L. Sacconi 6 50019 Sesto Fiorentino Italy
- Dipartimento di Chimica “Ugo Schiff” Università degli Studi di Firenze Via della Lastruccia 3 50019 Sesto Fiorentino Italy
| | - Claudio Luchinat
- Magnetic Resonance Center (CERM) and Interuniversity Consortium for Magnetic Resonance of Metallo Proteins (CIRMMP) Via L. Sacconi 6 50019 Sesto Fiorentino Italy
- Dipartimento di Chimica “Ugo Schiff” Università degli Studi di Firenze Via della Lastruccia 3 50019 Sesto Fiorentino Italy
| | - Luca Sgheri
- Istituto per le Applicazioni del Calcolo (CNR) sede di Firenze via Madonna del Piano 10 50019 Sesto Fiorentino Italy
| | - Enrico Ravera
- Magnetic Resonance Center (CERM) and Interuniversity Consortium for Magnetic Resonance of Metallo Proteins (CIRMMP) Via L. Sacconi 6 50019 Sesto Fiorentino Italy
- Dipartimento di Chimica “Ugo Schiff” Università degli Studi di Firenze Via della Lastruccia 3 50019 Sesto Fiorentino Italy
| |
Collapse
|
35
|
Monroe JI, Hatch HW, Mahynski NA, Shell MS, Shen VK. Extrapolation and interpolation strategies for efficiently estimating structural observables as a function of temperature and density. J Chem Phys 2020; 153:144101. [PMID: 33086808 DOI: 10.1063/5.0014282] [Citation(s) in RCA: 9] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/28/2022] Open
Abstract
Thermodynamic extrapolation has previously been used to predict arbitrary structural observables in molecular simulations at temperatures (or relative chemical potentials in open-system mixtures) different from those at which the simulation was performed. This greatly reduces the computational cost in mapping out phase and structural transitions. In this work, we explore the limitations and accuracy of thermodynamic extrapolation applied to water, where qualitative shifts from anomalous to simple-fluid-like behavior are manifested through shifts in the liquid structure that occur as a function of both temperature and density. We present formulas for extrapolating in volume for canonical ensembles and demonstrate that linear extrapolations of water's structural properties are only accurate over a limited density range. On the other hand, linear extrapolation in temperature can be accurate across the entire liquid state. We contrast these extrapolations with classical perturbation theory techniques, which are more conservative and slowly converging. Indeed, we show that such behavior is expected by demonstrating exact relationships between extrapolation of free energies and well-known techniques to predict free energy differences. An ideal gas in an external field is also studied to more clearly explain these results for a toy system with fully analytical solutions. We also present a recursive interpolation strategy for predicting arbitrary structural properties of molecular fluids over a predefined range of state conditions, demonstrating its success in mapping qualitative shifts in water structure with density.
Collapse
Affiliation(s)
- Jacob I Monroe
- National Institute of Standards and Technology, Gaithersburg, Maryland 20899, USA
| | - Harold W Hatch
- National Institute of Standards and Technology, Gaithersburg, Maryland 20899, USA
| | - Nathan A Mahynski
- National Institute of Standards and Technology, Gaithersburg, Maryland 20899, USA
| | - M Scott Shell
- University of California - Santa Barbara, Santa Barbara, California 93106, USA
| | - Vincent K Shen
- National Institute of Standards and Technology, Gaithersburg, Maryland 20899, USA
| |
Collapse
|
36
|
Sora V, Kumar M, Maiani E, Lambrughi M, Tiberti M, Papaleo E. Structure and Dynamics in the ATG8 Family From Experimental to Computational Techniques. Front Cell Dev Biol 2020; 8:420. [PMID: 32587856 PMCID: PMC7297954 DOI: 10.3389/fcell.2020.00420] [Citation(s) in RCA: 14] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/15/2020] [Accepted: 05/06/2020] [Indexed: 12/31/2022] Open
Abstract
Autophagy is a conserved and essential intracellular mechanism for the removal of damaged components. Since autophagy deregulation is linked to different kinds of pathologies, it is fundamental to gain knowledge on the fine molecular and structural details related to the core proteins of the autophagy machinery. Among these, the family of human ATG8 proteins plays a central role in recruiting other proteins to the different membrane structures involved in the autophagic pathway. Several experimental structures are available for the members of the ATG8 family alone or in complex with their different biological partners, including disordered regions of proteins containing a short linear motif called LC3 interacting motif. Recently, the first structural details of the interaction of ATG8 proteins with biological membranes came into light. The availability of structural data for human ATG8 proteins has been paving the way for studies on their structure-function-dynamic relationship using biomolecular simulations. Experimental and computational structural biology can help to address several outstanding questions on the mechanism of human ATG8 proteins, including their specificity toward different interactors, their association with membranes, the heterogeneity of their conformational ensemble, and their regulation by post-translational modifications. We here summarize the main results collected so far and discuss the future perspectives within the field and the knowledge gaps. Our review can serve as a roadmap for future structural and dynamics studies of the ATG8 family members in health and disease.
Collapse
Affiliation(s)
- Valentina Sora
- Computational Biology Laboratory, Center for Autophagy, Recycling and Disease, Danish Cancer Society Research Center, Copenhagen, Denmark
| | - Mukesh Kumar
- Computational Biology Laboratory, Center for Autophagy, Recycling and Disease, Danish Cancer Society Research Center, Copenhagen, Denmark
| | - Emiliano Maiani
- Computational Biology Laboratory, Center for Autophagy, Recycling and Disease, Danish Cancer Society Research Center, Copenhagen, Denmark
| | - Matteo Lambrughi
- Computational Biology Laboratory, Center for Autophagy, Recycling and Disease, Danish Cancer Society Research Center, Copenhagen, Denmark
| | - Matteo Tiberti
- Computational Biology Laboratory, Center for Autophagy, Recycling and Disease, Danish Cancer Society Research Center, Copenhagen, Denmark
| | - Elena Papaleo
- Computational Biology Laboratory, Center for Autophagy, Recycling and Disease, Danish Cancer Society Research Center, Copenhagen, Denmark
- Translational Disease System Biology, Faculty of Health and Medical Sciences, Novo Nordisk Foundation Center for Protein Research, University of Copenhagen, Copenhagen, Denmark
| |
Collapse
|
37
|
Nishimura Y, Nakai H. Hierarchical parallelization of divide‐and‐conquer density functional tight‐binding molecular dynamics and metadynamics simulations. J Comput Chem 2020; 41:1759-1772. [DOI: 10.1002/jcc.26217] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/19/2020] [Revised: 04/15/2020] [Accepted: 04/20/2020] [Indexed: 11/08/2022]
Affiliation(s)
- Yoshifumi Nishimura
- Waseda Research Institute for Science and Engineering Waseda University Tokyo Japan
| | - Hiromi Nakai
- Waseda Research Institute for Science and Engineering Waseda University Tokyo Japan
- Department of Chemistry and Biochemistry School of Advanced Science and Engineering, Waseda University Tokyo Japan
- Elements Strategy Initiative for Catalysts and Batteries Kyoto University Kyoto Japan
| |
Collapse
|
38
|
Horton JT, Allen AEA, Cole DJ. Modelling flexible protein–ligand binding in p38α MAP kinase using the QUBE force field. Chem Commun (Camb) 2020; 56:932-935. [DOI: 10.1039/c9cc08574b] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/08/2023]
Abstract
The accuracy of quantum mechanical bespoke (QUBE) force fields for protein–ligand binding free energy calculations are benchmarked against experiment.
Collapse
Affiliation(s)
- Joshua T. Horton
- School of Natural and Environmental Sciences
- Newcastle University
- Newcastle upon Tyne NE1 7RU
- UK
| | | | - Daniel J. Cole
- School of Natural and Environmental Sciences
- Newcastle University
- Newcastle upon Tyne NE1 7RU
- UK
| |
Collapse
|
39
|
Wang J, Alekseenko A, Kozakov D, Miao Y. Improved Modeling of Peptide-Protein Binding Through Global Docking and Accelerated Molecular Dynamics Simulations. Front Mol Biosci 2019; 6:112. [PMID: 31737642 PMCID: PMC6835073 DOI: 10.3389/fmolb.2019.00112] [Citation(s) in RCA: 43] [Impact Index Per Article: 8.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/22/2019] [Accepted: 10/09/2019] [Indexed: 01/31/2023] Open
Abstract
Peptides mediate up to 40% of known protein-protein interactions in higher eukaryotes and play a key role in cellular signaling, protein trafficking, immunology, and oncology. However, it is challenging to predict peptide-protein binding with conventional computational modeling approaches, due to slow dynamics and high peptide flexibility. Here, we present a prototype of the approach which combines global peptide docking using ClusPro PeptiDock and all-atom enhanced simulations using Gaussian accelerated molecular dynamics (GaMD). For three distinct model peptides, the lowest backbone root-mean-square deviations (RMSDs) of their bound conformations relative to X-ray structures obtained from PeptiDock were 3.3–4.8 Å, being medium quality predictions according to the Critical Assessment of PRediction of Interactions (CAPRI) criteria. GaMD simulations refined the peptide-protein complex structures with significantly reduced peptide backbone RMSDs of 0.6–2.7 Å, yielding two high quality (sub-angstrom) and one medium quality models. Furthermore, the GaMD simulations identified important low-energy conformational states and revealed the mechanism of peptide binding to the target proteins. Therefore, PeptiDock+GaMD is a promising approach for exploring peptide-protein interactions.
Collapse
Affiliation(s)
- Jinan Wang
- Center for Computational Biology and Department of Molecular Biosciences, University of Kansas, Lawrence, KS, United States
| | - Andrey Alekseenko
- Laufer Center for Physical and Quantitative Biology, Stony Brook University, Stony Brook, NY, United States.,Department of Applied Mathematics and Statistics, Stony Brook University, Stony Brook, NY, United States
| | - Dima Kozakov
- Laufer Center for Physical and Quantitative Biology, Stony Brook University, Stony Brook, NY, United States.,Department of Applied Mathematics and Statistics, Stony Brook University, Stony Brook, NY, United States
| | - Yinglong Miao
- Center for Computational Biology and Department of Molecular Biosciences, University of Kansas, Lawrence, KS, United States
| |
Collapse
|
40
|
Li Y, Wang M, Gao N, Li D, Lin J. The effect of dimerization on the activation and conformational dynamics of adenosine A 1 receptor. Phys Chem Chem Phys 2019; 21:22763-22773. [PMID: 31595279 DOI: 10.1039/c9cp04060a] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
The adenosine A1 receptor (A1R) is one of four adenosine receptors in humans, which are involved in the function of the cardiovascular, respiratory and central nervous systems. Experimental results indicate that A1R can form a homodimer and that the protomer-protomer interaction in the A1R dimer is related to certain pharmacological characteristics of A1R activation. In this work, we performed docking, metadynamics simulation, conventional molecular dynamics simulations, Gaussian-accelerated molecular dynamics simulations, potential of mean force calculations, dynamic cross-correlation motions analysis and community network analysis to study the binding mode of 5'-N-ethylcarboxamidoadenosine (NECA) to A1R and the effect of dimerization on the activation of A1R. Our results show that NECA binds to A1R in a similar mode to adenosine in the A1R crystal structure and NECA in the A2AR crystal structure. The A1R homodimer can be activated by one or two agonists with NECA occupying its orthosteric pockets in one (which we call the NECA-A1R system) or both protomers (which we call the dNECA-A1R system). In the NECA-A1R system, activation is predicated in the protomer without NECA bound. In the dNECA-A1R system, only one protomer achieves the active state. These findings suggest an asymmetrical activation mechanism of the homodimer and a negative cooperativity between the two protomers. We envision that our results may further facilitate the drug development of A1R.
Collapse
Affiliation(s)
- Yang Li
- State Key Laboratory of Medicinal Chemical Biology, College of Pharmacy and Tianjin Key Laboratory of Molecular Drug Research, Nankai University, Haihe Education Park, 38 Tongyan Road, Tianjin 300353, People's Republic of China.
| | | | | | | | | |
Collapse
|
41
|
Wang AH, Zhang ZC, Li GH. Advances in enhanced sampling molecular dynamics simulations for biomolecules. CHINESE J CHEM PHYS 2019. [DOI: 10.1063/1674-0068/cjcp1905091] [Citation(s) in RCA: 20] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/12/2022]
Affiliation(s)
- An-hui Wang
- Laboratory of Molecular Modeling and Design, State Key Laboratory of Molecular Reaction Dynamics, Dalian Institute of Chemical Physics, Chinese Academy of Sciences, Dalian 116023, China
- State Key Laboratory of Fine Chemicals, School of Chemistry, Dalian University of Technology, Dalian 116024, China
| | - Zhi-chao Zhang
- State Key Laboratory of Fine Chemicals, School of Chemistry, Dalian University of Technology, Dalian 116024, China
| | - Guo-hui Li
- Laboratory of Molecular Modeling and Design, State Key Laboratory of Molecular Reaction Dynamics, Dalian Institute of Chemical Physics, Chinese Academy of Sciences, Dalian 116023, China
| |
Collapse
|
42
|
Trapl D, Horvacanin I, Mareska V, Ozcelik F, Unal G, Spiwok V. Anncolvar: Approximation of Complex Collective Variables by Artificial Neural Networks for Analysis and Biasing of Molecular Simulations. Front Mol Biosci 2019; 6:25. [PMID: 31058167 PMCID: PMC6482212 DOI: 10.3389/fmolb.2019.00025] [Citation(s) in RCA: 18] [Impact Index Per Article: 3.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/22/2019] [Accepted: 04/01/2019] [Indexed: 11/23/2022] Open
Abstract
The state of a molecular system can be described in terms of collective variables. These low-dimensional descriptors of molecular structure can be used to monitor the state of the simulation, to calculate free energy profiles or to accelerate rare events by a bias potential or a bias force. Frequent calculation of some complex collective variables may slow down the simulation or analysis of trajectories. Moreover, many collective variables cannot be explicitly calculated for newly sampled structures. In order to address this problem, we developed a new package called anncolvar. This package makes it possible to build and train an artificial neural network model that approximates a collective variable. It can be used to generate an input for the open-source enhanced sampling simulation PLUMED package, so the collective variable can be monitored and biased by methods available in this program. The computational efficiency and the accuracy of anncolvar are demonstrated on selected molecular systems (cyclooctane derivative, Trp-cage miniprotein) and selected collective variables (Isomap, molecular surface area).
Collapse
Affiliation(s)
- Dalibor Trapl
- Department of Biochemistry and Microbiology, University of Chemistry and Technology in Prague, Prague, Czechia
| | - Izabela Horvacanin
- Department of Biochemistry and Microbiology, University of Chemistry and Technology in Prague, Prague, Czechia.,Faculty of Science, University of Zagreb, Zagreb, Croatia
| | - Vaclav Mareska
- Department of Biochemistry and Microbiology, University of Chemistry and Technology in Prague, Prague, Czechia
| | - Furkan Ozcelik
- Computer Engineering Department, Istanbul Technical University, Istanbul, Turkey
| | - Gozde Unal
- Computer Engineering Department, Istanbul Technical University, Istanbul, Turkey
| | - Vojtech Spiwok
- Department of Biochemistry and Microbiology, University of Chemistry and Technology in Prague, Prague, Czechia
| |
Collapse
|
43
|
Bignon E, Rizza S, Filomeni G, Papaleo E. Use of Computational Biochemistry for Elucidating Molecular Mechanisms of Nitric Oxide Synthase. Comput Struct Biotechnol J 2019; 17:415-429. [PMID: 30996821 PMCID: PMC6451115 DOI: 10.1016/j.csbj.2019.03.011] [Citation(s) in RCA: 20] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/21/2018] [Revised: 03/17/2019] [Accepted: 03/21/2019] [Indexed: 12/25/2022] Open
Abstract
Nitric oxide (NO) is an essential signaling molecule in the regulation of multiple cellular processes. It is endogenously synthesized by NO synthase (NOS) as the product of L-arginine oxidation to L-citrulline, requiring NADPH, molecular oxygen, and a pterin cofactor. Two NOS isoforms are constitutively present in cells, nNOS and eNOS, and a third is inducible (iNOS). Despite their biological relevance, the details of their complex structural features and reactivity mechanisms are still unclear. In this review, we summarized the contribution of computational biochemistry to research on NOS molecular mechanisms. We described in detail its use in studying aspects of structure, dynamics and reactivity. We also focus on the numerous outstanding questions in the field that could benefit from more extensive computational investigations.
Collapse
Affiliation(s)
- Emmanuelle Bignon
- Computational Biology Laboratory, Danish Cancer Society Research Center, Strandboulevarden 49, 2100 Copenhagen, Denmark
| | - Salvatore Rizza
- Redox Signaling and Oxidative Stress Group, Cell Stress and Survival Unit, Danish Cancer Society Research Center, Strandboulevarden 49, 2100 Copenhagen, Denmark
| | - Giuseppe Filomeni
- Redox Signaling and Oxidative Stress Group, Cell Stress and Survival Unit, Danish Cancer Society Research Center, Strandboulevarden 49, 2100 Copenhagen, Denmark.,Department of Biology, University of Rome Tor Vergata, Rome, Italy
| | - Elena Papaleo
- Computational Biology Laboratory, Danish Cancer Society Research Center, Strandboulevarden 49, 2100 Copenhagen, Denmark.,Translational Disease Systems Biology, Faculty of Health and Medical Sciences, Novo Nordisk Foundation Center for Protein Research University of Copenhagen, Copenhagen, Denmark
| |
Collapse
|
44
|
Enkavi G, Javanainen M, Kulig W, Róg T, Vattulainen I. Multiscale Simulations of Biological Membranes: The Challenge To Understand Biological Phenomena in a Living Substance. Chem Rev 2019; 119:5607-5774. [PMID: 30859819 PMCID: PMC6727218 DOI: 10.1021/acs.chemrev.8b00538] [Citation(s) in RCA: 173] [Impact Index Per Article: 34.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/23/2022]
Abstract
![]()
Biological
membranes are tricky to investigate. They are complex
in terms of molecular composition and structure, functional
over a wide range of time scales, and characterized
by nonequilibrium conditions. Because of all of these
features, simulations are a great technique to study biomembrane
behavior. A significant part of the functional processes
in biological membranes takes place at the molecular
level; thus computer simulations are the method of
choice to explore how their properties emerge from specific
molecular features and how the interplay among the numerous
molecules gives rise to function over spatial and
time scales larger than the molecular ones. In this
review, we focus on this broad theme. We discuss the current
state-of-the-art of biomembrane simulations that, until
now, have largely focused on a rather narrow picture
of the complexity of the membranes. Given this, we
also discuss the challenges that we should unravel in the
foreseeable future. Numerous features such as the actin-cytoskeleton
network, the glycocalyx network, and nonequilibrium
transport under ATP-driven conditions have so far
received very little attention; however, the potential
of simulations to solve them would be exceptionally high. A
major milestone for this research would be that one day
we could say that computer simulations genuinely research
biological membranes, not just lipid bilayers.
Collapse
Affiliation(s)
- Giray Enkavi
- Department of Physics , University of Helsinki , P.O. Box 64, FI-00014 Helsinki , Finland
| | - Matti Javanainen
- Department of Physics , University of Helsinki , P.O. Box 64, FI-00014 Helsinki , Finland.,Institute of Organic Chemistry and Biochemistry of the Czech Academy of Sciences , Flemingovo naḿesti 542/2 , 16610 Prague , Czech Republic.,Computational Physics Laboratory , Tampere University , P.O. Box 692, FI-33014 Tampere , Finland
| | - Waldemar Kulig
- Department of Physics , University of Helsinki , P.O. Box 64, FI-00014 Helsinki , Finland
| | - Tomasz Róg
- Department of Physics , University of Helsinki , P.O. Box 64, FI-00014 Helsinki , Finland.,Computational Physics Laboratory , Tampere University , P.O. Box 692, FI-33014 Tampere , Finland
| | - Ilpo Vattulainen
- Department of Physics , University of Helsinki , P.O. Box 64, FI-00014 Helsinki , Finland.,Computational Physics Laboratory , Tampere University , P.O. Box 692, FI-33014 Tampere , Finland.,MEMPHYS-Center for Biomembrane Physics
| |
Collapse
|
45
|
Wang J, Miao Y. Recent advances in computational studies of GPCR-G protein interactions. ADVANCES IN PROTEIN CHEMISTRY AND STRUCTURAL BIOLOGY 2019; 116:397-419. [PMID: 31036298 PMCID: PMC6986689 DOI: 10.1016/bs.apcsb.2018.11.011] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/14/2022]
Abstract
Protein-protein interactions are key in cellular signaling. G protein-coupled receptors (GPCRs), the largest superfamily of human membrane proteins, are able to transduce extracellular signals (e.g., hormones and neurotransmitters) to intracellular proteins, in particular the G proteins. Since GPCRs serve as primary targets of ~1/3 of currently marketed drugs, it is important to understand mechanisms of GPCR signaling in order to design selective and potent drug molecules. This chapter focuses on recent advances in computational studies of the GPCR-G protein interactions using bioinformatics, protein-protein docking and molecular dynamics simulation approaches.
Collapse
Affiliation(s)
- Jinan Wang
- Center for Computational Biology and Department of Molecular Biosciences, University of Kansas, Lawrence, KS, United States
| | - Yinglong Miao
- Center for Computational Biology and Department of Molecular Biosciences, University of Kansas, Lawrence, KS, United States.
| |
Collapse
|
46
|
Lambrughi M, Tiberti M, Allega MF, Sora V, Nygaard M, Toth A, Salamanca Viloria J, Bignon E, Papaleo E. Analyzing Biomolecular Ensembles. Methods Mol Biol 2019; 2022:415-451. [PMID: 31396914 DOI: 10.1007/978-1-4939-9608-7_18] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/10/2023]
Abstract
Several techniques are available to generate conformational ensembles of proteins and other biomolecules either experimentally or computationally. These methods produce a large amount of data that need to be analyzed to identify structure-dynamics-function relationship. In this chapter, we will cover different tools to unveil the information hidden in conformational ensemble data and to guide toward the rationalization of the data. We included routinely used approaches such as dimensionality reduction, as well as new methods inspired by high-order statistics and graph theory.
Collapse
Affiliation(s)
- Matteo Lambrughi
- Computational Biology Laboratory, Danish Cancer Society Research Center, Copenhagen, Denmark
| | - Matteo Tiberti
- Computational Biology Laboratory, Danish Cancer Society Research Center, Copenhagen, Denmark
| | - Maria Francesca Allega
- Computational Biology Laboratory, Danish Cancer Society Research Center, Copenhagen, Denmark
| | - Valentina Sora
- Computational Biology Laboratory, Danish Cancer Society Research Center, Copenhagen, Denmark
| | - Mads Nygaard
- Computational Biology Laboratory, Danish Cancer Society Research Center, Copenhagen, Denmark
| | - Agota Toth
- Computational Biology Laboratory, Danish Cancer Society Research Center, Copenhagen, Denmark
| | - Juan Salamanca Viloria
- Computational Biology Laboratory, Danish Cancer Society Research Center, Copenhagen, Denmark
| | - Emmanuelle Bignon
- Computational Biology Laboratory, Danish Cancer Society Research Center, Copenhagen, Denmark
| | - Elena Papaleo
- Computational Biology Laboratory, Danish Cancer Society Research Center, Copenhagen, Denmark.
| |
Collapse
|
47
|
Miao Y, Bhattarai A, Nguyen ATN, Christopoulos A, May LT. Structural Basis for Binding of Allosteric Drug Leads in the Adenosine A 1 Receptor. Sci Rep 2018; 8:16836. [PMID: 30442899 PMCID: PMC6237911 DOI: 10.1038/s41598-018-35266-x] [Citation(s) in RCA: 50] [Impact Index Per Article: 8.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/10/2018] [Accepted: 11/02/2018] [Indexed: 12/12/2022] Open
Abstract
Despite intense interest in designing positive allosteric modulators (PAMs) as selective drugs of the adenosine A1 receptor (A1AR), structural binding modes of the receptor PAMs remain unknown. Using the first X-ray structure of the A1AR, we have performed all-atom simulations using a robust Gaussian accelerated molecular dynamics (GaMD) technique to determine binding modes of the A1AR allosteric drug leads. Two prototypical PAMs, PD81723 and VCP171, were selected. Each PAM was initially placed at least 20 Å away from the receptor. Extensive GaMD simulations using the AMBER and NAMD simulation packages at different acceleration levels captured spontaneous binding of PAMs to the A1AR. The simulations allowed us to identify low-energy binding modes of the PAMs at an allosteric site formed by the receptor extracellular loop 2 (ECL2), which are highly consistent with mutagenesis experimental data. Furthermore, the PAMs stabilized agonist binding in the receptor. In the absence of PAMs at the ECL2 allosteric site, the agonist sampled a significantly larger conformational space and even dissociated from the A1AR alone. In summary, the GaMD simulations elucidated structural binding modes of the PAMs and provided important insights into allostery in the A1AR, which will greatly facilitate the receptor structure-based drug design.
Collapse
Affiliation(s)
- Yinglong Miao
- Center for Computational Biology and Department of Molecular Biosciences, University of Kansas, Lawrence, KS, 66047, USA.
| | - Apurba Bhattarai
- Center for Computational Biology and Department of Molecular Biosciences, University of Kansas, Lawrence, KS, 66047, USA
| | - Anh T N Nguyen
- Drug Discovery Biology, Monash Institute of Pharmaceutical Sciences and Department of Pharmacology, Monash University, Parkville, VIC, 3052, Australia
| | - Arthur Christopoulos
- Drug Discovery Biology, Monash Institute of Pharmaceutical Sciences and Department of Pharmacology, Monash University, Parkville, VIC, 3052, Australia
| | - Lauren T May
- Drug Discovery Biology, Monash Institute of Pharmaceutical Sciences and Department of Pharmacology, Monash University, Parkville, VIC, 3052, Australia
| |
Collapse
|
48
|
Abstract
INTRODUCTION Understanding pathways and mechanisms of drug binding to receptors is important for rational drug design. Remarkable advances in supercomputing and methodological developments have opened a new era for application of computer simulations in predicting drug-receptor interactions at an atomistic level. Gaussian accelerated molecular dynamics (GaMD) is a computational enhanced sampling technique that works by adding a harmonic boost potential to reduce energy barriers. GaMD enables free energy calculations without the requirement of predefined collective variables. GaMD has proven useful in biomolecular simulations, in particular, the prediction of drug-receptor interactions. Areas covered: Herein, the authors review recent GaMD simulation studies that elucidated pathways of drug binding to proteins including the G-protein-coupled receptors and HIV protease. Expert opinion: GaMD is advantageous for enhanced simulations of, amongst many biological processes, drug binding to target receptors. Compared with conventional molecular dynamics, GaMD speeds up biomolecular simulations by orders of magnitude. GaMD enables routine drug binding simulations using personal computers with GPUs or common computing clusters. GaMD and, more broadly, enhanced sampling simulations are expected to dramatically increase our capabilities to determine the mechanisms of drug binding to a wide range of receptors in the near future. This will greatly facilitate computer-aided drug design.
Collapse
Affiliation(s)
- Apurba Bhattarai
- Center for Computational Biology and Department of Molecular Biosciences, University of Kansas, Lawrence, KS 66047, USA,
| | - Yinglong Miao
- Center for Computational Biology and Department of Molecular Biosciences, University of Kansas, Lawrence, KS 66047, USA,
| |
Collapse
|
49
|
Miao Y. Acceleration of biomolecular kinetics in Gaussian accelerated molecular dynamics. J Chem Phys 2018; 149:072308. [PMID: 30134710 PMCID: PMC6901173 DOI: 10.1063/1.5024217] [Citation(s) in RCA: 30] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/04/2023] Open
Abstract
Recent studies demonstrated that Gaussian accelerated molecular dynamics (GaMD) is a robust computational technique, which provides simultaneous unconstrained enhanced sampling and free energy calculations of biomolecules. However, the exact acceleration of biomolecular dynamics or speedup of kinetic rates in GaMD simulations and, more broadly, in enhanced sampling methods, remains a challenging task to be determined. Here, the GaMD acceleration is examined using alanine dipeptide in explicit solvent as a biomolecular model system. Relative to long conventional molecular dynamics simulation, GaMD simulations exhibited ∼36-67 times speedup for sampling of the backbone dihedral transitions. The acceleration depended on level of the GaMD boost potential. Furthermore, Kramers' rate theory was applied to estimate GaMD acceleration using simulation-derived diffusion coefficients, curvatures and barriers of free energy profiles. In most cases, the calculations also showed significant speedup of dihedral transitions in GaMD, although the GaMD acceleration factors tended to be underestimated by ∼3-96 fold. Because greater boost potential can be applied in GaMD simulations of systems with increased sizes, which potentially leads to higher acceleration, it is subject to future studies on accelerating the dynamics and recovering kinetic rates of larger biomolecules such as proteins and protein-protein/nucleic acid complexes.
Collapse
Affiliation(s)
- Yinglong Miao
- Center for Computational Biology and Department of Molecular Biosciences, University of Kansas, Lawrence, Kansas 66047, USA
| |
Collapse
|
50
|
Ferreira AF, Comune M, Rai A, Ferreira L, Simões PN. Atomistic-Level Investigation of a LL37-Conjugated Gold Nanoparticle By Well-Tempered Metadynamics. J Phys Chem B 2018; 122:8359-8366. [PMID: 30106592 DOI: 10.1021/acs.jpcb.8b05717] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/06/2023]
Abstract
LL37 is a cathelicidin-derived antimicrobial peptide (AMP) with a broad spectrum of antimicrobial activity and wound-healing potential. The enhancement of these characteristics was recently demonstrated for a cysteine (CYS)-modified cathelicidin-derived LL37-SH conjugated with gold nanoparticles (AuNPs). Considering the potential of this peptide, we hereby report a computational study in which well-tempered metadynamics was applied to unveil the interaction of LL37-SH and LL37 with a AuNP with atomistic detail. A structural analysis combined with the free energy surface (FES) characterization allowed the assessment of the role of CYS residue during the formation of the conjugate, as well as to understand how the AuNP improves the antimicrobial activity of the peptide. It was found that CYS promotes a lower conformational entropy (before and after adsorption onto the AuNP) and a faster adsorption process when compared to the LL37 without CYS. The FES for LL37-SH is characterized by one global minimum, while for LL37 a potential metastable state was found. The presence of the AuNP leads to an elongation of the peptides along with the adsorption, which translates into the increase of the solvent-accessible surface area. This elongation combined with the greater availability of positively charged residues upon adsorption rationalizes the observed enhancement of the activity of the LL37-SH/AuNP conjugate.
Collapse
Affiliation(s)
- André F Ferreira
- CIEPQPF, Department of Chemical Engineering , University of Coimbra , 3030-790 Coimbra , Portugal
| | - Michela Comune
- CNC-Center for Neurosciences and Cell Biology , University of Coimbra , 3004-517 Coimbra , Portugal.,Biocant , Biotechnology Innovation Center , 3060-197 Cantanhede , Portugal
| | - Akhilesh Rai
- CNC-Center for Neurosciences and Cell Biology , University of Coimbra , 3004-517 Coimbra , Portugal.,Biocant , Biotechnology Innovation Center , 3060-197 Cantanhede , Portugal
| | - Lino Ferreira
- CNC-Center for Neurosciences and Cell Biology , University of Coimbra , 3004-517 Coimbra , Portugal.,Biocant , Biotechnology Innovation Center , 3060-197 Cantanhede , Portugal
| | - Pedro N Simões
- CIEPQPF, Department of Chemical Engineering , University of Coimbra , 3030-790 Coimbra , Portugal
| |
Collapse
|