Instantaneous Normal Modes and the Protein Glass Transition

In the instantaneous normal mode method, normal mode analysis is performed at instantaneous configurations of a condensed-phase system, leading to modes with negative eigenvalues. These negative modes provide a means of characterizing local anharmonicities of the potential energy surface. Here, we apply instantaneous normal mode to analyze temperature-dependent diffusive dynamics in molecular dynamics simulations of a small protein (a scorpion toxin). Those characteristics of the negative modes are determined that correlate with the dynamical (or glass) transition behavior of the protein, as manifested as an increase in the gradient with T of the average atomic mean-square displacement at ∼220 K. The number of negative eigenvalues shows no transition with temperature. Further, although filtering the negative modes to retain only those with eigenvectors corresponding to double-well potentials does reveal a transition in the hydration water, again, no transition in the protein is seen. However, additional filtering of the protein double-well modes, so as to retain only those that, on energy minimization, escape to different regions of configurational space, finally leads to clear protein dynamical transition behavior. Partial minimization of instantaneous configurations is also found to remove nondiffusive imaginary modes. In summary, examination of the form of negative instantaneous normal modes is shown to furnish a physical picture of local diffusive dynamics accompanying the protein glass transition.     

Coupling between Normal Modes Drives Protein Conformational Dynamics: Illustrations Using Allosteric Transitions in Myosin II

Structure-based elastic network models (ENMs) have been remarkably successful in describing conformational transitions in a variety of biological systems. Low-frequency normal modes are usually calculated from the ENM that characterizes elastic interactions between residues in contact in a given protein structure with a uniform force constant. To explore the dynamical effects of nonuniform elastic interactions, we calculate the robustness and coupling of the low-frequency modes in the presence of nonuniform variations in the ENM force constant. The variations in the elastic interactions, approximated here by Gaussian noise, approximately account for perturbation effects of heterogeneous residue-residue interactions or evolutionary sequence changes within a protein family. First-order perturbation theory provides an efficient and qualitatively correct estimate of the mode robustness and mode coupling for finite perturbations to the ENM force constant. The mode coupling analysis and the mode robustness analysis identify groups of strongly coupled modes that encode for protein functional motions. We illustrate the new concepts using myosin II motor protein as an example. The biological implications of mode coupling in tuning the allosteric couplings among the actin-binding site, the nucleotide-binding site, and the force-generating converter and lever arm in myosin isoforms are discussed. We evaluate the robustness of the correlation functions that quantify the allosteric couplings among these three key structural motifs.     

Thermodynamics of Conformational Transitions in a Disordered Protein Backbone Model

Conformational entropy is expected to contribute significantly to the thermodynamics of structural transitions in intrinsically disordered proteins or regions in response to protein/ligand binding, posttranslational modifications, and environmental changes. We calculated the backbone (dihedral) conformational entropy of oligoglycine $\left({\text{Gly}}_N\right)$, a protein backbone mimic and model intrinsically disordered region, as a function of chain length ($N=3$, 4, 5, 10, and 15) from simulations using three different approaches. The backbone conformational entropy scales linearly with chain length with a slope consistent with the entropy of folding of well-structured proteins. The entropic contributions of second-order dihedral correlations are predominantly through intraresidue ?-ψ pairs, suggesting that oligoglycine may be thermodynamically modeled as a system of independent glycine residues. We find the backbone conformational entropy to be largely independent of global structural parameters, like the end-to-end distance and radius of gyration. We introduce a framework referred to herein as “ensemble confinement” to estimate the loss (gain) of conformational free energy and its entropic component when individual residues are constrained to (released from) particular regions of the ?-ψ map. Quantitatively, we show that our protein backbone model resists ordering/folding with a significant, unfavorable ensemble confinement free energy because of the loss of a substantial portion of the absolute backbone entropy. Proteins can couple this free-energy reservoir to distal binding events as a regulatory mechanism to promote or suppress binding.     

A Comprehensive Multidimensional-Embedded, One-Dimensional Reaction Coordinate for Protein Unfolding/Folding

The goal of the Dynameomics project is to perform, store, and analyze molecular dynamics simulations of representative proteins, of all known globular folds, in their native state and along their unfolding pathways. To analyze unfolding simulations, the location of the protein along the unfolding reaction coordinate (RXN) must be determined. Properties such as the fraction of native contacts and radius of gyration are often used; however, there is an issue regarding degeneracy with these properties, as native and nonnative species can overlap. Here, we used 15 physical properties of the protein to construct a multidimensional-embedded, one-dimensional RXN coordinate that faithfully captures the complex nature of unfolding. The unfolding RXN coordinates for 188 proteins (1534 simulations and 22.9 μs in explicit water) were calculated. Native, transition, intermediate, and denatured states were readily identified with the use of this RXN coordinate. A global native ensemble based on the native-state properties of the 188 proteins was created. This ensemble was shown to be effective for calculating RXN coordinates for folds outside the initial 188 targets. These RXN coordinates enable, high-throughput assignment of conformational states, which represents an important step in comparing protein properties across fold space as well as characterizing the unfolding of individual proteins.     

Conformational Recognition of an Intrinsically Disordered Protein

There is a growing interest in understanding the properties of intrinsically disordered proteins (IDPs); however, the characterization of these states remains an open challenge. IDPs appear to have functional roles that diverge from those of folded proteins and revolve around their ability to act as hubs for protein-protein interactions. To gain a better understanding of the modes of binding of IDPs, we combined statistical mechanics, calorimetry, and NMR spectroscopy to investigate the recognition and binding of a fragment from the disordered protein Gab2 by the growth factor receptor-bound protein 2 (Grb2), a key interaction for normal cell signaling and cancer development. Structural ensemble refinement by NMR chemical shifts, thermodynamics measurements, and analysis of point mutations indicated that the population of preexisting bound conformations in the free-state ensemble of Gab2 is an essential determinant for recognition and binding by Grb2. A key role was found for transient polyproline II (PPII) structures and extended conformations. Our findings are likely to have very general implications for the biological behavior of IDPs in light of the evidence that a large fraction of these proteins possess a specific propensity to form PPII and to adopt conformations that are more extended than the typical random-coil states.     

Conformational Recognition of an Intrinsically Disordered Protein

There is a growing interest in understanding the properties of intrinsically disordered proteins (IDPs); however, the characterization of these states remains an open challenge. IDPs appear to have functional roles that diverge from those of folded proteins and revolve around their ability to act as hubs for protein-protein interactions. To gain a better understanding of the modes of binding of IDPs, we combined statistical mechanics, calorimetry, and NMR spectroscopy to investigate the recognition and binding of a fragment from the disordered protein Gab2 by the growth factor receptor-bound protein 2 (Grb2), a key interaction for normal cell signaling and cancer development. Structural ensemble refinement by NMR chemical shifts, thermodynamics measurements, and analysis of point mutations indicated that the population of preexisting bound conformations in the free-state ensemble of Gab2 is an essential determinant for recognition and binding by Grb2. A key role was found for transient polyproline II (PPII) structures and extended conformations. Our findings are likely to have very general implications for the biological behavior of IDPs in light of the evidence that a large fraction of these proteins possess a specific propensity to form PPII and to adopt conformations that are more extended than the typical random-coil states.     

Conformational Equilibrium of CDK/Cyclin Complexes by Molecular Dynamics with Excited Normal Modes

Cyclin-dependent kinases (CDKs) and their associated regulatory cyclins are central for timely regulation of cell-cycle progression. They constitute attractive pharmacological targets for development of anticancer therapeutics, since they are frequently deregulated in human cancers and contribute to sustained, uncontrolled tumor proliferation. Characterization of their structural/dynamic features is essential to gain in-depth insight into structure-activity relationships. In addition, the identification of druggable pockets or key intermediate conformations yields potential targets for the development of novel classes of inhibitors. Structural studies of CDK2/cyclin A have provided a wealth of information concerning monomeric/heterodimeric forms of this kinase. There is, however, much less structural information for other CDK/cyclin complexes, including CDK4/cyclin D1, which displays an alternative (open) position of the cyclin partner relative to CDK, contrasting with the closed CDK2/cyclin A conformation. In this study, we carried out normal-mode analysis and enhanced sampling simulations with our recently developed method, molecular dynamics with excited normal modes, to understand the conformational equilibrium on these complexes. Interestingly, the lowest-frequency normal mode computed for each complex described the transition between the open and closed conformations. Exploration of these motions with an explicit-solvent representation using molecular dynamics with excited normal modes confirmed that the closed conformation is the most stable for the CDK2/cyclin A complex, in agreement with their experimentally available structures. On the other hand, we clearly show that an open↔closed equilibrium may exist in CDK4/cyclin D1, with closed conformations resembling that captured for CDK2/cyclin A. Such conformational preferences may result from the distinct distributions of frustrated contacts in each complex. Using the same approach, the putative roles of the Thr¹⁶⁰ phosphoryl group and the T-loop conformation were investigated. These results provide a dynamic view of CDKs revealing intermediate conformations not yet characterized for CDK members other than CDK2, which will be useful for the design of inhibitors targeting critical conformational transitions.     

DNA and Protein Requirements for Substrate Conformational Changes Necessary for Human Flap Endonuclease-1-catalyzed Reaction

Human flap endonuclease-1 (hFEN1) catalyzes the essential removal of single-stranded flaps arising at DNA junctions during replication and repair processes. hFEN1 biological function must be precisely controlled, and consequently, the protein relies on a combination of protein and substrate conformational changes as a prerequisite for reaction. These include substrate bending at the duplex-duplex junction and transfer of unpaired reacting duplex end into the active site. When present, 5′-flaps are thought to thread under the helical cap, limiting reaction to flaps with free 5′-termini in vivo. Here we monitored DNA bending by FRET and DNA unpairing using 2-aminopurine exciton pair CD to determine the DNA and protein requirements for these substrate conformational changes. Binding of DNA to hFEN1 in a bent conformation occurred independently of 5′-flap accommodation and did not require active site metal ions or the presence of conserved active site residues. More stringent requirements exist for transfer of the substrate to the active site. Placement of the scissile phosphate diester in the active site required the presence of divalent metal ions, a free 5′-flap (if present), a Watson-Crick base pair at the terminus of the reacting duplex, and the intact secondary structure of the enzyme helical cap. Optimal positioning of the scissile phosphate additionally required active site conserved residues Tyr⁴⁰, Asp¹⁸¹, and Arg¹⁰⁰ and a reacting duplex 5′-phosphate. These studies suggest a FEN1 reaction mechanism where junctions are bound and 5′-flaps are threaded (when present), and finally the substrate is transferred onto active site metals initiating cleavage.     

Elastic Interactions of Photosynthetic Reaction Center Proteins Affecting Phase Transitions and Protein Distributions

Reaction-center proteins of Rhodopseudomonas Sphaeroides reconstituted into phosphatidylcholine vesicles shift and broaden the fluid-gel transition of the lipid bilayer. The amount of broadening and temperature shift of the transition depend both on protein concentration and on lipid chain length. In particular, the direction of the transition temperature shift is very sensitive to lipid chain length. Electron micrographs show homogeneous protein distribution on the fluid surface whereas the solid phase contains protein aggregates the type depending on chain length. The results can qualitatively be understood in the framework of a mattress model of lipid/protein interactions in membranes.     

Conformational Transitions in Yeast Chorismate Mutase Important for Allosteric Regulation as Identified by Nuclear Magnetic Resonance Spectroscopy

Proteins fluctuate between different conformations in solution, and these conformational fluctuations can be important for protein function and allosteric regulation. The chorismate mutase from Saccharomyces cerevisiae (ScCM), a key enzyme in the biosynthesis of aromatic amino acids, is allosterically activated and inhibited by tryptophan and tyrosine, respectively. It was initially proposed that in the absence of effector, ScCM fluctuates between activated R and inhibited T conformations according to the Monod-Wyman-Changeux (MWC) model, although a more complex regulation pattern was later suggested by mutagenesis and kinetic data. Here we used NMR relaxation dispersion experiments to understand the conformational fluctuations on the microsecond-to-millisecond timescale that occur in ScCM. In the absence of allosteric effectors, ScCM did not exclusively exchange between T and R conformations, suggesting that the two-state MWC model is insufficient to explain conformational dynamics. Addition of tyrosine led to the quenching of much of the motion on this timescale, while new motions were identified in the presence of tryptophan. These new motions are consistent with conformational fluctuations into an alternative conformation that may be important for enzyme activity.     

Conformational Transitions of Poly(dl-dC) in Aqueous Solution as Studied by Ultraviolet Resonance Raman Spectroscopy

Abstract

Poly(dI-dC) in aqueous solution can undergo different equilibrium geometries, which strongly depend on salt nature and concentrations. These equilibrium structures have been monitored by resonance Raman spectroscopy (RRS) measurements in the ultraviolet region, i. e. by using 257 and 281 nm laser excitation wavelengths which favor the resonance enhancement of the Raman contributions from inosine and cytosine residues of poly(dI-dC), respectively. Spectral changes depending on the NaCl concentration and on the presence of Ni²⁺ ions have been observed and interpreted in comparison with RRS results previously obtained for other alternating purine-pyrimidine polydeoxyribonucleotides, i.e. poly(dG-dC), poly(dA- dT) and poly(dA-dC). poly(dG-dT), which also showed B to Z conformational transitions in varying the salt concentrations. It is shown here that: i) the base stacking geometries are nearly the same in the high-salt form (5 M NaCl) of poly(dl-dC) as in the low-salt form (0.1M NaCl) of the polymer, ii) however, the high-salt structure yields important differences from a B-helix (obtained in low-salt solution) as regards the nucleoside conformations (sugar puckering and base-sugar orientation), and: iii) the addition of 9 mM NiCl₂ in the high-salt (5 M NaCl) solution of poly(dI-dC) induces the Z-conformation of the polymer.     

Structural Basis for Conformational Dynamics of GTP-bound Ras Protein

Ras family small GTPases assume two interconverting conformations, “inactive” state 1 and “active” state 2, in their GTP-bound forms. Here, to clarify the mechanism of state transition, we have carried out x-ray crystal structure analyses of a series of mutant H-Ras and M-Ras in complex with guanosine 5′-(β,γ-imido)triphosphate (GppNHp), representing various intermediate states of the transition. Crystallization of H-RasT35S-GppNHp enables us to solve the first complete tertiary structure of H-Ras state 1 possessing two surface pockets unseen in the state 2 or H-Ras-GDP structure. Moreover, determination of the two distinct crystal structures of H-RasT35S-GppNHp, showing prominent polysterism in the switch I and switch II regions, reveals a pivotal role of the guanine nucleotide-mediated interaction between the two switch regions and its rearrangement by a nucleotide positional change in the state 2 to state 1 transition. Furthermore, the ³¹P NMR spectra and crystal structures of the GppNHp-bound forms of M-Ras mutants, carrying various H-Ras-type amino acid substitutions, also reveal the existence of a surface pocket in state 1 and support a similar mechanism based on the nucleotide-mediated interaction and its rearrangement in the state 1 to state 2 transition. Intriguingly, the conformational changes accompanying the state transition mimic those that occurred upon GDP/GTP exchange, indicating a common mechanistic basis inherent in the high flexibility of the switch regions. Collectively, these results clarify the structural features distinguishing the two states and provide new insights into the molecular basis for the state transition of Ras protein.