Search

Metamorphic, or fold-switching, proteins feature different folds that are physiologically relevant. The human chemokine XCL1 (or Lymphotactin) is a metamorphic protein that features two native states, an $$\alpha -\beta$$ α - β and an all$$-\beta$$ - β fold, which have similar stability at physiological condition. Here, extended molecular dynamics (MD) simulations, principal component analysis of atomic fluctuations and thermodynamic modeling based on both the configurational volume and free energy landscape, are used to obtain a detailed characterization of the conformational thermodynamics of human Lymphotactin and of one of its ancestors (as was previously obtained by genetic reconstruction). Comparison of our computational results with the available experimental data show that the MD-based thermodynamics can explain the experimentally observed variation of the conformational equilibrium between the two proteins. In particular, our computational data provide an interpretation of the thermodynamic evolution in this protein, revealing the relevance of the configurational entropy and of the shape of the free energy landscape within the essential space (i.e., the space defined by the generalized internal coordinates providing the largest, typically non-Gaussian, structural fluctuations).

Background Essential Dynamics (ED) is a common application of principal component analysis (PCA) to extract biologically relevant motions from atomic trajectories of proteins. Covariance and correlation based PCA are two common approaches to determine PCA modes (eigenvectors) and their eigenvalues. Protein dynamics can be characterized in terms of Cartesian coordinates or internal distance pairs. In understanding protein dynamics, a comparison of trajectories taken from a set of proteins for similarity assessment provides insight into conserved mechanisms. Comprehensive software is needed to facilitate comparative-analysis with user-friendly features that are rooted in best practices from multivariate statistics. Results We developed a Java based Essential Dynamics toolkit called JED to compare the ED from multiple protein trajectories. Trajectories from different simulations and different proteins can be pooled for comparative studies. JED implements Cartesian-based coordinates (cPCA) and internal distance pair coordinates (dpPCA) as options to construct covariance (Q) or correlation (R) matrices. Statistical methods are implemented for treating outliers, benchmarking sampling adequacy, characterizing the precision of Q and R, and reporting partial correlations. JED output results as text files that include transformed coordinates for aligned structures, several metrics that quantify protein mobility, PCA modes with their eigenvalues, and displacement vector (DV) projections onto the top principal modes. Pymol scripts together with PDB files allow movies of individual Q- and R-cPCA modes to be visualized, and the essential dynamics occurring within user-selected time scales. Subspaces defined by the top eigenvectors are compared using several statistical metrics to quantify similarity/overlap of high dimensional vector spaces. Free energy landscapes can be generated for both cPCA and dpPCA. Conclusions JED offers a convenient toolkit that encourages best practices in applying multivariate statistics methods to perform comparative studies of essential dynamics over multiple proteins. For each protein, Cartesian coordinates or internal distance pairs can be employed over the entire structure or user-selected parts to quantify similarity/differences in mobility and correlations in dynamics to develop insight into protein structure/function relationships. Electronic supplementary material The online version of this article (doi:10.1186/s12859-017-1676-y) contains supplementary material, which is available to authorized users.

Molecular dynamics simulations of the wild type bovine ADP/ATP mitochondrial carrier, and of the single Ala113Pro and double Ala113Pro/Val180Met mutants, embedded in a lipid bilayer, have been carried out for 30 ns to shed light on the structural-dynamical changes induced by the Val180Met mutation restoring the carrier function in the Ala113Pro pathologic mutant. Principal component analysis indicates that, for the three systems, the protein dynamics is mainly characterized by the motion of the matrix loops and of the odd-numbered helices having a conserved proline in their central region. Analysis of the motions shows a different behaviour of single pathological mutant with respect of the other two systems. The single mutation induces a regularization and rigidity of the H3 helix, lost upon the introduction of the second mutation. This is directly correlated to the salt bridge distribution involving residues Arg79, Asp134 and Arg234, hypothesized to interact with the substrate. In fact, in the wild type simulation two stable inter-helices salt bridges, crucial for substrate binding, are present almost over all the simulation time. In line with the impaired ADP transport, one salt interaction is lost in the single mutant trajectory but reappears in the double mutant simulation, where a salt bridge network matching the wild type is restored. Other important structural-dynamical properties, such as the trans-membrane helices mobility, analyzed via the principal component analysis, are similar for the wild type and double mutant while are different for the single mutant, providing a mechanistic explanation for their different functional properties.

The main purpose of this work is to analyse, by means of molecular dynamics (MD) simulations both structural and mechanical-dynamical differences between Hepcidin-20 and Hepcidin-25 in both oxidized and reduced states in aqueous solution. Results indicate that the presence of disulfide bonds is essential, in both peptides, for maintaining their β-hairpin motif. As a matter of fact, the lack of this intra-peptide covalent interactions produces an almost immediate deviation from the oxidized, plausibly active, structure in both the systems. Interestingly, reduced Hepcidin-25 turns out to be characterized by a highly fluctuating structure which is found to rapidly span a large number of configurations at equilibrium. On the other hand, loss of disulfide bonds in the shorter peptide, results in a more compact and relatively rigid double-turn structure. Comparison of mechanical–dynamical properties and sidechains–sidechains interactions in oxidized Hepcidin-20 and Hepcidin-25 strongly suggest also the key role of N-terminus in the aggregation tendency of Hepcidin-25. © 2010 Wiley Periodicals, Inc. Biopolymers 93: 917–926, 2010.

Aim of this work is to assess the informativeness of protein dynamics in the detection of similarities among distant homologous proteins. To this end, an approach to perform large-scale comparisons of protein domain flexibilities is proposed. CONCOORD is confirmed as a reliable method for fast conformational sampling. The root mean square fluctuation of alpha carbon positions in the essential dynamics subspace is employed as a measure of local flexibility and a synthetic index of similarity is presented. The dynamics of a large collection of protein domains from ASTRAL/SCOP40 is analyzed and the possibility to identify relationships, at both the family and the superfamily levels, on the basis of the dynamical features is discussed. The obtained picture is in agreement with the SCOP classification, and furthermore suggests the presence of a distinguishable familiar trend in the flexibility profiles. The results support the complementarity of the dynamical and the structural information, suggesting that information from dynamics analysis can arise from functional similarities, often partially hidden by a static comparison. On the basis of this first test, flexibility annotation can be expected to help in automatically detecting functional similarities otherwise unrecoverable. © 2007 The Author(s).

In this work we report the study of a peptide, the Contryphan Vn produced by Conus ventricosus, a vermivorous cone snail living in the temperate Mediterranean sea. This cyclic peptide of nine residues is a ring closed by a Cys-Cys (Cys: cysteine) disulfide bond containing two proline (Pro) residues and two tryptophans (Trp), one of them being a D-Trp. We present a statistical mechanical characterization of the peptide, simulated in water for about 200 ns with classical molecular dynamics (MD). In recent years there has been a growing interest in the study of the mechanics and dynamics of biological molecules, and in particular for proteins and peptides, about the relationship between collective motions and the active conformations which exert the biological function. To this aim we used the essential dynamics analysis on the MD trajectory and extracted, from the total fluctuations of the molecule, the dominant dynamical modes responsible of the principal conformational transitions. The Contryphan Vn small size allowed us to investigate in details the all-atoms dynamics and the corresponding thermodynamics in conformational space defined by the most significant intramolecular motions. © 2004 Wiley Periodicals, Inc. Biopoly- mers 74: 448 - 456, 2004

The implementation of cutinase from Fusarium solani pisi as a fat-stain removing ingredient in laundry washing is hampered by its unfolding in the presence of anionic surfactants. In this work we present molecular dynamics (MD) computer simulations on cutinase and analysis procedures to distinguish the movements related to its functional behavior (e.g., substrate binding) from those related to the unfolding of the enzyme. Two kinds of MD-simulations were performed: a simulation mimicking the thermal motion at room temperature, and several simulations in which unfolding is induced either by high temperature or by using a modified water-protein interaction potential. Essential dynamics analyses (A. Amadei et al., Proteins 17:412–425, 1993) on the simulations identify distinct regions in the molecular structure of cutinase in which the motions occur for function and initial unfolding. The unfolding in various simulations starts in a similar way, suggesting that mutations in the regions involved might stabilize the enzyme without affecting its functionality. Proteins 33:253–264, 1998. © 1998 Wiley-Liss, Inc.

Model-free methods are introduced to determine quantities pertaining to protein domain motions from normal mode analyses and molecular dynamics simulations, For the normal mode analysis, the methods are based on the assumption that in low frequency modes, domain motions can be well approximated by modes of motion external to the domains, To analyze the molecular dynamics trajectory, a principal component analysis tailored specifically to analyze interdomain motions is applied, A method based on the curl of the atomic displacements is described, which yields a sharp discrimination of domains, and which defines a unique interdomain screw-axis, Hinge axes are defined and classified as twist or closure axes depending on their direction, The methods have been tested on lysozyme, A remarkable correspondence was found between the first normal mode axis and the first principal mode axis, with both axes passing within 3 Angstrom of the alpha-carbon atoms of residues 2, 39, and 56 of human lysozyme, and near the interdomain helix, The axes of the first modes are overwhelmingly closure axes, A lesser degree of correspondence is found for the second modes, but in both cases they are more twist axes than closure axes, Both analyses reveal that the interdomain connections allow only these two degrees of freedom, one more than provided by a pure mechanical hinge. (C) 1997 Wiley-Liss, Inc.

The Gal4p mediated transcriptional activation of GAL genes requires the interaction between Gal3p bound with ATP and galactose and Gal80p. Though numerous studies suggest that galactose and ATP activate Gal3p/Gal1p interaction with Gal80p, neither the mechanism of activation nor the interacting surface that binds to Gal80p is well understood. In this study we investigated the dynamics of Gal3p and Gal1p in the presence and absence of ligands ATP and galactose to understand the role played by dynamics in the function of these proteins through molecular dynamics simulation and protein–protein docking studies. We performed simulations totaling to 510 ns on both Gal1p and Gal3p proteins in the presence and absence of ligands ATP and galactose. We find that, while binding of ligands ATP and galactose to Gal3p/Gal1p do not affect the global conformation of proteins, some local conformational changes around upper-lip helix including insertion domain are observed. We observed that only in the presence of ATP and galactose, Gal3p displays opening and closing motion between the two domains. And because of this motion, a binding interface, which is largely hydrophobic, opens up on the surface of Gal3p and this surface can bind to Gal80p. From our simulation studies we infer probable docking sites for Gal80p on Gal3p/Gal1p, which were further ascertained by the docking of Gal80p on to ligand bound Gal1p and Gal3p proteins, and the residues at the interface between Gal3p and Gal80p are identified. Our results correlate quite well with the existing body of literature on functional and dynamical aspects of Gal1p and Gal3p proteins.

NMR structures are typically deposited in databases such as the PDB in the form of an ensemble of structures. Generally, each of the models in such an ensemble satisfies the experimental data and is equally valid. No unique solution can be calculated because the experimental NMR data is insufficient, in part because it reflects the conformational variability and dynamical behavior of the molecule in solution. Even for relatively rigid molecules, the limited number of structures that are typically deposited cannot completely encompass the structural diversity allowed by the observed NMR data, but they can be chosen to try and maximize its representation. We describe here the adaptation and application of techniques more commonly used to examine large ensembles from molecular dynamics simulations, to the analysis of NMR ensembles. The approach, which is based on principal component analysis, we call COCO ("Complementary Coordinates"). The COCO approach analyses the distribution of an NMR ensemble in conformational space, and generates a new ensemble that fills "gaps" in the distribution. The method is very rapid, and analysis of a 25-member ensemble and generation of a new 25 member ensemble typically takes 1-2 min on a conventional workstation. Applied to the 545 structures in the RECOORD database, we find that COCO generates new ensembles that are as structurally diverse-both from each other and from the original ensemble-as are the structures within the original ensemble. The COCO approach does not explicitly take into account the NMR restraint data, yet in tests on selected structures from the RECOORD database, the COCO ensembles are frequently good matches to this data, and certainly are structures that can be rapidly refined against the restraints to yield high-quality, novel solutions. COCO should therefore be a useful aid in NMR structure refinement and in other situations where a richer representation of conformational variability is desired-for example in docking studies. COCO is freely accessible via the website www.ccpb.ac.uk/COCO.

Periplasmic binding proteins from Gram-negative bacteria possess a common architecture, comprised of two domains linked by a hinge region, a fold which they share with the neurotransmitter-binding domains of ionotropic glutamate receptors (GluRs). Glutamine-binding protein (GlnBP) is one such protein, whose crystal structure has been solved in both open and closed forms. Multi-nanosecond molecular dynamics simulations have been used to explore motions about the hinge region and how they are altered by ligand binding. Glutamine binding is seen to significantly reduce inter-domain motions about the hinge region. Essential dynamics analysis of inter-domain motion revealed the presence of both hinge-bending and twisting motions, as has been reported for a related sugar-binding protein. Significantly, the influence of the ligand on GlnBP dynamics is similar to that previously observed in simulations of rat glutamate receptor (GluR2) ligand-binding domain. The essential dynamics analysis of GlnBP also revealed a third class of motion which suggests a mechanism for signal transmission in GluRs.

Human neutrophil collagenase (HNC,MMP-8) is one of the target enzymes for drug treatment of pathologic extracellular matrix degradation. Peptidomimetic inhibitors bind in the S'-side of the enzyme active site occupying the S1' primary specificity pocket by their large hydrophobic side-chains. The crystal structure of the complex between the catalytic domain of MMP-8 and Pro-Leu-TrpP(OH)2 (PLTP) showed that this phosphonate inhibitor binds in the S side of the active site. This finding was unexpected since it represents the first example of accomodation of the bulky Trp indolyl chain in the S1 rather than in the S1' subsite. Dynamical and structural factors favouring this uncommon mode of binding were therefore investigated. Molecular dynamics (MD) simulations performed on the uncomplexed enzyme show that its structure in aqueous solution is only slightly different from the crystal structure found in the complex with PLTP. Essential dynamics analysis of the MD simulations, performed on PLTP alternatively interacting with the S- or S'-side of the active site, shows that the enzyme fluctuation increases in both cases. The main contribution to the overall enzyme fluctuation is given by the loop 164-173. The fluctuation of this loop is spread over more degrees of freedom when PLTP interacts with the S-side. This dynamical factor can enhance the preference of PLTP for the S subsites of MMP-8. MD simulations also show that ligation of PLTP in the S subsites is further favoured by better zinc chelation, a cation-p interaction at the S3 subsite and unstrained binding conformations. The role of the S3, S3' and S1' subsites in determining the inhibitor binding is dicussed.

The active site of thermolysin-like proteases (TLPs) is located at the bottom of a cleft between the N- and C-terminal domains. Crystallographic studies have shown that the active-site cleft is more closed in Ligand-binding TLPs than in Ligand-free TLPs. Accordingly, it has been proposed that TLPs undergo a hinge-bending motion during catalysis resulting in "closure" and "opening" of the active-site cleft. Two hinge regions have been proposed. One is located around a conserved glycine 78; the second involves residues 135 and 136. The importance of conserved glycine residues in these hinge regions was studied experimentally by analyzing the effects of Gly --> Ala mutations on catalytic activity. Eight such mutations were made in the TLP of Bacillus stearothermophilus (TLP-ste) and their effects on activity toward casein and various peptide substrates were determined. Only the Gly78Ala, Gly136Ala, and Gly135Ala + Gly136Ala mutants decreased catalytic activity significantly. These mutants displayed a reduction in k(cat)/K-m for 3-(2-furylacryloyl)-L-glycyl-L-leucine amide of 73%, 62%, and 96%, respectively. Comparisons of effects on k(cat)/K-m for various substrates with effects on the K-i for phosphoramidon suggested that the mutation at position 78 primarily had an effect on substrate binding, whereas the mutations at positions 135 and 136 primarily influence k(cat). The apparent importance of conserved glycine residues in proposed hinge-bending regions for TLP activity supports the idea that hinge-bending is an essential part of catalysis.

We address the question how well proteins can be modelled on the basis of NMR data, when these data are incorporated into the protein model using distance restraints in a molecular dynamics simulation. We found, using HPr as a model protein, that distance restraining freezes the essential motion of proteins, as defined by Amadei et al. [Amadei, A., Linssen, A.B.M. and Berendsen, H.J.C. (1993) Protein Struct. Funct. Genet., 17, 412-425]. We discuss how modelling protocols can be improved in order to solve this problem.

The cellular retinol-binding protein (CRBP) is an intracellular retinol carrier protein belonging to a family of hydrophobic ligand-binding proteins, It transports retinol to specific locations in the cell where, for instance, it is esterified for storage, Recently solved crystallographic structures of CRBP homologues with and without bound ligand do not provide evidence for a ligand-induced conformational change, However, it has been shown that there is a difference in binding of holo-CRBP and apo-CRBP to lecithin-retinol acyltransferase, Moreover, proteolysis of holo-CRBP and apo-CRBP yields different products, indicating a difference in structure or dynamics between the two forms, Here, we present the results of molecular dynamics simulations of holo-CRBP and apo-CRBP, The simulations show a significant difference in conformation, in agreement with experimental results, The essential dynamics method was used to study differences in dynamics between the apo and hole forms of CRBP, and showed inhibition of essential motions upon ligand binding, It also revealed large correlated motions of retinol with regions of the protein, pointing to a possible retinol entry/exit site.