Large Concerted Motions Research Articles

Glycosaminoglycan Recognition of Neutrophil‐Activating Chemokines

Neutrophil‐activating chemokines are known for their fundamental and diverse role in recruiting immune cells to the site of inflammation/infection. Extracellular glycosaminoglycans (GAGs), especially heparan sulfate, sequester chemokines and engineer a concentration gradient to regulate the immune cell trafficking. Not much is known about the structural and dynamical aspects of chemokine – GAG interactions. We have studied the interaction of several human and mouse chemokines, including CXCL1, CXCL2, and CXCL5 with heparin/heparan sulfate oligosaccharides using computational (modeling and dynamics) and biophysical (NMR & ITC) techniques. These studies have led to fine atomistic level understanding of GAG recognition of chemokines. First, our studies show that even though these chemokines are structurally homologous, they exhibit striking differences in their recognition of GAGs. Diverse binding geometries have been observed in which a GAG chain may bind each monomer present in a chemokine dimer or a GAG chain may bind two non‐overlapping sites by spanning across the dimer interface. Second, molecular dynamics simulations show large concerted motions of the Lys/Arg side chains bound to sulfate/carboxylates of GAGs, which imply translational motion of a GAG chain over the binding interface leading to considerable plasticity in binding. At the same time, a GAG chain may have restricted motion arising from specific interactions with a group of core amino acid residues. Third, a large number of hydrogen bond contributing residues occur in chemokines, e.g., in N‐loop, 40s turn, β3 strand and C‐helix, which provide strong/weak interactions to stabilize the complex. Interestingly, GAG recognition of chemokines is also brought about by water‐mediated interactions between amino acid side chains and sulfates/carboxylates, which introduce additional stability to GAG's dynamical motion. Fourth, relative importance of individual amino acid residues correlate well with experimental observations, thus implying that computational trajectories can be reliably used to understand GAG–chemokine systems. In addition, the computational trajectories can lead to identification of preferred chain length, mode of binding, and dynamical motions of GAGs in recognition of neutrophil‐activating chemokines. Our work contributes fundamental understanding on the structural biology of GAG–chemokine interactions.Support or Funding InformationThis work was supported in part by grants from the NIH including HL107152, HL090586 and HL128639This abstract is from the Experimental Biology 2019 Meeting. There is no full text article associated with this abstract published in The FASEB Journal.

Insight derived from molecular dynamics simulation into dynamics and molecular motions of cuticle-degrading serine protease Ver112

Cuticle-degrading serine protease Ver112, which derived from a nematophagous fungus Lecanicillium psalliotae, has been exhibited to have high cuticle-degrading and nematicidal activities. We have performed molecular dynamics (MD) simulation based on the crystal structure of Ver112 to investigate its dynamic properties and large-scale concerted motions. The results indicate that the structural core of Ver112 shows a small fluctuation amplitude, whereas the substrate binding sites, and the regions close to and opposite the substrate binding sites experience significant conformational fluctuations. The large concerted motions obtained from essential dynamics (ED) analysis of MD trajectory can lead to open or close of the substrate binding sites, which are proposed to be linked to the functional properties of Ver112, such as substrate binding, orientation, catalytic, and release. The significant motion in the loop regions that is located opposite the binding sites are considered to play an important role in modulating the dynamics of the substrate binding sites. Furthermore, the bottom of free energy landscape (FEL) of Ver112 are rugged, which is mainly caused by the fluctuations of substrate binding regions and loops located opposite the binding site. In addition, the mechanism underlying the high flexibility and catalytic activity of Ver112 was also discussed. Our simulation study complements the biochemical and structural studies, and provides insight into the dynamics-function relationship of cuticle-degrading serine protease Ver112.

Insight derived from molecular dynamics simulations into molecular motions, thermodynamics and kinetics of HIV-1 gp120.

Although the crystal structures of the HIV-1 gp120 core bound and pre-bound by CD4 are known, the details of dynamics involved in conformational equilibrium and transition in relation to gp120 function have remained elusive. The homology models of gp120 comprising the N- and C-termini and loops V3 and V4 in the CD4-bound and CD4-unbound states were built and subjected to molecular dynamics (MD) simulations to investigate the differences in dynamic properties and molecular motions between them. The results indicate that the CD4-bound gp120 adopted a more compact and stable conformation than the unbound form during simulations. For both the unbound and bound gp120, the large concerted motions derived from essential dynamics (ED) analyses can influence the size/shape of the ligand-binding channel/cavity of gp120 and, therefore, were related to its functional properties. The differences in motion direction between certain structural components of these two forms of gp120 were related to the conformational interconversion between them. The free energy calculations based on the metadynamics simulations reveal a more rugged and complex free energy landscape (FEL) for the unbound than for the bound gp120, implying that gp120 has a richer conformational diversity in the unbound form. The estimated free energy difference of ∼−6.0 kJ/mol between the global minimum free energy states of the unbound and bound gp120 indicates that gp120 can transform spontaneously from the unbound to bound states, revealing that the bound state represents a high-probability “ground state” for gp120 and explaining why the unbound state resists crystallization. Our results provide insight into the dynamics-and-function relationship of gp120, and facilitate understandings of the thermodynamics, kinetics and conformational control mechanism of HIV-1 gp120.

153 Wonderful roles of the entropy in protein dynamics, binding and folding

The entropy, which is central to the second law of thermodynamics, determines that the thermal energy always flows spontaneously from regions of higher temperature to regions of lower temperature. In the protein–solvent thermodynamic system, the entropy is defined as a measure of how evenly the thermal energy would distribute over the entire system (Liu et al., 2012). Such tendency to distribute energy as evenly as possible will reduce the state of order of the initial system, and hence, the entropy can be regarded as an expression of the disorder, or randomness of the system (Yang et al., 2012). For a protein–solvent system under a constant solvent condition, the origin of entropy is the thermal energy stored in atoms, which makes atoms jostle around and bump onto one another, thus leading to vibrations of the covalent bonds connecting two atoms (occurring on the fs timescale) and the rotational and translational motions of amino acid side chain groups (occurring on ps timescale) and water molecules. These motions break the noncovalent bonds around structural regions that are weakly constrained thereby triggering the competitive interactions among residues or between residues and water molecules leading ultimately to the loop motions (occurring on ns timescale) around the protein surface. The loop motions can further transmit either through the water network around the protein surface or via specific structural components (such as the hinge-bending regions) over the entire protein molecule leading to large concerted motions (occurring on μs to s timescales) that are most relevant to protein functions (Amadei, Linssen & Berendsen, 1993; Tao, Rao & Liu, 2010). Thus, the multiple hierarchies of the protein dynamics on distinct timescales (Henzler-Wildman & Kern, 2007) are a consequence of the cascade amplification of the microscopic motions of atoms and groups for which the entropy originating from atomic thermal energy is most fundamental. In the case of protein–ligand binding, the importance of the entropy is embodied in the following aspects. (i) The release of the water molecule kinetic energy (which is a process of the solvent entropy maximization) will cause Brownian motions of individual water molecules which result in strong Brownian bombardments to solute molecules causing molecule wanders/diffusions and subsequent accident contacts/collisions between proteins and ligands. (ii) Such collisions will inevitably cause water molecule displacement and, if the contact interfaces are properly complementary, the requirement to increase the solvent entropy would further displace the water network around the binding interfaces thus leading to the formation the initial protein-ligand complex. (iii) In the initial complex, the loose association of the two partners provide the opportunity for protein to increase conformational entropy, thus triggering the conformational adjustments through competitive interaction between protein residues and ligand, leading ultimately to the formation of tightly associated complex (Liu et al., 2012). In the protein folding process, the first stage, i.e. the rapid hydrophobic collapse (Agashe, Shastry & Udgaonkar, 1995; Dill, 1985), is in fact driven by the effect of the solvent entropy maximization. Specifically, the requirement to maintain as many as possible the dynamic hydrogen bonds among the water molecules will squeeze/sequestrate the hydrophobic amino acid side chains into the interior of the folding intermediates and expose the polar/charged side chains onto the intermediate surface. This will minimize the solvent accessible surface area of the folding intermediates and as thus maximize the entropy of the solvent. The resulting molten globule states (Ohgushi & Wada, 1983) may contain a few secondary structural components and native tertiary contacts, while many native contacts, or close residue–residue interactions present in the native state have not yet formed. However, the nature to increase the protein conformational entropy can trigger a further conformational adjustment process, i.e. the conformational entropy increase breaks the transient secondary or tertiary contacts and triggers the competitive interactions among protein residues and between residues and water. This process may repeat many rounds until the negative enthalpy change resulting from the noncovalent formations can overcompensate for protein conformational entropy loss. In summary, we consider that the tendency to maximize the entropy of the protein–solvent system, which originates from the atomic thermal energy, is the most fundamental driving factor for protein folding, binding, and dynamics, whereas the enthalpy reduction, an opposing factor that tends to make the system become ordered, can compensate for the effect of entropy loss to ultimately allow the system to reach equilibrium at the free energy minima, either global or local.

Insight Derived from Molecular Dynamics Simulation into Substrate-Induced Changes in Protein Motions of Proteinase K

Because of the significant industrial, agricultural and biotechnological importance of serine protease proteinase K, it has been extensively investigated using experimental approaches such as X-ray crystallography, site-directed mutagenesis and kinetic measurement. However, detailed aspects of enzymatic mechanism such as substrate binding, release and relevant regulation remain unstudied. Molecular dynamics (MD) simulations of the proteinase K alone and in complex with the peptide substrate AAPA were performed to investigate the effect of substrate binding on the dynamics/molecular motions of proteinase K. The results indicate that during simulations the substrate-complexed proteinase K adopt a more compact and stable conformation than the substrate-free form. Further essential dynamics (ED) analysis reveals that the major internal motions are confined within a subspace of very small dimension. Upon substrate binding, the overall flexibility of the protease is reduced; and the noticeable displacements are observed not only in substrate-binding regions but also in regions opposite the substrate-binding groove/pockets. The dynamic pockets caused by the large concerted motions are proposed to be linked to the substrate recognition, binding, orientation and product release; and the significant displacements in regions opposite the binding groove/pockets are considered to play a role in modulating the dynamics of enzyme-substrate interaction. Our simulation results complement the biochemical and structural studies, highlighting the dynamic mechanism of the functional properties of proteinase K.

Insights derived from molecular dynamics simulation into the molecular motions of serine protease proteinase K

Serine protease proteinase K, a member of the subtilisin family of enzymes, is of significant industrial, agricultural and biotechnological importance. Despite the wealth of structural information about proteinase K provided by static X-ray structures, a full understanding of the enzymatic mechanism requires further insight into the dynamic properties of this enzyme. Molecular dynamics simulations and essential dynamics (ED) analysis were performed to investigate the molecular motions in proteinase K. The results indicate that the internal core of proteinase K is relatively rigid, whereas the surface-exposed loops, most notably the substrate-binding regions, exhibit considerable conformational fluctuations. Further ED analysis reveals that the large concerted motions in the substrate-binding regions cause opening/closing of the substrate-binding pockets, thus supporting the proposed induced-fit mechanism of substrate binding. The distinct electrostatic/hydrogen-bonding interactions between Asp39 and His69 and between His69 and Ser224 within the catalytic triad lead to different thermal motions and orientations of these three catalytic residues, which can be related to their different functional roles in the catalytic process. Statistical analyses of the geometrical/functional properties as well as evolutionary conservation of the glycines in proteinase K-like proteins reveal that glycines may play an important role in determining the folding architecture and structural flexibility of this class of enzymes. Our simulation study complements the biochemical and structural studies and provides new insights into the dynamic structural basis of the functional properties of this class of enzymes.

Conformational substates in different crystal forms of the photoactive yellow protein--correlation with theoretical and experimental flexibility.

The conformational changes during the photocycle of the photoactive yellow protein have been the subject of many recent studies. Spectroscopic measurements have shown that the photocycle also occurs in a crystalline environment, and this has been the basis for subsequent Laue diffraction and cryocrystallographic studies. These studies have shown that conformational changes during the photocycle are limited to the chromophore and its immediate environment. However, spectroscopic studies suggest the presence of large conformational changes in the protein. Here, we address this apparent discrepancy in two ways. First, we obtain a description of large concerted motions in the ground state of the yellow protein from NMR data and theoretical calculations. Second, we describe the high-resolution structure of the yellow protein crystallized in a different space group. The structure of the yellow protein differs significantly between the two crystal forms. We show that these differences can be used to obtain a description of the flexibility of the protein that is consistent with the motions observed in solution.

Protein Folding Pathways Revealed by Essential Dynamics Sampling.

The characterization of the protein folding process represents one of the major challenges in molecular biology. Here, a method to simulate the folding process of a protein to its native state is reported, the essential dynamics sampling (EDS) method, and is successfully applied to detecting the correct folding pathways of two small proteins, the all-β SH3 domain of Src tyrosine kinase transforming protein (SH3) and the α/β B1 domain of streptococcal protein G (GB1). The main idea of the method is that a subset of the natural modes of fluctuation in the native state is key in directing the folding process. A biased molecular dynamics simulation is performed, in which the restrained degrees of freedom are chosen among those obtained by a principal component, or essential dynamics, analysis of the positional fluctuations of the Cα atoms in the native state. Successful folding is obtained if the restraints are applied only to the eigenvectors with lowest eigenvalues, representing the most rigid quasi-constraint motions. If the essential eigenvectors, the ones accounting for most of the variance, are used, folding is not successful. These results clearly show that the eigenvectors with lowest eigenvalues contain the main mechanical information necessary to drive the folding process, while the essential eigenvectors represent the large concerted motions which can occur without folding/unfolding the protein.

Engineering Photocycle Dynamics: CRYSTAL STRUCTURES AND KINETICS OF THREE PHOTOACTIVE YELLOW PROTEIN HINGE-BENDING MUTANTS

Crystallographic and spectroscopic analyses of three hinge-bending mutants of the photoactive yellow protein are described. Previous studies have identified Gly(47) and Gly(51) as possible hinge points in the structure of the protein, allowing backbone segments around the chromophore to undergo large concerted motions. We have designed, crystallized, and solved the structures of three mutants: G47S, G51S, and G47S/G51S. The protein dynamics of these mutants are significantly affected. Transitions in the photocycle, measured with laser induced transient absorption spectroscopy, show rates up to 6-fold different from the wild type protein and show an additive effect in the double mutant. Compared with the native structure, no significant conformational differences were observed in the structures of the mutant proteins. We conclude that the structural and dynamic integrity of the region around these mutations is of crucial importance to the photocycle and suggest that the hinge-bending properties of Gly(51) may also play a role in PAS domain proteins where it is one of the few conserved residues.

Concerted motions in copper plastocyanin and azurin: an essential dynamics study

Essential dynamics analysis of molecular dynamics simulation trajectories (1.1 ns) of two copper containing electron transfer proteins, plastocyanin and azurin, has been performed. The protein essential modes have been analysed in order to identify large concerted motions which could be relevant for the electron transfer function exerted by these proteins. The analysis, conducted for temporal windows of different lengths along the protein trajectories, shows a rapid convergence and indicates that for both the proteins the predominant internal motions occur in a subspace of only a few degrees of freedom. Moreover, it is found that for both the proteins the likely binding sites (i.e. the hydrophobic and negative patches) with the reaction partners move in a concerted fashion with a few structural regions far from the active site. Such results are discussed in connection with the possible involvement of large concerted motions in the recognition and binding interaction with physiological electron transfer partners.

Functional concerted motions in the bovine serum retinol-binding protein.

The large concerted motions in the apo/holo bovine serum retinol-binding protein were studied using molecular dynamics simulation and 'essential dynamics' analysis. Initially, concerted motions were calculated from conformational differences between various crystal structures. The dynamic behaviour of the protein in the configurational space directions, described by these concerted motions, is analysed. This reveals that the large backbone dynamics of the protein is not influenced by the presence of retinol. Study of free retinol dynamics and retinol in the retinol binding site reveals that the protein binds retinol in a favourable conformation, as opposed to what has been previously described for the bovine cellular retinol-binding protein.

Essential spaces defined by NMR structure ensembles and molecular dynamics simulation show significant overlap

Large concerted motions of proteins which span its "essential space," are an important component of protein dynamics. We investigate to what extent structure ensembles generated with standard structure calculation techniques such as simulated annealing can capture these motions by comparing them to long-time molecular dynamics (MD) trajectories. The motions are analyzed by principal component analysis and compared using inner products of eigenvectors of the respective covariance matrices. Two very different systems are studied, the beta-spectrin PH domain and the single-stranded DNA binding protein (ssDBP) from the filamentous phage Pf3. A comparison of the ensembles from NMR and MD shows significant overlap of the essential spaces, which in the case of ssDBP is extraordinarily high. The influence of variations in the specifications of distance restraints is investigated. We also study the influence of the selection criterion for the final structure ensemble on the definition of mobility. The results suggest a modified criterion that improves conformational sampling in terms of amplitudes of correlated motion.

Essential spaces defined by NMR structure ensembles and molecular dynamics simulation show significant overlap

Large concerted motions of proteins which span its “essential space,” are an important component of protein dynamics. We investigate to what extent structure ensembles generated with standard structure calculation techniques such as simulated annealing can capture these motions by comparing them to long-time molecular dynamics (MD) trajectories. The motions are analyzed by principal component analysis and compared using inner products of eigenvectors of the respective covariance matrices. Two very different systems are studied, the β-spectrin PH domain and the single-stranded DNA binding protein (ssDBP) from the filamentous phage Pf3. A comparison of the ensembles from NMR and MD shows significant overlap of the essential spaces, which in the case of ssDBP is extraordinarily high. The influence of variations in the specifications of distance restraints is investigated. We also study the influence of the selection criterion for the final structure ensemble on the definition of mobility. The results suggest a modified criterion that improves conformational sampling in terms of amplitudes of correlated motion. Proteins 31:370–382, 1998. © 1998 Wiley-Liss, Inc.

Essential Dynamics from NMR Clusters: Dynamic Properties of the Myb DNA-Binding Domain and a Hinge-Bending Enhancing Variant

Application of the “essential dynamics” method to the NMR cluster of structures for the R2R3 DNA-binding domain of the mouse c-Myb transcriptional activator is described. Using this method, large concerted fluctuations of atoms are extracted showing a hinge-bending motion between the two (R2 and R3) Myb repeats on the basis of NMR data alone. Molecular dynamics simulation of the same protein allowed quantitative comparison of the large concerted motions calculated from experimental and theoretical data, showing a significant degree of similarity. Detailed inspection of the motions reveals a conserved proline that plays a key role in determining hinge flexibility. The proline-to-alanine mutation at this position, which has previously been characterized biochemically, was subjected to molecular dynamics and subsequent essential dynamics analysis. The hinge-bending motion between the two repeats was found to be enhanced for the mutant. The approach described should have general applications, predicting the effect of mutations on protein dynamic properties of other proteins.

Protein dynamics derived from clusters of crystal structures

A method is presented to mathematically extract concerted structural transitions in proteins from collections of crystal structures. The "essential dynamics" procedure is used to filter out small-amplitude fluctuations from such a set of structures; the remaining large conformational changes describe motions such as those important for the uptake/release of substrate/ligand and in catalytic reactions. The method is applied to sets of x-ray structures for a number of proteins, and the results are compared with the results from essential dynamics as applied to molecular dynamics simulations of those proteins. A significant degree of similarity is found, thereby providing a direct experimental basis for the application of such simulations to the description of large concerted motions in proteins.

A comparison of techniques for calculating protein essential dynamics

Recently the basic theory of essential dynamics, a method for extracting large concerted motions from protein molecular dynamics trajectories, was described. Here, we introduce and test new aspects. A method for diagonalizing large covariance matrices is presented. We show that it is possible to perform essential dynamics using different subsets of atoms and compare these to the basic C-α analysis. Essential dynamics analyses are also compared to the normal modes method. The stability of the essential space during a simulation is investigated by comparing the two halves of a trajectory. Apart from the analyses in Cartesian space, the essential dynamics in ϕ/ψ torsion angle space is discussed. © 1997 by John Wiley & Sons, Inc.

The consistency of large concerted motions in proteins in molecular dynamics simulations

A detailed investigation is presented into the effect of limited sampling time and small changes in the force field on molecular dynamics simulations of a protein. Thirteen independent simulations of the B1 IgG-binding domain of streptococcal protein G were performed, with small changes in the simulation parameters in each simulation. Parameters studied included temperature, bond constraints, cut-off radius for electrostatic interactions, and initial placement of hydrogen atoms. The essential dynamics technique was used to reveal dynamic differences between the simulations. Similar essential dynamics properties were found for all simulations, indicating that the large concerted motions found in the simulations are not particularly sensitive to small changes in the force field. A thorough investigation into the stability of the essential dynamics properties as derived from a molecular dynamics simulation of a few hundred picoseconds is provided. Although the definition of the essential modes of motion has not fully converged in these short simulations, the subspace in which these modes are confined is found to be reproducible.