1
|
Loschwitz J, Olubiyi OO, Hub JS, Strodel B, Poojari CS. Computer simulations of protein-membrane systems. PROGRESS IN MOLECULAR BIOLOGY AND TRANSLATIONAL SCIENCE 2020; 170:273-403. [PMID: 32145948 PMCID: PMC7109768 DOI: 10.1016/bs.pmbts.2020.01.001] [Citation(s) in RCA: 19] [Impact Index Per Article: 4.8] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Indexed: 12/12/2022]
Abstract
The interactions between proteins and membranes play critical roles in signal transduction, cell motility, and transport, and they are involved in many types of diseases. Molecular dynamics (MD) simulations have greatly contributed to our understanding of protein-membrane interactions, promoted by a dramatic development of MD-related software, increasingly accurate force fields, and available computer power. In this chapter, we present available methods for studying protein-membrane systems with MD simulations, including an overview about the various all-atom and coarse-grained force fields for lipids, and useful software for membrane simulation setup and analysis. A large set of case studies is discussed.
Collapse
Affiliation(s)
- Jennifer Loschwitz
- Institute of Theoretical and Computational Chemistry, Heinrich Heine University Düsseldorf, Düsseldorf, Germany; Institute of Biological Information Processing (IBI-7: Structural Biochemistry), Forschungszentrum Jülich, Jülich, Germany
| | - Olujide O Olubiyi
- Institute of Biological Information Processing (IBI-7: Structural Biochemistry), Forschungszentrum Jülich, Jülich, Germany; Department of Pharmaceutical Chemistry, Faculty of Pharmacy, Obafemi Awolowo University, Ile-Ife, Nigeria
| | - Jochen S Hub
- Theoretical Physics and Center for Biophysics, Saarland University, Saarbrücken, Germany
| | - Birgit Strodel
- Institute of Theoretical and Computational Chemistry, Heinrich Heine University Düsseldorf, Düsseldorf, Germany; Institute of Biological Information Processing (IBI-7: Structural Biochemistry), Forschungszentrum Jülich, Jülich, Germany
| | - Chetan S Poojari
- Theoretical Physics and Center for Biophysics, Saarland University, Saarbrücken, Germany.
| |
Collapse
|
2
|
Pluhackova K, Kirsch SA, Han J, Sun L, Jiang Z, Unruh T, Böckmann RA. A Critical Comparison of Biomembrane Force Fields: Structure and Dynamics of Model DMPC, POPC, and POPE Bilayers. J Phys Chem B 2016; 120:3888-903. [DOI: 10.1021/acs.jpcb.6b01870] [Citation(s) in RCA: 107] [Impact Index Per Article: 13.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/16/2022]
Affiliation(s)
- Kristyna Pluhackova
- Computational
Biology, Department of Biology, Friedrich-Alexander-University of Erlangen-Nürnberg, Staudtstr. 5, 91058 Erlangen, Germany
| | - Sonja A. Kirsch
- Computational
Biology, Department of Biology, Friedrich-Alexander-University of Erlangen-Nürnberg, Staudtstr. 5, 91058 Erlangen, Germany
| | - Jing Han
- Computational
Biology, Department of Biology, Friedrich-Alexander-University of Erlangen-Nürnberg, Staudtstr. 5, 91058 Erlangen, Germany
| | - Liping Sun
- Computational
Biology, Department of Biology, Friedrich-Alexander-University of Erlangen-Nürnberg, Staudtstr. 5, 91058 Erlangen, Germany
| | - Zhenyan Jiang
- Computational
Biology, Department of Biology, Friedrich-Alexander-University of Erlangen-Nürnberg, Staudtstr. 5, 91058 Erlangen, Germany
| | - Tobias Unruh
- Lehrstuhl
für Kristallografie und Strukturphysik, Department Physik, Friedrich-Alexander-University of Erlangen-Nürnberg, Staudtstr. 3, 91058 Erlangen, Germany
| | - Rainer A. Böckmann
- Computational
Biology, Department of Biology, Friedrich-Alexander-University of Erlangen-Nürnberg, Staudtstr. 5, 91058 Erlangen, Germany
| |
Collapse
|
3
|
Atomistic resolution structure and dynamics of lipid bilayers in simulations and experiments. BIOCHIMICA ET BIOPHYSICA ACTA-BIOMEMBRANES 2016; 1858:2512-2528. [PMID: 26809025 DOI: 10.1016/j.bbamem.2016.01.019] [Citation(s) in RCA: 56] [Impact Index Per Article: 7.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 11/20/2015] [Revised: 01/15/2016] [Accepted: 01/19/2016] [Indexed: 01/18/2023]
Abstract
Accurate details on the sampled atomistic resolution structures of lipid bilayers can be experimentally obtained by measuring C-H bond order parameters, spin relaxation rates and scattering form factors. These parameters can be also directly calculated from the classical atomistic resolution molecular dynamics simulations (MD) and compared to the experimentally achieved results. This comparison measures the simulation model quality with respect to 'reality'. If agreement is sufficient, the simulation model gives an atomistic structural interpretation of the acquired experimental data. Significant advance of MD models is made by jointly interpreting different experiments using the same structural model. Here we focus on phosphatidylcholine lipid bilayers, which out of all model membranes have been studied mostly by experiments and simulations, leading to the largest available dataset. From the applied comparisons we conclude that the acyl chain region structure and rotational dynamics are generally well described in simulation models. Also changes with temperature, dehydration and cholesterol concentration are qualitatively correctly reproduced. However, the quality of the underlying atomistic resolution structural changes is uncertain. Even worse, when focusing on the lipid bilayer properties at the interfacial region, e.g. glycerol backbone and choline structures, and cation binding, many simulation models produce an inaccurate description of experimental data. Thus extreme care must be applied when simulations are applied to understand phenomena where the interfacial region plays a significant role. This work is done by the NMRlipids Open Collaboration project running at https://nmrlipids.blogspot.fi and https://github.com/NMRLipids. This article is part of a Special Issue entitled: Biosimulations edited by Ilpo Vattulainen and Tomasz Róg.
Collapse
|
4
|
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
|
5
|
Botan A, Favela-Rosales F, Fuchs PFJ, Javanainen M, Kanduč M, Kulig W, Lamberg A, Loison C, Lyubartsev A, Miettinen MS, Monticelli L, Määttä J, Ollila OHS, Retegan M, Róg T, Santuz H, Tynkkynen J. Toward Atomistic Resolution Structure of Phosphatidylcholine Headgroup and Glycerol Backbone at Different Ambient Conditions. J Phys Chem B 2015; 119:15075-88. [PMID: 26509669 PMCID: PMC4677354 DOI: 10.1021/acs.jpcb.5b04878] [Citation(s) in RCA: 95] [Impact Index Per Article: 10.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/21/2015] [Revised: 10/19/2015] [Indexed: 11/28/2022]
Abstract
Phospholipids are essential building blocks of biological membranes. Despite a vast amount of very accurate experimental data, the atomistic resolution structures sampled by the glycerol backbone and choline headgroup in phoshatidylcholine bilayers are not known. Atomistic resolution molecular dynamics simulations have the potential to resolve the structures, and to give an arrestingly intuitive interpretation of the experimental data, but only if the simulations reproduce the data within experimental accuracy. In the present work, we simulated phosphatidylcholine (PC) lipid bilayers with 13 different atomistic models, and compared simulations with NMR experiments in terms of the highly structurally sensitive C-H bond vector order parameters. Focusing on the glycerol backbone and choline headgroups, we showed that the order parameter comparison can be used to judge the atomistic resolution structural accuracy of the models. Accurate models, in turn, allow molecular dynamics simulations to be used as an interpretation tool that translates these NMR data into a dynamic three-dimensional representation of biomolecules in biologically relevant conditions. In addition to lipid bilayers in fully hydrated conditions, we reviewed previous experimental data for dehydrated bilayers and cholesterol-containing bilayers, and interpreted them with simulations. Although none of the existing models reached experimental accuracy, by critically comparing them we were able to distill relevant chemical information: (1) increase of choline order parameters indicates the P-N vector tilting more parallel to the membrane, and (2) cholesterol induces only minor changes to the PC (glycerol backbone) structure. This work has been done as a fully open collaboration, using nmrlipids.blogspot.fi as a communication platform; all the scientific contributions were made publicly on this blog. During the open research process, the repository holding our simulation trajectories and files ( https://zenodo.org/collection/user-nmrlipids ) has become the most extensive publicly available collection of molecular dynamics simulation trajectories of lipid bilayers.
Collapse
Affiliation(s)
- Alexandru Botan
- Institut
Lumière Matière, UMR5306 Université
Lyon 1-CNRS, Université de Lyon, 69622 Villeurbanne, France
| | - Fernando Favela-Rosales
- Departamento
de Física, Centro de Investigación
y de Estudios Avanzados del IPN, Apartado, Postal 14-740, Mexico City, 07000 México
D.F., México
| | - Patrick F. J. Fuchs
- Institut
Jacques Monod, UMR 7592 CNRS, Université Paris
Diderot, Sorbonne, Paris Cité, F-75205 Paris, France
| | - Matti Javanainen
- Department
of Physics, Tampere University of Technology, Tampere, 33101 Finland
| | - Matej Kanduč
- Fachbereich
Physik, Freie Universität Berlin, Berlin, 14195 Germany
| | - Waldemar Kulig
- Department
of Physics, Tampere University of Technology, Tampere, 33101 Finland
| | - Antti Lamberg
- Department
of Chemical Engineering, Kyoto University, 615-8510 Kyoto, Japan
| | - Claire Loison
- Institut
Lumière Matière, UMR5306 Université
Lyon 1-CNRS, Université de Lyon, 69622 Villeurbanne, France
| | - Alexander Lyubartsev
- Division
of Physical Chemistry, Department of Materials and Environmental Chemistry, Stockholm University, S-106 91 Stockholm, Sweden
| | | | - Luca Monticelli
- Institut
de Biologie et Chimie des Protéines (IBCP), CNRS UMR 5086, Lyon 69 367, France
| | - Jukka Määttä
- Department of Chemistry, Aalto University, 00076 Aalto, Finland
| | - O. H. Samuli Ollila
- Department of Neuroscience and Biomedical Engineering, Aalto University, 00076 Aalto, Finland
| | - Marius Retegan
- Max Planck Institute
for Chemical Energy Conversion, Stiftstr. 34-38, 45470 Mülheim an der Ruhr, Germany
| | - Tomasz Róg
- Department
of Physics, Tampere University of Technology, Tampere, 33101 Finland
| | - Hubert Santuz
- INSERM, UMR_S 1134, DSIMB, Paris 75739, France
- Université
Paris Diderot, Sorbonne Paris Cité, UMR_S 1134, Paris, France
- Institut
National de la Transfusion Sanguine (INTS), Paris 75739, France
- Laboratoire d’Excellence GR-Ex, Paris 75015, France
| | - Joona Tynkkynen
- Department
of Physics, Tampere University of Technology, Tampere, 33101 Finland
| |
Collapse
|
6
|
Lockhart C, Klimov DK. Binding of Aβ peptide creates lipid density depression in DMPC bilayer. BIOCHIMICA ET BIOPHYSICA ACTA-BIOMEMBRANES 2014; 1838:2678-88. [DOI: 10.1016/j.bbamem.2014.07.010] [Citation(s) in RCA: 19] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Received: 05/09/2014] [Revised: 07/04/2014] [Accepted: 07/07/2014] [Indexed: 12/12/2022]
|
7
|
Structural dynamics of the cell wall precursor lipid II in the presence and absence of the lantibiotic nisin. BIOCHIMICA ET BIOPHYSICA ACTA-BIOMEMBRANES 2014; 1838:3061-8. [PMID: 25128154 DOI: 10.1016/j.bbamem.2014.07.024] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 06/24/2014] [Revised: 07/23/2014] [Accepted: 07/25/2014] [Indexed: 01/06/2023]
Abstract
Representing a physiological "Achilles' heel", the cell wall precursor lipid II (LII) is a prime target for various classes of antibiotics. Over the years LII-binding agents have been recognized as promising candidates and templates in the search for new antibacterial compounds to complement or replace existing drugs. To elucidate the molecular structural basis underlying LII functional mechanism and to better understand if and how lantibiotic binding alters the molecular behavior of LII, we performed molecular dynamics (MD) simulations of phospholipid membrane-embedded LII in the absence and presence of the LII-binding lantibiotic nisin. In a series of 2×4 independent, unbiased 100ns MD simulations we sampled the conformational dynamics of nine LII as well as nine LII-nisin complexes embedded in an aqueous 150mM NaCl/POPC phospholipid membrane environment. We found that nisin binding to LII induces a reduction of LII mobility and flexibility, an outward shift of the LII pentapeptide, an inward movement of the LII disaccharide section, and an overall deeper insertion of the LII tail group into the membrane. The latter effect might indicate an initial step in adopting a stabilizing, scaffold-like structure in the process of nisin-induced membrane leakage. At the same time nisin conformation and LII interaction remain similar to the 1WCO LII-nisin NMR solution structure.
Collapse
|
8
|
Dékány Fraňová M, Vattulainen I, Samuli Ollila O. Can pyrene probes be used to measure lateral pressure profiles of lipid membranes? Perspective through atomistic simulations. BIOCHIMICA ET BIOPHYSICA ACTA-BIOMEMBRANES 2014; 1838:1406-11. [DOI: 10.1016/j.bbamem.2014.01.030] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Received: 10/18/2013] [Revised: 01/07/2014] [Accepted: 01/30/2014] [Indexed: 10/25/2022]
|
9
|
Paloncýová M, DeVane R, Murch B, Berka K, Otyepka M. Amphiphilic Drug-Like Molecules Accumulate in a Membrane below the Head Group Region. J Phys Chem B 2014; 118:1030-9. [DOI: 10.1021/jp4112052] [Citation(s) in RCA: 80] [Impact Index Per Article: 8.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/21/2023]
Affiliation(s)
- Markéta Paloncýová
- Department
of Physical Chemistry, Faculty of Science, Regional Centre
of Advanced Technologies and Materials, Palacký University Olomouc, tř.
17. listopadu 12, 771 46 Olomouc, Czech Republic
| | - Russell DeVane
- Corporate Modeling & Simulation, Procter & Gamble, 8611 Beckett Road, West Chester, Ohio 45069, United States
| | - Bruce Murch
- Corporate Modeling & Simulation, Procter & Gamble, 8611 Beckett Road, West Chester, Ohio 45069, United States
| | - Karel Berka
- Department
of Physical Chemistry, Faculty of Science, Regional Centre
of Advanced Technologies and Materials, Palacký University Olomouc, tř.
17. listopadu 12, 771 46 Olomouc, Czech Republic
| | - Michal Otyepka
- Department
of Physical Chemistry, Faculty of Science, Regional Centre
of Advanced Technologies and Materials, Palacký University Olomouc, tř.
17. listopadu 12, 771 46 Olomouc, Czech Republic
| |
Collapse
|
10
|
Rabinovich AL, Lyubartsev AP. Computer simulation of lipid membranes: Methodology and achievements. POLYMER SCIENCE SERIES C 2013. [DOI: 10.1134/s1811238213070060] [Citation(s) in RCA: 23] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 01/30/2023]
|
11
|
Lama D, Modi V, Sankararamakrishnan R. Behavior of solvent-exposed hydrophobic groove in the anti-apoptotic Bcl-XL protein: clues for its ability to bind diverse BH3 ligands from MD simulations. PLoS One 2013; 8:e54397. [PMID: 23468841 PMCID: PMC3585337 DOI: 10.1371/journal.pone.0054397] [Citation(s) in RCA: 9] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/03/2012] [Accepted: 12/13/2012] [Indexed: 11/19/2022] Open
Abstract
Bcl-XL is a member of Bcl-2 family of proteins involved in the regulation of intrinsic pathway of apoptosis. Its overexpression in many human cancers makes it an important target for anti-cancer drugs. Bcl-XL interacts with the BH3 domain of several pro-apoptotic Bcl-2 partners. This helical bundle protein has a pronounced hydrophobic groove which acts as a binding region for the BH3 domains. Eight independent molecular dynamics simulations of the apo/holo forms of Bcl-XL were carried out to investigate the behavior of solvent-exposed hydrophobic groove. The simulations used either a twin-range cut-off or particle mesh Ewald (PME) scheme to treat long-range interactions. Destabilization of the BH3 domain-containing helix H2 was observed in all four twin-range cut-off simulations. Most of the other major helices remained stable. The unwinding of H2 can be related to the ability of Bcl-XL to bind diverse BH3 ligands. The loss of helical character can also be linked to the formation of homo- or hetero-dimers in Bcl-2 proteins. Several experimental studies have suggested that exposure of BH3 domain is a crucial event before they form dimers. Thus unwinding of H2 seems to be functionally very important. The four PME simulations, however, revealed a stable helix H2. It is possible that the H2 unfolding might occur in PME simulations at longer time scales. Hydrophobic residues in the hydrophobic groove are involved in stable interactions among themselves. The solvent accessible surface areas of bulky hydrophobic residues in the groove are significantly buried by the loop LB connecting the helix H2 and subsequent helix. These observations help to understand how the hydrophobic patch in Bcl-XL remains stable in the solvent-exposed state. We suggest that both the destabilization of helix H2 and the conformational heterogeneity of loop LB are important factors for binding of diverse ligands in the hydrophobic groove of Bcl-XL.
Collapse
Affiliation(s)
- Dilraj Lama
- Department of Biological Sciences & Bioengineering, Indian Institute of Technology Kanpur, Kanpur, India
| | - Vivek Modi
- Department of Biological Sciences & Bioengineering, Indian Institute of Technology Kanpur, Kanpur, India
| | | |
Collapse
|
12
|
Li J, Liu S, Lakshminarayanan R, Bai Y, Pervushin K, Verma C, Beuerman RW. Molecular simulations suggest how a branched antimicrobial peptide perturbs a bacterial membrane and enhances permeability. BIOCHIMICA ET BIOPHYSICA ACTA-BIOMEMBRANES 2012; 1828:1112-21. [PMID: 23274275 DOI: 10.1016/j.bbamem.2012.12.015] [Citation(s) in RCA: 55] [Impact Index Per Article: 4.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 06/28/2012] [Revised: 12/03/2012] [Accepted: 12/18/2012] [Indexed: 11/29/2022]
Abstract
A covalently, branched antimicrobial peptide (BAMP) B2088 demonstrating enhanced antimicrobial effects and without additional toxicity when compared to its linear counterpart, has been developed. Atomistic molecular dynamics simulations have been used to investigate the mode of interaction of B2088 with model bacterial and mammalian membranes. These simulations suggest that both long-range electrostatic interactions and short-range hydrogen bonding play important roles in steering B2088 toward the negatively charged membranes. The reason why B2088 is selective towards the bacterial membrane is postulated to be the greater density of negative charges on the bacterial membrane which enables rapid accumulation of B2088 on the bacterial membrane to a high surface concentration, stabilizing it through excess hydrogen bond formation. The majority of hydrogen bonds are seen between the side chains of the basic residues (Arg or Lys) with the PO4 groups of lipids. In particular, formation of the bidentate hydrogen bonds between the guanidinium group of Arg and PO4 groups are found to be more favorable, both geometrically and energetically. Moreover, the planar gaunidinium group and its hydrophobic character enable the Arg side chains to solvate into the hydrophobic membrane. Structural perturbation of the bacterial membrane is found to be concentration dependent and is significant at higher concentrations of B2088, resulting in a large number of water translocations across the bacterial membrane. These simulations enhance our understanding of the action mechanism of a covalently branched antimicrobial peptide with model membranes and provide practical guidance for the design of new antimicrobial peptides.
Collapse
Affiliation(s)
- Jianguo Li
- Singapore Eye Research Institute, Singapore, Singapore
| | | | | | | | | | | | | |
Collapse
|
13
|
Poger D, Mark AE. Lipid Bilayers: The Effect of Force Field on Ordering and Dynamics. J Chem Theory Comput 2012; 8:4807-17. [PMID: 26605633 DOI: 10.1021/ct300675z] [Citation(s) in RCA: 66] [Impact Index Per Article: 5.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/12/2022]
Abstract
The sensitivity of the structure and dynamics of a fully hydrated pure bilayer of 1,2-dipalmitoyl-sn-glycero-3-phosphocholine (DPPC) in molecular dynamics simulations to changes in force-field and simulation parameters has been assessed. Three related force fields (the Gromos 54A7 force field, a Gromos 53A6-derived parameter set and a variant of the Berger parameters) in combination with either particle-mesh Ewald (PME) or a reaction field (RF) were compared. Structural properties such as the area per lipid, carbon-deuterium order parameters, electron density profile and bilayer thicknesses, are reproduced by all the parameter sets within the uncertainty of the available experimental data. However, there are clear differences in the ordering of the glycerol backbone and choline headgroup, and the orientation of the headgroup dipole. In some cases, the degree of ordering was reminiscent of a liquid-ordered phase. It is also shown that, although the lateral diffusion of the lipids in the plane of the bilayer is often used to validate lipid force fields, because of the uncertainty in the experimental measurements and the fact that the lateral diffusion is dependent on the choice of the simulation conditions, it should not be employed as a measure of quality. Finally, the simulations show that the effect of small changes in force-field parameters on the structure and dynamics of a bilayer is more significant than the treatment of the long-range electrostatic interactions using RF or PME. Overall, the Gromos 54A7 best reproduced the range of experimental data examined.
Collapse
Affiliation(s)
- David Poger
- The University of Queensland, School of Chemistry and Molecular Biosciences, Brisbane QLD 4072, Australia
| | - Alan E Mark
- The University of Queensland, School of Chemistry and Molecular Biosciences, Brisbane QLD 4072, Australia.,The University of Queensland, Institute for Molecular Bioscience, Brisbane QLD 4072, Australia
| |
Collapse
|
14
|
Jämbeck JPM, Lyubartsev AP. An Extension and Further Validation of an All-Atomistic Force Field for Biological Membranes. J Chem Theory Comput 2012; 8:2938-48. [PMID: 26592132 DOI: 10.1021/ct300342n] [Citation(s) in RCA: 346] [Impact Index Per Article: 28.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/08/2023]
Abstract
Biological membranes are versatile in composition and host intriguing molecular processes. In order to be able to study these systems, an accurate model Hamiltonian or force field (FF) is a necessity. Here, we report the results of our extension of earlier developed all-atomistic FF parameters for fully saturated phospholipids that complements an earlier parameter set for saturated phosphatidylcholine lipids (J. Phys. Chem. B, 2012, 116, 3164-3179). The FF, coined Slipids (Stockholm lipids), now also includes parameters for unsaturated phosphatidylcholine and phosphatidylethanolamine lipids, e.g., POPC, DOPC, SOPC, POPE, and DOPE. As the extended set of parameters is derived with the same philosophy as previously applied, the resulting FF has been developed in a fully consistent manner. The capabilities of Slipids are demonstrated by performing long simulations without applying any surface tension and using the correct isothermal-isobaric (NPT) ensemble for a range of temperatures and carefully comparing a number of properties with experimental findings. Results show that several structural properties are very well reproduced, such as scattering form factors, NMR order parameters, thicknesses, and area per lipid. Thermal dependencies of different thicknesses and area per lipid are reproduced as well. Lipid diffusion is systematically slightly underestimated, whereas the normalized lipid diffusion follows the experimental trends. This is believed to be due to the lack of collective movement in the relatively small bilayer patches used. Furthermore, the compatibility with amino acid FFs from the AMBER family is tested in explicit transmembrane complexes of the WALP23 peptide with DLPC and DOPC bilayers, and this shows that Slipids can be used to study more complex and biologically relevant systems.
Collapse
Affiliation(s)
- Joakim P M Jämbeck
- Division of Physical Chemistry, Arrhenius Laboratory, Stockholm University, Stockholm, SE-10691, Sweden
| | - Alexander P Lyubartsev
- Division of Physical Chemistry, Arrhenius Laboratory, Stockholm University, Stockholm, SE-10691, Sweden
| |
Collapse
|
15
|
Jämbeck JPM, Lyubartsev AP. Derivation and systematic validation of a refined all-atom force field for phosphatidylcholine lipids. J Phys Chem B 2012; 116:3164-79. [PMID: 22352995 PMCID: PMC3320744 DOI: 10.1021/jp212503e] [Citation(s) in RCA: 409] [Impact Index Per Article: 34.1] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/27/2011] [Revised: 02/18/2012] [Indexed: 11/29/2022]
Abstract
An all-atomistic force field (FF) has been developed for fully saturated phospholipids. The parametrization has been largely based on high-level ab initio calculations in order to keep the empirical input to a minimum. Parameters for the lipid chains have been developed based on knowledge about bulk alkane liquids, for which thermodynamic and dynamic data are excellently reproduced. The FFs ability to simulate lipid bilayers in the liquid crystalline phase in a tensionless ensemble was tested in simulations of three lipids: 1,2-diauroyl-sn-glycero-3-phospocholine (DLPC), 1,2-dimyristoyl-sn-glycero-3-phosphocholine (DMPC), and 1,2-dipalmitoyl-sn-glycero-3-phospcholine (DPPC). Computed areas and volumes per lipid, and three different kinds of bilayer thicknesses, have been investigated. Most importantly NMR order parameters and scattering form factors agree in an excellent manner with experimental data under a range of temperatures. Further, the compatibility with the AMBER FF for biomolecules as well as the ability to simulate bilayers in gel phase was demonstrated. Overall, the FF presented here provides the important balance between the hydrophilic and hydrophobic forces present in lipid bilayers and therefore can be used for more complicated studies of realistic biological membranes with protein insertions.
Collapse
Affiliation(s)
- Joakim P. M. Jämbeck
- Division of Physical Chemistry, Arrhenius Laboratory, Stockholm University, Stockholm, SE-10691, Sweden
| | - Alexander P. Lyubartsev
- Division of Physical Chemistry, Arrhenius Laboratory, Stockholm University, Stockholm, SE-10691, Sweden
| |
Collapse
|
16
|
Intramolecular hydrogen bonding in articaine can be related to superior bone tissue penetration: a molecular dynamics study. Biophys Chem 2010; 154:18-25. [PMID: 21227568 DOI: 10.1016/j.bpc.2010.12.002] [Citation(s) in RCA: 29] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/15/2010] [Revised: 12/08/2010] [Accepted: 12/08/2010] [Indexed: 11/24/2022]
Abstract
Local anesthetics (LAs) are drugs that cause reversible loss of nociception during surgical procedures. Articaine is a commonly used LA in dentistry that has proven to be exceptionally effective in penetrating bone tissue and induce anesthesia on posterior teeth in maxilla and mandibula. In the present study, our aim was to gain a deeper understanding of the penetration of articaine through biological membranes by studying the interactions of articaine with a phospholipid membrane. Our approach involves Langmuir monolayer experiments combined with molecular dynamics simulations. Membrane permeability of LAs can be modulated by pH due to a titratable amine group with a pKa value close to physiological pH. A change in protonation state is thus known to act as a lipophilicity switch in LAs. Our study shows that articaine has an additional unique lipophilicity switch in its ability to form an intramolecular hydrogen bond. We suggest this intramolecular hydrogen bond as a novel and additional solvent-dependent mechanism for modulation of lipophilicity of articaine which may enhance its diffusion through membranes and connective tissue.
Collapse
|
17
|
Klauda JB, Venable RM, Freites JA, O’Connor JW, Tobias DJ, Mondragon-Ramirez C, Vorobyov I, MacKerell AD, Pastor RW. Update of the CHARMM all-atom additive force field for lipids: validation on six lipid types. J Phys Chem B 2010; 114:7830-43. [PMID: 20496934 PMCID: PMC2922408 DOI: 10.1021/jp101759q] [Citation(s) in RCA: 3115] [Impact Index Per Article: 222.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
A significant modification to the additive all-atom CHARMM lipid force field (FF) is developed and applied to phospholipid bilayers with both choline and ethanolamine containing head groups and with both saturated and unsaturated aliphatic chains. Motivated by the current CHARMM lipid FF (C27 and C27r) systematically yielding values of the surface area per lipid that are smaller than experimental estimates and gel-like structures of bilayers well above the gel transition temperature, selected torsional, Lennard-Jones and partial atomic charge parameters were modified by targeting both quantum mechanical (QM) and experimental data. QM calculations ranging from high-level ab initio calculations on small molecules to semiempirical QM studies on a 1,2-dipalmitoyl-sn-phosphatidylcholine (DPPC) bilayer in combination with experimental thermodynamic data were used as target data for parameter optimization. These changes were tested with simulations of pure bilayers at high hydration of the following six lipids: DPPC, 1,2-dimyristoyl-sn-phosphatidylcholine (DMPC), 1,2-dilauroyl-sn-phosphatidylcholine (DLPC), 1-palmitoyl-2-oleoyl-sn-phosphatidylcholine (POPC), 1,2-dioleoyl-sn-phosphatidylcholine (DOPC), and 1-palmitoyl-2-oleoyl-sn-phosphatidylethanolamine (POPE); simulations of a low hydration DOPC bilayer were also performed. Agreement with experimental surface area is on average within 2%, and the density profiles agree well with neutron and X-ray diffraction experiments. NMR deuterium order parameters (S(CD)) are well predicted with the new FF, including proper splitting of the S(CD) for the aliphatic carbon adjacent to the carbonyl for DPPC, POPE, and POPC bilayers. The area compressibility modulus and frequency dependence of (13)C NMR relaxation rates of DPPC and the water distribution of low hydration DOPC bilayers also agree well with experiment. Accordingly, the presented lipid FF, referred to as C36, allows for molecular dynamics simulations to be run in the tensionless ensemble (NPT), and is anticipated to be of utility for simulations of pure lipid systems as well as heterogeneous systems including membrane proteins.
Collapse
Affiliation(s)
- Jeffery B. Klauda
- Department of Chemical and Biomolecular Engineering, University of Maryland, College Park, MD 20742
| | - Richard M. Venable
- Laboratory of Computational Biology, National Heart, Lung, and Blood Institute, National Institutes of Health, Bethesda, MD 20892
| | - J. Alfredo Freites
- Department of Chemistry, University of California, Irvine, CA 92697-2025
| | - Joseph W. O’Connor
- Department of Chemical and Biomolecular Engineering, University of Maryland, College Park, MD 20742
| | - Douglas J. Tobias
- Department of Chemistry, University of California, Irvine, CA 92697-2025
| | - Carlos Mondragon-Ramirez
- Department of Pharmaceutical Sciences, 20 Penn Street HSF II, University of Maryland, Baltimore, Maryland 21201
| | - Igor Vorobyov
- Department of Pharmaceutical Sciences, 20 Penn Street HSF II, University of Maryland, Baltimore, Maryland 21201
| | - Alexander D. MacKerell
- Department of Pharmaceutical Sciences, 20 Penn Street HSF II, University of Maryland, Baltimore, Maryland 21201
| | - Richard W. Pastor
- Laboratory of Computational Biology, National Heart, Lung, and Blood Institute, National Institutes of Health, Bethesda, MD 20892
| |
Collapse
|