1
|
Hwang W, Austin SL, Blondel A, Boittier ED, Boresch S, Buck M, Buckner J, Caflisch A, Chang HT, Cheng X, Choi YK, Chu JW, Crowley MF, Cui Q, Damjanovic A, Deng Y, Devereux M, Ding X, Feig MF, Gao J, Glowacki DR, Gonzales JE, Hamaneh MB, Harder ED, Hayes RL, Huang J, Huang Y, Hudson PS, Im W, Islam SM, Jiang W, Jones MR, Käser S, Kearns FL, Kern NR, Klauda JB, Lazaridis T, Lee J, Lemkul JA, Liu X, Luo Y, MacKerell AD, Major DT, Meuwly M, Nam K, Nilsson L, Ovchinnikov V, Paci E, Park S, Pastor RW, Pittman AR, Post CB, Prasad S, Pu J, Qi Y, Rathinavelan T, Roe DR, Roux B, Rowley CN, Shen J, Simmonett AC, Sodt AJ, Töpfer K, Upadhyay M, van der Vaart A, Vazquez-Salazar LI, Venable RM, Warrensford LC, Woodcock HL, Wu Y, Brooks CL, Brooks BR, Karplus M. CHARMM at 45: Enhancements in Accessibility, Functionality, and Speed. J Phys Chem B 2024. [PMID: 39303207 DOI: 10.1021/acs.jpcb.4c04100] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 09/22/2024]
Abstract
Since its inception nearly a half century ago, CHARMM has been playing a central role in computational biochemistry and biophysics. Commensurate with the developments in experimental research and advances in computer hardware, the range of methods and applicability of CHARMM have also grown. This review summarizes major developments that occurred after 2009 when the last review of CHARMM was published. They include the following: new faster simulation engines, accessible user interfaces for convenient workflows, and a vast array of simulation and analysis methods that encompass quantum mechanical, atomistic, and coarse-grained levels, as well as extensive coverage of force fields. In addition to providing the current snapshot of the CHARMM development, this review may serve as a starting point for exploring relevant theories and computational methods for tackling contemporary and emerging problems in biomolecular systems. CHARMM is freely available for academic and nonprofit research at https://academiccharmm.org/program.
Collapse
Affiliation(s)
- Wonmuk Hwang
- Department of Biomedical Engineering, Texas A&M University, College Station, Texas 77843, United States
- Department of Materials Science and Engineering, Texas A&M University, College Station, Texas 77843, United States
- Department of Physics and Astronomy, Texas A&M University, College Station, Texas 77843, United States
- Center for AI and Natural Sciences, Korea Institute for Advanced Study, Seoul 02455, Republic of Korea
| | - Steven L Austin
- Department of Chemistry, University of South Florida, Tampa, Florida 33620, United States
| | - Arnaud Blondel
- Institut Pasteur, Université Paris Cité, CNRS UMR3825, Structural Bioinformatics Unit, 28 rue du Dr. Roux F-75015 Paris, France
| | - Eric D Boittier
- Department of Chemistry, University of Basel, Klingelbergstrasse 80, CH-4056 Basel, Switzerland
| | - Stefan Boresch
- Faculty of Chemistry, Department of Computational Biological Chemistry, University of Vienna, Wahringerstrasse 17, 1090 Vienna, Austria
| | - Matthias Buck
- Department of Physiology and Biophysics, Case Western Reserve University, School of Medicine, Cleveland, Ohio 44106, United States
| | - Joshua Buckner
- Department of Chemistry, University of Michigan, Ann Arbor, Michigan 48109, United States
| | - Amedeo Caflisch
- Department of Biochemistry, University of Zürich, CH-8057 Zürich, Switzerland
| | - Hao-Ting Chang
- Institute of Bioinformatics and Systems Biology, National Yang Ming Chiao Tung University, Hsinchu 30010, Taiwan, ROC
| | - Xi Cheng
- Shanghai Institute of Materia Medica, Chinese Academy of Sciences, Shanghai 201203, China
| | - Yeol Kyo Choi
- Department of Biological Sciences, Lehigh University, Bethlehem, Pennsylvania 18015, United States
| | - Jhih-Wei Chu
- Institute of Bioinformatics and Systems Biology, Department of Biological Science and Technology, Institute of Molecular Medicine and Bioengineering, and Center for Intelligent Drug Systems and Smart Bio-devices (IDS2B), National Yang Ming Chiao Tung University, Hsinchu 30010, Taiwan, ROC
| | - Michael F Crowley
- Renewable Resources and Enabling Sciences Center, National Renewable Energy Laboratory, Golden, Colorado 80401, United States
| | - Qiang Cui
- Department of Chemistry, Boston University, 590 Commonwealth Avenue, Boston, Massachusetts 02215, United States
- Department of Physics, Boston University, 590 Commonwealth Avenue, Boston, Massachusetts 02215, United States
- Department of Biomedical Engineering, Boston University, 44 Cummington Mall, Boston, Massachusetts 02215, United States
| | - Ana Damjanovic
- Department of Biophysics, Johns Hopkins University, Baltimore, Maryland 21218, United States
- Department of Physics and Astronomy, Johns Hopkins University, Baltimore, Maryland 21218, United States
- Laboratory of Computational Biology, National Heart Lung and Blood Institute, National Institutes of Health, Bethesda, Maryland 20892, United States
| | - Yuqing Deng
- Shanghai R&D Center, DP Technology, Ltd., Shanghai 201210, China
| | - Mike Devereux
- Department of Chemistry, University of Basel, Klingelbergstrasse 80, CH-4056 Basel, Switzerland
| | - Xinqiang Ding
- Department of Chemistry, Tufts University, Medford, Massachusetts 02155, United States
| | - Michael F Feig
- Department of Biochemistry and Molecular Biology, Michigan State University, East Lansing, Michigan 48824, United States
| | - Jiali Gao
- School of Chemical Biology & Biotechnology, Peking University Shenzhen Graduate School, Shenzhen, Guangdong 518055, China
- Institute of Systems and Physical Biology, Shenzhen Bay Laboratory, Shenzhen, Guangdong 518055, China
- Department of Chemistry and Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455, United States
| | - David R Glowacki
- CiTIUS Centro Singular de Investigación en Tecnoloxías Intelixentes da USC, 15705 Santiago de Compostela, Spain
| | - James E Gonzales
- Department of Biomedical Engineering, Texas A&M University, College Station, Texas 77843, United States
- Laboratory of Computational Biology, National Heart Lung and Blood Institute, National Institutes of Health, Bethesda, Maryland 20892, United States
| | - Mehdi Bagerhi Hamaneh
- Department of Physiology and Biophysics, Case Western Reserve University, School of Medicine, Cleveland, Ohio 44106, United States
| | | | - Ryan L Hayes
- Department of Chemical and Biomolecular Engineering, University of California, Irvine, Irvine, California 92697, United States
- Department of Pharmaceutical Sciences, University of California, Irvine, Irvine, California 92697, United States
| | - Jing Huang
- Key Laboratory of Structural Biology of Zhejiang Province, School of Life Sciences, Westlake University, Hangzhou, Zhejiang 310024, China
| | - Yandong Huang
- College of Computer Engineering, Jimei University, Xiamen 361021, China
| | - Phillip S Hudson
- Department of Chemistry, University of South Florida, Tampa, Florida 33620, United States
- Medicine Design, Pfizer Inc., Cambridge, Massachusetts 02139, United States
| | - Wonpil Im
- Department of Biological Sciences, Lehigh University, Bethlehem, Pennsylvania 18015, United States
| | - Shahidul M Islam
- Department of Chemistry, Delaware State University, Dover, Delaware 19901, United States
| | - Wei Jiang
- Computational Science Division, Argonne National Laboratory, Argonne, Illinois 60439, United States
| | - Michael R Jones
- Laboratory of Computational Biology, National Heart Lung and Blood Institute, National Institutes of Health, Bethesda, Maryland 20892, United States
| | - Silvan Käser
- Department of Chemistry, University of Basel, Klingelbergstrasse 80, CH-4056 Basel, Switzerland
| | - Fiona L Kearns
- Department of Chemistry, University of South Florida, Tampa, Florida 33620, United States
| | - Nathan R Kern
- Department of Biological Sciences, Lehigh University, Bethlehem, Pennsylvania 18015, United States
| | - Jeffery B Klauda
- Department of Chemical and Biomolecular Engineering, Institute for Physical Science and Technology, Biophysics Program, University of Maryland, College Park, Maryland 20742, United States
| | - Themis Lazaridis
- Department of Chemistry, City College of New York, New York, New York 10031, United States
| | - Jinhyuk Lee
- Disease Target Structure Research Center, Korea Research Institute of Bioscience and Biotechnology, Daejeon 34141, Republic of Korea
- Department of Bioinformatics, KRIBB School of Bioscience, University of Science and Technology, Daejeon 34141, Republic of Korea
| | - Justin A Lemkul
- Department of Biochemistry, Virginia Polytechnic Institute and State University, Blacksburg, Virginia 24061, United States
| | - Xiaorong Liu
- Department of Chemistry, University of Michigan, Ann Arbor, Michigan 48109, United States
| | - Yun Luo
- Department of Biotechnology and Pharmaceutical Sciences, College of Pharmacy, Western University of Health Sciences, Pomona, California 91766, United States
| | - Alexander D MacKerell
- Department of Pharmaceutical Sciences, University of Maryland School of Pharmacy, Baltimore, Maryland 21201, United States
| | - Dan T Major
- Department of Chemistry and Institute for Nanotechnology & Advanced Materials, Bar-Ilan University, Ramat-Gan 52900, Israel
| | - Markus Meuwly
- Department of Chemistry, University of Basel, Klingelbergstrasse 80, CH-4056 Basel, Switzerland
- Department of Chemistry, Brown University, Providence, Rhode Island 02912, United States
| | - Kwangho Nam
- Department of Chemistry and Biochemistry, University of Texas at Arlington, Arlington, Texas 76019, United States
| | - Lennart Nilsson
- Karolinska Institutet, Department of Biosciences and Nutrition, SE-14183 Huddinge, Sweden
| | - Victor Ovchinnikov
- Harvard University, Department of Chemistry and Chemical Biology, Cambridge, Massachusetts 02138, United States
| | - Emanuele Paci
- Dipartimento di Fisica e Astronomia, Universitá di Bologna, Bologna 40127, Italy
| | - Soohyung Park
- Department of Biological Sciences, Lehigh University, Bethlehem, Pennsylvania 18015, United States
| | - Richard W Pastor
- Laboratory of Computational Biology, National Heart Lung and Blood Institute, National Institutes of Health, Bethesda, Maryland 20892, United States
| | - Amanda R Pittman
- Department of Chemistry, University of South Florida, Tampa, Florida 33620, United States
| | - Carol Beth Post
- Borch Department of Medicinal Chemistry and Molecular Pharmacology, Purdue University, West Lafayette, Indiana 47907, United States
| | - Samarjeet Prasad
- Laboratory of Computational Biology, National Heart Lung and Blood Institute, National Institutes of Health, Bethesda, Maryland 20892, United States
| | - Jingzhi Pu
- Department of Chemistry and Chemical Biology, Indiana University Indianapolis, Indianapolis, Indiana 46202, United States
| | - Yifei Qi
- School of Pharmacy, Fudan University, Shanghai 201203, China
| | | | - Daniel R Roe
- Laboratory of Computational Biology, National Heart Lung and Blood Institute, National Institutes of Health, Bethesda, Maryland 20892, United States
| | - Benoit Roux
- Department of Chemistry, University of Chicago, Chicago, Illinois 60637, United States
| | | | - Jana Shen
- Department of Pharmaceutical Sciences, University of Maryland School of Pharmacy, Baltimore, Maryland 21201, United States
| | - Andrew C Simmonett
- Laboratory of Computational Biology, National Heart Lung and Blood Institute, National Institutes of Health, Bethesda, Maryland 20892, United States
| | - Alexander J Sodt
- Eunice Kennedy Shriver National Institute of Child Health and Human Development, National Institutes of Health, Bethesda, Maryland 20892, United States
| | - Kai Töpfer
- Department of Chemistry, University of Basel, Klingelbergstrasse 80, CH-4056 Basel, Switzerland
| | - Meenu Upadhyay
- Department of Chemistry, University of Basel, Klingelbergstrasse 80, CH-4056 Basel, Switzerland
| | - Arjan van der Vaart
- Department of Chemistry, University of South Florida, Tampa, Florida 33620, United States
| | | | - Richard M Venable
- Laboratory of Computational Biology, National Heart Lung and Blood Institute, National Institutes of Health, Bethesda, Maryland 20892, United States
| | - Luke C Warrensford
- Department of Chemistry, University of South Florida, Tampa, Florida 33620, United States
| | - H Lee Woodcock
- Department of Chemistry, University of South Florida, Tampa, Florida 33620, United States
| | - Yujin Wu
- Department of Chemistry, University of Michigan, Ann Arbor, Michigan 48109, United States
| | - Charles L Brooks
- Department of Chemistry, University of Michigan, Ann Arbor, Michigan 48109, United States
| | - Bernard R Brooks
- Laboratory of Computational Biology, National Heart Lung and Blood Institute, National Institutes of Health, Bethesda, Maryland 20892, United States
| | - Martin Karplus
- Harvard University, Department of Chemistry and Chemical Biology, Cambridge, Massachusetts 02138, United States
- Laboratoire de Chimie Biophysique, ISIS, Université de Strasbourg, 67000 Strasbourg, France
| |
Collapse
|
2
|
Su J, Wang P, Zhou W, Peydayesh M, Zhou J, Jin T, Donat F, Jin C, Xia L, Wang K, Ren F, Van der Meeren P, García de Arquer FP, Mezzenga R. Single-site iron-anchored amyloid hydrogels as catalytic platforms for alcohol detoxification. NATURE NANOTECHNOLOGY 2024; 19:1168-1177. [PMID: 38740933 PMCID: PMC11329373 DOI: 10.1038/s41565-024-01657-7] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 10/10/2023] [Accepted: 03/21/2024] [Indexed: 05/16/2024]
Abstract
Constructing effective antidotes to reduce global health impacts induced by alcohol prevalence is a challenging topic. Despite the positive effects observed with intravenous applications of natural enzyme complexes, their insufficient activities and complicated usage often result in the accumulation of toxic acetaldehyde, which raises important clinical concerns, highlighting the pressing need for stable oral strategies. Here we present an effective solution for alcohol detoxification by employing a biomimetic-nanozyme amyloid hydrogel as an orally administered catalytic platform. We exploit amyloid fibrils derived from β-lactoglobulin, a readily accessible milk protein that is rich in coordinable nitrogen atoms, as a nanocarrier to stabilize atomically dispersed iron (ferrous-dominated). By emulating the coordination structure of the horseradish peroxidase enzyme, the single-site iron nanozyme demonstrates the capability to selectively catalyse alcohol oxidation into acetic acid, as opposed to the more toxic acetaldehyde. Administering the gelatinous nanozyme to mice suffering from alcohol intoxication significantly reduced their blood-alcohol levels (decreased by 55.8% 300 min post-alcohol intake) without causing additional acetaldehyde build-up. Our hydrogel further demonstrates a protective effect on the liver, while simultaneously mitigating intestinal damage and dysbiosis associated with chronic alcohol consumption, introducing a promising strategy in effective alcohol detoxification.
Collapse
Affiliation(s)
- Jiaqi Su
- Department of Health Sciences and Technology, ETH Zurich, Zurich, Switzerland.
- Particle and Interfacial Technology Group, Faculty of Bioscience Engineering, Ghent University, Ghent, Belgium.
| | - Pengjie Wang
- Department of Nutrition and Health, Beijing Higher Institution Engineering Research Center of Animal Products, China Agricultural University, Beijing, China
| | - Wei Zhou
- Department of Chemistry and Applied Biosciences, ETH Zurich, Zurich, Switzerland
| | - Mohammad Peydayesh
- Department of Health Sciences and Technology, ETH Zurich, Zurich, Switzerland
| | - Jiangtao Zhou
- Department of Health Sciences and Technology, ETH Zurich, Zurich, Switzerland
| | - Tonghui Jin
- Department of Health Sciences and Technology, ETH Zurich, Zurich, Switzerland
| | - Felix Donat
- Institute of Energy and Process Engineering, Department of Mechanical and Process Engineering, ETH Zurich, Zurich, Switzerland
| | - Cuiyuan Jin
- Institute of Translational Medicine, Zhejiang Shuren University, Zhejiang, China
| | - Lu Xia
- ICFO-Institut de Ciències Fotòniques, The Barcelona Institute of Science and Technology, Barcelona, Spain
| | - Kaiwen Wang
- ICFO-Institut de Ciències Fotòniques, The Barcelona Institute of Science and Technology, Barcelona, Spain
| | - Fazheng Ren
- Department of Nutrition and Health, Beijing Higher Institution Engineering Research Center of Animal Products, China Agricultural University, Beijing, China
| | - Paul Van der Meeren
- Particle and Interfacial Technology Group, Faculty of Bioscience Engineering, Ghent University, Ghent, Belgium
| | - F Pelayo García de Arquer
- ICFO-Institut de Ciències Fotòniques, The Barcelona Institute of Science and Technology, Barcelona, Spain
| | - Raffaele Mezzenga
- Department of Health Sciences and Technology, ETH Zurich, Zurich, Switzerland.
- Department of Materials, ETH Zurich, Zurich, Switzerland.
| |
Collapse
|
3
|
Guvench O. Effect of Lipid Bilayer Anchoring on the Conformational Properties of the Cytochrome P450 2D6 Binding Site. J Phys Chem B 2024; 128:7188-7198. [PMID: 39016537 DOI: 10.1021/acs.jpcb.4c03097] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 07/18/2024]
Abstract
Human cytochrome P450 (CYP) proteins metabolize 75% of small-molecule pharmaceuticals, which makes structure-based modeling of CYP metabolism and inhibition, bolstered by the current availability of X-ray crystal structures of CYP globular catalytic domains, an attractive prospect. Accounting for this broad metabolic capacity is a combination of the existence of multiple different CYP proteins and the capacity of a single CYP protein to metabolize multiple different small molecules. It is thought that structural plasticity and flexibility contribute to this latter property; therefore, incorporating diverse conformational states of a particular CYP is likely an important consideration in structure-based CYP metabolism and inhibition modeling. While all-atom explicit-solvent molecular dynamics simulations can be used to generate conformational ensembles under biologically relevant conditions, existing CYP crystal structures are of the globular domain only, whereas human CYPs contain N-terminal transmembrane and linker peptides that anchor the globular catalytic domain to the endoplasmic reticulum. To determine whether this can cause significant differences in the sampled binding site conformations, microsecond scale all-atom explicit-solvent molecular dynamics simulations of the CYP2D6 globular domain in an aqueous environment were compared with those of the full-length protein anchored in a POPC lipid bilayer. While bilayer-anchoring damped some structural fluctuations in the globular domain relative to the aqueous simulations, none of the affected residues included binding site pocket residues. Furthermore, clustering of molecular dynamics snapshots based on either pairwise binding site pocket RMSD or volume differences demonstrated a lack of separation of snapshots from the two simulation conditions into different clusters. These results suggest the substantially simpler and computationally cheaper aqueous simulation approach can be used to generate a relevant conformational ensemble of the CYP2D6 binding site for structure-based metabolism and inhibition modeling.
Collapse
Affiliation(s)
- Olgun Guvench
- Department of Pharmaceutical Sciences and Administration, School of Pharmacy, Westbrook College of Health Professions, University of New England, 716 Stevens Ave, Portland, Maine 04103, United States
| |
Collapse
|
4
|
D'Agostino M, Simonetti A, Motta S, Wolff P, Romagnoli A, Piccinini A, Spinozzi F, Di Marino D, La Teana A, Ennifar E. Crystal structure of archaeal IF5A-DHS complex reveals insights into the hypusination mechanism. Structure 2024; 32:878-888.e4. [PMID: 38582076 DOI: 10.1016/j.str.2024.03.008] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/18/2023] [Revised: 02/12/2024] [Accepted: 03/12/2024] [Indexed: 04/08/2024]
Abstract
The translation factor IF5A is highly conserved in Eukarya and Archaea and undergoes a unique post-translational hypusine modification by the deoxyhypusine synthase (DHS) enzyme. DHS transfers the butylamine moiety from spermidine to IF5A using NAD as a cofactor, forming a deoxyhypusine intermediate. IF5A is a key player in protein synthesis, preventing ribosome stalling in proline-rich sequences during translation elongation and facilitating translation elongation and termination. Additionally, human eIF5A participates in various essential cellular processes and contributes to cancer metastasis, with inhibiting hypusination showing anti-proliferative effects. The hypusination pathway of IF5A is therefore an attractive new therapeutic target. We elucidated the 2.0 Å X-ray crystal structure of the archaeal DHS-IF5A complex, revealing hetero-octameric architecture and providing a detailed view of the complex active site including the hypusination loop. This structure, along with biophysical data and molecular dynamics simulations, provides new insights into the catalytic mechanism of the hypusination reaction.
Collapse
Affiliation(s)
- Mattia D'Agostino
- Department of Life and Environmental Sciences, Polytechnic University of Marche, Via Brecce Bianche, 60131 Ancona, Italy; Architecture et Réactivité de l'ARN, CNRS UPR 9002, Institut de Biologie Moléculaire et Cellulaire, Université de Strasbourg, Strasbourg, France
| | - Angelita Simonetti
- Architecture et Réactivité de l'ARN, CNRS UPR 9002, Institut de Biologie Moléculaire et Cellulaire, Université de Strasbourg, Strasbourg, France
| | - Stefano Motta
- Department of Earth and Environmental Sciences, University of Milano Bicocca, Piazza della Scienza 1, 20126 Milan, Italy
| | - Philippe Wolff
- Architecture et Réactivité de l'ARN, CNRS UPR 9002, Institut de Biologie Moléculaire et Cellulaire, Université de Strasbourg, Strasbourg, France
| | - Alice Romagnoli
- Department of Life and Environmental Sciences, Polytechnic University of Marche, Via Brecce Bianche, 60131 Ancona, Italy; New York-Marche Structural Biology Center (Ny-Masbic), Polytechnic University of Marche, Via Brecce Bianche, 60131 Ancona, Italy
| | - Astra Piccinini
- Department of Life and Environmental Sciences, Polytechnic University of Marche, Via Brecce Bianche, 60131 Ancona, Italy
| | - Francesco Spinozzi
- Department of Life and Environmental Sciences, Polytechnic University of Marche, Via Brecce Bianche, 60131 Ancona, Italy
| | - Daniele Di Marino
- Department of Life and Environmental Sciences, Polytechnic University of Marche, Via Brecce Bianche, 60131 Ancona, Italy; New York-Marche Structural Biology Center (Ny-Masbic), Polytechnic University of Marche, Via Brecce Bianche, 60131 Ancona, Italy; Neuronal Death and Neuroprotection Unit, Department of Neuroscience, Mario Negri Institute for Pharmacological Research-IRCCS, Via Mario Negri 2, 20156 Milano, Italy.
| | - Anna La Teana
- Department of Life and Environmental Sciences, Polytechnic University of Marche, Via Brecce Bianche, 60131 Ancona, Italy; New York-Marche Structural Biology Center (Ny-Masbic), Polytechnic University of Marche, Via Brecce Bianche, 60131 Ancona, Italy.
| | - Eric Ennifar
- Architecture et Réactivité de l'ARN, CNRS UPR 9002, Institut de Biologie Moléculaire et Cellulaire, Université de Strasbourg, Strasbourg, France.
| |
Collapse
|
5
|
Nan Y, MacKerell AD. Balancing Group I Monatomic Ion-Polar Compound Interactions for Condensed Phase Simulation in the Polarizable Drude Force Field. J Chem Theory Comput 2024; 20:3242-3257. [PMID: 38588064 PMCID: PMC11039353 DOI: 10.1021/acs.jctc.3c01380] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/10/2024]
Abstract
Molecular dynamics (MD) simulations are a commonly used method for investigating molecular behavior at the atomic level. Achieving reliable MD simulation results necessitates the use of an accurate force field. In the present work, we present a protocol to enhance the quality of group 1 monatomic ions (specifically Li+, Na+, K+, Rb+, and Cs+) with respect to their interactions with common polar model compounds in biomolecules in condensed phases in the context of the Drude polarizable force field. Instead of adjusting preexisting individual parameters for ions, model compounds, and water, we employ atom-pair specific Lennard-Jones (LJ) (known as NBFIX in CHARMM) and through-space Thole dipole screening (NBTHOLE) terms to fine-tune the balance of ion-model compound, ion-water, and model compound-water interactions. This involved establishing a protocol for the optimization of NBFIX and NBTHOLE parameters targeting the difference between molecular mechanical (MM) and quantum mechanical (QM) potential energy scans (PES). It is shown that targeting PES involving complexes that include multiple model compounds and/or ions as trimers and tetramers yields parameters that produce condensed phase properties in agreement with experimental data. Validation of this protocol involved the reproduction of experimental thermodynamic benchmarks, including solvation free energies of ions in methanol and N-methylacetamide, osmotic pressures, ionic conductivities, and diffusion coefficients within the condensed phase. These results show the importance of including more complex ion-model compound complexes beyond dimers in the QM target data to account for many-body effects during parameter fitting. The presented parameters represent a significant refinement of the Drude polarizable force field, which will lead to improved accuracy for modeling ion-biomolecular interactions.
Collapse
Affiliation(s)
- Yiling Nan
- University of Maryland Computer-Aided Drug Design Center, Department of Pharmaceutical Sciences, School of Pharmacy, University of Maryland, Baltimore, MD 21201 MD
| | - Alexander D. MacKerell
- University of Maryland Computer-Aided Drug Design Center, Department of Pharmaceutical Sciences, School of Pharmacy, University of Maryland, Baltimore, MD 21201 MD
| |
Collapse
|
6
|
Li S, Wu B, Luo YL, Han W. Simulations of Functional Motions of Super Large Biomolecules with a Mixed-Resolution Model. J Chem Theory Comput 2024; 20:2228-2245. [PMID: 38374639 PMCID: PMC10938502 DOI: 10.1021/acs.jctc.3c01046] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/22/2023] [Revised: 01/18/2024] [Accepted: 01/29/2024] [Indexed: 02/21/2024]
Abstract
Many large protein machines function through an interplay between large-scale movements and intricate conformational changes. Understanding functional motions of these proteins through simulations becomes challenging for both all-atom and coarse-grained (CG) modeling techniques because neither approach alone can readily capture the full details of these motions. In this study, we develop a multiscale model by employing the popular MARTINI CG model to represent a heterogeneous environment and structurally stable proteins and using the united-atom (UA) model PACE to describe proteins undergoing subtle conformational changes. PACE was previously developed to be compatible with the MARTINI solvent and membrane. Here, we couple the protein descriptions of the two models by directly mixing UA and CG interaction parameters to greatly simplify parameter determination. Through extensive validations with diverse protein systems in solution or membrane, we demonstrate that only additional parameter rescaling is needed to enable the resulting model to recover the stability of native structures of proteins under mixed representation. Moreover, we identify the optimal scaling factors that can be applied to various protein systems, rendering the model potentially transferable. To further demonstrate its applicability for realistic systems, we apply the model to a mechanosensitive ion channel Piezo1 that has peripheral arms for sensing membrane tension and a central pore for ion conductance. The model can reproduce the coupling between Piezo1's large-scale arm movement and subtle pore opening in response to membrane stress while consuming much less computational costs than all-atom models. Therefore, our model shows promise for studying functional motions of large protein machines.
Collapse
Affiliation(s)
- Shu Li
- Centre
for Artificial Intelligence Driven Drug Discovery, Faculty of Applied
Sciences, Macao Polytechnic University, Macao 999078, China
- State
Key Laboratory of Chemical Oncogenomics, Guangdong Provincial Key
Laboratory of Chemical Genomics, School of Chemical Biology and Biotechnology, Peking University Shenzhen Graduate School, Shenzhen 518055, China
| | - Bohua Wu
- State
Key Laboratory of Chemical Oncogenomics, Guangdong Provincial Key
Laboratory of Chemical Genomics, School of Chemical Biology and Biotechnology, Peking University Shenzhen Graduate School, Shenzhen 518055, China
| | - Yun Lyna Luo
- Department
of Biotechnology and Pharmaceutical Sciences, Western University of Health Sciences, Pomona, California 91766, United States
| | - Wei Han
- State
Key Laboratory of Chemical Oncogenomics, Guangdong Provincial Key
Laboratory of Chemical Genomics, School of Chemical Biology and Biotechnology, Peking University Shenzhen Graduate School, Shenzhen 518055, China
- Department
of Chemistry, Faculty of Science, Hong Kong
Baptist University, Hong Kong SAR 999077, China
- Shenzhen
Bay Laboratory, Institute of Chemical Biology, Shenzhen 518132, China
| |
Collapse
|
7
|
Xie M, Zhang M, Jin Z. Machine Learning-Based Interfacial Tension Equations for (H 2 + CO 2)-Water/Brine Systems over a Wide Range of Temperature and Pressure. LANGMUIR : THE ACS JOURNAL OF SURFACES AND COLLOIDS 2024; 40:5369-5377. [PMID: 38417158 DOI: 10.1021/acs.langmuir.3c03831] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 03/01/2024]
Abstract
Large-scale underground hydrogen storage (UHS) plays a vital role in energy transition. H2-brine interfacial tension (IFT) is a crucial parameter in structural trapping in underground geological locations and gas-water two-phase flow in subsurface porous media. On the other hand, cushion gas, such as CO2, is often co-injected with H2 to retain reservoir pressure. Therefore, it is imperative to accurately predict the (H2 + CO2)-water/brine IFT under UHS conditions. While there have been a number of experimental measurements on H2-water/brine and (H2 + CO2)-water/brine IFT, an accurate and efficient (H2 + CO2)-water/brine IFT model under UHS conditions is still lacking. In this work, we use molecular dynamics (MD) simulations to generate an extensive (H2 + CO2)-water/brine IFT databank (840 data points) over a wide range of temperature (from 298 to 373 K), pressure (from 50 to 400 bar), gas composition, and brine salinity (up to 3.15 mol/kg) for typical UHS conditions, which is used to develop an accurate and efficient machine learning (ML)-based IFT equation. Our ML-based IFT equation is validated by comparing to available experimental data and other IFT equations for various systems (H2-brine/water, CO2-brine/water, and (H2 + CO2)-brine/water), rendering generally good performance (with R2 = 0.902 against 601 experimental data points). The developed ML-based IFT equation can be readily applied and implemented in reservoir simulations and other UHS applications.
Collapse
Affiliation(s)
- Minjunshi Xie
- School of Mining and Petroleum Engineering, Department of Civil and Environmental Engineering, University of Alberta, Edmonton, AB T6G 1H9, Canada
| | - Mingshan Zhang
- Key Laboratory of Ministry of Education on Safe Mining of Deep Metal Mines, School of Resources and Civil Engineering, Northeastern University, Shenyang 110819, China
- Key Laboratory of Liaoning Province on Deep Engineering and Intelligent Technology, Northeastern University, Shenyang 110819, China
| | - Zhehui Jin
- School of Mining and Petroleum Engineering, Department of Civil and Environmental Engineering, University of Alberta, Edmonton, AB T6G 1H9, Canada
| |
Collapse
|
8
|
Kinose K, Shinoda K, Konishi T, Kawasaki H. Mutational analysis in Corynebacterium stationis MFS transporters for improving nucleotide bioproduction. Appl Microbiol Biotechnol 2024; 108:251. [PMID: 38436751 PMCID: PMC10912292 DOI: 10.1007/s00253-024-13080-y] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/04/2023] [Revised: 02/16/2024] [Accepted: 02/19/2024] [Indexed: 03/05/2024]
Abstract
Product secretion from an engineered cell can be advantageous for microbial cell factories. Extensive work on nucleotide manufacturing, one of the most successful microbial fermentation processes, has enabled Corynebacterium stationis to transport nucleotides outside the cell by random mutagenesis; however, the underlying mechanism has not been elucidated, hindering its applications in transporter engineering. Herein, we report the nucleotide-exporting major facilitator superfamily (MFS) transporter from the C. stationis genome and its hyperactive mutation at the G64 residue. Structural estimation and molecular dynamics simulations suggested that the activity of this transporter improved via two mechanisms: (1) enhancing interactions between transmembrane helices through the conserved "RxxQG" motif along with substrate binding and (2) trapping substrate-interacting residue for easier release from the cavity. Our results provide novel insights into how MFS transporters change their conformation from inward- to outward-facing states upon substrate binding to facilitate efflux and can contribute to the development of rational design approaches for efflux improvements in microbial cell factories. KEYPOINTS: • An MFS transporter from C. stationis genome and its mutation at residue G64 were assessed • It enhanced the transporter activity by strengthening transmembrane helix interactions and trapped substrate-interacting residues • Our results contribute to rational design approach development for efflux improvement.
Collapse
Affiliation(s)
- Keita Kinose
- Agro-Biotechnology Research Center, Graduate School of Agriculture and Life Sciences, The University of Tokyo, Tokyo, Japan
- Nagahama Institute for Biochemical Science, Oriental Yeast Co., Ltd., Nagahama, Shiga, Japan
| | - Keiko Shinoda
- Agro-Biotechnology Research Center, Graduate School of Agriculture and Life Sciences, The University of Tokyo, Tokyo, Japan
- Collaborative Research Institute for Innovative Microbiology, The University of Tokyo, Tokyo, Japan
- Research Organization of Information and Systems, The Institute of Statistical Mathematics, Tachikawa, Japan
| | - Tomoyuki Konishi
- Agro-Biotechnology Research Center, Graduate School of Agriculture and Life Sciences, The University of Tokyo, Tokyo, Japan
| | - Hisashi Kawasaki
- Agro-Biotechnology Research Center, Graduate School of Agriculture and Life Sciences, The University of Tokyo, Tokyo, Japan.
- Collaborative Research Institute for Innovative Microbiology, The University of Tokyo, Tokyo, Japan.
| |
Collapse
|
9
|
Paschek D, Busch J, Mock E, Ludwig R, Strate A. Computing the frequency-dependent NMR relaxation of 1H nuclei in liquid water. J Chem Phys 2024; 160:074102. [PMID: 38364003 DOI: 10.1063/5.0191052] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/11/2023] [Accepted: 01/23/2024] [Indexed: 02/18/2024] Open
Abstract
We present a computational framework for reliably determining the frequency-dependent intermolecular and intramolecular nuclear magnetic resonance (NMR) dipole-dipole relaxation rates of spin 1/2 nuclei from Molecular Dynamics (MD) simulations. This approach avoids the alterations caused by the well-known finite-size effects of translational diffusion. Moreover, a procedure is derived to control and correct for effects caused by fixed distance-sampling cutoffs and periodic boundary conditions. By construction, this approach is capable of accurately predicting the correct low-frequency scaling behavior of the intermolecular NMR dipole-dipole relaxation rate and thus allows for the reliable calculation of the frequency-dependent relaxation rate over many orders of magnitude. Our approach is based on the utilization of the theory of Hwang and Freed for the intermolecular dipole-dipole correlation function and its corresponding spectral density [L.-P. Hwang and J. H. Freed, J. Chem. Phys. 63, 4017-4025 (1975)] and its combination with data from MD simulations. The deviations from the Hwang and Freed theory caused by periodic boundary conditions and sampling distance cutoffs are quantified by means of random walker Monte Carlo simulations. An expression based on the Hwang and Freed theory is also suggested for correcting those effects. As a proof of principle, our approach is demonstrated by computing the frequency-dependent intermolecular and intramolecular dipolar NMR relaxation rates of 1H nuclei in liquid water at 273 and 298 K based on the simulations of the TIP4P/2005 model. Our calculations are suggesting that the intermolecular contribution to the 1H NMR relaxation rate of the TIP4P/2005 model in the extreme narrowing limit has previously been substantially underestimated.
Collapse
Affiliation(s)
- Dietmar Paschek
- Institut für Chemie, Abteilung Physikalische und Theoretische Chemie, Universität Rostock, Albert-Einstein-Str. 27, D-18059 Rostock, Germany
| | - Johanna Busch
- Institut für Chemie, Abteilung Physikalische und Theoretische Chemie, Universität Rostock, Albert-Einstein-Str. 27, D-18059 Rostock, Germany
| | - Eduard Mock
- Institut für Chemie, Abteilung Physikalische und Theoretische Chemie, Universität Rostock, Albert-Einstein-Str. 27, D-18059 Rostock, Germany
| | - Ralf Ludwig
- Institut für Chemie, Abteilung Physikalische und Theoretische Chemie, Universität Rostock, Albert-Einstein-Str. 27, D-18059 Rostock, Germany
- Leibniz-Institut für Katalyse an der Universität Rostock e.V., Albert-Einstein-Str. 29a, D-18059 Rostock, Germany
- Department Life, Light & Matter, Universität Rostock, Albert-Einstein-Str. 25, D-18059 Rostock, Germany
| | - Anne Strate
- Institut für Chemie, Abteilung Physikalische und Theoretische Chemie, Universität Rostock, Albert-Einstein-Str. 27, D-18059 Rostock, Germany
| |
Collapse
|
10
|
Busch J, Paschek D. Computing Accurate True Self-Diffusion Coefficients and Shear Viscosities Using the OrthoBoXY Approach. J Phys Chem B 2024; 128:1040-1052. [PMID: 38240259 DOI: 10.1021/acs.jpcb.3c07540] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/02/2024]
Abstract
In a recent paper [Busch, J.; Paschek, D. J. Phys. Chem. B 2023, 127, 7983-7987], we have shown that for molecular dynamics (MD) simulations using orthorhombic periodic boundary conditions with "magic" box length ratios of Lz/Lx = Lz/Ly = 2.7933596497, the self-diffusion coefficients Dx and Dy in x- and y-directions are independent of the system size. They both represent the true self-diffusion coefficient D0 = (Dx + Dy)/2, while the shear viscosity can be calculated from diffusion coefficients in x-, y-, and z-directions, using η = kBT·8.1711245653/[3πLz(Dx + Dy - 2Dz)]. In this contribution, we test this "OrthoBoXY" approach by its application to a variety of different systems: liquid water, dimethyl ether, methanol, triglyme, water/methanol mixtures, water/triglyme mixtures, and imidazolium-based ionic liquids. The chosen systems range from small-sized molecular liquids to complex mixtures and ionic liquids, while spanning a viscosity range of almost 3 orders of magnitude. We assess the efficiency of the method for computing true self-diffusion and viscosity data and provide simple formulas for estimating the required MD simulation lengths and sizes for delivering reliable data with targeted uncertainty levels. Our analysis of the system size dependence of statistical uncertainties for both the viscosity and the self-diffusion coefficient leads us to the conclusion that it is preferable to extend the simulation length instead of increasing the system size. MD simulations consisting of 768 molecules or ion pairs seem to be perfectly adequate.
Collapse
Affiliation(s)
- Johanna Busch
- Institut für Chemie, Abteilung Physikalische und Theoretische Chemie, Universität Rostock, Albert-Einstein-Str. 27, D-18059 Rostock, Germany
| | - Dietmar Paschek
- Institut für Chemie, Abteilung Physikalische und Theoretische Chemie, Universität Rostock, Albert-Einstein-Str. 27, D-18059 Rostock, Germany
| |
Collapse
|
11
|
Brooks CL, MacKerell AD, Post CB, Nilsson L. Biomolecular dynamics in the 21st century. Biochim Biophys Acta Gen Subj 2024; 1868:130534. [PMID: 38065235 PMCID: PMC10842176 DOI: 10.1016/j.bbagen.2023.130534] [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/26/2023] [Revised: 11/28/2023] [Accepted: 11/29/2023] [Indexed: 01/03/2024]
Abstract
The relevance of motions in biological macromolecules has been clear since the early structural analyses of proteins by X-ray crystallography. Computer simulations have been applied to provide a deeper understanding of the dynamics of biological macromolecules since 1976, and are now a standard tool in many labs working on the structure and function of biomolecules. In this mini-review we highlight some areas of current interest and active development for simulations, in particular all-atom molecular dynamics simulations.
Collapse
Affiliation(s)
- Charles L Brooks
- University of Michigan, Department of Chemistry, Ann Arbor, MI 48109, USA.
| | | | - Carol B Post
- Purdue University, Department of Medicinal Chemistry and Molecular Pharmacology, West Lafayette, IN 47907-2091, USA.
| | - Lennart Nilsson
- Karolinska Institutet, Department of Biosciences and Nutrition, SE-14183 Huddinge, Sweden.
| |
Collapse
|
12
|
Busch J, Paschek D. An OrthoBoXY-method for various alternative box geometries. Phys Chem Chem Phys 2024; 26:2907-2914. [PMID: 38086638 DOI: 10.1039/d3cp04916g] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/25/2024]
Abstract
We have shown in a recent contribution [Busch and Paschek, J. Phys. Chem. B, 2023 127, 7983-7987] that for molecular dynamics (MD) simulations of isotropic fluids based on orthorhombic periodic boundary conditions with "magical" box length ratios of Lz/Lx = Lz/Ly = 2.7933596497, the computed self-diffusion coefficients Dx and Dy in x- and y-direction become system size independent. They thus represent the true self-diffusion coefficient D0 = (Dx + Dy)/2, while the shear viscosity can be determined from diffusion coefficients in x-, y-, and z-direction, using the expression η = kBT·8.1711245653/[3πLz(Dx + Dy - 2Dz)]. Here we present a more generalized version of this "OrthoBoXY"-approach, which can be applied to any orthorhombic MD box of any shape. In particular, we would like to test, how the efficiency is affected by using a shape more akin to the cubic form, albeit with different box length ratios Lx/Lz ≠ Ly/Lz and Lx < Ly < Lz. We use NVT and NpT simulations of systems of 1536 TIP4P/2005 water molecules as a benchmark and explore different box geometries to determine the influence of the box shape on the computed statistical uncertainties for D0 and η. Moreover, another "magical" set of box length ratios is discovered with Ly/Lz = 0.57804765578 and Lx/Lz = 0.33413909235, where the self-diffusion coefficient in x-direction becomes system size independent, such that D0 = Dx.
Collapse
Affiliation(s)
- Johanna Busch
- Institut für Chemie, Abteilung Physikalische und Theoretische Chemie, Universität Rostock, Albert-Einstein-Straße 27, D-18059 Rostock, Germany.
| | - Dietmar Paschek
- Institut für Chemie, Abteilung Physikalische und Theoretische Chemie, Universität Rostock, Albert-Einstein-Straße 27, D-18059 Rostock, Germany.
| |
Collapse
|
13
|
Guvench O. Water Exchange from the Buried Binding Sites of Cytochrome P450 Enzymes 1A2, 2D6, and 3A4 Correlates with Conformational Fluctuations. Molecules 2024; 29:494. [PMID: 38276571 PMCID: PMC10820051 DOI: 10.3390/molecules29020494] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/19/2023] [Revised: 01/11/2024] [Accepted: 01/16/2024] [Indexed: 01/27/2024] Open
Abstract
Human cytochrome P450 enzymes (CYPs) are critical for the metabolism of small-molecule pharmaceuticals (drugs). As such, the prediction of drug metabolism by and drug inhibition of CYP activity is an important component of the drug discovery and design process. Relative to the availability of a wide range of experimental atomic-resolution CYP structures, the development of structure-based CYP activity models has been limited. To better characterize the role of CYP conformational fluctuations in CYP activity, we perform multiple microsecond-scale all-atom explicit-solvent molecular dynamics (MD) simulations on three CYP isoforms, 1A2, 2D6, and 3A4, which together account for the majority of CYP-mediated drug metabolism. The MD simulations employ a variety of positional restraints, ranging from keeping all CYP atoms close to their experimentally determined coordinates to allowing full flexibility. We find that, with full flexibility, large fluctuations in the CYP binding sites correlate with efficient water exchange from these buried binding sites. This is especially true for 1A2, which, when restrained to its crystallographic conformation, is unable to exchange water between the binding site and bulk solvent. These findings imply that, in addition to crystal structures, a representative ensemble of conformational states ought to be included when developing structure-based CYP activity models.
Collapse
Affiliation(s)
- Olgun Guvench
- Department of Pharmaceutical Sciences and Administration, School of Pharmacy, Westbrook College of Health Professions, University of New England, 716 Stevens Avenue, Portland, ME 04103, USA
| |
Collapse
|
14
|
Alamri SH, Haque S, Alghamdi BS, Tayeb HO, Azhari S, Farsi RM, Elmokadem A, Alamri TA, Harakeh S, Prakash A, Kumar V. Comprehensive mapping of mutations in TDP-43 and α-Synuclein that affect stability and binding. J Biomol Struct Dyn 2023:1-13. [PMID: 38126188 DOI: 10.1080/07391102.2023.2293258] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/24/2023] [Accepted: 11/11/2023] [Indexed: 12/23/2023]
Abstract
Abnormal aggregation and amyloid inclusions of TAR DNA-binding protein 43 (TDP-43) and α-Synuclein (α-Syn) are frequently co-observed in amyotrophic lateral sclerosis, Parkinson's disease, and Alzheimer's disease. Several reports showed TDP-43 C-terminal domain (CTD) and α-Syn interact with each other and the aggregates of these two proteins colocalized together in different cellular and animal models. Molecular dynamics simulation was conducted to elucidate the stability of the TDP-43 and Syn complex structure. The interfacial mutations in protein complexes changes the stability and binding affinity of the protein that may cause diseases. Here, we have utilized the computational saturation mutagenesis approach including structure-based stability and binding energy calculations to compute the systemic effects of missense mutations of TDP-43 CTD and α-Syn on protein stability and binding affinity. Most of the interfacial mutations of CTD and α-Syn were found to destabilize the protein and reduced the protein binding affinity. The results thus shed light on the functional consequences of missense mutations observed in TDP-43 associated proteinopathies and may provide the mechanisms of co-morbidities involving these two proteins.Communicated by Ramaswamy H. Sarma.
Collapse
Affiliation(s)
- Sultan H Alamri
- Department of Family Medicine, Faculty of Medicine, King Abdulaziz University, Jeddah, Saudi Arabia
| | - Shafiul Haque
- Research and Scientific Studies Unit, College of Nursing and Allied Health Sciences, Jazan University, Jazan, Saudi Arabia
- Gilbert and Rose-Marie Chagoury School of Medicine, Lebanese American University, Beirut, Lebanon
- Centre of Medical and Bio-Allied Health Sciences Research, Ajman University, Ajman, United Arab Emirates
| | - Badra S Alghamdi
- Department of Physiology, Neuroscience Unit, Faculty of Medicine, King Abdulaziz University, Jeddah, Saudi Arabia
- Pre-Clinical Research Unit, King Fahd Medical Research Center, King Abdulaziz University, Jeddah, Saudi Arabia
| | - Haythum O Tayeb
- The Mind and Brain Studies Initiative, Neuroscience Research Unit, Department of Neurology, Faculty of Medicine, King Abdulaziz University, Jeddah, Saudi Arabia
| | - Shereen Azhari
- Department of Biological Sciences, King Abdulaziz University, Jeddah, Saudi Arabia
| | - Reem M Farsi
- Department of Biological Sciences, King Abdulaziz University, Jeddah, Saudi Arabia
| | - Abear Elmokadem
- Department of Hematology/Pediatric Oncology, King Abdulaziz University Hospital, Faculty of Medicine, King Abdulaziz University, Jeddah, Saudi Arabia
| | - Turki A Alamri
- Family and Community Medicine Department, Faculty of Medicine in Rabigh, King Abdulaziz University, Jeddah, Saudi Arabia
| | - Steve Harakeh
- King Fahd Medical Research Center, King Abdulaziz University, Jeddah, Jeddah, Saudi Arabia
- Yousef Abdul Latif Jameel Scientific Chair of Prophetic Medicine Application, Faculty of Medicine, King Abdulaziz University, Jeddah, Saudi Arabia
| | - Amresh Prakash
- Amity Institute of Integrative Sciences and Health (AIISH), Amity University Haryana, Gurgaon, India
| | - Vijay Kumar
- Amity Institute of Neuropsychology & Neurosciences, Amity University, Noida, India
| |
Collapse
|
15
|
Guillien M, Mouhand A, Sagar A, Fournet A, Allemand F, Pereira GAN, Thureau A, Bernadó P, Banères JL, Sibille N. Phosphorylation motif dictates GPCR C-terminal domain conformation and arrestin interaction. Structure 2023; 31:1394-1406.e7. [PMID: 37669668 DOI: 10.1016/j.str.2023.08.011] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/21/2023] [Revised: 07/07/2023] [Accepted: 08/09/2023] [Indexed: 09/07/2023]
Abstract
Arrestin-dependent G protein-coupled receptor (GPCR) signaling pathway is regulated by the phosphorylation state of GPCR's C-terminal domain, but the molecular bases of arrestin:receptor interaction are to be further illuminated. Here we investigated the impact of phosphorylation on the conformational features of the C-terminal region from three rhodopsin-like GPCRs, the vasopressin V2 receptor (V2R), the growth hormone secretagogue or ghrelin receptor type 1a (GHSR), and the β2-adernergic receptor (β2AR). Using phosphomimetic variants, we identified pre-formed secondary structure elements, or short linear motifs (SLiMs), that undergo specific conformational transitions upon phosphorylation. Of importance, such conformational transitions appear to favor arrestin-2 binding. Hence, our results suggest a model in which the phosphorylation-dependent structuration of the GPCR C-terminal regions would modulate arrestin binding and therefore signaling outcomes in arrestin-dependent pathways.
Collapse
Affiliation(s)
- Myriam Guillien
- Centre de Biologie Structurale (CBS), CNRS, University Montpellier, Inserm, Montpellier, France
| | - Assia Mouhand
- Centre de Biologie Structurale (CBS), CNRS, University Montpellier, Inserm, Montpellier, France
| | - Amin Sagar
- Centre de Biologie Structurale (CBS), CNRS, University Montpellier, Inserm, Montpellier, France
| | - Aurélie Fournet
- Centre de Biologie Structurale (CBS), CNRS, University Montpellier, Inserm, Montpellier, France
| | - Frédéric Allemand
- Centre de Biologie Structurale (CBS), CNRS, University Montpellier, Inserm, Montpellier, France
| | - Glaécia A N Pereira
- Institut des Biomolécules Max Mousseron (IBMM), UMR-5247, University Montpellier, CNRS, ENSCM, Montpellier, France
| | - Aurélien Thureau
- HélioBio Section, Synchrotron SOLEIL, L'Orme des Merisiers, Saint-Aubin BP 48, 91190 Gif-sur-Yvette, France
| | - Pau Bernadó
- Centre de Biologie Structurale (CBS), CNRS, University Montpellier, Inserm, Montpellier, France
| | - Jean-Louis Banères
- Institut des Biomolécules Max Mousseron (IBMM), UMR-5247, University Montpellier, CNRS, ENSCM, Montpellier, France
| | - Nathalie Sibille
- Centre de Biologie Structurale (CBS), CNRS, University Montpellier, Inserm, Montpellier, France.
| |
Collapse
|
16
|
Busch J, Paschek D. OrthoBoXY: A Simple Way to Compute True Self-Diffusion Coefficients from MD Simulations with Periodic Boundary Conditions without Prior Knowledge of the Viscosity. J Phys Chem B 2023; 127:7983-7987. [PMID: 37683293 DOI: 10.1021/acs.jpcb.3c04492] [Citation(s) in RCA: 4] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 09/10/2023]
Abstract
Recently, an analytical expression for the system size dependence and direction-dependence of self-diffusion coefficients for neat liquids due to hydrodynamic interactions has been derived for molecular dynamics (MD) simulations using orthorhombic unit cells. Based on this description, we show that for systems with a "magic" box length ratio of Lz/Lx = Lz/Ly = 2.7933596497 the computed self-diffusion coefficients Dx and Dy in the x- and y-direction become system-size independent and represent the true self-diffusion coefficient D0 = (Dx + Dy)/2. Moreover, by using this particular box geometry, the viscosity can be determined with a reasonable degree of accuracy from the difference of components of the diffusion coefficients in x-, y-, and z-directions using the simple expression η = kBT × 8.1711245653/[3πLz(Dx + Dy - 2Dz)], where kB denotes Boltzmann's constant and T represents the temperature. MD simulations of TIP4P/2005 water for various system sizes using both orthorhombic and cubic box geometries are used to test the approach.
Collapse
Affiliation(s)
- Johanna Busch
- Institut für Chemie, Abteilung Physikalische und Theoretische Chemie, Universität Rostock, Albert-Einstein-Str. 27, D-18059 Rostock, Germany
| | - Dietmar Paschek
- Institut für Chemie, Abteilung Physikalische und Theoretische Chemie, Universität Rostock, Albert-Einstein-Str. 27, D-18059 Rostock, Germany
| |
Collapse
|
17
|
Karwounopoulos J, Kaupang Å, Wieder M, Boresch S. Calculations of Absolute Solvation Free Energies with Transformato─Application to the FreeSolv Database Using the CGenFF Force Field. J Chem Theory Comput 2023; 19:5988-5998. [PMID: 37616333 PMCID: PMC10500982 DOI: 10.1021/acs.jctc.3c00691] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/23/2023] [Indexed: 08/26/2023]
Abstract
We recently introduced transformato, an open-source Python package for the automated setup of large-scale calculations of relative solvation and binding free energy differences. Here, we extend the capabilities of transformato to the calculation of absolute solvation free energy differences. After careful validation against the literature results and reference calculations with the PERT module of CHARMM, we used transformato to compute absolute solvation free energies for most molecules in the FreeSolv database (621 out of 642). The force field parameters were obtained with the program cgenff (v2.5.1), which derives missing parameters from the CHARMM general force field (CGenFF v4.6). A long-range correction for the Lennard-Jones interactions was added to all computed solvation free energies. The mean absolute error compared to the experimental data is 1.12 kcal/mol. Our results allow a detailed comparison between the AMBER and CHARMM general force fields and provide a more in-depth understanding of the capabilities and limitations of the CGenFF small molecule parameters.
Collapse
Affiliation(s)
- Johannes Karwounopoulos
- Faculty
of Chemistry, Institute of Computational Biological Chemistry, University of Vienna, Währingerstr. 17, 1090 Vienna, Austria
- Vienna
Doctoral School of Chemistry (DoSChem), University of Vienna, Währingerstr. 42, 1090 Vienna, Austria
| | - Åsmund Kaupang
- Department
of Pharmacy, Section for Pharmaceutical Chemistry, University of Oslo, 0316 Oslo, Norway
| | - Marcus Wieder
- Department
of Pharmaceutical Sciences, Pharmaceutical Chemistry Division, University of Vienna, Althanstrasse 14, 1090 Vienna, Austria
| | - Stefan Boresch
- Faculty
of Chemistry, Institute of Computational Biological Chemistry, University of Vienna, Währingerstr. 17, 1090 Vienna, Austria
| |
Collapse
|
18
|
Zhang I, Rufa DA, Pulido I, Henry MM, Rosen LE, Hauser K, Singh S, Chodera JD. Identifying and overcoming the sampling challenges in relative binding free energy calculations of a model protein:protein complex. BIORXIV : THE PREPRINT SERVER FOR BIOLOGY 2023:2023.03.07.530278. [PMID: 36945557 PMCID: PMC10028896 DOI: 10.1101/2023.03.07.530278] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Indexed: 06/18/2023]
Abstract
Relative alchemical binding free energy calculations are routinely used in drug discovery projects to optimize the affinity of small molecules for their drug targets. Alchemical methods can also be used to estimate the impact of amino acid mutations on protein:protein binding affinities, but these calculations can involve sampling challenges due to the complex networks of protein and water interactions frequently present in protein:protein interfaces. We investigate these challenges by extending a GPU-accelerated open-source relative free energy calculation package (Perses) to predict the impact of amino acid mutations on protein:protein binding. Using the well-characterized model system barnase:barstar, we describe analyses for identifying and characterizing sampling problems in protein:protein relative free energy calculations. We find that mutations with sampling problems often involve charge-changes, and inadequate sampling can be attributed to slow degrees of freedom that are mutation-specific. We also explore the accuracy and efficiency of current state-of-the-art approaches-alchemical replica exchange and alchemical replica exchange with solute tempering-for overcoming relevant sampling problems. By employing sufficiently long simulations, we achieve accurate predictions (RMSE 1.61, 95% CI: [1.12, 2.11] kcal/mol), with 86% of estimates within 1 kcal/mol of the experimentally-determined relative binding free energies and 100% of predictions correctly classifying the sign of the changes in binding free energies. Ultimately, we provide a model workflow for applying protein mutation free energy calculations to protein:protein complexes, and importantly, catalog the sampling challenges associated with these types of alchemical transformations. Our free open-source package (Perses) is based on OpenMM and available at https://github.com/choderalab/perses .
Collapse
Affiliation(s)
- Ivy Zhang
- Computational and Systems Biology Program, Sloan Kettering Institute, Memorial Sloan Kettering Cancer Center, New York, NY 10065
- Tri-Institutional PhD Program in Computational Biology and Medicine, Weill Cornell Medical College, Cornell University, New York, NY 10065
| | - Dominic A. Rufa
- Computational and Systems Biology Program, Sloan Kettering Institute, Memorial Sloan Kettering Cancer Center, New York, NY 10065
- Tri-Institutional PhD Program in Chemical Biology, Weill Cornell Medical College, Cornell University, New York, NY 10065
| | - Iván Pulido
- Computational and Systems Biology Program, Sloan Kettering Institute, Memorial Sloan Kettering Cancer Center, New York, NY 10065
| | - Michael M. Henry
- Computational and Systems Biology Program, Sloan Kettering Institute, Memorial Sloan Kettering Cancer Center, New York, NY 10065
| | | | | | - Sukrit Singh
- Computational and Systems Biology Program, Sloan Kettering Institute, Memorial Sloan Kettering Cancer Center, New York, NY 10065
| | - John D. Chodera
- Computational and Systems Biology Program, Sloan Kettering Institute, Memorial Sloan Kettering Cancer Center, New York, NY 10065
| |
Collapse
|
19
|
Liu Q, Chen M, Sun L, Liu G, Xu R. Pore density effect on separations of water/ethanol and methanol/ethanol through graphene oxide membranes: A theoretical study. Sep Purif Technol 2023. [DOI: 10.1016/j.seppur.2022.122975] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/24/2022]
|
20
|
Vasudevan K, Udhaya Kumar S, Mithun A, Raghavendra B, George Priya Doss C. Structure-based virtual screening to identify potential lipase inhibitors to reduce lipid storage in Wolman disorder. ADVANCES IN PROTEIN CHEMISTRY AND STRUCTURAL BIOLOGY 2023; 133:351-363. [PMID: 36707205 DOI: 10.1016/bs.apcsb.2022.10.010] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/03/2022]
Abstract
Wolman disorder (WD) was first described in Iranian-Jewish (IJ) children, and it is caused by a deficiency of the lysosomal acid lipase (LAL). Newborns with WD are healthy and active at birth but soon develop severe malnutrition symptoms and often die before 1 year. In particular, spleens, livers, bone marrows, intestines, adrenal glands, and lymph nodes accumulate harmful amounts of lipids. G87V mutation in LIPA is responsible for Wolman disorder. Some reports suggest that δ-tocopherol can reduce lipid accumulation in cholesterol storage disorders. Hence, we used δ-tocopherol for the virtual screening process in this study. Initially, the lead compounds were docked with native and G87V mutant LIPA. Subsequently, the ADME and toxicity parameters for screened compounds were determined to ensure the safety profiles. Finally, the molecular dynamics simulations result indicated that dl-alpha-Tocopherol-13C3, a molecule obtained from the PubChem database, is identified as a potential and stable lead molecule that could be effective against the G87V mutant form of LIPA.
Collapse
Affiliation(s)
- Karthick Vasudevan
- Department of Biotechnology, School of Applied Sciences, REVA University, Bengaluru, Karnataka, India
| | - S Udhaya Kumar
- Laboratory of Integrative Genomics, Department of Integrative Biology, School of Bio Sciences and Technology, Vellore Institute of Technology, Vellore, Tamil Nadu, India
| | - A Mithun
- Department of Biotechnology, School of Applied Sciences, REVA University, Bengaluru, Karnataka, India
| | - B Raghavendra
- Department of Biotechnology, School of Applied Sciences, REVA University, Bengaluru, Karnataka, India
| | - C George Priya Doss
- Laboratory of Integrative Genomics, Department of Integrative Biology, School of Bio Sciences and Technology, Vellore Institute of Technology, Vellore, Tamil Nadu, India.
| |
Collapse
|
21
|
Nandi R, Bhowmik D, Srivastava R, Prakash A, Kumar D. Discovering potential inhibitors against SARS-CoV-2 by targeting Nsp13 Helicase. J Biomol Struct Dyn 2022; 40:12062-12074. [PMID: 34455933 DOI: 10.1080/07391102.2021.1970024] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/24/2022]
Abstract
The rise in the incidence of COVID-19 as a result of SARS-CoV-2 infection has threatened public health globally. Till now, there have been no proper prophylactics available to fight COVID-19, necessitating the advancement and evolution of effective curative against SARS-CoV-2. This study aimed at the nonstructural protein 13 (nsp13) helicase as a promising target for drug development against COVID-19. A unique collection of nucleoside analogs was screened against the SARS-CoV-2 helicase protein, for which a molecular docking experiment was executed to depict the selected ligand's binding affinity with the SARS-CoV-2 helicase proteins. Simultaneously, molecular dynamic simulations were performed to examine the protein's binding site's conformational stability, flexibility, and interaction with the ligands. Key nucleoside ligands were selected for pharmacokinetic analysis based on their docking scores. Selected ligands (cordycepin and pritelivir) showed excellent pharmacokinetics and were well stabilized at the proteins' binding site throughout the MD simulation. We have also performed binding free energy analysis or the binding characteristics of ligands with Nsp13 by using MM-PBSA and MM-GBSA. Free energy calculation by MM-PBSA and MM-GBSA analysis suggests that pritelivir may work as viable therapeutics for efficient drug advancement against SARS-CoV-2 Nsp13 helicase, potentially arresting the SARS-CoV-2 replication.Communicated by Ramaswamy H. Sarma.
Collapse
Affiliation(s)
- Rajat Nandi
- Department of Microbiology, Assam University, Silchar, Assam, India
| | - Deep Bhowmik
- Department of Microbiology, Assam University, Silchar, Assam, India
| | - Rakesh Srivastava
- School of Computational and Integrative Sciences, Jawaharlal Nehru University, New Delhi, Delhi, India
| | - Amresh Prakash
- Amity Institute of Integrative Sciences and Health, Amity University Haryana, Gurgaon, Haryana, India
| | - Diwakar Kumar
- Department of Microbiology, Assam University, Silchar, Assam, India
| |
Collapse
|
22
|
Liu Q, Wang X, Guo Y, Liu G, Zhou KG. Mechanism of ethanol/water reverse separation through a functional graphene membrane: a molecular simulation investigation. Front Chem Sci Eng 2022. [DOI: 10.1007/s11705-022-2246-z] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/23/2022]
|
23
|
Eistrikh-Heller PA, Rubinsky SV, Samygina VR, Lashkov AA. Calculation of Free Energy of Binding for Widely Specific Pyrimidine-Nucleoside Phosphorylase and Suspected Inhibitors. RUSSIAN JOURNAL OF BIOORGANIC CHEMISTRY 2022. [DOI: 10.1134/s1068162022060103] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/25/2022]
|
24
|
Zhang R, Xu Y, Lan J, Fan S, Huang J, Xu F. Structural Achievability of an NH-π Interaction between Gln and Phe in a Crystal Structure of a Collagen-like Peptide. Biomolecules 2022; 12:1433. [PMID: 36291642 PMCID: PMC9599227 DOI: 10.3390/biom12101433] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/30/2022] [Revised: 09/22/2022] [Accepted: 09/29/2022] [Indexed: 11/17/2022] Open
Abstract
NH-π interactions between polar and aromatic residues are well distributed in proteins whose stabilizing effects have been investigated in globular and fibrous proteins. In order to gain structural insights into side chain NH-π interactions, we solved a crystal structure of a collagen-like peptide containing Gln-Phe pairs. The Gln-Phe NH-π interactions were further characterized by quantum calculations, molecular simulations, and structural bioinformatics. The analyses indicated that the NH-π interactions are robust under various solvent conditions, can be distributed either on the protein surface or in its hydrophobic core and can form at a wide range of distances between residues. This study suggested that NH-π interactions can play a versatile role in protein design, including engineering hydrophobic cores, solvent accessible surfaces, and protein-protein interfaces.
Collapse
Affiliation(s)
- Ruixue Zhang
- Ministry of Education Key Laboratory of Industrial Biotechnology, School of Biotechnology, Jiangnan University, Wuxi 214122, China
| | - You Xu
- Key Laboratory of Structural Biology of Zhejiang Province, School of Life Sciences, Westlake University, Hangzhou 310024, China
- Westlake AI Therapeutics Lab, Westlake Laboratory of Life Sciences and Biomedicine, Hangzhou 310024, China
| | - Jun Lan
- Center for Structural Biology, School of Life Sciences, Tsinghua University, Beijing 100084, China
| | - Shilong Fan
- Center for Structural Biology, School of Life Sciences, Tsinghua University, Beijing 100084, China
| | - Jing Huang
- Key Laboratory of Structural Biology of Zhejiang Province, School of Life Sciences, Westlake University, Hangzhou 310024, China
- Westlake AI Therapeutics Lab, Westlake Laboratory of Life Sciences and Biomedicine, Hangzhou 310024, China
| | - Fei Xu
- Ministry of Education Key Laboratory of Industrial Biotechnology, School of Biotechnology, Jiangnan University, Wuxi 214122, China
| |
Collapse
|
25
|
Oliveira MP, Gonçalves YMH, Ol Gheta SK, Rieder SR, Horta BAC, Hünenberger PH. Comparison of the United- and All-Atom Representations of (Halo)alkanes Based on Two Condensed-Phase Force Fields Optimized against the Same Experimental Data Set. J Chem Theory Comput 2022; 18:6757-6778. [PMID: 36190354 DOI: 10.1021/acs.jctc.2c00524] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
The level of accuracy that can be achieved by a force field is influenced by choices made in the interaction-function representation and in the relevant simulation parameters. These choices, referred to here as functional-form variants (FFVs), include for example the model resolution, the charge-derivation procedure, the van der Waals combination rules, the cutoff distance, and the treatment of the long-range interactions. Ideally, assessing the effect of a given FFV on the intrinsic accuracy of the force-field representation requires that only the specific FFV is changed and that this change is performed at an optimal level of parametrization, a requirement that may prove extremely challenging to achieve in practice. Here, we present a first attempt at such a comparison for one specific FFV, namely the choice of a united-atom (UA) versus an all-atom (AA) resolution in a force field for saturated acyclic (halo)alkanes. Two force-field versions (UA vs AA) are optimized in an automated way using the CombiFF approach against 961 experimental values for the pure-liquid densities ρliq and vaporization enthalpies ΔHvap of 591 compounds. For the AA force field, the torsional and third-neighbor Lennard-Jones parameters are also refined based on quantum-mechanical rotational-energy profiles. The comparison between the UA and AA resolutions is also extended to properties that have not been included as parameterization targets, namely the surface-tension coefficient γ, the isothermal compressibility κT, the isobaric thermal-expansion coefficient αP, the isobaric heat capacity cP, the static relative dielectric permittivity ϵ, the self-diffusion coefficient D, the shear viscosity η, the hydration free energy ΔGwat, and the free energy of solvation ΔGche in cyclohexane. For the target properties ρliq and ΔHvap, the UA and AA resolutions reach very similar levels of accuracy after optimization. For the nine other properties, the AA representation leads to more accurate results in terms of η; comparably accurate results in terms of γ, κT, αP, ϵ, D, and ΔGche; and less accurate results in terms of cP and ΔGwat. This work also represents a first step toward the calibration of a GROMOS-compatible force field at the AA resolution.
Collapse
Affiliation(s)
- Marina P Oliveira
- Laboratorium für Physikalische Chemie, ETH Zürich, ETH-Hönggerberg, HCI, CH-8093 Zürich, Switzerland
| | - Yan M H Gonçalves
- Laboratorium für Physikalische Chemie, ETH Zürich, ETH-Hönggerberg, HCI, CH-8093 Zürich, Switzerland
| | - S Kashef Ol Gheta
- Laboratorium für Physikalische Chemie, ETH Zürich, ETH-Hönggerberg, HCI, CH-8093 Zürich, Switzerland
| | - Salomé R Rieder
- Laboratorium für Physikalische Chemie, ETH Zürich, ETH-Hönggerberg, HCI, CH-8093 Zürich, Switzerland
| | - Bruno A C Horta
- Laboratorium für Physikalische Chemie, ETH Zürich, ETH-Hönggerberg, HCI, CH-8093 Zürich, Switzerland
| | - Philippe H Hünenberger
- Laboratorium für Physikalische Chemie, ETH Zürich, ETH-Hönggerberg, HCI, CH-8093 Zürich, Switzerland
| |
Collapse
|
26
|
Yu H, Yang S, Chen Z, Xu Z, Quan X, Zhou J. Orientation and Conformation of Hydrophobin at the Oil-Water Interface: Insights from Molecular Dynamics Simulations. LANGMUIR : THE ACS JOURNAL OF SURFACES AND COLLOIDS 2022; 38:6191-6200. [PMID: 35508911 DOI: 10.1021/acs.langmuir.2c00614] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/14/2023]
Abstract
Hydrophobins, a new class of potential protein emulsifiers, have been extensively employed in the food, pharmaceutical, and chemical industries. However, the knowledge of the underlying molecular mechanism of protein adsorption at the oil-water interface remains elusive. In this study, all-atom molecular dynamics simulations were performed to probe the adsorption orientation and conformation change of class II hydrophobin HFBI at the cyclohexane-water interface. It was proposed that a hydrophobic dipole of the protein could be used to quantitatively predict the orientation of the adsorbed HFBI. Simulation results revealed that HFBI adsorbed at the interface with the patch-up orientation toward the oil phase, regardless of its initial orientations. HFBI's secondary structure was maintained to be intact in the course of simulations despite relatively significant variations in the tertiary structure observed, which could well preserve the bioactivity of HFBI. From the energy analysis, the driving force for interface adsorption was primarily determined by van der Waals interactions between HFBI and cyclohexane. Further analysis indicated that the adsorption orientation and conformation of HFBI at the oil-water interface were typically regulated by the hydrophobic patch and some key residues. This study provides some insights into the orientation, conformation, and adsorption mechanism of proteins at the oil-water interface and theoretical guidelines for the design and development of novel biological emulsifiers involved in the food, pharmaceutical, and chemical industries.
Collapse
Affiliation(s)
- Hai Yu
- School of Chemistry and Chemical Engineering, Guangdong Provincial Key Lab for Green Chemical Product Technology, South China University of Technology, Guangzhou 510640, P. R. China
| | - Shengjiang Yang
- School of Chemistry and Chemical Engineering, Guangdong Provincial Key Lab for Green Chemical Product Technology, South China University of Technology, Guangzhou 510640, P. R. China
| | - Zheng Chen
- School of Chemistry and Chemical Engineering, Guangdong Provincial Key Lab for Green Chemical Product Technology, South China University of Technology, Guangzhou 510640, P. R. China
| | - Zhiyong Xu
- School of Chemistry and Chemical Engineering, Guangdong Provincial Key Lab for Green Chemical Product Technology, South China University of Technology, Guangzhou 510640, P. R. China
| | - Xuebo Quan
- School of Chemistry and Chemical Engineering, Guangdong Provincial Key Lab for Green Chemical Product Technology, South China University of Technology, Guangzhou 510640, P. R. China
| | - Jian Zhou
- School of Chemistry and Chemical Engineering, Guangdong Provincial Key Lab for Green Chemical Product Technology, South China University of Technology, Guangzhou 510640, P. R. China
| |
Collapse
|
27
|
Li R, Singh R, Kashav T, Yang C, Sharma RD, Lynn AM, Prasad R, Prakash A, Kumar V. Computational Insights of Unfolding of N-Terminal Domain of TDP-43 Reveal the Conformational Heterogeneity in the Unfolding Pathway. Front Mol Neurosci 2022; 15:822863. [PMID: 35548668 PMCID: PMC9083116 DOI: 10.3389/fnmol.2022.822863] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/26/2021] [Accepted: 03/22/2022] [Indexed: 02/05/2023] Open
Abstract
TDP-43 proteinopathies is a disease hallmark that characterizes amyotrophic lateral sclerosis (ALS) and frontotemporal lobar degeneration (FTLD). The N-terminal domain of TDP-43 (NTD) is important to both TDP-43 physiology and TDP-43 proteinopathy. However, its folding and dimerization process is still poorly characterized. In the present study, we have investigated the folding/unfolding of NTD employing all-atom molecular dynamics (MD) simulations in 8 M dimethylsulfoxide (DMSO) at high temperatures. The MD results showed that the unfolding of the NTD at high temperature evolves through the formation of a number of conformational states differing in their stability and free energy. The presence of structurally heterogeneous population of intermediate ensembles was further characterized by the different extents of solvent exposure of Trp80 during unfolding. We suggest that these non-natives unfolded intermediate ensembles may facilitate NTD oligomerization and subsequently TDP-43 oligomerization, which might lead to the formation of irreversible pathological aggregates, characteristics of disease pathogenesis.
Collapse
Affiliation(s)
- Ruiting Li
- School of Engineering, Guangzhou College of Technology and Business, Guangzhou, China
| | - Ruhar Singh
- School of Computational & Integrative Sciences, Jawaharlal Nehru University, New Delhi, India
| | - Tara Kashav
- Department of Life Science, Central University of South Bihar, Gaya, India
| | - ChunMin Yang
- School of Engineering, Guangzhou College of Technology and Business, Guangzhou, China
| | - Ravi Datta Sharma
- Amity Institute of Biotechnology, Amity University Haryana, Gurgaon, India
| | - Andrew M. Lynn
- School of Computational & Integrative Sciences, Jawaharlal Nehru University, New Delhi, India
| | - Rajendra Prasad
- Amity Institute of Biotechnology, Amity University Haryana, Gurgaon, India
| | - Amresh Prakash
- Amity Institute of Integrative Sciences and Health (AIISH), Amity University Haryana, Gurgaon, India
- *Correspondence: Vijay Kumar Amresh Prakash
| | - Vijay Kumar
- Amity Institute of Neuropsychology & Neurosciences (AINN), Amity University, Noida, India
- *Correspondence: Vijay Kumar Amresh Prakash
| |
Collapse
|
28
|
Rizzuti B. Molecular simulations of proteins: From simplified physical interactions to complex biological phenomena. BIOCHIMICA ET BIOPHYSICA ACTA. PROTEINS AND PROTEOMICS 2022; 1870:140757. [PMID: 35051666 DOI: 10.1016/j.bbapap.2022.140757] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 10/27/2021] [Revised: 01/08/2022] [Accepted: 01/10/2022] [Indexed: 12/22/2022]
Abstract
Molecular dynamics simulation is the most popular computational technique for investigating the structural and dynamical behaviour of proteins, in search of the molecular basis of their function. Far from being a completely settled field of research, simulations are still evolving to best capture the essential features of the atomic interactions that govern a protein's inner motions. Modern force fields are becoming increasingly accurate in providing a physical description adequate to this purpose, and allow us to model complex biological systems under fairly realistic conditions. Furthermore, the use of accelerated sampling techniques is improving our access to the observation of progressively larger molecular structures, longer time scales, and more hidden functional events. In this review, the basic principles of molecular dynamics simulations and a number of key applications in the area of protein science are summarized, and some of the most important results are discussed. Examples include the study of the structure, dynamics and binding properties of 'difficult' targets, such as intrinsically disordered proteins and membrane receptors, and the investigation of challenging phenomena like hydration-driven processes and protein aggregation. The findings described provide an overall picture of the current state of this research field, and indicate new perspectives on the road ahead to the upcoming future of molecular simulations.
Collapse
Affiliation(s)
- Bruno Rizzuti
- CNR-NANOTEC, SS Rende (CS), Department of Physics, University of Calabria, 87036 Rende, Italy; Institute for Biocomputation and Physics of Complex Systems (BIFI), Joint Unit GBsC-CSIC-BIFI, University of Zaragoza, 50018 Zaragoza, Spain.
| |
Collapse
|
29
|
Tempra C, Ollila OHS, Javanainen M. Accurate Simulations of Lipid Monolayers Require a Water Model with Correct Surface Tension. J Chem Theory Comput 2022; 18:1862-1869. [PMID: 35133839 PMCID: PMC8908734 DOI: 10.1021/acs.jctc.1c00951] [Citation(s) in RCA: 15] [Impact Index Per Article: 7.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/28/2022]
Abstract
Lipid monolayers provide our lungs and eyes their functionality and serve as proxy systems in biomembrane research. Therefore, lipid monolayers have been studied intensively including using molecular dynamics simulations, which are able to probe their lateral structure and interactions with, e.g., pharmaceuticals or nanoparticles. However, such simulations have struggled in describing the forces at the air-water interface. Particularly, the surface tension of water and long-range van der Waals interactions have been considered critical, but their importance in monolayer simulations has been evaluated only separately. Here, we combine the recent C36/LJ-PME lipid force field that includes long-range van der Waals forces with water models that reproduce experimental surface tensions to elucidate the importance of these contributions in monolayer simulations. Our results suggest that a water model with correct surface tension is necessary to reproduce experimental surface pressure-area isotherms and monolayer phase behavior. The latter includes the liquid expanded and liquid condensed phases, their coexistence, and the opening of pores at the correct area per lipid upon expansion. Despite these improvements of the C36/LJ-PME with certain water models, the standard cutoff-based CHARMM36 lipid model with the 4-point OPC water model still provides the best agreement with experiments. Our results emphasize the importance of using high-quality water models in applications and parameter development in molecular dynamics simulations of biomolecules.
Collapse
Affiliation(s)
- Carmelo Tempra
- Institute of Organic Chemistry and Biochemistry, Czech Academy of Sciences, Flemingovo nám. 542/2, 160 00 Prague 6, Czech Republic
| | - O H Samuli Ollila
- Institute of Biotechnology, University of Helsinki, 00790 Helsinki, Finland
| | - Matti Javanainen
- Institute of Organic Chemistry and Biochemistry, Czech Academy of Sciences, Flemingovo nám. 542/2, 160 00 Prague 6, Czech Republic.,Institute of Biotechnology, University of Helsinki, 00790 Helsinki, Finland
| |
Collapse
|
30
|
Efficient separation of (C1–C2) alcohol solutions by graphyne membranes: A molecular simulation study. J Memb Sci 2022. [DOI: 10.1016/j.memsci.2021.120139] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/17/2022]
|
31
|
Hsieh MK, Yu Y, Klauda JB. All-Atom Modeling of Complex Cellular Membranes. LANGMUIR : THE ACS JOURNAL OF SURFACES AND COLLOIDS 2022; 38:3-17. [PMID: 34962814 DOI: 10.1021/acs.langmuir.1c02084] [Citation(s) in RCA: 7] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/14/2023]
Abstract
Cell membranes are composed of a variety of lipids and proteins where they interact with each other to fulfill their roles. The first step in modeling these interactions in molecular simulations is to have reliable mimetics of the membrane's lipid environment. This Feature Article presents our recent efforts to model complex cellular membranes using all-atom force fields. A short review of the CHARMM36 (C36) lipid force field and its recent update to incorporate the long-range dispersion is presented. Key examples of model membranes mimicking various species and organelles are given. These include single-celled organisms such as bacteria (E. coli., chlamydia, and P. aeruginosa) and yeast (plasma membrane, endoplasmic reticulum, and trans-Golgi network) and more advanced ones such as plants (soybean and Arabidopsis thaliana) and mammals (ocular lens, stratum corneum, and peripheral nerve myelin). Leaflet asymmetry in composition has also been applied to some of these models. With the increased lipid diversity in the C36 lipid FF, these complex models can better reflect the structural, mechanical, and dynamic properties of realistic membranes and open an opportunity to study biological processes involving other molecules.
Collapse
|
32
|
Yao J, Chen C, Zhang J, Zhang L, Zhang W, Shen JW, Liang L. Molecular understanding of charge effect on desalination performance in lamellar MoS 2 membranes. Phys Chem Chem Phys 2022; 24:26879-26889. [DOI: 10.1039/d2cp02145e] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/06/2022]
Abstract
The effect of atomic charge information on the desalination performance of lamellar MoS2 membranes was investigated at the molecular level.
Collapse
Affiliation(s)
- Junhui Yao
- Department of Chemistry, Key Laboratory of Surface & Interface Science of Polymer Materials of Zhejiang Province. Zhejiang Sci-Tech University, Hangzhou, 310018, P. R. China
| | - Chen Chen
- Department of Chemistry, Key Laboratory of Surface & Interface Science of Polymer Materials of Zhejiang Province. Zhejiang Sci-Tech University, Hangzhou, 310018, P. R. China
| | - Jing Zhang
- Department of Chemistry, Key Laboratory of Surface & Interface Science of Polymer Materials of Zhejiang Province. Zhejiang Sci-Tech University, Hangzhou, 310018, P. R. China
| | - Li Zhang
- Department of Chemistry, Key Laboratory of Surface & Interface Science of Polymer Materials of Zhejiang Province. Zhejiang Sci-Tech University, Hangzhou, 310018, P. R. China
| | - Wei Zhang
- Department of Chemistry, Key Laboratory of Surface & Interface Science of Polymer Materials of Zhejiang Province. Zhejiang Sci-Tech University, Hangzhou, 310018, P. R. China
| | - Jia-Wei Shen
- School of Pharmacy, Hangzhou Normal University, Hangzhou, 310016, P. R. China
| | - Lijun Liang
- College of Life Information Science and Instrument Engineering, Hangzhou Dianzi University, Hangzhou, 310018, P. R. China
| |
Collapse
|
33
|
Carlson S, Becker M, Brünig FN, Ataka K, Cruz R, Yu L, Tang P, Kanduč M, Haag R, Heberle J, Makki H, Netz RR. Hydrophobicity of Self-Assembled Monolayers of Alkanes: Fluorination, Density, Roughness, and Lennard-Jones Cutoffs. LANGMUIR : THE ACS JOURNAL OF SURFACES AND COLLOIDS 2021; 37:13846-13858. [PMID: 34787431 DOI: 10.1021/acs.langmuir.1c02187] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/13/2023]
Abstract
The interplay of fluorination and structure of alkane self-assembled monolayers and how these affect hydrophobicity are explored via molecular dynamics simulations, contact angle goniometry, and surface-enhanced infrared absorption spectroscopy. Wetting coefficients are found to grow linearly in the monolayer density for both alkane and perfluoroalkane monolayers. The larger contact angles of monolayers of perfluorinated alkanes are shown to be primarily caused by their larger molecular volume, which leads to a larger nearest-neighbor grafting distance and smaller tilt angle. Increasing the Lennard-Jones force cutoff in simulations is found to increase hydrophilicity. Specifically, wetting coefficients scale like the inverse square of the cutoff, and when extrapolated to the infinite cutoff limit, they yield contact angles that compare favorably to experimental values. Nanoscale roughness is also found to reliably increase monolayer hydrophobicity, mostly via the reduction of the entropic part of the work of adhesion. Analysis of depletion lengths shows that droplets on nanorough surfaces partially penetrate the surface, intermediate between Wenzel and Cassie-Baxter states.
Collapse
Affiliation(s)
- Shane Carlson
- Fachbereich Physik, Freie Universität Berlin, Arnimallee 14, 14195 Berlin, Germany
| | - Maximilian Becker
- Fachbereich Physik, Freie Universität Berlin, Arnimallee 14, 14195 Berlin, Germany
| | - Florian N Brünig
- Fachbereich Physik, Freie Universität Berlin, Arnimallee 14, 14195 Berlin, Germany
| | - Kenichi Ataka
- Fachbereich Physik, Freie Universität Berlin, Arnimallee 14, 14195 Berlin, Germany
| | - Rubén Cruz
- Fachbereich Physik, Freie Universität Berlin, Arnimallee 14, 14195 Berlin, Germany
| | - Leixiao Yu
- Institut für Chemie und Biochemie, Freie Universität Berlin, Takustraße 3, 14195 Berlin, Germany
| | - Peng Tang
- Institut für Chemie und Biochemie, Freie Universität Berlin, Takustraße 3, 14195 Berlin, Germany
| | - Matej Kanduč
- Department of Theoretical Physics, Jožef Stefan Institute, Jamova cesta 39, 1000 Ljubljana, Slovenia
| | - Rainer Haag
- Institut für Chemie und Biochemie, Freie Universität Berlin, Takustraße 3, 14195 Berlin, Germany
| | - Joachim Heberle
- Fachbereich Physik, Freie Universität Berlin, Arnimallee 14, 14195 Berlin, Germany
| | - Hesam Makki
- Polymer and Color Engineering, Amirkabir University of Technology, 424 Hafez Ave, Tehran 15875-4413, Iran
| | - Roland R Netz
- Fachbereich Physik, Freie Universität Berlin, Arnimallee 14, 14195 Berlin, Germany
| |
Collapse
|
34
|
Xu Y, Huang J. Validating the CHARMM36m protein force field with LJ-PME reveals altered hydrogen bonding dynamics under elevated pressures. Commun Chem 2021; 4:99. [PMID: 36697521 PMCID: PMC9814493 DOI: 10.1038/s42004-021-00537-8] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/11/2021] [Accepted: 06/08/2021] [Indexed: 01/28/2023] Open
Abstract
The pressure-temperature phase diagram is important to our understanding of the physics of biomolecules. Compared to studies on temperature effects, studies of the pressure dependence of protein dynamic are rather limited. Molecular dynamics (MD) simulations with fine-tuned force fields (FFs) offer a powerful tool to explore the influence of thermodynamic conditions on proteins. Here we evaluate the transferability of the CHARMM36m (C36m) protein force field at varied pressures compared with NMR data using ubiquitin as a model protein. The pressure dependences of J couplings for hydrogen bonds and order parameters for internal motion are in good agreement with experiment. We demonstrate that the C36m FF combined with the Lennard-Jones particle-mesh Ewald (LJ-PME) method is suitable for simulations in a wide range of temperature and pressure. As the ubiquitin remains stable up to 2500 bar, we identify the mobility and stability of different hydrogen bonds in response to pressure. Based on those results, C36m is expected to be applied to more proteins in the future to further investigate protein dynamics under elevated pressures.
Collapse
Affiliation(s)
- You Xu
- grid.494629.40000 0004 8008 9315Key Laboratory of Structural Biology of Zhejiang Province, School of Life Sciences, Westlake University, Hangzhou, Zhejiang China ,Westlake AI Therapeutics Lab, Westlake Laboratory of Life Sciences and Biomedicine, Hangzhou, Zhejiang China ,grid.494629.40000 0004 8008 9315Institute of Biology, Westlake Institute for Advanced Study, Hangzhou, Zhejiang China
| | - Jing Huang
- grid.494629.40000 0004 8008 9315Key Laboratory of Structural Biology of Zhejiang Province, School of Life Sciences, Westlake University, Hangzhou, Zhejiang China ,Westlake AI Therapeutics Lab, Westlake Laboratory of Life Sciences and Biomedicine, Hangzhou, Zhejiang China ,grid.494629.40000 0004 8008 9315Institute of Biology, Westlake Institute for Advanced Study, Hangzhou, Zhejiang China
| |
Collapse
|
35
|
Kashefolgheta S, Wang S, Acree WE, Hünenberger PH. Evaluation of nine condensed-phase force fields of the GROMOS, CHARMM, OPLS, AMBER, and OpenFF families against experimental cross-solvation free energies. Phys Chem Chem Phys 2021; 23:13055-13074. [PMID: 34105547 PMCID: PMC8207520 DOI: 10.1039/d1cp00215e] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/15/2021] [Accepted: 04/28/2021] [Indexed: 12/02/2022]
Abstract
Experimental solvation free energies are nowadays commonly included as target properties in the validation of condensed-phase force fields, sometimes even in their calibration. In a previous article [Kashefolgheta et al., J. Chem. Theory. Comput., 2020, 16, 7556-7580], we showed how the involved comparison between experimental and simulation results could be made more systematic by considering a full matrix of cross-solvation free energies . For a set of N molecules that are all in the liquid state under ambient conditions, such a matrix encompasses N×N entries for considering each of the N molecules either as solute (A) or as solvent (B). In the quoted study, a cross-solvation matrix of 25 × 25 experimental value was introduced, considering 25 small molecules representative for alkanes, chloroalkanes, ethers, ketones, esters, alcohols, amines, and amides. This experimental data was used to compare the relative accuracies of four popular condensed-phase force fields, namely GROMOS-2016H66, OPLS-AA, AMBER-GAFF, and CHARMM-CGenFF. In the present work, the comparison is extended to five additional force fields, namely GROMOS-54A7, GROMOS-ATB, OPLS-LBCC, AMBER-GAFF2, and OpenFF. Considering these nine force fields, the correlation coefficients between experimental values and simulation results range from 0.76 to 0.88, the root-mean-square errors (RMSEs) from 2.9 to 4.8 kJ mol-1, and average errors (AVEEs) from -1.5 to +1.0 kJ mol-1. In terms of RMSEs, GROMOS-2016H66 and OPLS-AA present the best accuracy (2.9 kJ mol-1), followed by OPLS-LBCC, AMBER-GAFF2, AMBER-GAFF, and OpenFF (3.3 to 3.6 kJ mol-1), and then by GROMOS-54A7, CHARM-CGenFF, and GROMOS-ATB (4.0 to 4.8 kJ mol-1). These differences are statistically significant but not very pronounced, and are distributed rather heterogeneously over the set of compounds within the different force fields.
Collapse
Affiliation(s)
- Sadra Kashefolgheta
- Laboratorium für Physikalische Chemie, ETH Zürich, ETH-Hönggerberg, HCICH-8093 ZürichSwitzerland+41 44 632 55 03
| | - Shuzhe Wang
- Laboratorium für Physikalische Chemie, ETH Zürich, ETH-Hönggerberg, HCICH-8093 ZürichSwitzerland+41 44 632 55 03
| | - William E. Acree
- Department of Chemistry, University of North Texas1155 Union Circle Drive #305070DentonTexas 76203USA
| | - Philippe H. Hünenberger
- Laboratorium für Physikalische Chemie, ETH Zürich, ETH-Hönggerberg, HCICH-8093 ZürichSwitzerland+41 44 632 55 03
| |
Collapse
|
36
|
Abstract
Molecular simulations of biological molecules require an accurate description of molecular interactions through a force field (FF). The focus of this Perspective is on all-atom lipid FFs. Recent additions to the CHARMM36 lipid FF continue to expand a researcher's ability to probe membrane structure and function with a wide variety of biologically important lipids. Currently, there is an effort to reduce the assumptions in all-atom lipid FFs. The inclusion of long-range dispersion interaction through particle-mesh Ewald is allowing for more accurate descriptions of lipid bilayer and monolayer properties without additional computational cost. Soon, simulations with lipid FFs will no longer depend on short-range cutoffs and will accurately represent long-range dispersion. This requires efficient FF parametrization with an automated approach due to FF complexity. In addition, polarizable FFs for lipids will be important for the next generation of simulations that accurately represent how molecule interactions respond to a varied environment.
Collapse
|
37
|
Busch J, Neumann J, Paschek D. An exact a posteriori correction for hydrogen bond population correlation functions and other reversible geminate recombinations obtained from simulations with periodic boundary conditions. Liquid water as a test case. J Chem Phys 2021; 154:214501. [PMID: 34240960 DOI: 10.1063/5.0053445] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
The kinetics of breaking and re-formation of hydrogen bonds (HBs) in liquid water is a prototype of reversible geminate recombination. HB population correlation functions (HBPCFs) are a means to study the HB kinetics. The long-time limiting behavior of HBPCFs is controlled by translatoric diffusion and shows a t-3/2 time-dependence, which can be described by analytical expressions based on the HB acceptor density and the donor-acceptor inter-diffusion coefficient. If the trajectories are not properly "unwrapped," the presence of periodic boundary conditions (PBCs) can perturb this long-time limiting behavior. Keeping the trajectories "wrapped," however, allows for a more efficient calculation of HBPCFs. We discuss the consequences of PBCs in combination with "wrapped" trajectories following from the approximations according to Luzar-Chandler and according to Starr, each deviating in a different fashion from the true long-time limiting behavior, but enveloping the unperturbed function. A simple expression is given for estimating the maximum time up to which the computed HBPCFs reliably describe the long-time limiting behavior. In addition, an exact a posteriori correction for systems with PBCs for "wrapped" trajectories is derived, which can be easily computed and which is able to fully recover the true t-3/2 long-time behavior. For comparison, HBPCFs are computed from MD simulations of TIP4P/2005 model water for varying system sizes and temperatures of 273 and 298 K using this newly introduced correction. Implications for the computations of HB lifetimes and the effect of the system-size are discussed.
Collapse
Affiliation(s)
- Johanna Busch
- Institut für Chemie, Abteilung Physikalische und Theoretische Chemie, Universität Rostock, Albert-Einstein-Str. 21, D-18059 Rostock, Germany
| | - Jan Neumann
- Institut für Chemie, Abteilung Physikalische und Theoretische Chemie, Universität Rostock, Albert-Einstein-Str. 21, D-18059 Rostock, Germany
| | - Dietmar Paschek
- Institut für Chemie, Abteilung Physikalische und Theoretische Chemie, Universität Rostock, Albert-Einstein-Str. 21, D-18059 Rostock, Germany
| |
Collapse
|
38
|
Raturi S, Nair AV, Shinoda K, Singh H, Bai B, Murakami S, Fujitani H, van Veen HW. Engineered MATE multidrug transporters reveal two functionally distinct ion-coupling pathways in NorM from Vibrio cholerae. Commun Biol 2021; 4:558. [PMID: 33976372 PMCID: PMC8113278 DOI: 10.1038/s42003-021-02081-6] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/25/2020] [Accepted: 04/01/2021] [Indexed: 11/09/2022] Open
Abstract
Multidrug and toxic compound extrusion (MATE) transport proteins confer multidrug resistance on pathogenic microorganisms and affect pharmacokinetics in mammals. Our understanding of how MATE transporters work, has mostly relied on protein structures and MD simulations. However, the energetics of drug transport has not been studied in detail. Many MATE transporters utilise the electrochemical H+ or Na+ gradient to drive substrate efflux, but NorM-VC from Vibrio cholerae can utilise both forms of metabolic energy. To dissect the localisation and organisation of H+ and Na+ translocation pathways in NorM-VC we engineered chimaeric proteins in which the N-lobe of H+-coupled NorM-PS from Pseudomonas stutzeri is fused to the C-lobe of NorM-VC, and vice versa. Our findings in drug binding and transport experiments with chimaeric, mutant and wildtype transporters highlight the versatile nature of energy coupling in NorM-VC, which enables adaptation to fluctuating salinity levels in the natural habitat of V. cholerae.
Collapse
Affiliation(s)
- Sagar Raturi
- Department of Pharmacology, University of Cambridge, Cambridge, UK
- University College Dublin Clinical Research Centre, St. Vincent's University Hospital, Dublin, Ireland
| | - Asha V Nair
- Department of Pharmacology, University of Cambridge, Cambridge, UK
| | - Keiko Shinoda
- Microbial Membrane Transport Engineering, Biotechnology Research Center, The University of Tokyo, Bunkyo-ku, Tokyo, Japan
| | - Himansha Singh
- Department of Pharmacology, University of Cambridge, Cambridge, UK
| | - Boyan Bai
- Department of Pharmacology, University of Cambridge, Cambridge, UK
| | - Satoshi Murakami
- Department of Life Science and Technology, Tokyo Institute of Technology, Nagatsuta, Midori-ku, Yokohama, Japan
| | - Hideaki Fujitani
- Laboratories for Systems Biology and Medicine, Research Center for Advanced Science and Technology, The University of Tokyo, Meguro-ku, Tokyo, Japan
| | | |
Collapse
|
39
|
Ploetz EA, Karunaweera S, Bentenitis N, Chen F, Dai S, Gee MB, Jiao Y, Kang M, Kariyawasam NL, Naleem N, Weerasinghe S, Smith PE. Kirkwood-Buff-Derived Force Field for Peptides and Proteins: Philosophy and Development of KBFF20. J Chem Theory Comput 2021; 17:2964-2990. [PMID: 33878263 DOI: 10.1021/acs.jctc.1c00075] [Citation(s) in RCA: 12] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022]
Abstract
A new classical nonpolarizable force field, KBFF20, for the simulation of peptides and proteins is presented. The force field relies heavily on the use of Kirkwood-Buff theory to provide a comparison of simulated and experimental Kirkwood-Buff integrals for solutes containing the functional groups common in proteins, thus ensuring intermolecular interactions that provide a good balance between the peptide-peptide, peptide-solvent, and solvent-solvent distributions observed in solution mixtures. In this way, it differs significantly from other biomolecular force fields. Further development and testing of the intermolecular potentials are presented here. Subsequently, rotational potentials for the ϕ/ψ and χ dihedral degrees of freedom are obtained by analysis of the Protein Data Bank, followed by small modifications to provide a reasonable balance between simulated and observed α and β percentages for small peptides. This, the first of two articles, describes in detail the philosophy and development behind KBFF20.
Collapse
Affiliation(s)
- Elizabeth A Ploetz
- Department of Chemistry, Kansas State University, 213 CBC Building, 1212 Mid-Campus Drive North, Manhattan, Kansas 66506, United States
| | - Sadish Karunaweera
- Department of Chemistry, Kansas State University, 213 CBC Building, 1212 Mid-Campus Drive North, Manhattan, Kansas 66506, United States
| | - Nikolaos Bentenitis
- Department of Chemistry, Kansas State University, 213 CBC Building, 1212 Mid-Campus Drive North, Manhattan, Kansas 66506, United States
| | - Feng Chen
- Department of Chemistry, Kansas State University, 213 CBC Building, 1212 Mid-Campus Drive North, Manhattan, Kansas 66506, United States
| | - Shu Dai
- Department of Chemistry, Kansas State University, 213 CBC Building, 1212 Mid-Campus Drive North, Manhattan, Kansas 66506, United States
| | - Moon B Gee
- Department of Chemistry, Kansas State University, 213 CBC Building, 1212 Mid-Campus Drive North, Manhattan, Kansas 66506, United States
| | - Yuanfang Jiao
- Department of Chemistry, Kansas State University, 213 CBC Building, 1212 Mid-Campus Drive North, Manhattan, Kansas 66506, United States
| | - Myungshim Kang
- Department of Chemistry, Kansas State University, 213 CBC Building, 1212 Mid-Campus Drive North, Manhattan, Kansas 66506, United States
| | - Nilusha L Kariyawasam
- Department of Chemistry, Kansas State University, 213 CBC Building, 1212 Mid-Campus Drive North, Manhattan, Kansas 66506, United States
| | - Nawavi Naleem
- Department of Chemistry, Kansas State University, 213 CBC Building, 1212 Mid-Campus Drive North, Manhattan, Kansas 66506, United States
| | | | - Paul E Smith
- Department of Chemistry, Kansas State University, 213 CBC Building, 1212 Mid-Campus Drive North, Manhattan, Kansas 66506, United States
| |
Collapse
|
40
|
Zhang C, Huang J. Interactions Between Nucleosomes: From Atomistic Simulation to Polymer Model. Front Mol Biosci 2021; 8:624679. [PMID: 33912585 PMCID: PMC8072053 DOI: 10.3389/fmolb.2021.624679] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/31/2020] [Accepted: 02/09/2021] [Indexed: 11/23/2022] Open
Abstract
The organization of genomes in space and time dimension plays an important role in gene expression and regulation. Chromatin folding occurs in a dynamic, structured way that is subject to biophysical rules and biological processes. Nucleosomes are the basic unit of chromatin in living cells, and here we report on the effective interactions between two nucleosomes in physiological conditions using explicit-solvent all-atom simulations. Free energy landscapes derived from umbrella sampling simulations agree well with recent experimental and simulation results. Our simulations reveal the atomistic details of the interactions between nucleosomes in solution and can be used for constructing the coarse-grained model for chromatin in a bottom-up manner.
Collapse
Affiliation(s)
- Chengwei Zhang
- College of Life Sciences, Zhejiang University, Hangzhou, China
- Key Laboratory of Structural Biology of Zhejiang Province, School of Life Sciences, Westlake University, Hangzhou, China
- Westlake Laboratory of Life Sciences and Biomedicine, Hangzhou, China
- Institute of Biology, Westlake Institute for Advanced Study, Hangzhou, China
| | - Jing Huang
- Key Laboratory of Structural Biology of Zhejiang Province, School of Life Sciences, Westlake University, Hangzhou, China
- Westlake Laboratory of Life Sciences and Biomedicine, Hangzhou, China
- Institute of Biology, Westlake Institute for Advanced Study, Hangzhou, China
| |
Collapse
|
41
|
Yu Y, Krämer A, Venable RM, Brooks BR, Klauda JB, Pastor RW. CHARMM36 Lipid Force Field with Explicit Treatment of Long-Range Dispersion: Parametrization and Validation for Phosphatidylethanolamine, Phosphatidylglycerol, and Ether Lipids. J Chem Theory Comput 2021; 17:1581-1595. [PMID: 33620194 PMCID: PMC8130185 DOI: 10.1021/acs.jctc.0c01327] [Citation(s) in RCA: 34] [Impact Index Per Article: 11.3] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/08/2023]
Abstract
Long-range Lennard-Jones (LJ) interactions have been incorporated into the CHARMM36 (C36) lipid force field (FF) using the LJ particle-mesh Ewald (LJ-PME) method in order to remove the inconsistency of bilayer and monolayer properties arising from the exclusion of long-range dispersion [Yu, Y.; Semi-automated Optimization of the CHARMM36 Lipid Force Field to Include Explicit Treatment of Long-Range Dispersion. J. Chem. Theory Comput. 2021, 10.1021/acs.jctc.0c01326. (preceding article in this issue)]. The new FF is denoted C36/LJ-PME. While the first optimization was based on three phosphatidylcholines (PCs), this work extends the validation and parametrization to more lipids including PC, phosphatidylethanolamine (PE), phosphatidylglycerol (PG), and ether lipids. The agreement with experimental structure data is excellent for PC, PE, and ether lipids. C36/LJ-PME also compares favorably with scattering data of PG bilayers but less so with NMR deuterium order parameters of 1,2-dimyristoyl-sn-glycero-3-phospho-(1'-rac-glycerol) (DMPG) at 303.15 K, indicating a need for future optimization regarding PG-specific parameters. Frequency dependence of NMR T1 spin-lattice relaxation times is well-described by C36/LJ-PME, and the overall agreement with experiment is comparable to C36. Lipid diffusion is slower than C36 due to the added long-range dispersion causing a higher viscosity, although it is still too fast compared to experiment after correction for periodic boundary conditions. When using a 10 Å real-space cutoff, the simulation speed of C36/LJ-PME is roughly equal to C36. While more lipids will be incorporated into the FF in the future, C36/LJ-PME can be readily used for common lipids and extends the capability of the CHARMM FF by supporting monolayers and eliminating the cutoff dependence.
Collapse
Affiliation(s)
- Yalun Yu
- Biophysics Graduate Program, University of Maryland, College Park, Maryland 20742, United States
- Laboratory of Computational Biology, National Heart, Lung and Blood Institute, National Institutes of Health, Bethesda, Maryland 20892, United States
| | - Andreas Krämer
- Laboratory of Computational Biology, National Heart, Lung and Blood Institute, National Institutes of Health, Bethesda, Maryland 20892, United States
| | - Richard M Venable
- Laboratory of Computational Biology, National Heart, Lung and Blood Institute, National Institutes of Health, Bethesda, Maryland 20892, United States
| | - Bernard R Brooks
- Laboratory of Computational Biology, National Heart, Lung and Blood Institute, National Institutes of Health, Bethesda, Maryland 20892, United States
| | - Jeffery B Klauda
- Biophysics Graduate Program, University of Maryland, College Park, Maryland 20742, United States
- Department of Chemical and Biomolecular Engineering, University of Maryland, College Park, Maryland 20742, United States
| | - Richard W Pastor
- Laboratory of Computational Biology, National Heart, Lung and Blood Institute, National Institutes of Health, Bethesda, Maryland 20892, United States
| |
Collapse
|
42
|
Yu Y, Krämer A, Venable RM, Simmonett AC, MacKerell AD, Klauda JB, Pastor RW, Brooks BR. Semi-automated Optimization of the CHARMM36 Lipid Force Field to Include Explicit Treatment of Long-Range Dispersion. J Chem Theory Comput 2021; 17:1562-1580. [PMID: 33620214 PMCID: PMC8059446 DOI: 10.1021/acs.jctc.0c01326] [Citation(s) in RCA: 37] [Impact Index Per Article: 12.3] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/31/2022]
Abstract
The development of the CHARMM lipid force field (FF) can be traced back to the early 1990s with its current version denoted CHARMM36 (C36). The parametrization of C36 utilized high-level quantum mechanical data and free energy calculations of model compounds before parameters were manually adjusted to yield agreement with experimental properties of lipid bilayers. While such manual fine-tuning of FF parameters is based on intuition and trial-and-error, automated methods can identify beneficial modifications of the parameters via their sensitivities and thereby guide the optimization process. This work introduces a semi-automated approach to reparametrize the CHARMM lipid FF with consistent inclusion of long-range dispersion through the Lennard-Jones particle-mesh Ewald (LJ-PME) approach. The optimization method is based on thermodynamic reweighting with regularization with respect to the C36 set. Two independent optimizations with different topology restrictions are presented. Targets of the optimizations are primarily liquid crystalline phase properties of lipid bilayers and the compression isotherm of monolayers. Pair correlation functions between water and lipid functional groups in aqueous solution are also included to address headgroup hydration. While the physics of the reweighting strategy itself is well-understood, applying it to heterogeneous, complex anisotropic systems poses additional challenges. These were overcome through careful selection of target properties and reweighting settings allowing for the successful incorporation of the explicit treatment of long-range dispersion, and we denote the newly optimized lipid force field as C36/LJ-PME. The current implementation of the optimization protocol will facilitate the future development of the CHARMM and related lipid force fields.
Collapse
Affiliation(s)
- Yalun Yu
- Biophysics Graduate Program, University of Maryland, College Park, Maryland 20742, United States
- Laboratory of Computational Biology, National Heart, Lung and Blood Institute, National Institutes of Health, Bethesda, Maryland 20892, United States
| | - Andreas Krämer
- Laboratory of Computational Biology, National Heart, Lung and Blood Institute, National Institutes of Health, Bethesda, Maryland 20892, United States
| | - Richard M Venable
- Laboratory of Computational Biology, National Heart, Lung and Blood Institute, National Institutes of Health, Bethesda, Maryland 20892, United States
| | - Andrew C Simmonett
- Laboratory of Computational Biology, National Heart, Lung and Blood Institute, National Institutes of Health, Bethesda, Maryland 20892, United States
| | - Alexander D MacKerell
- Department of Pharmaceutical Sciences, School of Pharmacy, University of Maryland, Baltimore, Maryland 21201, United States
| | - Jeffery B Klauda
- Biophysics Graduate Program, University of Maryland, College Park, Maryland 20742, United States
- Department of Chemical and Biomolecular Engineering, University of Maryland, College Park, Maryland 20742, United States
| | - Richard W Pastor
- Laboratory of Computational Biology, National Heart, Lung and Blood Institute, National Institutes of Health, Bethesda, Maryland 20892, United States
| | - Bernard R Brooks
- Laboratory of Computational Biology, National Heart, Lung and Blood Institute, National Institutes of Health, Bethesda, Maryland 20892, United States
| |
Collapse
|
43
|
Schellhammer KS, Cuniberti G, Ortmann F. Investigating a Combined Stochastic Nucleation and Molecular Dynamics-Based Equilibration Approach for Constructing Large-Scale Polycrystalline Films. J Chem Theory Comput 2021; 17:1266-1275. [PMID: 33434021 DOI: 10.1021/acs.jctc.0c01196] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
The morphology of small-molecule organic semiconducting materials can vary from single crystals via polycrystalline films with varying grain sizes to amorphous structures, depending on the process conditions. This structural variety affects the electronic properties and, thus, the performance of organic electronic devices. A nucleation-equilibration approach is investigated, whose focus is on the construction of morphologies with controlled variations in the average grain size. Its computational requirements are low because nucleation is purely based on geometrical considerations, thus allowing the construction of model systems of experimentally relevant sizes. Its application is demonstrated for C60 and pentacene by generating single-component films that vary from amorphous to crystalline structures. It is further generalized to two-component films and applied to C60: pentacene blends as well as dilute n-doped C60 structures. When combined with electronic structure calculations in the future, the nucleation-equilibration approach can offer insights into the impact of polycrystallinity on electronic and charge-transport properties in the absence of any knowledge about the growth mechanism and for a broad set of systems.
Collapse
Affiliation(s)
- K Sebastian Schellhammer
- Institute for Materials Science and Max Bergmann Center of Biomaterials, Technische Universität Dresden, 01062 Dresden, Germany.,Dresden Center for Computational Materials Science, Technische Universität Dresden, 01062 Dresden, Germany.,Center for Advancing Electronics Dresden, Technische Universität Dresden, 01062 Dresden, Germany
| | - Gianaurelio Cuniberti
- Institute for Materials Science and Max Bergmann Center of Biomaterials, Technische Universität Dresden, 01062 Dresden, Germany.,Dresden Center for Computational Materials Science, Technische Universität Dresden, 01062 Dresden, Germany.,Center for Advancing Electronics Dresden, Technische Universität Dresden, 01062 Dresden, Germany
| | - Frank Ortmann
- Dresden Center for Computational Materials Science, Technische Universität Dresden, 01062 Dresden, Germany.,Center for Advancing Electronics Dresden, Technische Universität Dresden, 01062 Dresden, Germany.,Department of Chemistry, Technische Universität München, Lichtenbergstr. 4, 85748 Garching, Germany
| |
Collapse
|
44
|
Man VH, Wu X, He X, Xie XQ, Brooks BR, Wang J. Determination of van der Waals Parameters Using a Double Exponential Potential for Nonbonded Divalent Metal Cations in TIP3P Solvent. J Chem Theory Comput 2021; 17:1086-1097. [PMID: 33503371 DOI: 10.1021/acs.jctc.0c01267] [Citation(s) in RCA: 15] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/16/2022]
Abstract
A double exponential (DE) functional form for Lennard-Jones (LJ) interactions, proposed in our previous study, has many advantages over LJ potentials including a natural softcore characteristic for the convenience of the pathway-based free-energy calculations, fast convergence, and flexibility in use. In this work, we put the first step on the application of the DE functional form by identifying a DE potential, coined DE-TIP3P, for molecular simulations using the TIP3P water model. The developed DE-TIP3 potential was better than LJ potential in reproducing the experimental water properties. Afterward, we developed the nonbonded models of 15 divalent metal ions, which frequently appear and play vital roles in biological systems, to be consistent with the DE-TIP3P potential and TIP3P water model. Our nonbonded models were as good as the complicated nonbonded dummy cationic models by Jiang et al. and the nonbonded 12-6-4 LJ models by Li and Merz in reproducing the experimental properties of those ions. Moreover, our nonbonded models achieved a better performance than the compromise (CM) LJ models and 12-6-4 LJ models, developed by Li and Merz, in reproducing the properties of MgCl2 in aqueous solution.
Collapse
Affiliation(s)
- Viet Hoang Man
- Department of Pharmaceutical Sciences and Computational Chemical Genomics Screening Center, School of Pharmacy, University of Pittsburgh, Pittsburgh, Pennsylvania 15261, United States
| | - Xiongwu Wu
- Laboratory of Computational Biology, National Heart, Lung, and Blood Institute, National Institute of Health, Bethesda, Maryland 20892, United States
| | - Xibing He
- Department of Pharmaceutical Sciences and Computational Chemical Genomics Screening Center, School of Pharmacy, University of Pittsburgh, Pittsburgh, Pennsylvania 15261, United States
| | - Xiang-Qun Xie
- Department of Pharmaceutical Sciences and Computational Chemical Genomics Screening Center, School of Pharmacy, University of Pittsburgh, Pittsburgh, Pennsylvania 15261, United States
| | - Bernard R Brooks
- Laboratory of Computational Biology, National Heart, Lung, and Blood Institute, National Institute of Health, Bethesda, Maryland 20892, United States
| | - Junmei Wang
- Department of Pharmaceutical Sciences and Computational Chemical Genomics Screening Center, School of Pharmacy, University of Pittsburgh, Pittsburgh, Pennsylvania 15261, United States
| |
Collapse
|
45
|
Simmonett AC, Brooks BR. A compression strategy for particle mesh Ewald theory. J Chem Phys 2021; 154:054112. [PMID: 33557541 DOI: 10.1063/5.0040966] [Citation(s) in RCA: 23] [Impact Index Per Article: 7.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/20/2022] Open
Abstract
Particle Mesh Ewald (PME) has become a standard method for treating long-range electrostatics in molecular simulations. Although the method has inferior asymptotic computational complexity to its linear scaling competitors, it remains enormously popular due to its high efficiency, which stems from the use of fast Fourier transforms (FFTs). This use of FFTs provides great challenges for scaling the method up to massively parallel systems, in large part because of the need to transfer large amounts of data. In this work, we demonstrate that this data transfer volume can be greatly reduced as a natural consequence of the structure of the PME equations. We also suggest an alternative algorithm that supplants the FFT with a linear algebra approach, which further decreases communication costs at the expense of increased asymptotic computational complexity. This linear algebra based approach is demonstrated to have great potential for latency hiding by interleaving communication and computation steps of the short- and long-range electrostatic terms.
Collapse
Affiliation(s)
- Andrew C Simmonett
- Laboratory of Computational Biology, National Heart, Lung and Blood Institute, National Institutes of Health, Bethesda, Maryland 20892, USA
| | - Bernard R Brooks
- Laboratory of Computational Biology, National Heart, Lung and Blood Institute, National Institutes of Health, Bethesda, Maryland 20892, USA
| |
Collapse
|
46
|
Oliveira MP, Andrey M, Rieder SR, Kern L, Hahn DF, Riniker S, Horta BAC, Hünenberger PH. Systematic Optimization of a Fragment-Based Force Field against Experimental Pure-Liquid Properties Considering Large Compound Families: Application to Saturated Haloalkanes. J Chem Theory Comput 2020; 16:7525-7555. [DOI: 10.1021/acs.jctc.0c00683] [Citation(s) in RCA: 13] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/19/2022]
Affiliation(s)
- Marina P. Oliveira
- Laboratorium für Physikalische Chemie, ETH Zürich, ETH-Honggerberg, HCI, CH-8093 Zürich, Switzerland
| | - Maurice Andrey
- Laboratorium für Physikalische Chemie, ETH Zürich, ETH-Honggerberg, HCI, CH-8093 Zürich, Switzerland
| | - Salomé R. Rieder
- Laboratorium für Physikalische Chemie, ETH Zürich, ETH-Honggerberg, HCI, CH-8093 Zürich, Switzerland
| | - Leyla Kern
- Laboratorium für Physikalische Chemie, ETH Zürich, ETH-Honggerberg, HCI, CH-8093 Zürich, Switzerland
| | - David F. Hahn
- Laboratorium für Physikalische Chemie, ETH Zürich, ETH-Honggerberg, HCI, CH-8093 Zürich, Switzerland
| | - Sereina Riniker
- Laboratorium für Physikalische Chemie, ETH Zürich, ETH-Honggerberg, HCI, CH-8093 Zürich, Switzerland
| | - Bruno A. C. Horta
- Instituto de Química, Universidade Federal do Rio de Janeiro, Rio de Janeiro 21941-909, Brazil
| | - Philippe H. Hünenberger
- Laboratorium für Physikalische Chemie, ETH Zürich, ETH-Honggerberg, HCI, CH-8093 Zürich, Switzerland
| |
Collapse
|
47
|
Kashefolgheta S, Oliveira MP, Rieder SR, Horta BAC, Acree WE, Hünenberger PH. Evaluating Classical Force Fields against Experimental Cross-Solvation Free Energies. J Chem Theory Comput 2020; 16:7556-7580. [DOI: 10.1021/acs.jctc.0c00688] [Citation(s) in RCA: 19] [Impact Index Per Article: 4.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/21/2022]
Affiliation(s)
- Sadra Kashefolgheta
- Laboratorium für Physikalische Chemie, ETH Zürich, ETH-Hönggerberg, HCI, CH-8093 Zürich, Switzerland
| | - Marina P. Oliveira
- Laboratorium für Physikalische Chemie, ETH Zürich, ETH-Hönggerberg, HCI, CH-8093 Zürich, Switzerland
| | - Salomé R. Rieder
- Laboratorium für Physikalische Chemie, ETH Zürich, ETH-Hönggerberg, HCI, CH-8093 Zürich, Switzerland
| | - Bruno A. C. Horta
- Instituto de Química, Universidade Federal do Rio de Janeiro, Rio de Janeiro 21941-909, Brazil
| | - William E. Acree
- Department of Chemistry, University of North Texas, 1155 Union Circle Drive #305070, Denton, Texas 76203, United States
| | - Philippe H. Hünenberger
- Laboratorium für Physikalische Chemie, ETH Zürich, ETH-Hönggerberg, HCI, CH-8093 Zürich, Switzerland
| |
Collapse
|
48
|
Liu Q, Cheng L, Liu G. Enhanced Selective Hydrogen Permeation through Graphdiyne Membrane: A Theoretical Study. MEMBRANES 2020; 10:E286. [PMID: 33076414 PMCID: PMC7650590 DOI: 10.3390/membranes10100286] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 09/29/2020] [Revised: 10/09/2020] [Accepted: 10/13/2020] [Indexed: 11/16/2022]
Abstract
Graphdiyne (GDY), with uniform pores and atomic thickness, is attracting widespread attention for application in H2 separation in recent years. However, the challenge lies in the rational design of GDYs for fast and selective H2 permeation. By MD and DFT calculations, several flexible GDYs were constructed to investigate the permeation properties of four pure gas (H2, N2, CO2, and CH4) and three equimolar binary mixtures (H2/N2, H2/CO2, and H2/CH4) in this study. When the pore size is smaller than 2.1 Å, the GDYs acted as an exceptional filter for H2 with an approximately infinite H2 selectivity. Beyond the size-sieving effect, in the separation process of binary mixtures, the blocking effect arising from the strong gas-membrane interaction was proven to greatly impede H2 permeation. After understanding the mechanism, the H2 permeance of the mixtures of H2/CO2 was further increased to 2.84 × 105 GPU by reducing the blocking effect with the addition of a tiny amount of surface charges, without sacrificing the selectivity. This theoretical study provides an additional atomic understanding of H2 permeation crossing GDYs, indicating that the GDY membrane could be a potential candidate for H2 purification.
Collapse
Affiliation(s)
- Quan Liu
- Analytical and Testing Center, Anhui University of Science and Technology, Huainan 232001, China
| | - Long Cheng
- State Key Laboratory of Materials-Oriented Chemical Engineering, College of Chemical Engineering, Nanjing Tech University, 30 Puzhu Road (S), Nanjing 211816, China;
| | - Gongping Liu
- State Key Laboratory of Materials-Oriented Chemical Engineering, College of Chemical Engineering, Nanjing Tech University, 30 Puzhu Road (S), Nanjing 211816, China;
| |
Collapse
|
49
|
Liu Q, Chen M, Mao Y, Liu G. Theoretical study on Janus graphene oxide membrane for water transport. Front Chem Sci Eng 2020. [DOI: 10.1007/s11705-020-1954-5] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/23/2022]
|
50
|
Yu Y, Klauda JB. Update of the CHARMM36 United Atom Chain Model for Hydrocarbons and Phospholipids. J Phys Chem B 2020; 124:6797-6812. [PMID: 32639155 DOI: 10.1021/acs.jpcb.0c04795] [Citation(s) in RCA: 9] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/22/2022]
Abstract
Accurate lipid force field (FF) parameters used in molecular dynamics (MD) simulations are crucial for understanding the properties of lipid-containing systems and biological processes related to lipids. The last update of the CHARMM36 united atom chain model (C36UA) was in 2013 [Lee, S. J. Phys. Chem. B 2014, 118, 547 556]; it utilized CHARMM36 (C36) lipid FF parameters for headgroups and OPLS-UA Lennard-Jones (LJ) parameters for tails. Simulations with the FF were able to reproduce many experimental observables of lipid bilayers accurately, but to be more applicable for a wide range of lipids, additional FF parameter optimization was needed. In this work, we present an update of the model, named C36UAr. The parameterization included the LJ parameters for hydrocarbons and related dihedrals. Bulk liquid properties (density, heat of vaporization, isothermal compressibility, and diffusion constant) of model compounds were used as targets for the LJ parameter fitting, and dihedrals were fit to either quantum mechanical (QM) or potential of mean force (PMF) calculations using C36. Thermodynamic reweighting was used to further improve the parameters. Bilayer simulations of various lipid headgroups (phosphatidylcholine, phosphatidylethanolamine, and phosphatidylglycerol) and tails (saturated, monounsaturated, and polyunsaturated) were performed to validate the model, and significant improvements were seen in bilayer properties, including surface area, membrane thicknesses, NMR deuterium order parameters, and density profiles. C36UAr was also compared to the hydrogen mass repartitioning (HMR) method. The high accuracy and competitive efficiency shown in this study make C36UAr one of the best choices for studies of membrane structure and membrane-associated proteins.
Collapse
|