Dynamic and structural analysis of isotropically distributed molecular ensembles.

An efficient new method is presented for the characterization of motional correlations derived from a set of protein structures without requiring the separation of overall and internal motion. In this method, termed isotropically distributed ensemble (IDE) analysis, each structure is represented by an ensemble of isotropically distributed replicas corresponding to the situation found in an isotropic protein solution. This leads to a covariance matrix of the cartesian atomic positions with elements proportional to the ensemble average of scalar products of the position vectors with respect to the center of mass. Diagonalization of the covariance matrix yields eigenmodes and amplitudes that describe concerted motions of atoms, including overall rotational and intramolecular dynamics. It is demonstrated that this covariance matrix naturally distinguishes between "rigid" and "mobile" parts without necessitating a priori selection of a reference structure and an atom set for the orientational alignment process. The method was applied to the analysis of a 5-ns molecular dynamics trajectory of native ubiquitin and a 40-ns trajectory of a partially folded state of ubiquitin. The results were compared with essential dynamics analysis. By taking advantage of the spherical symmetry of the IDE covariance matrix, more than a 10-fold speed up is achieved for the computation of eigenmodes and mode amplitudes. IDE analysis is particularly suitable for studying the correlated dynamics of flexible and large molecules.     

Essential dynamics of proteins

Analysis of extended molecular dynamics (MD) simulations of lysozyme in vacuo and in aqueous solution reveals that it is possible to separate the configurational space into two subspaces: (1) an “essential” subspace containing only a few degrees of freedom in which anharmonic motion occurs that comprises most of the positional fluctuations; and (2) the remaining space in which the motion has a narrow Gaussian distribution and which can be considered as “physically constrained.” If overall translation and rotation are eliminated, the two spaces can be constructed by a simple linear transformation in Cartesian coordinate space, which remains valid over several hundred picoseconds. The transformation follows from the covariance matrix of the positional deviations. The essential degrees of freedom seem to describe motions which are relevant for the function of the protein, while the physically constrained subspace merely describes irrelevant local fluctuations. The near-constraint behavior of the latter subspace allows the separation of equations of motion and promises the possibility of investigating independently the essential space and performing dynamic simulations only in this reduced space. © 1993 Wiley-Liss, Inc.     

Optimal Superpositioning of Flexible Molecule Ensembles

Analysis of the internal dynamics of a biological molecule requires the successful removal of overall translation and rotation. Particularly for flexible or intrinsically disordered peptides, this is a challenging task due to the absence of a well-defined reference structure that could be used for superpositioning. In this work, we started the analysis with a widely known formulation of an objective for the problem of superimposing a set of multiple molecules as variance minimization over an ensemble. A negative effect of this superpositioning method is the introduction of ambiguous rotations, where different rotation matrices may be applied to structurally similar molecules. We developed two algorithms to resolve the suboptimal rotations. The first approach minimizes the variance together with the distance of a structure to a preceding molecule in the ensemble. The second algorithm seeks for minimal variance together with the distance to the nearest neighbors of each structure. The newly developed methods were applied to molecular-dynamics trajectories and normal-mode ensembles of the Aβ peptide, RS peptide, and lysozyme. These new (to our knowledge) superpositioning methods combine the benefits of variance and distance between nearest-neighbor(s) minimization, providing a solution for the analysis of intrinsic motions of flexible molecules and resolving ambiguous rotations.     

Clustering algorithms for identifying core atom sets and for assessing the precision of protein structure ensembles

An important open question in the field of NMR-based biomolecular structure determination is how best to characterize the precision of the resulting ensemble of structures. Typically, the RMSD, as minimized in superimposing the ensemble of structures, is the preferred measure of precision. However, the presence of poorly determined atomic coordinates and multiple "RMSD-stable domains"--locally well-defined regions that are not aligned in global superimpositions--complicate RMSD calculations. In this paper, we present a method, based on a novel, structurally defined order parameter, for identifying a set of core atoms to use in determining superimpositions for RMSD calculations. In addition we present a method for deciding whether to partition that core atom set into "RMSD-stable domains" and, if so, how to determine partitioning of the core atom set. We demonstrate our algorithm and its application in calculating statistically sound RMSD values by applying it to a set of NMR-derived structural ensembles, superimposing each RMSD-stable domain (or the entire core atom set, where appropriate) found in each protein structure under consideration. A parameter calculated by our algorithm using a novel, kurtosis-based criterion, the epsilon-value, is a measure of precision of the superimposition that complements the RMSD. In addition, we compare our algorithm with previously described algorithms for determining core atom sets. The methods presented in this paper for biomolecular structure superimposition are quite general, and have application in many areas of structural bioinformatics and structural biology.     

Motion of groups of atoms in DNA studied by molecular dynamics simulation

The analysis of Molecular Dynamics simulations of two double stranded oligonucleotides is presented in terms of motions of quasi rigid subunits. First, a strategy is presented for grouping atoms submitted to concerted internal motions. The method is based on the analysis of the interatomic distance RMS matrix. It is found that each nucleotide can reasonably be decomposed into 3 or 4 rigid groups of atoms depending on the tolerance of the definition of a rigid body. In the second part, the different kinds of motions of the subunits (deformation, translation and rotation) are studied in terms of correlation using the canonical correlation analysis of data. It is shown that the residual deformation of any subunit does not influence the translational and rotational motions of the others, except perhaps for long time dynamics. Received: 3 October 1997 / Revised version: 22 December 1997 / Accepted: 23 December 1997     

Comparative methods with sampling error and within-species variation: contrasts revisited and revised

Comparative methods analyses have usually assumed that the species phenotypes are the true means for those species. In most analyses, the actual values used are means of samples of modest size. The covariances of contrasts then involve both the covariance of evolutionary changes and a fraction of the within-species phenotypic covariance, the fraction depending on the sample size for that species. Ives et al. have shown how to analyze data in this case when the within-species phenotypic covariances are known. The present model allows them to be unknown and to be estimated from the data. A multivariate normal statistical model is used for multiple characters in samples of finite size from species related by a known phylogeny, under the usual Brownian motion model of change and with equal within-species phenotypic covariances. Contrasts in each character can be obtained both between individuals within a species and between species. Each contrast can be taken for all of the characters. These sets of contrasts, each the same contrast taken for different characters, are independent. The within-set covariances are unequal and depend on the unknown true covariance matrices. An expectation-maximization algorithm is derived for making a reduced maximum likelihood estimate of the covariances of evolutionary change and the within-species phenotypic covariances. It is available in the Contrast program of the PHYLIP package. Computer simulations show that the covariances are biased when the finiteness of sample size is not taken into account and that using the present model corrects the bias. Sampling variation reduces the power of inference of covariation in evolution of different characters. An extension of this method to incorporate estimates of additive genetic covariances from a simple genetic experiment is also discussed.     

Normal mode analysis of macromolecular motions in a database framework: developing mode concentration as a useful classifying statistic

We investigated protein motions using normal modes within a database framework, determining on a large sample the degree to which normal modes anticipate the direction of the observed motion and were useful for motions classification. As a starting point for our analysis, we identified a large number of examples of protein flexibility from a comprehensive set of structural alignments of the proteins in the PDB. Each example consisted of a pair of proteins that were considerably different in structure given their sequence similarity. On each pair, we performed geometric comparisons and adiabatic-mapping interpolations in a high-throughput pipeline, arriving at a final list of 3,814 putative motions and standardized statistics for each. We then computed the normal modes of each motion in this list, determining the linear combination of modes that best approximated the direction of the observed motion. We integrated our new motions and normal mode calculations in the Macromolecular Motions Database, through a new ranking interface at http://molmovdb.org. Based on the normal mode calculations and the interpolations, we identified a new statistic, mode concentration, related to the mathematical concept of information content, which describes the degree to which the direction of the observed motion can be summarized by a few modes. Using this statistic, we were able to determine the fraction of the 3,814 motions where one could anticipate the direction of the actual motion from only a few modes. We also investigated mode concentration in comparison to related statistics on combinations of normal modes and correlated it with quantities characterizing protein flexibility (e.g., maximum backbone displacement or number of mobile atoms). Finally, we evaluated the ability of mode concentration to automatically classify motions into a variety of simple categories (e.g., whether or not they are "fragment-like"), in comparison to motion statistics. This involved the application of decision trees and feature selection (particular machine-learning techniques) to training and testing sets derived from merging the "list" of motions with manually classified ones.     

Fluorescence depolarization and rotational modes of tyrosine in bovine pancreatic trypsin inhibitor

The fluorescence of bovine pancreatic trypsin inhibitor (BPTI) is due to one or more of its four tyrosine residues. Observations of the stationary polarization of the fluorescence over a large range of temperatures and viscosities permit the demonstration of at least three modes of tyrosine rotation, and perhaps an ultrafast fourth one. The slowest mode is one of motion of the whole molecule; the second, a much faster motion limited to an amplitude of 11 degrees, is not changed by quenching of the fluorescence through addition of citrate and is therefore ascribed to the motion of internal tyrosines of BPTI. The third mode of motion is faster still; it has an amplitude similar to that of the second and, being sensitive to citrate quenching, is attributed to the rotation of the external tyrosine residue. A residual depolarization corresponding to a rotational amplitude of 22 degrees is deduced by comparison of the polarizations of BPTI and tyrosine dissolved in 80% glycerol-water at -40 degrees C. It is in accord in amplitude with the picosecond tyrosine rotations predicted by Karplus and collaborators from molecular dynamics computer simulations, but it could also originate, in whole or in part, from electronic energy transfer among the tyrosines.     

Comparison of tRNA motions in the free and ribosomal bound structures

A general method is presented that allows the separation of the rigid body motions from the nonrigid body motions of structural subunits when bound in a complex. The application presented considers the motions of the tRNAs: free, bound to the ribosome and to a synthase. We observe that both the rigid body and nonrigid body motions of the structural subunits are highly controlled by the large ribosomal assembly and are important for the functional motions of the assembly. For the intact ribosome, its major parts, the 30S and the 50S subunits, are found to have counterrotational motions in the first few slowest modes, which are consistent with the experimentally observed ratchet motion. The tRNAs are found to have on average approximately 72-75% rigid body motions and principally translational motions within the first 100 slow modes of the complex. Although the three tRNAs exhibit different apparent total motions, after the rigid body motions are removed, the remaining internal motions of all three tRNAs are essentially the same. The direction of the translational motions of the tRNAs are in the same direction as the requisite translocation step, especially in the first slowest mode. Surprisingly the small intrinsically flexible mRNA has all of its internal motions completely inhibited and shows mainly a rigid-body translation in the slow modes of the ribosome complex. On the other hand, the required nonrigid body motions of the tRNA during translocation reveal that the anticodon-stem-loop, as well as the acceptor arm, of the tRNA enjoy a large mobility but act as rigid structural units. In summary, the ribosome exerts its control by enforcing rigidity in the functional parts of the tRNAs as well as in the mRNA.     

Protein motions at zero-total angular momentum: the importance of long-range correlations

A constant-energy molecular dynamics simulation is used to monitor protein motion at zero-total angular momentum. With a simple protein model, it is shown that overall rotation is possible at zero-total angular momentum as a result of flexibility. Since the rotational motion is negligible on a time scale of 1000 reduced time units, the essentially rotation-free portion of the trajectory provides an unbiased test of the common approximate methods for separating overall rotation from internal motions by optimal superposition. Removing rotation by minimizing the root-mean-square deviation (RMSD) for the entire system is found to be more appropriate than using the RMSD for only the more rigid part of the system. The results verify the existence of positive cross-correlation in the motions of atoms separated by large distances.     

Systematic Study of Anharmonic Features in a Principal Component Analysis of Gramicidin A

We use principal component analysis (PCA) to detect functionally interesting collective motions in molecular-dynamics simulations of membrane-bound gramicidin A. We examine the statistical and structural properties of all PCA eigenvectors and eigenvalues for the backbone and side-chain atoms. All eigenvalue spectra show two distinct power-law scaling regimes, quantitatively separating large from small covariance motions. Time trajectories of the largest PCs converge to Gaussian distributions at long timescales, but groups of small-covariance PCs, which are usually ignored as noise, have subdiffusive distributions. These non-Gaussian distributions imply anharmonic motions on the free-energy surface. We characterize the anharmonic components of motion by analyzing the mean-square displacement for all PCs. The subdiffusive components reveal picosecond-scale oscillations in the mean-square displacement at frequencies consistent with infrared measurements. In this regime, the slowest backbone mode exhibits tilting of the peptide planes, which allows carbonyl oxygen atoms to provide surrogate solvation for water and cation transport in the channel lumen. Higher-frequency modes are also apparent, and we describe their vibrational spectra. Our findings expand the utility of PCA for quantifying the essential features of motion on the anharmonic free-energy surface made accessible by atomistic molecular-dynamics simulations.     

Mandibular kinematics represented by a non-orthogonal floating axis joint coordinate system

There are many methods used to represent joint kinematics (e.g., roll, pitch, and yaw angles; instantaneous center of rotation; kinematic center; helical axis). Often in biomechanics internal landmarks are inferred from external landmarks. This study represents mandibular kinematics using a non-orthogonal floating axis joint coordinate system based on 3-D geometric models with parameters that are "clinician friendly" and mathematically rigorous. Kinematics data for two controls were acquired from passive fiducial markers attached to a custom dental clutch. The geometric models were constructed from MRI data. The superior point along the arc of the long axis of the condyle was used to define the coordinate axes. The kinematic data and geometric models were registered through fiducial markers visible during both protocols. The mean absolute maxima across the subjects for sagittal rotation, coronal rotation, axial rotation, medial-lateral translation, anterior-posterior translation, and inferior-superior translation were 34.10 degrees, 1.82 degrees, 1.14 degrees, 2.31, 21.07, and 6.95 mm, respectively. All the parameters, except for one subject's axial rotation, were reproducible across two motion recording sessions. There was a linear correlation between sagittal rotation and translation, the dominant motion plane, with approximately 1.5 degrees of rotation per millimeter of translation. The novel approach of combining the floating axis system with geometric models succinctly described mandibular kinematics with reproducible and clinician friendly parameters.     

Large‐scale analysis of the dynamics of enzymes

Protein enzymes enable the cell to execute chemical reactions in short time by accelerating the rate of the reactions in a selective manner. The motions or dynamics of the enzymes are essential for their function. Comparison of the dynamics of a set of 1247 nonhomologous enzymes was performed. For each enzyme, the slowest modes of motion are calculated using the Gaussian network model (GNM) and they are globally aligned. Alignment is done using the dynamic programming algorithm of Needleman and Wunsch, commonly used for sequence alignment. Only 96 pairs of proteins were identified to have three similar GNM slow modes with 63 of them having a similar structure. The most frequent slowest mode of motion describes a two domains anticorrelated motion that characterizes at least 23% of the enzymes. Therefore, dynamics uniqueness cannot be accounted for by the slowest mode itself but rather by the combination of several slow modes. Different quaternary structure packing can restrain the motion of enzyme subunits differently and may serve as another mechanism that increases the dynamics uniqueness. Proteins 2013; 81:1910–1918. © 2013 Wiley Periodicals, Inc.     

A coarse-grained normal mode approach for macromolecules: an efficient implementation and application to Ca(2+)-ATPase

A block normal mode (BNM) algorithm, originally proposed by Tama et al., (Proteins Struct. Func. Genet. 41:1-7, 2000) was implemented into the simulation program CHARMM. The BNM approach projects the hessian matrix into local translation/rotation basis vectors and, therefore, dramatically reduces the size of the matrix involved in diagonalization. In the current work, by constructing the atomic hessian elements required in the projection operation on the fly, the memory requirement for the BNM approach has been significantly reduced from that of standard normal mode analysis and previous implementation of BNM. As a result, low frequency modes, which are of interest in large-scale conformational changes of large proteins or protein-nucleic acid complexes, can be readily obtained. Comparison of the BNM results with standard normal mode analysis for a number of small proteins and nucleic acids indicates that many properties dominated by low frequency motions are well reproduced by BNM; these include atomic fluctuations, the displacement covariance matrix, vibrational entropies, and involvement coefficients for conformational transitions. Preliminary application to a fairly large system, Ca(2+)-ATPase (994 residues), is described as an example. The structural flexibility of the cytoplasmic domains (especially domain N), correlated motions among residues on domain interfaces and displacement patterns for the transmembrane helices observed in the BNM results are discussed in relation to the function of Ca(2+)-ATPase. The current implementation of the BNM approach has paved the way for developing efficient sampling algorithms with molecular dynamics or Monte Carlo for studying long-time scale dynamics of macromolecules.     

Calculation of hydrodynamic properties of small nucleic acids from their atomic structure

Hydrodynamic properties (translational diffusion, sedimentation coefficients and correlation times) of short B-DNA oligonucleotides are calculated from the atomic-level structure using a bead modeling procedure in which each non-hydrogen atom is represented by a bead. Using available experimental data of hydrodynamic properties for several oligonucleotides, the best fit for the hydrodynamic radius of the atoms is found to be ~2.8 Å. Using this value, the predictions for the properties corresponding to translational motion and end-over-end rotation are accurate to within a few percent error. Analysis of NMR correlation times requires accounting for the internal flexibility of the double helix, and allows an estimation of ~0.85 for the Lipari–Szabo generalized order parameter. Also, the degree of hydration can be determined from hydrodynamics, with a result of ~0.3 g (water)/g (DNA). These numerical results are quite similar to those found for globular proteins. If the hydrodynamic model for the short DNA is simply a cylindrical rod, the predictions for overall translation and rotation are slightly worse, but the NMR correlation times and the degree of hydration, which depend more on the cross-sectional structure, are more severely affected.     

Thermal motion of whole protein molecules in protein solids

A preliminary analysis is presented of whether and to what extent thermal motion of protein molecules as rigid bodies contributes to the mobility found in X-ray and M?ssbauer studies. A simple theory is advanced enabling the estimation of mean-square amplitudes of translational and librational motion of the protein molecules in crystals and amorphous glasses from the experimental data on the elastic properties of these solids. The values calculated and their dependence on the crystal packing, temperature and hydration level were found to be in good accord with the data of X-ray analysis and M?ssbauer spectroscopy. The external modes were concluded to contribute significantly to the values of mean value of chi 2 measured by the last two methods meaning that the conventional amplitudes of internal motion in proteins were overestimated. The real average amplitude of thermal motion of atoms in the protein interior should be close to that in molecular crystals, in accordance with the crystal-like packing of atoms inside the protein globule and some other "crystal-like" physical characteristics such as Young's modulus, adiabatic compressibility and thermal expansion coefficient. Factors are discussed which determine the temperature dependence of the amplitudes of external and internal modes of protein motion.     

On selection criteria and estimation of parameters when the variance is heterogeneous

Summary Procedures for ranking candidates for selection and for estimating genetic and environmental parameters when variances are heterogeneous are discussed. The best linear unbiased predictor (BLUP) accounts automatically for heterogeneous variance provided that the covariance structure is known and that the assumptions of the model hold. Under multivariate normality BLUP allowing for heterogeneous variance maximizes expected genetic progress. Examples of application of BLUP to selection when residual or genetic variances are heterogeneous are given. Restricted maximum likelihood estimation of heterogeneous variances and covariances via the expectation-maximization algorithm is presented.     

Solution Structure of the Lipoyl Domain of the 2-Oxoglutarate Dehydrogenase Complex fromAzotobacter vinelandii

The three-dimensional solution structure of the lipoyl domain of the 2-oxoglutarate dehydrogenase complex fromAzotobacter vinelandiihas been determined from nuclear magnetic resonance data by using distance geometry and dynamical simulated annealing refinement. The structure determination is based on a total of 580 experimentally derived distance constraints and 65 dihedral angle constraints. The solution structure is represented by an ensemble of 25 structures with an average root-mean-square deviation between the individual structures of the ensemble and the mean coordinates of 0.71 Å for backbone atoms and 1.08 Å for all heavy atoms. The overall fold of the lipoyl domain is that of a β-barrel-sandwich hybrid. It consists of two almost parallel four-stranded anti-parallel β-sheets formed around a well-defined hydrophobic core, with a central position of the single tryptophan 21. The lipoylation site, lysine 42, is found in a β-turn at the far end of one of the sheets, and is close in space to a solvent-exposed loop comprising residues 7 to 15. The lipoyl domain displays a remarkable internal symmetry that projects one β-sheet onto the other β-sheet after rotation of approximately 180° about a 2-fold rotational symmetry axis. There is close structural similarity between the structure of this 2-oxoglutarate dehydrogenase complex lipoyl domain and the structures of the lipoyl domains of pyruvate dehydrogenase complexes fromBacillus stearothermophilusandEscherichia coli, and conformational differences occur primarily in a solvent-exposed loop close in space to the lipoylation site. The lipoyl domain structure is discussed in relation to the process of molecular recognition of lipoyl domains by their parent 2-oxo acid dehydrogenase.