Amyloidogenicity as a driving force for the formation of functional oligomers.

Insoluble amyloid fibrils formed by self-assembly of amyloidogenic regions of proteins have a cross-ß-structure. In this work, by using targeted molecular dynamics and rigid body simulation, we demonstrate that if a protein consists of an amyloidogenic region and a globular domain(s) and if the linker between them is short enough, such molecules cannot assemble into amyloid fibrils, instead, they form oligomers with a defined and limited number of ß-strands in the cross-ß core. We show that this blockage of the amyloid growth is due to the steric repulsion of the globular structures linked to amyloidogenic regions. Furthermore, we establish a relationship between the linker length and the number of monomers in such nanoparticles. We hypothesise that such oligomerisation can be a yet unrecognised way to form natural protein complexes involved in biological processes. Our results can also be used in protein engineering for designing soluble nanoparticles carrying different functional domains.

Coarse-Grained Prediction of Peptide Binding to G-Protein Coupled Receptors.

In this study, we used the Martini Coarse-Grained model with no applied restraints to predict the binding mode of some peptides to G-Protein Coupled Receptors (GPCRs). Both the Neurotensin-1 and the chemokine CXCR4 receptors were used as test cases. Their ligands, NTS8-13 and CVX15 peptides, respectively, were initially positioned in the surrounding water box. Using a protocol based on Replica Exchange Molecular Dynamics (REMD), both opening of the receptors and entry of the peptides into their dedicated pockets were observed on the µs time-scale. After clustering, the most statistically representative orientations were closely related to the X-ray structures of reference, sharing both RMSD lower than 3 Å and most of the native contacts. These results demonstrate that such a model, that does not require access to tremendous computational facilities, can be helpful in predicting peptide binding to GPCRs as well as some of the receptor's conformational changes required for this key step. We also discuss how such an approach can now help to predict, de novo, the interactions of GPCRs with other intra- or extracellular peptide/protein partners.

Structural and Molecular Determinants of Membrane Binding by the HIV-1 Matrix Protein.

Assembly of HIV-1 particles is initiated by the trafficking of viral Gag polyproteins from the cytoplasm to the plasma membrane, where they co-localize and bud to form immature particles. Membrane targeting is mediated by the N-terminally myristoylated matrix (MA) domain of Gag and is dependent on the plasma membrane marker phosphatidylinositol-4,5-bisphosphate [PI(4,5)P2]. Recent studies revealed that PI(4,5)P2 molecules containing truncated acyl chains [tr-PI(4,5)P2] are capable of binding MA in an "extended lipid" conformation and promoting myristoyl exposure. Here we report that tr-PI(4,5)P2 molecules also readily bind to non-membrane proteins, including HIV-1 capsid, which prompted us to re-examine MA-PI(4,5)P2 interactions using native lipids and membrane mimetic liposomes and bicelles. Liposome binding trends observed using a recently developed NMR approach paralleled results of flotation assays, although the affinities measured under the equilibrium conditions of NMR experiments were significantly higher. Native PI(4,5)P2 enhanced MA binding to liposomes designed to mimic non-raft-like regions of the membrane, suggesting the possibility that binding of the protein to disordered domains may precede Gag association with, or nucleation of, rafts. Studies with bicelles revealed a subset of surface and myr-associated MA residues that are sensitive to native PI(4,5)P2, but cleft residues that interact with the 2'-acyl chains of tr-PI(4,5)P2 molecules in aqueous solution were insensitive to native PI(4,5)P2 in bicelles. Our findings call to question extended-lipid MA:membrane binding models, and instead support a model put forward from coarse-grained simulations indicating that binding is mediated predominantly by dynamic, electrostatic interactions between conserved basic residues of MA and multiple PI(4,5)P2 and phosphatidylserine molecules.

Discrimination between olfactory receptor agonists and non-agonists.

A joint approach combining free-energy calculations and calcium-imaging assays on the broadly tuned human 1G1 olfactory receptor is reported. The free energy of binding of ten odorants was computed by means of molecular-dynamics simulations. This state function allows separating the experimentally determined eight agonists from the two non-agonists. This study constitutes a proof-of-principle for the computational deorphanization of olfactory receptors.

Coarse-grained simulations of the HIV-1 matrix protein anchoring: revisiting its assembly on membrane domains.

In the accepted model for human immunodeficiency virus preassembly in infected host cells, the anchoring to the intracellular leaflet of the membrane of the matrix domain (MA) that lies at the N-terminus of the viral Gag protein precursor appears to be one of the crucial steps for particle assembly. In this study, we simulated the membrane anchoring of human immunodeficiency virus-1 myristoylated MA protein using a coarse-grained representation of both the protein and the membrane. Our calculations first suggest that the myristoyl group could spontaneously release from its initial hydrophobic pocket before MA protein interacts with the lipid membrane. All-atom simulations confirmed this possibility with a related energy cost estimated to be ~5 kcal.mol(-1). The phosphatidylinositol (4,5) bisphosphate (PI(4,5)P2) head binds preferentially to the MA highly basic region as described in available NMR data, but interestingly without flipping of its 2' acyl chain into the MA protein. Moreover, MA was able to confine PI(4,5)P2 lipids all around its molecular surface after having found a stable orientation at the membrane surface. Our results suggest that this orientation is dependent on Myr anchoring and that this confinement induces a lateral segregation of PI(4,5)P2 in domains. This is consistent with a PI(4,5)P2 enrichment of the virus envelope as compared to the host cell membrane.

Probing the helical forms of Ca2+ -guluronan junction zones in alginate gels by molecular dynamics 1: Duplexes.

A molecular dynamics investigation of the helical forms adopted by (1→4)-α-L-guluronan in explicit water environment was carried out. Single chains and duplexes were modeled at 300 K starting both from 21 or 32 helical conformations and in the presence of a neutralizing amount of Ca(2+) ions. All systems were allowed full conformational freedom. The initial perfect helices with integral screw symmetries were lost at the very beginning of simulations and two distinct behaviors were observed: At equilibrium the 21 models mostly retained the 21 local helical conformations while exploring the 32 ones the rest of the time. In duplexes the two chains, which behaved similarly, were well extended and slightly twisted. By contrast, the chains in 32 duplex models were dissimilar and explored a much broader conformational space in which 21 and 32 local helical conformations were dominant and equally represented but the 31 and other conformations were also present. The wide variety of conformations revealed in this study is consistent with the general difficulty in obtaining crystals of Ca(2+)-guluronate with suitable lateral dimensions for crystallographic studies.

Molecular modelling of odorant/olfactory receptor complexes.

Providing a rationale that associates a chemical structure of an odorant to its induced perception has been sought for a long time. To achieve this, a detailed atomic structure of both the odorant and the olfactory receptor must be known. State-of-the-art techniques to model the 3D structure of an olfactory receptor in complex with various odorants are presented here. These range from sequence alignment with known structures to molecular dynamics simulations in a realistic environment.

Dissociation of membrane-anchored heterotrimeric G-protein induced by G(α) subunit binding to GTP.

Heterotrimeric G-proteins' activation on the intracellular side of the cell membrane is initiated by stimulation of the G-Protein Coupled Receptors (GPCRs) extra-cellular part. This two-step activation mechanism includes (1) an exchange between GDP and GTP molecules in the G(α) subunit and (2) a dissociation of the whole G(αßγ) complex into two membrane-anchored blocks, namely the isolated G(α) and G(ßγ) subunits. Although X-ray data are available for both inactive G(αßγ):GDP and active G(α):GTP complexes, intermediate steps involved in the molecular mechanism of the dissociation have not yet been addressed at the molecular level. In this study, we first built a membrane-anchored intermediate G(iαßγ):GTP complex. This model was then equilibrated by molecular dynamics simulations before the Targeted Molecular Dynamics (TMD) technique was used to force the G(α) subunit to evolve from its inactive (GDP-bound) to its active (GTP-bound) conformations, as described by available X-ray data. The TMD constraint was applied only to the G(α) subunit so that the resulting global rearrangements acting on the whole G(αßγ):GTP heterotrimer could be analyzed. We showed how these mainly local conformational changes of G(α) could initiate large domain:domain motions of the whole complex, the G(ßγ) behaving as an almost quasi-rigid block. This separation of the two G(α):GTP and G(ßγ) subunits required the loss of several interactions at the G(α):G(ßγ) interface that were reported. This study provided an atomistic view of the crucial intermediate step of the G-proteins activation, e.g., the dissociation, that could hardly be elucidated by the experiment.

How broadly tuned olfactory receptors equally recognize their agonists. Human OR1G1 as a test case.

The molecular features that dominate the binding mode of agonists by a broadly tuned olfactory receptor are analyzed through a joint approach combining cell biology, calcium imaging, and molecular modeling. The odorant/receptor affinities, estimated through statistics accrued during molecular dynamics simulations, are in accordance with the experimental ranking. Although in many systems receptors recognize their target through a network of oriented interactions, such as H-bonding, the binding by broadly tuned olfactory receptors is dominated by non-polar terms. We show how such a feature allows chemicals belonging to different chemical families to similarly activate the receptors through compensations of interactions within the binding site.

Molecular modeling of the structural and dynamical properties of secondary plant cell walls: influence of lignin chemistry.

A modeling of lignified secondary plant cell walls adapted to grass has been achieved, using molecular dynamics for time up to 180 ns, applied to systems composed of cellulose, xylan, water, and lignin. The overall model, which was 70 nm thick for a volume of 74.4 nm(3), consisted of two crystalline cellulose layers, each being two molecules deep, separated by an interlayer space where the three other components were located. Whereas the cellulose and xylan chemistry was fixed, 18 lignin systems were considered that varied not only in guaiacyl, syringyl, and p-hydroxyphenyl composition, but also in chain length, linkage types, and the presence or absence of coumaryl units. The stabilized models showed a well-defined interface between xylan and cellulose, but some interpenetration of xylan into the lignin part of the models. A survey of the 18 models showed that their lignin component was amorphous and that their density profile was very variable and essentially model dependent. This variability was also reflected in the co-orientation of the lignin phenyl rings with respect to the cellulose surfaces, some systems showing some orientation whereas others did not. The pattern of void distribution accessible to water varied from one system to the next, but the overall void volume was systematically established at around 3%, accepting around 200 water molecules. The estimated mobility of the water molecules interacting with lignin was 1.5 times greater than that interacting with carbohydrates.

Deciphering the selectivity of Bombyx mori pheromone binding protein for bombykol over bombykal: a theoretical approach.

In this article we report calculations dedicated to estimate the selectivity of the Bombyx mori pheromone binding protein towards the two closely related pheromonal components Bombykol and Bombykal. The selectivity is quantified by the binding free-energy difference, obtained either by the thermodynamic integration or by the MM-GBSA approach. In the latter, the selectivity is decomposed on a per-residue basis, which identifies the residues considered crucial for the selectivity of the protein for Bombykol over Bombykal. A discussion on the role of Bombyx mori pheromone binding protein is provided on the basis of these results.

Binding free energy prediction in strongly hydrophobic biomolecular systems.

We present a comparison of various computational approaches aiming at predicting the binding free energy in ligand-protein systems where the ligand is located within a highly hydrophobic cavity. The relative binding free energy between similar ligands is obtained by means of the thermodynamic integration (TI) method and compared to experimental data obtained through isothermal titration calorimetry measurements. The absolute free energy of binding prediction was obtained on a similar system (a pyrazine derivative bound to a lipocalin) by TI, potential of mean force (PMF) and also by means of the MMPBSA protocols. Although the TI protocol performs poorly either with an explicit or an implicit solvation scheme, the PMF calculation using an implicit solvation scheme leads to encouraging results, with a prediction of the binding affinity being 2 kcal mol(-1) lower than the experimental value. The use of an implicit solvation scheme appears to be well suited for the study of such hydrophobic systems, due to the lack of water molecules within the binding site.