Abstract
Biological membranes are complex environments consisting of different types of lipids and membrane proteins. The structure of a lipid bilayer is typically difficult to study because the membrane liquid crystalline state is made up of multiple disordered lipid molecules. This complicates the description of the lipid membrane properties by the conformation of any single lipid molecule. Molecular dynamics (MD) simulations have been used extensively to investigate properties of membrane lipids, lipid vesicles, and membrane protein systems. All-atom membrane models can elucidate detailed contacts between membrane proteins and its surrounding lipids, while united-atom and coarse-grained description have allowed larger models and longer timescales up to microsecond mark to be probed. Additionally, membrane models with mixed phospholipids and lipopolysaccharide content have made it possible to model improved views of biological membranes. Here, we present an overview of commonly used lipid force fields by the biosimulation community, useful tools for membrane MD simulations, and recent advances in membrane simulations.
Keywords
- all-atom (AT)
- coarse-grain (CG)
- force field
- lipid bilayer
- united-atom (UA)
1. Introduction
The biological membrane is a selective barrier delineating the boundaries of cells and organizing cellular organelles into compartments. One of the major functions of the membrane is to regulate movements of water, ions, and certain constituents in and out of the cell or organelle. Singer and Nicholson described the membrane as a “fluid mosaic” model in which lipids and proteins freely diffuse in the membrane plane [1]. While the biological membrane is actually made up of 1.5- to 4-fold proteins by weight, it is usually described as phospholipids arranged in a bilayer.
The lipid composition in membrane is highly specific with variation identified from membranes of different organisms, different cells of the same organism, and even different membranes of the same cell [2]. Generally, a lipid consists of a polar head group and hydrophobic tail region where a cell may have over a thousand different lipid species. Formation of bilayers and other lipid structures such as micelles occurs as the hydrophobic tails orientates away from the aqueous environment while the polar head groups interact favorably with water. Changes in the bilayer structures can affect the function of the membrane, membrane-embedded proteins and complexes. The surrounding environment can greatly affect the structure and functions of proteins, much like the effects of water on water-soluble proteins. The membrane produces a complex environment for embedded proteins, and membrane thickness, fluidity, charge, curvature, and phase have been demonstrated to play a role in protein structure and function determination [3].
The structure and properties of membranes are a complex matter and cannot be easily described by a single-lipid molecule alone. Molecular dynamics (MD) simulations offer a viable alternative to study the properties of membrane and lipid-forming structures such as vesicles. This review serves as an introduction to currently available membrane lipid force fields and recent advances in membrane simulations.
2. Force fields for lipid simulation
In general, all-atom (AT), united-atom (UA), and coarse-grained (CG) are the three-membrane lipid force fields. Figure 1 illustrates the AT, UA and CG force field of a lipid by spherical representation.
2.1. All-atom (AT) MD simulation: details and details
AT MD simulation represents every atom in the system as a single interaction site. Figure 2 shows the example of AT simulation system of
These major weaknesses of the CHARMM27r force field eventually motivated the development of CHARMM36. CHARMM36 corrected most of the existing problems with lipid force fields, most importantly allowing the simulation of lipid bilayer without the use of surface tension and improved the reproducibility of the experimental deuterium order parameters in the glycerol and upper chain regions of phospholipids [8]. A comparative study of lipid force fields by Piggot and colleagues found that CHARMM36 was the only force field that accurately reproduced the experimental order parameters of carbon 2 in both acyl chains of 1,2-dipalmitoyl-
AMBER force fields were generally less used for membrane protein simulation due to the lack of a specific parameter set for lipids. However, the general AMBER force field (GAFF) which was originally parameterized for the simulation of arbitrary organic molecules with pre-existing AMBER force fields have been shown to reproduce lipid bilayer parameters satisfactorily [10]. GAFF has been tested on a range of lipids, including 1,2-dimyristoyl-d54-
In addition, torsion parameters were also re-optimized using high-level quantum chemical data with “Paramfit.” This strategy has allowed a tensionless simulation of the lipids in the isothermal-isobaric (NPT) ensemble, achieving high level of agreement with experiment data for volume per lipid (VL) and thickness value within 5% of experiment [11]. However, the APL for POPE lipids is lower than expected and thickness value was overestimated by 10% [11]. The force field has also been shown to be able to reproduce the data on large lipid bilayers containing 288 lipids, with little changes in the APL, VL, and bilayer thickness [11]. The comparison with CHARMM27 and Berger force field also displayed GAFF's ability to reproduce the order asymmetry found at the beginning of lipid acyl chains [13]. Prior to the development of CHARMM36, GAFF was the only force field capable of capturing the carbon-2 deuterium order parameter splitting [13], although it has also produced lower surface APL and higher deuterium order parameters compared to experimental data [10, 12, 13]. Even so, compatibility of GAFF with previous AMBER force fields makes it a suitable choice for the simultaneous simulation of membrane, protein, and organic molecule simulation [13].
A recently developed Lipid14 force field with an updated Lipid11 head group and tail group charges and parameters enabled proper tensionless simulation of lipid bilayers in AMBER [14]. Lipid14 LJ and torsion parameters were modified to reproduce the experimental density,
Furthermore, the isothermal area compressibility modulus,
Recent analysis applying full AT force field have reached timescales of up to tens of nanoseconds, with ambitious simulations pushing the microsecond mark, e.g. to characterize the interaction of the multiple sclerosis synthetic biomarker CSF114(Glc) with the membrane bilayer [17], probing of the huntingtin Htt17 membrane anchor on a POPC bilayer [18], mixing of lipids [19], and membrane-binding mechanism of the yeast Osh4 peripheral membrane protein [20]. Even though microsecond simulations have been reported, the accessibility of this approach is limited as many biologically relevant phenomena may require sampling time beyond the microsecond timescale and investigators may not have access to intensive computational resources. Given that lipid reorganization is quite slow, a properly equilibrated membrane may require 20–40 ns of simulation time. Convergence of membrane protein simulations still remains a question, as 100 ns of simulation may still be insufficient to fully describe e.g. rhodopsin loop fluctuations in a membrane [21]. Indeed, since all atomic-level interactions are retained and time-steps for integration of Newtonian motions are in the femtosecond range, AT simulations could be a time consuming and computationally expensive practice. The issue of convergence has been rectified in part by running several long simulations and using non-equilibrium sampling method or steered MD to calculate the observables of interest. A series of extended simulations (~100 ns) on the voltage sensor domain of potassium channels revealed the importance of lipid phosphates in accommodating the significantly charged S4 helix [22]. Using umbrella sampling, potential of mean force has been calculated for permeation and effect of a potassium ion on the second ion passing through gramicidin A channel [23].
Regardless, AT simulations still provided the highest level of detail and reliability when it comes to quantitative prediction of properties such as motional timescales or interaction strengths. Such details recently provided insightful interactions between protein and membrane, such as the rearrangement of amino acid side chains and local bilayer deformation due to hydrophobic mismatch and how the hydrophobic membrane layer accommodates charged arginine side chain of outer membrane phospholipase A [24, 25].
2.2. United-atom (UA) lipid models: the best of both worlds?
The UA representation of lipids simplifies the carbon tails of the lipid by associating the aliphatic carbon and its hydrogen atoms into a single particle. Because the non-polar hydrogen atoms are treated implicitly, the number of interaction sites per lipid can be reduced by two third. The computational costs for simulations of such membrane systems become relatively cheap as the 60% of the pairwise interactions in the membrane is reduced. The model lipid DPPC can be represented by 50 particles in UA force field, but needed 130 interaction sites in an AT force field. Since limited physical information may be collected from explicit acyl chain hydrogens, it became desirable to utilize UA force fields for membrane with an AT protein force field for membrane protein simulation [26, 27].
The UA lipid models parameterized by Berger et al
A CHARMM UA representation of the acyl chains is also available, and compatible for simulations with AT CHARMM protein force fields [26]. The CHARMM UA model was derived from C27 phospholipids, where explicit hydrogen atoms of the acyl chains were replaced with a UA representation. The force field, called C27-UA, retains the accuracy of the AT counterpart and provides a practical alternative when used in simulations of proteins and other compounds described by C27. C27-UA was parameterized by fitting to experimental data and AT simulations of liquid phase model systems (pentadecane for saturated hydrocarbon chains,
Several GROMOS parameters are used in membrane simulations, such as 43A1-S3 [35], 53A6 [36], and Kukol's modification of the 53A6 parameters [37]. 43A1-S3 is an extension and modification of the 43A1 force field designed to improve the properties of lipid membranes in simulations. 43A1-S3 employs charges from Chiu et al. [29] while van der Waals and dihedral parameters were modified from 43A1 to improve hydrocarbon and choline head group dynamics [38]. GROMOS 43A1-S3 accurately reproduces many properties of DPPC bilayers, including APL, lipid diffusion, and the order parameter [39]. The 53A6 parametrization also uses Chiu et al.'s charges [29] and LJ parameters of the choline methyls and phosphate ester oxygen atoms, thus providing good agreement for the APL, density profile, and order parameter for saturated and unsaturated acyl chain PC lipids [40, 41]. Even with improvements introduced, the isothermal area compressibility modulus was still overestimated. Meanwhile, Kukol reparametrize 53A6 following reports that the force field failed to reproduce DPPC satisfactorily [37]. He used Chiu's charges and increased the van der Waals radii of the carbonyl carbons of DPPC, DMPC, POPC, and POPG. He also employed non-standard GROMOS dihedral parameters for the double bond in the unsaturated lipids. This results in fairly good agreement with experimentally determined properties with the exception of order parameters for the
With the advances in computational power, large-scale simulation projects using UA representation of the membrane region have been able to reach the microsecond mark. Access to such timescale allowed the observation of structural rearrangement and changes in hydrogen bonding pattern in an integral Kv1.2 channel embedded in a hyperpolarized membrane [42]. On a shorter timescale, CorA magnesium transport channel has been probed to undergo conformational changes to a putative open state in 110 ns [43].
2.3. Coarse-grained (CG): the need for speed
CG simulations are being widely used to investigate phenomenon occurring in timescales not accessible by AT simulation. In a CG simulation, 3–4 heavy atoms (non-H) are grouped together and represented by a single particle. For example, a DMPC lipid consisting of 130 atoms can be represented by 12 interaction sites [44]. The choline moiety is modeled with a single positively charged particle, the phosphate group with a negatively charged particle, the glycerol linkages with two nonpolar particles, while the lipid chains are modeled with 4 apolar particles each [45].
Early CG approaches were typically parameterized based on comparison to AT simulations by using inverted Monte Carlo schemes [46–48] or force matching [49], which aims to reproduce the structural details at a particular state for a particular system. Instead, the MARTINI (note: MARTINI is a CG force field developed by Marrink and coworkers [44]) CG force field calibrated the building blocks of biomolecules against thermodynamic data, particularly the oil/water partitioning coefficients so that there is no need to re-parameterize the model each time. A consistent, atomic level compatible CG force field will be beneficial for multi-scale applications. In MARTINI, an average of four heavy atoms were represented by a single interaction site, with the exception of ring structures which has 2 or 3 ring atoms mapped to a CG bead. The beads are then distinguished into four main types, i.e., polar, nonpolar, apolar, and charged. MARTINI was able to reproduce the properties of lipid bilayers on a semiquantitative level. These are including the APL, the distribution of groups across the membrane and the bending and area compression moduli.
MARTINI 2.0 further improved the stress profile across the lipid bilayer and its tendency to form pores. Additionally, the free energy of lipid desorption and lipid flip-flop across the bilayer agreed well with AT simulations and the condensing effect of cholesterol on the APL can be reproduced. MARTINI achieved 5- to 10-fold faster sampling of the configurational space of liquid hydrocarbons and the lipid tails inside a bilayer compared to AT force fields [44]. Among others, MARTINI has allowed applications in simulations of vesicle formation and fusion [50, 51], phase transition of lipid bilayers [52], and the structure and dynamics of membrane-protein assemblies [53, 54]. Currently, MARTINI has also been supplemented with an implicit solvent force field, named Dry MARTINI, which provided 1–2 order speed up and is expected to find application in simulations of large membranes containing millions of lipids [55].
The ELBA (an acronym for “electrostatic-based” by Orsi and Essex [56]) force field, so named because it is electrostatics based, offers another alternative for CG simulations of lipids and lipid bilayers. The ELBA model features two approaches in contrast to other CG methods. Notably, LJ interactions are treated using standard Lorentz-Berthelot mixing rules, similar to AT force fields. This is due to the explicit treatment of lipid electrostatics and water dipoles whereby a relative dielectric constant of unity (ɛr = 1) was used to model their interactions. In addition, a realistic dipolar water model based on a simpler soft sticky dipole potential, i.e. the Stockmayer potential is used in ELBA, which helped provide a correct diffusion coefficient of lipids in the liquid phase, a point often overestimated by other CG models. It also shows 15 and 200 times speed improvement over UA and AT models, respectively. Validation has been performed on DOPC, DOPE, and gel phase DSPC, whereby ELBA satisfactorily reproduced several fundamental experimental properties including APL, volume per lipid, curvature elastics constants, electrostatic potential distribution, electrostatic diffusion constant, and lipid diffusion coefficient. At the time of the development, ELBA was only available with the authors' in-house software, and it has since been implemented in LAMMPS [57]. However, there has yet to be a current parameterization for protein and other organic molecules, potentially limiting its use in mixed systems.
CG models for proteins are available and have been adopted into protein-bilayer and peptide-bilayer systems [58, 59]. The advantage in employing CG models is the improved speed and size. As CG force fields have a reduced degree of freedom and removed high-frequency motions such as hydrogen bond vibrations, integration time step could be pushed up to 20–40 fs and highly increases the accessible timescale by 100-fold. Therefore, it is possible to access phenomenon occurring at timescales not accessible to classical simulations, for example, membrane protein aggregation, demonstrated by the formation and domain-specific distribution of Ras proteins in plasma membranes [60]. CG simulation has also been used to self-assemble lipid bilayers around membrane proteins, providing information into protein positioning in the bilayer [61]. Comparison of membrane positioning for six proteins in the CG approach with available experimental data showed high qualitative similarity, even though there are discrepancies from different lipid composition used in the simulation and experiments [61].
Nevertheless, due to simplification of the representations, loss of details is inevitable. Two approaches have been employed to probe for detailed interactions while simulating in CG membrane models is by mixed AT-CG simulation and multiscaling. To elucidate small molecules, membrane peptides and proteins at high resolution, some approaches have used mixed CG-AT systems by simulating the ligand and protein atomistically in a CG membrane [62, 63]. Such mixed descriptions can be likened to quantum mechanical/molecular mechanical (QM/MM) simulations that have achieved considerable success in the past decade (see [64, 65] for reviews regarding usage of QM/MM with proteins). A major consideration in building a mixed AT-CG system is the parameterization and definition of the interactions at the AT-CG interface. As many CG force fields are derived from AT simulations, force-matching procedures can be used to derive effective pairwise CG force field from AT simulations. In AT-CG approach, all atoms simulation of the whole system is first carried out, the system will be subsequently divided into AT and CG parts whereby the most interesting part of the system will retain the atomistic details. The effective AT-CG force field is subsequently obtained by treating the AT and CG parts equally in the force matching procedure [63].
CHARMM-GUI PACE CG Builder [66] was developed for the modeling of large and complex biological system. It utilizes the mixed UA/CG where the PACE force field was used for protein in UA format and Martini CG force field for the water, ions and lipids. Thus, the number of atoms in a system can be reduced by the factor of 10. Analysis showed the PACE/MARTINI hybrid simulation has most of the proteins in the root means square deviation (RMSD) of less than 3 Å.
On the other hand, several methods have been developed for conversion of CG system to AT representation through fragment-based approach [67], simulated annealing [68], and force-matching [49]. These methods, termed the multiscale approach allows the use of CG simulations to explore membrane-lipid interactions, which, after a sufficient equilibration period is converted back to AT model simulation for detailed characterization [67].
3. Tools for membrane simulation setup and analysis
3.1. Automated setup of membrane simulation systems
The availability of many different force fields and parameters for a range of lipid molecules has made it easier to construct systems consisting of mixed lipids. Therefore, the lipid-converter tool is beneficial to easily adapt a system between force fields [69]. The lipid-converter tool can be used in command line by defining a PDB or Gromacs coordinate file and is also available as a web server [70]. The tool currently supports the Berger, GROMOS 43A1-S3, GROMOS 53A6, GROMOS 53A7, CHARMM36, OPLS-UA, and Lipid11. Thus, Stockholm lipids that are compatible with AMBER force field and used CHARMM nomenclature are also supported by this extension. Moreover, lipid converter may be useful to build non-conventional systems as it can generate asymmetric lipid distribution and even label leaflets in curved systems like vesicles.
Towards the end of automating the process of building heterogeneous membrane, the CHARMM-GUI Membrane Builder [71] is a useful tool to generate coordinates for membrane models and protein/membrane systems [72–74]. The membrane builder offers a selection of commonly used lipid models in addition to cholesterol which can be customized according to concentration, APL, hydration number, and thickness of the water layer [73].
Alternatively, a new web server MemGen is able to automatically set up lipid membrane simulation systems without restrictions in force fields, lipid types, or MD simulation software [75]. The user can upload one or more lipid structure files as well as amphiphilic molecules such as alcohol or detergent. A compact representation of each lipid aligned along the
iMembrane is another useful web-based tool which can predict the orientation of a membrane protein within the membrane [76]. Early approaches use a two-state membrane model or a simple hydrophobic slab to model the orientation of a membrane protein in the membrane. Scott and colleagues developed a CG MD to simulate membrane proteins in the presence of membrane lipids which self-assemble into a lipid bilayer [61]. Using the simulation results, iMembrane can predict the orientation of proteins of homologous structure or sequence. BLAST is first performed against the CGDB database for any input sequence or structure. Matches are subsequently realigned to the query using MUSCLE [77] for sequence realignment or “MAMMOTH” [78] for structure super-positioning. Residues in the query are then annotated as N (not in contact with the membrane), H (in contact with polar head group of the membrane lipids), or T (in contact with the lipid hydrophobic tails).
3.2. Tools for membrane MD simulation analysis
As MD simulations for membrane and membrane protein systems became widespread, many groups began developing tools to allow more efficient analysis of the MD trajectories. APL is an important indicator of the membrane phase and stability of the simulation. “GridMAT-MD” is a Perl program which can calculate the APL as well as thickness of a membrane [79]. For bilayer thickness calculation, the user can define a reference atom (such as the lipid phosphate, P atom) and the program first uses the upper leaflet as a reference and assign a paired lipid in the lower leaflet with the upper leaflet based on proximity in the
Mori and colleagues proposed a more sophisticated method for calculating the APL using Voronoi tessellation and Monte Carlo simulation [80]. Coordinates of center of mass for each lipid molecules and coordinates for protein atoms located between the maximum and minimum
Not only that, analysis of other membrane properties became easier with the development of “Membrainy,” an intelligent membrane analysis tool that can provide the calculation of various membrane-specific properties for planar bilayer trajectories [81]. This include APL, order parameter, head group orientation, lipid mixing/demixing entropy, time evolution of the transmembrane voltage, 2D surface map generation, gel percentage, membrane thickness, detection of lipid flip-flop and annular shell lipid analysis. While the program has been primarily designed for use with Gromacs MD package, it is also compatible with pdb trajectory from other MD packages. Currently, it is implemented with CHARMM36, Berger/GROMOS87 and Martini v2.0 force fields, but is also expandable to include other force fields and trajectory formats. Output graphs can be readable by the Grace plotting software.
Other than that, MEMBPLUGIN [82] is another tool to study the MD trajectories of membrane-protein and complex membrane structures. This is a plugin in Visual MD package to measure biophysical properties in the simulated membranes.
4. Towards realistic bilayer simulations
Beyond the improvements in computational power, force field developments, and CG methodologies, more accurate representations of the membrane continued to evolve. The biological membrane is a complex entity composed of numerous lipid species such as phosphatidylcholines, phosphatidylethanolamines, and phosphatidylserine. Other molecules such as cholesterols, sphingomyelins, and cardiolipins also play a role in regulating membrane structure and function. In the outer membrane of many species of Gram-negative bacteria, the presence of lipopolysaccharide in the upperleaflet modulates the insertion, folding, and dynamics of outer membrane proteins within the membrane. Available tools for generate mixed membrane bilayers are including CHARMM-GUI Membrane Builder [73] and MemBuilder [83]. CHARMM-GUI Membrane Builder and MemBuilder supports a total of 32 and 18 different lipids types, respectively.
Straatsma and Soares first reported the simulation of the outer membrane protein OprF in an asymmetric outer membrane with a lipopolysaccharide and phospholipids, describing the saccharides component using GLYCAM parameters [84]. Holdbrook et al
5. Conclusion
The continuous update and improvement of atomistic force fields expanded the types of lipid molecules which could be simulated and increased the accuracy to better match experimental data. Depending on the level of detail of the simulation, UA force fields are an excellent alternative to balance between accuracy and speed. By using CG force fields approaches, sampling and size limitations may be tackled efficiently. Ultimately, advances in computational power and hardware have improved the timescale and system size where MD can be employed. In membrane simulations, the microsecond mark has been reached and simulations are slowly becoming routine work to complement experimental results. In addition, various web server tools and useful analysis programs have been developed to aid membrane simulation analysis. Further advances in lipid force fields will make it possible to characterize membrane structures in greater time and physical scale.
Acknowledgments
This work was funded by the Malaysia Ministry of Higher Education Fundamental Research Grant Scheme (FRGS) (203/CIPPM/6711439) and the Higher Institutions Center of Excellence (HICoE) Grant. S. W. Leong would also like to thank the Malaysia Ministry of Higher Education for MyBrain Science scholarship.
References
- 1.
Singer SJ, Nicolson GL. The fluid mosaic model of the structure of cell membranes. Science. 1972; 175 :720–731. doi: 10.1126/science.175.4023.720 - 2.
van Meer G, Voelker DR, Feigenson GW. Membrane lipids: where they are and how they behave. Nat. Rev. Mol. Cell Biol. 2008; 9 :112–124. doi: 10.1038/nrm2330 - 3.
Lee AG. How lipids affect the activities of integral membrane proteins. Biochim. Biophys. Acta 2004; 1666 :62–87. doi: 10.1016/j.bbamem.2004.05.012 - 4.
Klauda JB, Brooks BR, MacKerell AD, Venable RM, Pastor RW. An ab Initio study on the torsional surface of alkanes and its effect on molecular simulations of alkanes and a DPPC bilayer. J. Phys. Chem. B. 2005;109 :5300–5311. doi: 10.1021/jp0468096 - 5.
Klauda JB, Pastor RW, Brooks BR. Adjacent gauche stabilization in linear alkanes: implications for polymer models and conformational analysis. J. Phys. Chem. B. 2005; 109 :15684–15686. doi: 10.1021/jp0527608 - 6.
Tien HT, Diana AL. Bimolecular lipid membranes: a review and a summary of some recent studies. Chem. Phys. Lipids 1968; 2 :55–101. doi: 10.1016/0009-3084(68)90035-2 - 7.
Marsh D. Renormalization of the tension and area expansion modulus in fluid membranes. Biophys. J. 1997; 73 :865–869. doi: 10.1016/S0006-3495(97)78119-0 - 8.
Klauda JB, Venable RM, Freites JA, O'Connor JW, Tobias DJ, Mondragon-Ramirez C, Vorobyov I, MacKerell AD, Pastor RW, Connor JWO. Update of the CHARMM all-atom additive force field for lipids: validation on six lipid types. J. Phys. Chem. B 2010; 114 :7830–7843. doi: 10.1021/jp101759q - 9.
Piggot TJ, Piñeiro Á, Khalid S. Molecular dynamics simulations of phosphatidylcholine membranes: a comparative force field study. J. Chem. Theory Comput. 2012; 8 :4593–4609. doi: 10.1021/ct3003157 - 10.
Jójárt B, Martinek TA, Jo ZS. Performance of the general AMBER force field in modeling aqueous POPC membrane bilayers. J. Comput. Chem. 2007; 28 :2051–2058. doi: 10.1002/jcc - 11.
Dickson CJ, Rosso L, Betz RM, Walker RC, Gould IR. GAFFlipid: a General AMBER force field for the accurate molecular dynamics simulation of phospholipid. Soft Matter. 2012; 8 :9617–9627. doi: 10.1039/c2sm26007g - 12.
Rosso L, Gould IR. Structure and dynamics of phospholipid bilayers using recently developed general all-atom force fields. J. Comput. Chem. 2007; 29 :24–37. doi: 10.1002/jcc.20675 - 13.
Siu SWI, Vácha R, Jungwirth P, Böckmann R. Biomolecular simulations of membranes: physical properties from different force fields. J. Chem. Phys. 2008; 128 :125103–125103. doi: 10.1063/1.2897760 - 14.
Dickson CJ, Madej BD, Skjevik ÅA, Betz RM, Teigen K, Gould IR, Walker RC. Lipid14: the AMBER lipid force field. J. Chem. Theory Comput. 2014; 10 :865–879. doi: 10.1021/ct4010307 - 15.
Rand RP, Parsegian VA. Hydration forces between phospholipid bilayers. BBA Rev. Biomembranes 1989; 988 :351–376. doi: 10.1016/0304-4157(89)90010-5 - 16.
Rappolt M, Hickel A, Bringezu F, Lohner K. Mechanism of the lamellar/inverse hexagonal phase transition examined by high resolution X-ray diffraction. Biophys. J. 2003; 84 :3111–3122. doi: 10.1016/S0006-3495(03)70036-8 - 17.
Bruno A, Scrima M, Novellino E, D'Errico G, D'Ursi AM, Limongelli V. The glycan role in the glycopeptide immunogenicity revealed by atomistic simulations and spectroscopic experiments on the multiple sclerosis biomarker CSF114(Glc). Sci. Rep. 2015; 5 : 9200. doi: 10.1038/srep09200 - 18.
Côté S, Binette V, Salnikov ES, Bechinger B, Mousseau N. Probing the huntingtin 1-17 membrane anchor on a phospholipid bilayer by using all-atom simulations. Biophys. J. 2015; 108 :1187–1198. doi: 10.1016/j.bpj.2015.02.001 - 19.
Hong C, Tieleman DP, Wang Y. Microsecond molecular dynamics simulations of lipid mixing. Langmuir 2014; 30 :11993–12001. doi: 10.1021/la502363b - 20.
Rogaski B, Klauda JB. Membrane-binding mechanism of a peripheral membrane protein through microsecond molecular dynamics simulations. J. Mol. Biol. 2012; 423 :847–861. doi: 10.1016/j.jmb.2012.08.015 - 21.
Grossfield A, Feller SE, Pitman MC. Convergence of molecular dynamics simulations of membrane proteins. Proteins 2007; 40 :31–40. doi: 10.1002/prot - 22.
Andersson M, Freites JA, Tobias DJ, White SH. Structural dynamics of the S4 voltage-sensor helix in lipid bilayers lacking phosphate groups. J. Phys. Chem. B 2011; 115 :8732–8738. doi: 10.1021/jp2001964 - 23.
Li Y, Andersen OS, Roux B. Energetic of double ion occupancy in the gramicidin A channel. J. Phys. Chem. B 2010; 114 :13881–13888. doi: 10.1021/jp105820u - 24.
Fleming PJ, Freites JA, Moon CP, Tobias DJ, Fleming KG. Outer membrane phospholipase A in phospholipid bilayers: a model system for concerted computational and experimental investigations of amino acid side chain partitioning into lipid bilayers. BBA Biomembranes 2012; 1818 :126–134. doi: 10.1016/j.bbamem.2011.07.016 - 25.
Sonntag Y, Musgaard M, Olesen C, Schiøtt B, Møller JV, Nissen P, Thøgersen L. Mutual adaptation of a membrane protein and its lipid bilayer during conformational changes. Nat. Commun. 2011; 2 :304. doi: 10.1038/ncomms1307 - 26.
Hénin J, Shinoda W, Klein ML. United-aom acyl chains for CHARMM phospholipids. J. Phys. Chem. B 2008; 112 :7008–7015. doi: 10.1021/jp800687p - 27.
Tieleman DP, Maccallum JL, Ash WL, Kandt C, Xu Z, Monticelli L. Membrane protein simulations with a united-atom lipid and all-atom protein model: lipid-protein interactions, side chain transfer free energies and model proteins. J. Phys. Condens. Matter. 2006; 18 :S1221–S1234. doi: 10.1088/0953-8984/18/28/s07 - 28.
Essex JW, Hann MM, Richards WG. Molecular dynamics simulation of a hydrated phospholipid bilayer. Phil. Trans. R. Soc. B 1994; 344 :239–260. doi: 10.2307/55918 - 29.
Chiu SW, Clark M, Balaji V, Subramaniam S, Scott HL, Jakobsson E. Incorporation of surface tension into molecular dynamics simulation of an interface: a fluid phase lipid bilayer membrane. Biophys. J. 1995; 69 :1230–1245. doi: 10.1016/S0006-3495(95)80005-6 - 30.
Berger O, Edholm O, Jähnig F. Molecular dynamics simulations of a fluid bilayer of dipalmitoylphosphatidylcholine at full hydration, constant pressure, and constant temperature. Biophys. J. 1997; 72 :2002–2013. doi: 10.1016/S0006-3495(97)78845-3 - 31.
Cordomí A, Caltabiano G, Pardo L. Membrane protein simulations using AMBER force field and Berger lipid parameters. J. Chem. Theory Comput. 2012; 8 :948–958. doi: 10.1021/ct200491c - 32.
Jerabek H, Pabst G, Rappolt M, Stockner T. Membrane-mediated effect on ion channels induced by the anesthetic drug ketamine. J. Am. Chem. Soc. 2010; 132 :7990–7997. doi: 10.1021/ja910843d - 33.
Lemkul JA, Bevan DR. Characterization of interactions between PilA from Pseudomonas aeruginosa strain K and a model membrane. J. Phys. Chem. B 2011;115 :8004–8008. doi: 10.1021/jp202217f - 34.
Lensink MF, Govaerts C, Ruysschaert J-M. Identification of specific lipid-binding sites in integral membrane proteins. J. Biol. Chem. 2010; 285 :10519–10526. doi: 10.1074/jbc.M109.068890 - 35.
Braun AR, Sachs JN, Nagle JF. Comparing simulations of lipid bilayers to scattering data: the GROMOS 43A1-S3 force field. J. Phys. Chem. B 2013; 117 :5065–5072. doi: 10.1021/jp401718k - 36.
Oostenbrink C, Villa A, Mark AE, Van Gunsteren WF. A biomolecular force field based on the free enthalpy of hydration and solvation: the GROMOS force-field parameter sets 53A5 and 53A6. J. Comput. Chem. 2004; 25 :1656–1676. doi: 10.1002/jcc.20090 - 37.
Kukol A. Lipid models for united-atom molecular dynamics simulations of proteins. J. Chem. Theory Comput. 2009; 5 :615–626. doi: 10.1021/ct8003468 - 38.
Scott WRP, Hünenberger PH, Tironi IG, Mark AE, Billeter SR, Fennen J, Torda AE, Huber T, Krüger P, van Gunsteren WF. The GROMOS biomolecular simulation program package. J. Phys. Chem. A 1999; 103 :3596–3607. doi: 10.1021/jp984217f - 39.
Chiu SW, Pandit SA, Scott HL, Jakobsson E. An improved united atom force field for simulation of mixed lipid bilayers. J. Phys. Chem. B 2009; 113 :2748–2763. doi: 10.1021/jp807056c - 40.
Poger D, Van Gunsteren WF, Mark AE. A new force field for simulating phosphatidylcholine bilayers. J. Comput. Chem. 2010; 31 :1117–1125. doi: 10.1002/jcc.21396 - 41.
Poger D, Mark AE. On the validation of molecular dynamics simulations of saturated and cis-monounsaturated phosphatidylcholine lipid bilayers: a comparison with experiment. J. Chem. Theory Comput. 2010; 6 :325–336. doi: 10.1021/ct900487a - 42.
Bjelkmar P, Niemelä PS, Vattulainen I, Lindahl E. Conformational changes and slow dynamics through microsecond polarized atomistic molecular simulation of an integral Kv1.2 ion channel. PLoS Comput. Biol. 2009; 5 :e1000289. doi: 10.1371/journal.pcbi.1000289 - 43.
Chakrabarti N, Neale C, Payandeh J, Pai EF, Pomès R. An iris-like mechanism of pore dilation in the CorA magnesium transport system. Biophys. J. 2010; 98 :784–792. doi: 10.1016/j.bpj.2009.11.009 - 44.
Marrink SJ, Risselada HJ, Yefimov S, Tieleman DP, de Vries AH. The MARTINI force field: coarse grained model for biomolecular simulations. J. Phys. Chem. B 2007; 111 :7812–7824. doi: 10.1021/jp071097f - 45.
Marrink SJ, de Vries AH, Mark AE. Coarse grained model for semiquantitative lipid simulations. J. Phys. Chem. B 2004; 108 :750–760. doi: 10.1021/jp036508g - 46.
Shelley JC, Shelley MY, Reeder RC, Bandyopadhyay S, Klein ML. A coarse grain model for phospholipid simulations. J. Phys. Chem. B 2001; 105 :4464–4470. doi: 10.1021/jp010238p - 47.
Lyubartsev AP. Multiscale modeling of lipids and lipid bilayers. Eur. Biophys. J. 2005; 35 :53–61. doi: 10.1007/s00249-005-0005-y - 48.
Elezgaray J, Laguerre M. A systematic method to derive force fields for coarse-grained simulations of phospholipids. Comput. Phys. Commun. 2006; 175 :264–268. doi: 10.1016/j.cpc.2006.01.009 - 49.
Izvekov S, Voth GA. A multiscale coarse-graining method for biomolecular systems. J. Phys. Chem. B 2005; 109 :2469–2473. doi: 10.1021/jp044629q - 50.
Marrink SJ, Mark AE. The mechanism of vesicle fusion as revealed by molecular dynamics simulations. J. Am. Chem. Soc. 2003; 125 :11144–11145. doi: 10.1021/ja036138+ - 51.
Marrink SJ, Mark AE. Molecular dynamics simulation of the formation, structure, and dynamics of small phospholipid vesicles. J. Am. Chem. Soc. 2003; 125 :15233–15242. doi: 10.1021/ja0352092 - 52.
Prates Ramalho JP, Gkeka P, Sarkisov L. Structure and phase transformations of DPPC lipid bilayers in the presence of nanoparticles: insights from coarse-grained molecular dynamics simulations. Langmuir 2011; 27 :3723–3730. doi: 10.1021/la200236d - 53.
Sansom MSP, Scott KA, Bond PJ. Coarse-grained simulation: a high-throughput computational approach to membrane proteins. Biochem. Soc. Trans. 2008; 36 :27–32. doi: 10.1042/bst0360027 - 54.
Bucher D, Hsu Y-H, Mouchlis VD, Dennis EA, McCammon JA. Insertion of the Ca2+-independent phospholipase A2 into a phospholipid bilayer via coarse-grained and atomistic molecular dynamics simulations. PLoS Comput. Biol. 2013; 9 :e1003156. doi: 10.1371/journal.pcbi.1003156 - 55.
Arnarez C, Uusitalo JJ, Masman MF, Ingólfsson HI, de Jong DH, Melo MN, Periole X, de Vries AH, Marrink SJ. Dry Martini, a coarse-grained force field for lipid membrane simulations with implicit solvent. J. Chem. Theory Comput. 2015; 11 :260–275. doi: 10.1021/ct500477k - 56.
Orsi M, Essex JW. The ELBA force field for coarse-grain modeling of lipid membranes. PLoS One 2011; 6 :e28637. doi: 10.1371/journal.pone.0028637 - 57.
Plimpton S. Fast parallel algorithms for short-range molecular dynamics. J. Comput. Phys. 1995; 117 :1–19. doi: 10.1006/jcph.1995.1039 - 58.
Monticelli L, Kandasamy SK, Periole X, Larson RG, Tieleman DP, Marrink SJ. The MARTINI coarse-grained force field: extension to proteins. J. Chem. Theory Comput. 2008; 4 :819–834. doi: 10.1021/ct700324x - 59.
Spijker P, van Hoof B, Debertrand M, Markvoort AJ, Vaidehi N, Hilbers PJ. Coarse grained molecular dynamics simulations of transmembrane protein-lipid systems. Int. J. Mol. Sci. 2010; 11 :2393–2420. doi: 10.3390/ijms11062393 - 60.
Janosi L, Li Z, Hancock JF, Gorfe AA. Organization, dynamics, and segregation of Ras nanoclusters in membrane domains. Proc. Natl. Acad. Sci. USA 2012; 109 :8097–8102. doi: 10.1073/pnas.1200773109 - 61.
Scott KA, Bond PJ, Ivetac A, Chetwynd AP, Khalid S, Sansom MSP. Coarse-grained MD simulations of membrane protein-bilayer self-assembly. Structure 2008; 16 :621–30. doi: 10.1016/j.str.2008.01.014 - 62.
Orsi M, Essex JW. Permeability of drugs and hormones through a lipid bilayer: insights from dual-resolution molecular dynamics. Soft Matter. 2010; 6 :3797–3797. doi: 10.1039/c0sm00136h - 63.
Shi Q, Izvekov S, Voth G. Mixed atomistic and coarse-grained molecular dynamics: simulation of a membrane-bound ion channel. J. Phys. Chem. B 2006; 110 :15045–15048. doi: 10.1021/jp062700h - 64.
Murphy RB, Philipp DM, Friesner RA. A mixed quantum mechanics/molecular mechanics (QM/MM) method for large-scale modeling of chemistry in protein environments. J. Comput. Chem. 2000; 21 :1442–1457. doi: 10.1002/1096-987X(200012)21:16<1442::AID-JCC3>3.0.CO;2-O - 65.
Senn HM, Thiel W. QM/MM studies of enzymes. Curr. Op. Chem. Biol. 2007; 11 :182–187. doi: 10.1016/j.cbpa.2007.01.684 - 66.
Qi Y, Cheng X, Han W, Jo S, Schulten K, Im W. CHARMM-GUI PACE CG Builder for solution, micelle, and bilayer coarse-grained simulations. J. Chem. Inf. Model. 2014; 54 :1003–1009. doi: 10.1021/ci500007n - 67.
Stansfeld PJ, Sansom MSPP. From coarse grained to atomistic: a serial multiscale approach to membrane protein simulations. J. Chem. Theory Comput. 2011; 7 :1157–1166. doi: 10.1021/ct100569y - 68.
Rzepiela AJ, Schäfer LV, Goga N, Risselada HJ, De Vries AH, Marrink SJ. Reconstruction of atomistic details from coarse-grained structures. J. Comput. Chem. 2010; 31 :1333–1343. doi: 10.1002/jcc.21415 - 69.
Larsson P, Kasson PM. Lipid converter, a framework for lipid manipulations in molecular dynamics simulations. J. Membr. Biol. 2014; 247 :1137–1140. doi: 10.1007/s00232-014-9705-5 - 70.
Larsson P, Kasson P. Lipid-converter [Internet]. 2014. Available from: http://lipid-converter.appspot.com/ [accessed: 2016-02-16]. - 71.
Jo S, Lim JB, Klauda JB, Im W. CHARMM-GUI Membrane Builder [Internet]. 2009. Available from: http://www.charmm-gui.org/input/membrane [Accessed: 2016-02-16] - 72.
Jo S, Kim T, Im W. Automated builder and database of protein/membrane complexes for molecular dynamics simulations. PLoS One 2007; 2 :e880. doi: 10.1371/journal.pone.0000880 - 73.
Jo S, Lim JB, Klauda JB, Im W. CHARMM-GUI membrane builder for mixed bilayers and its application to yeast membranes. Biophys. J. 2009; 97 :50–58. doi: 10.1016/j.bpj.2009.04.013 - 74.
Wu EL, Cheng X, Jo S, Rui H, Song KC, Davila-Contreras EM, Qi Y, Lee J, Monje-Galvan V, Venable RM, Klauda JB, Im W. CHARMM-GUI membrane builder toward realistic biological membrane simulations. J. Comput. Chem. 2014; 35 :1997–2004. doi: 10.1002/jcc.23702 - 75.
Knight CJ, Hub JS. MemGen: a general web server for the setup of lipid membrane simulation systems. Bioinformatics 2015; 31 :2897–2899. doi: 10.1093/bioinformatics/btv292 - 76.
Kelm S, Shi J, Deane CM. iMembrane: homology-based membrane-insertion of proteins. Bioinformatics 2009; 25 :1086–1088. doi: 10.1093/bioinformatics/btp102 - 77.
Edgar RC. MUSCLE: multiple sequence alignment with high accuracy and high throughput. Nucleic Acids Res. 2004; 32 :1792–1797. doi: 10.1093/nar/gkh340 - 78.
Ortiz AR, Strauss CEM, Olmea O. MAMMOTH (matching molecular models obtained from theory): an automated method for model comparison. Protein Sci. 2002; 11 :2606–2621. doi: 10.1110/ps.0215902 - 79.
Allen WJ, Lemkul JA, Bevan DR. GridMAT-MD: a grid-based membrane analysis tool for use with molecular dynamics. J. Comput. Chem. 2009; 30 :1952–1958. doi: 10.1002/jcc.21172 - 80.
Mori T, Ogushi F, Sugita Y. Analysis of lipid surface area in protein–membrane systems combining voronoi tessellation and monte carlo integration methods. J. Comput. Chem. 2012; 33 :286–293. doi: 10.1002/jcc.21973 - 81.
Carr M, MacPhee CE. Membrainy: a ‘smart’, unified membrane analysis tool. Source Code Biol. Med. 2015; 10 :1–10. doi: 10.1186/s13029-015-0033-7 - 82.
Guixa-Gonzalez R, Rodriguez-Espigares I, Ramirez-Anguita JM, Carrio-Gaspar P, Martinez-Seara H, Giorgino T, Selent J. MEMBPLUGIN: studying membrane complexity in VMD. Bioinformatics 2014; 30 :1478–1480. doi: 10.1093/bioinformatics/btu037 - 83.
Ghahremanpour MM, Arab SS, Aghazadeh SB, Zhang J, van der Spoel D. MemBuilder: a web-based graphical interface to build heterogeneously mixed membrane bilayers for the GROMACS biomolecular simulation program. Bioinformatics 2013; 30 :439–441. doi: 10.1093/bioinformatics/btt680 - 84.
Straatsma TP, Soares TA. Characterization of the outer membrane protein OprF of Pseudomonas aeruginosa in a lipopolysaccharide membrane by computer simulation. Proteins 2009;74 :475–488. doi: 10.1002/prot.22165 - 85.
Holdbrook DA, Piggot TJ, Sansom MSP, Khalid S. Stability and membrane interactions of an autotransport protein: MD simulations of the Hia translocator domain in a complex membrane environment. Biochim. Biophys. Acta 2013; 1828 :715–723. doi: 10.1016/j.bbamem.2012.09.002