Lactobacillus demonstrate thiol-independent metabolism of methylglyoxal: Implications toward browning prevention in Parmesan cheese.

Endogenous production of α-dicarbonyls by lactic acid bacteria can influence the quality and consistency of fermented foods and beverages. Methylglyoxal (MG) in Parmesan cheese can contribute toward undesired browning during low temperature ripening and storage conditions, leading to the economic depreciation of affected cheeses. We demonstrate the effects of exogenously added MG on browning and volatile formation using a Parmesan cheese extract (PCE). To determine the influence of Lactobacillus on α-dicarbonyls, strains were screened for their ability to modulate concentrations of MG, glyoxal, and diacetyl in PCE. It was found that a major metabolic pathway of MG in Lactobacillus is a thiol-independent reduction, whereby MG is partially or fully reduced to acetol and 1,2-propanediol, respectively. The majority of lactobacilli grown in PCE accumulated the intermediate acetol, whereas Lactobacillus brevis 367 formed exclusively 1,2-propanediol and Lactobacillus fermentum 14931 formed both metabolites. In addition, we determined the inherent tolerance to bacteriostatic concentrations of MG among lactobacilli grown in rich media. It was found that L. brevis 367 reduces MG exclusively to 1,2-propanediol, which correlates to both its ability to significantly decrease MG concentrations in PCE, as well as its significantly higher tolerance to MG, in comparison to other lactobacilli screened. These findings have broader implications toward lactobacilli as a viable solution for reducing MG-mediated browning of Parmesan cheese.

Natural J-coupling analysis: interpretation of scalar J-couplings in terms of natural bond orbitals.

The natural J-coupling (NJC) method presented here analyzes the Fermi contact portion of J-coupling in the framework of finite perturbation theory applied to ab initio/density function theory (DFT) wave functions, to compute individual and pairwise orbital contributions to the net J-coupling. The approach is based on the concepts and formalisms of natural bond orbital (NBO) methods. Computed coupling contributions can be classified as Lewis (individual orbital contributions corresponding to the natural Lewis structure of the molecule), delocalization (resulting from pairwise donor-acceptor interactions), and residual repolarization (corresponding to correlation-like interactions). This approach is illustrated by an analysis of the angular and distance dependences of the contributions to vicinal (3)J(HH) couplings in ethane and to the long-range (6)J(HH) couplings in pentane. The results indicate that approximately 70% or more of the net J-coupling is propagated by steric exchange antisymmetry interactions between Lewis orbitals (predominantly sigma bonding orbitals). Hyperconjugative sigma to sigma delocalization interactions account for the remainder of the coupling. Calculated pairwise-steric and hyperconjugative-delocalization energies provide a means for relating coupling mechanisms to molecular energetics. In this way, J-coupling contributions can be related directly to the localized features of the molecular electronic structure in order to explain measured J-coupling patterns and to predict J-coupling trends that have yet to be measured.

Monomeric solution structure of the prototypical 'C' chemokine lymphotactin.

Lymphotactin, the sole identified member of the C class of chemokines, specifically attracts T lymphocytes and natural killer cells. This 93-residue protein lacks 2 of the 4 conserved cysteine residues characteristic of the other 3 classes of chemokines and possesses an extended carboxyl terminus, which is required for chemotactic activity. We have determined the three-dimensional solution structure of recombinant human lymphotactin by NMR spectroscopy. Under the conditions used for the structure determination, lymphotactin was predominantly monomeric; however, pulsed field gradient NMR self-diffusion measurements and analytical ultracentrifugation revealed evidence of dimer formation. Sequence-specific chemical shift assignments were determined through analysis of two- and three-dimensional NMR spectra of (15)N- and (13)C/(15)N-enriched protein samples. Input for the torsion angle dynamics calculations used in determining the structure included 1258 unique NOE-derived distance constraints and 60 dihedral angle constraints obtained from chemical-shift-based searching of a protein conformational database. The ensemble of 20 structures chosen to represent the structure had backbone and heavy atom rms deviations of 0.46 +/- 0.11 and 1.02 +/- 0.14 A, respectively. The results revealed that human lymphotactin adopts the conserved chemokine fold, which is characterized by a three-stranded antiparallel beta-sheet and a C-terminal alpha-helix. Two regions are dynamically disordered as evidenced by (1)H and (13)C chemical shifts and [(15)N]-(1)H NOEs: residues 1-9 of the amino terminus and residues 69-93 of the C-terminal extension. A functional role for the C-terminal extension, which is unique to lymphotactin, remains to be elucidated.

NMR chemical shift mapping of the binding site of a protein proteinase inhibitor: changes in the (1)H, (13)C and (15)N NMR chemical shifts of turkey ovomucoid third domain upon binding to bovine chymotrypsin A(alpha).

The substrate-like inhibition of serine proteinases by avian ovomucoid domains has provided an excellent model for protein inhibitor-proteinase interactions of the standard type. 1H,15N and 13C NMR studies have been undertaken on complexes formed between turkey ovomucoid third domain (OMTKY3)2 and chymotrypsin A(alpha) (Ctr) in order to characterize structural changes occurring in the Ctr binding site of OMTKY3. 15N and 13C were incorporated uniformly into OMTKY3, allowing backbone resonances to be assigned for OMTKY3 in both its free and complex states. Chemical shift perturbation mapping indicates that the two regions, K13-P22 and N33-A40, are the primary sites in OMTKY3 involved in Ctr binding, in full agreement with the 12 consensus proteinase-contact residues of OMTKY3 defined previously on the basis of X-ray crystallographic and mutational analysis. Smaller chemical shift perturbations in selected other regions may result from minor structural changes on binding. Through-bond 15N-13C correlations between P1-13C' and P1'-15N in two-dimensional H(N)CO and HN(CO) NMR spectra of selectively labeled OMTKY3 complexed with Ctr indicate that the scissile peptide bond between L18 and E19 of the inhibitor is intact in the complex. The chemical shifts of the reactive site peptide bond indicate that it is predominantly trigonal, although the data are not inconsistent with a slight perturbation of the hybridization of the peptide bond toward the first tetrahedral state along the reaction coordinate.

Ovine caveolin-1: cDNA cloning, E. coli expression, and association with endothelial nitric oxide synthase.

Caveolin-1 (Cav-1), the principal coat protein of caveolae, plays an obligatory role in regulating the activity of endothelial nitric oxide (NO) synthase (eNOS). We propose that Cav-1 may be critical to eNOS-NO mediated uterine vasodilatation during pregnancy and estrogen replacement therapy. To test this hypothesis in the sheep model, we isolated the full-length cDNA of ovine Cav-1 (oCav-1) from a Lambda ZAP cDNA library of ovine placental artery endothelial cells. Thirty-two positive oCav-1 clones were recognized by a partial oCav-1 cDNA from this library, of which eight were sequenced. Restriction digestion of these clones revealed that the cDNAs of oCav-1 ranged from approximately 2.1 to 2.7 kb. Northern analysis of Cav-1 mRNAs in ovine uterine artery endothelial cells (UAEC) showed two transcripts of approximately 2.1 and 2.7 kb, respectively. Immunoreactive Cav-1 protein, but not caveolin-2 or caveolin-3, was detected in UAEC. Sequence analysis revealed that in addition to a 537-bp open reading frame encoding a 178 amino acid oCav-1 protein, full-length oCav-1 cDNAs apparently possess a approximately 1.6-2.1 kb 3'-untranslated region. Database searches with oCav-1 cDNA revealed that the coding region of mammalian Cav-1 genes is highly conserved. We prepared a recombinant full-length oCav-1 protein in which six consecutive histidine residues were tagged at the end of its COOH-terminus and developed a [His]6-tagged oCav-1 'pull-down assay' for studying the association of eNOS with Cav-1. Incubation of exogenous [His]6-tagged oCav-1 with resting UAEC extracts led to the formation of a [His]6-tagged oCav-1-eNOS complex. In the presence of a synthetic caveolin-scaffolding domain (CSD, aa 82-101) peptide, but not a mutated CSD peptide, [His]6-tagged oCav-1 associated eNOS was dose (0-10 microM)-dependently inhibited. eNOS association with Cav-1 in UAEC was further confirmed by the facts that eNOS co-immunoprecipitated with Cav-1 and vice versa, and that eNOS co-existed with Cav-1 during the isolation of caveolae membranes. Because dissociation of eNOS from Cav-1 is required for the activation of eNOS, eNOS association with Cav-1 in UAEC suggests an important role of Cav-1 in regulating UAEC production of NO and possibly NO-mediated uterine vasodilatation.

Solution structure of the toluene 4-monooxygenase effector protein (T4moD).

Toluene 4-monooxygenase (T4MO) from Pseudomonas mendocina catalyzes the NADH- and O(2)-dependent hydroxylation of toluene to form p-cresol. The complex consists of an NADH oxidoreductase (T4moF), a Rieske ferredoxin (T4moC), a diiron hydroxylase [T4moH, with (alphabetagamma)(2) quaternary structure], and a catalytic effector protein (T4moD). The solution structure of the 102-amino acid T4moD effector protein has been determined from 2D and 3D (1)H, (13)C, and (15)N NMR spectroscopic data. The structural model was refined through simulated annealing by molecular dynamics in torsion angle space (DYANA software) with input from 1467 experimental constraints, comprising 1259 distance constraints obtained from NOEs, 128 dihedral angle constraints from J-couplings, and 80 hydrogen bond constraints. Of 60 conformers that met the acceptance criteria, the 20 that best satisfied the input constraints were selected to represent the solution structure. With exclusion of the ill-defined N- and C-terminal segments (Ser1-Asn11 and Asp99-Met102), the atomic root-mean-square deviation for the 20 conformers with respect to the mean coordinates was 0.71 A for the backbone and 1.24 A for all non-hydrogen atoms. The secondary structure of T4moD consists of three alpha-helices and seven beta-strands arranged in an N-terminal betaalphabetabeta and a C-terminal betaalphaalphabetabetabeta domain topology. Although the published NMR structures of the methane monooxygenase effector proteins from Methylosinus trichosporium OB3b and Methylococcus capsulatus (Bath) have a similar secondary structure topology, their three-dimensional structures differ from that of T4moD. The major differences in the structures of the three effector proteins are in the relative orientations of the two beta-sheets and the interactions between the alpha-helices in the two domains. The structure of T4moD is closer to that of the methane monooxygenase effector protein from M. capsulatus (Bath) than that from M. trichosporium OB3b. The specificity of T4moD as an effector protein was investigated by replacing it in reconstituted T4MO complexes with effector proteins from monooxygenases from other bacterial species: Pseudomonas pickettii PKO1 (TbuV, toluene 3-monooxygenase); Pseudomonas species JS150 (TbmC, toluene 2-monooxygenase); and Burkeholderia cepacia G4 (S1, toluene 2-monooxygenase). The results showed that the closely related TbuV effector protein (55% sequence identity) provided partial activation of the complex, whereas the more distantly related TbmC (34% sequence identity) and S1 (29% sequence identity) did not. The (1)H NMR chemical shifts of the side-chain amide protons of Asn34, a conserved, structurally relevant amino acid, were found to be similar in spectra of effector proteins T4moD and TbuV but not in the spectrum of TbmC. This suggests that the region around Asn34 may be involved in structural aspects contributing to functional specificity.

Functional recombinant rabbit muscle phosphoglucomutase from Escherichia coli.

The gene coding for phosphoglucomutase (PGM) from Oryctolagus cuniculus (rabbit) has been expressed in Escherichia coli under a T7 expression system with a His-tag. About half of the expressed PGM protein was present in inclusion bodies, but this protein was inactive when solubilized. The protein in the soluble cell fraction was isolated and purified in one step on a Ni-NTA column. The eluate from this column was adjusted to 95% saturated ammonium sulfate, and the resulting protein precipitate was resuspended in sodium phosphate buffer and dialyzed against 2.5 M ammonium sulfate. The final yield of protein was about 10 mg per liter of LB medium. The protein was judged to be greater than 90% pure on the basis of gel electrophoresis and activity measurements (128 U per milligram). Our motivation for developing this bacterial production system for PGM has been to prepare sufficient quantities of stable-isotope-labeled protein for experiments that utilize recently developed NMR technologies suitable for proteins the size of PGM (61.6 kDa). Preliminary NMR studies indicate that the current level of purity is adequate for this work. The construct described here was designed to incorporate an N-terminal His-tag for ease of isolation. Although PGM is a metalloprotein, the His-tag does not appear to interfere with activity. The presence of the His-tag should not pose a problem for proposed (31)P NMR investigations of the protein and its complexes in aqueous solution or incorporated into reverse micelles. However, we plan to design a cleavable His-tag for later (1)H, (13)C, (15)N studies of the active site, which includes essential histidine residues.

Efficient production of recombinant brazzein, a small, heat-stable, sweet-tasting protein of plant origin.

Brazzein is a 54-amino-acid sweet-tasting protein first isolated from the fruit of Pentadiplandra brazzeana Baillon found in West Africa. Brazzein, as isolated from the fruit, is 500 times sweeter than sucrose on a weight basis (9500 times sweeter on a per-molecule basis). A minor component of brazzein from fruit, des-pGlu1-brazzein, has 53 amino acid residues and has twice the sweetness of the parent protein. We have designed a gene for des-pGlu1- brazzein that incorporates codons that are optimal for protein production in Escherichia coli. Production of brazzein from the chemically synthesized gene resulted in recombinant protein with sweetness similar to that of brazzein isolated from the original source. The best yields were achieved by producing brazzein as a fusion with staphylococcal nuclease with a designed cyanogen bromide cleavage site. Because of its intense sweetness and stability at high pH and temperature, brazzein is an ideal system for investigating the chemical and structural requirements involved in sweet-taste properties. This efficient protein production system for brazzein will facilitate such investigations.

Sweetness determinant sites of brazzein, a small, heat-stable, sweet-tasting protein.

Brazzein, originally isolated from the fruit of the African plant Pentadiplandra brazzeana Baillon, is the smallest, most heat-stable and pH-stable member of the set of proteins known to have intrinsic sweetness. These properties make brazzein an ideal system for investigating the chemical and structural requirements of a sweet-tasting protein. We have used the three-dimensional structure of the protein (J. E. Caldwell et al. (1998) Nat. Struct. Biol. 5, 427-431) as a guide in designing 15 synthetic genes in expression constructs aimed at delineating the sweetness determinants of brazzein. Protein was produced heterologously in Escherichia coli, isolated, and purified as described in the companion paper (Assadi-Porter, F. M., Aceti, D., Cheng, H., and Markley, J. L., this issue). Analysis by one-dimensional (1)H NMR spectroscopy indicated that all but one of these variants had folded properly under the conditions used. A taste panel compared the gustatory properties of solutions of these proteins to those of sucrose and brazzein isolated from fruit. Of the 14 mutations in the des-pGlu1-brazzein background, four exhibited almost no sweetness, six had significantly reduced sweetness, two had taste properties equivalent to des-pGlu1-brazzein (two times as sweet as the major form of brazzein isolated from fruit which contains pGlu1), and two were about twice as sweet as des-pGlu1-brazzein. Overall, the results suggest that two regions of the protein are critical for the sweetness of brazzein: a region that includes the N- and C-termini of the protein, which are located close to one another, and a region that includes the flexible loop around Arg43.

Electron-nuclear interactions in two prototypical [2Fe-2S] proteins: selective (chiral) deuteration and analysis of (1)H and (2)H NMR signals from the alpha and beta hydrogens of cysteinyl residues that ligate the iron in the active sites of human ferredoxin and Anabaena 7120 vegetative ferredoxin.

A vertebrate ferredoxin (human ferredoxin) and a plant-type ferredoxin (the ferredoxin from the vegetative form of Anabaena 7120) were labeled selectively with deuterium at their active site cysteines. The recombinant proteins were produced in Escherichia coli and labeled by replacing natural abundance cysteine in the defined culture medium with (2)H(alpha)-cysteine, (2)H(beta2), (2)H(beta3)-cysteine, or (2)H(beta2)-cystine. The chiral labeled cystine ((2)H(beta2)-cystine) was prepared by selective hydrogen exchange catalyzed by cystathionine gamma-synthase. NMR spectra of these samples in their oxidized and reduced states support unambiguous identifications by atom type of (1)H and (2)H NMR signals from the cysteine alpha and beta hydrogens. These signals lie outside the normal diamagnetic spectral region as a result of interaction of the hydrogens with unpaired electron density from the iron-sulfur cluster, and their chemical shifts are highly dependent on local conformation at the active site. The very different chemical properties of the iron centers of plant-type and vertebrate ferredoxins reflect relatively small differences in the conformation of the iron-sulfur cluster ligands.

Oligomerization of the EK18 mutant of the trp repressor of Escherichia coli as observed by NMR spectroscopy.

The regulation of the trp repressor system of Escherichia coli is frequently modeled by a single equilibrium, that between the aporepressor (TR) and the corepressor, l-tryptophan (Trp), at their intracellular concentrations. The actual mechanism, which is much more complex and more finely tuned, involves multiple equilibria: TR and Trp association, TR oligomerization, specific and nonspecific binding of various states of TR to DNA, and interactions between these various species and ions. TR in isolation exists primarily as a homodimer, but the state of oligomerization increases as the TR concentration goes up and/or the salt concentration goes down, leading to species with lower affinity for DNA. We have used multinuclear, multidimensional NMR spectroscopy to investigate structural changes that accompany the oligomerization of TR. For these investigations, the superrepressor mutant EK18 (TR with Glu 18 replaced by Lys) was chosen because it exhibits less severe oligomerization at higher protein concentration than other known variants; this made it possible to study the dimer to tetramer oligomerization step by NMR. The NMR results suggest that the interaction between TR dimers is structurally linked to folding of the DNA binding domain and that it likely involves direct contacts between the C-terminal residues of the C-helix of one dimer with the next dimer. This implies that oligomerization can compete with DNA binding and thus serves as a factor in the fine-tuning of gene expression.

Thermodynamic criterion for the conformation of P1 residues of substrates and of inhibitors in complexes with serine proteinases.

Eglin c, turkey ovomucoid third domain, and bovine pancreatic trypsin inhibitor (Kunitz) are all standard mechanism, canonical protein inhibitors of serine proteinases. Each of the three belongs to a different inhibitor family. Therefore, all three have the same canonical conformation in their combining loops but differ in their scaffoldings. Eglin c (Leu45 at P1) binds to chymotrypsin much better than its Ala45 variant (the difference in standard free energy changes on binding is -5.00 kcal/mol). Similarly, turkey ovomucoid third domain (Leu18 at P1) binds to chymotrypsin much better than its Ala18 variant (the difference in standard free energy changes on binding is -4.70 kcal/mol). As these two differences are within the +/-400 cal/mol bandwidth (expected from the experimental error), one can conclude that the system is additive. On the basis that isoenergetic is isostructural, we expect that within both the P1 Ala pair and the P1 Leu pair, the conformation of the inhibitor's P1 side chain and of the enzyme's specificity pocket will be identical. This is confirmed, within the experimental error, by the available X-ray structures of complexes of bovine chymotrypsin Aalpha with eglin c () and with turkey ovomucoid third domain (). A comparison can also be made between the structures of P1 (Lys+)15 of bovine pancreatic trypsin inhibitor (Kunitz) ( and ) and of the P1 (Lys+)18 variant of turkey ovomucoid third domain (), both interacting with chymotrypsin. In this case, the conformation of the side chains is strikingly different. Bovine pancreatic trypsin inhibitor with (Lys+)15 at P1 binds to chymotrypsin more strongly than its Ala15 variant (the difference in standard free energy changes on binding is -1.90 kcal/mol). In contrast, turkey ovomucoid third domain variant with (Lys+)18 at P1 binds to chymotrypsin less strongly than its Ala18 variant (the difference in standard free energies of association is 0.95 kcal/mol). In this case, P1 Lys+ is neither isostructural nor isoenergetic. Thus, a thermodynamic criterion for whether the conformation of a P1 side chain in the complex matches that of an already determined one is at hand. Such a criterion may be useful in reducing the number of required X-ray crystallographic structure determinations. More importantly, the criterion can be applied to situations where direct determination of the structure is extremely difficult. Here, we apply it to determine the conformation of the Lys+ side chain in the transition state complex of a substrate with chymotrypsin. On the basis of kcat/KM measurements, the difference in free energies of activation for Suc-AAPX-pna when X is Lys+ and X is Ala is 1.29 kcal/mol. This is in good agreement with the corresponding difference for turkey ovomucoid third domain variants but in sharp contrast to the bovine pancreatic trypsin inhibitor (Kunitz) data. Therefore, we expect that in the transition state complex of this substrate with chymotrypsin, the P1 Lys+ side chain is deeply inserted into the enzyme's specificity pocket as it is in the (Lys+)18 turkey ovomucoid third domain complex with chymotrypsin.

NMR studies of retinoid-protein interactions: the conformation of [13C]-beta-ionones bound to beta-lactoglobulin B.

PURPOSE: Vitamin A (retinol) and its metabolites comprise the natural retinoids. While the biological action of these molecules are thought to be primarily mediated by ca. 55 kDa nuclear retinoic acid receptors, a number of structurally similar 15-20 kDa proteins are involved in the transport, and possibly metabolism, of these compounds. The milk protein beta-lactoglobulin B (beta-LG) is an 18 kDa protein which binds retinol and may be involved in oral delivery of retinol to neonates. beta-LG also binds drugs and other natural products and is of potential interest as a protective delivery vehicle. METHODS: To examine the conformation of the model retinoid beta-ionone both in solution and when bound to beta-LG, NMR and computational methods have been employed. RESULTS: Taken together, NMR studies of beta-ionone in solution measuring scalar and dipolar coupling, as well as CHARMm calculations, suggest beta-ionone prefers a slightly twisted 6-s-cis conformation. Isotope-edited NMR studies of 13C-labeled beta-ionones bound to beta-LG, primarily employing the HMQC-NOE experiment, suggest beta-ionone also binds to beta-LG in its 6-s-cis conformation. CONCLUSIONS: The methods employed here allow estimates of protein-bound ligand conformation. However, additional sites of ligand labeling will be necessary to aid in binding site localization.

Structure of the insect cytokine peptide plasmatocyte-spreading peptide 1 from Pseudoplusia includens.

The structure of the recently identified plasmatocyte spreading peptide from the moth Pseudoplusia includens (PSP1) has been determined by NMR spectroscopy. This novel insect cytokine consists of 23 amino acid residues and a single disulfide bond. Torsion angle dynamics calculations utilizing a total of 337 distance constraints yielded an ensemble of 30 structures with an average backbone root mean square deviation for residues 7-22 of 0.18 A from the mean structure. The structure consists of a disordered N-terminal region and a well defined core that is stabilized by numerous hydrophobic interactions and a short beta-hairpin. Structural comparisons confirm that PSP1 adopts an epidermal growth factor (EGF)-like fold with close similarity to the C-terminal subdomain of EGF-like module 5 of human thrombomodulin. The combination of the three-dimensional structure of PSP1 and the extensive literature on EGF-receptor interactions should accelerate the process of identifying the specific residues responsible for receptor binding activity of this family of immunoregulatory peptides.

Properties of the His57-Asp102 dyad of rat trypsin D189S in the zymogen, activated enzyme, and alpha1-proteinase inhibitor complexed forms.

Structural and biochemical studies suggest that serpins induce structural rearrangements in their target serine-proteinases. Previous NMR studies of the complex between a serpin, alpha1-proteinase inhibitor, and a mutant of recombinant rat trypsin (the Asp189 to Ser mutant, D189S, which is much more stable than wild-type rat trypsin against autoproteolysis) provided information about the state of catalytic residues in this complex: the hydrogen bond between Asp102 and His57 remains intact in the complex, and spectral properties of His57 are more like those of the zymogen than of the activated enzyme (G. Kaslik, et al., 1997, Biochemistry 36, 5455-5464). Here we report the protonation and exchange behavior of His57 of recombinant rat trypsin D189S in three states: the zymogen, the active enzyme, and the complex with human alpha1-proteinase inhibitor and compare these with analogous behavior of His57 of bovine chymotrypsinogen and alpha-chymotrypsin. In these studies the pKa of His57 has been determined from the pH dependence of the 1H NMR signal from the Hdelta1 proton of histidine in the Asp102-His57 dyad, and a measure of the accessibility of this part of the active site has been obtained from the rate of appearance of this signal following its selective saturation. The activation of rat trypsinogen D189S (zymogen, pKa = 7.8 +/- 0.1; Hill coefficient = 0. 86 +/- 0.05) decreased the pKa of His57 by 1.1 unit and made the protonation process cooperative (active enzyme, pKa = 6.7 +/- 0.1; Hill coefficient = 1.37 +/- 0.08). The binding of alpha1-proteinase inhibitor to trypsin D189S led to an increase in the pKa value of His57 to a value higher than that of the zymogen and led to negative cooperativity in the protonation process (complex, pKa = 8.1 +/- 0. 1; Hill coefficient = 0.70 +/- 0.08), as was observed for the zymogen. In spite of these differences in the pKa of His57 in the zymogen, active enzyme, and alpha1-proteinase inhibitor complex, the solvent exchange lifetime of the His57 Hdelta1 proton was the same, within experimental error, in all three states (lifetime = 2 to 12.5 ms). The linewidth of the 1H NMR signal from the Hdelta1 proton of His57 was relatively sharp, at temperatures between 5 and 20 degrees C at both low pH (5.2) and high pH (10.0), in spectra of bovine alpha-chymotrypsin, recombinant rat trypsin D189S, and the complex between rat trypsin D189S and human alpha1-proteinase inhibitor; however, in spectra of the complex between alpha-chymotrypsin and human alpha1-proteinase inhibitor, the peak was broader and could be well-resolved only at the lower temperature (5 degrees C).

Detection and classification of hyperfine-shifted 1H, 2H, and 15N resonances of the Rieske ferredoxin component of toluene 4-monooxygenase.

T4MOC is a 12.3 kDa soluble Rieske ferredoxin that is obligately required for electron transfer between the oxidoreductase and diiron hydroxylase components of toluene 4-monooxygenase from Pseudomonas mendocina KR1. Our preliminary 1H NMR studies of oxidized and reduced T4MOC [Markley, J. L., Xia, B., Chae, Y. K., Cheng, H., Westler, W. M., Pikus, J. D., and Fox, B. G. (1996) in Protein Structure Function Relationships (Zaidi, Z., and Smith, D., Eds.) pp 135-146, Plenum Press, London] revealed the presence of hyperfine-shifted 1H resonances whose short relaxation times made it impractical to use nuclear Overhauser effect (NOE) measurements for assignment purposes. We report here the use of selective isotopic labeling to analyze the hyperfine-shifted 1H, 2H, and 15N signals from T4MOC. Selective deuteration led to identification of signals from the four Hbeta atoms of cluster ligands C45 and C64 in the oxidized and reduced forms of T4MOC. In the reduced state, the Curie temperature dependence of the Hbeta protons corresponded to that predicted from the simple vector spin-coupling model for nuclei associated with the localized ferric site. The signal at 25.5 ppm in the 1H spectrum of reduced T4MOC was assigned on the basis of selective 2H labeling to the His Hepsilon1 atom of one of the cluster ligands (H47 or H67). This assignment was corroborated by a one bond 1H-13C correlation (at 25.39 ppm 1H and 136.11 ppm 13C) observed in spectra of [U-13C]T4MOC with a 1H-13C coupling constant of approximately 192 Hz. The carbon chemical shift and one bond coupling constant are those expected for 1Hepsilon1-13Cepsilon1 in the imidazolium ring of histidine and are inconsistent with values expected for cysteine 1Halpha-13Calpha. The His Hepsilon1 proton exhibited weak Curie temperature dependence from 283 to 303 K, contrary to the anti-Curie temperature dependence predicted from the spin coupling model for nuclei associated with the localized ferrous site. A 1H peak at -12.3 ppm was observed in spectra of reduced T4MOC; this signal was found to correspond to a hydrogen (probably in an H-bond to the cluster) that exchanged with solvent with a half-time of about 2 days in the oxidized state but with a much longer (undetectable) half-time in the reduced state. These results with T4MOC call into question certain 1H assignments recently reported on the basis of NOE measurements for the comparable Rieske ferredoxin component of an evolutionarily related alkene monooxygenase from Xanthobacter sp. Py2 [Holz, R. C., Small, F. J., and Ensign, S. A, (1997) Biochemistry 36, 14690-14696]. Selective 15N labeling was used to identify hyperfine-shifted 15N NMR signals from the backbone nitrogens of all four cluster ligands (C45, H47, C64, and H67), from the Nepsilon2 atoms of the two histidine ligands (H47 and H67), and from nonligand Gln and Ala residues (Q48 and A66) present in the cluster-binding motif of T4MOC in the oxidized and reduced states. The results indicate that the Ndelta1 of each of the two ligand histidines of T4MOC are ligated to an iron atom and reveal a pattern of H-bonding to the Rieske [2Fe-2S] center involving four (H47, Q48, A66, and H67 of T4MOC) of the five backbone amide H-bonds expected on the basis of comparison with the crystal structures of other related Rieske proteins; the fifth backbone amide (I50 of T4MOC) failed to exhibit a hyperfine shift. This anomaly may arise from the lack of an associated disulfide in T4MOC, a fundamental structural difference between the three types of Rieske proteins that may be related to functional diversity in this protein family.

Fractionation factors and activation energies for exchange of the low barrier hydrogen bonding proton in peptidyl trifluoromethyl ketone complexes of chymotrypsin.

NMR investigations have been carried out of complexes between bovine chymotrypsin Aalpha and a series of four peptidyl trifluoromethyl ketones, listed here in order of increasing affinity for chymotrypsin: N-Acetyl-L-Phe-CF3, N-Acetyl-Gly-L-Phe-CF3, N-Acetyl-L-Val-L-Phe-CF3, and N-Acetyl-L-Leu-L-Phe-CF3. The D/H fractionation factors (phi) for the hydrogen in the H-bond between His 57 and Asp 102 (His 57-Hdelta1) in these four complexes at 5 degreesC were in the range phi = 0.32-0.43, expected for a low-barrier hydrogen bond. For this series of complexes, measurements also were made of the chemical shifts of His 57-Hepsilon1 (delta2,2-dimethylsilapentane-5-sulfonic acid 8.97-9. 18), the exchange rate of the His 57-Hdelta1 proton with bulk water protons (284-12.4 s-1), and the activation enthalpies for this hydrogen exchange (14.7-19.4 kcal.mol-1). It was found that the previously noted correlations between the inhibition constants (Ki 170-1.2 microM) and the chemical shifts of His 57-Hdelta1 (delta2, 2-dimethylsilapentane-5-sulfonic acid 18.61-18.95) for this series of peptidyl trifluoromethyl ketones with chymotrypsin [Lin, J., Cassidy, C. S. & Frey, P. A. (1998) Biochemistry 37, 11940-11948] could be extended to include the fractionation factors, hydrogen exchange rates, and hydrogen exchange activation enthalpies. The results support the proposal of low barrier hydrogen bond-facilitated general base catalysis in the addition of Ser 195 to the peptidyl carbonyl group of substrates in the mechanism of chymotrypsin-catalyzed peptide hydrolysis. Trends in the enthalpies for hydrogen exchange and the fractionation factors are consistent with a strong, double-minimum or single-well potential hydrogen bond in the strongest complexes. The lifetimes of His 57-Hdelta1, which is solvent shielded in these complexes, track the strength of the hydrogen bond. Because these lifetimes are orders of magnitude shorter than those of the complexes themselves, the enzyme must have a pathway for hydrogen exchange at this site that is independent of dissociation of the complexes.

Practical model fitting approaches to the direct extraction of NMR parameters simultaneously from all dimensions of multidimensional NMR spectra.

A maximum likelihood (ML)-based approach has been established for the direct extraction of NMR parameters (e.g., frequency, amplitude, phase, and decay rate) simultaneously from all dimensions of a D-dimensional NMR spectrum. The approach, referred to here as HTFD-ML (hybrid time frequency domain maximum likelihood), constructs a time-domain model composed of a sum of exponentially-decaying sinusoidal signals. The apodized Fourier transform of this time-domain signal is a model spectrum that represents the 'best-fit' to the equivalent frequency-domain data spectrum. The desired amplitude and frequency parameters can be extracted directly from the signal model constructed by the HTFD-ML algorithm. The HTFD-ML approach presented here, as embodied in the software package CHIFIT, is designed to meet the challenges posed by model fitting of D-dimensional NMR data sets, where each consists of many data points (10(8) is not uncommon) encoding information about numerous signals (up to 10(5) for a protein of moderate size) that exhibit spectral overlap. The suitability of the approach is demonstrated by its application to the concerted analysis of a series of ten 2D 1H-15N HSQC experiments measuring 15N T1 relaxation. In addition to demonstrating the practicality of performing maximum likelihood analysis on large, multidimensional NMR spectra, the results demonstrate that this parametric model-fitting approach provides more accurate amplitude and frequency estimates than those obtained from conventional peak-based analysis of the FT spectrum. The improved performance of the model fitting approach derives from its ability to take into account the simultaneous contributions of all signals in a crowded spectral region (deconvolution) as well as to incorporate prior knowledge in constructing models to fit the data.