PubMed Health⌕ Search

Biomedical subjects

J C Shillcock

Publications and source records attributed to J C Shillcock.

5 recordsLinked to original sources

Two-component membrane material properties and domain formation from dissipative particle dynamics.

The material parameters (area stretch modulus and bending rigidity) of two-component amphiphilic membranes are determined from dissipative particle dynamics simulations. The preferred area per molecule for each species is varied so as to produce homogeneous mixtures or nonhomogeneous mixtures that form domains. If the latter mixtures are composed of amphiphiles with the same tail length, but different preferred areas per molecule, their material parameters increase monotonically as a function of composition. By contrast, mixtures of amphiphiles that differ in both tail length and preferred area per molecule form both homogeneous and nonhomogeneous mixtures that both exhibit smaller values of their material properties compared to the corresponding pure systems. When the same nonhomogeneous mixtures of amphiphiles are assembled into planar membrane patches and vesicles, the resulting domain shapes are different when the bending rigidities of the domains are sufficiently different. Additionally, both bilayer and monolayer domains are observed in vesicles. We conclude that the evolution of the domain shapes is influenced by the high curvature of the vesicles in the simulation, a result that may be relevant for biological vesicle membranes.

Computer Simulation↗

Effect of chain length and asymmetry on material properties of bilayer membranes.

Dissipative particle dynamics is used to extract the material parameters (bending and area stretch moduli) of a bilayer membrane patch. Some experiments indicate that the area stretch modulus of lipid vesicles varies little as the chain length of the lipids composing the bilayer increases. Here we show that making the interactions between the hydrophilic head groups of the model amphiphiles proportional to the hydrophobic tail length reproduces the above result for the area stretch modulus. We also show that the area stretch modulus of bilayers composed of amphiphiles with the same number of tail beads but with asymmetric chains is less than that of bilayers with symmetric chains. The effects on the bilayer density and lateral stress profiles of changes to the amphiphile architecture are also presented.

Journal Article↗

Lateral and transverse diffusion in two-component bilayer membranes.

We study the lateral and transverse diffusion of amphiphiles in two-component bilayer membranes, using a coarse-grained model for amphiphilic molecules and combined Monte Carlo-Molecular Dynamics simulations. Membrane structural properties, such as the mean thickness, are also measured. The dependence of such properties on membrane composition, inter-molecular interactions, and amphiphile stiffness is determined. In particular, we show that addition of shorter amphiphiles drives the model membrane towards a more fluid state, with increased amphiphile lateral diffusion rates. These results can be understood in the framework of a simple free-volume model. Furthermore, we observe an increase in the trans-membrane diffusion when the interaction energy of amphiphiles with their neighboring molecules is decreased.

Biological Transport↗

Thermally induced proliferation of pores in a model fluid membrane.

The growth of thermally induced pores in a two-dimensional model fluid membrane is investigated by Monte Carlo simulation. Holes appear in the membrane via an activated process, and their subsequent growth is controlled by an edge energy per unit length or line tension. The barrier height and line tension, together with a lateral tension, are the independent parameters of the model. In the resulting phase diagram, a rupture transition separates an intact membrane from a disintegrated state. The approach to the ruptured state shows distinct regimes. Reducing the barrier height at large line tension produces multiple, quasi-independent, small holes whose behavior is dominated by their edge energy, whereas at lower line tensions shape fluctuations of the holes facilitate their coalescence into a single large hole. At a small value of line tension and large barrier height, a single hole spontaneously permeabilizes the membrane in an entropically driven phase transition. Entropy dominates pore growth for line tensions not far below those measured for artificial vesicles. Permeabilization of lipid bilayers by certain peptides involves perturbing lipid-lipid cohesive energies, and our simulations show that at small line tensions the entropy of hole shape fluctuations destroys the model membrane's stability.

Algorithms↗

Entropy-driven instability and rupture of fluid membranes.

A computer simulation is used to investigate hole formation in a model membrane. The model parameters are the stress applied to the membrane, and the edge energy per unit length along the hole boundary (edge tension). Even at zero stress, the membrane has an entropically driven instability against hole formation. Within the model, the minimum edge tension required for the stability of a typical biological membrane is in the region of 1 x 10(-11) J/m, which is similar to the edge tension obtained in many measurements of biomembranes. At the zero-stress instability threshold, the hole shape is the same as a self-avoiding ring, but under compression, the hole shape assumes a branched polymer form. In the presence of large holes at zero stress, the membrane itself behaves like a branched polymer. The boundaries of the phase diagram for membrane stability are obtained, and general features of the rate of membrane rupture under stress are investigated. A model in which the entropy of hole formation is proportional to the hole perimeter is used to interpret the simulation results at small stress near the instability threshold.

Biomechanical Phenomena↗