Waldemar Kulig

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.

Cholesterol oxidation products and their biological importance.

The main biological cause of oxysterols is the oxidation of cholesterol. They differ from cholesterol by the presence of additional polar groups that are typically hydroxyl, keto, hydroperoxy, epoxy, or carboxyl moieties. Under typical conditions, oxysterol concentration is maintained at a very low and precisely regulated level, with an excess of cholesterol. Like cholesterol, many oxysterols are hydrophobic and hence confined to cell membranes. However, small chemical differences between the sterols can significantly affect how they interact with other membrane components, and this in turn can have a substantial effect on membrane properties. In this spirit, this review describes the biological importance and the roles of oxysterols in the human body. We focus primarily on the effect of oxysterols on lipid membranes, but we also consider other issues such as enzymatic and nonenzymatic synthesis processes of oxysterols as well as pathological conditions induced by oxysterols.

Cis and trans unsaturated phosphatidylcholine bilayers: A molecular dynamics simulation study.

Trans unsaturated lipids are uncommon in nature. In the human diet, they occur as natural products of ruminal bacteria or from industrial food processing like hydrogenation of vegetable oils. Consumption of trans unsaturated lipids has been shown to have a negative influence on human health; in particular, the risk of cardiovascular disease is higher when the amount of trans unsaturated lipids in the diet is elevated. In this study, we first performed quantum mechanical calculations to specifically and accurately parameterize cis and trans mono-unsaturated lipids and subsequently validated the newly derived parameter set. Then, we carried out molecular dynamics (MD) simulations of lipid bilayers composed of cis or trans unsaturated lipids with and without cholesterol. Our results show that trans mono-unsaturated chains are more flexible than cis mono-unsaturated chains due to lower barriers for rotation around the single bonds next to the trans double bond than those next to the cis double bond. In effect, interactions between cholesterol and trans unsaturated chains are stronger than cis unsaturated chains, which results in a higher ordering effect of cholesterol in trans unsaturated bilayers.

A sodium atom in a large water cluster: Electron delocalization and infrared spectra

Ab initio molecular dynamics simulations modeling low-energy collisions of a sodium atom with a cluster with more than 30 water molecules are presented. We follow the dynamics of the atom-cluster interaction and the delocalization of the valence electron of sodium together with the changes in the electron binding energy. This electron tends to be shared by the nascent sodium cation and the water cluster. IR spectra of the sodium-water cluster are both computationally and experimentally obtained, with a good agreement between the two approaches.

Cholesterol under oxidative stress-How lipid membranes sense oxidation as cholesterol is being replaced by oxysterols.

The behavior of oxysterols in phospholipid membranes and their effects on membrane properties were investigated by means of dynamic light scattering, fluorescence spectroscopy, NMR, and extensive atomistic simulations. Two families of oxysterols were scrutinized-tail-oxidized sterols, which are mostly produced by enzymatic processes, and ring-oxidized sterols, formed mostly via reactions with free radicals. The former family of sterols was found to behave similar to cholesterol in terms of molecular orientation, roughly parallel to the bilayer normal, leading to increasing membrane stiffness and suppression of its membrane permeability. In contrast, ring-oxidized sterols behave quantitatively differently from cholesterol. They acquire tilted orientations and therefore disrupt the bilayer structure with potential implications for signaling and other biochemical processes in the membranes.

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

Topologies, structures and parameter files for lipid simulations in GROMACS with the OPLS-aa force field: DPPC, POPC, DOPC, PEPC, and cholesterol

In this data article we provide topologies and force field parameters files for molecular dynamics simulations of lipids in the OPLS-aa force field using the GROMACS package. This is the first systematic parameterization of lipid molecules in this force field. Topologies are provided for four phosphatidylcholines: saturated DPPC, mono-cis unsaturated POPC and DOPC, and mono-trans unsaturated PEPC. Parameterization of the phosphatidylcholines was achieved in two steps: first, we supplemented the OPLS force field parameters for DPPC with new parameters for torsion angles and van der Waals parameters for the carbon and hydrogen atoms in the acyl chains, as well as new partial atomic charges and parameters for torsion angles in the phosphatidylcholine and glycerol moieties [1]. Next, we derived parameters for the cis and trans double bonds and the neighboring them single bonds [2]. Additionally, we provide GROMACS input files with parameters describing simulation conditions (md.mdp), which are strongly recommended to be used with these lipids models. The data are associated with the research article “Cis and trans unsaturated phosphatidylcholine bilayers: a molecular dynamics simulation study” [2] and provided as supporting materials.

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.

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.

Steric and electronic effects in the host-guest hydrogen bonding in clathrate hydrates.

Clathrate hydrates with polar guest molecules (dimethyl ether, ethylene oxide, trimethylene oxide, tetrahydrofuran, and tetrahydropyran) were studied by means of the density functional theory. A model of a large cage of structure-I clathrate was employed. Optimal configurations of encaged guests were investigated with a focus on the host-guest hydrogen bond formation. Weak hydrogen bonds were found to be formed by each guest, while for THP a strong hydrogen bond and formation of L-defect was also observed. This is in accord with previous computational and experimental studies. Steric factors were shown to play a key role for the strength of the hydrogen bond formed. Interestingly, the host-guest binding is influenced not only by the size of a guest molecule but also by its shape. This work demonstrates that both electronic and steric properties of a polar guest should be considered for a full description of clathrate systems.