Search

In this work we conduct a systematic analysis of the influence of curvature on various properties of a realistic model of mammalian plasma membrane with asymmetric lipid content of monolayers and a realistic concentration of cholesterol. In order to do this we developed new technique for simulating membranes with the global membrane curvature restricted to any desirable value while keeping free lateral diffusion of the lipids and without introducing artifacts or perturbing the membrane structure. We show that the thickness of the hydrophobic core of the concave monolayer decreases by approximately 1.3 Å in comparison to that of the flat membrane, while the thickness of the convex monolayer does not change. The order parameter of the lipid tails decreases significantly in the certain layers of the curved membrane. The area per lipid increases in the convex monolayer and decreases in the concave monolayer. The cholesterol inclination angle decreases when the curvature of a particular monolayer changes from concave to convex. The amount of cholesterol in the minor fraction located between the membrane leaflets is zero in the membrane with positive curvature and increases to 1.7% in the flat membrane and to 2.5% in the membrane with negative curvature.

Asymmetric lipid composition of the cell membranes plays an important role in the multitude of important biological functions. Much less is known, however, about the distribution and dynamics of cholesterol in asymmetric biological membranes. In this work we show how this issue could be addressed computationally by molecular dynamics simulations. The influence of the lipid head group charge, acyl chain saturation, spontaneous membrane curvature and the surface tension of the membrane on cholesterol distribution in asymmetric lipid bilayers is investigated. Four asymmetric bilayers containing DOPC, DOPS, DSPC, or DSPS lipids, were simulated on the time scale extended to tens of microseconds. We show that cholesterol strongly prefers anionic lipids to neutral and saturated lipid tails to unsaturated with distribution ratio ~0.4-0.6. Multiple flip-flop transitions of cholesterol were observed directly and their mean times range from 350 to 2,000 ns. It was shown that the distribution of cholesterol in the asymmetric bilayer depends not only on the type of lipids but also on the local membrane curvature and the surface tension. The geometric shape of spontaneously curved asymmetric bilayer changes dramatically in the presence of cholesterol. The membrane curvature becomes less homogeneous with large patches of flattened regions interleaved by rather sharp bends.

Transmembrane translocation of C60 fullerenes functionalized by the single amino-derivative in neutral and charged forms was studies by extensive all-atom molecular dynamics simulations. It is shown that these complexes exhibit very strong affinity to the membrane core, but their spontaneous translocation through the membrane is not possible at practical time scale. In contrast, free amino derivatives translocate through the membrane much easier than their complexes with fullerenes, but do not have pronounced affinity to the membrane interior. Our results suggest that monofunctionalized C60 could be extremely efficient membrane targeting agents, which facilitate accumulation of the water-soluble compounds in the hydrophobic core of lipid bilayer.

The Hierarchical Domain-Wise Alignment (HDWA) technique of domain identification in proteins is presented. HDWA is designed to identify hierarchically organized dynamic domains in proteins using the MD trajectories by eliminating systematic motions from MD trajectories recursively in a model-free manner. The method is tested on the proteins from different structural classes.

The major challenge in setting up membrane protein simulations is embedding the protein into the pre-equilibrated lipid bilayer. Several techniques were proposed to achieve optimal packing of the lipid molecules around the protein. However, all of them possess serious disadvantages, which limit their applicability and discourage the users of simulation packages from using them. In the present work, we analyzed existing approaches and proposed a new procedure of protein insertion into the lipid bilayer, which is implemented in the ProtSqueeze software. The advantages of ProtSqueeze are as follows: (1) the insertion algorithm is simple, understandable, and controllable; (2) the software can work with virtually any simulation package on virtually any platform; (3) no modification of the source code of the simulation package is needed; (4) the procedure of insertion is as automated as possible; (5) ProtSqueeze is distributed for free under a general public license. In this work, we present the architecture and the algorithm of ProtSqueeze and demonstrate its usage in case studies.