PubMed Health⌕ Search

Biomedical subjects

Patrick McMurtry

Publications and source records attributed to Patrick McMurtry.

3 recordsLinked to original sources

Coupling field theory with continuum mechanics: a simulation of domain formation in giant unilamellar vesicles.

Domain formation is modeled on the surface of giant unilamellar vesicles using a Landau field theory model for phase coexistence coupled to elastic deformation mechanics (e.g., membrane curvature). Smooth particle applied mechanics, a form of smoothed particle continuum mechanics, is used to solve either the time-dependent Landau-Ginzburg or Cahn-Hilliard free-energy models for the composition dynamics. At the same time, the underlying elastic membrane is modeled using smooth particle applied mechanics, resulting in a unified computational scheme capable of treating the response of the composition fields to arbitrary deformations of the vesicle and vice versa. The results indicate that curvature coupling, along with the field theory model for composition free energy, gives domain formations that are correlated with surface defects on the vesicle. In the case that external deformations are included, the domain structures are seen to respond to such deformations. The present simulation capability provides a significant step forward toward the simulation of realistic cellular membrane processes.

Biomechanical Phenomena↗

Interfacing molecular dynamics and macro-scale simulations for lipid bilayer vesicles.

A continuum-level model for a giant unilamellar vesicle (GUV) is bridged to a corresponding atomistic model of a dimyristoylphosphatidylcholine (DMPC) bilayer at various cholesterol concentrations via computation of the bulk modulus. The bulk modulus and other microscopically determined parameters are passed to a continuum-level model operating in time- and length-scales orders of magnitude beyond that which is accessible by atomistic-level simulation. The continuum-level simulation method used is the material point method (MPM), and the particular variation used here takes advantage of the spherical nature of many GUVs. An osmotic pressure gradient due to a solvent concentration change is incorporated into the continuum-level simulation, resulting in osmotic swelling of the vesicle. The model is then extended to treat mixtures of DMPC and cholesterol, where small domains of different composition are considered.

Algorithms↗

Calculating the bulk modulus for a lipid bilayer with nonequilibrium molecular dynamics simulation.

Nonequilibrium molecular dynamics (NEMD) computer simulations are used to calculated the bulk modulus for a dimyristoylphosphatidylcholine bilayer. A methodology is developed whereby NEMD can be effectively used to calculate material properties for complex systems that undergo long time-scale conformational changes. It is found that the bulk modulus upon expansion from a zero stress state agrees well with experimental estimates. However, it is also found that the modulus upon contraction from a zero stress state is larger. From a molecular perspective, it is possible to explain this phenomena by examining the molecular origins of the pressure response. The finding that the two moduli are not equal upon compression and expansion is in apparent contradiction to osmotic stress experiments where the area modulus was found to be the same upon expansion and contraction. This issue is addressed.

Algorithms↗