Dataset Information

Acceleration of Lateral Equilibration in Mixed Lipid Bilayers Using Replica Exchange with Solute Tempering.

ABSTRACT: The lateral heterogeneity of cellular membranes plays an important role in many biological functions such as signaling and regulating membrane proteins. This heterogeneity can result from preferential interactions between membrane components or interactions with membrane proteins. One major difficulty in molecular dynamics simulations aimed at studying the membrane heterogeneity is that lipids diffuse slowly and collectively in bilayers, and therefore, it is difficult to reach equilibrium in lateral organization in bilayer mixtures. Here, we propose the use of the replica exchange with solute tempering (REST) approach to accelerate lateral relaxation in heterogeneous bilayers. REST is based on the replica exchange method but tempers only the solute, leaving the temperature of the solvent fixed. Since the number of replicas in REST scales approximately only with the degrees of freedom in the solute, REST enables us to enhance the configuration sampling of lipid bilayers with fewer replicas, in comparison with the temperature replica exchange molecular dynamics simulation (T-REMD) where the number of replicas scales with the degrees of freedom of the entire system. We apply the REST method to a cholesterol and 1,2-dipalmitoyl-sn-glycero-3-phosphocholine (DPPC) bilayer mixture and find that the lateral distribution functions of all molecular pair types converge much faster than in the standard MD simulation. The relative diffusion rate between molecules in REST is, on average, an order of magnitude faster than in the standard MD simulation. Although REST was initially proposed to study protein folding and its efficiency in protein folding is still under debate, we find a unique application of REST to accelerate lateral equilibration in mixed lipid membranes and suggest a promising way to probe membrane lateral heterogeneity through molecular dynamics simulation.

SUBMITTER: Huang K

PROVIDER: S-EPMC4196747 | biostudies-literature | 2014 Oct

REPOSITORIES: biostudies-literature

ACCESS DATA

Publications

Acceleration of Lateral Equilibration in Mixed Lipid Bilayers Using Replica Exchange with Solute Tempering.

Huang Kun K García Angel E AE

Journal of chemical theory and computation 20140912 10

The lateral heterogeneity of cellular membranes plays an important role in many biological functions such as signaling and regulating membrane proteins. This heterogeneity can result from preferential interactions between membrane components or interactions with membrane proteins. One major difficulty in molecular dynamics simulations aimed at studying the membrane heterogeneity is that lipids diffuse slowly and collectively in bilayers, and therefore, it is difficult to reach equilibrium in lat ...[more]

PMID: 25328493

Similar Datasets

Project description:Time-resolved fluorescence and differential scanning calorimetry (DSC) were used to examine how two amino acids, L-phenylalanine (L-PA) and N-acetyl-DL-tryptophan (NAT), affect the temperature-dependent membrane affinity of two structurally similar coumarin solutes for 1,2-dipalmitoyl-sn-glycero-3-phosphocholine (DPPC) vesicles. The 7-aminocoumarin solutes, coumarin 151 (C151) and coumarin 152 (C152), differ in their substitution at amine position-C151 is a primary amine, and C152 is a tertiary amine-and both solutes show different tendencies to associate with lipid bilayers consistent with differences in their respective log-P-values. Adding L-PA to the DPPC vesicle solution did not change C151's propensity to remain freely solvated in aqueous solution, but C152 showed a greater tendency to partition into the hydrophobic bilayer interior at temperatures below DPPC's gel-liquid crystalline transition temperature (Tgel-lc). This finding is consistent with L-PA's ability to enhance membrane permeability by disrupting chain-chain interactions. Adding NAT to DPPC-vesicle-containing solutions changed C151 and C152 affinity for the DPPC membranes in unexpected ways. DSC data show that NAT interacts strongly with the lipid bilayer, lowering Tgel-lc by up to 2°C at concentrations of 10 mM. These effects disappear when either C151 or C152 is added to solution at concentrations below 10 μM, and Tgel-lc returns to a value consistent with unperturbed DPPC bilayers. Together with DSC results, fluorescence data imply that NAT promotes coumarin adsorption to the vesicle bilayer surface. NAT's effects diminish above Tgel-lc and imply that unlike L-PA, NAT does not penetrate into the bilayer but instead remains adsorbed to the bilayer's exterior. Taken in their entirety, these discoveries suggest that amino acids-and by inference, polypeptides and proteins-change solute affinity for lipid bilayers with specific effects that depend on individualized amino-acid-lipid-bilayer interactions.

Project description:Polarizable force fields for lipid and solvent environments are used for molecular dynamics simulations of a fully hydrated dipalmitoylphosphatidylcholine (DPPC) bilayer and gramicidin A (gA) dimer embedded in a dimyristoylphosphatidylcholine (DMPC) bilayer. The lipid bilayer is modeled using the CHARMM charge equilibration (CHEQ) polarizable force field for lipids and the TIP4P-FQ force field to represent solvent. For the DPPC bilayer system, results are compared to the same system simulated using the nonpolarizable CHARMM27r (C27r) force field and TIP3P water. Calculated atomic and electron density profiles, head group orientations as measured by the phosphorus-nitrogen vector orientation, and deuterium order parameters are found to be consistent with previous simulations and with experiment. The CHEQ model exhibits greater water penetration into the bilayer interior, as demonstrated by the potential of mean force calculated from the water density profile. This is a result of the variation of the water molecular dipole from 2.55 D in the bulk to 1.88 D in the interior. We discuss this finding in the context of previous studies (both simulation and experiment) that have investigated the extent of penetration of water into DPPC bilayers. We also discuss the effects of including explicit polarization on the water dipole moment variation as a function of distance from the bilayer. We show distributions of atomic charges over the course of the simulation since the CHEQ model allows the charges to fluctuate. We have calculated the interfacial dipole potential, which the CHEQ model predicts to be 0.95 V compared to 0.86 V as predicted by the C27r model. We also discuss dielectric permittivity profiles and the differences arising between the two models. We obtain bulk values of 72.77 for the CHEQ model (TIP4P-FQ water) and 91.22 for C27r (TIP3P), and values approaching unity in the membrane interior. Finally, we present results of simulations of gA embedded in a DMPC bilayer using the CHEQ model and discuss structural properties.

Dataset Information

Acceleration of Lateral Equilibration in Mixed Lipid Bilayers Using Replica Exchange with Solute Tempering.

Publications

Acceleration of Lateral Equilibration in Mixed Lipid Bilayers Using Replica Exchange with Solute Tempering.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets