Induced fit DNA recognition by a minor groove binding analogue of Hoechst 33258: fluctuations in DNA A tract structure investigated by NMR and molecular dynamics simulations

NMR analysis and molecular dynamics simulations of d(GGTAATTACC)2 and its complex with a tetrahydropyrimidinium analogue of Hoechst 33258 suggest that DNA minor groove recognition in solution involves a combination of conformational selection and induced fit, rather than binding to a preorganised site. Analysis of structural fluctuations in the bound and unbound states suggests that the degree of induced fit observed is primarily a consequence of optimising van der Waals contacts with the walls of the minor groove resulting in groove narrowing through: (i) changes in base step parameters, including increased helical twist and propeller twist; (ii) changes to the sugar–phosphate backbone conformation to engulf the bound ligand; (iii) suppression of bending modes at the TpA steps. In contrast, the geometrical arrangement of hydrogen bond acceptors on the groove floor appears to be relatively insensitive to DNA conformation (helical twist and propeller twist). We suggest that effective recognition of DNA sequences (in this case an A tract structure) appears to depend to a significant extent on the sequence being flexible enough to be able to adopt the geometrically optimal conformation compatible with the various binding interactions, rather than involving ‘lock and key’ recognition.

NMR structure of the calreticulin P-domain

The NMR structure of the rat calreticulin P-domain, comprising residues 189–288, CRT(189–288), shows a hairpin fold that involves the entire polypeptide chain, has the two chain ends in close spatial proximity, and does not fold back on itself. This globally extended structure is stabilized by three antiparallel β-sheets, with the β-strands comprising the residues 189–192 and 276–279, 206–209 and 262–265, and 223–226 and 248–251, respectively. The hairpin loop of residues 227–247 and the two connecting regions between the β-sheets contain a hydrophobic cluster, where each of the three clusters includes two highly conserved tryptophyl residues, one from each strand of the hairpin. The three β-sheets and the three hydrophobic clusters form a repeating pattern of interactions across the hairpin that reflects the periodicity of the amino acid sequence, which consists of three 17-residue repeats followed by three 14-residue repeats. Within the global hairpin fold there are two well-ordered subdomains comprising the residues 219–258, and 189–209 and 262–284, respectively. These are separated by a poorly ordered linker region, so that the relative orientation of the two subdomains cannot be precisely described. The structure type observed for CRT(189–288) provides an additional basis for functional studies of the abundant endoplasmic reticulum chaperone calreticulin.

Dissection of the ion-induced folding of the hammerhead ribozyme using 19F NMR

We have used 19F NMR to analyze the metal ion-induced folding of the hammerhead ribozyme by selective incorporation of 5fluorouridine. We have studied the chemical shift and linewidths of 19F resonances of 5-fluorouridine at the 4 and 7 positions in the ribozyme core as a function of added Mg2+. The data fit well to a simple two-state model whereby the formation of domain 1 is induced by the noncooperative binding of Mg2+ with an association constant in the range of 100 to 500 M−1, depending on the concentration of monovalent ions present. The results are in excellent agreement with data reporting on changes in the global shape of the ribozyme. However, the NMR experiments exploit reporters located in the center of the RNA sections undergoing the folding transitions, thereby allowing the assignment of specific nucleotides to the separate stages. The results define the folding pathway at high resolution and provide a time scale for the first transition in the millisecond range.

New methods of structure refinement for macromolecular structure determination by NMR

Recent advances in multidimensional NMR methodology have permitted solution structures of proteins in excess of 250 residues to be solved. In this paper, we discuss several methods of structure refinement that promise to increase the accuracy of macromolecular structures determined by NMR. These methods include the use of a conformational database potential and direct refinement against three-bond coupling constants, secondary 13C shifts, 1H shifts, T1/T2 ratios, and residual dipolar couplings. The latter two measurements provide long range restraints that are not accessible by other solution NMR parameters.

Initiation sites of protein folding by NMR analysis.

Detailed characterization of denatured states of proteins is necessary to understand the interactions that funnel the large number of possible conformations along fast routes for folding. Nuclear magnetic resonance experiments based on the nuclear Overhauser effect (NOE) detect hydrogen atoms close in space and provide information about local structure. Here we present an NMR procedure that detects almost all sequential NOEs between amide hydrogen atoms (HN-HN NOE), including those in random coil regions in a protein, barnase, in urea solutions. A semi-quantitative analysis of these HN-HN NOEs identified partly structured regions that are in remarkable agreement with those found to form early on the reaction pathway. Our results strongly suggest that the folding of barnase initiates at the first helix and the beta-turn between the third and the fourth strands. This strategy of defining residual structure has also worked for cold-denatured barstar and guanidinium hydrochloride-denatured chymotrypsin inhibitor 2 and so should be generally applicable.

NMR evidence for the participation of a low-barrier hydrogen bond in the mechanism of delta 5-3-ketosteroid isomerase.

Delta 5-3-Ketosteroid isomerase (EC 5.3.3.1) promotes an allylic rearrangement involving intramolecular proton transfer via a dienolic intermediate. This enzyme enhances the catalytic rate by a factor of 10(10). Two residues, Tyr-14, the general acid that polarizes the steroid 3-carbonyl group and facilitates enolization, and Asp-38 the general base that abstracts and transfers the 4 beta-proton to the 6 beta-position, contribute 10(4.7) and 10(5.6) to the rate increase, respectively. A major mechanistic enigma is the huge disparity between the pKa values of the catalytic groups and their targets. Upon binding of an analog of the dienolate intermediate to isomerase, proton NMR detects a highly deshielded resonance at 18.15 ppm in proximity to aromatic protons, and with a 3-fold preference for protium over deuterium (fractionation factor, phi = 0.34), consistent with formation of a short, strong (low-barrier) hydrogen bond to Tyr-14. The strength of this hydrogen bond is estimated to be at least 7.1 kcal/mol. This bond is relatively inaccessible to bulk solvent and is pH insensitive. Low-barrier hydrogen bonding of Tyr-14 to the intermediate, in conjunction with the previously demonstrated tunneling contribution to the proton transfer by Asp-38, provide a plausible and quantitative explanation for the high catalytic power of this isomerase.

Increased tricarboxylic acid cycle flux in rat brain during forepaw stimulation detected with 1H[13C]NMR.

NMR spectroscopy was used to test recent proposals that the additional energy required for brain activation is provided through nonoxidative glycolysis. Using localized NMR spectroscopic methods, the rate of C4-glutamate isotopic turnover from infused [1-(13)C]glucose was measured in the somatosensory cortex of rat brain both at rest and during forepaw stimulation. Analysis of the glutamate turnover data using a mathematical model of cerebral glucose metabolism showed that the tricarboxylic acid cycle flux [(V(TCA)] increased from 0.49 +/- 0.03 at rest to 1.48 +/- 0.82 micromol/g/min during stimulation (P < 0.01). The minimum fraction of C4-glutamate derived from C1-glucose was approximately 75%, and this fraction was found in both the resting and stimulated rats. Hence, the percentage increase in oxidative cerebral metabolic rate of glucose use (CMRglc) equals the percentage increases in V(TCA) and cerebral metabolic rate of oxygen consumption (CMRO2). Comparison with previous work for the same rat model, which measured total CMRglc [Ueki, M., Linn, F. & Hossman, K. A. (1988) J. Cereb. Blood Flow Metab. 8, 486-4941, indicates that oxidative CMRglc supplies the majority of energy during sustained brain activation.

Natural abundance solid-state carbon NMR studies of photosynthetic reaction centers with photoinduced polarization.

Solid-state NMR spectra of natural abundance 13C in reaction centers from photosynthetic bacteria Rhodobacter sphaeroides R-26 was measured. When the quinone acceptors were removed and continuous visible illumination of the sample was provided, exceptionally strong nuclear spin polarization was observed in NMR lines with chemical shifts resembling those of the aromatic carbons in bacteriochlorophyll and bacteriopheophytin. The observation of spin polarized 15N nuclei in bacteriochlorophyll and bacteriopheophytin was previously demonstrated with nonspecifically 15N-labeled reaction centers. Both the carbon and the nitrogen NMR studies indicate that the polarization is developed on species that carry unpaired electrons in the early electron transfer steps, including the bacteriochlorophyll dimer donor P860 and probably the bacteriopheophytin acceptor. I. Both enhanced-absorptive and emissive polarization were seen in the carbon spectrum; most lines were absorptive but the methine carbons of the porphyrin ring (alpha, beta, gamma, ) exhibited emissive polarization. The change in the sign of the hyperfine coupling at these sites indicates the existence of nodes in the spin density distribution on the tetrapyrrole cofactors flanking each methine carbon bridge.

31P NMR magnetization transfer study of the control of ATP turnover in Saccharomyces cerevisiae.

31P NMR magnetization transfer measurements have been used to measure the steady state flux between Pi and ATP in yeast cells genetically modified to overexpress an adenine nucleotide translocase isoform. An increase in Pi -> ATP flux and apparent ratio of moles of ATP synthesized/atoms of oxygen consumed (P:O ratio), when these cells were incubated with glucose, demonstrated that the reactions catalyzed by the translocase and F1F0 ATP synthase were readily reversible in vivo. However, when the same cells were incubated with ethanol alone, translocase overexpression had no effect on the measured Pi -> ATP flux or apparent P:O ratio, suggesting that the synthase was now operating irreversibly. This change was accompanied by an increase in the intracellular ADP concentration. These observations are consistent with a model proposed for the kinetic control of mitochondrial ATP synthesis, which was based on isotope exchange measurements with isolated mammalian mitochondria [LaNoue, K. F., Jeffries, F. M. H. & Radda, G. K. (1986) Biochemistry 25, 7667-7675].

NMR studies of muscle glycogen synthesis in insulin-resistant offspring of parents with non-insulin-dependent diabetes mellitus immediately after glycogen-depleting exercise.

To examine the impact of insulin resistance on the insulin-dependent and insulin-independent portions of muscle glycogen synthesis during recovery from exercise, we studied eight young, lean, normoglycemic insulin-resistant (IR) offspring of individuals with non-insulin-dependent diabetes mellitus and eight age-weight matched control (CON) subjects after plantar flexion exercise that lowered muscle glycogen to approximately 25% of resting concentration. After approximately 20 min of exercise, intramuscular glucose 6-phosphate and glycogen were simultaneously monitored with 31P and 13C NMR spectroscopies. The postexercise rate of glycogen resynthesis was nonlinear. Glycogen synthesis rates during the initial insulin independent portion (0-1 hr of recovery) were similar in the two groups (IR, 15.5 +/- 1.3 mM/hr and CON, 15.8 +/- 1.7 mM/hr); however, over the next 4 hr, insulin-dependent glycogen synthesis was significantly reduced in the IR group [IR, 0.1 +/- 0.5 mM/hr and CON, 2.9 +/- 0.2 mM/hr; (P < or = 0.001)]. After exercise there was an initial rise in glucose 6-phosphate concentrations that returned to baseline after the first hour of recovery in both groups. In summary, we found that following muscle glycogen-depleting exercise, IR offspring of parents with non-insulin-dependent diabetes mellitus had (i) normal rates of muscle glycogen synthesis during the insulin-independent phase of recovery from exercise and (ii) severely diminished rates of muscle glycogen synthesis during the subsequent recovery period (2-5 hr), which has previously been shown to be insulin-dependent in normal CON subjects. These data provide evidence that exercise and insulin stimulate muscle glycogen synthesis in humans by different mechanisms and that in the IR subjects the early response to stimulation by exercise is normal.

The catalytic mechanism of beta-lactamases: NMR titration of an active-site lysine residue of the TEM-1 enzyme.

Beta-Lactamases are widespread in the bacterial world, where they are responsible for resistance to penicillins, cephalosporins, and related compounds, currently the most widely used antibacterial agents. Detailed structural and mechanistic understanding of these enzymes can be expected to guide the design of new antibacterial compounds resistant to their action. A number of high-resolution structures are available for class A beta-lactamases, whose catalytic mechanism involves the acylation of a serine residue at the active site. The identity of the general base which participates in the activation of this serine residue during catalysis has been the subject of controversy, both a lysine residue and a glutamic acid residue having been proposed as candidates for this role. We have used the pH dependence of chemical modification of epsilon-amino groups by 2,4,6,-trinitrobenzenesulfonate and the pH dependence of the epsilon-methylene 1H and 13C chemical shifts (in enzyme selectively labeled with [epsilon-13C]lysine) to estimate the pKa of the relevant lysine residue, lysine-73, of TEM-1 beta-lactamase. Both methods show that the pKa of this residue is > 10, making it very unlikely that this residue could act as a proton acceptor in catalysis. An alternative mechanism in which this role is performed by glutamate-166 through an intervening water molecule is described.

Observation via one-dimensional 13Calpha NMR of local conformational substates in thermal unfolding equilibria of a synthetic analog of the GCN4 leucine zipper.

Synthesis of a 33-residue, capped leucine zipper analogous to that in GCN4 is reported. Histidine and arginine residues are mutated to lysine to reduce the unfolding temperature. CD and ultracentrifugation studies indicate that the molecule is a two-stranded coiled coil under benign conditions. Versions of the same peptide are made with 99% 13Calpha at selected sites. One-dimensional 13C NMR spectra are assigned by inspection and used to study thermal unfolding equilibria over the entire transition from 8 to 73 degrees C. Spectra at the temperature extremes establish the approximate chemical shifts for folded and unfolded forms at each labeled site. Resonances for each amino acid appear at both locations at intermediate T, indicating that folded and unfolded forms interconvert slowly (> >2 ms) on the NMR time scale. Moreover, near room temperature, the structured form's resonance is double at several, but not all, sites, indicating at least two slowly interconverting, structured, local conformational substates. Analysis of the dynamics is possible. For example, near room temperature at the Val-9, Ala-24, and Gly-31 positions, the equilibrium constant for interconversion of the two structured forms is near unity and the time scale is > or= 10-20 ms.

Determining RNA solution structure by segmental isotopic labeling and NMR: application to Caenorhabditis elegans spliced leader RNA 1.

Recent developments in multidimensional heteronuclear NMR spectroscopy and large-scale synthesis of uniformly 13C- and 15N-labeled oligonucleotides have greatly improved the prospects for determination of the solution structure of RNA. However, there are circumstances in which it may be advantageous to label only a segment of the entire RNA chain. For example, in a larger RNA molecule the structural question of interest may reside in a localized domain. Labeling only the corresponding nucleotides simplifies the spectrum and resonance assignments because one can filter proton spectra for coupling to 13C and 15N. Another example is in resolving alternative secondary structure models that are indistinguishable in imino proton connectivities. Here we report a general method for enzymatic synthesis of quantities of segmentally labeled RNA molecules required for NMR spectroscopy. We use the method to distinguish definitively two competing secondary structure models for the 5' half of Caenorhabditis elegans spliced leader RNA by comparison of the two-dimensional [15N] 1H heteronuclear multiple quantum correlation spectrum of the uniformly labeled sample with that of a segmentally labeled sample. The method requires relatively small samples; solutions in the 200-300 microM concentration range, with a total of 30 nmol or approximately 40 micrograms of RNA in approximately 150 microliters, give strong NMR signals in a short accumulation time. The method can be adapted to label an internal segment of a larger RNA chain for study of localized structural problems. This definitive approach provides an alternative to the more common enzymatic and chemical footprinting methods for determination of RNA secondary structure.

Stopped-flow NMR spectroscopy: real-time unfolding studies of 6-19F-tryptophan-labeled Escherichia coli dihydrofolate reductase.

Escherichia coli dihydrofolate reductase (DHFR; EC 1.5.1.3) contains five tryptophan residues that have been replaced with 6-19F-tryptophan. The 19F NMR assignments are known in the native, unliganded form and the unfolded form. We have used these assignments with stopped-flow 19F NMR spectroscopy to investigate the behavior of specific regions of the protein in real time during urea-induced unfolding. The NMR data show that within 1.5 sec most of the intensities of the native 19F resonances of the protein are lost but only a fraction (approximately 20%) of the intensities of the unfolded resonances appears. We postulate that the early disappearance of the native resonances indicates that most of the protein rapidly forms an intermediate in which the side chains have considerable mobility. Stopped-flow far-UV circular dichroism measurements indicate that this intermediate retains native-like secondary structure. Eighty percent of the intensities of the NMR resonances assigned to the individual tryptophans in the unfolded state appear with similar rate constants (k approximately 0.14 sec-1), consistent with the major phase of unfolding observed by stopped-flow circular dichroism (representing 80% of total amplitude). These data imply that after formation of the intermediate, which appears to represent an expanded structural form, all regions of the protein unfold at the same rate. Stopped-flow measurements of the fluorescence and circular dichroism changes associated with the urea-induced unfolding show a fast phase (half-time of about 1 sec) representing 20% of the total amplitude in addition to the slow phase mentioned above. The NMR data show that approximately 20% of the total intensity for each of the unfolded tryptophan resonances is present at 1.5 sec, indicating that these two phases may represent the complete unfolding of the two different populations of the native protein.