Matti Javanainen

Anomalous diffusion of phospholipids and cholesterols in a lipid bilayer and its origins.

Combining extensive molecular dynamics simulations of lipid bilayer systems of varying chemical compositions with single-trajectory analyses, we systematically elucidate the stochastic nature of the lipid motion. We observe subdiffusion over more than 4 orders of magnitude in time, clearly stretching into the submicrosecond domain. The lipid motion depends on the lipid chemistry, the lipid phase, and especially the presence of cholesterol. We demonstrate that fractional Langevin equation motion universally describes the lipid motion in all phases, including the gel phase, and in the presence of cholesterol. The results underline the relevance of anomalous diffusion in lipid bilayers and the strong effects of the membrane composition.

Toward Atomistic Resolution Structure of Phosphatidylcholine Headgroup and Glycerol Backbone at Different Ambient Conditions

Phospholipids are essential building blocks of biological membranes. Despite a vast amount of very accurate experimental data, the atomistic resolution structures sampled by the glycerol backbone and choline headgroup in phoshatidylcholine bilayers are not known. Atomistic resolution molecular dynamics simulations have the potential to resolve the structures, and to give an arrestingly intuitive interpretation of the experimental data, but only if the simulations reproduce the data within experimental accuracy. In the present work, we simulated phosphatidylcholine (PC) lipid bilayers with 13 different atomistic models, and compared simulations with NMR experiments in terms of the highly structurally sensitive C–H bond vector order parameters. Focusing on the glycerol backbone and choline headgroups, we showed that the order parameter comparison can be used to judge the atomistic resolution structural accuracy of the models. Accurate models, in turn, allow molecular dynamics simulations to be used as an interpretation tool that translates these NMR data into a dynamic three-dimensional representation of biomolecules in biologically relevant conditions. In addition to lipid bilayers in fully hydrated conditions, we reviewed previous experimental data for dehydrated bilayers and cholesterol-containing bilayers, and interpreted them with simulations. Although none of the existing models reached experimental accuracy, by critically comparing them we were able to distill relevant chemical information: (1) increase of choline order parameters indicates the P–N vector tilting more parallel to the membrane, and (2) cholesterol induces only minor changes to the PC (glycerol backbone) structure. This work has been done as a fully open collaboration, using nmrlipids.blogspot.fi as a communication platform; all the scientific contributions were made publicly on this blog. During the open research process, the repository holding our simulation trajectories and files (https://zenodo.org/collection/user-nmrlipids) has become the most extensive publicly available collection of molecular dynamics simulation trajectories of lipid bilayers.

Membrane omega-3 fatty acids modulate the oligomerisation kinetics of adenosine A2A and dopamine D2 receptors

Membrane levels of docosahexaenoic acid (DHA), an essential omega-3 polyunsaturated fatty acid (ω-3 PUFA), are decreased in common neuropsychiatric disorders. DHA modulates key cell membrane properties like fluidity, thereby affecting the behaviour of transmembrane proteins like G protein-coupled receptors (GPCRs). These receptors, which have special relevance for major neuropsychiatric disorders have recently been shown to form dimers or higher order oligomers, and evidence suggests that DHA levels affect GPCR function by modulating oligomerisation. In this study, we assessed the effect of membrane DHA content on the formation of a class of protein complexes with particular relevance for brain disease: adenosine A2A and dopamine D2 receptor oligomers. Using extensive multiscale computer modelling, we find a marked propensity of DHA for interaction with both A2A and D2 receptors, which leads to an increased rate of receptor oligomerisation. Bioluminescence resonance energy transfer (BRET) experiments performed on living cells suggest that this DHA effect on the oligomerisation of A2A and D2 receptors is purely kinetic. This work reveals for the first time that membrane ω-3 PUFAs play a key role in GPCR oligomerisation kinetics, which may have important implications for neuropsychiatric conditions like schizophrenia or Parkinson’s disease.

Mechanism for translocation of fluoroquinolones across lipid membranes

Classical atom-scale molecular dynamics simulations, constrained free energy calculations, and quantum mechanical (QM) calculations are employed to study the diffusive translocation of ciprofloxacin (CPFX) across lipid membranes. CPFX is considered here as a representative of the fluoroquinolone antibiotics class. Neutral and zwitterionic CPFX coexist at physiological pH, with the latter being predominant. Simulations reveal that only the neutral form permeates the bilayer, and it does so through a novel mechanism that involves dissolution of concerted stacks of zwitterionic ciprofloxacins. Subsequent QM analysis of the observed molecular stacking shows the important role of partial charge neutralization in the stacks, highlighting how the zwitterionic form of the drug is neutralized for translocation. The findings propose a translocation mechanism in which zwitterionic CPFX molecules approach the membrane in stacks, but they diffuse through the membrane as neutral CPFX monomers due to intermolecular transfer of protons favored by partial solvation loss. The mechanism is expected to be of importance in the permeation and translocation of a variety of ampholitic drugs with stacking tendencies.

Enzymatic oxidation of cholesterol: properties and functional effects of cholestenone in cell membranes.

Bacterial cholesterol oxidase is commonly used as an experimental tool to reduce cellular cholesterol content. That the treatment also generates the poorly degradable metabolite 4-cholesten-3-one (cholestenone) has received less attention. Here, we investigated the membrane partitioning of cholestenone using simulations and cell biological experiments and assessed the functional effects of cholestenone in human cells. Atomistic simulations predicted that cholestenone reduces membrane order, undergoes faster flip-flop and desorbs more readily from membranes than cholesterol. In primary human fibroblasts, cholestenone was released from membranes to physiological extracellular acceptors more avidly than cholesterol, but without acceptors it remained in cells over a day. To address the functional effects of cholestenone, we studied fibroblast migration during wound healing. When cells were either cholesterol oxidase treated or part of cellular cholesterol was exchanged for cholestenone with cyclodextrin, cell migration during 22 h was markedly inhibited. Instead, when a similar fraction of cholesterol was removed using cyclodextrin, cells replenished their cholesterol content in 3 h and migrated similarly to control cells. Thus, cholesterol oxidation produces long-term functional effects in cells and these are in part due to the generated membrane active cholestenone.

Mechanism of allosteric regulation of β2-adrenergic receptor by cholesterol

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. DOI: http://dx.doi.org/10.7554/eLife.18432.001

How well does cholesteryl hemisuccinate mimic cholesterol in saturated phospholipid bilayers

Cholesteryl hemisuccinate is a detergent that is often used to replace cholesterol in crystallization of membrane proteins. Here we employ atomistic molecular dynamics simulations to characterize how well the properties of cholesteryl hemisuccinate actually match those of cholesterol in saturated protein-free lipid membranes. We show that the protonated form of cholesteryl hemisuccinate mimics many of the membrane properties of cholesterol quite well, while the deprotonated form of cholesteryl hemisuccinate is less convincing in this respect. Based on the results, we suggest that cholesteryl hemisuccinate in its protonated form is a quite faithful mimic of cholesterol for membrane protein crystallization, if specific cholesterol-protein interactions (not investigated here) are not playing a crucial role.

Universal Method for Embedding Proteins into Complex Lipid Bilayers for Molecular Dynamics Simulations

The growing interest toward membrane protein simulations calls for a universal and efficient protocol for embedding multiple proteins simultaneously into asymmetric many-component lipid membranes. To this end, here, we present a new and simple technique based on pushing proteins into a lipid membrane from its side by applying a high lateral pressure on the system. This approach is compatible with most (if not all) simulation software packages, and it is independent of external program codes. Most importantly, our protocol does not alter the lipid composition or the transmembrane lipid distribution of the host bilayer membrane. It can be employed with both atomistic and coarse-grained models, and it allows multiple proteins to be embedded into a membrane at the same time. It is shown that the protein structure is unaffected by the pressure applied to the system during the procedure, and the simulation resources required for protein insertion are shown to be modest, regarding both atomistic and coarse-grained simulation models.

Experimental determination and computational interpretation of biophysical properties of lipid bilayers enriched by cholesteryl hemisuccinate

Cholesteryl hemisuccinate (CHS) is one of the cholesterol-mimicking detergents not observed in nature. It is, however, widely used in protein crystallography, in biochemical studies of proteins, and in pharmacology. Here, we performed an extensive experimental and theoretical study on the behavior of CHS in lipid membranes rich in unsaturated phospholipids. We found that the deprotonated form of CHS (that is the predominant form under physiological conditions) does not mimic cholesterol very well. The protonated form of CHS does better in this regard, but also its ability to mimic the physical effects of cholesterol on lipid membranes is limited. Overall, although ordering and condensing effects characteristic to cholesterol are present in systems containing any form of CHS, their strength is appreciably weaker compared to cholesterol. Based on the considerable amount of experimental and atomistic simulation data, we conclude that these differences originate from the fact that the ester group of CHS does not anchor it in an optimal position at the water-membrane interface. The implications of these findings for considerations of protein-cholesterol interactions are briefly discussed.

Excessive aggregation of membrane proteins in the Martini model

The coarse-grained Martini model is employed extensively to study membrane protein oligomerization. While this approach is exceptionally promising given its computational efficiency, it is alarming that a significant fraction of these studies demonstrate unrealistic protein clusters, whose formation is essentially an irreversible process. This suggests that the protein–protein interactions are exaggerated in the Martini model. If this held true, then it would limit the applicability of Martini to study multi-protein complexes, as the rapidly clustering proteins would not be able to properly sample the correct dimerization conformations. In this work we first demonstrate the excessive protein aggregation by comparing the dimerization free energies of helical transmembrane peptides obtained with the Martini model to those determined from FRET experiments. Second, we show that the predictions provided by the Martini model for the structures of transmembrane domain dimers are in poor agreement with the corresponding structures resolved using NMR. Next, we demonstrate that the first issue can be overcome by slightly scaling down the Martini protein–protein interactions in a manner, which does not interfere with the other Martini interaction parameters. By preventing excessive, irreversible, and non-selective aggregation of membrane proteins, this approach renders the consideration of lateral dynamics and protein–lipid interactions in crowded membranes by the Martini model more realistic. However, this adjusted model does not lead to an improvement in the predicted dimer structures. This implicates that the poor agreement between the Martini model and NMR structures cannot be cured by simply uniformly reducing the interactions between all protein beads. Instead, a careful amino-acid specific adjustment of the protein–protein interactions is likely required.