A molecular view on the organizational complexity of proteins in membranes

Alle levende wezens bestaan uit cellen, en elke cel is omgeven door een celmembraan. Tot voorkort dachten we dat de celmembraan niks meer was dan een vliesje dat de inhoud van de cel afscheid van alles daarbuiten, maar meer en meer ontdekken we dat dit vliesje, dat voornamelijk bestaat uit lipiden, cholesterol en eiwitten, nog vele andere functies heeft. In mijn onderzoek heb ik computer simulaties gebruikt om de wisselwerking tussen lipiden en eiwitten in de membraan te bestuderen. Door het kleine formaat van celmembranen en de hele korte tijdschaal waarin processen plaatsvinden, zijn de wetenschappelijke vraagstukken lastig met experimentele technieken te bestuderen. Met behulp van computer modellen, zoals het in Groningen ontwikkelde “Martini” model, kunnen we deze processen op bijna atomistische resolutie observeren. Het onderzoek richtte zich op twee aspecten. Ten eerste is het Martini krachtenveld verder getest en ontwikkeld, wat de mogelijkheid geeft tot nog betrouwbaardere simulaties. Het tweede aspect is het simuleren van kleine membraaneiwitten in een (lipide)membraan met behulp van het Martini krachtenveld. Deze simulaties hebben meer inzicht verschaft in hoe de eiwitten zich bewegen tussen verschillende domeinen binnen het membraan en de belangrijke rol die een specifiek lipide hierin speelt: het GM1 ganglioside lipide

Characterization of thylakoid lipid membranes from cyanobacteria and higher plants by molecular dynamics simulations

AbstractThe thylakoid membrane is mainly composed of non-common lipids, so called galactolipids. Despite the importance of these lipids for the function of the photosynthetic reaction centers, the molecular organization of these membranes is largely unexplored. Here we use multiscale molecular dynamics simulations to characterize the thylakoid membrane of both cyanobacteria and higher plants. We consider mixtures of up to five different galactolipids plus phosphatidylglycerol to represent these complex membranes. We find that the different lipids generally mix well, although nanoscale heterogeneities are observed especially in case of the plant membrane. The fluidity of the cyanobacterial membrane is markedly reduced compared to the plant membrane, even considering elevated temperatures at which thermophilic cyanobacteria are found. We also find that the plant membrane more readily undergoes a phase transformation to an inverted hexagonal phase. We furthermore characterized the conformation and dynamics of the cofactors plastoquinone and plastoquinol, revealing of the fast flip-flop rates for the non-reduced form. Together, our results provide a molecular view on the dynamical organization of the thylakoid membrane

The power of coarse graining in biomolecular simulations

Computational modeling of biological systems is challenging because of the multitude of spatial and temporal scales involved. Replacing atomistic detail with lower resolution, coarse grained (CG), beads has opened the way to simulate large-scale biomolecular processes on time scales inaccessible to all-atom models. We provide an overview of some of the more popular CG models used in biomolecular applications to date, focusing on models that retain chemical specificity. A few state-of-the-art examples of protein folding, membrane protein gating and self-assembly, DNA hybridization, and modeling of carbohydrate fibers are used to illustrate the power and diversity of current CG modeling

Phytochemicals Perturb Membranes and Promiscuously Alter Protein Function

A wide variety of phytochemicals are consumed for their perceived health benefits. Many of these phytochemicals have been found to alter numerous cell functions, but the mechanisms underlying their biological activity tend to be poorly understood. Phenolic phytochemicals are particularly promiscuous modifiers of membrane protein function, suggesting that some of their actions may be due to a common, membrane bilayer-mediated mechanism. To test whether bilayer perturbation may underlie this diversity of actions, we examined five bioactive phenols reported to have medicinal value: capsaicin from chili peppers, curcumin from turmeric, EGCG from green tea, genistein from soybeans, and resveratrol from grapes. We find that each of these widely consumed phytochemicals alters lipid bilayer properties and the function of diverse membrane proteins. Molecular dynamics simulations show that these phytochemicals modify bilayer properties by localizing to the bilayer/solution interface. Bilayer-modifying propensity was verified using a gramicidin-based assay, and indiscriminate modulation of membrane protein function was demonstrated using four proteins: membrane-anchored metalloproteases, mechanosensitive ion channels, and voltage-dependent potassium and sodium channels. Each protein exhibited similar responses to multiple phytochemicals, consistent with a common, bilayer-mediated mechanism. Our results suggest that many effects of amphiphilic phytochemicals are due to cell membrane perturbations, rather than specific protein binding

Dimerization of Amino Acid Side Chains: Lessons from the Comparison of Different Force Fields

The interactions between amino acid side chains govern protein secondary, tertiary, and quaternary structure formation. For molecular modeling approaches to be able to realistically describe these phenomena, the underlying force fields have to represent these interactions as accurately as possible. Here, we compare the side chain−side chain interactions for a number of commonly used force fields, namely the all-atom OPLS, the united-atom GROMOS, and the coarse-grain MARTINI force field. We do so by calculating the dimerization free energies between selected pairs of side chains and structural characterization of their binding modes. To mimic both polar and nonpolar environments, the simulations are performed in water, n-octanol, and decane. In general, reasonable correlations are found between all three force fields, with deviations on the order of 1 kT in aqueous solvent. In apolar solvent, however, significantly larger differences are found, especially for charged amino acid pairs between the OPLS and GROMOS force fields, and for polar interactions in the MARTINI force field in comparison to the higher resolution models. Interestingly, even in cases where the dimerization free energies are similar, the binding mode may differ substantially between the force fields. This was found to be especially the case for aromatic residues. In addition to the inter-force-field comparison, we compared the various force fields to a knowledge-based potential. The two independent approaches show good correlation in aqueous solvent with an exception of aromatic residues for which the interaction strength is lower in the knowledge-based potentials.

Molecular view on protein sorting into liquid-ordered membrane domains mediated by gangliosides and lipid anchors

We present results from coarse grain molecular dynamics simulations of mixed model membranes consisting of saturated and unsaturated lipids together with cholesterol, in which lipid-anchored membrane proteins are embedded. The membrane proteins studied are the peripherally bound H-Ras, N-Ras, and Hedgehog, and the transmembrane peptides WALP and LAT. We provide a molecular view on how the presence and nature of these lipid anchors affects partitioning of the proteins between liquid-ordered and liquid-disordered domains. In addition, we probed the role of the ganglioside lipid GM1 on the protein sorting, showing formation of GM1-protein nano-domains that act as shuttles between the differently ordered membrane regions.

Simulation of polyethylene glycol and calcium-mediated membrane fusion

We report on the mechanism of membrane fusion mediated by polyethylene glycol (PEG) and Ca2+ by means of a coarse-grained molecular dynamics simulation approach. Our data provide a detailed view on the role of cations and polymer in modulating the interaction between negatively charged apposed membranes. The PEG chains cause a reduction of the inter-lamellar distance and cause an increase in concentration of divalent cations. When thermally driven fluctuations bring the membranes at close contact, a switch from cis to trans Ca2+-lipid complexes stabilizes a focal contact acting as a nucleation site for further expansion of the adhesion region. Flipping of lipid tails induces subsequent stalk formation. Together, our results provide a molecular explanation for the synergistic effect of Ca2+ and PEG on membrane fusion. (c) 2014 AIP Publishing LLC

Gaussian curvature elasticity determined from global shape transformations and local stress distributions: a comparative study using the MARTINI model

We calculate the Gaussian curvature modulus ¯κ of a systematically coarse-grained (CG) one-component lipid membrane by applying a recently proposed method to the MARTINI representation of 1,2-dimyristoyl-sn-glycero-3-phosphocholine (DMPC). We find the value ¯κ/κ = −1.04 ± 0.03 for the elastic ratio between the Gaussian and the mean curvature modulus and deduce ¯κm/κm ≈ −0.98 ± 0.09 for the monolayer elastic ratio, where the latter is based on plausible assumptions for the distance z0 of the monolayer neutral surface from the bilayer midplane and the spontaneous lipid curvature K0m. By also analyzing the lateral stress profile σ0(z) of our system, two other lipid types and pertinent data from the literature, we show that determining K0m and ¯κ through the first and second moment of σ0(z) gives rise to physically implausible values for these observables. This discrepancy, which we previously observed for a much simpler CG model, suggests that the moment conditions derived from simple continuum assumptions miss the effect of physically important correlations in the lipid bilayer.