Evgenii L. Kovrigin

The mechanism of rate-limiting motions in enzyme function

The ability to use conformational flexibility is a hallmark of enzyme function. Here we show that protein motions and catalytic activity in a RNase are coupled and display identical solvent isotope effects. Solution NMR relaxation experiments identify a cluster of residues, some distant from the active site, that are integral to this motion. These studies implicate a single residue, histidine-48, as the key modulator in coupling protein motion with enzyme function. Mutation of H48 to alanine results in loss of protein motion in the isotope-sensitive region of the enzyme. In addition, kcat decreases for this mutant and the kinetic solvent isotope effect on kcat, which was 2.0 in WT, is near unity in H48A. Despite being located 18 Å from the enzyme active site, H48 is essential in coordinating the motions involved in the rate-limiting enzymatic step. These studies have identified, of ≈160 potential exchangeable protons, a single site that is integral in the rate-limiting step in RNase A enzyme function.

Global conformational dynamics in ras.

Ras and its homologues are central to regulation of a multitude of cellular processes. Ras in complex with GTP binds and activates its downstream signaling partners. (31)P NMR studies indicated that the Ras-GTP conformation is heterogeneous on a millisecond time scale, but details of its conformational dynamics remain unknown. Here we present evidence that the conformational exchange process in human H-Ras complexed with GTP mimic GppNHp is global, encompassing most of the GTPase catalytic domain. The correlated character of conformational dynamics in Ras opens opportunities for understanding allosteric effects in Ras function.

NMR line shapes and multi-state binding equilibria

Biological function of proteins relies on conformational transitions and binding of specific ligands. Protein–ligand interactions are thermodynamically and kinetically coupled to conformational changes in protein structures as conceptualized by the models of pre-existing equilibria and induced fit. NMR spectroscopy is particularly sensitive to complex ligand-binding modes—NMR line-shape analysis can provide for thermodynamic and kinetic constants of ligand-binding equilibria with the site-specific resolution. However, broad use of line shape analysis is hampered by complexity of NMR line shapes in multi-state systems. To facilitate interpretation of such spectral patterns, I computationally explored systems where isomerization or dimerization of a protein (receptor) molecule is coupled to binding of a ligand. Through an extensive analysis of multiple exchange regimes for a family of three-state models, I identified signature features to guide an NMR experimentalist in recognizing specific interaction mechanisms. Results show that distinct multi-state models may produce very similar spectral patterns. I also discussed aggregation of a receptor as a possible source of spurious three-state line shapes and provided specific suggestions for complementary experiments that can ensure reliable mechanistic insight.

On the stabilizing action of protein denaturants: acetonitrile effect on stability of lysozyme in aqueous solutions.

Stability of hen lysozyme in the presence of acetonitrile (MeCN) at different pH values of the medium was studied by scanning microcalorimetry with a special emphasis on determination of reliable values of the denaturational heat capacity change. It was found that the temperature of denaturation decreases on addition of MeCN. However, the free energy extrapolation showed that below room temperature the thermodynamic stability increases at low concentrations of MeCN in spite of the general destabilizing effect at higher concentrations and temperatures. Charge-induced contribution to this stabilization was shown to be negligible (no pH-dependence was found); therefore, the most probable cause for the phenomenon is an increase of hydrophobic interactions at low temperatures in aqueous solutions containing small amounts of the organic additive. The difference in preferential solvation of native and denatured states of lysozyme was calculated from the stabilization free energy data. It was found that the change in preferential solvation strongly depends on the temperature in the water-rich region. At the higher MeCN content this dependence decreases until, at 0.06 mole fractions of MeCN, the difference in the preferential solvation between native and denatured lysozyme becomes independent of the temperature over a range of 60 K. The importance of taking into account non-ideality of a mixed solution, when analyzing preferential solvation phenomena was emphasized.

Alteration of hydrogen bonding in the vicinity of histidine 48 disrupts millisecond motions in RNase A.

The motion of amino acid residues on the millisecond (ms) time scale is involved in the tight regulation of catalytic function in numerous enzyme systems. Using a combination of mutational, enzymological, and relaxation-compensated (15)N Carr-Purcell-Meiboom-Gill (CPMG) methods, we have previously established the conformational significance of the distant His48 residue and the neighboring loop 1 in RNase A function. These studies suggested that RNase A relies on an intricate network of hydrogen bonding interactions involved in propagating functionally relevant, long-range ms motions to the catalytic site of the enzyme. To further investigate the dynamic importance of this H-bonding network, this study focuses on the individual replacement of Thr17 and Thr82 with alanine, effectively altering the key H-bonding interactions that connect loop 1 and His48 to the rest of the protein. (15)N CPMG dispersion studies, nuclear magnetic resonance (NMR) chemical shift analysis, and NMR line shape analysis of point mutants T17A and T82A demonstrate that the evolutionarily conserved single H-bond linking His48 to Thr82 is essential for propagating ms motions from His48 to the active site of RNase A on the time scale of catalytic turnover, whereas the T17A mutation increases the off rate and conformational exchange motions in loop 1. Accumulating evidence from our mutational studies indicates that residues experiencing conformational exchange in RNase A can be grouped into two separate clusters displaying distinct dynamical features, which appear to be independently affected by mutation. Overall, this study illuminates how tightly controlled and finely tuned ms motions are in RNase A, suggesting that designed modulation of protein motions may be possible.

Folding Under Inequilibrium Conditions as a Possible Reason for Partial Irreversibility of Heat-Denatured Proteins: Computer Simulation Study

Using computer simulations we have studied possible effects of heating and cooling at different scan rates on unfolding and refolding of macromolecules. We have shown that even the simplest two-state reversible transition can behave irreversibly when an unfavorable combination of cooling rate, relaxation time and activation energy of refolding occurs. On the basis of this finding we suppose that apparent irreversibility of some proteins denatured by heat may result from slow relaxation on cooling rather than thermodynamic instability and/or irreversible alterations of the polypeptide chain. Using this kinetic reversible two-state model, we estimated the effects of the scan rate and kinetic parameters of the macromolecule on its unfolding-refolding process. A few recommendations are suggested on how to reach maximal possible recovery after denaturation if refolding appears to be under kinetic control.

Conservation of Flexible Residue Clusters among Structural and Functional Enzyme Homologues

Background: It remains unclear whether structural homologues rely on similar concerted motions to promote enzyme function. Results: Ribonuclease homologues display similar, contiguous clustering motions that can be modulated by mutagenesis. Conclusion: Conformational flexibility can be conserved between distant structural homologues. Significance: Controlling dynamics to modulate function has broad implications in protein engineering and allosteric drug design. Conformational flexibility between structural ensembles is an essential component of enzyme function. Although the broad dynamical landscape of proteins is known to promote a number of functional events on multiple time scales, it is yet unknown whether structural and functional enzyme homologues rely on the same concerted residue motions to perform their catalytic function. It is hypothesized that networks of contiguous and flexible residue motions occurring on the biologically relevant millisecond time scale evolved to promote and/or preserve optimal enzyme catalysis. In this study, we use a combination of NMR relaxation dispersion, model-free analysis, and ligand titration experiments to successfully capture and compare the role of conformational flexibility between two structural homologues of the pancreatic ribonuclease family: RNase A and eosinophil cationic protein (or RNase 3). In addition to conserving the same catalytic residues and structural fold, both homologues show similar yet functionally distinct clusters of millisecond dynamics, suggesting that conformational flexibility can be conserved among analogous protein folds displaying low sequence identity. Our work shows that the reduced conformational flexibility of eosinophil cationic protein can be dynamically and functionally reproduced in the RNase A scaffold upon creation of a chimeric hybrid between the two proteins. These results support the hypothesis that conformational flexibility is partly required for catalytic function in homologous enzyme folds, further highlighting the importance of dynamic residue sectors in the structural organization of proteins.

Complete Thermodynamic and Kinetic Characterization of the Isomer-Specific Interaction between Pin1-WW Domain and the Amyloid Precursor Protein Cytoplasmic Tail Phosphorylated at Thr668

Peptidyl prolyl cis-trans isomerization acts as an effective molecular timer that plays significant roles in biological and pathological processes. Enzymes such as Pin1 catalyze cis-trans isomerization, accelerating the otherwise slow isomerization rate into time scales relevant for cellular signaling. Here we have combined NMR line shape analysis, fluorescence spectroscopy, and isothermal titration calorimetry to determine the kinetic and thermodynamic parameters describing the trans-specific interaction between the binding domain of Pin1 (WW domain) and a key cis-trans molecular switch in the amyloid precursor protein cytoplasmic tail. A three-state model, in which the cis-trans isomerization equilibrium is coupled to the binding equilibrium through the trans isomer, was found to fit the data well. The trans isomer binds the WW domain with ∼22 μM affinity via very fast association (approaching the diffusion limit) and dissociation rates. The common structural and electrostatic characteristics of Pin1 substrates, which contain a phosphorylated serine/threonine-proline motif, suggest that very rapid binding kinetics are a general feature of Pin1 interactions with other substrates. The fast binding kinetics of the WW domain allows rapid response of Pin1 to the dynamic events of phosphorylation and dephosphorylation in the cell that alter the relative populations of diverse Pin1 substrates. Furthermore, our results also highlight the vastly different rates at which slow uncatalyzed cis-trans isomerization and fast isomer-specific binding events occur. These results, along with the experimental methods presented herein, should guide future experiments aimed at the thermodynamic and kinetic characterization of cis-trans molecular switches and isomer-specific interactions involved in various biological processes.

The Ras G Domain Lacks the Intrinsic Propensity to Form Dimers

Ras GTPase is a molecular switch controlling a number of cellular pathways including growth, proliferation, differentiation, and apoptosis. Recent reports indicated that Ras undergoes dimerization at the membrane surface through protein-protein interactions. If firmly established this property of Ras would require profound reassessment of a large amount of published data and modification of the Ras signaling paradigm. One proposed mechanism of dimerization involves formation of salt bridges between the two GTPase domains (G domains) leading to formation of a compact dimer as observed in Ras crystal structures. In this work, we interrogated the intrinsic ability of Ras to self-associate in solution by creating conditions of high local concentration through irreversibly tethering the two G domains together at their unstructured C-terminal tails. We evaluated possible self-association in this inverted tandem conjugate via analysis of the time-domain fluorescence anisotropy and NMR chemical shift perturbations. We did not observe the increased rotational correlation time expected for the G domain dimer. Variation of the ionic strength (to modulate stability of the salt bridges) did not affect the rotational correlation time in the tandem further supporting independent rotational diffusion of two G domains. In a parallel line of experiments to detect and map weak self-association of the G domains, we analyzed NMR chemical shifts perturbations at a number of sites near the crystallographic dimer interface. The nearly complete lack of chemical shift perturbations in the tandem construct supported a simple model with the independent G domains repelled from each other by their overall negative charge. These results lead us to the conclusion that self-association of the G domains cannot be responsible for homodimerization of Ras reported in the literature.

Characterization of the second ion-binding site in the G domain of H-Ras.

Ras is a small monomeric GTPase acting as molecular switch in multiple cellular processes. The N-terminal G domain of Ras binds GTP or GDP accompanied by a magnesium ion, which is strictly required for GTPase activity and performs a structural role. Another ion-binding site on the opposite face of the G domain has been recently observed to specifically associate with calcium acetate in the crystal [Buhrman, G., et al. (2010) Proc. Natl. Aacd. Sci. U.S.A. 107, 4931-4936]. In this article, we report thermodynamic measurements of the affinity and specificity of the remote ion-binding site in H-Ras as observed in solution. Using (15)N-(1)H nuclear magnetic resonance spectroscopy, we determined that, in contrast to the crystalline state, the remote site in solution is specific for a divalent cation, binding both calcium and magnesium with anions playing a minimal role. The affinity of the remote site for divalent cations is in the low millimolar range and remarkably different for GDP- and GppNHp-bound forms of the G domain, indicating that the GTP-binding pocket and the remote site are allosterically coupled through the distance of more than 25 Å. Considering that the remote site is oriented toward the membrane surface in vivo, we hypothesize that its cognate biological ligand might be a positively charged group extending from a lipid or an integral membrane protein.