Indira Chandrasekhar

A consistent potential energy parameter set for lipids: dipalmitoylphosphatidylcholine as a benchmark of the GROMOS96 45A3 force field

The performance of the GROMOS96 parameter set 45A3 developed for aliphatic alkanes is tested on a bilayer of dipalmitoylphosphatidylcholine (DPPC) in water in the liquid-crystalline Lα phase. Variants of the force-field parameter set as well as different sets of simulation conditions or simulation parameter sets are evaluated. In the case of the force-field parameters, the van der Waals constants for the non-bonded interaction of the ester carbonyl carbon and the partial charges and charge group definition of the phosphatidylcholine head group are examined. On the methodological side, different cut-off distances for the non-bonded interactions, use of a reaction-field force due to long-range electrostatic interactions, the frequency of removal of the centre of mass motion and the strength of the coupling of the pressure of the system to the pressure bath are tested. The area per lipid, as a measure of structure, the order parameters of the chain carbons, as a measure of membrane fluidity, and the translational diffusion of the lipids in the plane of the bilayer are calculated and compared with experimental values. An optimal set of simulation parameters for which the GROMOS96 parameter set 45A3 yields a head group area, chain order parameters and a lateral diffusion coefficient in accordance with the experimental data is listed.

Interaction of the Disaccharide Trehalose with a Phospholipid Bilayer: A Molecular Dynamics Study

The disaccharide trehalose is well known for its bioprotective properties. Produced in large amounts during stress periods in the life of organisms able to survive potentially damaging conditions, trehalose plays its protective role by stabilizing biostructures such as proteins and lipid membranes. In this study, molecular dynamics simulations are used to investigate the interaction of trehalose with a phospholipid bilayer at atomistic resolution. Simulations of the bilayer in the absence and in the presence of trehalose at two different concentrations (1 or 2 molal) are carried out at 325 K and 475 K. The results show that trehalose is able to minimize the disruptive effect of the elevated temperature and stabilize the bilayer structure. At both temperature, trehalose is found to interact directly with the bilayer through hydrogen bonds. However, the water molecules at the bilayer surface are not completely replaced. At high temperature, the protective effect of trehalose is correlated with a significant increase in the number of trehalose-bilayer hydrogen bonds, predominantly through an increase in the number of trehalose molecules bridging three or more lipid molecules.

A molecular dynamics study of the bee venom melittin in aqueous solution, in methanol, and inserted in a phospholipid bilayer

The structural properties of melittin, a small amphipathic peptide found in the bee venom, are investigated in three different environments by molecular dynamics simulation. Long simulations have been performed for monomeric melittin solvated in water, in methanol, and shorter ones for melittin inserted in a dimyristoylphosphatidylcholine bilayer. The resulting trajectories were analysed in terms of structural properties of the peptide and compared to the available NMR data. While in water and methanol solution melittin is observed to partly unfold, the peptide retains its structure when embedded in a lipid bilayer. The latter simulation shows good agreement with the experimentally derived 3J-coupling constants. Generally, it appears that higher the stability of the helical conformation of melittin, lower is the dielectric permittivity of the environment. In addition, peptide-lipid interactions were investigated showing that the C-terminus of the peptide provides an anchor to the lipid bilayer by forming hydrogen bonds with the lipid head groups.

A comparison of the potential energy parameters of aliphatic alkanes: molecular dynamics simulations of triacylglycerols in the alpha phase

Abstract. A comparison of six GROMOS parameter sets for aliphatic alkanes is presented using simulations of a triglyceride, trioctanoin, in the gel or α phase. It is found that the parameter set 43A2 results in a system that is rigid and close packed, forcing a collective tilt of the molecules, known experimentally to be disallowed. The parameter set 45A3-45 results in an expanded system where the molecules undergo a conformational change from the chair to tuning-fork form. The lamellae are fragmented, with the fragments shifted in the direction parallel to the layer normal. The newest parameter set 45A3, when used with either of two sets of partial charges for the glycerol ester group, retains the characteristics of the gel phase. The last parameter set tested, 45A3-45×12, also performs well when combined with either set of partial charges, although the system expands in the plane of the bilayer.

Simulating the Physiological Phase of Hydrated DPPC Bilayers: The Ester Moiety

Abstract The charge on the ester oxygen of the sn2 group of the dipalmitoylphosphatidylcholine (DPPC) has a remarkable effect on the square area per lipid in simulations of a hydrated bilayer. This is in contrast to simulations of nonpolar, neutral lipids, where it has been found to have little effect. The charges associated with the GROMOS96 45A3 and 45A4 biomolecular force fields have been previously shown to cause significant membrane shrinkage. We find that the use of larger charges at the ester groups alone (as opposed to on all the polar moieties in the head group) remedies the shrinkage. The source of this effect in DPPC lies in the fact that the charge distribution of this polar group profoundly influences its free energy of hydration and, correspondingly, the water distribution around it. In an attempt to rationally tune the ester parameters, the repulsive Lennard–Jones parameters that represent the van der Waals interaction have been refined to reproduce the experimental density and heat of vaporization, and the charges of the ester groups have been tuned to reproduce the experimental free energies of hydration of a series of alkane esters. The new parameters form part of the GROMOS96 53A5 and 53A6 force fields. However, with the new force‐field parameters, the area per lipid in simulations of hydrated DPPC bilayers lies below that of the physiological liquid‐crystalline phase, the implications of which are discussed.

Molecular dynamics simulation of lipid bilayers with GROMOS96: Application of surface tension

The GROMOS96 force fields 45A3 and 53A5, when applied to dipalmitoylphosphatidylcholine (DPPC) membranes, have a tendency to result in a reduced area per lipid in constant pressure simulations. The application of surface tension is effective in increasing the area per lipid, a measure of the phase of the membrane, but only if the area is already close to the experimental range. Therefore the surface tension cannot compensate for strong inadequacies in the force-field parameters. The behaviour of the 45A3 force field from long NP n γT simulations of tens of nanoseconds is analysed over a range of different surface tensions. Comparisons are made with the corresponding NP n AT simulations.