|
|||
| 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). |