A coarse-grained Langevin molecular dynamics approach to de novo protein structure prediction

De novo prediction of protein structures, the prediction of structures from amino acid sequences which are not similar to those of hitherto resolved structures, has been one of the major challenges in molecular biophysics. In this paper, we develop a new method of de novo prediction, which combines the fragment assembly method and the simulation of physical folding process: structures which have consistently assembled fragments are dynamically searched by Langevin molecular dynamics of conformational change. The benchmarking test shows that the prediction is improved when the candidate structures are cross-checked by an empirically derived score function.     

Conformational transition of an α-helix studied by molecular dynamics

Molecular dynamics simulations were performed on a 20-residue polyalanine helix and a spontaneous transition from a kinked to a straight conformation was observed. The kinetics of the transition was analyzed within the framework of the Kramers model for chemical reactions and within a random walk model. The Kramers model which is based on diffusion along a one-dimensional reaction pathway and the crossing of an energy barrier was found to be inadequate. Instead, a random walk model based on diffusion in the high-dimensional phase space of the system was found to be compatible with the data. The high dimensionality of the phase space permits the system to circumvent high energy barriers and diffuse rapidly at about constant energy, but decelerates the reaction since in the labyrinth of pathways the transition state is reached rarely.     

Residual motion of hemoglobin-bound spin labels and protein dynamics: viscosity dependence of the rotational correlation times

The residual motion of spin labels bound to cysteine 93 and to lysines of methemoglobin has been studied by electron paramagnetic resonance spectroscopy. To separate the influences of the solvent and the protein environment of the label fluctuations, the correlation times, , were analyzed as a function of temperature for fixed solvent viscosities, . Results show that over a wide range of viscosity the dependence of on may be empirically described by a power law ^k . The exponent k depends strongly on the location of the label on the protein surface. If one regards the spin labels as artificial amino acid side chains, characteristic values of correlation times and amplitudes of the rotational motion at the surface can be given. For =1 cP and T=297 K the correlation time of the labels bound to lysines is found to be =9 · 10^–10 s and the rotational diffusion is nearly isotropic. The spin label bound to cysteine 93 occupies a protein pocket, its rotational motion is therefore restricted. The correlation time of the label motion within a limited motion cone of semi angle =30° ± 3° is found to be =1.3 · 10^–9 s for =1 cP and T=297 K.     

Proteome dynamics in complex organisms: using stable isotopes to monitor individual protein turnover rates

The complete definition of changes in a proteome requires information about dynamics and specifically the rate at which the individual proteins are turned over intracellularly. Whilst this can be achieved in single-cell culture using stable isotope precursors, it is more challenging to develop methods for intact animals. In this study, we show how dietary administration of stable isotope-labelled amino acids can obtain information on the relative rates of synthesis and degradation of individual proteins in a proteome. The pattern of stable isotope-labelling in tryptic peptides can be deconstructed to yield a highly reliable measure of the isotope abundance of the precursor pool, a parameter that is often difficult to acquire. We demonstrate this approach using chickens fed a semisynthetic diet containing [(2)H(8)]valine at a calculated relative isotope abundance (RIA) of 0.5. When the labelling pattern of gel-resolved muscle proteins was analyzed, the intracellular precursor isotope abundance was 0.35, consistent with dilution of the amino acid precursor pool with unlabelled amino acids derived from degradation of pre-existing proteins. However, the RIA was stable over an extended labelling window, and permitted calculation of the rates of synthesis and degradation of individual proteins isolated by gel electrophoresis. For the first time, it is feasible to contemplate the analysis of turnover of individual proteins in intact animals.     

Influence of protein dynamics on the metal-sites of ovotransferrin

Using the perturbed angular correlations (PAC) technique, the formation of hafnium-ovotransferrin complexes has been studied. Two binding configurations at each of the two specific binding-sites of the protein have been observed. They are characterized by well-defined electric quadrupole frequencies. Information about the dynamics of the protein was derived from temperature dependent measurements of the relaxation constant. The well-resolved spectra taken with fast BaF₂-detectors allow a precise determination of the relaxation behaviour of the protein. The results are compared with the predictions from a hydrodynamic model for the reorientation of macromolecules. Thus the hydrodynamic volume of ovotransferrin and its N-terminal half-molecule were determined. The ovotransferrin volume is in agreement with a value derived for human serum transferrin from small angle neutron scattering. From experiments with immobilized protein material there is evidence for internal protein dynamics which is probed by the Hf-ion bound to the specific metal-sites.Abbreviations PAC perturbed angular correlations technique - TF serum transferrin - LF lactoferrin - OTF ovotransferrin - OTF/2N N-terminal half-molecule of ovotransferrin - NQR nuclear quadrupole resonance - EFG electric field gradient - NQI nuclear quadrupole interaction - NTA nitrilotriacetate - MES 2-(N-morpholino)ethanesulphonic acid - HEPPS N-2-hydroxyethyl-piperazine-N-3-propanesulphonic acid - TRIS tris(hydroxymethyl)aminomethane Correspondence to: F. J. Schwab.     

Integrating protein structural dynamics and evolutionary analysis with Bio3D

Background

Popular bioinformatics approaches for studying protein functional dynamics include comparisons of crystallographic structures, molecular dynamics simulations and normal mode analysis. However, determining how observed displacements and predicted motions from these traditionally separate analyses relate to each other, as well as to the evolution of sequence, structure and function within large protein families, remains a considerable challenge. This is in part due to the general lack of tools that integrate information of molecular structure, dynamics and evolution.

Results

Here, we describe the integration of new methodologies for evolutionary sequence, structure and simulation analysis into the Bio3D package. This major update includes unique high-throughput normal mode analysis for examining and contrasting the dynamics of related proteins with non-identical sequences and structures, as well as new methods for quantifying dynamical couplings and their residue-wise dissection from correlation network analysis. These new methodologies are integrated with major biomolecular databases as well as established methods for evolutionary sequence and comparative structural analysis. New functionality for directly comparing results derived from normal modes, molecular dynamics and principal component analysis of heterogeneous experimental structure distributions is also included. We demonstrate these integrated capabilities with example applications to dihydrofolate reductase and heterotrimeric G-protein families along with a discussion of the mechanistic insight provided in each case.

Conclusions

The integration of structural dynamics and evolutionary analysis in Bio3D enables researchers to go beyond a prediction of single protein dynamics to investigate dynamical features across large protein families. The Bio3D package is distributed with full source code and extensive documentation as a platform independent R package under a GPL2 license from http://thegrantlab.org/bio3d/.

Electronic supplementary material

The online version of this article (doi:10.1186/s12859-014-0399-6) contains supplementary material, which is available to authorized users.     

Exploring the energy landscape of protein folding using replica-exchange and conventional molecular dynamics simulations

Two independent replica-exchange molecular dynamics (REMD) simulations with an explicit water model were performed of the Trp-cage mini-protein. In the first REMD simulation, the replicas started from the native conformation, while in the second they started from a nonnative conformation. Initially, the first simulation yielded results qualitatively similar to those of two previously published REMD simulations: the protein appeared to be over-stabilized, with the predicted melting temperature 50-150K higher than the experimental value of 315K. However, as the first REMD simulation progressed, the protein unfolded at all temperatures. In our second REMD simulation, which starts from a nonnative conformation, there was no evidence of significant folding. Transitions from the unfolded to the folded state did not occur on the timescale of these simulations, despite the expected improvement in sampling of REMD over conventional molecular dynamics (MD) simulations. The combined 1.42 micros of simulation time was insufficient for REMD simulations with different starting structures to converge. Conventional MD simulations at a range of temperatures were also performed. In contrast to REMD, the conventional MD simulations provide an estimate of Tm in good agreement with experiment. Furthermore, the conventional MD is a fraction of the cost of REMD and continuous, realistic pathways of the unfolding process at atomic resolution are obtained.     

Water dynamics clue to key residues in protein folding

A computational method independent of experimental protein structure information is proposed to recognize key residues in protein folding, from the study of hydration water dynamics. Based on all-atom molecular dynamics simulation, two key residues are recognized with distinct water dynamical behavior in a folding process of the Trp-cage protein. The identified key residues are shown to play an essential role in both 3D structure and hydrophobic-induced collapse. With observations on hydration water dynamics around key residues, a dynamical pathway of folding can be interpreted.     

Kinetic analysis of biphasic protein modification reactions cooperative effects

A mathematical treatment of protein modification reactions is presented, and it is shown thai in these cases protein modification is described by a summation of exponential functions of reaction time, the number of exponentials being equal to the number of modified protein species. It is shown that in cases of protein modification cooperativity, there is a strict dependence of the coefficients of the multiexponential modification equation on the constants of the same equation. The conditions necessary for a reduction of a multiexponential protein modification equation to one of a summation of two exponentials only are examined. The possible formulae for the coefficients of a two-exponential-summation equation, used to describe the modification of protein models with two, three or four modifiable residues (as well as some aspects of models with five and six modifiable residues) per protein molecule are derived. It is seen that the number of such coefficients is severely limited. The most frequently obtained formula for the lower stoichiomelric coefficient of a 'wo-exponential-summation equation is Ak_a/(k_a-k_b). where k_b and k_b are the constants of the two exponentials of the equation, and A is a constant. The value most frequently arrived at for A is (n?1)/n, where n is the number of modifiable residues per protein molecule, while values such as 1/n, or a/n (where a is an integer, and also where a < n) are also possible. In most of the cooperative protein modification models worked out, k_a is identical with k_n, viz., k_a is identical with the rate constant for the first stoichiometric protein modification.     

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     

Role of methyl groups in dynamics and evolution of biomolecules

Recent studies have discovered strong differences between the dynamics of nucleic acids (RNA and DNA) and proteins, especially at low hydration and low temperatures. This difference is caused primarily by dynamics of methyl groups that are abundant in proteins, but are absent or very rare in RNA and DNA. In this paper, we present a hypothesis regarding the role of methyl groups as intrinsic plasticizers in proteins and their evolutionary selection to facilitate protein dynamics and activity. We demonstrate the profound effect methyl groups have on protein dynamics relative to nucleic acid dynamics, and note the apparent correlation of methyl group content in protein classes and their need for molecular flexibility. Moreover, we note the fastest methyl groups of some enzymes appear around dynamical centers such as hinges or active sites. Methyl groups are also of tremendous importance from a hydrophobicity/folding/entropy perspective. These significant roles, however, complement our hypothesis rather than preclude the recognition of methyl groups in the dynamics and evolution of biomolecules.     

Ion-specific diffusion rates through transmembrane protein channels. A molecular dynamics study

The migration of different alkali metal cations through a transmembrane model channel is simulated by means of the molecular dynamics technique. The parameters of the model are chosen in close relation to the gramicidin A channel. Coulomb- and van der Waals-type potentials between the ions and flexible carbonyl groups of the pore-forming molecule are used to describe the ion channel interaction. The diffusion properties of the ions are obtained from three-dimensional trajectory calculations. The diffusion rates for the different ions Li+, Na+, K+ and Rb+ are affected not only by the mass of the particles but also very strongly by their size. The latter effect is more pronounced for rigid channels, i.e., for binding vibrational frequencies of the CO groups with v greater than 400 cm-1. In this range the selectivity sequence for the diffusion rates is the inverse of that expected from normal rate theory but agrees with that found in experiments for gramicidin A.     

Viscosity B-coefficients and standard partial molar volumes of amino acids, and their roles in interpreting the protein (enzyme) stabilization

This review systematically surveys the viscosity B-coefficients and standard partial molar volumes of amino acids at various temperatures as these data are quite important for interpreting the hydration and other properties of peptides and proteins. The effect of organic solutes and various ions on the viscometric and volumetric properties of amino acids has also been discussed in terms of their kosmotropic ('structure-making') effects on the hydration of amino acids. The comparison of these effects on the amino acid hydration enables us to have a better understanding of the influence of organic solute and salt on the protein stabilization. In addition, the viscometric and volumetric behaviors of amino acid ions (cations and anions) are also summarized because these ions have recently been incorporated as part of novel ionic liquids, which have wide applications in biocatalysis and protein stabilization.     

Kinesins and protein kinases: key players in the regulation of microtubule dynamics and organization

Microtubule dynamics is controlled and amplified in vivo by complex sets of regulators. Among these regulatory proteins, molecular motors from the kinesin superfamily are taking an increasing importance. Here we review how microtubule disassembly or assembly into interphase microtubules, mitotic spindle or cilia may involve kinesins and how protein kinases may participate in these kinesin-dependent regulations.     

Yeast proteomics and protein microarrays

Our understanding of biological processes as well as human diseases has improved greatly thanks to studies on model organisms such as yeast. The power of scientific approaches with yeast lies in its relatively simple genome, its facile classical and molecular genetics, as well as the evolutionary conservation of many basic biological mechanisms. However, even in this simple model organism, systems biology studies, especially proteomic studies had been an intimidating task. During the past decade, powerful high-throughput technologies in proteomic research have been developed for yeast including protein microarray technology. The protein microarray technology allows the interrogation of protein–protein, protein–DNA, protein–small molecule interaction networks as well as post-translational modification networks in a large-scale, high-throughput manner. With this technology, many groundbreaking findings have been established in studies with the budding yeast Saccharomyces cerevisiae, most of which could have been unachievable with traditional approaches. Discovery of these networks has profound impact on explicating biological processes with a proteomic point of view, which may lead to a better understanding of normal biological phenomena as well as various human diseases.     

Semi-automated microplate monitoring of protein polymerization and aggregation

Static light scattering (SLS) is a commonly used technique for monitoring dynamics of high molecular weight protein complexes such as protein oligomers or aggregates. However, traditional methods are limited to testing a single condition and typically require large amounts of protein and specialized equipment. We show that a standard microplate reader can be used to characterize the molecular dynamics of different types of protein complexes, with the multiple advantages of microscale experimental volumes, semi-automated protocols and highly parallel processing.