ProMode: a database of normal mode analyses on protein molecules with a full-atom model

MOTIVATION: Although information from protein dynamics simulation is important to understand principles of architecture of a protein structure and its function, simulations such as molecular dynamics and Monte Carlo are very CPU-intensive. Although the ability of normal mode analysis (NMA) is limited because of the need for a harmonic approximation on which NMA is based, NMA is adequate to carry out routine analyses on many proteins to compute aspects of the collective motions essential to protein dynamics and function. Furthermore, it is hoped that realistic animations of the protein dynamics can be observed easily without expensive software and hardware, and that the dynamic properties for various proteins can be compared with each other. RESULTS: ProMode, a database collecting NMA results on protein molecules, was constructed. The NMA calculations are performed with a full-atom model, by using dihedral angles as independent variables, faster and more efficiently than the calculations using Cartesian coordinates. In ProMode, an animation of the normal mode vibration is played with a free plug-in, Chime (MDL Information Systems, Inc.). With the full-atom model, the realistic three-dimensional motions at an atomic level are displayed with Chime. The dynamic domains and their mutual screw motions defined from the NMA results are also displayed. Properties for each normal mode vibration and their time averages, e.g. fluctuations of atom positions, fluctuations of dihedral angles and correlations between the atomic motions, are also presented graphically for characterizing the collective motions in more detail. AVAILABILITY: http://promode.socs.waseda.ac.jp     

600 ps Molecular dynamics reveals stable substructures and flexible hinge points in cAMP dependent protein kinase

Molecular dynamics simulations of the catalytic subunit of cAMP dependent protein kinase (cAPK) have been performed in an aqueous environment. The relations among the protein hydrogen‐bonding network, secondary structural elements, and the internal motions of rigid domains were examined. The values of fluctuations of protein dihedral angles during dynamics show quite distinct maxima in the regions of loops and minima in the regions of α‐helices and β‐strands. Analyses of conformation snapshots throughout the run show stable subdomains and indicate that these rigid domains are constrained during the dynamics by a stable network of hydrogen bonds. The most stable subdomain during the dynamics was in the small lobe including part of the carboxy‐terminal tail. The most significant flexible region was the highly conserved glycine‐rich loop between β strands 1 and 2 in the small lobe. Many of the main chain dihedral angle changes measured in a comparison of the crystallographic structures of “open” and “closed” conformations of cAPK correspond to the highly flexible residues found during dynamics. © 1999 John Wiley & Sons, Inc. Biopoly 50: 513–524, 1999     

Studying backbone torsional dynamics of intrinsically disordered proteins using fluorescence depolarization kinetics

Intrinsically disordered proteins (IDPs) do not autonomously adopt a stable unique 3D structure and exist as an ensemble of rapidly interconverting structures. They are characterized by significant conformational plasticity and are associated with several biological functions and dysfunctions. The rapid conformational fluctuation is governed by the backbone segmental dynamics arising due to the dihedral angle fluctuation on the Ramachandran ?–ψ conformational space. We discovered that the intrinsic backbone torsional mobility can be monitored by a sensitive fluorescence readout, namely fluorescence depolarization kinetics, of tryptophan in an archetypal IDP such as α-synuclein. This methodology allows us to map the site-specific torsional mobility in the dihedral space within picosecond-nanosecond time range at a low protein concentration under the native condition. The characteristic timescale of ~?1.4 ns, independent of residue position, represents collective torsional dynamics of dihedral angles (? and ψ) of several residues from tryptophan and is independent of overall global tumbling of the protein. We believe that fluorescence depolarization kinetics methodology will find broad application to study both short-range and long-range correlated motions, internal friction, binding-induced folding, disorder-to-order transition, misfolding and aggregation of IDPs.     

VADAR: a web server for quantitative evaluation of protein structure quality

VADAR (Volume Area Dihedral Angle Reporter) is a comprehensive web server for quantitative protein structure evaluation. It accepts Protein Data Bank (PDB) formatted files or PDB accession numbers as input and calculates, identifies, graphs, reports and/or evaluates a large number (>30) of key structural parameters both for individual residues and for the entire protein. These include excluded volume, accessible surface area, backbone and side chain dihedral angles, secondary structure, hydrogen bonding partners, hydrogen bond energies, steric quality, solvation free energy as well as local and overall fold quality. These derived parameters can be used to rapidly identify both general and residue-specific problems within newly determined protein structures. The VADAR web server is freely accessible at http://redpoll.pharmacy.ualberta.ca/vadar.     

Dynamics of transfer RNAs analyzed by normal mode calculation.

Normal mode calculation is applied to tRNAPhe and tRNAAsp, and their structural and vibrational aspects are analyzed. Dihedral angles along the phosphate-ribose backbone (alpha, beta, gamma, epsilon, zeta) and dihedral angles of glycosyl bonds (chi) are selected as movable parameters. The calculated displacement of each atom agrees with experimental data. In modes with frequencies higher than 130 cm-1, the motions are localized around each stem and the elbow region of the L-shape. On the other hand, collective motions such as bending or twisting of arms are seen in modes with lower frequencies. Hinge axes and bend angles are calculated without prior knowledge. Movements in modes with very low frequencies are combinations of hinge bending motions with various hinge axes and bend angles. The thermal fluctuations of dihedral angles well reflect the structural characters of transfer RNAs. There are some dihedral angles of nucleotides located around the elbow region of L-shape, which fluctuate about five to six times more than the average value. Nucleotides in the position seem to be influential in the dynamics of the entire structure. The normal mode calculation seems to provide much information for the study of conformational changes of transfer RNAs induced by aminoacyl-tRNA synthetase or codon during molecular recognition.     

Mass-weighted molecular dynamics simulation and conformational analysis of polypeptide.

Atomic motions in protein molecules have been studied by molecular dynamics (MD) simulations; dynamics simulation methods have also been employed in conformational studies of polypeptide molecules. It was found that when atomic masses are weighted, the molecular dynamics method can significantly increase the sampling of dihedral conformation space in such studies, compared to a conventional MD simulation of the same total simulation time length. Herein the theoretical study of molecular conformation sampling by the molecular dynamics-based simulation method in which atomic masses are weighted is reported in detail; moreover, a numerical scheme for analyzing the extensive conformational sampling in the simulation of a tetrapeptide amide molecule is presented. From numerical analyses of the mass-weighted molecular dynamics trajectories of backbone dihedral angles, low-resolution structures covering the entire backbone dihedral conformation space of the molecule were determined, and the distribution of rotationally stable conformations in this space were analyzed quantitatively. The theoretical analyses based on the computer simulation and numerical analytical methods suggest that distinctive regimes in the conformational space of the peptide molecule can be identified.     

WEBnm@ v2.0: Web server and services for comparing protein flexibility

Background

Normal mode analysis (NMA) using elastic network models is a reliable and cost-effective computational method to characterise protein flexibility and by extension, their dynamics. Further insight into the dynamics–function relationship can be gained by comparing protein motions between protein homologs and functional classifications. This can be achieved by comparing normal modes obtained from sets of evolutionary related proteins.

Results

We have developed an automated tool for comparative NMA of a set of pre-aligned protein structures. The user can submit a sequence alignment in the FASTA format and the corresponding coordinate files in the Protein Data Bank (PDB) format. The computed normalised squared atomic fluctuations and atomic deformation energies of the submitted structures can be easily compared on graphs provided by the web user interface. The web server provides pairwise comparison of the dynamics of all proteins included in the submitted set using two measures: the Root Mean Squared Inner Product and the Bhattacharyya Coefficient. The Comparative Analysis has been implemented on our web server for NMA, WEBnm@, which also provides recently upgraded functionality for NMA of single protein structures. This includes new visualisations of protein motion, visualisation of inter-residue correlations and the analysis of conformational change using the overlap analysis. In addition, programmatic access to WEBnm@ is now available through a SOAP-based web service. Webnm@ is available at http://apps.cbu.uib.no/webnma.

Conclusion

WEBnm@ v2.0 is an online tool offering unique capability for comparative NMA on multiple protein structures. Along with a convenient web interface, powerful computing resources, and several methods for mode analyses, WEBnm@ facilitates the assessment of protein flexibility within protein families and superfamilies. These analyses can give a good view of how the structures move and how the flexibility is conserved over the different structures.

Electronic supplementary material

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

Ligand-induced conformational change of a protein reproduced by a linear combination of displacement vectors obtained from normal mode analysis

The conformational change of a protein upon ligand binding was examined by normal mode analysis (NMA) based on an elastic-network model (ENM) for a full-atom system using dihedral angles as independent variables. Specifically, we investigated the extent to which conformational change vectors of atoms from an apo form to a holo form of a protein can be represented by a linear combination of the displacement vectors of atoms in the apo form calculated for the lowest-frequency m normal modes (m=1, 2,…, 20). In this analysis, the latter vectors were best fitted to the former ones by the least-squares method. Twenty-two paired proteins in the holo and apo forms, including three dimer pairs, were examined. The results showed that, in most cases, the conformational change vectors were reproduced well by a linear combination of the displacement vectors of a small number of low-frequency normal modes. The conformational change around an active site was reproduced as well as the entire conformational change, except for some proteins that only undergo significant conformational changes around active sites. The weighting factors for 20 normal modes optimized by the least-squares fitting characterize the conformational changes upon ligand binding for these proteins. The conformational changes sampled around the apo form of a protein by the linear combination of the displacement vectors obtained by ENM-based NMA may help solve the flexible-docking problem of a protein with another molecule because the results presented herein suggest that they have a relatively high probability of being involved in an actual conformational change.     

Structure and internal mobility of proteins: a molecular dynamics study of hen egg white lysozyme

The structure and internal motions of the protein hen egg white lysozyme are studied by analysis of simulation and experimental data. A molecular dynamics simulation and an energy minimization of the protein in vacuum have been made and the results compared with high-resolution structures and temperature factors of hen egg white lysozyme in two different crystal forms and of the homologous protein human lysozyme. The structures obtained from molecular dynamics and energy minimization have root-mean-square deviations for backbone atoms of 2.3 Å and 1.1–1.3 Å, respectively, relative to the crystal structures; the different crystal structures have root-mean-square deviations of 0.73–0.81 Å for the backbone atoms. In comparing the backbone dihedral angles, the difference between the dynamics and the crystal structure on which it is based is the same as that between any two crystal structures. The internal fluctuations of atomic positions calculated from the molecular dynamics trajectory agree well with the temperature factors from the three structures. Simulation and crystal results both show that there are large motions for residues involved in exposed turns of the backbone chain, relatively smaller motions for residues involved in the middle of helices or β-sheet structures, and relatively small motions of residues near disulfide bridges. Also, both the simulation and crystal data show that side-chain atoms have larger fluctuations than main-chain atoms. Moreover, the regions that have large deviations among the x-ray crystal structures, which indicates flexibility, are found to have large fluctuations in the simulation.     

Equilibrium transitions between side‐chain conformations in leucine and isoleucine

Despite recent improvements in computational methods for protein design, we still lack a quantitative, predictive understanding of the intrinsic probabilities for amino acids to adopt particular side‐chain conformations. Surprisingly, this question has remained unsettled for many years, in part because of inconsistent results from different experimental approaches. To explicitly determine the relative populations of different side‐chain dihedral angles, we performed all‐atom hard‐sphere Langevin Dynamics simulations of leucine (Leu) and isoleucine (Ile) dipeptide mimetics with stereo‐chemical constraints and repulsive‐only steric interactions between non‐bonded atoms. We determine the relative populations of the different χ₁ and χ₂ dihedral angle combinations as a function of the backbone dihedral angles ? and ψ. We also propose, and test, a mechanism for inter‐conversion between the different side‐chain conformations. Specifically, we discover that some of the transitions between side‐chain dihedral angle combinations are very frequent, whereas others are orders of magnitude less frequent, because they require rare coordinated motions to avoid steric clashes. For example, to transition between different values of χ₂, the Leu side‐chain bond angles κ₁ and κ₂ must increase, whereas to transition in χ₁, the Ile bond angles λ₁ and λ₂ must increase. These results emphasize the importance of computational approaches in stimulating further experimental studies of the conformations of side‐chains in proteins. Moreover, our studies emphasize the power of simple steric models to inform our understanding of protein structure, dynamics, and design. Proteins 2015; 83:1488–1499. © 2015 Wiley Periodicals, Inc.     

A new spin on protein dynamics

Site-directed spin labeling is a general method for investigating structure and conformational switching in soluble and membrane proteins. It will also be an important tool for exploring protein backbone dynamics. A semi-empirical analysis of nitroxide sidechain dynamics in spin-labeled proteins reveals contributions from fluctuations in backbone dihedral angles and rigid-body (collective) motions of alpha helices. Quantitative analysis of sidechain dynamics is sometimes possible, and contributions from backbone modes can be expressed in terms of relative order parameters and rates. Dynamic sequences identified by site-directed spin labeling correlate with functional domains, and so nitroxide scanning could provide an efficient strategy for identifying such domains in high-molecular weight proteins, supramolecular complexes and membrane proteins.     

Comparison of normal mode analyses on a small globular protein in dihedral angle space and Cartesian coordinate space

Normal mode analyses on the protein, bovine pancreatic trypsin inhibitor, in dihedral angle space and Cartesian coordinate space are compared. In Cartesian coordinate space it is found that modes of frequencies lower than 30 cm(-1) contribute 80% of the total mean-square fluctuation and are represented almost completely by motions in the dihedral angles. Bond angle and length fluctuations dominate in modes above 200 cm(-1), but contribute less than 2% to the total mean-square fluctuation. In the low-frequency modes a good correspondence between patterns of atomic displacements was found, but on average the root-mean-square fluctuations of the Cartesian coordinate modes are 13% greater than their dihedral angle counterparts. The main effect of fluctuations in the bond angles and lengths, therefore, is to allow the dihedral angles to become more flexible. As the important subspaces determined from the two methods overlap considerably, dihedral angle space analysis can be applied to proteins too large for Cartesian coordinate space analysis.     

Statistical and conformational analysis of the electron density of protein side chains

Protein side chains make most of the specific contacts between proteins and other molecules, and their conformational properties have been studied for many years. These properties have been analyzed primarily in the form of rotamer libraries, which cluster the observed conformations into groups and provide frequencies and average dihedral angles for these groups. In recent years, these libraries have improved with higher resolution structures and using various criteria such as high thermal factors to eliminate side chains that may be misplaced within the crystallographic model coordinates. Many of these side chains have highly non-rotameric dihedral angles. The origin of side chains with high B-factors and/or with non-rotameric dihedral angles is of interest in the determination of protein structures and in assessing the prediction of side chain conformations. In this paper, using a statistical analysis of the electron density of a large set of proteins, it is shown that: (1) most non-rotameric side chains have low electron density compared to rotameric side chains; (2) up to 15% of chi1 non-rotameric side chains in PDB models can clearly be fit to density at a single rotameric conformation and in some cases multiple rotameric conformations; (3) a further 47% of non-rotameric side chains have highly dispersed electron density, indicating potentially interconverting rotameric conformations; (4) the entropy of these side chains is close to that of side chains annotated as having more than one chi(1) rotamer in the crystallographic model; (5) many rotameric side chains with high entropy clearly show multiple conformations that are not annotated in the crystallographic model. These results indicate that modeling of side chains alternating between rotamers in the electron density is important and needs further improvement, both in structure determination and in structure prediction.     

High-throughput modeling and analysis of protein structural dynamics

Protein function is a dynamic property closely related to the conformational mechanisms of protein structure in its physiological environment. To understand and control the function of target proteins, it becomes increasingly important to develop methods and tools for predicting collective motions at the molecular level. In this article, we review computational methods for predicting conformational dynamics and discuss software tools for data analysis. In particular, we discuss a high-throughput, web-based system called iGNM for protein structural dynamics. iGNM contains a database of protein motions for more than 20 000 PDB structures and supports online calculations for newly deposited PDB structures or user-modified structures. iGNM allows dynamics analysis of protein structures ranging from enzymes to large complexes and assemblies, and enables the exploration of protein sequence-structure-dynamics-function relations.     

Coarse-grained modeling of the actin filament derived from atomistic-scale simulations

A coarse-grained (CG) procedure that incorporates the information obtained from all-atom molecular dynamics (MD) simulations is presented and applied to actin filaments (F-actin). This procedure matches the averaged values and fluctuations of the effective internal coordinates that are used to define a CG model to the values extracted from atomistic MD simulations. The fluctuations of effective internal coordinates in a CG model are computed via normal-mode analysis (NMA), and the computed fluctuations are matched with the atomistic MD results in a self-consistent manner. Each actin monomer (G-actin) is coarse-grained into four sites, and each site corresponds to one of the subdomains of G-actin. The potential energy of a CG G-actin contains three bonds, two angles, and one dihedral angle; effective harmonic bonds are used to describe the intermonomer interactions in a CG F-actin. The persistence length of a CG F-actin was found to be sensitive to the cut-off distance of assigning intermonomer bonds. Effective harmonic bonds for a monomer with its third nearest neighboring monomers are found to be necessary to reproduce the values of persistence length obtained from all-atom MD simulations. Compared to the elastic network model, incorporating the information of internal coordinate fluctuations enhances the accuracy and robustness for a CG model to describe the shapes of low-frequency vibrational modes. Combining the fluctuation-matching CG procedure and NMA, the achievable time- and length scales of modeling actin filaments can be greatly enhanced. In particular, a method is described to compute the force-extension curve using the CG model developed in this work and NMA. It was found that F-actin is easily buckled under compressive deformation, and a writhing mode is developed as a result. In addition to the bending and twisting modes, this novel writhing mode of F-actin could also play important roles in the interactions of F-actin with actin-binding proteins and in the force-generation process via polymerization.     

Mollib: a molecular and NMR data analysis software

Mollib is a software framework for the analysis of molecular structures, properties and data with an emphasis on data collected by NMR. It uses an open source model and a plugin framework to promote community-driven development of new and enhanced features. Mollib includes tools for the automatic retrieval and caching of protein databank (PDB) structures, the hydrogenation of biomolecules, the analysis of backbone dihedral angles and hydrogen bonds, and the fitting of residual dipolar coupling (RDC) and residual anisotropic chemical shift (RACS) data. In this article, we release version 1.0 of mollib and demonstrate its application to common molecular and NMR data analyses.     

Structure and dynamics of an amphiphilic peptide in a lipid bilayer: a molecular dynamics study

A molecular dynamics simulation of a simple model membrane system composed of a single amphiphilic helical peptide (ace-K2GL16K2A-amide) in a fully hydrated 1,2-dimyristoyl-sn-glycero-3-phosphocholine bilayer was performed for a total of 1060 ps. The secondary structure of the peptide and its stability were described in terms of average dihedral angles, phi and psi, and the C alpha torsion angles formed by backbone atoms; by the average translation per residue along the helix axis; and by the intramolecular peptide hydrogen bonds. The results indicated that residues 6 through 15 remain in a stable right-handed alpha-helical conformation, whereas both termini exhibit substantial fluctuations. A change in the backbone dihedral angles for residues 16 and 17 is accompanied by the loss of two intramolecular hydrogen bonds, leading to a local but long-lived disruption of the helix. The dynamics of the peptide was characterized in terms of local and global helix motions. The local motions of the N-H bond angles were described in terms of the autocorrelation functions of P2[cos thetaNH(t, t + tau)] and reflected the different degrees of local peptide order as well as a variation in time scale for local motions. The chi1 and chi2 dihedral angles of the leucine side chains underwent frequent transitions between potential minima. No connection between the side-chain positions and their mobility was observed, however. In contrast, the lysine side chains displayed little mobility during the simulation. The global peptide motions were characterized by the tilting and bending motions of the helix. Although the peptide was initially aligned parallel to the bilayer normal, during the simulation it was observed to tilt away from the normal, reaching an angle of approximately 25 degrees by the end of the simulation. In addition, a slight bend of the helix was detected. Finally, the solvation of the peptide backbone and side-chain atoms was also investigated.     

Vibrational approach to the dynamics of an α-helix

The dynamics of a finite α-helix have been studied in the harmonic approximation by a vibrational analysis of the atomic motions about their equilibrium positions. The system were represented by an empirical potential energy function, and all degrees of freedom (bond lengths, bond angles, and torsional angles) were allowed to vary. The complete results were compared with a more restrictive model in which the peptide dihedral angle was kept rigid; also, a model potential excluding hydrogen bonds was examined. Thermal fluctuations in the backbone dihedral angles ? and ψ are 12° to 15°. The fluctuations of adjacent dihedral angles are highly correlated, and the correlation pattern is affected by the flexibility of the peptide dihedral angle. Time-dependent autocorrelations in the motion of ? and ψ appear to decay due to dephasing in less than 1 psec, while the motions of the carbonyl oxygen and amide hydrogens out of the peptide plane are more harmonic. Length fluctuations have been evaluated and exhibit a strong end effect; the calculated elastic modulus is in agreement with other values. Rigid and adiabatic total energy surfaces corresponding to dihedral angle rotations in the middle of the helix have been obtained and compared with the quadratic approximation to those surfaces. The magnitudes and correlations between the fluctuations obtained by averaging over the adiabatic energy surface most closely resemble the vibrational results. Of particular interest is the fact that hydrogen bonds play a relatively small role in the local dihedral angle fluctuations, though the hydrogen bonds are important in the energy of overall length changes.     

Relationships between amino acid sequence and backbone torsion angle preferences

Statistical averages and correlations for backbone torsion angles of chymotrypsin inhibitor 2 are calculated by using the Rotational Isomeric States model of chain statistics. Statistical weights of torsional states of phipsi pairs, needed for the statistics of the full chain, are obtained in two different ways: 1) by using knowledge-based pairwise dependent phipsi energy maps from Protein Data Bank (PDB) and 2) by collecting torsion angle data from a large number of random coil configurations of an all-atom protein model with volume exclusion. Results obtained by using PDB data show strong correlations between adjacent torsion angle pairs belonging to both the same and different residues. These correlations favor the choice of the native-state torsion angles, and they are strongly context dependent, determined by the specific amino acid sequence of the protein. Excluded volume or steric clashes, only, do not introduce context-dependent phipsi correlations into the chain that would affect the choice of native-state torsional angles.     

Conformational motion correlations in the formation of polypeptide secondary structure in a viscous medium

The Langevin dynamics method and statistical correlation analysis were used to study the α-helical structure folding dynamics of the (Ala)₅₀, (AlaGly)₂₅, and (AlaGly)₇₅ polypeptides depending on the viscosity of the medium. Friction forces that arise when the effective viscosity of the medium is similar to the viscosity of water were found to result in strong correlations between the backbone torsion angles. The polypeptides under study folded mainly to produce α-helical structures. A structure of two contacting α-helices that were approximately equal in length and had a loop between them was observed for a longer chain of 150 residues. A method to visualize the correlation matrix of the dihedral angles of a polypeptide chain was developed for analyzing the effects of the dynamic correlation of conformational degrees of freedom. The analysis of the dynamics of the correlation matrix showed that rotations involving angles of the same type (φ–φ and ψ–ψ) occur predominantly in one direction. Rotations invoving different angles (φ–ψ) occur predominantly in opposite directions, so that the total macromolecule does not rotate. A significant reduction in the effective viscosity of the medium disrupts the correlation and makes the rotations stochastic, thus distorting the formation of the regular (helical) structure. The effects of correlated conformational motions are consequences of viscous friction forces. This conclusion agrees with our previous results that outlined the principle of the minimum rate of energy dissipation and the equipartition of energy dissipation rate between conformational degrees of freedom.