Tensorial elastic network model for protein dynamics: Integration of the anisotropic network model with bond‐bending and twist elasticities

We present a tensorial elastic network model (TNM) to describe the equilibrium fluctuations of proteins near their native fold structure. The model combines the anisotropic network model (ANM), bond bending elasticity, and backbone twist elasticity, and can predict both the isotropic fluctuations, similar to the Gaussian network model (GNM), and anisotropic fluctuations, similar to the ANM. TNM performs equally well for B‐factor predictions as GNM and predicts the anisotropy of B‐factors better than ANM. The model also outperforms the ANM in its predictability of the complete anisotropic displacement parameters. Proteins 2012; © 2012 Wiley Periodicals, Inc.     

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.     

Sampling of the native conformational ensemble of myoglobin via structures in different crystalline environments

Proteins sample multiple conformational substates in their native environment, but the process of crystallization selects the conformers that allow for close packing. The population of conformers can be shifted by varying the environment through a range of crystallization conditions, often resulting in different space groups and changes in the packing arrangements. Three high resolution structures of myoglobin (Mb) in different crystal space groups are presented, including one in a new space group P6(1)22 and two structures in space groups P2(1)2(1)2(1) and P6. We compare coordinates and anisotropic displacement parameters (ADPs) from these three structures plus an existing structure in space group P2(1). While the overall changes are small, there is substantial variation in several external regions with varying patterns of crystal contacts across the space group packing arrangements. The structural ensemble containing four different crystal forms displays greater conformational variance (Calpha rmsd of 0.54-0.79 A) in comparison to a collection of four Mb structures with different ligands and mutations in the same crystal form (Calpha rmsd values of 0.28-0.37 A). The high resolution of the data enables comparison of both the magnitudes and directions of ADPs, which are found to be suppressed by crystal contacts. A composite dynamic profile of Mb structural variation from the four structures was compared with an independent structural ensemble developed from NMR refinement. Despite the limitations and biases of each method, the ADPs of the crystallographic ensemble closely match the positional variance from the solution NMR ensemble with linear correlation of 0.8. This suggests that crystal packing selects conformers representative of the solution ensemble, and several different crystal forms give a more complete view of the plasticity of a protein structure.     

Anisotropic network model: systematic evaluation and a new web interface

MOTIVATION: The Anisotropic Network Model (ANM) is a simple yet powerful model for normal mode analysis of proteins. Despite its broad use for exploring biomolecular collective motions, ANM has not been systematically evaluated to date. A lack of a convenient interface has been an additional obstacle for easy usage. RESULTS: ANM has been evaluated on a large set of proteins to establish the optimal model parameters that achieve the highest correlation with experimental data and its limits of accuracy and applicability. Residue fluctuations in globular proteins are shown to be more accurately predicted than those in nonglobular proteins, and core residues are more accurately described than solvent-exposed ones. Significant improvement in agreement with experiments is observed with increase in the resolution of the examined structure. A new server for ANM calculations is presented, which offers flexible options for controlling model parameters and output formats, interactive animation of collective modes and advanced graphical features. AVAILABILITY: ANM server (http://www.ccbb.pitt.edu/anm)     

Normal mode refinement: crystallographic refinement of protein dynamic structure. II. Application to human lysozyme.

The dynamic structure of a protein, human lysozyme, is determined by the normal mode refinement of X-ray crystal structure. This method uses the normal modes of both internal and external motions to distinguish the real internal dynamics from the external terms such as lattice disorder, and gives an anisotropic and concerted picture of atomic fluctuations. The refinement is carried out with diffraction data of 5.0 to 1.8 A resolution, which are collected on an imaging plate. The results of the refinement show: (1) Debye-Waller factor consists of two parts, highly anisotropic internal fluctuations and almost isotropic external terms. The former is smaller than the latter by a factor of 0.72 in the scale of B-factor. Therefore, the internal dynamics cannot be recognized directly from the apparent electron density distribution. (2) The internal fluctuations show basically similar features as those predicted by the normal mode analysis, with almost the same amplitude and a similar level of anisotropy. (3) Correlations of fluctuations are detected between two lobes forming the active site cleft, which move simultaneously in opposite directions. This corresponds to the hinge-bending motion of lysozyme.     

Dynamics of proteins predicted by molecular dynamics simulations and analytical approaches: application to alpha-amylase inhibitor

The dynamics of alpha-amylase inhibitors has been investigated using molecular dynamics (MD) simulations and two analytical approaches, the Gaussian network model (GNM) and anisotropic network model (ANM). MD simulations use a full atomic approach with empirical force fields, while the analytical approaches are based on a coarse-grained single-site-per-residue model with a single-parameter harmonic potential between sufficiently close (r 相似文献   

Anharmonic Normal Mode Analysis of Elastic Network Model Improves the Modeling of Atomic Fluctuations in Protein Crystal Structures

Protein conformational dynamics, despite its significant anharmonicity, has been widely explored by normal mode analysis (NMA) based on atomic or coarse-grained potential functions. To account for the anharmonic aspects of protein dynamics, this study proposes, and has performed, an anharmonic NMA (ANMA) based on the C_α-only elastic network models, which assume elastic interactions between pairs of residues whose C_α atoms or heavy atoms are within a cutoff distance. The key step of ANMA is to sample an anharmonic potential function along the directions of eigenvectors of the lowest normal modes to determine the mean-squared fluctuations along these directions. ANMA was evaluated based on the modeling of anisotropic displacement parameters (ADPs) from a list of 83 high-resolution protein crystal structures. Significant improvement was found in the modeling of ADPs by ANMA compared with standard NMA. Further improvement in the modeling of ADPs is attained if the interactions between a protein and its crystalline environment are taken into account. In addition, this study has determined the optimal cutoff distances for ADP modeling based on elastic network models, and these agree well with the peaks of the statistical distributions of distances between C_α atoms or heavy atoms derived from a large set of protein crystal structures.     

An analysis of the influence of protein intrinsic dynamical properties on its thermal unfolding behavior

The influence of the protein topology-encoded dynamical properties on its thermal unfolding motions was studied in the present work. The intrinsic dynamics of protein topology was obtained by the anisotropic network model (ANM). The ANM has been largely used to investigate protein collective functional motions, but it is not well elucidated if this model can also reveal the preferred large-scale motions during protein unfolding. A small protein barnase is used as a typical case study to explore the relationship between protein topology-encoded dynamics and its unfolding motions. Three thermal unfolding simulations at 500 K were performed for barnase and the entire unfolding trajectories were sampled and partitioned into several windows. For each window, the preferred unfolding motions were investigated by essential dynamics analysis, and then associated with the intrinsic dynamical properties of the starting conformation in this window, which is detected by ANM. The results show that only a few slow normal modes imposed by protein structure are sufficient to give a significant overlap with the preferred unfolding motions. Especially, the large amplitude unfolding movements, which imply that the protein jumps out of a local energy basin, can be well described by a single or several ANM slow modes. Besides the global motions, it is also found that the local residual fluctuations encoded in protein structure are highly correlated with those in the protein unfolding process. Furthermore, we also investigated the relationship between protein intrinsic flexibility and its unfolding events. The results show that the intrinsic flexible regions tend to unfold early. Several early unfolding events can be predicted by analysis of protein structural flexibility. These results imply that protein structure-encoded dynamical properties have significant influences on protein unfolding motions.     

Concordance of residual dipolar couplings, backbone order parameters and crystallographic B-factors for a small alpha/beta protein: a unified picture of high probability, fast atomic motions in proteins

Using ensemble refinement of the third immunoglobulin binding domain (GB3) of streptococcal protein G (a small alpha/beta protein of 56 residues), we demonstrate that backbone (N-H, N-C', Calpha-Halpha, Calpha-C') residual dipolar coupling data in five independent alignment media, generalized order parameters from 15N relaxation data, and B-factors from a high-resolution (1.1A), room temperature crystal structure are entirely consistent with one another within experimental error. The optimal ensemble size representation is between four and eight, as assessed by complete cross-validation of the residual dipolar couplings. Thus, in the case of GB3, all three observables reflect the same low-amplitude anisotropic motions arising from fluctuations in backbone phi/psi torsion angles in the picosecond to nanosecond regime in both solution and crystalline environments, yielding a unified picture of fast, high-probability atomic motions in proteins. An understanding of these motions is crucial for understanding the impact of protein dynamics on protein function, since they provide part of the driving force for triggered conformational changes that occur, for example, upon ligand binding, signal transduction and enzyme catalysis.     

An atomic resolution structure for human fibroblast growth factor 1

A 1.10-A atomic resolution X-ray structure of human fibroblast growth factor 1 (FGF-1), a member of the beta-trefoil superfold, has been determined. The beta-trefoil is one of 10 fundamental protein superfolds and is the only superfold to exhibit 3-fold structural symmetry (comprising 3 "trefoil" units). The quality of the diffraction data permits unambiguous assignment of Asn, Gln, and His rotamers, Pro ring pucker, as well as refinement of atomic anisotropic displacement parameters (ADPs). The FGF-1 structure exhibits numerous core-packing defects, detectable using a 1.0-A probe radius. In addition to contributing to the relatively low thermal stability of FGF-1, these defects may also permit domain motions within the structure. The availability of refined ADPs allows a translation/libration/screw (TLS) analysis of putative rigid body domains. The TLS analysis shows that beta-strands 6-12 together form a rigid body, and there is a clear demarcation in TLS motions between the adjacent carboxyl- and amino-termini. Although separate from beta-strands 6-12, the individual beta-strands 1-5 do not exhibit correlated motions; thus, this region appears to be comparatively flexible. The heparin-binding contacts of FGF-1 are located within beta-strands 6-12; conversely, a significant portion of the receptor-binding contacts are located within beta-strands 1-5. Thus, the observed rigid body motion in FGF-1 appears related to the ligand-binding functionalities.     

Collective dynamics of EcoRI-DNA complex by elastic network model and molecular dynamics simulations

Anisotropic network model (ANM) is used to analyze the collective motions of restriction enzyme EcoRI in free form and in complex with DNA. For comparison, three independent molecular dynamics (MD) simulations, each of 1.5 ns duration, are also performed for the EcoRI-DNA complex in explicit water. Although high mobility (equilibrium fluctuations) of inner and outer loops that surround the DNA is consistent in both methods and experiments, MD runs sample different conformational subspaces from which reliable collective dynamics cannot be extracted. However, ANM employed on different conformations from MD simulations indicates very similar collective motions. The stems of the inner loops are quite immobile even in the free enzyme and form a large, almost fixed, pocket for DNA binding. As a result, the residues that make specific and non-specific interactions with the DNA exhibit very low fluctuations in the free enzyme. The vibrational entropy difference between the EcoRI complex and free protein + unkinked DNA is positive (favorable), which may partially counteract the unfavorable enthalpy difference of DNA kink formation. Dynamic domains in EcoRI complex and cross-correlations between residue fluctuations indicate possible means of communication between the distal active sites.     

The effects of rigid motions on elastic network model force constants

Elastic network models provide an efficient way to quickly calculate protein global dynamics from experimentally determined structures. The model's single parameter, its force constant, determines the physical extent of equilibrium fluctuations. The values of force constants can be calculated by fitting to experimental data, but the results depend on the type of experimental data used. Here, we investigate the differences between calculated values of force constants and data from NMR and X-ray structures. We find that X-ray B factors carry the signature of rigid-body motions, to the extent that B factors can be almost entirely accounted for by rigid motions alone. When fitting to more refined anisotropic temperature factors, the contributions of rigid motions are significantly reduced, indicating that the large contribution of rigid motions to B factors is a result of over-fitting. No correlation is found between force constants fit to NMR data and those fit to X-ray data, possibly due to the inability of NMR data to accurately capture protein dynamics.     

Nucleic-acid structural deformability deduced from anisotropic displacement parameters

The growing numbers of very well resolved nucleic-acid crystal structures with anisotropic displacement parameters provide an unprecedented opportunity to learn about the natural motions of DNA and RNA. Here we report a new Monte-Carlo approach that takes direct account of this information to extract the distortions of covalent structure, base pairing, and dinucleotide geometry intrinsic to regularly organized double-helical molecules. We present new methods to test the validity of the anisotropic parameters and examine the apparent deformability of a variety of structures, including several A, B, and Z DNA duplexes, an AB helical intermediate, an RNA, a ligand-DNA complex, and an enzyme-bound DNA. The rigid-body parameters characterizing the positions of the bases in the structures mirror the mean parameters found when atomic motion is taken into account. The base-pair fluctuations intrinsic to a single structure, however, differ from those extracted from collections of nucleic-acid structures, although selected base-pair steps undergo conformational excursions along routes suggested by the ensembles. The computations reveal surprising new molecular insights, such as the stiffening of DNA and concomitant separation of motions of contacted nucleotides on opposite strands by the binding of Escherichia coli endonuclease VIII, which suggest how the protein may direct enzymatic action.     

Normal mode refinement: crystallographic refinement of protein dynamic structure. I. Theory and test by simulated diffraction data.

A dynamic structure refinement method for X-ray crystallography, referred to as the normal mode refinement, is proposed. The Debye-Waller factor is expanded in terms of the low-frequency normal modes whose amplitudes and eigenvectors are experimentally optimized in the process of the crystallographic refinement. In this model, the atomic fluctuations are treated as anisotropic and concerted. The normal modes of the external motion (TLS model) are also introduced to cover the factors other than the internal fluctuations, such as the lattice disorder and diffusion. A program for the normal mode refinement (NM-REF) has been developed. The method has first been tested against simulated diffraction data for human lysozyme calculated by a Monte Carlo simulation. Applications of the method have demonstrated that the normal mode refinement has: (1) improved the fitting to the diffraction data, even with fewer adjustable parameters; (2) distinguished internal fluctuations from external ones; (3) determined anisotropic thermal factors; and (4) identified concerted fluctuations in the protein molecule.     

Anisotropic dynamics of the JE-2147-HIV protease complex: drug resistance and thermodynamic binding mode examined in a 1.09 A structure

The structure of HIV protease (HIV Pr) bound to JE-2147 (also named AG1776 or KNI-764) is determined here to 1.09 A resolution. This highest-resolution structure for HIV Pr allows refinement of anisotropic displacement parameters (ADPs) for all atoms. Clustering based on the directional information in ADPs defines two sets of subdomains such that within each set, subdomains undergo similar anisotropic motion. These sets are (a) the core of monomer A grouped with both substrate-binding flaps and (b) the core of monomer B coupled to both catalytic aspartates (25A/B). The four-stranded beta-sheet (1-4 A/B and 95-99 A/B) that forms a significant part of the dimer interface exhibits large anisotropic amplitudes that differ from those of the other sets of subdomains. JE-2147 is shown here to be a picomolar inhibitor (K(i) = 41 +/- 18 pM). The structure is used to interpret the mechanism of association of JE-2147, a second-generation inhibitor for which binding is enthalpically driven, with respect to first-generation inhibitors for which binding is predominantly entropically driven [Velazquez-Campoy, A., et al. (2001) Arch. Biochem. Biophys. 390, 169-175]. Relative to the entropically driven inhibitor complexes, the JE-2147-HIV Pr complex exhibits an approximately 0.5 A movement of the substrate flaps in toward the substrate, suggesting a more compatible enthalpically driven association. Domains of the protease identified by clustering of ADPs also suggest a model of enthalpy-entropy compensation for all HIV Pr inhibitors in which dynamic coupling of the flaps is offset by an increased level of motion of the beta-sheet domain of the dimer interface (1-4 A/B and 95-99 A/B).     

Normal mode refinement: crystallographic refinement of protein dynamic structure applied to human lysozyme.

A new method of dynamic structure refinement of protein x-ray crystallography, normal mode refinement, is developed. In this method the Debye-Waller factor is expanded in terms of the low-frequency normal modes and external normal modes, whose amplitudes and couplings are optimized in the process of crystallographic refinement. By this method, internal and external contributions to the atomic fluctuations can be separated. Also, anisotropic atomic fluctuations and their interatomic correlations can be determined experimentally even with a relatively small number of adjustable parameters. The method is applied to the analysis of experimental data of human lysozyme to reveal its dynamic structure.     

Global ribosome motions revealed with elastic network model

The motions of large systems such as the ribosome are not fully accessible with conventional molecular simulations. A coarse-grained, less-than-atomic-detail model such as the anisotropic network model (ANM) is a convenient informative tool to study the cooperative motions of the ribosome. The motions of the small 30S subunit, the larger 50S subunit, and the entire 70S assembly of the two subunits have been analyzed using ANM. The lowest frequency collective modes predicted by ANM show that the 50S subunit and 30S subunit are strongly anti-correlated in the motion of the 70S assembly. A ratchet-like motion is observed that corresponds well to the experimentally reported ratchet motion. Other slow modes are also examined because of their potential links to the translocation steps in the ribosome. We identify several modes that may facilitate the E-tRNA exiting from the assembly. The A-site t-RNA and P-site t-RNA are found to be strongly coupled and positively correlated in these slow modes, suggesting that the translocations of these two t-RNAs occur simultaneously, while the motions of the E-site t-RNA are less correlated, and thus less likely to occur simultaneously. Overall the t-RNAs exhibit relatively large deformations. Animations of these slow modes of motion can be viewed at.     

Application of Elastic Network Models to Proteins in the Crystalline State

Normal mode analysis using elastic network models has grown popular for probing the low-frequency collective dynamics of proteins and other biomolecular assemblies. In most previous studies, these models were validated by comparing calculated atomic fluctuations for isolated proteins with experimental temperature factors determined in the crystalline state, although there were also hints that including crystal contacts in the calculations has an impact on the comparison. In this study, a set of 83 ultra-high resolution crystal structures with experimentally determined anisotropic displacement parameters is used to evaluate several C_α-based elastic network models that either ignore or treat the crystal environment in different ways; the latter include using periodic boundary conditions defined with respect to the asymmetric unit or the primitive unit cell as well as using the Born-von Kármán boundary condition that accounts for lattice vibrations. For all elastic network models, treating the crystal environment leads to better agreement with experimental anisotropic displacement parameters with the Born-von Kármán boundary condition giving the best agreement. Atomic correlations over the entire protein are clearly affected by the presence of the crystal contacts and fairly sensitive to the way that the crystal environment is treated. These observations highlight the importance of properly treating the protein system in an environment consistent with experiment when either evaluating approximate protein models or using approximate dynamic models in structural refinement application types. Finally, investigation of the scaling behaviors of the cumulative density of states and the heat capacity indicates that there are still gaps between simplified elastic models and all-atom models for proteins.     

Normal-mode refinement of anisotropic thermal parameters for potassium channel KcsA at 3.2 A crystallographic resolution

We report a normal-mode method for anisotropic refinement of membrane-protein structures, based on a hypothesis that the global near-native-state disordering of membrane proteins in crystals follows low-frequency normal modes. Thus, a small set of modes is sufficient to represent the anisotropic thermal motions in X-ray crystallographic refinement. By applying the method to potassium channel KcsA at 3.2 A, we obtained a structural model with an improved fit with the diffraction data. Moreover, the improved electron density maps allowed for large structural adjustments for 12 residues in each subunit, including the rebuilding of 3 missing side chains. Overall, the anisotropic KcsA structure at 3.2 A was systematically closer to a 2.0 A KcsA structure, especially in the selectivity filter. Furthermore, the anisotropic thermal ellipsoids from the refinement revealed functionally relevant structural flexibility. We expect this method to be a valuable tool for structural refinement of many membrane proteins with moderate-resolution diffraction data.     

A physical picture of atomic motions within the Dickerson DNA dodecamer in solution derived from joint ensemble refinement against NMR and large-angle X-ray scattering data

The structure and dynamics of the Dickerson DNA dodecamer [5'd(CGCGAATTCGCG)2] in solution have been investigated by joint simulated annealing refinement against NMR and large-angle X-ray scattering data (extending from 0.25 to 3 A-1). The NMR data comprise an extensive set of hetero- and homonuclear residual dipolar coupling and 31P chemical shift anisotropy restraints in two alignment media, supplemented by NOE and 3J coupling data. The NMR and X-ray scattering data cannot be fully ascribed to a single structure representation, indicating the presence of anisotropic motions that impact the experimental observables in different ways. Refinement with ensemble sizes (Ne) of >or=2 to represent the atomic motions reconciles all the experimental data within measurement error. Cross validation against both the dipolar coupling and X-ray scattering data suggests that the optimal ensemble size required to account for the current data is 4. The resulting ensembles permit one to obtain a detailed view of the conformational space sampled by the dodecamer in solution and permit one to analyze fluctuations in helicoidal parameters, sugar puckers, and BI-BII backbone transitions and to obtain quantitative metrics of atomic motion such as generalized order parameters and thermal B factors. The calculated order parameters are in good agreement with experimental order parameters obtained from 13C relaxation measurements. Although DNA behaves as a relatively rigid rod with a persistence length of approximately 150 bp, dynamic conformational heterogeneity at the base pair level is functionally important since it readily permits optimization of intermolecular protein-DNA interactions.