Modeling of the structure in aqueous solution of the exopolysaccharide produced by Lactobacillus helveticus 766.

A method is described for constructing a conformational model in water of a heteropolysaccharide built up from repeating units, and is applied to the exopolysaccharide produced by Lactobacillus helveticus 766. The molecular modeling method is based on energy minima, obtained from molecular mechanics calculations of each of the constituting disaccharide fragments of the repeating unit in vacuo, as starting points. Subsequently, adaptive umbrella sampling of the potential of mean force is applied to extract rotamer populations of glycosidic dihedral angles of oligosaccharide fragments in solution. From these analyses, the most probable conformations are constructed for the hexasaccharide-repeating unit of the polysaccharide. After exploring the conformational space of each of the individual structures by molecular dynamics simulations, the different repeating unit conformations are used as building blocks for the generation of oligo- and polysaccharide models, by using a polysaccharide building program. The created models of the exopolysaccharide produced by L. helveticus 766 exhibit a flexible twisted secondary structure and tend to adopt a random coil conformation as tertiary structure.     

Role of glycosylation in structure and stability of Erythrina corallodendron lectin (EcorL): a molecular dynamics study

The effect of glycosylation on structure and stability of glycoproteins has been a topic of considerable interest. In this work, we have investigated the solution conformation of the oligosaccharide and its effect on the structure and stability of the glycoprotein by carrying out a series of long Molecular dynamics (MD) simulations on glycosylated Erythrina corallodendron lectin (EcorL) and nonglycosylated recombinant Erythrina corallodendron lectin (rEcorL). Our results indicate that, despite the similarity in overall three dimensional structures, glycosylated EcorL has lesser nonpolar solvent accessible surface area compared to nonglycosylated EcorL. This might explain the experimental observation of higher thermodynamic stability for glycosylated EcorL compared to nonglycosylated EcorL. Analysis of the simulation results indicates that, dynamic view of interactions between protein residues and oligosaccharide is entirely different from the static picture seen in the crystal structure. The oligosaccharide moiety had dynamically stable interactions with Lys 55 and Tyr 53, both of which are separated in sequence from the site of glycosylation, Asn 17. It is possible that glycosylation helps in forming long-range contacts between amino acids, which are separated in sequence and thus provides a folding nucleus. Thus our simulations not only reveal the conformations sampled by the oligosaccharide, but also provide novel insights into possible molecular mechanisms by which glycosylation can help in folding of the glycoprotein by formation of folding nucleus involving specific contacts with the oligosaccharide moiety.     

Conformational study of cyclolinopeptide A. A distance geometry and molecular dynamics approach

The conformation of cyclolinopeptide A, c(Pro-Pro-Phe-Phe-Leu-Ile-Ile-Leu-Val), a naturally occurring peptide with remarkable cytoprotective activity, has been investigated by means of distance geometry calculations and molecular dynamics simulations. The starting points for all the calculations were an X-ray structure and other structures obtained from distance geometry calculations based on NMR data. Restrained and unrestrained molecular dynamics simulations are reported in vacuo and in CCl4. Structural and dynamic properties are investigated and compared with those experimentally determined. The conformation obtained from the MD simulations which best reproduces the NMR parameters is at the same time one of the most stable ones and is also fairly similar to the crystal structure. An explanation for the occurrence of multiple conformations in solution at room temperature is given.     

The MUMO (minimal under-restraining minimal over-restraining) method for the determination of native state ensembles of proteins

While reliable procedures for determining the conformations of proteins are available, methods for generating ensembles of structures that also reflect their flexibility are much less well established. Here we present a systematic assessment of the ability of ensemble-averaged molecular dynamics simulations with ensemble-averaged NMR restraints to simultaneously reproduce the average structure of proteins and their associated dynamics. We discuss the effects that under-restraining (overfitting) and over-restraining (underfitting) have on the structures generated in ensemble-averaged molecular simulations. We then introduce the MUMO (minimal under-restraining minimal over-restraining) method, a procedure in which different observables are averaged over a different number of molecules. As both over-restraining and under-restraining are significantly reduced in the MUMO method, it is possible to generate ensembles of conformations that accurately characterize both the structure and the dynamics of native states of proteins. The application of the MUMO method to the protein ubiquitin yields a high-resolution structural ensemble with an RDC Q-factor of 0.19.     

A molecular dynamics study of the 41-56 beta-hairpin from B1 domain of protein G.

The structural and dynamical behavior of the 41-56 beta-hairpin from the protein G B1 domain (GB1) has been studied at different temperatures using molecular dynamics (MD) simulations in an aqueous environment. The purpose of these simulations is to establish the stability of this hairpin in view of its possible role as a nucleation site for protein folding. The conformation of the peptide in the crystallographic structure of the protein GB1 (native conformation) was lost in all simulations. The new equilibrium conformations are stable for several nanoseconds at 300K (>10 ns), 350 K (>6.5 ns), and even at 450 K (up to 2.5 ns). The new structures have very similar hairpin-like conformations with properties in agreement with available experimental nuclear Overhauser effect (NOE) data. The stability of the structure in the hydrophobic core region during the simulations is consistent with the experimental data and provides further evidence for the role played by hydrophobic interactions in hairpin structures. Essential dynamics analysis shows that the dynamics of the peptide at different temperatures spans basically the same essential subspace. The main equilibrium motions in this subspace involve large fluctuations of the residues in the turn and ends regions. Of the six interchain hydrogen bonds, the inner four remain stable during the simulations. The space spanned by the first two eigenvectors, as sampled at 450 K, includes almost all of the 47 different hairpin structures found in the database. Finally, analysis of the hydration of the 300 K average conformations shows that the hydration sites observed in the native conformation are still well hydrated in the equilibrium MD ensemble.     

Solvation Effects on the Conformational Behaviour of Gellan and Calcium Ion Binding to Gellan Double Helices

The contribution of the presence of solvent to the conformations adopted by disaccharide fragments within the repeat unit of gellan have been studied by molecular modelling techniques. Initial conformational energy searches, using a dielectric continuum to represent the solvent, provided starting geometries for a series of molecular dynamics simulations. The solution behaviour from these simulations was subsequently compared to fibre diffraction data of the potassium gellan salt. The present calculations indicate considerable flexibility of the glycosidic linkages, and this is discussed in relation to its effect on gel formation. One of the fragments was solvated with explicit water molecules. These calculations showed the same conformational behaviour as those simulations conducted in implicit solvent.Finally, a series of molecular dynamics (MD) simulations were performed to study the calcium binding to gellan. The results from this clearly showed a well defined binding site for this ion.Abbreviations MM molecular mechanics - MD molecular dynamics     

Effects of phosphorylation on the intrinsic propensity of backbone conformations of serine/threonine

Each amino acid has its intrinsic propensity for certain local backbone conformations, which can be further modulated by the physicochemical environment and post-translational modifications. In this work, we study the effects of phosphorylation on the intrinsic propensity for different local backbone conformations of serine/threonine by molecular dynamics simulations. We showed that phosphorylation has very different effects on the intrinsic propensity for certain local backbone conformations for the serine and threonine. The phosphorylation of serine increases the propensity of forming polyproline II, whereas that of threonine has the opposite effect. Detailed analysis showed that such different responses to phosphorylation mainly arise from their different perturbations to the backbone hydration and the geometrical constraints by forming side-chain–backbone hydrogen bonds due to phosphorylation. Such an effect of phosphorylation on backbone conformations can be crucial for understanding the molecular mechanism of phosphorylation-regulated protein structures/dynamics and functions.     

Development and testing of protocols for computer-aided design of peptide drugs, using oxytocin

Considerable clinical interest in neuropeptides and peptide hormones has stimulated recent research and development of peptide-based drugs. This process differs from most classical drug discovery procedures because peptide molecules have considerable inherent flexibility. In the present paper, to identify lowest energy and metastable conformers for drug design, and to develop protocols for such studies, conformational search algorithms, incorporating empirical energy calculations, have been applied in the analysis of the peptide oxytocin. Energy minimization in torsion angle space was carried out from a variety of starting conformations, including published structures, in all-atom mode and all with distance constraints for disulphide bond formation. The energy-minimized conformations have been further optimized by a mapping method. Complementary simulations have been performed in united-atom mode and a model representing the effects of water using dummy sites has been developed and tested for this representation. Several of the preferred conformers together with de novo conformations have been used as starting points in molecular dynamics simulations; 28 low potential energy conformations were located at a temperature of 4 K. Conformations are analysed to identify hydrogen bonds, phi-psi angle distributions and the RMS values relative to the X-ray structure of deamino-oxytocin. The modelled structure of lowest energy in the molecular mechanics calculations was also that of least RMS deviation from the crystal structure; whilst structures of lower energy but larger deviation were identified by molecular dynamics techniques. A metastable structure has been identified which satisfies existing criteria for the "active form", and this model is tested by a theoretical residue-substitution technique, to provide clues on the agonist/antagonist relationship at the atomic level.     

Membrane-bound ARF1 peptide: interpretation of neutron diffraction data by molecular dynamics simulation methods

Adenosine diphosphate ribosylation factor-1 (ARF1) is activated by cell membrane binding of a self-folding N-terminal domain. We have previously presented four possible conformations of the membrane bound, human ARF1 N-terminal peptide in planar lipid bilayers of DOPC and DOPG (7:3 molar ratio), determined from lamellar neutron diffraction and circular dichroism data. In this paper we analyse the four possible conformations by molecular dynamics simulations. The aim of these simulations was to use MD to distinguish which of the four possible membrane bound structures was the most likely. The most likely conformation was determined according to the following criteria: (a) location of label positions on the peptide in relation to the bilayer, (b) lowest mean square displacement from the initial structure, (c) lowest system energy, (d) most peptide-lipid headgroup hydrogen bonding, (e) analysis of phi/psi angles of the peptide. These findings demonstrate the application of molecular dynamics simulations to explore neutron diffraction data.     

Similarity Measures for Protein Ensembles

Analyses of similarities and changes in protein conformation can provide important information regarding protein function and evolution. Many scores, including the commonly used root mean square deviation, have therefore been developed to quantify the similarities of different protein conformations. However, instead of examining individual conformations it is in many cases more relevant to analyse ensembles of conformations that have been obtained either through experiments or from methods such as molecular dynamics simulations. We here present three approaches that can be used to compare conformational ensembles in the same way as the root mean square deviation is used to compare individual pairs of structures. The methods are based on the estimation of the probability distributions underlying the ensembles and subsequent comparison of these distributions. We first validate the methods using a synthetic example from molecular dynamics simulations. We then apply the algorithms to revisit the problem of ensemble averaging during structure determination of proteins, and find that an ensemble refinement method is able to recover the correct distribution of conformations better than standard single-molecule refinement.     

Homogeneous and heterogeneous tertiary structure ensembles of amyloid-β peptides

The interplay of modern molecular simulation and high-quality nuclear magnetic resonance (NMR) experiments has reached a fruitful stage for quantitative characterization of structural ensembles of disordered peptides. Amyloid-β 1-42 (Aβ42), the primary peptide associated with Alzheimer's disease, and fragments such as Aβ21-30 are both classified as intrinsically disordered peptides (IDPs). We use a variety of NMR observables to validate de novo molecular dynamics simulations in explicit water to characterize the tertiary structure ensemble of Aβ42 and Aβ21-30 from the perspective of their classification as IDPs. Unlike the Aβ21-30 fragment that conforms to expectations of an IDP that is primarily extended, we find that Aβ42 samples conformations reflecting all possible secondary structure categories and spans the range of IDP classifications from collapsed structured states to highly extended conformations, making it an IDP with a far more heterogeneous tertiary ensemble.     

Reversible scaling of dihedral angle barriers during molecular dynamics to improve structure prediction of cyclic peptides.

Peptide cyclization or chemical cross-linking has frequently been used to restrict the conformational freedom of a peptide, for example, to enhance its capacity for selective binding to a target receptor molecule. Structure prediction of cyclic peptides is important to evaluate possible conformations prior to synthesis. Because of the conformational constraints imposed by cyclization low energy conformations of cyclic peptides can be separated by large energy barriers. In order to improve the conformational search properties of molecular dynamics (MD) simulations a potential scaling method has been designed. The approach consists of several consecutive MD simulations with a specific lowering of dihedral energy barriers and reduced nonbonded interactions between atoms separated by three atoms followed by gradually scaling the potential until the original barriers are reached. Application to four cyclic penta- and hexa-peptide test cases and a protein loop of known structure indicates that the potential scaling method is more efficient and faster in locating low energy conformations than standard MD simulations. Combined with a generalized Born implicit solvation model the low energy cyclic peptide conformations and the loop structure are in good agreement with experiment. Applications in the presence of explicit water molecules during the simulations showed also improved convergence to structures close to experiment compared with regular MD.     

Conformational dynamics of sialyl Lewisx in aqueous solution and its interaction with selectinE. A study by molecular dynamics

Three dimensional structures of sialyl Lewis(x) (SLe(x)) in aqueous solution and bound to selectinE are described based on an exhaustive conformational analysis and several long molecular dynamics simulations using different glycosidic regions as starting conformations. It appears from this study that when the oligosaccharide is free in solution the NeuNAcalpha(2-3)Gal segment favors glycosidic conformation in three different regions in the (Phi,Psi) plane with propensity of populations in the ratio 1:8:1. Each one of these structures is characteristically stabilized by specific hydrogen bonding interaction between NeuNAc and Gal. On the other hand, the Gal-GlcNAc-Fuc segment can exist in four different conformational states. Based on the topology of SLe(x) we are able to predict that out of all the allowed conformations in solution only two of these structures possess a geometry that would fit without steric clashes into the binding location of selectinE. In both of these binding modes, segment Gal-GlcNAc-Fuc adopts a unique conformation. The only difference between the two SLe(x) conformers that can successfully bind to selectinE is given by two possible regions in glycosidic space in the fragment NeuNAcalpha(2-3)Gal. A large conformational departure from the crystallographic data is observed for two lysine residues at the binding site of selectinE. These two residues play an important role when SLe(x) binds selectinE in aqueous solution. These findings help reconcile the X-ray data, in which these residues appear to be 1 nm away from SLe(x), with recent liquid NMR data reporting couplings between these protein residues and the sugar.     

Membrane-bound ARF1 peptide: interpretation of neutron diffraction data by molecular dynamics simulation methods

Adenosine diphosphate ribosylation factor-1 (ARF1) is activated by cell membrane binding of a self-folding N-terminal domain. We have previously presented four possible conformations of the membrane bound, human ARF1 N-terminal peptide in planar lipid bilayers of DOPC and DOPG (7:3 molar ratio), determined from lamellar neutron diffraction and circular dichroism data. In this paper we analyse the four possible conformations by molecular dynamics simulations. The aim of these simulations was to use MD to distinguish which of the four possible membrane bound structures was the most likely. The most likely conformation was determined according to the following criteria: (a) location of label positions on the peptide in relation to the bilayer, (b) lowest mean square displacement from the initial structure, (c) lowest system energy, (d) most peptide-lipid headgroup hydrogen bonding, (e) analysis of phi/psi angles of the peptide. These findings demonstrate the application of molecular dynamics simulations to explore neutron diffraction data.     

Discrete analyses of protein dynamics

Abstract

Protein structures are highly dynamic macromolecules. This dynamics is often analysed through experimental and/or computational methods only for an isolated or a limited number of proteins. Here, we explore large-scale protein dynamics simulation to observe dynamics of local protein conformations using different perspectives. We analysed molecular dynamics to investigate protein flexibility locally, using classical approaches such as RMSf, solvent accessibility, but also innovative approaches such as local entropy. First, we focussed on classical secondary structures and analysed specifically how β-strand, β–turns, and bends evolve during molecular simulations. We underlined interesting specific bias between β–turns and bends, which are considered as the same category, while their dynamics show differences. Second, we used a structural alphabet that is able to approximate every part of the protein structures conformations, namely protein blocks (PBs) to analyse (i) how each initial local protein conformations evolve during dynamics and (ii) if some exchange can exist among these PBs. Interestingly, the results are largely complex than simple regular/rigid and coil/flexible exchange. Abbreviations N_eq number of equivalent

PB Protein Blocks

PDB Protein DataBank

RMSf root mean square fluctuations

Communicated by Ramaswamy H. Sarma     

Simulation of the structure and dynamics of nonhelical RNA motifs

Computer simulation methods are increasingly being used to study possible conformations and dynamics of structural motifs in RNA. Recent results of molecular dynamics simulations and continum solvent studies of RNA structures and RNA-ligand complexes show promising agreement with experimental data. Combined with the ongoing progress in the experimental characterization of RNA structure and thermodynamics, these computational approaches can help to better understand the mechanism of RNA structure formation and the binding of ligands.     

Flexible fitting of atomic structures into electron microscopy maps using molecular dynamics

A novel method to flexibly fit atomic structures into electron microscopy (EM) maps using molecular dynamics simulations is presented. The simulations incorporate the EM data as an external potential added to the molecular dynamics force field, allowing all internal features present in the EM map to be used in the fitting process, while the model remains fully flexible and stereochemically correct. The molecular dynamics flexible fitting (MDFF) method is validated for available crystal structures of protein and RNA in different conformations; measures to assess and monitor the fitting process are introduced. The MDFF method is then used to obtain high-resolution structures of the E. coli ribosome in different functional states imaged by cryo-EM.     

Restricted N-glycan Conformational Space in the PDB and Its Implication in Glycan Structure Modeling

Understanding glycan structure and dynamics is central to understanding protein-carbohydrate recognition and its role in protein-protein interactions. Given the difficulties in obtaining the glycan''s crystal structure in glycoconjugates due to its flexibility and heterogeneity, computational modeling could play an important role in providing glycosylated protein structure models. To address if glycan structures available in the PDB can be used as templates or fragments for glycan modeling, we present a survey of the N-glycan structures of 35 different sequences in the PDB. Our statistical analysis shows that the N-glycan structures found on homologous glycoproteins are significantly conserved compared to the random background, suggesting that N-glycan chains can be confidently modeled with template glycan structures whose parent glycoproteins share sequence similarity. On the other hand, N-glycan structures found on non-homologous glycoproteins do not show significant global structural similarity. Nonetheless, the internal substructures of these N-glycans, particularly, the substructures that are closer to the protein, show significantly similar structures, suggesting that such substructures can be used as fragments in glycan modeling. Increased interactions with protein might be responsible for the restricted conformational space of N-glycan chains. Our results suggest that structure prediction/modeling of N-glycans of glycoconjugates using structure database could be effective and different modeling approaches would be needed depending on the availability of template structures.     

The Influence of 150-Cavity Binders on the Dynamics of Influenza A Neuraminidases as Revealed by Molecular Dynamics Simulations and Combined Clustering

Neuraminidase inhibitors are the main pharmaceutical agents employed for treatments of influenza infections. The neuraminidase structures typically exhibit a 150-cavity, an exposed pocket that is adjacent to the catalytic site. This site offers promising additional contact points for improving potency of existing pharmaceuticals, as well as generating entirely new candidate inhibitors. Several inhibitors based on known compounds and designed to interact with 150-cavity residues have been reported. However, the dynamics of any of these inhibitors remains unstudied and their viability remains unknown. This work reports the outcome of long-term, all-atom molecular dynamics simulations of four such inhibitors, along with three standard inhibitors for comparison. Each is studied in complex with four representative neuraminidase structures, which are also simulated in the absence of ligands for comparison, resulting in a total simulation time of 9.6µs. Our results demonstrate that standard inhibitors characteristically reduce the mobility of these dynamic proteins, while the 150-binders do not, instead giving rise to many unique conformations. We further describe an improved RMSD-based clustering technique that isolates these conformations – the structures of which are provided to facilitate future molecular docking studies – and reveals their interdependence. We find that this approach confers many advantages over previously described techniques, and the implications for rational drug design are discussed.     

Interaction of anesthetics with open and closed conformations of a potassium channel studied via molecular dynamics and normal mode analysis

A variety of experiments suggest that membrane proteins are important targets of anesthetic molecules, and that ion channels interact differently with anesthetics in their open and closed conformations. The availability of an open and a closed structural model for the KirBac1.1 potassium channel has made it possible to perform a comparative analysis of the interactions of anesthetics with the same channel in its open and closed states. To this end, all-atom molecular dynamics simulations supplemented by normal mode analysis have been employed to probe the interactions of the inhalational anesthetic halothane with both an open and closed conformer of KirBac1.1 embedded in a lipid bilayer. Normal mode analysis on the closed and open channel, in the presence and absence of halothane, reveals that the anesthetic modulates the global as well as the local dynamics of both conformations differently. In the case of the open channel, the observed reduction of flexibility of residues in the inner helices suggests a functional modification action of anesthetics on ion channels. In this context, preferential quenching of the aromatic residue motion and modulation of global dynamics by halothane may be seen as steps toward potentiating or favoring open state conformations. These molecular dynamics simulations provide the first insights into possible specific interactions between anesthetic molecules and ion channels in different conformations.