Cell membranes contain hundreds of different proteins and lipids in an asymmetric arrangement. Our current understanding of the detailed organization of cell membranes remains rather elusive, because of the challenge to study fluctuating nanoscale assemblies of lipids and proteins with the required spatiotemporal resolution. Here, we use molecular dynamics simulations to characterize the lipid environment of 10 different membrane proteins. To provide a realistic lipid environment, the proteins are embedded in a model plasma membrane, where more than 60 lipid species are represented, asymmetrically distributed between the leaflets. The simulations detail how each protein modulates its local lipid environment in a unique way, through enrichment or depletion of specific lipid components, resulting in thickness and curvature gradients. Our results provide a molecular glimpse of the complexity of lipid–protein interactions, with potentially far-reaching implications for our understanding of the overall organization of real cell membranes.

No abstract

Our results provide a molecular glimpse of the complexity of lipid-protein interactions, with po-27 tentially far reaching implications for the overall organization of the cell membrane.

Gangliosides are glycolipids in which an oligosaccharide headgroup containing one or more sialic acids is connected to a ceramide. Gangliosides reside in the outer leaflet of the plasma membrane and play a crucial role in various physiological processes such as cell signal transduction and neuronal differentiation by modulating structures and functions of membrane proteins. Because the detailed behavior of gangliosides and protein-ganglioside interactions are poorly known, we investigated the interactions between the gangliosides GM1 and GM3 and the proteins aquaporin (AQP1) and WALP23 using equilibrium molecular dynamics simulations and potential of mean force calculations at both coarse-grained (CG) and atomistic levels. In atomistic simulations, on the basis of the GROMOS force field, ganglioside aggregation appears to be a result of the balance between hydrogen bond interactions and steric hindrance of the headgroups. GM3 clusters are slightly larger and more ordered than GM1 clusters due to the smaller headgroup of GM3. The different structures of GM1 and GM3 clusters from atomistic simulations are not observed at the CG level based on the Martini model, implying a difference in driving forces for ganglioside interactions in atomistic and CG simulations. For protein-ganglioside interactions, in the atomistic simulations, GM1 lipids bind to specific sites on the AQP1 surface, whereas they are depleted from WALP23. In the CG simulations, the ganglioside binding sites on the AQP1 surface are similar, but ganglioside aggregation and protein-ganglioside interactions are more prevalent than in the atomistic simulations. Using the polarizable Martini water model, results were closer to the atomistic simulations. Although experimental data for validation is lacking, we proposed modified Martini parameters for gangliosides to more closely mimic the sizes and structures of ganglioside clusters observed at the atomistic level.

Understanding the lateral organization in plasma membranes remains an open problem despite a large body of research. Model membranes with coexisting micrometer-size domains are routinely employed as simplified models of plasma membranes. Many molecular dynamics simulations have investigated phase separation in model membranes at the coarse-grained level, but atomistic simulations remain computationally challenging. We simulate DPPC:DOPC and DPPC:DOPC:cholesterol lipid bilayers to investigate phase transitions at temperatures from 310 to 270 K. In this temperature range, the binary mixture forms a liquid phase (Lα) and a coexistence of Lα and either gel or ripple phases. The ternary mixture forms a liquid disordered (Ld) phase and a coexistence of liquid ordered (Lo) and either Ld or gel phases. We quantify the coexisting phases and discuss their properties against the background of experimental results. We observe partial registration of growing domains in both mixtures. We characterize specific cholesterol−cholesterol and cholesterol−phospholipid interaction geometries underlying its increased partitioning and the smoothed phase transition in the ternary mixture compared to the binary mixture. By comparing coexisting domains with homogeneous bilayers of the same composition, we demonstrate how domain coexistence affects their properties. Our simulations provide important insights into the lipid−lipid interactions in model lipid bilayers and improve our understanding of the lateral organization in plasma membranes with higher compositional complexity.

The pharmacological activity of a series of novel amide derivatives was characterized on several nicotinic acetylcholine receptors (AChRs). Ca(2+) influx results indicate that these compounds are not agonists of the human (h) α4β2, α3β4, α7, and α1β1γδ AChRs; compounds 2-4 are specific positive allosteric modulators (PAMs) of hα7 AChRs, whereas compounds 1-4, 7, and 12 are noncompetitive antagonists of the other AChRs. Radioligand binding results indicate that PAMs do not inhibit binding of [(3)H]methyllycaconitine but enhance binding of [(3)H]epibatidine to hα7 AChRs, indicating that these compounds do not directly, but allosterically, interact with the hα7 agonist sites. Additional competition binding results indicate that the antagonistic action mediated by these compounds is produced by direct interaction with neither the phencyclidine site in the Torpedo AChR ion channel nor the imipramine and the agonist sites in the hα4β2 and hα3β4 AChRs. Molecular dynamics and docking results suggest that the binding site for compounds 2-4 is mainly located in the inner β-sheet of the hα7-α7 interface, ∼12 Å from the agonist locus. Hydrogen bond interactions between the amide group of the PAMs and the hα7 AChR binding site are found to be critical for their activity. The dual PAM and antagonistic activities elicited by compounds 2-4 might be therapeutically important.

Cholesterol is the most abundant molecule in the plasma membrane of mammals. Its distribution across the two membrane leaflets is critical for understanding how cells work. Cholesterol trans-bilayer motion (flip-flop) is a key process influencing its distribution in membranes. Despite extensive investigations, the rate of cholesterol flip-flop and its dependence on the lateral heterogeneity of membranes remain uncertain. In this work, we used atomistic molecular dynamics simulations to sample spontaneous cholesterol flip-flop events in a DPPC:DOPC:cholesterol mixture with heterogeneous lateral distribution of lipids. In addition to an overall flip-flop rate at the time scale of submilliseconds, we identified a significant impact of local environment on flip-flop rate. We discuss the atomistic details of the flip-flop events observed in our simulations.

Alzheimer's disease (AD) is one of the most common dementia. The aggregation and deposition of the amyloid-β peptide (Aβ) in neural tissue is its characteristic symptom. To destabilize and dissolve Aβ fibrils, a number of candidate molecules have been proposed. wgx-50 is a compound extracted from Sichuan pepper (Zanthoxylum bungeanum) and a potential candidate drug for treating AD. Our early experiments show it is effective in disassembling Aβ42 aggregations. A series of molecular dynamics simulations were performed in this work to explain the molecular mechanism of the destabilization of Aβ42 protofibril by wgx-50. It is found that there were three possible stable binding sites including two sites in hydrophobic grooves on surface of Aβ protofibril that made no significant changes in Aβ structures and one site in the interior that caused destabilization of the protofibril. In this site, wgx-50 was packed against the side chains of I32 and L34, disrupted the D23-K28 salt bridges, and partially opened the tightly compacted two β-sheets. The results were confirmed by simulations at 320 K, where deeper insertion of wgx-50 into the whole protofibril was observed. The molecular mechanism of this novel drug candidate wgx-50 to disaggregate Aβ protofibril may provide some insight into the strategy of structure-based drug design for AD.