1
|
Duboué-Dijon E, Javanainen M, Delcroix P, Jungwirth P, Martinez-Seara H. A practical guide to biologically relevant molecular simulations with charge scaling for electronic polarization. J Chem Phys 2021; 153:050901. [PMID: 32770904 DOI: 10.1063/5.0017775] [Citation(s) in RCA: 51] [Impact Index Per Article: 17.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/23/2022] Open
Abstract
Molecular simulations can elucidate atomistic-level mechanisms of key biological processes, which are often hardly accessible to experiment. However, the results of the simulations can only be as trustworthy as the underlying simulation model. In many of these processes, interactions between charged moieties play a critical role. Current empirical force fields tend to overestimate such interactions, often in a dramatic way, when polyvalent ions are involved. The source of this shortcoming is the missing electronic polarization in these models. Given the importance of such biomolecular systems, there is great interest in fixing this deficiency in a computationally inexpensive way without employing explicitly polarizable force fields. Here, we review the electronic continuum correction approach, which accounts for electronic polarization in a mean-field way, focusing on its charge scaling variant. We show that by pragmatically scaling only the charged molecular groups, we qualitatively improve the charge-charge interactions without extra computational costs and benefit from decades of force field development on biomolecular force fields.
Collapse
Affiliation(s)
- E Duboué-Dijon
- CNRS, Université de Paris, UPR 9080, Laboratoire de Biochimie Théorique, 13 rue Pierre et Marie Curie, 75005 Paris, France
| | - M Javanainen
- Institute of Organic Chemistry and Biochemistry, Czech Academy of Sciences, Flemingovo nam. 2, Prague 6 166 10, Czech Republic
| | - P Delcroix
- Institute of Organic Chemistry and Biochemistry, Czech Academy of Sciences, Flemingovo nam. 2, Prague 6 166 10, Czech Republic
| | - P Jungwirth
- Institute of Organic Chemistry and Biochemistry, Czech Academy of Sciences, Flemingovo nam. 2, Prague 6 166 10, Czech Republic
| | - H Martinez-Seara
- Institute of Organic Chemistry and Biochemistry, Czech Academy of Sciences, Flemingovo nam. 2, Prague 6 166 10, Czech Republic
| |
Collapse
|
2
|
Fornasier F, Souza LMP, Souza FR, Reynaud F, Pimentel AS. Lipophilicity of Coarse-Grained Cholesterol Models. J Chem Inf Model 2020; 60:569-577. [DOI: 10.1021/acs.jcim.9b00830] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Franccesca Fornasier
- Departamento de Quı́mica, Pontifı́cia Universidade Católica do Rio de Janeiro Rio de Janeiro, RJ 22453-900, Brazil
| | - Lucas M. P. Souza
- Departamento de Quı́mica, Pontifı́cia Universidade Católica do Rio de Janeiro Rio de Janeiro, RJ 22453-900, Brazil
| | - Felipe R. Souza
- Departamento de Quı́mica, Pontifı́cia Universidade Católica do Rio de Janeiro Rio de Janeiro, RJ 22453-900, Brazil
| | - Franceline Reynaud
- Institut Galien Paris-Sud, CNRS 8612, Université Paris-Saclay, 92290 Châtenay, Malabry, France
- Faculdade de Farmácia, Universidade Federal do Rio de Janeiro, Rio de Janeiro, RJ 21941-902, Brazil
| | - Andre S. Pimentel
- Departamento de Quı́mica, Pontifı́cia Universidade Católica do Rio de Janeiro Rio de Janeiro, RJ 22453-900, Brazil
| |
Collapse
|
3
|
Chu H, Peng X, Li Y, Zhang Y, Min H, Li G. Polarizable atomic multipole-based force field for DOPC and POPE membrane lipids. Mol Phys 2018. [DOI: 10.1080/00268976.2018.1436201] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/18/2022]
Affiliation(s)
- Huiying Chu
- Laboratory of Molecular Modeling and Design, State Key Laboratory of Molecular Reaction Dynamics, Dalian Institute of Chemical Physics, Chinese Academy of Science, Liaoning, China
| | - Xiangda Peng
- Laboratory of Molecular Modeling and Design, State Key Laboratory of Molecular Reaction Dynamics, Dalian Institute of Chemical Physics, Chinese Academy of Science, Liaoning, China
- University of Chinese Academy of Sciences, Chinese Academy of Sciences, Beijing, China
| | - Yan Li
- Laboratory of Molecular Modeling and Design, State Key Laboratory of Molecular Reaction Dynamics, Dalian Institute of Chemical Physics, Chinese Academy of Science, Liaoning, China
| | - Yuebin Zhang
- Laboratory of Molecular Modeling and Design, State Key Laboratory of Molecular Reaction Dynamics, Dalian Institute of Chemical Physics, Chinese Academy of Science, Liaoning, China
| | - Hanyi Min
- Chinese Academy of Medical Science & Peking Union Medical College Hospital, Ophthalmology, Beijing, China
| | - Guohui Li
- Laboratory of Molecular Modeling and Design, State Key Laboratory of Molecular Reaction Dynamics, Dalian Institute of Chemical Physics, Chinese Academy of Science, Liaoning, China
| |
Collapse
|
4
|
A Polarizable Atomic Multipole-Based Force Field for Molecular Dynamics Simulations of Anionic Lipids. Molecules 2017; 23:molecules23010077. [PMID: 29301229 PMCID: PMC6017617 DOI: 10.3390/molecules23010077] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/31/2017] [Revised: 12/26/2017] [Accepted: 12/28/2017] [Indexed: 11/18/2022] Open
Abstract
In all of the classical force fields, electrostatic interaction is simply treated and explicit electronic polarizability is neglected. The condensed-phase polarization, relative to the gas-phase charge distributions, is commonly accounted for in an average way by increasing the atomic charges, which remain fixed throughout simulations. Based on the lipid polarizable force field DMPC and following the same framework as Atomic Multipole Optimized Energetics for BiomoleculAr (AMOEBA) simulation, the present effort expands the force field to new anionic lipid models, in which the new lipids contain DMPG and POPS. The parameters are compatible with the AMOEBA force field, which includes water, ions, proteins, etc. The charge distribution of each atom is represented by the permanent atomic monopole, dipole and quadrupole moments, which are derived from the ab initio gas phase calculations. Many-body polarization including the inter- and intramolecular polarization is modeled in a consistent manner with distributed atomic polarizabilities. Molecular dynamics simulations of the two aqueous DMPG and POPS membrane bilayer systems, consisting of 72 lipids with water molecules, were then carried out to validate the force field parameters. Membrane width, area per lipid, volume per lipid, deuterium order parameters, electron density profile, electrostatic potential difference between the center of the bilayer and water are all calculated, and compared with limited experimental data.
Collapse
|
5
|
Chu H, Cao L, Peng X, Li G. Polarizable force field development for lipids and their efficient applications in membrane proteins. WILEY INTERDISCIPLINARY REVIEWS-COMPUTATIONAL MOLECULAR SCIENCE 2017. [DOI: 10.1002/wcms.1312] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 01/12/2023]
Affiliation(s)
- Huiying Chu
- Laboratory of Molecular Modeling and Design, State Key Laboratory of Molecular Reaction Dynamics; Dalian Institute of Chemical Physics, Chinese Academy of Science; Dalian China
| | - Liaoran Cao
- Laboratory of Molecular Modeling and Design, State Key Laboratory of Molecular Reaction Dynamics; Dalian Institute of Chemical Physics, Chinese Academy of Science; Dalian China
| | - Xiangda Peng
- Laboratory of Molecular Modeling and Design, State Key Laboratory of Molecular Reaction Dynamics; Dalian Institute of Chemical Physics, Chinese Academy of Science; Dalian China
| | - Guohui Li
- Laboratory of Molecular Modeling and Design, State Key Laboratory of Molecular Reaction Dynamics; Dalian Institute of Chemical Physics, Chinese Academy of Science; Dalian China
| |
Collapse
|
6
|
Witzke S, List NH, Olsen JMH, Steinmann C, Petersen M, Beerepoot MTP, Kongsted J. An averaged polarizable potential for multiscale modeling in phospholipid membranes. J Comput Chem 2017; 38:601-611. [PMID: 28160294 DOI: 10.1002/jcc.24718] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/11/2015] [Revised: 12/07/2016] [Accepted: 12/09/2016] [Indexed: 01/28/2023]
Abstract
A set of average atom-centered charges and polarizabilities has been developed for three types of phospholipids for use in polarizable embedding calculations. The lipids investigated are 1,2-dimyristoyl-sn-glycero-3-phosphocholine, 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine, and 1-palmitoyl-2-oleoyl-sn-glycerol-3-phospho-L-serine given their common use both in experimental and computational studies. The charges, and to a lesser extent the polarizabilities, are found to depend strongly on the molecular conformation of the lipids. Furthermore, the importance of explicit polarization is underlined for the description of larger assemblies of lipids, that is, membranes. In conclusion, we find that specially developed polarizable parameters are needed for embedding calculations in membranes, while common non-polarizable point-charge force fields usually perform well enough for structural and dynamical studies. © 2017 Wiley Periodicals, Inc.
Collapse
Affiliation(s)
- Sarah Witzke
- Department of Physics, Chemistry and Pharmacy, University of Southern Denmark, Odense M, DK-5230, Denmark
| | - Nanna Holmgaard List
- Division of Theoretical Chemistry and Biology, School of Biotechnology, KTH Royal Institute of Technology, Stockholm, SE-106 91, Sweden
| | | | - Casper Steinmann
- Centre for Computational Chemistry, School of Chemistry, University of Bristol, Bristol, BS8 1TS, United Kingdom
| | - Michael Petersen
- Department of Physics, Chemistry and Pharmacy, University of Southern Denmark, Odense M, DK-5230, Denmark
| | - Maarten T P Beerepoot
- Centre for Theoretical and Computational Chemistry, Department of Chemistry, University of Tromsø-The Arctic University of Norway, Tromsø, N-9037, Norway
| | - Jacob Kongsted
- Department of Physics, Chemistry and Pharmacy, University of Southern Denmark, Odense M, DK-5230, Denmark
| |
Collapse
|
7
|
Ghumman A, Robinson D. Modulation of the Turkey β1-Adrenergic Receptor by Membrane Rafts - Insight from Molecular Dynamics. ChemistrySelect 2016. [DOI: 10.1002/slct.201600849] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/07/2022]
Affiliation(s)
- Afshan Ghumman
- School of Chemistry; University Of Nottingham; University Park Nottingham NG7 2RD UK
| | - David Robinson
- School of Chemistry; University Of Nottingham; University Park Nottingham NG7 2RD UK
| |
Collapse
|
8
|
Fast calculation of molecular total energy with ABEEMσπ/MM method – For some series of organic molecules and peptides. Chem Phys 2016. [DOI: 10.1016/j.chemphys.2016.03.004] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/23/2022]
|
9
|
Guo W, Cheng L, Chu H, Cao L, Zhang D, Liu J, Xu P, Zheng X, Li G. Some polarisable force fields for molecular dynamics simulations of lipids, and bilayers. MOLECULAR SIMULATION 2016. [DOI: 10.1080/08927022.2016.1161190] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/21/2022]
|
10
|
Lyubartsev AP, Rabinovich AL. Force Field Development for Lipid Membrane Simulations. BIOCHIMICA ET BIOPHYSICA ACTA-BIOMEMBRANES 2016; 1858:2483-2497. [PMID: 26766518 DOI: 10.1016/j.bbamem.2015.12.033] [Citation(s) in RCA: 72] [Impact Index Per Article: 9.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 11/11/2015] [Revised: 12/21/2015] [Accepted: 12/23/2015] [Indexed: 02/04/2023]
Abstract
With the rapid development of computer power and wide availability of modelling software computer simulations of realistic models of lipid membranes, including their interactions with various molecular species, polypeptides and membrane proteins have become feasible for many research groups. The crucial issue of the reliability of such simulations is the quality of the force field, and many efforts, especially in the latest several years, have been devoted to parametrization and optimization of the force fields for biomembrane modelling. In this review, we give account of the recent development in this area, covering different classes of force fields, principles of the force field parametrization, comparison of the force fields, and their experimental validation. This article is part of a Special Issue entitled: Biosimulations edited by Ilpo Vattulainen and Tomasz Róg.
Collapse
Affiliation(s)
- Alexander P Lyubartsev
- Department of Materials and Environmental Chemistry, Stockholm University, SE 106 91, Stockholm, Sweden.
| | - Alexander L Rabinovich
- Institute of Biology, Karelian Research Center, Russian Academy of Sciences, Pushkinskaya 11, Petrozavodsk, 185910, Russian Federation.
| |
Collapse
|
11
|
Awoonor-Williams E, Rowley CN. Molecular simulation of nonfacilitated membrane permeation. BIOCHIMICA ET BIOPHYSICA ACTA-BIOMEMBRANES 2015; 1858:1672-87. [PMID: 26706099 DOI: 10.1016/j.bbamem.2015.12.014] [Citation(s) in RCA: 83] [Impact Index Per Article: 9.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 10/13/2015] [Revised: 12/05/2015] [Accepted: 12/09/2015] [Indexed: 12/29/2022]
Abstract
This is a review. Non-electrolytic compounds typically cross cell membranes by passive diffusion. The rate of permeation is dependent on the chemical properties of the solute and the composition of the lipid bilayer membrane. Predicting the permeability coefficient of a solute is important in pharmaceutical chemistry and toxicology. Molecular simulation has proven to be a valuable tool for modeling permeation of solutes through a lipid bilayer. In particular, the solubility-diffusion model has allowed for the quantitative calculation of permeability coefficients. The underlying theory and computational methods used to calculate membrane permeability are reviewed. We also discuss applications of these methods to examine the permeability of solutes and the effect of membrane composition on permeability. The application of coarse grain and polarizable models is discussed. This article is part of a Special Issue entitled: Membrane Proteins edited by J.C. Gumbart and Sergei Noskov.
Collapse
Affiliation(s)
- Ernest Awoonor-Williams
- Department of Chemistry, Memorial University of Newfoundland, St. John's, NL, A1B 3X7 Canada
| | - Christopher N Rowley
- Department of Chemistry, Memorial University of Newfoundland, St. John's, NL, A1B 3X7 Canada.
| |
Collapse
|
12
|
Szklarczyk OM, Arvaniti E, van Gunsteren WF. Polarizable coarse-grained models for molecular dynamics simulation of liquid cyclohexane. J Comput Chem 2015; 36:1311-21. [DOI: 10.1002/jcc.23929] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/17/2014] [Revised: 03/09/2015] [Accepted: 03/12/2015] [Indexed: 11/08/2022]
Affiliation(s)
- Oliwia M. Szklarczyk
- Laboratory of Physical Chemistry, Department of Chemistry and Applied Biosciences; Swiss Federal Institute of Technology ETH; 8093 Zürich Switzerland
| | - Eirini Arvaniti
- Laboratory of Physical Chemistry, Department of Chemistry and Applied Biosciences; Swiss Federal Institute of Technology ETH; 8093 Zürich Switzerland
| | - Wilfred F. van Gunsteren
- Laboratory of Physical Chemistry, Department of Chemistry and Applied Biosciences; Swiss Federal Institute of Technology ETH; 8093 Zürich Switzerland
| |
Collapse
|
13
|
Baker CM. Polarizable force fields for molecular dynamics simulations of biomolecules. WILEY INTERDISCIPLINARY REVIEWS-COMPUTATIONAL MOLECULAR SCIENCE 2015. [DOI: 10.1002/wcms.1215] [Citation(s) in RCA: 100] [Impact Index Per Article: 11.1] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/17/2022]
|
14
|
Yasuda T, Kinoshita M, Murata M, Matsumori N. Detailed comparison of deuterium quadrupole profiles between sphingomyelin and phosphatidylcholine bilayers. Biophys J 2014; 106:631-8. [PMID: 24507603 DOI: 10.1016/j.bpj.2013.12.034] [Citation(s) in RCA: 51] [Impact Index Per Article: 5.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/18/2013] [Revised: 11/25/2013] [Accepted: 12/26/2013] [Indexed: 11/25/2022] Open
Abstract
Lipid rafts are microdomains rich in sphingomyelin (SM) and cholesterol (Chol). The essential question is why natural lipid rafts prefer SM rather than saturated diacyl glycerophosphocholine, although both form ordered membranes with Chol in model systems. Hence in this study, we synthesized site-specifically deuterated 1-palmitoyl-2-stearoyl-sn-glycero-3-phosphocholines that match the acyl chain length of stearoyl-SM (SSM), and compared their deuterium quadrupole coupling profiles in detail. The results suggest a deeper distribution of Chol in the SSM membranes, a lower entropic penalty upon accommodation of Chol in SSM membranes, and a higher thermal stability of acyl-chain orders in the SSM-Chol bilayers than in the 1-palmitoyl-2-stearoyl-sn-glycero-3-phosphocholine-Chol system at various Chol concentrations. The entropy effect and thermal stability should render SM a more preferred raft constituent than saturated diacyl glycerophosphocholine. Our data also demonstrate that the selective and comprehensive deuteration strategy is indispensable for accurate comparison of order profiles.
Collapse
Affiliation(s)
- Tomokazu Yasuda
- Department of Chemistry, Graduate School of Science, Osaka University, Toyonaka, Osaka, Japan; Japan Science and Technology Agency, ERATO, Lipid Active Structure Project, Toyonaka, Osaka, Japan
| | - Masanao Kinoshita
- Japan Science and Technology Agency, ERATO, Lipid Active Structure Project, Toyonaka, Osaka, Japan
| | - Michio Murata
- Department of Chemistry, Graduate School of Science, Osaka University, Toyonaka, Osaka, Japan; Japan Science and Technology Agency, ERATO, Lipid Active Structure Project, Toyonaka, Osaka, Japan
| | - Nobuaki Matsumori
- Department of Chemistry, Graduate School of Science, Osaka University, Toyonaka, Osaka, Japan.
| |
Collapse
|
15
|
Blake HL, Robinson D. QM/MM studies of contemporary and novel membrane raft fluorescent probes. Molecules 2014; 19:10230-41. [PMID: 25029071 PMCID: PMC6271554 DOI: 10.3390/molecules190710230] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/16/2014] [Revised: 07/09/2014] [Accepted: 07/10/2014] [Indexed: 11/22/2022] Open
Abstract
We have studied a number of contemporary and novel membrane probes, selected for their structural similarity to membrane raft components, in order to properly anchor themselves within a sphingolipid/cholesterol rich region. A QM/MM approach was adopted in order to understand the structural and electrostatic influences of fluorescence emission shifts of the probes in different lipid and solvation environments. The proposed modifications to the membrane probes have shown encouraging data relating not only to emission shifts within the membrane, but also their ability to anchor within a membrane raft domain and the stability to internalization within a membrane system.
Collapse
Affiliation(s)
- Hannah L Blake
- School of Chemistry, University of Nottingham, University Park, Nottingham NG7 2RD, UK.
| | - David Robinson
- School of Chemistry, University of Nottingham, University Park, Nottingham NG7 2RD, UK.
| |
Collapse
|
16
|
Yang ZZ, Wang JJ, Zhao DX. Valence state parameters of all transition metal atoms in metalloproteins-development of ABEEMσπ fluctuating charge force field. J Comput Chem 2014; 35:1690-706. [DOI: 10.1002/jcc.23676] [Citation(s) in RCA: 27] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/11/2014] [Revised: 06/01/2014] [Accepted: 06/18/2014] [Indexed: 01/08/2023]
Affiliation(s)
- Zhong-Zhi Yang
- School of Chemistry and Chemical Engineering; Department of Chemistry, Liaoning Normal University; Dalian China 116029
| | - Jian-Jiang Wang
- School of Chemistry and Chemical Engineering; Department of Chemistry, Liaoning Normal University; Dalian China 116029
| | - Dong-Xia Zhao
- School of Chemistry and Chemical Engineering; Department of Chemistry, Liaoning Normal University; Dalian China 116029
| |
Collapse
|