1 / 44

Frank L. H. Brown Department of Chemistry and Biochemistry, UCSB

Simple models for biomembrane dynamics and structure. Frank L. H. Brown Department of Chemistry and Biochemistry, UCSB. Limitations of Fully Atomic Molecular Dynamics Simulation. A recent “large” membrane simulation (Pitman et. al., JACS, 127 , 4576 (2005))

malory
Download Presentation

Frank L. H. Brown Department of Chemistry and Biochemistry, UCSB

An Image/Link below is provided (as is) to download presentation Download Policy: Content on the Website is provided to you AS IS for your information and personal use and may not be sold / licensed / shared on other websites without getting consent from its author. Content is provided to you AS IS for your information and personal use only. Download presentation by click this link. While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server. During download, if you can't get a presentation, the file might be deleted by the publisher.

E N D

Presentation Transcript

  1. Simple models for biomembrane dynamics and structure Frank L. H. Brown Department of Chemistry and Biochemistry, UCSB

  2. Limitations of Fully Atomic Molecular Dynamics Simulation • A recent “large” membrane simulation • (Pitman et. al., JACS, 127, 4576 (2005)) • 1 rhodopsin, 99 lipids, 24 cholesterols, 7400 waters • (43,222 atoms total) • 5.5 x 7.7 x 10.3 nm periodic box for 118 ns duration • Length/time scales relevant • to cellular biology • ms, m (and longer) • A 1.0 x 1.0 x 0.1 m simulation for 1 ms • would be approximately 2 x 109 more • expensive than our abilities in 2005 • Moore’s law: this might be possible in 46 yrs.

  3. Outline • Elastic simulations for membrane dynamics • Protein motion on the surface of the red blood cell • Supported bilayers and “active” membranes • Molecular membrane model • Flexible lipid, solvent-free • Correspondence with experimental physical properties and stress profile • Protein induced deformation of the bilayer • Elastic model for peristaltic bilayer deformations • Fluctuation spectra • Protein deformations

  4. Helfrich bending free energy: Linear response for normal modes: Ornstein-Uhlenbeck process for each mode: Linear response, curvature elasticity model L h(r) Kc: Bending modulus L: Linear dimension T: Temperature : Cytoplasm viscosity

  5. Relaxation frequencies Non-inertial Navier-Stokes eq: Nonlocal Langevin equation: Oseen tensor (for infinite medium) Bending force Solve for relaxation of membrane modes coupled to a fluid in the overdamped limit: R. Granek, J. Phys. II France, 7, 1761-1788 (1997).

  6.   Membrane near an Impermeable wall Different viscoscities on both sides of membrane Membrane near a semi-permeable wall Generalization to other fluid environments is possible And various combinations Seifert PRE 94, Safran and Gov PRE 04, Lin and Brown (In press)

  7. Membrane Dynamics

  8. Helfrich bending free energy + additional interactions: Overdamped dynamics: (or generalized expressions) • Solve via Brownian dynamics • Handle bulk of calculation in Fourier Space (FSBD) • Efficient handling of hydrodynamics • Natural way to coarse grain over short length scales L. Lin and F. Brown, Phys. Rev. Lett., 93, 256001 (2004). L. Lin and F. Brown, Phys. Rev. E, 72, 011910 (2005). Extension to non-harmonic systems

  9. Protein motion on the surface of red blood cells

  10. S. Liu et al., J. Cell. Biol., 104, 527 (1987).

  11. S. Liu et al., J. Cell. Biol., 104, 527 (1987). • Spectrin “corrals” protein diffusion • Dmicro= 5x10-9 cm2/s (motion inside corral) • Dmacro= 7x10-11 cm2/s (hops between corrals)

  12. Proposed Models

  13. Why membrane undulations? • The mechanism has not been considered before. • “Thermal” membrane fluctuations are known to exist in red blood cells (flicker effect). • Physical constants necessary for a simple model are available from unrelated experiments.

  14. Dynamic undulation model Dmicro Kc=2x10-13 ergs =0.06 poise L=140 nm T=37oC Dmicro=0.53 m2/s h0=6 nm F. Brown, Biophys. J., 84, 842 (2003).

  15. Explicit Cytoskeletal Interactions • Harmonic anchoring of spectrin cytoskeleton to the bilayer • Additional repulsive interaction along the edges of the corral to mimic spectrin L. Lin and F. Brown, Biophys. J., 86, 764 (2004). L. Lin and F. Brown, Phys. Rev. Lett., 93, 256001 (2004).

  16. Dynamics with repulsive spectrin

  17. Escape rate for protein from a corral • Macroscopic diffusion constant on cell surface • (experimentally measured) Information extracted from the simulation • Probability that thermal bilayer fluctuation exceeds h0=6nm at equilibrium (intracellular domain size) • Probability that such a fluctuation persists longer than t0=23s (time to diffuse over spectrin)

  18. Median experimental value Calculated Dmacro • Used experimental median value of corral size L=110 nm

  19. Fluctuations of supported bilayers Y. Kaizuka and J. T. Groves, Biophys. J., 86, 905 (2004). Lin and Brown (submitted).

  20. Way off! Timescales consistent with experiment Fluctuations of supported bilayers (dynamics) Impermeable wall (different boundary conditions) No wall x10-5

  21. Fluctuations of active membranes (experiments) koff J.-B. Manneville, P. Bassereau, D. Levy and J. Prost, PRL, 4356, 1999. kon Off Push down koff kon Off Push up Fluctuations of “active” bilayers L. Lin, N. Gov and F.L.H. Brown, JCP, 124, 074903 (2006).

  22. Summary (elastic modeling) • Elastic models for membrane undulations can be extended to complex geometries and potentials via Brownian dynamics simulation. • “Thermal” undulations appear to be able to promote protein mobility on the RBC. • Other biophysical and biochemical systems are well suited to this approach.

  23. Molecular membrane models with implicit solvent

  24. Why? • Study the basic (minimal) requirements for bilayer stability and elasticity. • A particle based method is more versatile than elastic models. • Incorporate proteins, cytoskeleton, etc. • Non “planar” geometries • Computational efficiency.

  25. A range of resolutions... W. Helfrich, Z. Natuforsch, 28c, 693 (1973) J.M. Drouffe, A.C. Maggs, and S. Leibler, Science 254, 1353 (1991) Helmut Heller, Michael Schaefer, and Klaus Schulten, J. Phys. Chem 1993, 8343 (1993) Y. Kantor, M. Kardar, and D.R. Nelson, Phys. Rev. A 35, 3056 (1987) R. Goetz and R. Lipowsky, J. Chem. Phys. 108, 7397 (1998)

  26. Motivation D. Marsh, Biochim. Biophys. Acta, 1286, 183-223 (1996).

  27. A solvent free model G. Brannigan, P. F. Philips and F. L. H. Brown, Phys. Rev. E, 72, 011915 (2005).

  28. Self-Assembly (128 lipids)

  29. Phase Behavior Fluid (kBT = 0.9) “Gel” (kBT = 0.5)

  30. Flexible (and tunable) kc: 1 - 8 * 10-20 J kA: 40-220 mJ/m2

  31. Stress Profile(-surface tension vs. height) R. Goetz and R.Lipowsky, JCP 108 p. 7397, (1998) Explicit Solvent Implicit Solvent

  32. Protein Inclusions Single protein inclusion in the bilayer sheet Inclusion is composed of rigid “lipid molecules”

  33. Deformation of the bilayer data fit to elastic r H. Aranda-Espinoza et. al., Biophys. J., 71, 648 (1996). G. Brannigan and FLH Brown, Biophys. J., 90, 1501 (2006).

  34. Direct from simulation Thermal fluctuations Infer from stress Profile & fluctuations The fit is meaningful Using physical constant inferred from homogeneous bilayer simulations we can predict the “fit” values. Deformation profile has three independent constants formed from combinations of elastic properties

  35. A consistent elastic model protrusion bending • Begin with standard treatment for surfactant monolayers (e.g. Safran) • Couple monolayers by requiring volume conservation of hydrophobic tails and equal local area/lipid between leaflets • Include microscopic protrusions (harmonically bound to z fields & include interfacial tension)

  36. Undulations (bending & protrusion) Peristaltic bending Peristaltic protrusions (and coupling to bending) A consistent elastic model

  37. A consistent elastic model Predicts thermal fluctuation amplitudes for both undulation and peristaltic modes Continuous as function of q Spontaneous curvature included Elastic constants are interrelated 5 fit constants over two data sets Predicts response to hydrophobic mismatch Same elastic constants as fluctuations Protrusion modes explicitly considered

  38. Fits to fluctuation spectra Lindahl &Edholm Marrink & Mark Scott & coworkers

  39. These molecules are happy because they’re in their preferred curvature. These molecules are unhappy, but, there are relatively few of them. Effect of spontaneous curvature Positive spontaneous curvature & conservation of volume favors modulated thickness

  40. Fluctuation spectra: spontaneous curvature Of the available simulation data, this effect is most pronounced for DPPC: Lindahl, E. and O. Edholm. Biophysical Journal79:426(2000)

  41. gramicidin A channel lifetimes Rate vs. monolayer thickness

  42. In progress… Free energy calculations to establish potential of mean force.

  43. Summary • Implicit solvent models can reproduce an array of bilayer properties and behaviors. • Computation is significantly reduced relative to explicit solvent models, enabling otherwise difficult simulations. • A single elastic model can explain peristaltic fluctuations, height fluctuations and response to a protein deformations (if you include spontaneous curvature and protrusions)

  44. Acknowledgements Lawrence Lin Grace Brannigan Adele Tamboli Peter Phillips Jay Groves, Larry Scott, Jay Mashl Siewert-Jan Marrink Olle Edholm NSF, ACS-PRF, Sloan Foundation,UCSB

More Related