No Arabic abstract
We consider the hydrodynamics of lipid bilayers containing transmembrane proteins of arbitrary shape. This biologically-motivated problem is relevant to the cell membrane, whose fluctuating dynamics play a key role in phenomena ranging from cell migration, intercellular transport, and cell communication. Using Onsagers variational principle, we derive the equations that govern the relaxation dynamics of the membrane shape, of the mass densities of the bilayer leaflets, and of the diffusing proteins concentration. With our generic formalism, we obtain several results on membrane dynamics. We find that proteins that span the bilayer increase the intermonolayer friction coefficient. The renormalization, which can be significant, is in inverse proportion to the proteins mobility. Second, we find that asymmetric proteins couple to the membrane curvature and to the difference in monolayer densities. For practically all accessible membrane tensions ($sigma> 10^{-8}$ N/m) we show that the protein density is the slowest relaxing variable. Furthermore, its relaxation rate decreases at small wavelengths due to the coupling to curvature. We apply our formalism to the large-scale diffusion of a concentrated protein patch. We find that the diffusion profile is not self-similar, owing to the wavevector dependence of the effective diffusion coefficient.
We carry out a coarse-grained molecular dynamics simulation of phospholipid vesicles with transmembrane proteins. We measure the mean and Gaussian curvatures of our protein-embedded vesicles and quantitatively show how protein clusters change the shapes of their host vesicles. The effects of depletion force and vesiculation on protein clustering are also investigated. By increasing the protein concentration, clusters are fragmented to smaller bundles, which are then redistributed to form more symmetric structures corresponding to lower bending energies. Big clusters and highly aspherical vesicles cannot be formed when the fraction of protein to lipid molecules is large.
Nearly a quarter of genomic sequences and almost half of all receptors that are likely to be targets for drug design are integral membrane proteins. Understanding the detailed mechanisms of the folding of membrane proteins is a largely unsolved, key problem in structural biology. Here, we introduce a general model and use computer simulations to study the equilibrium properties and the folding kinetics of a $C_{alpha}$-based two helix bundle fragment (comprised of 66 amino-acids) of Bacteriorhodopsin. Various intermediates are identified and their free energy are calculated toghether with the free energy barrier between them. In 40% of folding trajectories, the folding rate is considerably increased by the presence of non-obligatory intermediates acting as traps. In all cases, a substantial portion of the helices is rapidly formed. This initial stage is followed by a long period of consolidation of the helices accompanied by their correct packing within the membrane. Our results provide the framework for understanding the variety of folding pathways of helical transmembrane proteins.
Protein aggregation in cell membrane is vital for the majority of biological functions. Recent experimental results suggest that transmembrane domains of proteins such as $alpha$-helices and $beta$-sheets have different structural rigidities. We use molecular dynamics simulation of a coarse-grained model of protein-embedded lipid membranes to investigate the mechanisms of protein clustering. For a variety of protein concentrations, our simulations under thermal equilibrium conditions reveal that the structural rigidity of transmembrane domains dramatically affects interactions and changes the shape of the cluster. We have observed stable large aggregates even in the absence of hydrophobic mismatch which has been previously proposed as the mechanism of protein aggregation. According to our results, semi-flexible proteins aggregate to form two-dimensional clusters while rigid proteins, by contrast, form one-dimensional string-like structures. By assuming two probable scenarios for the formation of a two-dimensional triangular structure, we calculate the lipid density around protein clusters and find that the difference in lipid distribution around rigid and semiflexible proteins determines the one- or two-dimensional nature of aggregates. It is found that lipids move faster around semiflexible proteins than rigid ones. The aggregation mechanism suggested in this paper can be tested by current state-of-the-art experimental facilities.
Bilayer membranes self-assembled from amphiphilic molecules such as lipids, surfactants and block copolymers are ubiquitous in biological and physiochemical systems. The shape and structure of bilayer membranes depend crucially on their mechanical properties such surface tension, bending moduli and line tension. Understanding how the molecular property of the amphiphiles determine the structure and mechanics of the self-assembled bilayers requires a molecularly detailed theoretical framework. The self-consistent field theory provides such a theoretical framework, which is capable of accurately predicting mechanical parameters of self-assembled bilayer membranes. In this mini review we summarize the formulation of the self-consistent field theory, as exemplified by a model system composed of flexible amphiphilic chains dissolved in hydrophilic polymeric solvents, and its application to the study of self-assembled bilayer membranes.
We use coarse grained molecular dynamics simulations to investigate diffusion properties of sheared lipid membranes with embedded transmembrane proteins. In membranes without proteins, we find normal in-plane diffusion of lipids in all flow conditions. Protein embedded membranes behave quite differently: by imposing a simple shear flow and sliding the monolayers of the membrane over each other, the motion of protein clusters becomes strongly superdiffusive in the shear direction. In such a circumstance, subdiffusion regime is predominant perpendicular to the flow. We show that superdiffusion is a result of accelerated chaotic motions of protein--lipid complexes within the membrane voids, which are generated by hydrophobic mismatch or the transport of lipids by proteins.