1
|
Pyzer-Knapp EO, Curioni A. Advancing biomolecular simulation through exascale HPC, AI and quantum computing. Curr Opin Struct Biol 2024; 87:102826. [PMID: 38733863 DOI: 10.1016/j.sbi.2024.102826] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/15/2023] [Revised: 04/11/2024] [Accepted: 04/11/2024] [Indexed: 05/13/2024]
Abstract
Biomolecular simulation can act as both a digital microscope and a crystal ball; offering the potential for a deeper understanding of experimental observations whilst also presenting a forward-looking avenue for the in silico design and evaluation of hitherto unsynthesized compounds. Indeed, as the intricacy of our scientific inquiries has grown, so too has the computational prowess we seek to deploy in our pursuit of answers. As we enter the Exascale era, this mini-review surveys the computational landscape from both the point of view of the development of new and ever more powerful systems, and the simulations that are run on them. Moreover, as we stand on the cusp of a transformative phase in computational biology, this article offers a contemplative glance into the future, speculating on the profound implications of artificial intelligence and quantum computing for large-scale biomolecular simulations.
Collapse
|
2
|
Asgharpour S, Chi LA, Spehr M, Carloni P, Alfonso-Prieto M. Fluoride Transport and Inhibition Across CLC Transporters. Handb Exp Pharmacol 2024; 283:81-100. [PMID: 36042142 DOI: 10.1007/164_2022_593] [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: 12/29/2022]
Abstract
The Chloride Channel (CLC) family includes proton-coupled chloride and fluoride transporters. Despite their similar protein architecture, the former exchange two chloride ions for each proton and are inhibited by fluoride, whereas the latter efficiently transport one fluoride in exchange for one proton. The combination of structural, mutagenesis, and functional experiments with molecular simulations has pinpointed several amino acid changes in the permeation pathway that capitalize on the different chemical properties of chloride and fluoride to fine-tune protein function. Here we summarize recent findings on fluoride inhibition and transport in the two prototypical members of the CLC family, the chloride/proton transporter from Escherichia coli (CLC-ec1) and the fluoride/proton transporter from Enterococcus casseliflavus (CLCF-eca).
Collapse
Affiliation(s)
- Somayeh Asgharpour
- Institute for Advanced Simulations IAS-5 and Institute of Neuroscience and Medicine INM-9, Computational Biomedicine, Forschungszentrum Jülich, Jülich, Germany
- Research Training Group 2416 MultiSenses-MultiScales, Institute for Biology II, RWTH Aachen University, Aachen, Germany
| | - L América Chi
- Laboratory for the Design and Development of New Drugs and Biotechnological Innovation, Escuela Superior de Medicina, Instituto Politécnico Nacional, Plan de San Luis y Díaz Mirón, Ciudad de México, Mexico
| | - Marc Spehr
- Research Training Group 2416 MultiSenses-MultiScales, Institute for Biology II, RWTH Aachen University, Aachen, Germany
- Department of Chemosensation, Institute for Biology II, RWTH Aachen University, Aachen, Germany
| | - Paolo Carloni
- Institute for Advanced Simulations IAS-5 and Institute of Neuroscience and Medicine INM-9, Computational Biomedicine, Forschungszentrum Jülich, Jülich, Germany.
- Research Training Group 2416 MultiSenses-MultiScales, Institute for Biology II, RWTH Aachen University, Aachen, Germany.
- Department of Physics, RWTH Aachen University, Aachen, Germany.
- JARA Institute Molecular Neuroscience and Neuroimaging (INM-11), Forschungszentrum Jülich, Jülich, Germany.
- JARA-HPC, Forschungszentrum Jülich, Jülich, Germany.
| | - Mercedes Alfonso-Prieto
- Institute for Advanced Simulations IAS-5 and Institute of Neuroscience and Medicine INM-9, Computational Biomedicine, Forschungszentrum Jülich, Jülich, Germany.
- Medical Faculty, Cécile and Oskar Vogt Institute for Brain Research, University Hospital Düsseldorf, Heinrich Heine University Düsseldorf, Düsseldorf, Germany.
| |
Collapse
|
3
|
Sabei A, Prentiss M, Prévost C. Modeling the Homologous Recombination Process: Methods, Successes and Challenges. Int J Mol Sci 2023; 24:14896. [PMID: 37834348 PMCID: PMC10573387 DOI: 10.3390/ijms241914896] [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: 09/04/2023] [Revised: 09/24/2023] [Accepted: 09/25/2023] [Indexed: 10/15/2023] Open
Abstract
Homologous recombination (HR) is a fundamental process common to all species. HR aims to faithfully repair DNA double strand breaks. HR involves the formation of nucleoprotein filaments on DNA single strands (ssDNA) resected from the break. The nucleoprotein filaments search for homologous regions in the genome and promote strand exchange with the ssDNA homologous region in an unbroken copy of the genome. HR has been the object of intensive studies for decades. Because multi-scale dynamics is a fundamental aspect of this process, studying HR is highly challenging, both experimentally and using computational approaches. Nevertheless, knowledge has built up over the years and has recently progressed at an accelerated pace, borne by increasingly focused investigations using new techniques such as single molecule approaches. Linking this knowledge to the atomic structure of the nucleoprotein filament systems and the succession of unstable, transient intermediate steps that takes place during the HR process remains a challenge; modeling retains a very strong role in bridging the gap between structures that are stable enough to be observed and in exploring transition paths between these structures. However, working on ever-changing long filament systems submitted to kinetic processes is full of pitfalls. This review presents the modeling tools that are used in such studies, their possibilities and limitations, and reviews the advances in the knowledge of the HR process that have been obtained through modeling. Notably, we will emphasize how cooperative behavior in the HR nucleoprotein filament enables modeling to produce reliable information.
Collapse
Affiliation(s)
- Afra Sabei
- CNRS, UPR 9080, Laboratoire de Biochimie Théorique, Université de Paris, 13 Rue Pierre et Marie Curie, F-75005 Paris, France;
- Institut de Biologie Physico-Chimique-Fondation Edmond de Rotschild, PSL Research University, F-75005 Paris, France
| | - Mara Prentiss
- Department of Physics, Harvard University, Cambridge, MA02138, USA;
| | - Chantal Prévost
- CNRS, UPR 9080, Laboratoire de Biochimie Théorique, Université de Paris, 13 Rue Pierre et Marie Curie, F-75005 Paris, France;
- Institut de Biologie Physico-Chimique-Fondation Edmond de Rotschild, PSL Research University, F-75005 Paris, France
| |
Collapse
|
4
|
Durumeric AEP, Voth GA. Using classifiers to understand coarse-grained models and their fidelity with the underlying all-atom systems. J Chem Phys 2023; 158:234101. [PMID: 37318166 DOI: 10.1063/5.0146812] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/16/2023] [Accepted: 05/26/2023] [Indexed: 06/16/2023] Open
Abstract
Bottom-up coarse-grained (CG) molecular dynamics models are parameterized using complex effective Hamiltonians. These models are typically optimized to approximate high dimensional data from atomistic simulations. However, human validation of these models is often limited to low dimensional statistics that do not necessarily differentiate between the CG model and said atomistic simulations. We propose that classification can be used to variationally estimate high dimensional error and that explainable machine learning can help convey this information to scientists. This approach is demonstrated using Shapley additive explanations and two CG protein models. This framework may also be valuable for ascertaining whether allosteric effects at the atomistic level are accurately propagated to a CG model.
Collapse
Affiliation(s)
- Aleksander E P Durumeric
- Department of Chemistry, Chicago Center for Theoretical Chemistry, James Franck Institute, and Institute for Biophysical Dynamics, The University of Chicago, 5735 S. Ellis Ave., Chicago, Illinois 60637, USA
| | - Gregory A Voth
- Department of Chemistry, Chicago Center for Theoretical Chemistry, James Franck Institute, and Institute for Biophysical Dynamics, The University of Chicago, 5735 S. Ellis Ave., Chicago, Illinois 60637, USA
| |
Collapse
|
5
|
Zhang J, Wang Z, Zhuang W, Rabiee H, Zhu C, Deng J, Ge L, Ying H. Amphiphilic Nanointerface: Inducing the Interfacial Activation for Lipase. ACS APPLIED MATERIALS & INTERFACES 2022; 14:39622-39636. [PMID: 35980131 DOI: 10.1021/acsami.2c11500] [Citation(s) in RCA: 4] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/15/2023]
Abstract
Graphene-based materials are widely used in the field of immobilized enzymes due to their easily tunable interfacial properties. We designed amphiphilic nanobiological interfaces between graphene oxide (GO) and lipase TL (Thermomyces lanuginosus) with tunable reduction degrees through molecular dynamics simulations and a facile chemical modulation, thus revealing the optimal interface for the interfacial activation of lipase TL and addressing the weakness of lipase TL, which exhibits weak catalytic activity due to an inconspicuous active site lid. It was demonstrated that the reduced graphene oxide (rGO) after 4 h of ascorbic acid reduction could boost the relative enzyme activity of lipase TL to reach 208%, which was 48% higher than the pristine GO and 120% higher than the rGO after 48 h of reduction. Moreover, TL-GO-4 h's tolerance against heat, organic solvent, and long-term storage environment was higher than that of free TL. The drawbacks of strong hydrophobic nanomaterials on lipase production were explored in depth with the help of molecular dynamics simulations, which explained the mechanism of enzyme activity enhancement. We demonstrated that nanomaterials with certain hydrophilicity could facilitate the lipase to undergo interfacial activation and improve its stability and protein loading rate, displaying the potential of the extensive application.
Collapse
Affiliation(s)
- Jihang Zhang
- College of Biotechnology and Pharmaceutical Engineering, Nanjing Tech University, No. 30, Puzhu South Road, Nanjing 211816, China
| | - Zhaoxin Wang
- School of Chemical Engineering, Zhengzhou University, Zhengzhou 450001, China
| | - Wei Zhuang
- College of Biotechnology and Pharmaceutical Engineering, Nanjing Tech University, No. 30, Puzhu South Road, Nanjing 211816, China
- State Key Laboratory of Materials-Oriented Chemical Engineering, National Engineering Technique Research Center for Biotechnology, Nanjing Tech University, No. 30, Puzhu South Road, Nanjing 211816, China
- School of Chemical Engineering, Zhengzhou University, Zhengzhou 450001, China
| | - Hesamoddin Rabiee
- Centre for Future Materials, University of Southern Queensland, Springfield, QLD 4300, Australia
- School of Chemical Engineering, The University of Queensland, Brisbane, QLD 4072, Australia
| | - Chenjie Zhu
- College of Biotechnology and Pharmaceutical Engineering, Nanjing Tech University, No. 30, Puzhu South Road, Nanjing 211816, China
- State Key Laboratory of Materials-Oriented Chemical Engineering, National Engineering Technique Research Center for Biotechnology, Nanjing Tech University, No. 30, Puzhu South Road, Nanjing 211816, China
| | - Jiawei Deng
- School of Chemical Engineering, Zhengzhou University, Zhengzhou 450001, China
| | - Lei Ge
- Centre for Future Materials, University of Southern Queensland, Springfield, QLD 4300, Australia
- School of Chemical Engineering, The University of Queensland, Brisbane, QLD 4072, Australia
| | - Hanjie Ying
- State Key Laboratory of Materials-Oriented Chemical Engineering, National Engineering Technique Research Center for Biotechnology, Nanjing Tech University, No. 30, Puzhu South Road, Nanjing 211816, China
- School of Chemical Engineering, Zhengzhou University, Zhengzhou 450001, China
| |
Collapse
|
6
|
Ion Mobility Mass Spectrometry for Structural Elucidation of Petroleum Compounds. Trends Analyt Chem 2022. [DOI: 10.1016/j.trac.2022.116597] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/19/2022]
|
7
|
Padhi AK, Rath SL, Tripathi T. Accelerating COVID-19 Research Using Molecular Dynamics Simulation. J Phys Chem B 2021; 125:9078-9091. [PMID: 34319118 PMCID: PMC8340580 DOI: 10.1021/acs.jpcb.1c04556] [Citation(s) in RCA: 39] [Impact Index Per Article: 13.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/24/2021] [Revised: 07/12/2021] [Indexed: 12/14/2022]
Abstract
The COVID-19 pandemic has emerged as a global medico-socio-economic disaster. Given the lack of effective therapeutics against SARS-CoV-2, scientists are racing to disseminate suggestions for rapidly deployable therapeutic options, including drug repurposing and repositioning strategies. Molecular dynamics (MD) simulations have provided the opportunity to make rational scientific breakthroughs in a time of crisis. Advancements in these technologies in recent years have become an indispensable tool for scientists studying protein structure, function, dynamics, interactions, and drug discovery. Integrating the structural data obtained from high-resolution methods with MD simulations has helped in comprehending the process of infection and pathogenesis, as well as the SARS-CoV-2 maturation in host cells, in a short duration of time. It has also guided us to identify and prioritize drug targets and new chemical entities, and to repurpose drugs. Here, we discuss how MD simulation has been explored by the scientific community to accelerate and guide translational research on SARS-CoV-2 in the past year. We have also considered future research directions for researchers, where MD simulations can help fill the existing gaps in COVID-19 research.
Collapse
Affiliation(s)
- Aditya K. Padhi
- Laboratory for Structural Bioinformatics, Center for
Biosystems Dynamics Research, RIKEN, 1-7-22 Suehiro, Tsurumi,
Yokohama, Kanagawa 230-0045, Japan
| | - Soumya Lipsa Rath
- Department of Biotechnology, National
Institute of Technology, Warangal, Telangana 506004,
India
| | - Timir Tripathi
- Molecular and Structural Biophysics Laboratory,
Department of Biochemistry, North-Eastern Hill University,
Shillong 793022, India
| |
Collapse
|
8
|
Schlick T, Portillo-Ledesma S, Myers CG, Beljak L, Chen J, Dakhel S, Darling D, Ghosh S, Hall J, Jan M, Liang E, Saju S, Vohr M, Wu C, Xu Y, Xue E. Biomolecular Modeling and Simulation: A Prospering Multidisciplinary Field. Annu Rev Biophys 2021; 50:267-301. [PMID: 33606945 PMCID: PMC8105287 DOI: 10.1146/annurev-biophys-091720-102019] [Citation(s) in RCA: 24] [Impact Index Per Article: 8.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022]
Abstract
We reassess progress in the field of biomolecular modeling and simulation, following up on our perspective published in 2011. By reviewing metrics for the field's productivity and providing examples of success, we underscore the productive phase of the field, whose short-term expectations were overestimated and long-term effects underestimated. Such successes include prediction of structures and mechanisms; generation of new insights into biomolecular activity; and thriving collaborations between modeling and experimentation, including experiments driven by modeling. We also discuss the impact of field exercises and web games on the field's progress. Overall, we note tremendous success by the biomolecular modeling community in utilization of computer power; improvement in force fields; and development and application of new algorithms, notably machine learning and artificial intelligence. The combined advances are enhancing the accuracy andscope of modeling and simulation, establishing an exemplary discipline where experiment and theory or simulations are full partners.
Collapse
Affiliation(s)
- Tamar Schlick
- Department of Chemistry, New York University, New York, New York 10003, USA;
- Courant Institute of Mathematical Sciences, New York University, New York, New York 10012, USA
- New York University-East China Normal University Center for Computational Chemistry, New York University Shanghai, Shanghai 200122, China
| | | | - Christopher G Myers
- Department of Chemistry, New York University, New York, New York 10003, USA;
| | - Lauren Beljak
- College of Arts and Science, New York University, New York, New York 10003, USA
| | - Justin Chen
- College of Arts and Science, New York University, New York, New York 10003, USA
| | - Sami Dakhel
- College of Arts and Science, New York University, New York, New York 10003, USA
| | - Daniel Darling
- College of Arts and Science, New York University, New York, New York 10003, USA
| | - Sayak Ghosh
- College of Arts and Science, New York University, New York, New York 10003, USA
| | - Joseph Hall
- College of Arts and Science, New York University, New York, New York 10003, USA
| | - Mikaeel Jan
- College of Arts and Science, New York University, New York, New York 10003, USA
| | - Emily Liang
- College of Arts and Science, New York University, New York, New York 10003, USA
| | - Sera Saju
- College of Arts and Science, New York University, New York, New York 10003, USA
| | - Mackenzie Vohr
- College of Arts and Science, New York University, New York, New York 10003, USA
| | - Chris Wu
- College of Arts and Science, New York University, New York, New York 10003, USA
| | - Yifan Xu
- College of Arts and Science, New York University, New York, New York 10003, USA
| | - Eva Xue
- College of Arts and Science, New York University, New York, New York 10003, USA
| |
Collapse
|
9
|
Fossépré M, Tuvi-Arad I, Beljonne D, Richeter S, Clément S, Surin M. Binding Mode Multiplicity and Multiscale Chirality in the Supramolecular Assembly of DNA and a π-Conjugated Polymer. Chemphyschem 2020; 21:2543-2552. [PMID: 32910539 DOI: 10.1002/cphc.202000630] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/17/2020] [Revised: 09/04/2020] [Indexed: 02/04/2023]
Abstract
Water-soluble π-conjugated polymers are increasingly considered for DNA biosensing. However, the conformational rearrangement, supramolecular organization and dynamics upon interaction with DNA have been overlooked, which prevents the rational design of such detection tools. To elucidate the binding of a cationic polythiophene (CPT) to DNA with atomistic resolution, we performed molecular simulations of their supramolecular assembly. Comparison of replicated simulations show a multiplicity of CPT binding geometries that contribute to the wrapping of CPT around DNA. The different binding geometries are stabilized by both electrostatic interactions between CPT lateral cations and DNA phosphodiesters and van der Waals interactions between the CPT backbone and the DNA grooves. Simulated circular dichroism (CD) spectra show that the induced CD signal stems from a conserved geometrical feature across the replicated simulations, i. e. the presence of segments of syn configurations between thiophene units along the CPT chain. At the macromolecular scale, we inspected the different shapes related to the CPT binding modes around the DNA through symmetry metrics. Altogether, molecular dynamics (MD) simulations, model Hamiltonian calculations of the CD spectra, and symmetry indices provide insights into the origin of induced chirality from the atomic to the macromolecular scale. Our multidisciplinary approach points out the hierarchical aspect of CPT chiral organization induced by DNA.
Collapse
Affiliation(s)
- Mathieu Fossépré
- Laboratory for Chemistry of Novel Materials, Centre of Innovation and Research in Materials and Polymers (CIRMAP), University of Mons - UMONS, 20 Place du Parc, Mons, 7000, Belgium
| | - Inbal Tuvi-Arad
- Department of Natural Sciences, The Open University of Israel, Raanana, Israel
| | - David Beljonne
- Laboratory for Chemistry of Novel Materials, Centre of Innovation and Research in Materials and Polymers (CIRMAP), University of Mons - UMONS, 20 Place du Parc, Mons, 7000, Belgium
| | | | | | - Mathieu Surin
- Laboratory for Chemistry of Novel Materials, Centre of Innovation and Research in Materials and Polymers (CIRMAP), University of Mons - UMONS, 20 Place du Parc, Mons, 7000, Belgium
| |
Collapse
|
10
|
Matsunaga Y, Sugita Y. Use of single-molecule time-series data for refining conformational dynamics in molecular simulations. Curr Opin Struct Biol 2020; 61:153-159. [DOI: 10.1016/j.sbi.2019.12.022] [Citation(s) in RCA: 9] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/30/2019] [Revised: 12/24/2019] [Accepted: 12/27/2019] [Indexed: 12/18/2022]
|
11
|
Durumeric AEP, Voth GA. Adversarial-residual-coarse-graining: Applying machine learning theory to systematic molecular coarse-graining. J Chem Phys 2019; 151:124110. [PMID: 31575201 DOI: 10.1063/1.5097559] [Citation(s) in RCA: 18] [Impact Index Per Article: 3.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/06/2023] Open
Abstract
We utilize connections between molecular coarse-graining (CG) approaches and implicit generative models in machine learning to describe a new framework for systematic molecular CG. Focus is placed on the formalism encompassing generative adversarial networks. The resulting method enables a variety of model parameterization strategies, some of which show similarity to previous CG methods. We demonstrate that the resulting framework can rigorously parameterize CG models containing CG sites with no prescribed connection to the reference atomistic system (termed virtual sites); however, this advantage is offset by the lack of a closed-form expression for the CG Hamiltonian at the resolution obtained after integration over the virtual CG sites. Computational examples are provided for cases in which these methods ideally return identical parameters as relative entropy minimization CG but where traditional relative entropy minimization CG optimization equations are not applicable.
Collapse
Affiliation(s)
- Aleksander E P Durumeric
- Department of Chemistry, James Franck Institute, Institute for Biophysical Dynamics, and Computation Institute, The University of Chicago, Chicago, Illinois 60637, USA
| | - Gregory A Voth
- Department of Chemistry, James Franck Institute, Institute for Biophysical Dynamics, and Computation Institute, The University of Chicago, Chicago, Illinois 60637, USA
| |
Collapse
|
12
|
Lim D, Park Y, Chang R, Ahmed A, Kim S. Application of molecular dynamics simulation to improve the theoretical prediction for collisional cross section of aromatic compounds with long alkyl chains in crude oils. RAPID COMMUNICATIONS IN MASS SPECTROMETRY : RCM 2019; 33:650-656. [PMID: 30710409 DOI: 10.1002/rcm.8400] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 10/15/2018] [Revised: 01/21/2019] [Accepted: 01/23/2019] [Indexed: 06/09/2023]
Abstract
RATIONALE Molecular dynamics (MD) simulations with finite temperature were performed to improve the theoretical prediction of collisional cross section (CCS) values, especially for aromatic compounds containing long alkyl chains. METHODS In this study, the CCS values of 11 aromatic compounds with long alkyl chains were calculated by MD simulations while considering internal energy at 300, 500, and 700 K, and the results were compared with experimentally determined values. RESULTS The CCS values calculated at higher energies showed better agreement with the experimental values. Polycyclic aromatic hydrocarbons (PAHs) such as pentacene and benz[b]anthracene were also investigated, and better agreement between the theoretical and experimental results was observed when higher temperature (or higher internal energy) was considered. CONCLUSIONS The data presented in this study show that the internal degrees of freedom of ions must be considered to accurately predict the CCS values of aromatic compounds with a flexible structure measured by ion mobility mass spectrometry.
Collapse
Affiliation(s)
- Dongwan Lim
- Department of Chemistry, Kyungpook National University, Daegu, 41566, Republic of Korea
| | - Yunjae Park
- Department of Chemistry, Kwangwoon University, Seoul, 01897, Republic of Korea
| | - Rakwoo Chang
- Department of Chemistry, Kwangwoon University, Seoul, 01897, Republic of Korea
| | - Arif Ahmed
- Department of Chemistry, Kyungpook National University, Daegu, 41566, Republic of Korea
| | - Sunghwan Kim
- Department of Chemistry, Kyungpook National University, Daegu, 41566, Republic of Korea
- Green Nano Center, Department of Chemistry, Daegu, 41566, Republic of Korea
| |
Collapse
|
13
|
Nilsson L, Villa A. Modeling and Simulation of Oligonucleotide Hybrids: Outlining a Strategy. Methods Mol Biol 2019; 2036:113-126. [PMID: 31410793 DOI: 10.1007/978-1-4939-9670-4_6] [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
Molecular dynamics simulations with a state-of-the-art force field provide an atomistic detailed description of the structural and thermodynamic features of biomolecules. Effects of chemical modifications and of the environment such as sequence, solvent, and ionic strength can explicitly be taken into account. Molecular simulation techniques can also provide insight in change in binding affinity, in protonation (pKa shift) and tautomeric propensity due to changes in the environment or in the molecular system. The quality and reliability of a simulation depend on the quality of the force field and on the reproducibility of the data, and validation depends on the availability of suitable experimental data. Here, we describe the workflow to investigate oligonucleotide hybrids using molecular simulation including hardware and software information.
Collapse
Affiliation(s)
- Lennart Nilsson
- Department of Biosciences and Nutrition, Karolinska Institutet, Huddinge, Sweden
| | - Alessandra Villa
- Department of Biosciences and Nutrition, Karolinska Institutet, Huddinge, Sweden.
| |
Collapse
|
14
|
Todorova N, Yarovsky I. The Enigma of Amyloid Forming Proteins: Insights From Molecular Simulations. Aust J Chem 2019. [DOI: 10.1071/ch19059] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022]
Abstract
Molecular level insight into the interplay between protein sequence, structure, and conformational dynamics is crucial for the comprehensive understanding of protein folding, misfolding, and aggregation phenomena that are pertinent to the formation of amyloid fibrils implicated in several degenerative diseases. Computational modelling provides insight into protein behaviour at spatial and temporal resolution still largely outside the reach of experiments. Herein we present an account of our theoretical modelling research conducted in collaboration with several experimental groups where we explored the effects of local environment on the structure and aggregation propensity of several types of amyloidogenic peptides and proteins, including apolipoprotein C-II, insulin, amylin, and amyloid-β using a variety of computational approaches.
Collapse
|
15
|
Dutta R, Brotzakis ZF, Mira A. Bayesian calibration of force-fields from experimental data: TIP4P water. J Chem Phys 2018; 149:154110. [DOI: 10.1063/1.5030950] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/30/2023] Open
Affiliation(s)
- Ritabrata Dutta
- Institute of Computational Science, Università della Svizzera Italiana, Lugano, Switzerland
| | - Zacharias Faidon Brotzakis
- Institute of Computational Science, Università della Svizzera Italiana, Lugano, Switzerland
- Department of Chemistry and Applied Bioscience, ETH Zürich, Zürich, Switzerland
| | - Antonietta Mira
- Institute of Computational Science, Università della Svizzera Italiana, Lugano, Switzerland
- Department of Science and High Technology, Università degli Studi dell’Insubria, Varese, Italy
| |
Collapse
|
16
|
Petrović D, Wang X, Strodel B. How accurately do force fields represent protein side chain ensembles? Proteins 2018; 86:935-944. [DOI: 10.1002/prot.25525] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/18/2017] [Revised: 05/02/2018] [Accepted: 05/14/2018] [Indexed: 12/22/2022]
Affiliation(s)
- Dušan Petrović
- Institute of Complex Systems, Structural Biochemistry, Forschungszentrum Jülich; Jülich, 52425 Germany
- Department of Cell and Molecular Biology; Uppsala University, BMC Box 596; Uppsala, 751 24 Sweden
| | - Xue Wang
- Institute of Complex Systems, Structural Biochemistry, Forschungszentrum Jülich; Jülich, 52425 Germany
- Institute of Theoretical and Computational Chemistry; Heinrich Heine University Düsseldorf, Universitätsstraße 1; Düsseldorf, 40225 Germany
| | - Birgit Strodel
- Institute of Complex Systems, Structural Biochemistry, Forschungszentrum Jülich; Jülich, 52425 Germany
- Institute of Theoretical and Computational Chemistry; Heinrich Heine University Düsseldorf, Universitätsstraße 1; Düsseldorf, 40225 Germany
| |
Collapse
|
17
|
Baz J, Gebhardt J, Kraus H, Markthaler D, Hansen N. Insights into Noncovalent Binding Obtained from Molecular Dynamics Simulations. CHEM-ING-TECH 2018. [DOI: 10.1002/cite.201800050] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/11/2022]
Affiliation(s)
- Jörg Baz
- University of Stuttgart; Institute of Thermodynamics and Thermal Process Engineering; Pfaffenwaldring 9 70569 Stuttgart Germany
| | - Julia Gebhardt
- University of Stuttgart; Institute of Thermodynamics and Thermal Process Engineering; Pfaffenwaldring 9 70569 Stuttgart Germany
| | - Hamzeh Kraus
- University of Stuttgart; Institute of Thermodynamics and Thermal Process Engineering; Pfaffenwaldring 9 70569 Stuttgart Germany
| | - Daniel Markthaler
- University of Stuttgart; Institute of Thermodynamics and Thermal Process Engineering; Pfaffenwaldring 9 70569 Stuttgart Germany
| | - Niels Hansen
- University of Stuttgart; Institute of Thermodynamics and Thermal Process Engineering; Pfaffenwaldring 9 70569 Stuttgart Germany
| |
Collapse
|
18
|
Setiawan D, Brender J, Zhang Y. Recent advances in automated protein design and its future challenges. Expert Opin Drug Discov 2018; 13:587-604. [PMID: 29695210 DOI: 10.1080/17460441.2018.1465922] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/07/2023]
Abstract
INTRODUCTION Protein function is determined by protein structure which is in turn determined by the corresponding protein sequence. If the rules that cause a protein to adopt a particular structure are understood, it should be possible to refine or even redefine the function of a protein by working backwards from the desired structure to the sequence. Automated protein design attempts to calculate the effects of mutations computationally with the goal of more radical or complex transformations than are accessible by experimental techniques. Areas covered: The authors give a brief overview of the recent methodological advances in computer-aided protein design, showing how methodological choices affect final design and how automated protein design can be used to address problems considered beyond traditional protein engineering, including the creation of novel protein scaffolds for drug development. Also, the authors address specifically the future challenges in the development of automated protein design. Expert opinion: Automated protein design holds potential as a protein engineering technique, particularly in cases where screening by combinatorial mutagenesis is problematic. Considering solubility and immunogenicity issues, automated protein design is initially more likely to make an impact as a research tool for exploring basic biology in drug discovery than in the design of protein biologics.
Collapse
Affiliation(s)
- Dani Setiawan
- a Department of Computational Medicine and Bioinformatics , University of Michigan , Ann Arbor , MI , USA
| | - Jeffrey Brender
- b Radiation Biology Branch , Center for Cancer Research, National Cancer Institute - NIH , Bethesda , MD , USA
| | - Yang Zhang
- a Department of Computational Medicine and Bioinformatics , University of Michigan , Ann Arbor , MI , USA.,c Department of Biological Chemistry , University of Michigan , Ann Arbor , MI , USA
| |
Collapse
|
19
|
Deplazes E. Molecular simulations of venom peptide-membrane interactions: Progress and challenges. Pept Sci (Hoboken) 2018. [DOI: 10.1002/pep2.24060] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/21/2022]
Affiliation(s)
- Evelyne Deplazes
- School of Pharmacy and Biomedical Sciences; Curtin Health Innovation Research Institute, Curtin Institute for Computation, Curtin University; Bentley, Perth WA 6102 Australia
| |
Collapse
|
20
|
Vargiu AV, Ramaswamy VK, Malvacio I, Malloci G, Kleinekathöfer U, Ruggerone P. Water-mediated interactions enable smooth substrate transport in a bacterial efflux pump. Biochim Biophys Acta Gen Subj 2018; 1862:836-845. [PMID: 29339082 DOI: 10.1016/j.bbagen.2018.01.010] [Citation(s) in RCA: 36] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/21/2017] [Revised: 11/28/2017] [Accepted: 01/11/2018] [Indexed: 01/12/2023]
Abstract
BACKGROUND Efflux pumps of the Resistance-Nodulation-cell Division superfamily confer multi-drug resistance to Gram-negative bacteria. The most-studied polyspecific transporter belonging to this class is the inner-membrane trimeric antiporter AcrB of Escherichia coli. In previous studies, a functional rotation mechanism was proposed for its functioning, according to which the three monomers undergo concerted conformational changes facilitating the extrusion of substrates. However, the molecular determinants and the energetics of this mechanism still remain unknown, so its feasibility must be proven mechanistically. METHODS A computational protocol able to mimic the functional rotation mechanism in AcrB was developed. By using multi-bias molecular dynamics simulations we characterized the translocation of the substrate doxorubicin driven by conformational changes of the protein. In addition, we estimated for the first time the free energy profile associated to this process. RESULTS We provided a molecular view of the process in agreement with experimental data. Moreover, we showed that the conformational changes occurring in AcrB enable the formation of a layer of structured waters on the internal surface of the transport channel. This water layer, in turn, allows for a fairly constant hydration of the substrate, facilitating its diffusion over a smooth free energy profile. CONCLUSIONS Our findings reveal a new molecular mechanism of polyspecific transport whereby water contributes by screening potentially strong substrate-protein interactions. GENERAL SIGNIFICANCE We provided a mechanistic understanding of a fundamental process related to multi-drug transport. Our results can help rationalizing the behavior of other polyspecific transporters and designing compounds avoiding extrusion or inhibitors of efflux pumps.
Collapse
Affiliation(s)
- Attilio Vittorio Vargiu
- Department of Physics, University of Cagliari, s.p. 8, Cittadella Universitaria, 09042 Monserrato (CA), Italy.
| | - Venkata Krishnan Ramaswamy
- Department of Physics, University of Cagliari, s.p. 8, Cittadella Universitaria, 09042 Monserrato (CA), Italy
| | - Ivana Malvacio
- Department of Physics, University of Cagliari, s.p. 8, Cittadella Universitaria, 09042 Monserrato (CA), Italy
| | - Giuliano Malloci
- Department of Physics, University of Cagliari, s.p. 8, Cittadella Universitaria, 09042 Monserrato (CA), Italy
| | - Ulrich Kleinekathöfer
- Department of Physics & Earth Sciences, Jacobs University Bremen, Campus Ring 1, 28759 Bremen, Germany
| | - Paolo Ruggerone
- Department of Physics, University of Cagliari, s.p. 8, Cittadella Universitaria, 09042 Monserrato (CA), Italy
| |
Collapse
|
21
|
Cacciotto P, Ramaswamy VK, Malloci G, Ruggerone P, Vargiu AV. Molecular Modeling of Multidrug Properties of Resistance Nodulation Division (RND) Transporters. Methods Mol Biol 2018; 1700:179-219. [PMID: 29177832 DOI: 10.1007/978-1-4939-7454-2_11] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/07/2023]
Abstract
Efflux pumps of the resistance nodulation division (RND) superfamily are among the major contributors to intrinsic and acquired multidrug resistance in Gram-negative bacteria. Structural information on AcrAB-TolC and MexAB-OprM, major efflux pumps of Escherichia coli and Pseudomonas aeruginosa respectively, boosted intensive research aimed at understanding the molecular mechanisms ruling the active extrusion processes. In particular, several studies were devoted to the understanding of the determinants behind the extraordinary broad specificity of the RND transporters AcrB and MexB. In this chapter, we discuss the ever-growing role computational methods have been playing in deciphering key structural and dynamical features of these transporters and of their interaction with substrates and inhibitors. We further discuss and illustrate examples from our lab of how molecular docking, homology modeling, all-atom molecular dynamics simulations and in silico free energy estimations can all together give precious insights into the processes of recognition and extrusion of substrates, as well as on the possible inhibition strategies.
Collapse
Affiliation(s)
- Pierpaolo Cacciotto
- Department of Physics, University of Cagliari, s.p. 8, 09042, Monserrato, CA, Italy
| | - Venkata K Ramaswamy
- Department of Physics, University of Cagliari, s.p. 8, 09042, Monserrato, CA, Italy
| | - Giuliano Malloci
- Department of Physics, University of Cagliari, s.p. 8, 09042, Monserrato, CA, Italy
| | - Paolo Ruggerone
- Department of Physics, University of Cagliari, s.p. 8, 09042, Monserrato, CA, Italy
| | - Attilio V Vargiu
- Department of Physics, University of Cagliari, s.p. 8, 09042, Monserrato, CA, Italy.
| |
Collapse
|
22
|
Garai A, Mogurampelly S, Bag S, Maiti PK. Overstretching of B-DNA with various pulling protocols: Appearance of structural polymorphism and S-DNA. J Chem Phys 2017; 147:225102. [DOI: 10.1063/1.4991862] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/26/2022] Open
Affiliation(s)
- Ashok Garai
- Centre for Condensed Matter Theory, Department of Physics, Indian Institute of Science, Bangalore 560012, India
- Department of Physics, The LNM Institute of Information Technology, Jamdoli, Jaipur 302031, India
| | - Santosh Mogurampelly
- Centre for Condensed Matter Theory, Department of Physics, Indian Institute of Science, Bangalore 560012, India
| | - Saientan Bag
- Centre for Condensed Matter Theory, Department of Physics, Indian Institute of Science, Bangalore 560012, India
| | - Prabal K. Maiti
- Centre for Condensed Matter Theory, Department of Physics, Indian Institute of Science, Bangalore 560012, India
| |
Collapse
|
23
|
Pabon-Martinez YV, Xu Y, Villa A, Lundin KE, Geny S, Nguyen CH, Pedersen EB, Jørgensen PT, Wengel J, Nilsson L, Smith CIE, Zain R. LNA effects on DNA binding and conformation: from single strand to duplex and triplex structures. Sci Rep 2017; 7:11043. [PMID: 28887512 PMCID: PMC5591256 DOI: 10.1038/s41598-017-09147-8] [Citation(s) in RCA: 27] [Impact Index Per Article: 3.9] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/08/2017] [Accepted: 07/20/2017] [Indexed: 12/19/2022] Open
Abstract
The anti-gene strategy is based on sequence-specific recognition of double-strand DNA by triplex forming (TFOs) or DNA strand invading oligonucleotides to modulate gene expression. To be efficient, the oligonucleotides (ONs) should target DNA selectively, with high affinity. Here we combined hybridization analysis and electrophoretic mobility shift assay with molecular dynamics (MD) simulations to better understand the underlying structural features of modified ONs in stabilizing duplex- and triplex structures. Particularly, we investigated the role played by the position and number of locked nucleic acid (LNA) substitutions in the ON when targeting a c-MYC or FXN (Frataxin) sequence. We found that LNA-containing single strand TFOs are conformationally pre-organized for major groove binding. Reduced content of LNA at consecutive positions at the 3'-end of a TFO destabilizes the triplex structure, whereas the presence of Twisted Intercalating Nucleic Acid (TINA) at the 3'-end of the TFO increases the rate and extent of triplex formation. A triplex-specific intercalating benzoquinoquinoxaline (BQQ) compound highly stabilizes LNA-containing triplex structures. Moreover, LNA-substitution in the duplex pyrimidine strand alters the double helix structure, affecting x-displacement, slide and twist favoring triplex formation through enhanced TFO major groove accommodation. Collectively, these findings should facilitate the design of potent anti-gene ONs.
Collapse
Affiliation(s)
- Y Vladimir Pabon-Martinez
- Department of Laboratory Medicine, Clinical Research Center, Karolinska Institutet, SE-141 86, Huddinge, Stockholm, Sweden
| | - You Xu
- Department of Biosciences and Nutrition, Karolinska Institutet, SE-141 83, Huddinge, Sweden
| | - Alessandra Villa
- Department of Biosciences and Nutrition, Karolinska Institutet, SE-141 83, Huddinge, Sweden
| | - Karin E Lundin
- Department of Laboratory Medicine, Clinical Research Center, Karolinska Institutet, SE-141 86, Huddinge, Stockholm, Sweden
| | - Sylvain Geny
- Department of Laboratory Medicine, Clinical Research Center, Karolinska Institutet, SE-141 86, Huddinge, Stockholm, Sweden
| | - Chi-Hung Nguyen
- Institut Curie, PSL Research University, UMR 9187-U 1196, CNRS-Institut Curie, INSERM, Centre Universitaire, Orsay, France
| | - Erik B Pedersen
- Department of Physics, Chemistry and Pharmacy, Nucleic Acid Center, University of Southern Denmark, DK-5230, Odense M, Denmark
| | - Per T Jørgensen
- Department of Physics, Chemistry and Pharmacy, Nucleic Acid Center, University of Southern Denmark, DK-5230, Odense M, Denmark
| | - Jesper Wengel
- Department of Physics, Chemistry and Pharmacy, Nucleic Acid Center, University of Southern Denmark, DK-5230, Odense M, Denmark
| | - Lennart Nilsson
- Department of Biosciences and Nutrition, Karolinska Institutet, SE-141 83, Huddinge, Sweden
| | - C I Edvard Smith
- Department of Laboratory Medicine, Clinical Research Center, Karolinska Institutet, SE-141 86, Huddinge, Stockholm, Sweden
| | - Rula Zain
- Department of Laboratory Medicine, Clinical Research Center, Karolinska Institutet, SE-141 86, Huddinge, Stockholm, Sweden.
- Department of Clinical Genetics, Centre for Rare Diseases, Karolinska University Hospital, SE-171 76, Stockholm, Sweden.
| |
Collapse
|
24
|
Demharter S, Knapp B, Deane CM, Minary P. Modeling Functional Motions of Biological Systems by Customized Natural Moves. Biophys J 2017; 111:710-721. [PMID: 27558715 PMCID: PMC5002067 DOI: 10.1016/j.bpj.2016.06.028] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/01/2016] [Revised: 06/20/2016] [Accepted: 06/22/2016] [Indexed: 11/30/2022] Open
Abstract
Simulating the functional motions of biomolecular systems requires large computational resources. We introduce a computationally inexpensive protocol for the systematic testing of hypotheses regarding the dynamic behavior of proteins and nucleic acids. The protocol is based on natural move Monte Carlo, a highly efficient conformational sampling method with built-in customization capabilities that allows researchers to design and perform a large number of simulations to investigate functional motions in biological systems. We demonstrate the use of this protocol on both a protein and a DNA case study. Firstly, we investigate the plasticity of a class II major histocompatibility complex in the absence of a bound peptide. Secondly, we study the effects of the epigenetic mark 5-hydroxymethyl on cytosine on the structure of the Dickerson-Drew dodecamer. We show how our customized natural moves protocol can be used to investigate causal relationships of functional motions in biological systems.
Collapse
Affiliation(s)
- Samuel Demharter
- Department of Computer Science, University of Oxford, Oxford, UK
| | - Bernhard Knapp
- Department of Statistics, University of Oxford, Oxford, UK
| | | | - Peter Minary
- Department of Computer Science, University of Oxford, Oxford, UK.
| |
Collapse
|
25
|
Pasi M, Angelov D, Bednar J, Dimitrov S, Lavery R. Extra views on structure and dynamics of DNA loops on nucleosomes studied with molecular simulations. Nucleus 2016; 7:554-559. [PMID: 27874316 PMCID: PMC5214536 DOI: 10.1080/19491034.2016.1260800] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/10/2016] [Accepted: 11/10/2016] [Indexed: 02/07/2023] Open
Abstract
It has been shown experimentally that the action of the RSC chromatin remodeler leads to the formation of an irregular, partially remodeled nucleosome, termed a remosome. The remosome contains an extra 30-40 base pairs of DNA compared to a canonical nucleosome. Large-scale molecular simulations have provided information on the probable structure of remosomes and have explained why they remain stable in the absence of RSC. Here we explain how these simulations were carried out and what the resulting remosome models imply in terms of the mechanism of action of RSC. We notably show that local kinks within DNA are key in explaining how extra DNA can be in added to nucleosomes without unduly disturbing DNA-histone binding.
Collapse
Affiliation(s)
- Marco Pasi
- Université de Lyon/CNRS UMR 5086, MMSB, Institut de Biologie et Chimie des Protéines, Lyon, France
| | - Dimitar Angelov
- Université de Lyon/CNRS/INRA UMR 5239, LBMC, Ecole Normale Supérieure de Lyon, Lyon, France
| | - Jan Bednar
- Institut Albert Bonniot, Université de Grenoble Alpes/INSERM U1209/CNRS UMR 5309, Grenoble, France
| | - Stefan Dimitrov
- Institut Albert Bonniot, Université de Grenoble Alpes/INSERM U1209/CNRS UMR 5309, Grenoble, France
| | - Richard Lavery
- Université de Lyon/CNRS UMR 5086, MMSB, Institut de Biologie et Chimie des Protéines, Lyon, France
| |
Collapse
|
26
|
Wang Y, Papaleo E, Lindorff-Larsen K. Mapping transiently formed and sparsely populated conformations on a complex energy landscape. eLife 2016; 5. [PMID: 27552057 PMCID: PMC5050026 DOI: 10.7554/elife.17505] [Citation(s) in RCA: 46] [Impact Index Per Article: 5.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/04/2016] [Accepted: 08/22/2016] [Indexed: 12/11/2022] Open
Abstract
Determining the structures, kinetics, thermodynamics and mechanisms that underlie conformational exchange processes in proteins remains extremely difficult. Only in favourable cases is it possible to provide atomic-level descriptions of sparsely populated and transiently formed alternative conformations. Here we benchmark the ability of enhanced-sampling molecular dynamics simulations to determine the free energy landscape of the L99A cavity mutant of T4 lysozyme. We find that the simulations capture key properties previously measured by NMR relaxation dispersion methods including the structure of a minor conformation, the kinetics and thermodynamics of conformational exchange, and the effect of mutations. We discover a new tunnel that involves the transient exposure towards the solvent of an internal cavity, and show it to be relevant for ligand escape. Together, our results provide a comprehensive view of the structural landscape of a protein, and point forward to studies of conformational exchange in systems that are less characterized experimentally. DOI:http://dx.doi.org/10.7554/eLife.17505.001 Proteins are the workhorses of cells, where they perform a wide range of roles. To do so, they adopt specific three-dimensional structures that enable them to interact with other molecules as necessary. Often a protein needs to be able to shift between different states with distinct structures as it goes about its job. To fully understand how a protein works, it is important to be able to characterize these different structures and how the protein changes between them. Many of the experimental techniques used to study protein structure rely on isolating the individual structural forms of a protein. Since many structures only exist briefly, this can be very difficult. To complement experimental results, computer simulations allow researchers to model how atoms behave within a molecule. However, a number of factors limit how well these models represent what happens experimentally, such as the accuracy of the physical description used for the modeling. Wang et al. set out to test and benchmark how well computer simulations could model changes in structure for a protein called T4 lysozyme, which has been studied extensively using experimental techniques. T4 lysozyme exists in two different states that have distinct structures. By comparing existing detailed experimental measurements with the results of their simulations, Wang et al. found that the simulations could capture key aspects of how T4 lysozyme changes its shape. The simulations described the structure of the protein in both states and accurately determined the relative proportion of molecules that are found in each state. They could also determine how long it takes for a molecule to change its shape from one state to the other. The findings allowed Wang et al. to describe in fine detail – down to the level of individual atoms – how the protein changes its shape and how mutations in the protein affect its ability to do so. A key question for future studies is whether these insights can be extended to other proteins that are less well characterized experimentally than T4 lysozyme. DOI:http://dx.doi.org/10.7554/eLife.17505.002
Collapse
Affiliation(s)
- Yong Wang
- Structural Biology and NMR Laboratory, Linderstrøm-Lang Centre for Protein Science, Department of Biology, University of Copenhagen, Copenhagen, Denmark
| | - Elena Papaleo
- Structural Biology and NMR Laboratory, Linderstrøm-Lang Centre for Protein Science, Department of Biology, University of Copenhagen, Copenhagen, Denmark
| | - Kresten Lindorff-Larsen
- Structural Biology and NMR Laboratory, Linderstrøm-Lang Centre for Protein Science, Department of Biology, University of Copenhagen, Copenhagen, Denmark
| |
Collapse
|
27
|
Santiago G, de Salas F, Lucas MF, Monza E, Acebes S, Martinez ÁT, Camarero S, Guallar V. Computer-Aided Laccase Engineering: Toward Biological Oxidation of Arylamines. ACS Catal 2016. [DOI: 10.1021/acscatal.6b01460] [Citation(s) in RCA: 44] [Impact Index Per Article: 5.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/06/2023]
Affiliation(s)
- Gerard Santiago
- Joint
BSC-CRG-IRB Research Program in Computational Biology, Barcelona Supercomputing Center, Jordi Girona 29, E-08034 Barcelona, Spain
| | - Felipe de Salas
- Centro de Investigaciones Biológicas, CSIC, Ramiro de Maeztu 9, E-28040 Madrid, Spain
| | - M. Fátima Lucas
- Joint
BSC-CRG-IRB Research Program in Computational Biology, Barcelona Supercomputing Center, Jordi Girona 29, E-08034 Barcelona, Spain
- Anaxomics Biotech, Balmes 89, E-08008 Barcelona, Spain
| | - Emanuele Monza
- Joint
BSC-CRG-IRB Research Program in Computational Biology, Barcelona Supercomputing Center, Jordi Girona 29, E-08034 Barcelona, Spain
| | - Sandra Acebes
- Joint
BSC-CRG-IRB Research Program in Computational Biology, Barcelona Supercomputing Center, Jordi Girona 29, E-08034 Barcelona, Spain
| | - Ángel T. Martinez
- Centro de Investigaciones Biológicas, CSIC, Ramiro de Maeztu 9, E-28040 Madrid, Spain
| | - Susana Camarero
- Centro de Investigaciones Biológicas, CSIC, Ramiro de Maeztu 9, E-28040 Madrid, Spain
| | - Víctor Guallar
- Joint
BSC-CRG-IRB Research Program in Computational Biology, Barcelona Supercomputing Center, Jordi Girona 29, E-08034 Barcelona, Spain
- ICREA, Passeig Lluís Companys 23, E-08010 Barcelona, Spain
| |
Collapse
|
28
|
Structural effects of modified ribonucleotides and magnesium in transfer RNAs. Bioorg Med Chem 2016; 24:4826-4834. [PMID: 27364608 DOI: 10.1016/j.bmc.2016.06.037] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Key Words] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/26/2016] [Revised: 06/16/2016] [Accepted: 06/17/2016] [Indexed: 11/20/2022]
Abstract
Modified nucleotides are ubiquitous and important to tRNA structure and function. To understand their effect on tRNA conformation, we performed a series of molecular dynamics simulations on yeast tRNAPhe and tRNAinit, Escherichia coli tRNAinit and HIV tRNALys. Simulations were performed with the wild type modified nucleotides, using the recently developed CHARMM compatible force field parameter set for modified nucleotides (J. Comput. Chem.2016, 37, 896), or with the corresponding unmodified nucleotides, and in the presence or absence of Mg2+. Results showed a stabilizing effect associated with the presence of the modifications and Mg2+ for some important positions, such as modified guanosine in position 37 and dihydrouridines in 16/17 including both structural properties and base interactions. Some other modifications were also found to make subtle contributions to the structural properties of local domains. While we were not able to investigate the effect of adenosine 37 in tRNAinit and limitations were observed in the conformation of E. coli tRNAinit, the presence of the modified nucleotides and of Mg2+ better maintained the structural features and base interactions of the tRNA systems than in their absence indicating the utility of incorporating the modified nucleotides in simulations of tRNA and other RNAs.
Collapse
|
29
|
Xu Y, Vanommeslaeghe K, Aleksandrov A, MacKerell AD, Nilsson L. Additive CHARMM force field for naturally occurring modified ribonucleotides. J Comput Chem 2016; 37:896-912. [PMID: 26841080 PMCID: PMC4801715 DOI: 10.1002/jcc.24307] [Citation(s) in RCA: 50] [Impact Index Per Article: 6.3] [Reference Citation Analysis] [Abstract] [Key Words] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/29/2015] [Revised: 11/16/2015] [Accepted: 01/06/2016] [Indexed: 01/13/2023]
Abstract
More than 100 naturally occurring modified nucleotides have been found in RNA molecules, in particular in tRNAs. We have determined molecular mechanics force field parameters compatible with the CHARMM36 all‐atom additive force field for all these modifications using the CHARMM force field parametrization strategy. Emphasis was placed on fine tuning of the partial atomic charges and torsion angle parameters. Quantum mechanics calculations on model compounds provided the initial set of target data, and extensive molecular dynamics simulations of nucleotides and oligonucleotides in aqueous solutions were used for further refinement against experimental data. The presented parameters will allow for computational studies of a wide range of RNAs containing modified nucleotides, including the ribosome and transfer RNAs. © 2016 The Authors. Journal of Computational Chemistry Published by Wiley Periodicals, Inc.
Collapse
Affiliation(s)
- You Xu
- Department of Biosciences and Nutrition, Karolinska Institutet, HUDDINGE, SE-141 83, Sweden
| | - Kenno Vanommeslaeghe
- Department of Pharmaceutical Sciences, School of Pharmacy, University of Maryland, 20 Penn Street, Baltimore, Maryland, 21201.,Department of Analytical Chemistry and Pharmaceutical Technology (FABI), Center for Pharmaceutical Research (CePhaR), Vrije Universiteit Brussel (VUB), Laarbeeklaan 103, Brussels, B-1090, Belgium
| | - Alexey Aleksandrov
- Department of Biology, Ecole Polytechnique, Laboratoire De Biochimie (CNRS UMR7654), Palaiseau, F-91128, France
| | - Alexander D MacKerell
- Department of Pharmaceutical Sciences, School of Pharmacy, University of Maryland, 20 Penn Street, Baltimore, Maryland, 21201
| | - Lennart Nilsson
- Department of Biosciences and Nutrition, Karolinska Institutet, HUDDINGE, SE-141 83, Sweden
| |
Collapse
|
30
|
Extracting Conformational Ensembles of Small Molecules from Molecular Dynamics Simulations: Ampicillin as a Test Case. COMPUTATION 2016. [DOI: 10.3390/computation4010005] [Citation(s) in RCA: 15] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/16/2022]
|
31
|
Machado MR, Pantano S. SIRAH tools: mapping, backmapping and visualization of coarse-grained models. Bioinformatics 2016; 32:1568-70. [DOI: 10.1093/bioinformatics/btw020] [Citation(s) in RCA: 69] [Impact Index Per Article: 8.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/09/2015] [Accepted: 01/08/2016] [Indexed: 11/13/2022] Open
|
32
|
Ley K, Christofferson A, Penna M, Winkler D, Maclaughlin S, Yarovsky I. Surface-water Interface Induces Conformational Changes Critical for Protein Adsorption: Implications for Monolayer Formation of EAS Hydrophobin. Front Mol Biosci 2015; 2:64. [PMID: 26636091 PMCID: PMC4644811 DOI: 10.3389/fmolb.2015.00064] [Citation(s) in RCA: 24] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/09/2015] [Accepted: 10/29/2015] [Indexed: 11/13/2022] Open
Abstract
The class I hydrophobin EAS is part of a family of small, amphiphilic fungal proteins best known for their ability to self-assemble into stable monolayers that modify the hydrophobicity of a surface to facilitate further microbial growth. These proteins have attracted increasing attention for industrial and biomedical applications, with the aim of designing surfaces that have the potential to maintain their clean state by resisting non-specific protein binding. To gain a better understanding of this process, we have employed all-atom molecular dynamics to study initial stages of the spontaneous adsorption of monomeric EAS hydrophobin on fully hydroxylated silica, a commonly used industrial and biomedical substrate. Particular interest has been paid to the Cys3-Cys4 loop, which has been shown to exhibit disruptive behavior in solution, and the Cys7-Cys8 loop, which is believed to be involved in the aggregation of EAS hydrophobin at interfaces. Specific and water mediated interactions with the surface were also analyzed. We have identified two possible binding motifs, one which allows unfolding of the Cys7-Cys8 loop due to the surfactant-like behavior of the Cys3-Cys4 loop, and another which has limited unfolding due to the Cys3-Cys4 loop remaining disordered in solution. We have also identified intermittent interactions with water which mediate the protein adsorption to the surface, as well as longer lasting interactions which control the diffusion of water around the adsorption site. These results have shown that EAS behaves in a similar way at the air-water and surface-water interfaces, and have also highlighted the need for hydrophilic ligand functionalization of the silica surface in order to prevent the adsorption of EAS hydrophobin.
Collapse
Affiliation(s)
- Kamron Ley
- Health Innovations Research Institute and School of Aerospace, Mechanical and Manufacturing Engineering, RMIT University Melbourne, VIC, Australia
| | - Andrew Christofferson
- Health Innovations Research Institute and School of Aerospace, Mechanical and Manufacturing Engineering, RMIT University Melbourne, VIC, Australia
| | - Matthew Penna
- Health Innovations Research Institute and School of Aerospace, Mechanical and Manufacturing Engineering, RMIT University Melbourne, VIC, Australia ; Australian Research Council (ARC) Research Hub for Australian Steel Manufacturing Wollongong, NSW, Australia
| | - Dave Winkler
- CSIRO, Manufacturing Flagship Clayton, VIC, Australia ; Institute of Pharmaceutical Science, Monash University Parkville, VIC, Australia ; Institute for Molecular Science, Latrobe University Bundoora, VIC, Australia
| | - Shane Maclaughlin
- Australian Research Council (ARC) Research Hub for Australian Steel Manufacturing Wollongong, NSW, Australia ; BlueScope Steel Research Laboratories Port Kembla, NSW, Australia
| | - Irene Yarovsky
- Health Innovations Research Institute and School of Aerospace, Mechanical and Manufacturing Engineering, RMIT University Melbourne, VIC, Australia ; Australian Research Council (ARC) Research Hub for Australian Steel Manufacturing Wollongong, NSW, Australia
| |
Collapse
|
33
|
Malloci G, Vargiu AV, Serra G, Bosin A, Ruggerone P, Ceccarelli M. A Database of Force-Field Parameters, Dynamics, and Properties of Antimicrobial Compounds. Molecules 2015; 20:13997-4021. [PMID: 26247924 PMCID: PMC6332394 DOI: 10.3390/molecules200813997] [Citation(s) in RCA: 44] [Impact Index Per Article: 4.9] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/11/2015] [Accepted: 07/28/2015] [Indexed: 02/01/2023] Open
Abstract
We present an on-line database of all-atom force-field parameters and molecular properties of compounds with antimicrobial activity (mostly antibiotics and some beta-lactamase inhibitors). For each compound, we provide the General Amber Force Field parameters for the major species at physiological pH, together with an analysis of properties of interest as extracted from µs-long molecular dynamics simulations in explicit water solution. The properties include number and population of structural clusters, molecular flexibility, hydrophobic and hydrophilic molecular surfaces, the statistics of intraand inter-molecular H-bonds, as well as structural and dynamical properties of solvent molecules within first and second solvation shells. In addition, the database contains several key molecular parameters, such as energy of the frontier molecular orbitals, vibrational properties, rotational constants, atomic partial charges and electric dipole moment, computed by Density Functional Theory. The present database (to our knowledge the first extensive one including dynamical properties) is part of a wider project aiming to build-up a database containing structural, physico-chemical and dynamical properties of medicinal compounds using different force-field parameters with increasing level of complexity and reliability. The database is freely accessible at http://www.dsf.unica.it/translocation/db/.
Collapse
Affiliation(s)
- Giuliano Malloci
- Dipartimento di Fisica, Università degli studi di Cagliari, Cittadella Universitaria, I-09042 Monserrato (Cagliari), Italy.
| | - Attilio Vittorio Vargiu
- Dipartimento di Fisica, Università degli studi di Cagliari, Cittadella Universitaria, I-09042 Monserrato (Cagliari), Italy.
| | - Giovanni Serra
- Dipartimento di Fisica, Università degli studi di Cagliari, Cittadella Universitaria, I-09042 Monserrato (Cagliari), Italy.
| | - Andrea Bosin
- Dipartimento di Fisica, Università degli studi di Cagliari, Cittadella Universitaria, I-09042 Monserrato (Cagliari), Italy.
| | - Paolo Ruggerone
- Dipartimento di Fisica, Università degli studi di Cagliari, Cittadella Universitaria, I-09042 Monserrato (Cagliari), Italy.
| | - Matteo Ceccarelli
- Dipartimento di Fisica, Università degli studi di Cagliari, Cittadella Universitaria, I-09042 Monserrato (Cagliari), Italy.
| |
Collapse
|