Molecular Dynamics of Photosystem II Embedded in the Thylakoid Membrane.

Photosystem II (PSII) is one of the key protein complexes in photosynthesis. We introduce a coarse grained model of PSII and present the analysis of 60 µs molecular dynamics simulations of PSII in both monomeric and dimeric form, embedded in a thylakoid membrane model that reflects its native lipid composition. We describe in detail the setup of the protein complex and the many natural cofactors and characterize their mobility. Overall we find that the protein subunits and cofactors are more flexible toward the periphery of the complex as well as near the PLQ exchange cavity and at the dimer interface. Of all cofactors, ß-carotenes show the highest mobility. Some of the ß-carotenes diffuse in and out of the protein complex via the thylakoid membrane. In contrast with the PSII dimer, the monomeric form adopts a tilted conformation in the membrane, with strong interactions between the soluble PsbO subunit and the glycolipid headgroups. Interestingly, the tilted conformation causes buckling of the membrane. Together, our results provide an unprecedented view of PSII dynamics on a microsecond time scale. Our data may be used as basis for the interpretation of experimental data as well as for theoretical models describing exciton energy transfer.

Imidazolium-Based Lipid Analogues and Their Interaction with Phosphatidylcholine Membranes.

4,5-Dialkylated imidazolium lipid salts are a new class of lipid analogues showing distinct biological activities. The potential effects of the imidazolium lipids on artificial lipid membranes and the corresponding membrane interactions was analyzed. Therefore, 1,2-dipalmitoyl-sn-glycero-3-phosphocholine (DPPC) was employed to create an established lipid monolayer model and a bilayer membrane. Mixed monolayers of DPPC and 4,5-dialkylimidazolium lipids differing by their alkyl chain length (C7, C11, and C15) were characterized by surface pressure-area (π-A) isotherms using a Wilhelmy film balance in combination with epifluorescence microscopy. Monolayer hysteresis for binary mixtures was examined by recording triplicate consecutive compression-expansion cycles. The lipid miscibility and membrane stability of DPPC/imidazolium lipids were subsequently evaluated by the excess mean molecular area (ΔAex) and the excess Gibbs free energy (ΔGex) of mixing. Furthermore, the thermotropic behavior of mixed liposomes of DPPC/imidazolium lipids was investigated by differential scanning calorimetry (DSC). The C15-imidazolium lipid (C15-IMe·HI) forms a thermodynamically favored and kinetically reversible Langmuir monolayer with DPPC and exhibits a rigidification effect on both DPPC monolayer and bilayer structures at low molar fractions (X ≤ 0.3). However, the incorporation of the C11-imidazolium lipid (C11-IMe·HI) causes the formation of an unstable and irreversible Langmuir-Gibbs monolayer with DPPC and disordered DPPC liposomes. The C7-imidazolium lipid (C7-IMe·HI) displays negligible membrane activity. To better understand these results on a molecular level, all-atom molecular dynamics (MD) simulations were performed. The simulations yield two opposing molecular mechanisms governing the different behavior of the three imidazolium lipids: a lateral ordering effect and a free volume/stretching effect. Overall, our study provides the first evidence that the membrane interaction of the C15 and C11 derivatives modulates the structural organization of lipid membranes. On the contrary, for the C7 derivative its membrane activity is too low to contribute to its earlier reported potent cytotoxicity.

P3HT:DiPBI bulk heterojunction solar cells: morphology and electronic structure probed by multiscale simulation and UV/vis spectroscopy.

Coarse grained molecular dynamics simulations are performed for a mixture of poly(3-hexylthiophene) (P3HT) and diperylene bisimide (DiPBI). The effect of different annealing and cooling protocols on the morphology is investigated and the resulting domain structures are analyzed. In particular, π-stacked clusters of DiPBI molecules are observed whose size decreases with increasing temperature. Domain structure and diffusivity data suggest that the DiPBI subsystem undergoes an order → disorder phase transition between 700 and 900 K. Electronic structure calculations based on density functional theory are carried out after backmapping the coarse grained model onto an atomistic force field representation built upon first principles. UV/vis absorption spectra of the P3HT:DiPBI mixture are computed using time-dependent density functional linear response theory and recorded experimentally for a spin-coated thin film. It is demonstrated that the absorption spectrum depends sensitively on the details of the amorphous structure, thus providing valuable insight into the morphology. In particular, the results show that the tempering procedure has a significant influence on the material's electronic properties. This knowledge may help to develop effective processing routines to enhance the performance of bulk heterojunction solar cells.

Dry Martini, a coarse-grained force field for lipid membrane simulations with implicit solvent.

Coarse-grained (CG) models allow simulation of larger systems for longer times by decreasing the number of degrees of freedom compared with all-atom models. Here we introduce an implicit-solvent version of the popular CG Martini model, nicknamed "Dry" Martini. To account for the omitted solvent degrees of freedom, the nonbonded interaction matrix underlying the Martini force field was reparametrized. The Dry Martini force field reproduces relatively well a variety of lipid membrane properties such as area per lipid, bilayer thickness, bending modulus, and coexistence of liquid-ordered and disordered domains. Furthermore, we show that the new model can be applied to study membrane fusion and tether formation, with results similar to those of the standard Martini model. Membrane proteins can also be included, but less quantitative results are obtained. The absence of water in Dry Martini leads to a significant speedup for large systems, opening the way to the study of complex multicomponent membranes containing millions of lipids.

Structure of P3HT crystals, thin films, and solutions by UV/Vis spectral analysis.

Optical absorption spectra of poly(3-hexylthiophene) (P3HT) are calculated in solution, spin-coated thin films, and the bulk crystal using a multiscale simulation approach. The structure of the amorphous thin film is obtained from coarse grained molecular dynamics (MD) simulations and subsequent back-mapping onto an atomistic force field representation. The absorption spectra are computed using TDDFT by statistically averaging over an ensemble of molecules taken from the MD simulations. Experimental UV/Vis spectra of spin-coated thin films and solutions are recorded with varying ratios of 'good' versus 'poor' solvent. The theoretical approach is able to faithfully predict the spectral position in the various phases and offers fundamental insight into the cause of any spectral shifts. The position of the main absorption peak is found to be chiefly determined by the level of torsion between the thiophene rings inside each molecule, while intermolecular effects are less important. Hence, optical absorption spectra hold valuable clues about the microscopic structure of disordered P3HT phases.

Atomistic and Coarse Grain Topologies for the Cofactors Associated with the Photosystem II Core Complex.

Electron transfers within and between protein complexes are core processes of the electron transport chains occurring in thylakoid (chloroplast), mitochondrial, and bacterial membranes. These electron transfers involve a number of cofactors. Here we describe the derivation of molecular mechanics parameters for the cofactors associated with the function of the photosystem II core complex: plastoquinone, plastoquinol, heme b, chlorophyll A, pheophytin, and ß-carotene. Parameters were also obtained for ubiquinol and ubiquinone, related cofactors involved in the respiratory chain. Parameters were derived at both atomistic and coarse grain (CG) resolutions, compatible with the building blocks of the GROMOS united-atom and Martini CG force fields, respectively. Structural and thermodynamic properties of the cofactors were compared to experimental values when available. The topologies were further tested in molecular dynamics simulations of the cofactors in their physiological environment, e.g., either in a lipid membrane environment or in complex with the heme binding protein bacterioferritin.

Optimized co-solute paramagnetic relaxation enhancement for the rapid NMR analysis of a highly fibrillogenic peptide.

Co-solute paramagnetic relaxation enhancement (PRE) is an attractive way to speed up data acquisition in NMR spectroscopy by shortening the T 1 relaxation time of the nucleus of interest and thus the necessary recycle delay. Here, we present the rationale to utilize high-spin iron(III) as the optimal transition metal for this purpose and characterize the properties of its neutral chelate form Fe(DO3A) as a suitable PRE agent. Fe(DO3A) effectively reduces the T 1 values across the entire sequence of the intrinsically disordered protein α-synuclein with negligible impact on line width. The agent is better suited than currently used alternatives, shows no specific interaction with the polypeptide chain and, due to its high relaxivity, is effective at low concentrations and in 'proton-less' NMR experiments. By using Fe(DO3A) we were able to complete the backbone resonance assignment of a highly fibrillogenic peptide from α1-antitrypsin by acquiring the necessary suite of multidimensional NMR datasets in 3 h.

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

The 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.

Simulation of polyethylene glycol and calcium-mediated membrane fusion.

We report on the mechanism of membrane fusion mediated by polyethylene glycol (PEG) and Ca(2+) 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 Ca(2+)-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 Ca(2+) and PEG on membrane fusion.

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.

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 kappa of a systematically coarse-grained (CG) one-component lipid membrane by applying the method recently proposed by Hu et al. [Biophys. J., 2012, 102, 1403] to the MARTINI representation of 1,2-dimyristoyl-sn-glycero-3-phosphocholine (DMPC). We find the value kappa/kappa = -1.04 +/- 0.03 for the elastic ratio between the Gaussian and the mean curvature modulus and deduce kappa(m)/kappa(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 K(0m). By also analyzing the lateral stress profile sigma0(z) of our system, two other lipid types and pertinent data from the literature, we show that determining K(0m) and kappa through the first and second moment of sigma0(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.

Improved Parameters for the Martini Coarse-Grained Protein Force Field.

The Martini coarse-grained force field has been successfully used for simulating a wide range of (bio)molecular systems. Recent progress in our ability to test the model against fully atomistic force fields, however, has revealed some shortcomings. Most notable, phenylalanine and proline were too hydrophobic, and dimers formed by polar residues in apolar solvents did not bind strongly enough. Here, we reparametrize these residues either through reassignment of particle types or by introducing embedded charges. The new parameters are tested with respect to partitioning across a lipid bilayer, membrane binding of Wimley-White peptides, and dimerization free energy in solvents of different polarity. In addition, we improve some of the bonded terms in the Martini protein force field that lead to a more realistic length of α-helices and to improved numerical stability for polyalanine and glycine repeats. The new parameter set is denoted Martini version 2.2.

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.

Determining equilibrium constants for dimerization reactions from molecular dynamics simulations.

With today's available computer power, free energy calculations from equilibrium molecular dynamics simulations "via counting" become feasible for an increasing number of reactions. An example is the dimerization reaction of transmembrane alpha-helices. If an extended simulation of the two helices covers sufficiently many dimerization and dissociation events, their binding free energy is readily derived from the fraction of time during which the two helices are observed in dimeric form. Exactly how the correct value for the free energy is to be calculated, however, is unclear, and indeed several different and contradictory approaches have been used. In particular, results obtained via Boltzmann statistics differ from those determined via the law of mass action. Here, we develop a theory that resolves this discrepancy. We show that for simulation systems containing two molecules, the dimerization free energy is given by a formula of the form ΔG ∝ ln(P(1) /P(0) ). Our theory is also applicable to high concentrations that typically have to be used in molecular dynamics simulations to keep the simulation system small, where the textbook dilute approximations fail. It also covers simulations with an arbitrary number of monomers and dimers and provides rigorous error estimates. Comparison with test simulations of a simple Lennard Jones system with various particle numbers as well as with reference free energy values obtained from radial distribution functions show full agreement for both binding free energies and dimerization statistics.

Lipid packing drives the segregation of transmembrane helices into disordered lipid domains in model membranes.

Cell membranes are comprised of multicomponent lipid and protein mixtures that exhibit a complex partitioning behavior. Regions of structural and compositional heterogeneity play a major role in the sorting and self-assembly of proteins, and their clustering into higher-order oligomers. Here, we use computer simulations and optical microscopy to study the sorting of transmembrane helices into the liquid-disordered domains of phase-separated model membranes, irrespective of peptide-lipid hydrophobic mismatch. Free energy calculations show that the enthalpic contribution due to the packing of the lipids drives the lateral sorting of the helices. Hydrophobic mismatch regulates the clustering into either small dynamic or large static aggregates. These results reveal important molecular driving forces for the lateral organization and self-assembly of transmembrane helices in heterogeneous model membranes, with implications for the formation of functional protein complexes in real cells.

Cholesterol in bilayers with PUFA chains: doping with DMPC or POPC results in sterol reorientation and membrane-domain formation.

Using neutron diffraction Harroun et al. [(2006) Biochemistry 45, 1227-1233; (2008) Biochemistry 47, 7090-7096] carried out studies that unequivocally demonstrated cholesterol preferentially sequestering in the middle of bilayers (i.e., flat orientation) made of lipids with polyunsaturated fatty acids (PUFA), in contrast to its "usual" position where its hydroxyl group locates near the lipid/water interface (i.e., upright orientation). Here we clearly show, using neutron diffraction, cholesterol's orientational preference in different lipid bilayers. For example, although it requires 50 mol % POPC (16:0-18:1 PC) in DAPC (di20:4 PC) bilayers to cause cholesterol to revert to its upright orientation, only 5 mol % DMPC (di14:0 PC) is needed to achieve the same effect. This result demonstrates not only cholesterol's affinity for saturated hydrocarbon chains, but also its aversion for PUFAs. Molecular dynamics (MD) simulations performed on similar systems show that in high PUFA content bilayers cholesterol is simultaneously capable of assuming different orientations within a bilayer. Although this result is known from previous MD studies by Marrink et al. [(2008) J. Am. Chem. Soc. 130, 10-11], it has yet to be confirmed experimentally. Importantly, MD simulations predict the formation of DMPC-rich domains, data corroborated by experiment (i.e., 10 mol % DMPC-doped DAPC bilayers), where cholesterol preferentially locates in its upright orientation, while in DMPC-depleted domains cholesterol is found mostly in the bilayer center (i.e., flat orientation). These results lend credence to DMPC's aversion for PUFAs, supporting the notion that domain formation is primarily driven by lipids.