Biomedical subjects

Siewert J Marrink

Publications and source records attributed to Siewert J Marrink.

8 recordsLinked to original sources

PubMed2006

Does isoprene protect plant membranes from thermal shock? A molecular dynamics study.

The question of why plants release isoprene when heat stressed has been continuously debated for more than half a century. In this work we use molecular dynamics simulation techniques to directly investigate the interaction between isoprene and a model phospholipid membrane in atomic detail. It is found that isoprene partitions preferentially in the center of the membrane and in a dose dependent manner enhances the order within the membrane without significantly changing the dynamical properties of the system. At a concentration of 20 mol% isoprene (16 isoprene molecules per 64 lipid molecules) the effect of the addition of isoprene on the membrane order is equivalent to a reduction in temperature of 10 K, rising to a reduction of 30 K at 43 mol% isoprene. The significance of the work is that it provides for the first time direct evidence that isoprene stabilizes lipid membranes and reduces the likelihood of a phospholipid membrane undergoing a heat induced phase transition. Furthermore it provides a clear mechanistic picture as to why plants specifically utilize isoprene for this purpose.

Butadienes↗

PubMed2006

Antimicrobial peptides in action.

Molecular dynamics simulations of the magainin MG-H2 peptide interacting with a model phospholipid membrane have been used to investigate the mechanism by which antimicrobial peptides act. Multiple copies of the peptide were randomly placed in solution close to the membrane. The peptide readily bound to the membrane, and above a certain concentration, the peptide was observed to cooperatively induce the formation of a nanometer-sized, toroidally shaped pore in the bilayer. In sharp contrast with the commonly accepted model of a toroidal pore, only one peptide was typically found near the center of the pore. The remaining peptides lay close to the edge of the pore, maintaining a predominantly parallel orientation with respect to the membrane.

1,2-Dipalmitoylphosphatidylcholine↗

PubMed2005

Molecular structure of the lecithin ripple phase.

Molecular dynamics simulations of lecithin lipid bilayers in water as they are cooled from the liquid crystalline phase show the spontaneous formation of rippled bilayers. The ripple consists of two domains of different length and orientation, connected by a kink. The organization of the lipids in one domain of the ripple is found to be that of a splayed gel; in the other domain the lipids are gel-like and fully interdigitated. In the concave part of the kink region between the domains the lipids are disordered. The results are consistent with the experimental information available and provide an atomic-level model that may be tested by further experiments.

1,2-Dipalmitoylphosphatidylcholine↗

PubMed2005

Simulation of gel phase formation and melting in lipid bilayers using a coarse grained model.

The transformation between a gel and a fluid phase in dipalmitoyl-phosphatidylcholine (DPPC) bilayers has been simulated using a coarse grained (CG) model by cooling bilayer patches composed of up to 8000 lipids. The critical step in the transformation process is the nucleation of a gel cluster consisting of 20-80 lipids, spanning both monolayers. After the formation of the critical cluster, a fast growth regime is entered. Growth slows when multiple gel domains start interacting, forming a percolating network. Long-lived fluid domains remain trapped and can be metastable on a microsecond time scale. From the temperature dependence of the rate of cluster growth, the line tension of the fluid-gel interface was estimated to be 3+/-2 pN. The reverse process is observed when heating the gel phase. No evidence is found for a hexatic phase as an intermediate stage of melting. The hysteresis observed in the freezing and melting transformation is found to depend both on the system size and on the time scale of the simulation. Extrapolating to macroscopic length and time scales, the transition temperature for heating and cooling converges to 295+/-5 K, in semi-quantitative agreement with the experimental value for DPPC (315 K). The phase transformation is associated with a drop in lateral mobility of the lipids by two orders of magnitude, and an increase in the rotational correlation time of the same order of magnitude. The lipid headgroups, however, remain fluid. These observations are in agreement with experimental findings, and show that the nature of the ordered phase obtained with the CG model is indeed a gel rather than a crystalline phase. Simulations performed at different levels of hydration furthermore show that the gel phase is stabilized at low hydration. A simulation of a small DPPC vesicle reveals that curvature has the opposite effect.

1,2-Dipalmitoylphosphatidylcholine↗

PubMed2004

Molecular dynamics simulation of the spontaneous formation of a small DPPC vesicle in water in atomistic detail.

Molecular dynamics simulations have been used to study the spontaneous aggregation of a concentrated solution of dipalmitoylphosphatidylcholine (DPPC) molecules in water into a small vesicle. The molecules were represented in atomistic detail. Starting from a DPPC solution in water, an oblong vesicle with a long axis of 15 nm and short axes of 10 nm was formed spontaneously. After 90 ns of simulation, the vesicle contained a number of water pores. Water pores were shown to facilitate exchange of lipids between inner and outer leaflets. Lipid tails were shown to be less ordered in the inner leaflet of the vesicle, as compared to those in the outer leaflet of the vesicle and an equilibrated lamellar bilayer.

1,2-Dipalmitoylphosphatidylcholine↗

PubMed2004

Molecular dynamics simulations of hydrophilic pores in lipid bilayers.

Hydrophilic pores are formed in peptide free lipid bilayers under mechanical stress. It has been proposed that the transport of ionic species across such membranes is largely determined by the existence of such meta-stable hydrophilic pores. To study the properties of these structures and understand the mechanism by which pore expansion leads to membrane rupture, a series of molecular dynamics simulations of a dipalmitoylphosphatidylcholine (DPPC) bilayer have been conducted. The system was simulated in two different states; first, as a bilayer containing a meta-stable pore and second, as an equilibrated bilayer without a pore. Surface tension in both cases was applied to study the formation and stability of hydrophilic pores inside the bilayers. It is observed that below a critical threshold tension of approximately 38 mN/m the pores are stabilized. The minimum radius at which a pore can be stabilized is 0.7 nm. Based on the critical threshold tension the line tension of the bilayer was estimated to be approximately 3 x 10(-11) N, in good agreement with experimental measurements. The flux of water molecules through these stabilized pores was analyzed, and the structure and size of the pores characterized. When the lateral pressure exceeds the threshold tension, the pores become unstable and start to expand causing the rupture of the membrane. In the simulations the mechanical threshold tension necessary to cause rupture of the membrane on a nanosecond timescale is much higher in the case of the equilibrated bilayers, as compared with membranes containing preexisting pores.

1,2-Dipalmitoylphosphatidylcholine↗

PubMed2003

Molecular dynamics simulation of the formation, structure, and dynamics of small phospholipid vesicles.

Here, we use coarse grained molecular dynamics (MD) simulations to study the spontaneous aggregation of dipalmitoylphosphatidylcholine (DPPC) lipids into small unilamellar vesicles. We show that the aggregation process occurs on a nanosecond time scale, with bicelles and cuplike vesicles formed at intermediate stages. Formation of hemifused vesicles is also observed at higher lipid concentration. With either 25% dipalmitoylphosphatidylethanolamine (DPPE) or lysoPC mixed into the system, the final stages of the aggregation process occur significantly faster. The structure of the spontaneously formed vesicles is analyzed in detail. Microsecond simulations of isolated vesicles reveal significant differences in the packing of the lipids between the inner and outer monolayers, and between PC, PE, and lysoPC. Due to the small size of the vesicles they remain almost perfectly spherical, undergoing very limited shape fluctuations or bilayer undulations. The lipid lateral diffusion rate is found to be faster in the outer than in the inner monolayer. The water permeability coefficient of the pure DPPC vesicles is of the order of 10(-)(3) cm s(-)(1), in agreement with experimental measurements.

1,2-Dipalmitoylphosphatidylcholine↗

PubMed2003

The mechanism of vesicle fusion as revealed by molecular dynamics simulations.

We describe molecular dynamics simulations elucidating the molecular details of the process of fusion for small lipid vesicles. The simulations are based on a coarse grained (CG) lipid model that accurately represents the lamellar state of a variety of phospholipids and enables us to observe intermediate stages during fusion at near atomic detail. Simulations were conducted on a variety of systems containing common phospholipids such as phosphatidylcholine (PC), phosphatidylethanolamine (PE), lysoPC, and mixtures of the above. The fusion intermediates found are in general agreement with the stalk-pore mechanism. Transient pores sometimes form adjacent to the stalk, however, resulting in the mixing of lipids from the outer and inner monolayers. The speed of stalk formation and the opening of the fusion pore can be modulated by altering the lipid composition in qualitative agreement with experimental observations.

Cytoplasmic Vesicles↗