Motivation
Many biological processes of interest occur on length and time scales inaccessible to atomistic and molecular dynamics simulations. For these particular applications, we must use simplified models of biological membranes to have any hope of quantitative analysis. We have developed an algorithm for the time evolution of a membrane in arbitrary potentials. In this scheme, harmonic interactions can be handled analytically while simulations can be used to deal with general forces.

An example of a process that can be studied with these methods is the diffusive motion of proteins on the surface of the red blood cells.  Single particle tracking experiments have shown that proteins are trapped in corrals (formed by the cytoskeleton beneath the membrane) with sizes of about a hundred nanometers for periods of about a hundred microseconds.  Here, our method is suitable for analyzing the role of thermal undulations of the membrane as a mechanism for carrying the protein over the cytoskeleton (see figure).


Model
Continuous sheets with elastic free energy are a natural model for studying the behavior of lipid bilayers on longer length and time scales. In the past, continuum models with curvature elasticity have been used to explain the various shapes of the red blood cell under different external conditions, as well as the repulsive behavior of interacting membranes. When thermal fluctuations are added to these elastic models, they provide an explanation for the observed flickering of erythrocytes when viewed under a light microscope.

Our membrane model consists of a curvature free energy plus additional arbitrary potentials, hydrodynamic coupling, and thermal fluctuations through a Langevin random force. All of the equations are most easily handled in Fourier space due to either decoupling of modes for analytical work or for reasons of computational efficiency for simulations. Our Fourier Space Brownian Dynamics (FSBD) algorithm is an efficient method for evolving membrane configurations in simulations that can be used for general systems where membrane undulations are important.

The elastic membrane model is, to some extent, analytically tractable for harmonic potentials. Specifically, we study localized harmonic interactions that pin the membrane to the underlying cytoskeleton and calculate the probability that a thermal undulation is large enough and persists long enough for the passage of the protein out of the corral (see figures). This escape probability is completely characterized by the height autocorrelation function which can be computed without simulation. A diffusion constant can then be estimated for the protein diffusing in the plane of the membrane. The results suggest that thermal fluctuations are a possible mechanism for promoting the passage of proteins over the cytoskeletal barrier.

To extend our model to include interactions between the cytoskeleton and the membrane, we can add a repulsive potential at the edges of the corrals (see figures). This interaction is nonharmonic and therefore requires simulation over small time steps. The results for the diffusion constant are modified somewhat from the previous case, but still suggest that thermal fluctuations may be a mechanism that contributes to the mobility of proteins on the surface of the red blood cell. More detail is presented in the following papers.


Related Publications
Dynamic Simulations of Membranes with Cytoskeletal Interactions, Lawrence C.-L. Lin  and F. L. H. Brown, Physical Review E, 72, 011910 (2005).
Brownian Dynamics in Fourier Space: Membrane Simulations Over Long Length and Time Scales, Lawrence C.-L. Lin  and F. L. H. Brown, Physical Review Letters, 93, 256001 (2004).
Dynamics of Pinned Membranes with Application to Protein Diffusion on the Surface of Red Blood Cells, Lawrence C.-L. Lin  and F. L. H. Brown, Biophysical Journal, 86, 764-780 (2004).
Regulation of Protein Mobility Via Thermal Membrane Undulations, F. L. H. Brown, Biophysical Journal, 84, 842-853 (2003).
Lateral diffusion of membrane proteins in the presence of static and dynamic corrals - suggestions for appropriate observables. F.L. H. Brown, D. M. Leitner, J. A. McCammon and K. R. Wilson, Biophysical Journal, 78, 2257-2269 (2000).
Regulation of protein mobility in cell membranes: A dynamic corral model, D. M. Leitner, F. L. H. Brown, and K. R. Wilson, Biophysical Journal, 78, 125-135 (2000).