1
|
Pan X, Van R, Epifanovsky E, Liu J, Pu J, Nam K, Shao Y. Accelerating Ab Initio Quantum Mechanical and Molecular Mechanical (QM/MM) Molecular Dynamics Simulations with Multiple Time Step Integration and a Recalibrated Semiempirical QM/MM Hamiltonian. J Phys Chem B 2022; 126:10.1021/acs.jpcb.2c02262. [PMID: 35653199 PMCID: PMC9715852 DOI: 10.1021/acs.jpcb.2c02262] [Citation(s) in RCA: 4] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Molecular dynamics (MD) simulations employing ab initio quantum mechanical and molecular mechanical (ai-QM/MM) potentials are considered to be the state of the art, but the high computational cost associated with the ai-QM calculations remains a theoretical challenge for their routine application. Here, we present a modified protocol of the multiple time step (MTS) method for accelerating ai-QM/MM MD simulations of condensed-phase reactions. Within a previous MTS protocol [Nam J. Chem. Theory Comput. 2014, 10, 4175], reference forces are evaluated using a low-level (semiempirical QM/MM) Hamiltonian and employed at inner time steps to propagate the nuclear motions. Correction forces, which arise from the force differences between high-level (ai-QM/MM) and low-level Hamiltonians, are applied at outer time steps, where the MTS algorithm allows the time-reversible integration of the correction forces. To increase the outer step size, which is bound by the highest-frequency component in the correction forces, the semiempirical QM Hamiltonian is recalibrated in this work to minimize the magnitude of the correction forces. The remaining high-frequency modes, which are mainly bond stretches involving hydrogen atoms, are then removed from the correction forces. When combined with a Langevin or SIN(R) thermostat, the modified MTS-QM/MM scheme remains robust with an up to 8 (with Langevin) or 10 fs (with SIN(R)) outer time step (with 1 fs inner time steps) for the chorismate mutase system. This leads to an over 5-fold speedup over standard ai-QM/MM simulations, without sacrificing the accuracy in the predicted free energy profile of the reaction.
Collapse
Affiliation(s)
- Xiaoliang Pan
- Department of Chemistry and Biochemistry, University of Oklahoma, 101 Stephenson Parkway, Norman, Oklahoma 73019-5251, United States
| | - Richard Van
- Department of Chemistry and Biochemistry, University of Oklahoma, 101 Stephenson Parkway, Norman, Oklahoma 73019-5251, United States
| | - Evgeny Epifanovsky
- Q-Chem, Inc., 6601 Owens Drive, Suite 105, Pleasanton, California 94588, United States
| | - Jian Liu
- Beijing National Laboratory for Molecular Sciences, Institute of Theoretical and Computational Chemistry, College of Chemistry and Molecular Engineering, Peking University, Beijing 100871, China
| | - Jingzhi Pu
- Department of Chemistry and Chemical Biology, Indiana University-Purdue University Indianapolis, 402 N Blackford St., LD326, Indianapolis, Indiana 46202, United States
| | - Kwangho Nam
- Department of Chemistry and Biochemistry, University of Texas at Arlington, Arlington, Texas 76019, United States
| | - Yihan Shao
- Department of Chemistry and Biochemistry, University of Oklahoma, 101 Stephenson Parkway, Norman, Oklahoma 73019-5251, United States
| |
Collapse
|
2
|
Ngo VA, Jha RK. Identifying key determinants and dynamics of SARS-CoV-2/ACE2 tight interaction. PLoS One 2021; 16:e0257905. [PMID: 34582502 PMCID: PMC8478251 DOI: 10.1371/journal.pone.0257905] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/27/2021] [Accepted: 09/13/2021] [Indexed: 12/15/2022] Open
Abstract
SARS-CoV-2 virus, the causative agent of Covid-19, has fired up a global pandemic. The virus interacts with the human receptor angiotensin-converting enzyme 2 (ACE2) for an invasion via receptor binding domain (RBD) on its spike protein. To provide a deeper understanding of this interaction, we performed microsecond simulations of the RBD-ACE2 complex for SARS-CoV-2 and compared it with the closely related SARS-CoV discovered in 2003. We show residues in the RBD of SARS-CoV-2 that were mutated from SARS-CoV, collectively help make the RBD anchor much stronger to the N-terminal part of ACE2 than the corresponding residues on RBD of SARS-CoV. This would result in a reduced dissociation rate of SARS-CoV-2 from human receptor protein compared to SARS-CoV. The phenomenon was consistently observed in simulations beyond 500 ns and was reproducible across different force fields. Altogether, our study adds more insight into the critical dynamics of the key residues at the virus spike and human receptor binding interface and potentially aids the development of diagnostics and therapeutics to combat the pandemic efficiently.
Collapse
Affiliation(s)
- Van A. Ngo
- Center for Nonlinear Studies, Los Alamos National Laboratory, Los Alamos, New Mexico, United States of America
| | - Ramesh K. Jha
- Bioscience Division, Los Alamos National Laboratory, Los Alamos, New Mexico, United States of America
| |
Collapse
|
3
|
In Silico Prediction of the Binding, Folding, Insertion, and Overall Stability of Membrane-Active Peptides. Methods Mol Biol 2021; 2315:161-182. [PMID: 34302676 DOI: 10.1007/978-1-0716-1468-6_10] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/18/2023]
Abstract
Membrane-active peptides (MAPs) are short-length peptides used for potential biomedical applications in diagnostic imaging of tissues, targeted drug delivery, gene delivery, and antimicrobials and antibiotics. The broad appeal of MAPs is that they are infinitely variable, relatively low cost, and biocompatible. However, experimentally characterizing the specific properties of a MAP or its many variants is a low-resolution and potentially time-consuming endeavor; molecular dynamics (MD) simulations have emerged as an invaluable tool in identifying the biophysical interactions that are fundamental to the function of MAPs. In this chapter, a step-by-step approach to discreetly model the binding, folding, and insertion of a membrane-active peptide to a model lipid bilayer using MD simulations is described. Detailed discussion is devoted to the critical aspects of running these types of simulations: prior knowledge of the system, understanding the strengths and weaknesses of molecular mechanics force fields, proper construction and equilibration of the system, realistically estimating both experimental and computational timescales, and leveraging analysis to make direct comparisons to experimental results as often as possible.
Collapse
|
4
|
Bureau HR, Quirk S, Hernandez R. The relative stability of trpzip1 and its mutants determined by computation and experiment. RSC Adv 2020; 10:6520-6535. [PMID: 35495997 PMCID: PMC9049704 DOI: 10.1039/d0ra00920b] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/13/2019] [Accepted: 02/04/2020] [Indexed: 11/21/2022] Open
Abstract
The single-point mutations of tprzip1 are indicated at left, and their relative energetics are compared at right.
Collapse
|
5
|
Nguyen TH, Ngo V, Castro Zerba JP, Noskov S, Minh DDL. Nonequilibrium path-ensemble averages for symmetric protocols. J Chem Phys 2019; 151:194103. [DOI: 10.1063/1.5121306] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Trung Hai Nguyen
- Laboratory of Theoretical and Computational Biophysics, Ton Duc Thang University, Ho Chi Minh City, Vietnam
- Faculty of Applied Sciences, Ton Duc Thang University, Ho Chi Minh City, Vietnam
| | - Van Ngo
- Center for Nonlinear Studies, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, USA
| | - João Paulo Castro Zerba
- Department of Chemistry, Illinois Institute of Technology, Chicago, Illinois 60616, USA
- Institute of Physics “Gleb Wataghin”, University of Campinas, Campinas, SP 13083-859, Brazil
| | - Sergei Noskov
- SYN: Centre for Molecular Simulation, Department of Biological Sciences, University of Calgary, Calgary, Alberta T2N 1N4, Canada
| | - David D. L. Minh
- Department of Chemistry, Illinois Institute of Technology, Chicago, Illinois 60616, USA
| |
Collapse
|
6
|
Célerse F, Lagardère L, Derat E, Piquemal JP. Massively Parallel Implementation of Steered Molecular Dynamics in Tinker-HP: Comparisons of Polarizable and Non-Polarizable Simulations of Realistic Systems. J Chem Theory Comput 2019; 15:3694-3709. [DOI: 10.1021/acs.jctc.9b00199] [Citation(s) in RCA: 16] [Impact Index Per Article: 3.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Frédéric Célerse
- Laboratoire de Chimie Théorique, UMR 7616 CNRS, Sorbonne Université, 75005 Paris, France
- Institut Parisien de Chimie Moléculaire, UMR 8232 CNRS, Sorbonne Université, 75005 Paris, France
| | - Louis Lagardère
- Institut des Sciences du Calcul et des Données, Sorbonne Université, 75005 Paris, France
- Institut Parisien de Chimie Physique et Théorique, FR 2622 CNRS, Sorbonne Université, 75005 Paris, France
- Laboratoire de Chimie théorique, UMR 7616 CNRS, Sorbonne Université, 75005 Paris, France
| | - Etienne Derat
- Institut Parisien de Chimie Moléculaire, UMR 8232 CNRS, Sorbonne Université, 75005 Paris, France
| | - Jean-Philip Piquemal
- Laboratoire de Chimie Théorique, UMR 7616 CNRS, Sorbonne Université, 75005 Paris, France
- Department of Biomedical Engineering, The University of Texas at Austin, Austin, Texas 78712, United States
- Institut Universitaire de France, 75005 Paris, France
| |
Collapse
|
7
|
Wang M, Mei Y, Ryde U. Predicting Relative Binding Affinity Using Nonequilibrium QM/MM Simulations. J Chem Theory Comput 2018; 14:6613-6622. [DOI: 10.1021/acs.jctc.8b00685] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/16/2022]
Affiliation(s)
- Meiting Wang
- State Key Laboratory of Precision Spectroscopy, School of Physics and Materials Science, East China Normal University, Shanghai 200062, China
- Department of Theoretical Chemistry, Lund University, Chemical Centre, P.O. Box 124, SE-221 00 Lund, Sweden
| | - Ye Mei
- State Key Laboratory of Precision Spectroscopy, School of Physics and Materials Science, East China Normal University, Shanghai 200062, China
- NYU−ECNU Center for Computational Chemistry at NYU Shanghai, Shanghai 200062, China
| | - Ulf Ryde
- Department of Theoretical Chemistry, Lund University, Chemical Centre, P.O. Box 124, SE-221 00 Lund, Sweden
| |
Collapse
|
8
|
Plata CA, Scholl ZN, Marszalek PE, Prados A. Relevance of the Speed and Direction of Pulling in Simple Modular Proteins. J Chem Theory Comput 2018; 14:2910-2918. [DOI: 10.1021/acs.jctc.8b00347] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Carlos A. Plata
- Física Teórica, Universidad de Sevilla, Apdo. de Correos 1065, Sevilla 41080, Spain
| | - Zackary N. Scholl
- Department of Physics, University of Alberta, Edmonton, Alberta T6G 2E1, Canada
| | - Piotr E. Marszalek
- Department of Mechanical Engineering and Materials Science, Duke University, Durham 27708, North Carolina, United States
| | - Antonio Prados
- Física Teórica, Universidad de Sevilla, Apdo. de Correos 1065, Sevilla 41080, Spain
| |
Collapse
|
9
|
Li PF, Jia XY, Wang MT, Mei Y. Comparison of Accuracy and Convergence Rate between Equilibrium and Nonequilibrium Alchemical Transformations for Calculation of Relative Binding Free Energy. CHINESE J CHEM PHYS 2017. [DOI: 10.1063/1674-0068/30/cjcp1711204] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/18/2023]
Affiliation(s)
- Peng-fei Li
- State Key Laboratory of Precision Spectroscopy, School of Physics and Materials Science, East China Normal University, Shanghai 200062, China
| | - Xiang-yu Jia
- State Key Laboratory of Precision Spectroscopy, School of Physics and Materials Science, East China Normal University, Shanghai 200062, China
- NYU-ECNU Center for Computational Chemistry at NYU Shanghai, Shanghai 200062, China
| | - Mei-ting Wang
- State Key Laboratory of Precision Spectroscopy, School of Physics and Materials Science, East China Normal University, Shanghai 200062, China
| | - Ye Mei
- State Key Laboratory of Precision Spectroscopy, School of Physics and Materials Science, East China Normal University, Shanghai 200062, China
- NYU-ECNU Center for Computational Chemistry at NYU Shanghai, Shanghai 200062, China
- Department of Chemistry and Biochemistry, University of Oklahoma, Norman OK 73019, USA
| |
Collapse
|
10
|
Miranda WE, Ngo VA, Perissinotti LL, Noskov SY. Computational membrane biophysics: From ion channel interactions with drugs to cellular function. BIOCHIMICA ET BIOPHYSICA ACTA. PROTEINS AND PROTEOMICS 2017; 1865:1643-1653. [PMID: 28847523 PMCID: PMC5764198 DOI: 10.1016/j.bbapap.2017.08.008] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 05/07/2017] [Revised: 08/16/2017] [Accepted: 08/16/2017] [Indexed: 12/16/2022]
Abstract
The rapid development of experimental and computational techniques has changed fundamentally our understanding of cellular-membrane transport. The advent of powerful computers and refined force-fields for proteins, ions, and lipids has expanded the applicability of Molecular Dynamics (MD) simulations. A myriad of cellular responses is modulated through the binding of endogenous and exogenous ligands (e.g. neurotransmitters and drugs, respectively) to ion channels. Deciphering the thermodynamics and kinetics of the ligand binding processes to these membrane proteins is at the heart of modern drug development. The ever-increasing computational power has already provided insightful data on the thermodynamics and kinetics of drug-target interactions, free energies of solvation, and partitioning into lipid bilayers for drugs. This review aims to provide a brief summary about modeling approaches to map out crucial binding pathways with intermediate conformations and free-energy surfaces for drug-ion channel binding mechanisms that are responsible for multiple effects on cellular functions. We will discuss post-processing analysis of simulation-generated data, which are then transformed to kinetic models to better understand the molecular underpinning of the experimental observables under the influence of drugs or mutations in ion channels. This review highlights crucial mathematical frameworks and perspectives on bridging different well-established computational techniques to connect the dynamics and timescales from all-atom MD and free energy simulations of ion channels to the physiology of action potentials in cellular models. This article is part of a Special Issue entitled: Biophysics in Canada, edited by Lewis Kay, John Baenziger, Albert Berghuis and Peter Tieleman.
Collapse
Affiliation(s)
- Williams E Miranda
- Centre for Molecular Simulations, Department of Biological Sciences, University of Calgary, Calgary, AB, Canada
| | - Van A Ngo
- Centre for Molecular Simulations, Department of Biological Sciences, University of Calgary, Calgary, AB, Canada
| | - Laura L Perissinotti
- Centre for Molecular Simulations, Department of Biological Sciences, University of Calgary, Calgary, AB, Canada
| | - Sergei Yu Noskov
- Centre for Molecular Simulations, Department of Biological Sciences, University of Calgary, Calgary, AB, Canada.
| |
Collapse
|
11
|
Wang M, Li P, Jia X, Liu W, Shao Y, Hu W, Zheng J, Brooks BR, Mei Y. Efficient Strategy for the Calculation of Solvation Free Energies in Water and Chloroform at the Quantum Mechanical/Molecular Mechanical Level. J Chem Inf Model 2017; 57:2476-2489. [DOI: 10.1021/acs.jcim.7b00001] [Citation(s) in RCA: 24] [Impact Index Per Article: 3.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Meiting Wang
- State Key Laboratory
of Precision Spectroscopy, School of Physics and Materials Science, East China Normal University, Shanghai 200062, China
| | - Pengfei Li
- State Key Laboratory
of Precision Spectroscopy, School of Physics and Materials Science, East China Normal University, Shanghai 200062, China
| | - Xiangyu Jia
- State Key Laboratory
of Precision Spectroscopy, School of Physics and Materials Science, East China Normal University, Shanghai 200062, China
| | - Wei Liu
- State Key Laboratory
of Precision Spectroscopy, School of Physics and Materials Science, East China Normal University, Shanghai 200062, China
| | - Yihan Shao
- Q-Chem Inc., 6601 Owens Drive, Suite 105, Pleasanton, California 94588, United States
- Department of Chemistry and Biochemistry, University of Oklahoma, Norman Oklahoma 73019, United States
| | - Wenxin Hu
- The Computer
Center, School of Computer Science and Software Engineering, East China Normal University, Shanghai 200062, China
| | - Jun Zheng
- The Computer
Center, School of Computer Science and Software Engineering, East China Normal University, Shanghai 200062, China
| | - Bernard R. Brooks
- Laboratory of Computational Biology, National Institutes of Health, National Heart, Lung and Blood Institute, 5635 Fishers Lane, T-900 Suite, Rockville, Maryland 20852, United States
| | - Ye Mei
- State Key Laboratory
of Precision Spectroscopy, School of Physics and Materials Science, East China Normal University, Shanghai 200062, China
- NYU-ECNU Center for Computational Chemistry at NYU Shanghai, Shanghai 200062, China
| |
Collapse
|
12
|
Mapping Ryanodine Binding Sites in the Pore Cavity of Ryanodine Receptors. Biophys J 2017; 112:1645-1653. [PMID: 28445755 DOI: 10.1016/j.bpj.2017.03.014] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/07/2016] [Revised: 03/10/2017] [Accepted: 03/21/2017] [Indexed: 02/03/2023] Open
Abstract
Ryanodine (Ryd) irreversibly targets ryanodine receptors (RyRs), a family of intracellular calcium release channels essential for many cellular processes ranging from muscle contraction to learning and memory. Little is known of the atomistic details about how Ryd binds to RyRs. In this study, we used all-atom molecular dynamics simulations with both enhanced and bidirectional sampling to gain direct insights into how Ryd interacts with major residues in RyRs that were experimentally determined to be critical for its binding. We found that the pyrrolic ring of Ryd displays preference for the R4892AGGG-F4921 residues in the cavity of RyR1, which explain the effects of the corresponding mutations in RyR2 in experiments. Particularly, the mutant Q4933A (or Q4863A in RyR2) critical for both the gating and Ryd binding not only has significantly less interaction with Ryd than the wild-type, but also yields more space for Ryd and water molecules in the cavity. These results describe clear binding modes of Ryd in the RyR cavity and offer structural mechanisms explaining functional data collected on RyR blockade.
Collapse
|
13
|
Zhekova HR, Ngo V, da Silva MC, Salahub D, Noskov S. Selective ion binding and transport by membrane proteins – A computational perspective. Coord Chem Rev 2017. [DOI: 10.1016/j.ccr.2017.03.019] [Citation(s) in RCA: 19] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/19/2023]
|
14
|
Edler E, Schulze E, Stein M. Membrane localization and dynamics of geranylgeranylated Rab5 hypervariable region. BIOCHIMICA ET BIOPHYSICA ACTA-BIOMEMBRANES 2017; 1859:1335-1349. [PMID: 28455099 DOI: 10.1016/j.bbamem.2017.04.021] [Citation(s) in RCA: 15] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 02/01/2017] [Revised: 04/18/2017] [Accepted: 04/23/2017] [Indexed: 12/11/2022]
Abstract
The small GTPase Rab5 is a key regulator of endosomal trafficking processes and a marker for the early endosome. The C-terminal hypervariable region (HVR) of Rab5 is post-translationally modified at residues Cys212 and Cys213 to accommodate two geranylgeranyl anchors (C20 carbon chain length) in order to associate Rab5 with the membrane. The structural role of the HVR regarding protein-early endosome membrane recruitment is not resolved due to its high degree of flexibility and lack of crystallographic information. Here, full-atomistic and coarse-grained molecular dynamics simulations of the truncated Rab5 HVR206-215 in three model membranes of increasing complexity (pure phospholipid bilayer, ternary membrane with cholesterol, six-component early endosome) were performed. Specific electrostatic interactions between the HVR206-215 Arg209 residue and the phosphate group of the inositol ring of PI(3)P were detected. This shows that PI(3)P acts as a first contact site of protein recruitment to the early endosome. The free energy change of HVR206-215 extraction from the bilayer was largest for the physiological negatively charged membrane. 5μs coarse-grained simulations revealed an active recruitment of PI(3)P to the HVR206-215 supporting the formation of Rab5- and PI(3)P enriched signaling platforms.
Collapse
Affiliation(s)
- Eileen Edler
- Molecular Simulations and Design Group, Max Planck Institute for Dynamics of Complex Technical Systems, Sandtorstrasse 1, 39106 Magdeburg, Germany
| | - Eric Schulze
- Molecular Simulations and Design Group, Max Planck Institute for Dynamics of Complex Technical Systems, Sandtorstrasse 1, 39106 Magdeburg, Germany; International Max Planck Research School (IMPRS) for Advanced Methods in Process and Systems Engineering, Magdeburg, Germany
| | - Matthias Stein
- Molecular Simulations and Design Group, Max Planck Institute for Dynamics of Complex Technical Systems, Sandtorstrasse 1, 39106 Magdeburg, Germany.
| |
Collapse
|
15
|
Bavi N, Bavi O, Vossoughi M, Naghdabadi R, Hill AP, Martinac B, Jamali Y. Nanomechanical properties of MscL α helices: A steered molecular dynamics study. Channels (Austin) 2016; 11:209-223. [PMID: 27753526 DOI: 10.1080/19336950.2016.1249077] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/20/2022] Open
Abstract
Gating of mechanosensitive (MS) channels is driven by a hierarchical cascade of movements and deformations of transmembrane helices in response to bilayer tension. Determining the intrinsic mechanical properties of the individual transmembrane helices is therefore central to understanding the intricacies of the gating mechanism of MS channels. We used a constant-force steered molecular dynamics (SMD) approach to perform unidirectional pulling tests on all the helices of MscL in M. tuberculosis and E. coli homologs. Using this method, we could overcome the issues encountered with the commonly used constant-velocity SMD simulations, such as low mechanical stability of the helix during stretching and high dependency of the elastic properties on the pulling rate. We estimated Young's moduli of the α-helices of MscL to vary between 0.2 and 12.5 GPa with TM2 helix being the stiffest. We also studied the effect of water on the properties of the pore-lining TM1 helix. In the absence of water, this helix exhibited a much stiffer response. By monitoring the number of hydrogen bonds, it appears that water acts like a 'lubricant' (softener) during TM1 helix elongation. These data shed light on another physical aspect underlying hydrophobic gating of MS channels, in particular MscL.
Collapse
Affiliation(s)
- N Bavi
- a Division of Molecular Cardiology and Biophysics , Victor Chang Cardiac Research Institute , Darlinghurst , NSW , Australia.,b St Vincent's Clinical School, Faculty of Medicine , University of New South Wales , Darlinghurst , NSW , Australia
| | - O Bavi
- c Institute for Nanoscience and Nanotechnology, Sharif University of Technology , Tehran , Iran
| | - M Vossoughi
- c Institute for Nanoscience and Nanotechnology, Sharif University of Technology , Tehran , Iran.,d Biochemical & Bioenvironmental Research Center (BBRC) , Tehran , Iran
| | - R Naghdabadi
- c Institute for Nanoscience and Nanotechnology, Sharif University of Technology , Tehran , Iran.,e Department of Mechanical Engineering , Sharif University of Technology , Tehran , Iran
| | - A P Hill
- a Division of Molecular Cardiology and Biophysics , Victor Chang Cardiac Research Institute , Darlinghurst , NSW , Australia
| | - B Martinac
- a Division of Molecular Cardiology and Biophysics , Victor Chang Cardiac Research Institute , Darlinghurst , NSW , Australia.,b St Vincent's Clinical School, Faculty of Medicine , University of New South Wales , Darlinghurst , NSW , Australia
| | - Y Jamali
- f Department of Mathematics , Tarbiat Modares University , Tehran , Iran.,g Computational Physical Sciences Research Laboratory , School of Nanoscience, Institute for Research in Fundamental Sciences (IPM) , Tehran , Iran
| |
Collapse
|
16
|
Soto-Delgado J, Tapia RA, Torras J. Multiscale Treatment for the Molecular Mechanism of a Diels-Alder Reaction in Solution: A QM/MM-MD Study. J Chem Theory Comput 2016; 12:4735-4742. [PMID: 27640742 DOI: 10.1021/acs.jctc.6b00772] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/12/2023]
Abstract
Thermodynamics and the solvent role in the acceleration of the Diels-Alder reaction between cyclopentadiene (CPD) and methyl vinyl ketone (MVK) have been revisited. In this work we use an ab initio hybrid QM/MM-MD scheme combined with multiple steered molecular dynamics to extract the free energy pofile in water and methanol using the bidirectional Minh-Adib estimator. We obtain 18.7 kcal mol-1 and 20.8 kcal mol-1 free energy barrier for the reaction in water and methanol, respectively. This methodology reproduces experimental values with an absolute error of about 0.8 kcal mol-1. The experimental difference between the activation free-energy barriers of water and methanol is also reproduced with an absolute error of about 0.1 kcal mol-1. We explore the charge transfer evolution along reaction coordinates to characterize the electronic behavior for this reaction. It is shown that the solvent molecules around the reaction system produce a global polarization along the reaction coordinate which is consistent with the solvent polarity. The results highlight the role of hydrogen bonding formed in the transition state to stabilize the system charge reorganization in the reaction process.
Collapse
Affiliation(s)
- Jorge Soto-Delgado
- Departamento de Ciencias Químicas, Facultad de Ciencias Exactas, Universidad Andrés Bello , Quillota 980, Viña del Mar, Chile
| | - Ricardo A Tapia
- Departamento de Química Orgánica, Pontificia Universidad Católica de Chile , Vicuña Mackenna, 4860 Santiago, Chile
| | - Juan Torras
- Departament d'Enginyeria Química, Escola d'Enginyeria de Barcelona Est (EEBE), Universitat Politècnica de Catalunya , C. Eduard Maristany 10-14, 08019 Barcelona, Spain
| |
Collapse
|
17
|
Noskov SY, Rostovtseva TK, Chamberlin AC, Teijido O, Jiang W, Bezrukov SM. Current state of theoretical and experimental studies of the voltage-dependent anion channel (VDAC). BIOCHIMICA ET BIOPHYSICA ACTA 2016; 1858:1778-90. [PMID: 26940625 PMCID: PMC4877207 DOI: 10.1016/j.bbamem.2016.02.026] [Citation(s) in RCA: 57] [Impact Index Per Article: 7.1] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 12/19/2015] [Revised: 02/09/2016] [Accepted: 02/10/2016] [Indexed: 01/04/2023]
Abstract
Voltage-dependent anion channel (VDAC), the major channel of the mitochondrial outer membrane provides a controlled pathway for respiratory metabolites in and out of the mitochondria. In spite of the wealth of experimental data from structural, biochemical, and biophysical investigations, the exact mechanisms governing selective ion and metabolite transport, especially the role of titratable charged residues and interactions with soluble cytosolic proteins, remain hotly debated in the field. The computational advances hold a promise to provide a much sought-after solution to many of the scientific disputes around solute and ion transport through VDAC and hence, across the mitochondrial outer membrane. In this review, we examine how Molecular Dynamics, Free Energy, and Brownian Dynamics simulations of the large β-barrel channel, VDAC, advanced our understanding. We will provide a short overview of non-conventional techniques and also discuss examples of how the modeling excursions into VDAC biophysics prospectively aid experimental efforts. This article is part of a Special Issue entitled: Membrane Proteins edited by J.C. Gumbart and Sergei Noskov.
Collapse
Affiliation(s)
- Sergei Yu Noskov
- Department of Biological Sciences and Centre for Molecular Simulation, University of Calgary, 2500 University Drive N.W., Calgary, Alberta T2N1N4, Canada.
| | - Tatiana K Rostovtseva
- Section on Molecular Transport, Eunice Kennedy Shriver National Institute of Child Health and Human Development, National Institutes of Health, Bethesda, MD 20892, USA.
| | | | - Oscar Teijido
- Section on Molecular Transport, Eunice Kennedy Shriver National Institute of Child Health and Human Development, National Institutes of Health, Bethesda, MD 20892, USA; Department of Medical Epigenetics, Institute of Medical Sciences and Genomic Medicine, EuroEspes Sta. Marta de Babío S/N, 15165 Bergondo, A Coruña, Spain
| | - Wei Jiang
- Leadership Computing Facility, Argonne National Laboratory, 9700S Cass Avenue, Lemont, IL 60439, USA
| | - Sergey M Bezrukov
- Section on Molecular Transport, Eunice Kennedy Shriver National Institute of Child Health and Human Development, National Institutes of Health, Bethesda, MD 20892, USA.
| |
Collapse
|