Structure of the type III secretion and substrate-binding domain of Yersinia YopH phosphatase.

Pathogenic strains of Yersinia deploy a type III secretion system to inject the potent tyrosine phosphatase YopH into host cells, where it dephosphorylates focal adhesion-associated substrates. The amino-terminal, non-catalytic domain of YopH is bifunctional; it is essential for the secretion and binding of the specific chaperone SycH, but also targets the catalytic domain to substrates in the infected cell. We describe the 2.2 A resolution crystal structure of residues 1-129 of YopH from Yersinia pseudotuberculosis. The amino-terminal alpha-helix (2-17), comprising the secretion signal, and beta-strand (24-28) of one molecule exchange with another molecule to form a domain-swapped dimer. Nuclear magnetic resonance (NMR) and gel filtration experiments demonstrated that YopH(1-129) could exist as a monomer and/or a dimer in solution. The topology of the dimer and the dynamics of a monomeric form in solution observed by NMR imply that YopH has the propensity to unfold partially. The dimer is probably not important physiologically, but may mimic how SycH binds to the exposed non-polar surfaces of a partially unfolded YopH. Phosphopeptide-induced perturbations in NMR chemical shifts define a substrate-binding surface on YopH(1-129) that includes residues previously shown by mutagenesis to be essential for YopH function.

Delineation of the allosteric mechanism of a cytidylyltransferase exhibiting negative cooperativity.

The dimeric enzyme CTP:glycerol-3-phosphate cytidylyltransferase (GCT) displays strong negative cooperativity between the first and second binding of its substrate, CTP. Using NMR to study the allosteric mechanism of this enzyme, we observe widespread chemical shift changes for the individual CTP binding steps. Mapping these changes onto the molecular structure allowed the formulation of a detailed model of allosteric conformational change. Upon the second step of ligand binding, NMR experiments indicate an extensive loss of conformational exchange broadening of the backbone resonances of GCT. This suggests that a fraction of the free energy of negative cooperativity is entropic in origin.

Mapping the interactions between flavodoxin and its physiological partners flavodoxin reductase and cobalamin-dependent methionine synthase.

Flavodoxins are electron-transfer proteins that contain the prosthetic group flavin mononucleotide. In Escherichia coli, flavodoxin is reduced by the FAD-containing protein NADPH:ferredoxin (flavodoxin) oxidoreductase; flavodoxins serve as electron donors in the reductive activation of anaerobic ribonucleotide reductase, biotin synthase, pyruvate formate lyase, and cobalamin-dependent methionine synthase. In addition, domains homologous to flavodoxin are components of the multidomain flavoproteins cytochrome P450 reductase, nitric oxide synthase, and methionine synthase reductase. Although three-dimensional structures are known for many of these proteins and domains, very little is known about the structural aspects of their interactions. We address this issue by using NMR chemical shift mapping to identify the surfaces on flavodoxin that bind flavodoxin reductase and methionine synthase. We find that these physiological partners bind to unique overlapping sites on flavodoxin, precluding the formation of ternary complexes. We infer that the flavodoxin-like domains of the cytochrome P450 reductase family form mutually exclusive complexes with their electron-donating and -accepting partners, complexes that require conformational changes for interconversion.

Functional dynamics in the active site of the ribonuclease binase.

Binase, a member of a family of microbial guanyl-specific ribonucleases, catalyzes the endonucleotic cleavage of single-stranded RNA. It shares 82% amino acid identity with the well-studied protein barnase. We used NMR spectroscopy to study the millisecond dynamics of this small enzyme, using several methods including the measurement of residual dipolar couplings in solution. Our data show that the active site of binase is flanked by loops that are flexible at the 300-micros time scale. One of the catalytic residues, His-101, is located on such a flexible loop. In contrast, the other catalytic residue, Glu-72, is located on a beta-sheet, and is static. The residues Phe-55, part of the guanine base recognition site, and Tyr-102, stabilizing the base, are the most dynamic. Our findings suggest that binase possesses an active site that has a well-defined bottom, but which has sides that are flexible to facilitate substrate access/egress, and to deliver one of the catalytic residues. The motion in these loops does not change on complexation with the inhibitor d(CGAG) and compares well with the maximum k(cat) (1,500 s(-1)) of these ribonucleases. This observation indicates that the NMR-measured loop motions reflect the opening necessary for product release, which is apparently rate limiting for the overall turnover.

An iterative fitting procedure for the determination of longitudinal NMR cross-correlation rates.

We present a method to measure (15)N-(1)H dipolar/(15)N CSA longitudinal cross-correlation rates in protonated proteins. The method depends on the measurement of four observables: the cumulative proton-proton cross relaxation rates, the (15)N R(1) relaxation rate, the multiexponential decay of 2N(Z)H(N)(Z) spin-order, and multiexponential buildup of 2N(Z)H(N)(Z) spin-order. The (15)N-(1)H dipolar/(15)N CSA longitudinal cross-correlation rate is extracted from these measurements by an iterative fitting procedure to the solution of differential equations describing the coupled relaxation dynamics of the z-magnetization of the (15)N nucleus, the two-spin-order 2N(Z)H(N)(Z), and a two-spin-order term 2N(Z)H(Q)(Z) describing the interaction with remote protons. The method is applied to the microbial ribonuclease binase. The method can also extract longitudinal cross-correlation rates for those amide protons that are involved in rapid solvent exchange. The experiment that serves for extracting proton-proton cross-relaxation rates is a modification of 3D (15)N-resolved NOESY-HSQC. The experiment restores the solvent magnetization to its equilibrium state during data detection for all phase cycling steps and all values of NOE mixing times and is recommended for use in standard applications as well.

Magnetization transfer via residual dipolar couplings: application to proton-proton correlations in partially aligned proteins.

A novel three-dimensional NMR experiment is reported that allows the observation of correlations between amide and other protons via residual dipolar couplings in partially oriented proteins. The experiment is designed to permit quantitative measurement of the magnitude of proton-proton residual dipolar couplings in larger molecules and at higher degree of alignments. The observed couplings contain data valuable for protein resonance assignment, local protein structure refinement, and determination of low-resolution protein folds.

Structural insights into substrate binding by the molecular chaperone DnaK.

How substrate affinity is modulated by nucleotide binding remains a fundamental, unanswered question in the study of 70 kDa heat shock protein (Hsp70) molecular chaperones. We find here that the Escherichia coli Hsp70, DnaK, lacking the entire alpha-helical domain, DnaK(1-507), retains the ability to support lambda phage replication in vivo and to pass information from the nucleotide binding domain to the substrate binding domain, and vice versa, in vitro. We determined the NMR solution structure of the corresponding substrate binding domain, DnaK(393-507), without substrate, and assessed the impact of substrate binding. Without bound substrate, loop L3,4 and strand beta3 are in significantly different conformations than observed in previous structures of the bound DnaK substrate binding domain, leading to occlusion of the substrate binding site. Upon substrate binding, the beta-domain shifts towards the structure seen in earlier X-ray and NMR structures. Taken together, our results suggest that conformational changes in the beta-domain itself contribute to the mechanism by which nucleotide binding modulates substrate binding affinity.

Secondary structure and fold homology of the ArsC protein from the Escherichia coli arsenic resistance plasmid R773.

Resistance to several toxic anions in Escherichia coli is conferred by the ars operon carried on plasmid R773. The gene products of this operon catalyze extrusion of antimonials and arsenicals from cells. In this paper, we report the determination of the overall fold for ArsC, a 16 kDa protein of the ars operon involved in the reduction of arsenate to arsenite, using multidimensional, multinuclear NMR. The protein is found to contain large regions of extensive mobility, particularly in the active site. A model fold, computed on the basis of a preliminary set of NOEs, was found to be structurally homologous to E. coli glutaredoxin, thiol transferases, and glutathione S-transferase. Some kinship to the structure of low molecular weight tyrosine phosphatases, based on rough topological similarity but more so on the basis of a common anion-binding-loop motif H-CX(n)R, was also detected. Although functional, secondary, and tertiary structural homology is observed with these molecules, no significant homology in primary structure was detected. The mobilities of the active site of ArsC and of other enzymes are discussed.

High-resolution four-dimensional HMQC-NOESY-HSQC spectroscopy.

Practical optimization of the 4D [(1)H, (13)C, (13)C, (1)H] HMQC-NOESY-HSQC experiment in terms of distribution of resolution over the indirect dimensions is analyzed in detail. Recommendations for an optimal experiment are based on computer simulations assessing the effective resolution of the experiment, defined as the percentage of all possible NOE cross peaks that can be assigned unambiguously on the basis of the spectral data alone. Using actual (13)C-(1)H spectra of an 18-kDa chaperone protein, the analysis shows that experiments with the best effective resolution are also among the most sensitive ones. When combined with an efficient aliasing scheme that reduces indirect spectral space 124-fold, a 4D experiment that yields unambiguous assignments for 41% of all possible NOE cross peaks can be recorded in 28 h. A high-resolution experiment, which can be recorded in 8 days, yields 61% unambiguous assignments and can be analyzed more easily using standard NMR display software. The predictions are verified with experimental 4D spectra from which 1850 NOEs (914 long-range) were extracted for the 18-kDa chaperone protein.

High-resolution solution structure of the 18 kDa substrate-binding domain of the mammalian chaperone protein Hsc70.

The three-dimensional structure for the substrate-binding domain of the mammalian chaperone protein Hsc70 of the 70 kDa heat shock class (HSP70) is presented. This domain includes residues 383-540 (18 kDa) and is necessary for the binding of the chaperone with substrate proteins and peptides. The high-resolution NMR solution structure is based on 4150 experimental distance constraints leading to an average root-mean-square precision of 0.38 A for the backbone atoms and 0.76 A for all atoms in the beta-sandwich sub-domain. The protein is observed to bind residue Leu539 in its hydrophobic substrate-binding groove by intramolecular interaction. The position of a helical latch differs dramatically from what is observed in the crystal and solution structures of the homologous prokaryotic chaperone DnaK. In the Hsc70 structure, the helix lies in a hydrophobic groove and is anchored by a buried salt-bridge. Residues involved in this salt-bridge appear to be important for the allosteric functioning of the protein. A mechanism for interdomain allosteric modulation of substrate-binding is proposed. It involves large-scale movements of the helical domain, redefining the location of the hinge area that enables such motions.

Line narrowing in spectra of proteins dissolved in a dilute liquid crystalline phase by band-selective adiabatic decoupling: application to 1HN-15N residual dipolar coupling measurements.

Residual heteronuclear dipolar couplings obtained from partially oriented protein samples can provide unique NMR constraints for protein structure determination. However, partial orientation of protein samples also causes severe 1H line broadening resulting from residual 1H-1H dipolar couplings. In this communication we show that band-selective 1H homonuclear decoupling during data acquisition is an efficient way to suppress residual 1H-1H dipolar couplings, resulting in spectra that are still amenable to solution NMR analysis, even with high degrees of alignment. As an example, we present a novel experiment with improved sensitivity for the measurement of one-bond 1HN-15N residual dipolar couplings in a protein sample dissolved in magnetically aligned liquid crystalline bicelles.

Characterizing semilocal motions in proteins by NMR relaxation studies.

The understanding of protein function is incomplete without the study of protein dynamics. NMR spectroscopy is valuable for probing nanosecond and picosecond dynamics via relaxation studies. The use of 15N relaxation to study backbone dynamics has become virtually standard. Here, we propose to measure the relaxation of additional nuclei on each peptide plane allowing for the observation of anisotropic local motions. This allows the nature of local motions to be characterized in proteins. As an example, semilocal rotational motion was detected for part of a helix of the protein Escherichia coli flavodoxin.

High-resolution detection of five frequencies in a single 3D spectrum: HNHCACO--a bidirectional coherence transfer experiment.

A new triple-resonance pulse sequence, 3D HNHCACO, is introduced and discussed, which identifies sequential correlations of the backbone nuclei (H alpha (i-1), C alpha (i-1), C(i-1), NH(i), N(i)) of doubly labeled proteins in H2O. The three-dimensional (3D) method utilizes a recording of 15N and 13C resonances in a single indirect time domain, the 13C' resonance in another indirect time domain, and detects both NH and H alpha protons. A bidirectional coherence transfer (NH(i) <--> N(i) <--> C(i-1) <--> C alpha (i-1) <--> H alpha (i-1)) is effectuated, resulting in a single high-resolution 3D spectrum that contains the frequencies of all five backbone nuclei. The experiment was applied to the 12.3 kDa ribonuclease from Bacillus intermedius (Binase).

NMR solution structure of the 21 kDa chaperone protein DnaK substrate binding domain: a preview of chaperone-protein interaction.

The solution structure of the 21 kDa substrate-binding domain of the Escherichia coli Hsp70-chaperone protein DnaK (DnaK 386-561) has been determined to a precision of 1.00 A (backbone of the beta-domain) from 1075 experimental restraints obtained from multinuclear, multidimensional NMR experiments. The domain is observed to bind to its own C-terminus and offers a preview of the interaction of this chaperone with other proteins. The bound protein region is tightly held at a single amino acid position (a leucyl residue) that is buried in a deep pocket lined with conserved hydrophobic residues. A second hydrophobic binding site was identified using paramagnetically labeled peptides. It is located in a region close to the N-terminus of the domain and may constitute the allosteric region that links substrate-binding affinity with nucleotide binding in the Hsp70 chaperones.

Structural analysis of the P10/11-P12 RNA domain of yeast RNase P RNA and its interaction with magnesium.

The P10/11-P12 RNA domain of yeast RNase P contains several highly conserved nucleotides within a conserved secondary structure. This RNA domain is essential for enzyme function in vivo, where it has a demonstrated role in divalent cation utilization. To better understand the function of this domain, its structure and alterations in response to magnesium have been investigated in vitro. A secondary structure model of the P10/11-P12 RNA domain had been previously developed by phylogenetic analysis. Computer modeling and energy minimization were applied to the Saccharomyces cerevisiae P10/11-P12 domain to explore alternatives and additional interactions not predicted by the phylogenetic consensus. The working secondary structure models were challenged with data obtained from 1H NMR and in vitro chemical and enzymatic probing experiments. The solution structure of the isolated domain was found to conform to the phylogenetic prediction within the context of the holoenzyme. Structure probing data also discriminated among additional base contacts predicted by energy minimization. The withdrawal of magnesium does not appear to cause gross refolding or rearrangement of the RNA domain structure. Instead, subtle changes occur in the solution accessibility of specific nucleotide positions. Most of the conserved nucleotides reported to be involved in magnesium utilization in vivo also display magnesium-dependent changes in vitro.

Identification of functional conserved residues of CTP:glycerol-3-phosphate cytidylyltransferase. Role of histidines in the conserved HXGH in catalysis.

The CTP:glycerol-3-phosphate cytidylyltransferase (GCT) of Bacillus subtilis has been shown to be similar in primary structure to the CTP:phosphocholine cytidylyltransferases of several organisms. To identify the residues of this cytidylyltransferase family that function in catalysis, the conserved hydrophilic amino acid residues plus a conserved tryptophan of the GCT were mutated to alanine. The most dramatic losses in activity occurred with H14A and H17A; these histidine residues are part of an HXGH sequence similar to that found in class I aminoacyl-tRNA synthetases. The kcat values for H14A and H17A were decreased by factors of 5 x 10(-5) and 4 x 10(-4), respectively, with no significant change in Km values. Asp-11, which is found near the HXGH sequence in the cytidylyltransferases but not aminoacyl-tRNA synthetases, was also important for activity, with the D11A mutation decreasing activity by a factor of 2 x 10(-3). Several residues found in the sequence RTEGISTT, a signature sequence for this cytidylyltransferase family, as well as other isolated residues were also shown to be important for activity, with kcat values decreasing by factors of 0.14-4 x 10(-4). The Km values of three mutant enzymes, D38A, W74A, and D94A, for both CTP and glycerol-3-phosphate were 6-130-fold higher than that of the wild-type enzyme. Mutant enzymes were analyzed by two-dimensional NMR to determine if the overall structures of the enzymes were intact. One of the mutant enzymes, D66A, was defective in overall structure, but several of the others, including H14A and H17A, were not. These results indicate that His-14 and His-17 play a role in catalysis and suggest that their role is similar to the role of the His residues in the HXGH sequence in class I aminoacyl-tRNA synthetases, i.e. to stabilize a pentacoordinate transition state.

Protein heteronuclear NMR assignments using mean-field simulated annealing.

A computational method for the assignment of the NMR spectra of larger (21 kDa) proteins using a set of six of the most sensitive heteronuclear multidimensional nuclear magnetic resonance experiments is described. Connectivity data obtained from HNC alpha, HN(CO)C alpha, HN(C alpha)H alpha, and H alpha (C alpha CO)NH and spin-system identification data obtained from CP-(H)CCH-TOCSY and CP-(H)C(C alpha CO)NH-TOCSY were used to perform sequence-specific assignments using a mean-field formalism and simulated annealing. This mean-field method reports the resonance assignments in a probabilistic fashion, displaying the certainty of assignments in an unambiguous and quantitative manner. This technique was applied to the NMR data of the 172-residue peptide-binding domain of the E. coli heat-shock protein, DnaK. The method is demonstrated to be robust to significant amounts of missing, spurious, noisy, extraneous, and erroneous data.