Coarse-Grained Parameterization of Nucleotide Cofactors and Metabolites: Protonation Constants, Partition Coefficients, and Model Topologies.

Nucleotides are structural units relevant not only in nucleic acids but also as substrates or cofactors in key biochemical reactions. The size- and timescales of such nucleotide-protein interactions fall well within the scope of coarse-grained molecular dynamics, which holds promise of important mechanistic insight. However, the lack of specific parameters has prevented accurate coarse-grained simulations of protein interactions with most nucleotide compounds. In this work, we comprehensively develop coarse-grained parameters for key metabolites/cofactors (FAD, FMN, riboflavin, NAD, NADP, ATP, ADP, AMP, and thiamine pyrophosphate) in different oxidation and protonation states as well as for smaller molecules derived from them (among others, nicotinamide, adenosine, adenine, ribose, thiamine, and lumiflavin), summing up a total of 79 different molecules. In line with the Martini parameterization methodology, parameters were tuned to reproduce octanol-water partition coefficients. Given the lack of existing data, we set out to experimentally determine these partition coefficients, developing two methodological approaches, based on 31P-NMR and fluorescence spectroscopy, specifically tailored to the strong hydrophilicity of most of the parameterized compounds. To distinguish the partition of each relevant protonation species, we further potentiometrically characterized the protonation constants of key molecules. This work successfully builds a comprehensive and relevant set of computational models that will boost the biochemical application of coarse-grained simulations. It does so based on the measurement of partition and acid-base physicochemical data that, in turn, covers important gaps in nucleotide characterization.

Computational 'microscopy' of cellular membranes.

Computational 'microscopy' refers to the use of computational resources to simulate the dynamics of a molecular system. Tuned to cell membranes, this computational 'microscopy' technique is able to capture the interplay between lipids and proteins at a spatio-temporal resolution that is unmatched by other methods. Recent advances allow us to zoom out from individual atoms and molecules to supramolecular complexes and subcellular compartments that contain tens of millions of particles, and to capture the complexity of the crowded environment of real cell membranes. This Commentary gives an overview of the main concepts of computational 'microscopy' and describes the state-of-the-art methods used to model cell membrane processes. We illustrate the power of computational modelling approaches by providing a few in-depth examples of large-scale simulations that move up from molecular descriptions into the subcellular arena. We end with an outlook towards modelling a complete cell in silico.

Identification of Two New Cholesterol Interaction Sites on the A2A Adenosine Receptor.

By mole, cholesterol is the most abundant component of animal cell plasma membranes. Many membrane proteins have been shown to be functionally dependent on cholesterol, several of which have also been shown to bind cholesterol at well-defined locations on their membrane-facing surface. In this work, a combination of coarse-grained "Martini" and all-atom simulations are used to identify two, to our knowledge, new cholesterol-binding sites on the A2A adenosine receptor, a G-protein-coupled receptor that is a target for the treatment of Parkinson's disease. One of the sites is also observed to bind cholesterol in several recent, high-resolution crystal structures of the protein, and in the simulations, interacts with cholesterol only when bound to the inverse agonist ZM241385. Cataloguing cholesterol-binding sites is a vital step in the effort to understand cholesterol-dependent function of membrane proteins. Given that cholesterol content in plasma membranes varies with cell type and on administration of widely prescribed pharmaceuticals, such as statins, understanding cholesterol-dependent function is an important step toward exploiting membrane compositional variation for therapeutic purposes.

High-Throughput Simulations Reveal Membrane-Mediated Effects of Alcohols on MscL Gating.

The mechanosensitive channels of large conductance (MscL) are bacterial membrane proteins that serve as last resort emergency release valves in case of severe osmotic downshock. Sensing bilayer tension, MscL channels are sensitive to changes in the bilayer environment and are, therefore, an ideal test case for exploring membrane protein coupling. Here, we use high-throughput coarse-grained molecular dynamics simulations to characterize MscL gating kinetics in different bilayer environments under the influence of alcohols. We performed over five hundred simulations to obtain sufficient statistics to reveal the subtle effects of changes in the membrane environment on MscL gating. MscL opening times were found to increase with the addition of the straight-chain alcohols ethanol, octanol, and to some extent dodecanol but not with hexadecanol. Increasing concentration of octanol increased the impeding effect, but only up to 10-20 mol %. Our in silico predictions were experimentally confirmed using reconstituted MscL in a liposomal fluorescent efflux assay. Our combined data reveal that the effect of alcohols on MscL gating arises not through specific binding sites but through a combination of the alcohol-induced changes to a number of bilayer properties and their alteration of the MscL-bilayer interface. Our work provides a key example of how extensive molecular simulations can be used to predict the functional modification of membrane proteins by subtle changes in their bilayer environment.

The activation mode of the mechanosensitive ion channel, MscL, by lysophosphatidylcholine differs from tension-induced gating.

One of the best-studied mechanosensitive channels is the mechanosensitive channel of large conductance (MscL). MscL senses tension in the membrane evoked by an osmotic down shock and directly couples it to large conformational changes leading to the opening of the channel. Spectroscopic techniques offer unique possibilities to monitor these conformational changes if it were possible to generate tension in the lipid bilayer, the native environment of MscL, during the measurements. To this end, asymmetric insertion of l-α-lysophosphatidylcholine (LPC) into the lipid bilayer has been effective; however, how LPC activates MscL is not fully understood. Here, the effects of LPC on tension-sensitive mutants of a bacterial MscL and on MscL homologs with different tension sensitivities are reported, leading to the conclusion that the mode of action of LPC is different from that of applied tension. Our results imply that LPC shifts the free energy of gating by interfering with MscL-membrane coupling. Furthermore, we demonstrate that the fine-tuned addition of LPC can be used for controlled activation of MscL in spectroscopic studies.

Lipid organization of the plasma membrane.

The detailed organization of cellular membranes remains rather elusive. Based on large-scale molecular dynamics simulations, we provide a high-resolution view of the lipid organization of a plasma membrane at an unprecedented level of complexity. Our plasma membrane model consists of 63 different lipid species, combining 14 types of headgroups and 11 types of tails asymmetrically distributed across the two leaflets, closely mimicking an idealized mammalian plasma membrane. We observe an enrichment of cholesterol in the outer leaflet and a general non-ideal lateral mixing of the different lipid species. Transient domains with liquid-ordered character form and disappear on the microsecond time scale. These domains are coupled across the two membrane leaflets. In the outer leaflet, distinct nanodomains consisting of gangliosides are observed. Phosphoinositides show preferential clustering in the inner leaflet. Our data provide a key view on the lateral organization of lipids in one of life's fundamental structures, the cell membrane.

Evidence for cardiolipin binding sites on the membrane-exposed surface of the cytochrome bc1.

The respiratory chain is located in the inner membrane of mitochondria and produces the major part of the ATP used by a cell. Cardiolipin (CL), a double charged phospholipid composing ~10-20% of the mitochondrial membrane, plays an important role in the function and supramolecular organization of the respiratory chain complexes. We present an extensive set of coarse-grain molecular dynamics (CGMD) simulations aiming at the determination of the preferential interfaces of CLs on the respiratory chain complex III (cytochrome bc(1), CIII). Six CL binding sites are identified, including the CL binding sites known from earlier structural studies and buried into protein cavities. The simulations revealed the importance of two subunits of CIII (G and K in bovine heart) for the structural integrity of these internal CL binding sites. In addition, new binding sites are found on the membrane-exposed protein surface. The reproducibility of these binding sites over two species (bovine heart and yeast mitochondria) points to an important role for the function of the respiratory chain. Interestingly the membrane-exposed CL binding sites are located on the matrix side of CIII in the inner membrane and thus may provide localized sources of proton ready for uptake by CIII. Furthermore, we found that CLs bound to those membrane-exposed sites bridge the proteins during their assembly into supercomplexes by sharing the binding sites.

Stability and structure of protein-lipoamino acid colloidal particles: toward nasal delivery of pharmaceutically active proteins.

To circumvent the painful intravenous injection of proteins in the treatment of children with growth deficiency, anemia, and calcium insufficiency, we investigated the stability and structure of protein-lipoamino acid complexes that could be nasally sprayed. Preparations that ensure a colloidal and structural stability of recombinant human growth hormone (rhGH), recombinant human erythropoietin (rhEPO), and salmon calcitonin (sCT) mixed with lauroyl proline (LP) were established. Protein structure was controlled by circular dichroism, and very small sizes of ca. 5 nm were determined by dynamic light scattering. The colloidal preparations could be sprayed with a droplet size of 20-30 µm. The molecular structure of aggregates was investigated by all-atom molecular dynamics. Whereas a lauroyl proline capping of globular proteins rhGH and rhEPO with preservation of their active structure was observed, a mixed micelle of sCT and lipoamino acids was formed. In the latter, aggregated LP constitutes the inner core and the surface is covered with calcitonins that acquire a marked α-helix character. Hydrophobic/philic interaction balance between proteins and LP drives the particles' stability. Passage through nasal cells grown at confluence was markedly increased by the colloidal preparations and could reach a 20 times increase in the case of EPO. Biological implications of such colloidal preparations are discussed in terms of furtiveness.

The flitting of electrons in complex I: a stochastic approach.

A stochastic approach based on the Gillespie algorithm is particularly well adapted to describe the time course of the redox reactions that occur inside the respiratory chain complexes because they involve the motion of single electrons between the individual unique redox centres of a given complex. We use this approach to describe the molecular functioning of the peripheral arm of complex I based on its known crystallographic structure and the rate constants of electron tunnelling derived from the Moser and Dutton phenomenological equations. There are several possible electrons pathways but we show that most of them take the route defined by the successive sites and redox centres: NADH+ site-FMN-N3-N1b-N4-N5-N6a-N6b-N2-Q site. However, the electrons do not go directly from NADH towards the ubiquinone molecule. They frequently jump back and forth between neighbouring redox centres with the result that the net flux of electrons through complex I (i.e. net number of electrons reducing a ubiquinone) is far smaller than the number of redox reactions which actually occur. While most of the redox centres are reduced in our simulations the degree of reduction can vary according to the individual midpoint potentials. The high turnover number observed in our simulation seems to indicate that, in the whole complex I, one or several slower step(s) follow(s) the redox reactions involved in the peripheral arm. It also appears that the residence time of FMNH* and SQ* (possible producers of ROS) is low (around 4% and between 1.6% and 5% respectively according to the values of the midpoint potentials). We did not find any evidence for a role of N7 which remains mainly reduced in our simulations. The role of N1a is complex and depends upon its midpoint potential. In all cases its presence slightly decreases the life time of the flavosemiquinone species. These simulations demonstrate the interest of this type of model which links the molecular physico-chemistry of the individual redox reactions to the more global level of the reaction, as is observed experimentally.

Curvature Energetics Determined by Alchemical Simulation on Four Topologically Distinct Lipid Phases.

The relative curvature energetics of two lipids are tested using thermodynamic integration (TI) on four topologically distinct lipid phases. Simulations use TI to switch between choline headgroup lipids (POPC; that prefers to be flat) and ethanolamine headgroup lipids (POPE; that prefer, for example, the inner monolayer of vesicles). The thermodynamical moving of the lipids between planar, inverse hexagonal (HII), cubic (QII; Pn3m space group), and vesicle topologies reveals differences in material parameters that were previously challenging to access. The methodology allows for predictions of two important lipid material properties: the difference in POPC/POPE monolayer intrinsic curvature (ΔJ0) and the difference in POPC/POPE monolayer Gaussian curvature modulus (Δκ̅m), both of which are connected to the energetics of topological variation. Analysis of the TI data indicates that, consistent with previous experiment and simulation, the J0 of POPE is more negative than POPC (ΔJ0 = -0.018 ± 0.001 Å-1). The theoretical framework extracts significant differences in κ̅m of which POPE is less negative than POPC by 2.0 to 4.0 kcal/mol. The range of these values is determined by considering subsets of the simulations, and disagreement between these subsets suggests separate mechanical parameters at very high curvature. Finally, the fit of the TI data to the model indicates that the position of the pivotal plane of curvature is not constant across topologies at high curvature. Overall, the results offer insights into lipid material properties, the limits of a single HC model, and how to test them using simulation.

Hysteresis and the Cholesterol Dependent Phase Transition in Binary Lipid Mixtures with the Martini Model.

Extensive Martini simulation data, totaling 5 ms, is presented for binary mixtures of dipalmitoylphosphatidylcholine (DPPC) and cholesterol. Using simulation initiated from both gel (so) and liquid-disordered (Ld) phases, significant and strongly cholesterol-dependent hysteresis in the enthalpy as a function of temperature is observed for cholesterol mole fractions from 0 to 20 mol %. Although the precise phase transition temperature cannot be determined due to the hysteresis, the data are consistent with a first order gel to fluid transition, which increases in temperature with cholesterol. At 30 mol % cholesterol, no hysteresis is observed, and there is no evidence for a continuous transition, in either structural parameters like the area per lipid or in the heat capacity as a function of temperature. The results are consistent with a single uniform phase above a critical cholesterol composition between 20 and 30 mol % in Martini, while highlighting the importance and difficulty of obtaining the equilibrium averages to locate phase boundaries precisely in computational models of lipid bilayers.

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.

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.

Helfrich model of membrane bending: from Gibbs theory of liquid interfaces to membranes as thick anisotropic elastic layers.

Helfrich model of membrane bending elasticity has been most influential in establishment and development of Soft-Matter Physics of lipid bilayers and biological membranes. Recently, Helfrich theory has been extensively used in Cell Biology to understand the phenomena of shaping, fusion and fission of cellular membranes. The general background of Helfrich theory on the one hand, and the ways of specifying the model parameters on the other, are important for quantitative treatment of particular biologically relevant membrane phenomena. Here we present the origin of Helfrich model within the context of the general Gibbs theory of capillary interfaces, and review the strategies of computing the membrane elastic moduli based on considering a lipid monolayer as a three-dimensional thick layer characterized by trans-monolayer profiles of elastic parameters. We present the results of original computations of these profiles by a state-of-the-art numerical approach.

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.