Propagating Subterahertz Longitudinal Phonon-like Modes in a Lipid Bilayer in Aqueous Medium

It has been revealed recently that the subterahertz/terahertz vibrational motions in enzymes and DNA immersed in aqueous solutions can be underdamped. Importantly, these motions are associated with coherent delocalized modes that control functional processes. Analogous propagating phonon-like modes have been found in free hydrated lipid bilayers. In the present work, subterahertz (frequencies of the order of tens and hundreds of gigahertz) longitudinal acoustic oscillations in a bilayer lipid membrane immersed in aqueous medium are investigated theoretically. We consider driven oscillations excited by tangential mechanical tensions at the bilayer surfaces and thermally induced phonon modes. The analysis is based on: (i) a generalized hydrodynamic model of two-dimensional lipid bilayer in aqueous medium; (ii) known estimates of frequencies and lifetimes of longitudinal acoustic phonons in free hydrated lipid bilayer and in water, which were obtained in the experiments on non-elastic X-ray scattering and the molecular dynamics simulations. We show that the membrane phonon-like excitations are underdamped for the typical values of the system parameters, and the contribution of aqueous medium to the membrane mode damping is small compared to the contribution of the lipid bilayer. The obtained results suggest the possibility of realization of thermally induced longitudinal membrane phonons in physiological conditions, as well as the possibility of resonance amplification of the impact of subnanosecond electric impulses and impulses of subterahertz electromagnetic radiation on membrane dynamics.     

Cholesterol-induced protein sorting: an analysis of energetic feasibility

The mechanism(s) underlying the sorting of integral membrane proteins between the Golgi complex and the plasma membrane remain uncertain because no specific Golgi retention signal has been found. Moreover one can alter a protein's eventual localization simply by altering the length of its transmembrane domain (TMD). M. S. Bretscher and S. Munro (SCIENCE: 261:1280-1281, 1993) therefore proposed a physical sorting mechanism based on the hydrophobic match between the proteins' TMD and the bilayer thickness, in which cholesterol would regulate protein sorting by increasing the lipid bilayer thickness. In this model, Golgi proteins with short TMDs would be excluded from cholesterol-enriched domains (lipid rafts) that are incorporated into transport vesicles destined for the plasma membrane. Although attractive, this model remains unproven. We therefore evaluated the energetic feasibility of a cholesterol-dependent sorting process using the theory of elastic liquid crystal deformations. We show that the distribution of proteins between cholesterol-enriched and cholesterol-poor bilayer domains can be regulated by cholesterol-induced changes in the bilayer physical properties. Changes in bilayer thickness per se, however, have only a modest effect on sorting; the major effect arises because cholesterol changes also the bilayer material properties, which augments the energetic penalty for incorporating short TMDs into cholesterol-enriched domains. We conclude that cholesterol-induced changes in the bilayer physical properties allow for effective and accurate sorting which will be important generally for protein partitioning between different membrane domains.     

Calculation of the solution properties of flexible macromolecules: methods and applications

While the prediction of hydrodynamic properties of rigid particles is nowadays feasible using simple and efficient computer programs, the calculation of such properties and, in general, the dynamic behavior of flexible macromolecules has not reached a similar situation. Although the theories are available, usually the computational work is done using solutions specific for each problem. We intend to develop computer programs that would greatly facilitate the task of predicting solution behavior of flexible macromolecules. In this paper, we first present an overview of the two approaches that are most practical: the Monte Carlo rigid-body treatment, and the Brownian dynamics simulation technique. The Monte Carlo procedure is based on the calculation of properties for instantaneous conformations of the macromolecule that are regarded as if they were instantaneously rigid. We describe how a Monte Carlo program can be interfaced to the programs in the HYDRO suite for rigid particles, and provide an example of such calculation, for a hypothetical particle: a protein with two domains connected by a flexible linker. We also describe briefly the essentials of Brownian dynamics, and propose a general mechanical model that includes several kinds of intramolecular interactions, such as bending, internal rotation, excluded volume effects, etc. We provide an example of the application of this methodology to the dynamics of a semiflexible, wormlike DNA.     

Intermonolayer friction and surface shear viscosity of lipid bilayer membranes

The flow behavior of lipid bilayer membranes is characterized by a surface viscosity for in-plane shear deformations, and an intermonolayer friction coefficient for slip between the two leaflets of the bilayer. Both properties have been studied for a variety of coarse-grained double-tailed model lipids, using equilibrium and nonequilibrium molecular dynamics simulations. For lipids with two identical tails, the surface shear viscosity rises rapidly with tail length, while the intermonolayer friction coefficient is less sensitive to the tail length. Interdigitation of lipid tails across the bilayer midsurface, as observed for lipids with two distinct tails, strongly enhances the intermonolayer friction coefficient, but hardly affects the surface shear viscosity. The simulation results are compared against the available experimental data.     

Hydrodynamics of segmentally flexible macromolecules

Segmentally flexible macromolecules are composed of a few rigid subunits linked by joints which are more or less flexible. The dynamics in solution of this type of macromolecule present special aspects that are reviewed here. Three alternative approaches are described. One is the rigid-body treatment, which is shown to be valid for overall dynamic properties such as translational diffusion and intrinsic viscosity. Another approach is the Harvey-Wegener treatment, which is particularly suited for rotational diffusion. The simplest version of this treatment, which ignores hydrodynamic interaction (HI) effects, is found to be quite accurate when compared to a more rigorous version including HI. A third approach is the Brownian dynamics simulation that, albeit at some computational cost, might describe rigorously cases of arbitrary complexity. This technique has been used to test the approximations in the rigid-body and Harvey-Wegener treatments, thus allowing a better understanding of their validity. Brownian trajectories of simplified models such as the trumbbell and the broken rod have been simulated. The comparison of the decay rates of some correlation functions with the predictions of the two treatments leads to a general conclusion: the Harvey-Wegener treatment determines the initial rate, while the long-time behavior is dominated by the rigid-body relaxation time. As an example of application to a specific biological macromolecule, we present a simulation of an immunoglobulin molecule, showing how Brownian Dynamics can be used to predict rotational and internal dynamics. Another typical example is myosin. Literature data of hydrodynamic properties of whole myosin and the myosin rod are compared with predictions from the Harvey-Wegener and rigid-body treatments. The present situation of the problem on myosin flexibility is analyzed, and some indications are given for future experimental and simulation work.     

The Fluid–Mosaic model of cell membranes: A brief introduction,historical features,some general principles,and its adaptation to current information

The Fluid–Mosaic Membrane (FMM) model was originally proposed as a general, nanometer-scale representation of cell membranes (Singer and Nicolson, 1972). The FMM model was based on some general principles, such as thermodynamic considerations, intercalation of globular proteins into a lipid bilayer, independent protein and lipid dynamics, cooperativity and other characteristics. Other models had trimolecular structures or membrane globular lipoprotein units. These latter models were flawed, because they did not allow autonomous lipids, membrane domains or discrete lateral dynamics. The FMM model was also consistent with membrane asymmetry, cis- and trans-membrane linkages and associations of membrane components into multi-molecular complexes and domains. It has remained useful for explaining the basic organizational principles and properties of various biological membranes. New information has been added, such as membrane-associated cytoskeletal assemblies, extracellular matrix interactions, transmembrane controls, specialized lipid-protein domains that differ in compositions, rotational and lateral mobilities, lifetimes, functions, and other characteristics. The presence of dense, structured membrane domains has reduced significantly the extent of fluid-lipid membrane areas, and the FMM model is now considered to be more mosaic and dense than the original proposal.     

Molecular dynamics study of the lauryl alcohol-laurate model bilayer

Molecular dynamics (MD) calculation of the fluid phase lauryl alcohol-laurate bilayer has been executed based on Berendsen's surface-constrained model. Structure and dynamics of the bilayer have been investigated by analyzing the trajectories of the chain configurations. Newly defined correlation functions as well as the conventional ones showed that the tilt and bend of the chain play an important role in the bilayer structure, including behavior of the order parameter. Interpenetration of the layers as well as formation of collectively ordered small domains was also found. The calculated lateral diffusion coefficient was in satisfactory agreement with the experimental one. Successive jumps of the head group, rather than the hydrodynamic continuous motion, were observed. Between the jumps, the molecule librated in a local site. Time-dependent autocorrelation functions showed evidence of several different modes of the chain motion, whose time constant ranged from a few tenths of picoseconds to several tens of picoseconds.     

Setting up and running molecular dynamics simulations of membrane proteins

Molecular dynamics simulations have become a popular and powerful technique to study lipids and membrane proteins. We present some general questions and issues that should be considered prior to embarking on molecular dynamics simulation studies of membrane proteins and review common simulation methods. We suggest a practical approach to setting up and running simulations of membrane proteins, and introduce two new (related) methods to embed a protein in a lipid bilayer. Both methods rely on placing lipids and the protein(s) on a widely spaced grid and then 'shrinking' the grid until the bilayer with the protein has the desired density, with lipids neatly packed around the protein. When starting from a grid based on a single lipid structure, or several potentially different lipid structures (method 1), the bilayer will start well-packed but requires more equilibration. When starting from a pre-equilibrated bilayer, either pure or mixed, most of the structure of the bilayer stays intact, reducing equilibration time (method 2). The main advantages of these methods are that they minimize equilibration time and can be almost completely automated, nearly eliminating one time consuming step in MD simulations of membrane proteins.     

Molecular and hydrodynamic properties of human epidermal growth factor receptor HER2 extracellular domain and its homodimer: Experiments and multi-scale simulations

BackgroundIn a broad range of human carcinomas gene amplification leads to HER2 overexpression, which has been proposed to cause spontaneous dimerization and activation in the absence of ligand. This makes HER2 attractive as a therapeutic target. However, the HER2 homodimerization mechanism remains unexplored. It has been suggested that the “back-to-back” homodimer does not form in solution. Notwithstanding, very recently the crystal structure of the HER2 extracellular domain homodimer formed with a “back-to-head” interaction has been resolved. We intend to explore the existence of such interactions.MethodsA combination of experiments, molecular dynamics and hydrodynamic modeling were used to monitor the transport properties of HER2 in solution.Results & conclusionsWe have detected the HER2 extracellular domain homodimer in solution. The results show a high degree of molecular flexibility, which ultimately leads to quite higher values of the intrinsic viscosity and lower values of diffusion coefficient than those corresponding to globular proteins. This flexibility obeys to the open conformation of the receptor and to the large fluctuations of the different domains. We also report that for obtaining the correct hydrodynamic constants from the modeling one must consider the glycosylation of the systems.General significanceConformational features of epidermal growth factor receptors regulate their hydrodynamic properties and control their activity. It is essential to understand the dynamics of these systems and the role of the specific domains involved. To find biophysical correlations between dynamics and macroscopic transport properties is of general interest for researches working in this area. This article is part of a Special Issue entitled “Biochemistry of Synthetic Biology - Recent Developments” Guest Editor: Dr. Ilka Heinemann and Dr. Patrick O'Donoghue.     

Membrane structural perturbations caused by anesthetics and nonimmobilizers: a molecular dynamics investigation.

The structural perturbations of the fully hydrated dimyristoyl-phosphatidylcholine bilayer induced by the presence of hexafluoroethane C(2F6), a "nonimmobilizer," have been examined by molecular dynamics simulations and compared with the effects produced by halothane CF3CHBrCl, an "anesthetic," on a similar bilayer (DPPC) (Koubi et al., Biophys. J. 2000. 78:800). We find that the overall structure of the lipid bilayer and the zwitterionic head-group dipole orientation undergo only a slight modification compared with the pure lipid bilayer, with virtually no change in the potential across the interface. This is in contrast to the anesthetic case in which the presence of the molecule led to a large perturbation of the electrostatic potential across to the membrane interface. Similarly, the analysis of the structural and dynamical properties of the lipid core are unchanged in the presence of the nonimmobilizer although there is a substantial increase in the microscopic viscosity for the system containing the anesthetic. These contrasting perturbations of the lipid membrane caused by those quite similarly sized molecules may explain the difference in their physiological effects as anesthetics and nonimmobilizers, respectively.     

A hybrid model for erythrocyte membrane: a single unit of protein network coupled with lipid bilayer

To investigate the nanomechanics of the erythrocyte membrane we developed a hybrid model that couples the actin-spectrin network to the lipid bilayer. This model features a Fourier space Brownian dynamics model of the bilayer, a Brownian dynamics model of the actin protofilament, and a modified wormlike-chain model of the spectrin (including a cable-dynamics model to predict the oscillation in tension). This model enables us to predict the nanomechanics of single or multiple units of the protein network, the lipid bilayer, and the effect of their interactions. The present work is focused on the attitude of the actin protofilament at the equilibrium states coupled with the elevations of the lipid bilayer through their primary linkage at the suspension complex in deformations. Two different actin-spectrin junctions are considered at the junctional complex. With a point-attachment junction, large pitch angles and bifurcation of yaw angles are predicted. Thermal fluctuations at bifurcation may lead to mode-switching, which may affect the network and the physiological performance of the membrane. In contrast, with a wrap-around junction, pitch angles remain small, and the occurrence of bifurcation is greatly reduced. These simulations suggest the importance of three-dimensional molecular junctions and the lipid bilayer/protein network coupling on cell membrane mechanics.     

Finding a Needle in a Haystack: The Role of Electrostatics in Target Lipid Recognition by PH Domains

Interactions between protein domains and lipid molecules play key roles in controlling cell membrane signalling and trafficking. The pleckstrin homology (PH) domain is one of the most widespread, binding specifically to phosphatidylinositol phosphates (PIPs) in cell membranes. PH domains must locate specific PIPs in the presence of a background of approximately 20% anionic lipids within the cytoplasmic leaflet of the plasma membrane. We investigate the mechanism of such recognition via a multiscale procedure combining Brownian dynamics (BD) and molecular dynamics (MD) simulations of the GRP1 PH domain interacting with phosphatidylinositol (3,4,5)-trisphosphate (PI(3,4,5)P(3)). The interaction of GRP1-PH with PI(3,4,5)P(3) in a zwitterionic bilayer is compared with the interaction in bilayers containing different levels of anionic 'decoy' lipids. BD simulations reveal both translational and orientational electrostatic steering of the PH domain towards the PI(3,4,5)P(3)-containing anionic bilayer surface. There is a payoff between non-PIP anionic lipids attracting the PH domain to the bilayer surface in a favourable orientation and their role as 'decoys', disrupting the interaction of GRP1-PH with the PI(3,4,5)P(3) molecule. Significantly, approximately 20% anionic lipid in the cytoplasmic leaflet of the bilayer is nearly optimal to both enhance orientational steering and to localise GRP1-PH proximal to the surface of the membrane without sacrificing its ability to locate PI(3,4,5)P(3) within the bilayer plane. Subsequent MD simulations reveal binding to PI(3,4,5)P(3), forming protein-phosphate contacts comparable to those in X-ray structures. These studies demonstrate a computational framework which addresses lipid recognition within a cell membrane environment, offering a link between structural and cell biological characterisation.     

Immobilization and activity assay of cytochrome P450 on patterned lipid membranes

We report on a methodology for immobilizing cytochrome P450 on the surface of micropatterned lipid bilayer membranes and measuring the enzymatic activity. The patterned bilayer comprised a matrix of polymeric lipid bilayers and embedded fluid lipid bilayers. The polymeric lipid bilayer domains act as a barrier to confine fluid lipid bilayers in defined areas and as a framework to stabilize embedded membranes. The fluid bilayer domains, on the other hand, can contain lipid compositions that facilitate the fusion between lipid membranes, and are intended to be used as the binding agent of microsomes containing rat CYP1A1. By optimizing the membrane compositions of the fluid bilayers, we could selectively immobilize microsomal membranes on these domains. The enzymatic activity was significantly higher on lipid bilayer substrates compared with direct adsorption on glass. Furthermore, competitive assay experiment between two fluorogenic substrates demonstrated the feasibility of bioassays based on immobilized P450s.     

Flicker spectroscopy of thermal lipid bilayer domain boundary fluctuations

This contribution describes measurements of lipid bilayer domain line tension based on two-dimensional thermal undulations of membranes with liquid ordered/liquid disordered phase coexistence and near-critical composition at room temperature. Lateral inhomogeneity of lipid and protein composition is currently a subject of avid research aimed at determining both fundamental properties and biological relevance of membrane domains. Line tension at fluid lipid bilayer membrane domain boundaries controls the kinetics of domain growth and therefore regulates the size of compositional heterogeneities. High line tension promotes membrane domain budding and fission. Line tension could therefore be an important control parameter regulating functional aspects of biological membranes. Here the established method of fluid domain flicker spectroscopy is applied to examine thermal domain wall fluctuations of phase-separated bilayer membranes. We find a Gaussian probability distribution for the first few excited mode amplitudes, which permits an analysis by means of appropriately specialized capillary wave theory. Time autocorrelation functions are found to decay exponentially, and relaxation times are fitted by means of a hydrodynamic theory relating line tensions and excited mode relaxation kinetics. Line tensions below 1 pN are obtained, with these two approaches yielding similar results. We examine experimental artifacts that perturb the Fourier spectrum of domain traces and discuss ways to identify the number of modes that yield reliable line tension information.     

Electrostriction and dynamics of solid supported lipid films

This review reports the significance of bilayer lipid membranes on a solid support (sBLM) for the construction of biosensors. The methods of formation of lipid membranes on different solid supports including different metals (silver, gold, stainless steel), agar and conducting polymers are presented. Several examples of the application of electrostriction and dielectric relaxation methods for the study of mechanical properties and dynamics of solid supported bilayers have been shown. We demonstrated that these methods are useful for determination of the binding of enzymes and antibodies to sBLM, for the study of hybridization of nucleic acids on membrane surfaces and for the study of physical properties of modified supported membranes.     

Protein shape change has a major effect on the gating energy of a mechanosensitive channel

Increasing experimental evidence has shown that membrane protein functionality depends on molecular composition of cell membranes. However, the origin of this dependence is not fully understood. It is reasonable to assume that specific lipid-protein interactions are important, yet more generic effects due to mechanical properties of lipid bilayers likely play a significant role too. Previously it has been demonstrated using models for elastic properties of membranes and lateral pressure profiles of lipid bilayers that the mechanical properties of a lipid bilayer can contribute as much as ∼10 k_BT to the free energy difference associated with a change in protein conformational state. Here, we extend those previous approaches to a more realistic model for a large mechanosensitive channel (MscL). We use molecular dynamics together with the MARTINI model to simulate the open and closed states of MscL embedded in a DOPC bilayer. We introduce a procedure to calculate the mechanical energy change in the channel gating using a three-dimensional pressure distribution inside a membrane, computed from the molecular dynamics simulations. We decompose the mechanical energy to terms associated with area dilation and shape contribution. Our results highlight that the lateral pressure profile of a lipid bilayer together with the shape change in gating can induce a contribution of ∼30 k_BT on the gating energy of MscL. This contribution arises largely from the interfacial tension between hydrophobic and hydrophilic regions in a lipid bilayer.     

Crossover from picosecond collective to single particle dynamics defines the mechanism of lateral lipid diffusion

It has been widely accepted that the thermally excited motions of the molecules in a cell membrane is the prerequisite for a cell to carry its biological functions. On the other hand, the detailed mapping of the ultrafast picosecond single-molecule and the collective lipid dynamics in a cell membrane remains rather elusive. Here, we report all-atom molecular dynamics simulations of a 1,2-dipalmitoyl-sn-glycero-3-phosphocholine bilayer over a wide range of temperature. We elucidate a molecular mechanism underlying the lateral lipid diffusion in a cell membrane across the gel, rippled, and liquid phases using an analysis of the longitudinal and transverse current correlation spectra, the velocity auto-correlation functions, and the molecules mean square displacements. The molecular mechanism is based on the anomalous ultrafast vibrational properties of lipid molecules at the viscous-to-elastic crossover. The macroscopic lipid diffusion coefficients predicted by the proposed diffusion model are in a good agreement with experimentally observed values. Furthermore, we unveil the role of water confined at the water-lipid interface in triggering collective vibrations in a lipid bilayer.     

Molecular mechanisms of spontaneous curvature and softening in complex lipid bilayer mixtures

Membrane reshaping is an essential biological process. The chemical composition of lipid membranes determines their mechanical properties and thus the energetics of their shape. Hundreds of distinct lipid species make up native bilayers, and this diversity complicates efforts to uncover what compositional factors drive membrane stability in cells. Simplifying assumptions, therefore, are used to generate quantitative predictions of bilayer dynamics based on lipid composition. One assumption commonly used is that “per lipid” mechanical properties are both additive and constant—that they are an intrinsic property of lipids independent of the surrounding composition. Related to this is the assumption that lipid bulkiness, or “shape,” determines its curvature preference, independently of context. In this study, all-atom molecular dynamics simulations on three separate multilipid systems were used to explicitly test these assumptions, applying methodology recently developed to isolate properties of single lipids or nanometer-scale patches of lipids. The curvature preference experienced by populations of lipid conformations were inferred from their redistribution on a dynamically fluctuating bilayer. Representative populations were extracted by both structural similarity and semi-automated hidden Markov model analysis. The curvature preferences of lipid dimers were then determined and compared with an additive model that combines the monomer curvature preference of both the individual lipids. In all three systems, we identified conformational subpopulations of lipid dimers that showed non-additive curvature preference, in each case mediated by a special chemical interaction (e.g., hydrogen bonding). Our study highlights the importance of specific chemical interactions between lipids in multicomponent bilayers and the impact of interactions on bilayer stiffness. We identify two mechanisms of bilayer softening: diffusional softening, driven by the dynamic coupling between lipid distributions and membrane undulations, and conformational softening, driven by the inter-conversion between distinct dimeric conformations.     

Lipid Bilayer Domain Fluctuations as a Probe of Membrane Viscosity

We argue that membrane viscosity, η_m, plays a prominent role in the thermal fluctuation dynamics of micron-scale lipid domains. A theoretical expression is presented for the timescales of domain shape relaxation, which reduces to the well-known η_m = 0 result of Stone and McConnell in the limit of large domain sizes. Experimental measurements of domain dynamics on the surface of ternary phospholipid and cholesterol vesicles confirm the theoretical results and suggest domain flicker spectroscopy as a convenient means to simultaneously measure both the line tension, σ, and the membrane viscosity, η_m, governing the behavior of individual lipid domains.     

Structure and dynamics of sphingomyelin bilayer: insight gained through systematic comparison to phosphatidylcholine

Sphingomyelin, one of the main lipid components of biological membranes, is actively involved in various cellular processes such as protein trafficking and signal transduction. In particular, specific lateral domains enriched in sphingomyelin and cholesterol have been proposed to play an important functional role in biomembranes, although their precise characteristics have remained unclear. A thorough understanding of the functional role of membranes requires detailed knowledge of their individual lipid components. Here, we employ molecular dynamics simulations to conduct a systematic comparison of a palmitoylsphingomyelin (PSM, 16:0-SM) bilayer with a membrane that comprises dipalmitoylphosphatidylcholine (DPPC) above the main phase transition temperature. We clarify atomic-scale properties that are specific to sphingomyelin due to its sphingosine moiety, and further discuss their implications for SM-rich membranes. We find that PSM bilayers, and in particular the dynamics of PSM systems, are distinctly different from those of a DPPC bilayer. When compared with DPPC, the strong hydrogen bonding properties characteristic to PSM are observed to lead to considerable structural changes in the polar headgroup and interface regions. The strong ordering of PSM acyl chains and specific ordering effects in the vicinity of a PSM-water interface reflect this issue clearly. The sphingosine moiety and related hydrogen bonding further play a crucial role in the dynamics of PSM bilayers, as most dynamic properties, such as lateral and rotational diffusion, are strongly suppressed. This is most evident in the rotational motion characterized by spin-lattice relaxation times and the decay of hydrogen bond autocorrelation functions that are expected to be important in complexation of SM with other lipids in many-component bilayers. A thorough understanding of SM bilayers would greatly benefit from nuclear magnetic resonance experiments for acyl chain ordering and dynamics, allowing full comparison of these simulations to experiments.