Probing the elastic response of lipid bilayers and nanovesicles to leaflet tensions via volume per lipid.

Biological and biomimetic membranes are based on lipid bilayers, consisting of two monolayers or leaflets. One important but challenging physical parameter of these membranes is their tension. For a long time, this tension was explicitly or implicitly taken to be the bilayer tension, acting on the whole bilayer membrane. More recently, it has been realized that it is useful to decompose the bilayer tension into two leaflet tensions and that these tensions are accessible to molecular dynamics simulations. To divide the bilayer up into two leaflets, it is necessary to introduce a midsurface that defines the spatial extent of the two leaflets. In previous studies, this midsurface was obtained from the density profiles across the bilayer and was then used to compute the molecular area per lipid. Here, we develop an alternative approach based on three-dimensional Voronoi tessellation and molecular volume per lipid. Using this volume-based approach, we determine the reference states with tensionless leaflets as well as the optimal volumes and areas per lipid. The optimal lipid volumes have practically the same value in both leaflets, irrespective of the size and curvature of the nanovesicles, whereas the optimal lipid areas are different for the two leaflets and depend on the vesicle size. In addition, we introduce lateral volume compressibilities to describe the elastic response of the lipid volume to the leaflet tensions. We show that the outer leaflet of a nanovesicle is more densely packed and less compressible than the inner leaflet and that this difference becomes more pronounced for smaller vesicles.

Dimeric states of transmembrane domains of insulin and IGF-1R receptors: Structures and possible role in activation.

Despite the biological significance of insulin signaling, the molecular mechanisms of activation of the insulin receptor (IR) and other proteins from its family remain elusive. Current hypothesis on signal transduction suggests ligand-triggered structural changes in the extracellular domain followed by transmembrane (TM) domains closure and dimerization leading to trans-autophosphorylation and kinase activity in intracellular segments of the receptor. Using NMR spectroscopy, we detected dimerization of isolated TM segments of IR in different membrane-mimicking environments and observed multiple signals of NH groups of protein backbone possibly corresponding to several dimer conformations. Taking available experimental data as constraints, several atomistic models of dimeric TM domains of IR and insulin-like growth factor 1 (IGF-1R) receptors were elaborated. Molecular dynamics simulations of IR ectodomain revealed noticeable collective movements potentially responsible for closure of the C-termini of FnIII-3 domains and spatial approaching of TM helices upon insulin-induced receptor activation. In addition, we demonstrated that the intracellular part of the receptor does not impose restrictions on the positioning of TM helices in the membrane. Finally, we used two independent structure prediction methods to generate a series of dimer conformations followed by their cluster analysis and dimerization free energy estimation to select the best dimer models. Biological relevance of the later was further tested via comparison of the hydrophobic organization of TM helices for both wild-type receptors and their mutants. Based on these data, the ability of several segments from other proteins to functionally replace IR and/or IGF-1R TM domains was explained.