Thermal stability of glucokinase (GK) as influenced by the substrate glucose, an allosteric glucokinase activator drug (GKA) and the osmolytes glycerol and urea.

We investigated how glycerol, urea, glucose and a GKA influence kinetics and stability of wild-type and mutant GK. Glycerol and glucose stabilized GK additively. Glycerol barely affected the TF spectra of all GKs but decreased k(cat), glucose S(0.5) and K(D) values and ATP K(M) while leaving cooperativity unchanged. Glycerol sensitized all GKs to GKA as shown by TF. Glucose increased TF of GKs without influence of glycerol on the effect. Glycerol and GKA affected kinetics and binding additively. The activation energies for thermal denaturation of GK were a function of glucose with K(D)s of 3 and 1mM without and with glycerol, respectively. High urea denatured wild type GK reversibly at 20 and 60°C and urea treatment of irreversibly heat denatured GK allowed refolding as demonstrated by TF including glucose response. We concluded: Glycerol stabilizes GK indirectly without changing the folding structure of the apoenzyme, by restructuring the surface water of the protein, whereas glucose stabilizes GK directly by binding to its substrate site and inducing a compact conformation. Glucose or glycerol (alone or combined) is unable to prevent irreversible heat denaturation above 40°C. However, urea denatures GK reversibly even at 60°C by binding to the protein backbone and directly interacting with hydrophobic side chains. It prevents irreversible aggregation allowing complete refolding when urea is removed. This study establishes the foundation for exploring numerous instability mutants among the more than 600 variant GKs causing diabetes in animals and humans.

Main chain and side chain dynamics of oxidized flavodoxin from Cyanobacterium anabaena.

Oxidized flavodoxin from Cyanobacterium anabaena PCC 7119 is used as a model system to investigate the fast internal dynamics of a flavin-bearing protein. Virtually complete backbone and side chain resonance NMR assignments of an oxidized flavodoxin point mutant (C55A) have been determined. Backbone and side chain dynamics in flavodoxin (C55A) were investigated using (15)N amide and deuterium methyl NMR relaxation methods. The squared generalized order parameters (S(NH)(2)) for backbone amide N-H bonds are found to be uniformly high ( approximately 0.923 over 109 residues in regular secondary structure), indicating considerable restriction of motion in the backbone of the protein. In contrast, methyl-bearing side chains are considerably heterogeneous in their amplitude of motion, as indicated by obtained symmetry axis squared generalized order parameters (S(axis)(2)). However, in comparison to nonprosthetic group-bearing proteins studied with these NMR relaxation methods, the side chains of oxidized flavodoxin are unusually rigid.

Oligomerization of the integrin alphaIIbbeta3: roles of the transmembrane and cytoplasmic domains.

Integrins are a family of alpha/beta heterodimeric membrane proteins, which mediate cell-cell and cell-matrix interactions. The molecular mechanisms by which integrins are activated and cluster are currently poorly understood. One hypothesis posits that the cytoplasmic tails of the alpha and beta subunits interact strongly with one another in a 1:1 interaction, and that this interaction is modulated in the course of the activation of alphaIIbbeta3 [Hughes, P. E., et al. (1996) J. Biol. Chem. 271, 6571-6574]. To examine the structural basis for this interaction, protein fragments encompassing the transmembrane helix plus cytoplasmic tails of the alpha and beta subunits of alphaIIbbeta3 were expressed and studied in phospholipid micelles at physiological salt concentrations. Analyses of these fragments by analytical ultracentrifugation, NMR, circular dichroism, and electrophoresis indicated that they had very little or no tendency to interact with one another. Instead, they formed homomeric interactions, with the alpha- and beta-fragments forming dimers and trimers, respectively. Thus, these regions of the protein structure may contribute to the clustering of integrins that accompanies cellular adhesion.

Water magnetic relaxation dispersion in biological systems: the contribution of proton exchange and implications for the noninvasive detection of cartilage degradation.

Magnetic relaxation has been used extensively to study and characterize biological tissues. In particular, spin-lattice relaxation in the rotating frame (T(1rho)) of water in protein solutions has been demonstrated to be sensitive to macromolecular weight and composition. However, the nature of the contribution from low frequency processes to water relaxation remains unclear. We have examined this problem by studying the water T(1rho) dispersion in peptide solutions ((14)N- and (15)N-labeled), glycosaminoglycan solutions, and samples of bovine articular cartilage before and after proteoglycan degradation. We find in model systems and tissue that hydrogen exchange from NH and OH groups to water dominates the low frequency water T(1rho) dispersion, in the context of the model used to interpret the relaxation data. Further, low frequency dispersion changes are correlated with loss of proteoglycan from the extra-cellular matrix of articular cartilage. This finding has significance for the noninvasive detection of matrix degradation.

Dynamic activation of protein function: a view emerging from NMR spectroscopy.

Recent developments in solution NMR methods have allowed for an unprecedented view of protein dynamics. Current insights into the nature of protein dynamics and their potential influence on protein structure, stability and function are reviewed. Particular emphasis is placed on the potential of fast side chain motion to report on the residual conformational entropy of proteins and how this entropy can enter into both the thermodynamic and kinetic aspects of protein function.

Dynamics of a de novo designed three-helix bundle protein studied by 15N, 13C, and 2H NMR relaxation methods.

Understanding how the amino acid sequence of a polypeptide chain specifies a unique, functional three-dimensional structure remains an important goal, especially in the context of the emerging discipline of de novo protein design. Alpha3D is a single chain protein of 73 amino acids resulting from a de novo design effort. Previous solution nuclear magnetic resonance studies of alpha3D confirm that the protein adopts the designed structure of a three-helix bundle. Furthermore, alpha3D has been previously shown to possess all of the major thermodynamic and structural characteristics of natural proteins, though it shares no sequence homology to any protein sequence in the database. In this work, the backbone and side-chain dynamics of alpha3D were investigated using 15N, 13C, and 2H nuclear magnetic resonance relaxation methods with the aim of assessing the character of the internal motions of this native-like protein of de novo design. At the backbone level, both 15N and 13C(alpha) relaxation studies indicate highly restrictive motion on the picosecond to nanosecond time scale in the alpha-helical regions of alpha3D, with increasing mobility at the ends of the alpha-helices and in the two loop regions. This is largely consistent with what is seen in proteins of natural origin. Overall, the view provided by both 2H and 13C methyl relaxation methods suggest that the side chains of alpha3D are more dynamic compared to natural proteins. Regions of relative flexibility bound clusters of rigid methyl-bearing side-chain groups that are interspersed with aromatic and beta-branched amino acids. The time scale of motions associated with methyl-bearing side chains of alpha3D are significantly longer than that seen in natural proteins. These results indicate that the strategies underlying the design of alpha3D have largely, but not completely, captured both the structural and dynamic character of natural proteins.

Hydrophobic modulation of heme properties in heme protein maquettes.

We have investigated the properties of the two hemes bound to histidine in the H10 positions of the uniquely structured apo form of the heme binding four-helix bundle protein maquette [H10H24-L6I,L13F](2), here called [I(6)F(13)H(24)](2) for the amino acids at positions 6 (I), 13 (F) and 24 (H), respectively. The primary structure of each alpha-helix, alpha-SH, in [I(6)F(13)H(24)](2) is Ac-CGGGEI(6)WKL.H(10)EEF(13)LKK.FEELLKL.H(24)EERLKK.L-CONH(2). In our nomenclature, [I(6)F(13)H(24)] represents the disulfide-bridged di-alpha-helical homodimer of this sequence, i.e., (alpha-SS-alpha), and [I(6)F(13)H(24)](2) represents the dimeric four helix bundle composed of two di-alpha-helical subunits, i.e., (alpha-SS-alpha)(2). We replaced the histidines at positions H24 in [I(6)F(13)H(24)](2) with hydrophobic amino acids incompetent for heme ligation. These maquette variants, [I(6)F(13)I(24)](2), [I(6)F(13)A(24)](2), and [I(6)F(13)F(24)](2), are distinguished from the tetraheme binding parent peptide, [I(6)F(13)H(24)](2), by a reduction in the heme:four-helix bundle stoichiometry from 4:1 to 2:1. Iterative redesign has identified phenylalanine as the optimal amino acid replacement for H24 in the context of apo state conformational specificity. Furthermore, the novel second generation diheme [I(6)F(13)F(24)](2) maquette was related to the first generation diheme [H10A24](2) prototype, [L(6)L(13)A(24)](2) in the present nomenclature, via a sequential path in sequence space to evaluate the effects of conservative hydrophobic amino acid changes on heme properties. Each of the disulfide-linked dipeptides studied was highly helical (>77% as determined from circular dichroism spectroscopy), self-associates in solution to form a dimer (as determined by size exclusion chromatography), is thermodynamically stable (-DeltaG(H)2(O) >18 kcal/mol), and possesses conformational specificity that NMR data indicate can vary from multistructured to single structured. Each peptide binds one heme with a dissociation constant, K(d1) value, tighter than 65 nM forming a series of monoheme maquettes. Addition of a second equivalent of heme results in heme binding with a K(d2) in the range of 35-800 nM forming the diheme maquette state. Single conservative amino acid changes between peptide sequences are responsible for up to 10-fold changes in K(d) values. The equilibrium reduction midpoint potential (E(m7.5)) determined in the monoheme state ranges from -156 to -210 mV vs SHE and in the diheme state ranges from -144 to -288 mV. An observed heme-heme electrostatic interaction (>70 mV) in the diheme state indicates a syn global topology of the di-alpha-helical monomers. The heme affinity and electrochemistry of the three H24 variants studied identify the tight binding sites (K(d1) and K(d2) values <200 nM) having the lower reduction midpoint potentials (E(m7.5) values of -155 and -260 mV) with the H10 bound hemes in the parent tetraheme state of [H10H24-L6I,L13F](2), here called [I(6)F(13)H(24)](2). The results of this study illustrate that conservative hydrophobic amino acid changes near the heme binding site can modulate the E(m) by up to +/-50 mV and the K(d) by an order of magnitude. Furthermore, the effects of multiple single amino acid changes on E(m) and K(d) do not appear to be additive.

Main chain and side chain dynamics of a heme protein: 15N and 2H NMR relaxation studies of R. capsulatus ferrocytochrome c2.

A detailed characterization of the main chain and side chain dynamics in R. capsulatus ferrocytochrome c(2) derived from (2)H NMR relaxation of methyl group resonances is presented. (15)N relaxation measurements confirm earlier results indicating that R. capsulatus ferrocytochrome c(2) exhibits minor rotational anisotropy in solution. The current study is focused on the use of deuterium relaxation in side chain methyl groups, which has been shown to provide a detailed and accurate measure of internal dynamics. Results obtained indicate that the side chains of ferrocytochrome c(2) exhibit a wide range of motional amplitudes, but are more rigid than generally found in the interior of nonprosthetic group bearing globular proteins. This unusual rigidity is ascribed to the interactions of the protein with the large heme prosthetic group. This observation has significant implications for the potential of the heme-protein interface to modulate the redox properties of the protein and also points to the need for great precision in the design and engineering of heme proteins.

Microscopic origins of entropy, heat capacity and the glass transition in proteins.

Internal motion is central to protein folding, to protein stability through the resulting residual entropy, and to protein function. Despite its importance, the precise nature of the internal motions of protein macromolecules remains a mystery. Here we report a survey of the temperature dependence of the fast dynamics of methyl-bearing side chains in a calmodulin-peptide complex using site-specific deuterium NMR relaxation methods. The amplitudes of motion had a markedly heterogeneous spectrum and segregated into three largely distinct classes. Other proteins studied at single temperatures tend to segregate similarly. Furthermore, a large variability in the degree of thermal activation of the dynamics in the calmodulin complex indicates a heterogeneous distribution of residual entropy and hence reveals the microscopic origins of heat capacity in proteins. These observations also point to an unexpected explanation for the low-temperature 'glass transition' of proteins. It is this transition that has been ascribed to the creation of protein motional modes that are responsible for biological activity.

Redistribution and loss of side chain entropy upon formation of a calmodulin-peptide complex.

The response of the internal dynamics of calcium-saturated calmodulin to the formation of a complex with a peptide model of the calmodulin-binding domain of the smooth muscle myosin light chain kinase has been studied using NMR relaxation methods. The backbone of calmodulin is found to be unaffected by the binding of the domain, whereas the dynamics of side chains are significantly perturbed. The changes in dynamics are interpreted in terms of a heterogeneous partitioning between structure (enthalpy) and dynamics (entropy). These data provide a microscopic view of the residual entropy of a protein in two functional states and suggest extensive enthalpy/entropy exchange during the formation of a protein-protein interface.

The binding site for UCH-L3 on ubiquitin: mutagenesis and NMR studies on the complex between ubiquitin and UCH-L3.

The ubiquitin fold is a versatile and widely used targeting signal that is added post-translationally to a variety of proteins. Covalent attachment of one or more ubiquitin domains results in localization of the target protein to the proteasome, the nucleus, the cytoskeleton or the endocytotic machinery. Recognition of the ubiquitin domain by a variety of enzymes and receptors is vital to the targeting function of ubiquitin. Several parallel pathways exist and these must be able to distinguish among ubiquitin, several different types of polymeric ubiquitin, and the various ubiquitin-like domains. Here we report the first molecular description of the binding site on ubiquitin for ubiquitin C-terminal hydrolase L3 (UCH-L3). The site on ubiquitin was experimentally determined using solution NMR, and site-directed mutagenesis. The site on UCH-L3 was modeled based on X-ray crystallography, multiple sequence alignments, and computer-aided docking. Basic residues located on ubiquitin (K6, K11, R72, and R74) are postulated to contact acidic residues on UCH-L3 (E10, E14, D33, E219). These putative interactions are testable and fully explain the selectivity of ubiquitin domain binding to this enzyme.

De novo proteins as models of radical enzymes.

Catalytically essential side-chain radicals have been recognized in a growing number of redox enzymes. Here we present a novel approach to study this class of redox cofactors. Our aim is to construct a de novo protein, a radical maquette, that will provide a protein framework in which to investigate how side-chain radicals are generated, controlled, and directed toward catalysis. A tryptophan and a tyrosine radical maquette, denoted alpha(3)W(1) and alpha(3)Y(1), respectively, have been synthesized. alpha(3)W(1) and alpha(3)Y(1) contain 65 residues each and have molecular masses of 7.4 kDa. The proteins differ only in residue 32, which is the position of their single aromatic side chain. Structural characterization reveals that the proteins fold in water solution into thermodynamically stable, alpha-helical conformations with well-defined tertiary structures. The proteins are resistant to pH changes and remain stable through the physiological pH range. The aromatic residues are shown to be located within the protein interior and shielded from the bulk phase, as designed. Differential pulse voltammetry was used to examine the reduction potentials of the aromatic side chains in alpha(3)W(1) and alpha(3)Y(1) and compare them to the potentials of tryptophan and tyrosine when dissolved in water. The tryptophan and tyrosine potentials were raised considerably when moved from a solution environment to a well-ordered protein milieu. We propose that the increase in reduction potential of the aromatic residues buried within the protein, relative to the solution potentials, is due to a lack of an effective protonic contact between the aromatic residues and the bulk solution.

Preparation of encapsulated proteins dissolved in low viscosity fluids.

The majority of proteins are too large to be comprehensively examined by solution NMR methods, primarily because they tumble too slowly in solution. One potential approach to making the NMR relaxation properties of large proteins amenable to modern solution NMR techniques is to encapsulate them in a reverse micelle which is dissolved in a low viscosity fluid. Unfortunately, promising low viscosity fluids such as the short chain alkanes, supercritical carbon dioxide, and various halocarbon refrigerants all require the application of significant pressure to be kept liquefied at room temperature. Here we describe the design and use of a simple cost effective NMR tube suitable for the preparation of solutions of proteins encapsulated in reverse micelles dissolved in such fluids.

Assessing potential bias in the determination of rotational correlation times of proteins by NMR relaxation.

The various factors that influence the reliable and efficient determination of the correlation time describing molecular reorientation of proteins by NMR relaxation methods are examined. Nuclear Overhauser effects, spin-lattice, and spin-spin relaxation parameters of 15N NMR relaxation in ubiquitin have been determined at 17.6, 14.1, 11.7 and 9.4 Tesla. This unusually broad set of relaxation parameters has allowed the examination of the influence of chemical shift anisotropy, the functional form of the model-free spectral density, and the reliability of determined spin-spin relaxation parameters on the characterization of global tumbling of the protein. Treating the 15N chemical shift anisotropy (CSA) as an adjustable parameter, a consensus value of -170 +/- 15 ppm for the breadth of the chemical shift tensor and a global isotropic correlation time of 4.1 ns are found when using the model-free spectral density to fit T1 and NOE data from all fields. The inclusion of T2 relaxation parameters in the determination of the global correlation time results in its increase to 4.6 ns. This apparent inconsistency may explain a large portion of the discrepancy often found between NMR- and fluorescence-derived tau m values for proteins. The near identity of observed T2 and T1 rho values suggests that contributions from slow motions are not the origin of the apparent inconsistency with obtained T1 and NOE data. Various considerations suggest that the origin of this apparent discrepancy may reside in a contribution to the spectral density at zero frequency that is not represented by the simple model-free formalism in addition to the usual experimental difficulties associated with the measurement of these relaxation parameters. Finally, an axially symmetric diffusion tensor for ubiquitin is obtained using exclusively T1 and NOE data. A recommendation is reached on the types and combinations of relaxation data that can be used to reliably determine tau m values. It is also noted that the reliable determination of tau m values from 15N T1 and NOE relaxation parameters will become increasingly difficult as tau m increases.

High-resolution NMR of encapsulated proteins dissolved in low-viscosity fluids.

The majority of known proteins are too large to be comprehensively examined by solution NMR methods, primarily because they tumble too slowly in solution. Here we introduce an approach to making the NMR relaxation properties of large proteins amenable to modern solution NMR techniques. The encapsulation of a protein in a reverse micelle dissolved in a low-viscosity fluid allows it to tumble as fast as a much smaller protein. The approach is demonstrated and validated with the protein ubiquitin encapsulated in reverse micelles prepared in a variety of alkane solvents.

Local stability and dynamics of apocytochrome b562 examined by the dependence of hydrogen exchange on hydrostatic pressure.

Hydrostatic pressure is used to perturb the manifold of states available to apocytochrome b562 and to examine the energetics and dynamics of the protein using hydrogen exchange monitored in real-time by heteronuclear spectroscopy at pressures ranging up to 1. 1 kbar. An analytical framework for interpreting the effects of hydrostatic pressure on the physical events leading to protein hydrogen exchange is presented. The protein is found to have three regions of subglobal cooperative stability. The most stable region, or core, is composed of the central two helices of the bundle. The dependence of the global unfolding free energy upon pressure is first order and associated with a negative volume change of -102 mL mol-1. Two additional regions of cooperative structure are identified, and both also have negative volume changes associated with their unfolding at high pressure. Surprisingly, one of the subglobal unfolding units shows a significant positive volume change at low pressures (<200 bar) suggesting the presence of a highly mispacked open state at ambient pressure. The three regions of cooperative stability are the same as identified by perturbation with chemical denaturant. The implications of these results for issues in protein folding and the form of the energy landscape of globular proteins are discussed.