Search

The human Niemann-Pick C1 (NPC1) gene encoding a 1278 amino acid protein is very heterogeneous. While some variants represent benign polymorphisms, NPC disease carriers and patients may possess rare variants, whose functional importance remains unknown. An NPC1 cDNA construct known as NPC1 wild-type variant (WT-V), distributed between laboratories and used as a WT control in several studies, also contains changes regarding specific amino acids compared to the NPC1 Genbank reference sequence. To improve the dissection of subtle functional differences, we generated human cells stably expressing NPC1 variants from the AAVS1 safe-harbor locus on an NPC1-null background engineered by CRISPR/Cas9 editing. We then employed high-content imaging with automated image analysis to quantitatively assess LDL-induced, time-dependent changes in lysosomal cholesterol content and lipid droplet formation. Our results indicate that the L472P change present in NPC1 WT-V compromises NPC1 functionality in lysosomal cholesterol export. All-atom molecular dynamics simulations suggest that the L472P change alters the relative position of the NPC1 middle and the C-terminal luminal domains, disrupting the recently characterized cholesterol efflux tunnel. These results reveal functional defects in NPC1 WT-V and highlight the strength of simulations and quantitative imaging upon stable protein expression in elucidating subtle differences in protein function. publishedVersion

HFD-bound and HFD-free simulation systems presented in the paper; see table_of_simulations.pdf Contents of the simulation_archive: _FORCEFIELD: Force field parameters and topologies in GROMACS format amber-ff14SB.ff: The required force field parameters compiled for all systems.*** actin: chicken actin | topology (mol.itp)***, structure (mol.gro) ncap: mouse HFD domain (CAP) | topology (mol.itp)***, structure (mol.gro) ADP: ADP | topology (mol.itp)***, structure (mol.gro) MMG: multisite Mg2+ | topology (mol.itp)***, structure (mol.gro) cofilin: rabbit cofilin | topology (mol.itp)***, structure (mol.gro) dock: the "platform" region | topology (mol.itp), structure (mol.gro), the position restraints (position_restraints.itp) _parametrize_HIC: Amber format parameters for methyl-histidine (Nτ-Methyl-L-histidine, HIC) to be used together with amber-ff14SB force field. Charges were calculated based on the protocol by Duan, Y. et al. A point-charge force field for molecular mechanics simulations of proteins based on condensed-phase quantum mechanical calculations. J. Comput. Chem. 24, 1999–2012 (2003). Missing bonded parameters were adopted from GAFF2. *** Converted to GROMACS from amber format using the ParmEd tool (https://github.com/ParmEd/ParmEd) HFD-bound: contains the directories for each simulation repeat listed in the table of simulations. 1, 2, and 3: simulations that were started by replacing the terminal actins with HFD-bound actins. 1' and 2': simulations that were started by docking HFD after some simulation without HFD (see below HFD-free simulations 1 and 2) HFD-free: contains the directories for each simulation repeat listed in the table of simulations. 1, 2, and 3: simulations of the cofilin-decorated actin filament. 2': simulations that was started by removing HFD from the HFD-bound system (see above HFD-bound simulation 2) Each subdirectory (under 1 2 3, etc.) contains prod.mdp: simulation setting file topol.top: GROMACS topology file t0.pdb: initial coordinates for simulation obtained after a set of restrained NVT and NpT equilibration steps prod0.tpr: the GROMACS tpr file index.ndx: index file nW.pdb: t0.pdb with all water molecules removed nW.xtc: the simulation trajectory with all water molecules removed. _FORCEFIELD: link to ../../_FORCEFIELD, {"references":["Kotila, T., Wioland, H., Enkavi, G. et al. Mechanism of synergistic actin filament pointed end depolymerization by cyclase-associated protein and cofilin. Nat Commun 10, 5320 (2019) doi:10.1038/s41467-019-13213-2"]}

Lipid nanodiscs are macromolecular assemblies, where a scaffold protein is wrapped around a nanosized disc of a lipid bilayer, thus protecting the hydrocarbon chains at the disc edges from unfavorable interactions with water. These nanostructures have numerous applications in, e.g., nanotechnology and pharmaceutics, and in investigations of membrane proteins. Here, we present results based on atomistic molecular dynamics simulations combined with electron paramagnetic spectroscopy measurements on the structure and dynamics of lipids in single-component nanodiscs. Our data highlight the existence of three distinctly different lipid fractions: central lipids residing in the center of a nanodisc, boundary lipids in direct contact with a scaffold protein, and intermediate lipids between these two regions. The central lipids are highly ordered and characterized by slow diffusion. In this part of the nanodisc, the membrane is the thickest and characterized by a gel-like or liquid-ordered phase, having features common to cholesterol-rich membranes. The boundary lipids in direct contact with the scaffold protein turned out to be less ordered and characterized by faster diffusion, and they remained in the liquid-disordered phase even at temperatures that were somewhat below the main phase transition temperature (Tm). The enthalpies associated with the central-boundary and central-intermediate transitions were similar to those observed for lipids going through the main phase transition. Overall, the study reveals lipid nanodiscs to be characterized by a complex internal structure, which is expected to influence membrane proteins placed in nanodiscs. publishedVersion

Cholesterol renders mammalian cell membranes more compact by reducing the amount of voids in the membrane structure. Because of this, cholesterol is known to regulate the ability of cell membranes to prevent the permeation of water and water-soluble molecules through the membranes. Meanwhile, it is also known that even seemingly tiny modifications in the chemical structure of cholesterol can lead to notable changes in membrane properties. The question is, how significantly do these small changes in cholesterol structure affect the permeability barrier function of cell membranes? In this work, we applied fluorescence methods as well as atomistic molecular dynamics simulations to characterize changes in lipid membrane permeability induced by cholesterol oxidation. The studied 7β-hydroxycholesterol (7β-OH-chol) and 27-hydroxycholesterol (27-OH-chol) represent two distinct groups of oxysterols, namely, ring- and tail-oxidized cholesterols, respectively. Our previous research showed that the oxidation of the cholesterol tail has only a marginal effect on the structure of a lipid bilayer; however, oxidation was found to disturb membrane dynamics by introducing a mechanism that allows sterol molecules to move rapidly back and forth across the membrane-bobbing. Herein, we show that bobbing of 27-OH-chol accelerates fluorescence quenching of NBD-lipid probes in the inner leaflet of liposomes by dithionite added to the liposomal suspension. Systematic experiments using fluorescence quenching spectroscopy and microscopy led to the conclusion that the presence of 27-OH-chol increases membrane permeability to the dithionite anion. Atomistic molecular dynamics simulations demonstrated that 27-OH-chol also facilitates water transport across the membrane. The results support the view that oxysterol bobbing gives rise to successive perturbations to the hydrophobic core of the membrane, and these perturbations promote the permeation of water and small water-soluble molecules through a lipid bilayer. The observed impairment of permeability can have important consequences for eukaryotic organisms. The effects described for 27-OH-chol were not observed for 7β-OH-chol which represents ring-oxidized sterols. publishedVersion

Synaptic neurotransmission has recently been proposed to function via either a membrane-independent or a membrane-dependent mechanism, depending on the neurotransmitter type. In the membrane-dependent mechanism, amphipathic neurotransmitters first partition to the lipid headgroup region and then diffuse along the membrane plane to their membrane-buried receptors. However, to date, this mechanism has not been demonstrated for any neurotransmitter-receptor complex. Here, we combined isothermal calorimetry measurements with a diverse set of molecular dynamics simulation methods to investigate the partitioning of an amphipathic neurotransmitter (dopamine) and the mechanism of its entry into the ligand-binding site. Our results show that the binding of dopamine to its receptor is consistent with the membrane-dependent binding and entry mechanism. Both experimental and simulation results showed that dopamine favors binding to lipid membranes especially in the headgroup region. Moreover, our simulations revealed a ligand-entry pathway from the membrane to the binding site. This pathway passes through a lateral gate between transmembrane alpha-helices 5 and 6 on the membrane-facing side of the protein. All in all, our results demonstrate that dopamine binds to its receptor by a membrane-dependent mechanism, and this is complemented by the more traditional binding mechanism directly through the aqueous phase. The results suggest that the membrane-dependent mechanism is common in other synaptic receptors, too. publishedVersion

Myelin protein P2 is a peripheral membrane protein of the fatty acid-binding protein family that functions in the formation and maintenance of the peripheral nerve myelin sheath. Several P2 gene mutations cause human Charcot-Marie-Tooth neuropathy, but the mature myelin sheath assembly mechanism is unclear. Here, cryo-EM of myelin-like proteolipid multilayers revealed an ordered three-dimensional (3D) lattice of P2 molecules between stacked lipid bilayers, visualizing supramolecular assembly at the myelin major dense line. The data disclosed that a single P2 layer is inserted between two bilayers in a tight intermembrane space of ∼3 nm, implying direct interactions between P2 and two membrane surfaces. X-ray diffraction from P2-stacked bicelle multilayers revealed lateral protein organization, and surface mutagenesis of P2 coupled with structure-function experiments revealed a role for both the portal region of P2 and its opposite face in membrane interactions. Atomistic molecular dynamics simulations of P2 on model membrane surfaces suggested that Arg-88 is critical for P2-membrane interactions, in addition to the helical lid domain. Negatively charged lipid headgroups stably anchored P2 on the myelin-like bilayer surface. Membrane binding may be accompanied by opening of the P2 β-barrel structure and ligand exchange with the apposing bilayer. Our results provide an unprecedented view into an ordered, multilayered biomolecular membrane system induced by the presence of a peripheral membrane protein from human myelin. This is an important step toward deciphering the 3D assembly of a mature myelin sheath at the molecular level. publishedVersion

FGF2 is a tumor cell survival factor that is exported from cells by an ER/Golgi-independent secretory pathway. This unconventional mechanism of protein secretion is based on direct translocation of FGF2 across the plasma membrane. The Na,K-ATPase has previously been shown to play a role in this process, however, the underlying mechanism has remained elusive. Here, we define structural elements that are critical for a direct physical interaction between FGF2 and the α1 subunit of the Na,K-ATPase. In intact cells, corresponding FGF2 mutant forms were impaired regarding both recruitment at the inner plasma membrane leaflet and secretion. Ouabain, a drug that inhibits both the Na,K-ATPase and FGF2 secretion, was found to impair the interaction of FGF2 with the Na,K-ATPase in cells. Our findings reveal the Na,K-ATPase as the initial recruitment factor for FGF2 at the inner plasma membrane leaflet being required for efficient membrane translocation of FGF2 to cell surfaces., Legrand et al. identify two lysine residues on molecular surface of Fibroblast Growth Factor 2 (FGF2) essential for its interaction with α1 subunit of the Na,K-ATPase. They further conclude that this interaction precedes interaction of the FGF2 with PI(4,5)P2 and facilitates its unconventional secretion across the membrane, which is impaired by Ouabain, an Na,K-ATPase inhibitor.

The cryo-EM structure of the complete insulin receptor ectodomain saturated with four insulin ligands reveals a T-like conformation with converging membrane-proximal domains and structural evidence of insulin binding sites 2/2′., Glucose homeostasis and growth essentially depend on the hormone insulin engaging its receptor. Despite biochemical and structural advances, a fundamental contradiction has persisted in the current understanding of insulin ligand–receptor interactions. While biochemistry predicts two distinct insulin binding sites, 1 and 2, recent structural analyses have resolved only site 1. Using a combined approach of cryo-EM and atomistic molecular dynamics simulation, we present the structure of the entire dimeric insulin receptor ectodomain saturated with four insulin molecules. Complementing the previously described insulin–site 1 interaction, we present the first view of insulin bound to the discrete insulin receptor site 2. Insulin binding stabilizes the receptor ectodomain in a T-shaped conformation wherein the membrane-proximal domains converge and contact each other. These findings expand the current models of insulin binding to its receptor and of its regulation. In summary, we provide the structural basis for a comprehensive description of ligand–receptor interactions that ultimately will inform new approaches to structure-based drug design., Graphical Abstract

Surfactant protein B (SP-B) is essential in transferring surface-active phospholipids from membrane-based surfactant complexes into the alveolar air-liquid interface. This allows maintaining the mechanical stability of the surfactant film under high pressure at the end of expiration; therefore, SP-B is crucial in lung function. Despite its necessity, the structure and the mechanism of lipid transfer by SP-B have remained poorly characterized. Earlier, we proposed higher-order oligomerization of SP-B into ring-like supramolecular assemblies. In the present work, we used coarse-grained molecular dynamics simulations to elucidate how the ring-like oligomeric structure of SP-B determines its membrane binding and lipid transfer. In particular, we explored how SP-B interacts with specific surfactant lipids, and how consequently SP-B reorganizes its lipid environment to modulate the pulmonary surfactant structure and function. Based on these studies, there are specific lipid-protein interactions leading to perturbation and reorganization of pulmonary surfactant layers. Especially, we found compelling evidence that anionic phospholipids and cholesterol are needed or even crucial in the membrane binding and lipid transfer function of SP-B. Also, on the basis of the simulations, larger oligomers of SP-B catalyze lipid transfer between adjacent surfactant layers. Better understanding of the molecular mechanism of SP-B will help in the design of therapeutic SP-B-based preparations and novel treatments for fatal respiratory complications, such as the acute respiratory distress syndrome. (C) 2020 The Author(s). Published by Elsevier Ltd.

Toll-like receptor 4 (TLR4) is activated by bacterial lipopolysaccharide (LPS), which drives the production of proinflammatory cytokines. Earlier studies have indicated that cholesterol- and glycosphingolipid-rich subregions of the plasma membrane (lipid domains) are important for TLR4-mediated signaling. We report that inhibition of glucosylceramide (GluCer) synthase, which resulted in decreased concentrations of the glycosphingolipid GluCer in lipid domains, reduced the LPS-induced inflammatory response in both mouse and human macrophages. Atomistic molecular dynamics simulations of the TLR4 dimer complex (with and without LPS in its MD-2 binding pockets) in membranes (in the presence and absence of GluCer) showed that: (1) LPS induced a tilted orientation of TLR4 and increased dimer integrity; (2) GluCer did not affect the integrity of the LPS/TLR4 dimer but reduced the LPS-induced tilt; and (3) GluCer increased electrostatic interactions between the membrane and the TLR4 extracellular domain, which could potentially modulate the tilt. We also showed that GCS inhibition reduced the interaction between TLR4 and the intracellular adaptor protein Mal. We conclude that the GluCer-induced effects on LPS/TLR4 orientation may influence the signaling capabilities of the LPS/TLR4 complex by affecting its interaction with downstream signaling proteins. publishedVersion

Accurately calculating rate constants of macroscopic chemical processes from molecular dynamics simulations is a long-sought but elusive goal. The problem is particularly relevant for processes occurring in biological systems, as is the case for ligand-protein and ligand-membrane interactions. Several formalisms to determine rate constants from easily accessible free-energy profiles [Δ Go( z)] of a molecule along a coordinate of interest have been proposed. However, their applicability for molecular interactions in condensed media has not been critically evaluated or validated. This work presents such evaluation and validation and introduces improved methodology. As a case study, we have characterized quantitatively the rate of translocation of cholesterol across 1-palmitoyl-2-oleoyl- sn-glycero-3-phosphocholine bilayers. Translocation across lipid bilayers is the rate-limiting step in the permeation of most drugs through biomembranes. We use coarse-grained molecular dynamics simulations and different kinetic formalisms to calculate this rate constant. A self-consistent test of the applicability of various available formalisms is provided by comparing their predictions with the translocation rates obtained from actual events observed in long unrestrained simulations. To this effect, a novel procedure was used to obtain the effective rate constant, based on an analysis of time intervals between transitions among different states along the reaction coordinate. While most tested formalisms lead to results in reasonable agreement (within a factor of 5) with this effective rate constant, the most adequate one is based on the explicit relaxation frequencies from the transition state in the forward and backward directions along the reaction coordinate. acceptedVersion

Folding and packing of membrane proteins are highly influenced by the lipidic component of the membrane. Here, we explore how the hydrophobic mismatch (the difference between the hydrophobic span of a transmembrane protein region and the hydrophobic thickness of the lipid membrane around the protein) influences transmembrane helix packing in a cellular environment. Using a ToxRED assay in Escherichia coli and a Bimolecular Fluorescent Complementation approach in human-derived cells complemented by atomistic molecular dynamics simulations we analyzed the dimerization of Glycophorin A derived transmembrane segments. We concluded that, biological membranes can accommodate transmembrane homo-dimers with a wide range of hydrophobic lengths. Hydrophobic mismatch and its effects on dimerization are found to be considerably weaker than those previously observed in model membranes, or under in vitro conditions, indicating that biological membranes (particularly eukaryotic membranes) can adapt to structural deformations through compensatory mechanisms that emerge from their complex structure and composition to alleviate membrane stress. Results based on atomistic simulations support this view, as they revealed that Glycophorin A dimers remain stable, despite of poor hydrophobic match, using mechanisms based on dimer tilting or local membrane thickness perturbations. Furthermore, hetero-dimers with large length disparity between their monomers are also tolerated in cells, and the conclusions that one can draw are essentially similar to those found with homo-dimers. However, large differences between transmembrane helices length hinder the monomer/dimer equilibrium, confirming that, the hydrophobic mismatch has, nonetheless, biologically relevant effects on helix packing in vivo. publishedVersion publishedVersion

Biological membranes are tricky to investigate. They are complex in terms of molecular composition and structure, functional over a wide range of time scales, and characterized by nonequilibrium conditions. Because of all of these features, simulations are a great technique to study biomembrane behavior. A significant part of the functional processes in biological membranes takes place at the molecular level; thus computer simulations are the method of choice to explore how their properties emerge from specific molecular features and how the interplay among the numerous molecules gives rise to function over spatial and time scales larger than the molecular ones. In this review, we focus on this broad theme. We discuss the current state-of-the-art of biomembrane simulations that, until now, have largely focused on a rather narrow picture of the complexity of the membranes. Given this, we also discuss the challenges that we should unravel in the foreseeable future. Numerous features such as the actin-cytoskeleton network, the glycocalyx network, and nonequilibrium transport under ATP-driven conditions have so far received very little attention; however, the potential of simulations to solve them would be exceptionally high. A major milestone for this research would be that one day we could say that computer simulations genuinely research biological membranes, not just lipid bilayers. publishedVersion

Understanding the molecular mechanisms governing nanoparticle–membrane interactions is of prime importance for drug delivery and biomedical applications. Neutron reflectometry (NR) experiments are combined with atomistic and coarse‐grained molecular dynamics (MD) simulations to study the interaction between cationic gold nanoparticles (AuNPs) and model lipid membranes composed of a mixture of zwitterionic di‐stearoyl‐phosphatidylcholine (DSPC) and anionic di‐stearoyl‐phosphatidylglycerol (DSPG). MD simulations show that the interaction between AuNPs and a pure DSPC lipid bilayer is modulated by a free energy barrier. This can be overcome by increasing temperature, which promotes an irreversible AuNP incorporation into the lipid bilayer. NR experiments confirm the encapsulation of the AuNPs within the lipid bilayer at temperatures around 55 °C. In contrast, the AuNP adsorption is weak and impaired by heating for a DSPC–DSPG (3:1) lipid bilayer. These results demonstrate that both the lipid charge and the temperature play pivotal roles in AuNP–membrane interactions. Furthermore, NR experiments indicate that the (negative) DSPG lipids are associated with lipid extraction upon AuNP adsorption, which is confirmed by coarse‐grained MD simulations as a lipid‐crawling effect driving further AuNP aggregation. Overall, the obtained detailed molecular view of the interaction mechanisms sheds light on AuNP incorporation and membrane destabilization. © 2019 The Authors. Published by WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim This is an open access article under the terms of the Creative Commons Attribution‐NonCommercial‐NoDerivs License, which permits use and distribution in any medium, provided the original work is properly cited, the use is non‐commercial and no modifications or adaptations are made.

A mutated KRAS protein is frequently observed in human cancers. Traditionally, the oncogenic properties of KRAS missense mutants at position 12 (G12X) have been considered as equal. Here, by assessing the probabilities of occurrence of all KRAS G12X mutations and KRAS dynamics we show that this assumption does not hold true. Instead, our findings revealed an outstanding mutational bias. We conducted a thorough mutational analysis of KRAS G12X mutations and assessed to what extent the observed mutation frequencies follow a random distribution. Unique tissue-specific frequencies are displayed with specific mutations, especially with G12R, which cannot be explained by random probabilities. To clarify the underlying causes for the nonrandom probabilities, we conducted extensive atomistic molecular dynamics simulations (170 μs) to study the differences of G12X mutations on a molecular level. The simulations revealed an allosteric hydrophobic signaling network in KRAS, and that protein dynamics is altered among the G12X mutants and as such differs from the wild-type and is mutation-specific. The shift in long-timescale conformational dynamics was confirmed with Markov state modeling. A G12X mutation was found to modify KRAS dynamics in an allosteric way, which is especially manifested in the switch regions that are responsible for the effector protein binding. The findings provide a basis to understand better the oncogenic properties of KRAS G12X mutants and the consequences of the observed nonrandom frequencies of specific G12X mutations., Author summary The oncogene KRAS is frequently mutated in various cancers. When the amino acid glycine 12 is mutated, KRAS protein acquires oncogenic properties that result in tumor cell-growth and cancer progression. These mutations prevail especially in the pancreatic ductal adenocarcinoma, which is a cancer with an exceptionally dismal prognosis. To date, there is a limited understanding of the different mutations at the position 12, also regarding whether the different mutations would have different consequences. These discrepancies could have major implications for the future drug therapies targeting KRAS mutant harboring tumors. In this study, we made a critical assessment of the observed frequency of KRAS G12X mutations and the underlying causes for these frequencies. We also assessed KRAS G12X mutant discrepancies on an atomistic level by utilizing state-of-the-art molecular dynamics simulations. We found that the dynamics of the mutants does not only differ from the wild-type protein, but there is also a profound difference among the different mutants. These results emphasize that the different KRAS G12X mutations are not equal, and thereby they suggest that the future research related to mutant KRAS biology should account for these observations.

Cardiolipins have multiple vital functions within biological cell membranes, most notably in the energy metabolism associated with the inner mitochondrial membrane. Considering their essential role, peroxidation of cardiolipins may plausibly have significant effects, as peroxidation is known to alter the functionality of lipid molecules. We used atomistic molecular dynamics simulations to study how peroxidation of cardiolipin affects the properties of the inner mitochondrial membrane. To this end, we explored what happens when varying fractions of fatty acid chains of cardiolipin are replaced by its four different oxidized products in systems modeling the inner mitochondrial membrane. We found that the oxidation of cardiolipin leads to a conformational change both in the backbone/head group and in chain regions of oxidized cardiolipin molecules. The oxidized groups were observed to shift closer to the membrane-water interface region, where they formed hydrogen bonds with several other groups. Additionally, the conformational change turned out to decrease bilayer thickness, and to increase the area per lipid chain, though these changes were minor. The acyl chain conformational order of unoxidized lipids exposed to interactions with oxidized cardiolipins was increased in carbons 3–5 and decreased in carbons 13–17 due to the structural reorganization of the cardiolipin molecules. Overall, the results bring up that the conformation of cardiolipin is altered upon oxidation, suggesting that its oxidation may interfere its interactions with mitochondrial proteins and thereby affect cardiolipin-dependent cellular processes such as electron and proton transport. publishedVersion

Lung surfactant and a tear film lipid layer are examples of biologically relevant macromolecular structures found at the air–water interface. Because of their complexity, they are often studied in terms of simplified lipid layers, the simplest example being a Langmuir monolayer. Given the profound biological significance of these lipid assemblies, there is a need to understand their structure and dynamics on the nanoscale, yet there are not many techniques able to provide this information. Atomistic molecular dynamics simulations would be a tool fit for this purpose; however, the simulation models suggested until now have been qualitative instead of quantitative. This limitation has mainly stemmed from the challenge to correctly describe the surface tension of water with simulation parameters compatible with other biomolecules. In this work, we show that this limitation can be overcome by using the recently introduced four-point OPC water model, whose surface tension for water is demonstrated to be quantitatively consistent with experimental data and which is also shown to be compatible with the commonly employed lipid models. We further establish that the approach of combining the OPC four-point water model with the CHARMM36 lipid force field provides nearly quantitative agreement with experiments for the surface pressure–area isotherm for POPC and DPPC monolayers, also including the experimentally observed phase coexistence in a DPPC monolayer. The simulation models reported in this work pave the way for nearly quantitative atomistic studies of lipid-rich biological structures at air–water interfaces. acceptedVersion

FGF2 is secreted from cells by an unconventional secretory pathway. This process is mediated by direct translocation across the plasma membrane. Here, we define the minimal molecular machinery required for FGF2 membrane translocation in a fully reconstituted inside-out vesicle system. FGF2 membrane translocation is thermodynamically driven by PI(4,5)P2-induced membrane insertion of FGF2 oligomers. The latter serve as dynamic translocation intermediates of FGF2 with a subunit number in the range of 8-12 FGF2 molecules. Vectorial translocation of FGF2 across the membrane is governed by sequential and mutually exclusive interactions with PI(4,5)P2 and heparan sulfates on opposing sides of the membrane. Based on atomistic molecular dynamics simulations, we propose a mechanism that drives PI(4,5)P2 dependent oligomerization of FGF2. Our combined findings establish a novel type of self-sustained protein translocation across membranes revealing the molecular basis of the unconventional secretory pathway of FGF2. publishedVersion publishedVersion

Hyaluronan is a polyanionic, megadalton-scale polysaccharide, which initiates cell signaling by interacting with several receptor proteins including CD44 involved in cell-cell interactions and cell adhesion. Previous studies of the CD44 hyaluronan binding domain have identified multiple widespread residues to be responsible for its recognition capacity. In contrast, the X-ray structural characterization of CD44 has revealed a single binding mode associated with interactions that involve just a fraction of these residues. In this study, we show through atomistic molecular dynamics simulations that hyaluronan can bind CD44 with three topographically different binding modes that in unison define an interaction fingerprint, thus providing a plausible explanation for the disagreement between the earlier studies. Our results confirm that the known crystallographic mode is the strongest of the three binding modes. The other two modes represent metastable configurations that are readily available in the initial stages of the binding, and they are also the most frequently observed modes in our unbiased simulations. We further discuss how CD44, fostered by the weaker binding modes, diffuses along HA when attached. This 1D diffusion combined with the constrained relative orientation of the diffusing proteins is likely to influence the aggregation kinetics of CD44. Importantly, CD44 aggregation has been suggested to be a possible mechanism in CD44-mediated signaling., Author summary Hyaluronan is a natural sugar polymer in our bodies. Besides acting as a space-filling agent for example in multiple connective tissues, it also functions as a cellular cue in cancer and inflammation. Our tissues sense hyaluronan through receptors—proteins that sit at the surface of cells and grab the molecules they are expected to recognize. Although the knowledge associated with hyaluronan and its receptors is constantly accumulating, the molecular-level insight is largely missing or incomplete due to the lack of techniques able to probe the dynamics of protein–carbohydrate interactions with sufficiently high resolution. In this work, we characterize the binding of hyaluronan to its receptor CD44 with atomistic precision. We achieve this level of precision by employing atomistic molecular dynamics simulations. This computational technique allows one to follow the movement of atoms of a virtual system at scales beyond the resolution of any experimental technique. Our work specifically focuses on the different stages of hyaluronan–CD44 binding, and we observe the process to involve three different binding modes, making it more versatile than previously thought. Our insights, therefore, promote the understanding of the interplay between hyaluronan and HA, thereby fostering development of new drugs or inhibitors to malignancies, such as cancer metastasis.

Biological membranes generate specific functions through compartmentalized regions such as cholesterol-enriched membrane nanodomains that host selected proteins. Despite the biological significance of nanodomains, details on their structure remain elusive. They cannot be observed via microscopic experimental techniques due to their small size, yet there is also a lack of atomistic simulation models able to describe spontaneous nanodomain formation in sufficiently simple but biologically relevant complex membranes. Here we use atomistic simulations to consider a binary mixture of saturated dipalmitoylphosphatidylcholine and cholesterol — the “minimal standard” for nanodomain formation. The simulations reveal how cholesterol drives the formation of fluid cholesterol-rich nanodomains hosting hexagonally packed cholesterol-poor lipid nanoclusters, both of which show registration between the membrane leaflets. The complex nanodomain substructure forms when cholesterol positions itself in the domain boundary region. Here cholesterol can also readily flip–flop across the membrane. Most importantly, replacing cholesterol with a sterol characterized by a less asymmetric ring region impairs the emergence of nanodomains. The model considered explains a plethora of controversial experimental results and provides an excellent basis for further computational studies on nanodomains. Furthermore, the results highlight the role of cholesterol as a key player in the modulation of nanodomains for membrane protein function.

Proteins embedded in the plasma membrane mediate interactions with the cell environment and play decisive roles in many signaling events. For cell-cell recognition molecules, it is highly likely that their structures and behavior have been optimized in ways that overcome the limitations of membrane tethering. In particular, the ligand binding regions of these proteins likely need to be maximally exposed. Here we show by means of atomistic simulations of membrane-bound CD2, a small cell adhesion receptor expressed by human T-cells and natural killer cells, that the presentation of its ectodomain is highly dependent on membrane lipids and receptor glycosylation acting in apparent unison. Detailed analysis shows that the underlying mechanism is based on electrostatic interactions complemented by steric interactions between glycans in the protein and the membrane surface. The findings are significant for understanding the factors that render membrane receptors accessible for binding and signaling. publishedVersion

This article has an erratum: Doi 10.1038/s41598-017-18751-7 Charcot-Marie-Tooth (CMT) disease is one of the most common inherited neuropathies. Recently, three CMT1-associated point mutations (I43N, T51P, and I52T) were discovered in the abundant peripheral myelin protein P2. These mutations trigger abnormal myelin structure, leading to reduced nerve conduction velocity, muscle weakness, and distal limb atrophy. P2 is a myelin-specific protein expressed by Schwann cells that binds to fatty acids and membranes, contributing to peripheral myelin lipid homeostasis. We studied the molecular basis of the P2 patient mutations. None of the CMT1-associated mutations alter the overall folding of P2 in the crystal state. P2 disease variants show increased aggregation tendency and remarkably reduced stability, T51P being most severe. In addition, P2 disease mutations affect protein dynamics. Both fatty acid binding by P2 and the kinetics of its membrane interactions are affected by the mutations. Experiments and simulations suggest opening of the beta barrel in T51P, possibly representing a general mechanism in fatty acid-binding proteins. Our findings demonstrate that altered biophysical properties and functional dynamics of P2 may cause myelin defects in CMT1 patients. At the molecular level, a few malformed hydrogen bonds lead to structural instability and misregulation of conformational changes related to ligand exchange and membrane binding.

Niemann-Pick Protein C2 (npc2) is a small soluble protein critical for cholesterol transport within and from the lysosome and the late endosome. Intriguingly, npc2-mediated cholesterol transport has been shown to be modulated by lipids, yet the molecular mechanism of npc2-membrane interactions has remained elusive. Here, based on an extensive set of atomistic simulations and free energy calculations, we clarify the mechanism and energetics of npc2-membrane binding and characterize the roles of physiologically relevant key lipids associated with the binding process. Our results capture in atomistic detail two competitively favorable membrane binding orientations of npc2 with a low interconversion barrier. The first binding mode (Prone) places the cholesterol binding pocket in direct contact with the membrane and is characterized by membrane insertion of a loop (V59-M60-G61-I62-P63-V64-P65). This mode is associated with cholesterol uptake and release. On the other hand, the second mode (Supine) places the cholesterol binding pocket away from the membrane surface, but has overall higher membrane binding affinity. We determined that bis(monoacylglycero)phosphate (bmp) is specifically required for strong membrane binding in Prone mode, and that it cannot be substituted by other anionic lipids. Meanwhile, sphingomyelin counteracts bmp by hindering Prone mode without affecting Supine mode. Our results provide concrete evidence that lipids modulate npc2-mediated cholesterol transport either by favoring or disfavoring Prone mode and that they impose this by modulating the accessibility of bmp for interacting with npc2. Overall, we provide a mechanism by which npc2-mediated cholesterol transport is controlled by the membrane composition and how npc2-lipid interactions can regulate the transport rate., Author summary Cholesterol plays essential structural and functional roles in all vertebrate cells. Abnormalities in cholesterol metabolism are associated with severe and wide-spread diseases. Cholesterol efflux is one of the essential steps in its metabolism and is mediated by Niemann-Pick Protein C2 (npc2) in endosomes and lysosomes. Mutations in npc2 result in the fatal genetic disease, called Niemann-Pick C disease, characterized by neuronal degeneration resulting in early death. We performed an extensive set of atomistic molecular dynamics simulations employing enhanced sampling approaches to investigate how npc2 mediates cholesterol transport and the roles of membrane lipids in the process. We show that the binding of npc2 to endosomal/lysosomal membranes and the consequent trafficking of cholesterol are concertedly regulated by bis(monoacylglycero)phosphate and sphingomyelin. Moreover, our results explain the molecular mechanism of experimentally observed affect of lipids on cholesterol transport rate. In essence, our findings provide key insights into the regulation of the lysosomal and endosomal cholesterol trafficking by npc2 and reveal key lipid-protein interactions.

Driven by interactions between lipids and proteins, biological membranes display lateral heterogeneity that manifests itself in a mosaic of liquid-ordered (Lo) or raft, and liquid-disordered (Ld) or non-raft domains with a wide range of different properties and compositions. In giant plasma membrane vesicles and giant unilamellar vesicles, specific binding of Cholera Toxin (CTxB) to GM1 glycolipids is a commonly used strategy to label raft domains or Lo membrane environments. However, these studies often use acyl-chain labeled bodipy-GM1 (bdGM1), whose headgroup accessibility and membrane order or phase partitioning may differ from those of GM1, rendering the interpretation of CTxB binding data quite problematic. To unravel the molecular basis of CTxB binding to GM1 and bdGM1, we explored the partitioning and the headgroup presentation of these gangliosides in the Lo and Ld phases using atomistic molecular dynamics simulations complemented by CTxB binding experiments. The conformation of both GM1 and bdGM1 was shown to be largely similar in the Lo and Ld phases. However, bdGM1 showed reduction in receptor availability when reconstituted into synthetic bilayer mixtures, highlighting that membrane phase partitioning of the gangliosides plays a considerable role in CTxB binding. Our results suggest that the CTxB binding is predominately modulated by the partitioning of the receptor to an appropriate membrane phase. Further, given that the Lo and Ld partitioning of bdGM1 differs from those of GM1, usage of bdGM1 for studying GM1 behavior in cells can lead to invalid interpretation of experimental data. publishedVersion

Molecular oxygen (O2) is key to all life on earth, as it is constantly cycled via photosynthesis and cellular respiration. Substantial scientific effort has been devoted to understanding every part of this cycle. Classical molecular dynamics (MD) simulations have been used to study some of the key processes involved in cellular respiration: O2 permeation through alveolar monolayers and cellular membranes, its binding to hemoglobin during transport in the bloodstream, as well as its transport along optimal pathways toward its reduction sites in proteins. Moreover, MD simulations can help interpret the results of several imaging techniques in which O2 is used because of its paramagnetic nature. However, despite the widespread use of computational models for the O2 molecule, their performances have never been systematically evaluated. In this paper, we assess the performances of 14 different models of O2 available in the literature by calculating four thermodynamic properties: density, heat of vaporization, free energy of hydration, and free energy of solvation in hexadecane. For each property, reliable experimental data are available. Most models perform reasonably well in predicting the correct trends, but they fail to reproduce the experimental data quantitatively. We then develop new models for O2, with and without a quadrupole moment, and compare their behavior with the behavior of previously published models. The new models show significant improvement in terms of density, heat of vaporization, and free energy of hydration. However, quantitative agreement with water-oil partitioning is not reached due to discrepancies between the calculated and measured free energies of solvation in hexadecane. We suggest that classical pairwise-additive models may be inadequate to properly describe the thermodynamics of solvation of apolar species, such as O2, in apolar solvents. acceptedVersion

There is evidence that lipids can be allosteric regulators of membrane protein structure and activation. However, there are no data showing how exactly the regulation emerges from specific lipid-protein interactions. Here we show in atomistic detail how the human β2-adrenergic receptor (β2AR) – a prototypical G protein-coupled receptor – is modulated by cholesterol in an allosteric fashion. Extensive atomistic simulations show that cholesterol regulates β2AR by limiting its conformational variability. The mechanism of action is based on the binding of cholesterol at specific high-affinity sites located near the transmembrane helices 5–7 of the receptor. The alternative mechanism, where the β2AR conformation would be modulated by membrane-mediated interactions, plays only a minor role. Cholesterol analogues also bind to cholesterol binding sites and impede the structural flexibility of β2AR, however cholesterol generates the strongest effect. The results highlight the capacity of lipids to regulate the conformation of membrane receptors through specific interactions., eLife, 5, ISSN:2050-084X

Extracellular and cytosolic leaflets in cellular membranes are distinctly different in lipid composition, yet they contribute together to signaling across the membranes. Here we consider a mechanism based on long-chain gangliosides for coupling the extracellular and cytosolic membrane leaflets together. Based on atomistic molecular dynamics simulations, we find that long-chain GM1 in the extracellular leaflet exhibits a strong tendency to protrude into the opposing bilayer leaflet. This interdigitation modulates the order in the cytosolic monolayer and thereby strengthens the interaction and coupling across a membrane. Coarse-grained simulations probing longer time scales in large membrane systems indicate that GM1 in the extracellular leaflet modulates the phase behavior in the cytosolic monolayer. While short-chain GM1 maintains phase-symmetric bilayers with a strong membrane registration effect, the situation is altered with long-chain GM1. Here, the significant interdigitation induced by long-chain GM1 modulates the behavior in the cytosolic GM1-free leaflet, weakening and slowing down the membrane registration process. The observed physical interaction mechanism provides a possible means to mediate or foster transmembrane communication associated with signal transduction. (C) 2017 Elsevier B.V. All rights reserved.

1,6-Diphenyl-1,3,5-hexatriene (DPH) is one of the most commonly used fluorescent probes to study dynamical and structural properties of lipid bilayers and cellular membranes via measuring steady-state or time-resolved fluorescence anisotropy. In this study, we present a limitation in the use of DPH to predict the order of lipid acyl chains when the lipid bilayer is doped with itraconazole (ITZ), an antifungal drug. Our steady-state fluorescence anisotropy measurements showed a significant decrease in fluorescence anisotropy of DPH embedded in the ITZ-containing membrane, suggesting a substantial increase in membrane fluidity, which indirectly indicates a decrease in the order of the hydrocarbon chains. This result or its interpretation is in disagreement with the fluorescence recovery after photobleaching measurements and molecular dynamics (MD) simulation data. The results of these experiments and calculations indicate an increase in the hydrocarbon chain order. The MD simulations of the bilayer containing both ITZ and DPH provide explanations for these observations. Apparently, in the presence of the drug, the DPH molecules are pushed deeper into the hydrophobic membrane core below the lipid double bonds, and the probe predominately adopts the orientation of the ITZ molecules that is parallel to the membrane surface, instead of orienting parallel to the lipid acyl chains. For this reason, DPH anisotropy provides information related to the less ordered central region of the membrane rather than reporting the properties of the upper segments of the lipid acyl chains. publishedVersion

In this Data in Brief article we provide a data package of GROMACS input files for atomistic molecular dynamics simulations of multicomponent, asymmetric lipid bilayers using the OPLS-AA force field. These data include 14 model bilayers composed of 8 different lipid molecules. The lipids present in these models are: cholesterol (CHOL), 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphatidylcholine (POPC), 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphatidylethanolamine (POPE), 1-stearoyl-2-oleoyl-sn-glycero-3-phosphatidyl-ethanolamine (SOPE), 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphatidylserine (POPS), 1-stearoyl-2-oleoyl-sn-glycero-3-phosphatidylserine (SOPS), N-palmitoyl-D-erythro-sphingosyl-phosphatidylcholine (SM16), and N-lignoceroyl-D-erythro-sphingosyl-phosphatidylcholine (SM24). The bilayers׳ compositions are based on lipidomic studies of PC-3 prostate cancer cells and exosomes discussed in Llorente et al. (2013) [1], showing an increase in the section of long-tail lipid species (SOPS, SOPE, and SM24) in the exosomes. Former knowledge about lipid asymmetry in cell membranes was accounted for in the models, meaning that the model of the inner leaflet is composed of a mixture of PC, PS, PE, and cholesterol, while the extracellular leaflet is composed of SM, PC and cholesterol discussed in Van Meer et al. (2008) [2]. The provided data include lipids׳ topologies, equilibrated structures of asymmetric bilayers, all force field parameters, and input files with parameters describing simulation conditions (md.mdp). The data is associated with the research article “Interdigitation of Long-Chain Sphingomyelin Induces Coupling of Membrane Leaflets in a Cholesterol Dependent Manner” (Róg et al., 2016) [3]., Graphical abstract

Atomistic molecular dynamics simulations were performed with 13 non-peptidic neurotransmitters (NTs) in three different membrane environments. The results provide compelling evidence that NTs are divided into membrane-binding and membrane-nonbinding molecules. NTs adhere to the postsynaptic membrane surface whenever the ligand-binding sites of their synaptic receptors are buried in the lipid bilayer. In contrast, NTs that have extracellular ligand-binding sites do not have a similar tendency to adhere to the membrane surface. This finding is a seemingly simple yet important addition to the paradigm of neurotransmission, essentially dividing it into membrane-independent and membrane-dependent mechanisms. Moreover, the simulations also indicate that the lipid composition especially in terms of charged lipids can affect the membrane partitioning of NTs. The revised paradigm, highlighting the importance of cell membrane and specific lipids for neurotransmission, should to be of interest to neuroscientists, drug industry and the general public alike.

Despite the vast amount of experimental and theoretical studies on the binding affinity of cations - especially the biologically relevant Na(+) and Ca(2+) - for phospholipid bilayers, there is no consensus in the literature. Here we show that by interpreting changes in the choline headgroup order parameters according to the 'molecular electrometer' concept [Seelig et al., Biochemistry, 1987, 26, 7535], one can directly compare the ion binding affinities between simulations and experiments. Our findings strongly support the view that in contrast to Ca(2+) and other multivalent ions, Na(+) and other monovalent ions (except Li(+)) do not specifically bind to phosphatidylcholine lipid bilayers at sub-molar concentrations. However, the Na(+) binding affinity was overestimated by several molecular dynamics simulation models, resulting in artificially positively charged bilayers and exaggerated structural effects in the lipid headgroups. While qualitatively correct headgroup order parameter response was observed with Ca(2+) binding in all the tested models, no model had sufficient quantitative accuracy to interpret the Ca(2+):lipid stoichiometry or the induced atomistic resolution structural changes. All scientific contributions to this open collaboration work were made publicly, using nmrlipids.blogspot.fi as the main communication platform. acceptedVersion

BackgroundIn lipid bilayers, cholesterol facilitates the formation of the liquid-ordered phase and enables the formation of laterally ordered structures such as lipid rafts. While these domains have an important role in a variety of cellular processes, the precise atomic-level mechanisms responsible for cholesterol's specific ordering and packing capability have remained unresolved.Methodology/principal findingsOur atomic-scale molecular dynamics simulations reveal that this ordering and the associated packing effects in membranes largely result from cholesterol's molecular structure, which differentiates cholesterol from other sterols. We find that cholesterol molecules prefer to be located in the second coordination shell, avoiding direct cholesterol-cholesterol contacts, and form a three-fold symmetric arrangement with proximal cholesterol molecules. At larger distances, the lateral three-fold organization is broken by thermal fluctuations. For other sterols having less structural asymmetry, the three-fold arrangement is considerably lost.Conclusions/significanceWe conclude that cholesterol molecules act collectively in lipid membranes. This is the main reason why the liquid-ordered phase only emerges for Chol concentrations well above 10 mol% where the collective self-organization of Chol molecules emerges spontaneously. The collective ordering process requires specific molecular-scale features that explain why different sterols have very different membrane ordering properties: the three-fold symmetry in the Chol-Chol organization arises from the cholesterol off-plane methyl groups allowing the identification of raft-promoting sterols from those that do not promote rafts.

We study the structure and dynamics of spherical high density lipoprotein (HDL) particles through coarse-grained multi-microsecond molecular dynamics simulations. We simulate both a lipid droplet without the apolipoprotein A-I (apoA-I) and the full HDL particle including two apoA-I molecules surrounding the lipid compartment. The present models are the first ones among computational studies where the size and lipid composition of HDL are realistic, corresponding to human serum HDL. We focus on the role of lipids in HDL structure and dynamics. Particular attention is paid to the assembly of lipids and the influence of lipid-protein interactions on HDL properties. We find that the properties of lipids depend significantly on their location in the particle (core, intermediate region, surface). Unlike the hydrophobic core, the intermediate and surface regions are characterized by prominent conformational lipid order. Yet, not only the conformations but also the dynamics of lipids are found to be distinctly different in the different regions of HDL, highlighting the importance of dynamics in considering the functionalization of HDL. The structure of the lipid droplet close to the HDL-water interface is altered by the presence of apoA-Is, with most prominent changes being observed for cholesterol and polar lipids. For cholesterol, slow trafficking between the surface layer and the regimes underneath is observed. The lipid-protein interactions are strongest for cholesterol, in particular its interaction with hydrophobic residues of apoA-I. Our results reveal that not only hydrophobicity but also conformational entropy of the molecules are the driving forces in the formation of HDL structure. The results provide the first detailed structural model for HDL and its dynamics with and without apoA-I, and indicate how the interplay and competition between entropy and detailed interactions may be used in nanoparticle and drug design through self-assembly., Author Summary Cardiovascular diseases are the primary cause of death in western countries. One of the main causes is lipid accumulation and plaque formation on arterial walls, called atherosclerosis. The risk of being exposed to this condition is reduced by high levels of high density lipoprotein (HDL). The functionality of HDL has remained elusive, and even its structure is not well understood. Through extensive coarse-grained simulations, we have clarified the structure of the lipid droplet in HDL and elucidated its interactions with the apolipoprotein A-I (apoA-I) that surrounds the droplet. We have found that the structural and dynamic properties of lipids depend significantly on their location in the particle (core, intermediate region, surface). As for apoA-I, we have observed it alter the overall structure of the lipid droplet close to the HDL-water interface, with prominent changes taking place for cholesterol and other polar lipids. The nature of lipid-protein interactions is most favorable for cholesterol. Our results reveal that not only hydrophobicity but also conformational entropy are the driving forces in the formation of HDL structure, suggesting how the interplay and competition between entropy and detailed interactions may be used in nanoparticle and drug design through self-assembly.