High-resolution experimental phases for tryptophanyl-tRNA synthetase (TrpRS) complexed with tryptophanyl-5'AMP.

Native data, anomalous data at three wavelengths and an independent peak-wavelength data set for SeMet-substituted protein have been collected from cryoprotected crystals of the TrpRS-adenylate product (TAM) complex to a resolution limit of 1.7 A. Independent phase sets were developed using SHARP and improved by solvent flipping with SOLOMON using molecular envelopes derived from experimental densities for, respectively, peak-wavelength SAD data from four different crystals, MAD data and their M(S)IRAS combinations with native data. Hendrickson-Lattman phase-probability coefficients from each phase set were used in BUSTER to drive maximum-likelihood refinements of well defined parts of the previously refined room-temperature 2.9 A structure. Maximum-entropy completion followed by manual rebuilding was then used to generate a model for the missing segments, bound ligand and solvent molecules. Surprisingly, peak-wavelength SAD experiments produced the smallest phase errors relative to the refined structures. Selenomethionylated models deviate from one another by 0.25 A and from the native model by 0.38 A, but all have r.m.s. deviations of approximately 1.0 A from the 2.9 A model. Difference Fourier calculations between amplitudes from the 300 K experiment and the new amplitudes at 100 K using 1.7 A model phases show no significant structural changes arising from temperature variation or addition of cryoprotectant. The main differences between low- and high-resolution structures arise from correcting side-chain rotamers in the core of the protein as well as on the surface. These changes improve various structure-validation criteria.

Four-body potentials reveal protein-specific correlations to stability changes caused by hydrophobic core mutations.

Mutational experiments show how changes in the hydrophobic cores of proteins affect their stabilities. Here, we estimate these effects computationally, using four-body likelihood potentials obtained by simplicial neighborhood analysis of protein packing (SNAPP). In this procedure, the volume of a known protein structure is tiled with tetrahedra having the center of mass of one amino acid side-chain at each vertex. Log-likelihoods are computed for the 8855 possible tetrahedra with equivalent compositions from structural databases and amino acid frequencies. The sum of these four-body potentials for tetrahedra present in a given protein yields the SNAPP score. Mutations change this sum by changing the compositions of tetrahedra containing the mutated residue and their related potentials. Linear correlation coefficients between experimental mutational stability changes, Delta(DeltaG(unfold)), and those based on SNAPP scoring range from 0.70 to 0.94 for hydrophobic core mutations in five different proteins. Accurate predictions for the effects of hydrophobic core mutations can therefore be obtained by virtual mutagenesis, based on changes to the total SNAPP likelihood potential. Significantly, slopes of the relation between Delta(DeltaG(unfold)) and DeltaSNAPP for different proteins are statistically distinct, and we show that these protein-specific effects can be estimated using the average SNAPP score per residue, which is readily derived from the analysis itself. This result enhances the predictive value of statistical potentials and supports previous suggestions that "comparable" mutations in different proteins may lead to different Delta(DeltaG(unfold)) values because of differences in their flexibility and/or conformational entropy.

Phase determination via Sayre-type equations with anomalous-scattering data.

The necessary background for the analysis of complex-valued electron-density maps is established. Various systems of structure-factor equations of convolutional type akin to Sayre's squaring method equations are tested for agreement on the real and imaginary parts of the electron density as well as approximations thereof. A system of convolutional structure-factor equations holding in a complex-valued electron density generated by two atom types is developed. The scope of application of these equations is determined and it is shown that the equations provide a method of extrapolating high-resolution phases from a low-resolution base phase set without introducing further model bias. Additional applications to phase refinement are explored.

Ancient adaptation of the active site of tryptophanyl-tRNA synthetase for tryptophan binding.

The amino acid binding domains of the tryptophanyl (TrpRS)- and tyrosyl-tRNA synthetases (TyrRS) of Bacillus stearothermophilus are highly homologous. These similarities suggest that conserved residues in TrpRS may be responsible for both determining tryptophan recognition and discrimination against tyrosine. This was investigated by the systematic mutation of TrpRS residues based upon the identity of homologous positions in TyrRS. Of the four residues which interact directly with the aromatic side chain of tryptophan (Phe5, Met129, Asp132, and Val141) replacements of Asp132 led to significant changes in the catalytic efficiency of Trp aminoacylation (200-1250-fold reduction in k(cat)/K(M)) and substitution of Val141 by the larger Glu side chain reduced k(cat)/K(M) by 300-fold. Mutation of Pro127, which determines the position of active-site residues, did not significantly effect Trp binding. Of the mutants tested, D132N TrpRS also showed a significant reduction in discrimination against Tyr, with Tyr acting as a competitive inhibitor but not a substrate. The analogous residue in B. stearothermophilusTyrRS (Asp176) has also been implicated as a determinant of amino acid specificity in earlier studies [de Prat Gay, G., Duckworth, H. W., and Fersht, A. R. (1993) FEBS Lett. 318, 167-171]. This striking similarity in the function of a highly conserved residue found in both TrpRS and TyrRS provides mechanistic support for a common origin of the two enzymes.

2.9 A crystal structure of ligand-free tryptophanyl-tRNA synthetase: domain movements fragment the adenine nucleotide binding site.

The crystal structure of ligand-free tryptophanyl-tRNA synthetase (TrpRS) was solved at 2.9 A using a combination of molecular replacement and maximum-entropy map/phase improvement. The dimeric structure (R = 23.7, Rfree = 26.2) is asymmetric, unlike that of the TrpRS tryptophanyl-5'AMP complex (TAM; Doublié S, Bricogne G, Gilmore CJ, Carter CW Jr, 1995, Structure 3:17-31). In agreement with small-angle solution X-ray scattering experiments, unliganded TrpRS has a conformation in which both monomers open, leaving only the tryptophan-binding regions of their active sites intact. The amino terminal alphaA-helix, TIGN, and KMSKS signature sequences, and the distal helical domain rotate as a single rigid body away from the dinucleotide-binding fold domain, opening the AMP binding site, seen in the TAM complex, into two halves. Comparison of side-chain packing in ligand-free TrpRS and the TAM complex, using identification of nonpolar nuclei (Ilyin VA, 1994, Protein Eng 7:1189-1195), shows that significant repacking occurs between three relatively stable core regions, one of which acts as a bearing between the other two. These domain rearrangements provide a new structural paradigm that is consistent in detail with the "induced-fit" mechanism proposed for TyrRS by Fersht et al. (Fersht AR, Knill-Jones JW, Beduelle H, Winter G, 1988, Biochemistry 27:1581-1587). Coupling of ATP binding determinants associated with the two catalytic signature sequences to the helical domain containing the presumptive anticodon-binding site provides a mechanism to coordinate active-site chemistry with relocation of the major tRNA binding determinants.

Patterned library analysis: a method for the quantitative assessment of hypotheses concerning the determinants of protein structure.

Site-directed mutagenesis and combinatorial libraries are powerful tools for providing information about the relationship between protein sequence and structure. Here we report two extensions that expand the utility of combinatorial mutagenesis for the quantitative assessment of hypotheses about the determinants of protein structure. First, we show that resin-splitting technology, which allows the construction of arbitrarily complex libraries of degenerate oligonucleotides, can be used to construct more complex protein libraries for hypothesis testing than can be constructed from oligonucleotides limited to degenerate codons. Second, using eglin c as a model protein, we show that regression analysis of activity scores from library data can be used to assess the relative contributions to the specific activity of the amino acids that were varied in the library. The regression parameters derived from the analysis of a 455-member sample from a library wherein four solvent-exposed sites in an alpha-helix can contain any of nine different amino acids are highly correlated (P < 0.0001, R(2) = 0. 97) to the relative helix propensities for those amino acids, as estimated by a variety of biophysical and computational techniques.

Cytidine deaminases from B. subtilis and E. coli: compensating effects of changing zinc coordination and quaternary structure.

Cytidine deaminase from E. coli is a dimer of identical subunits (M(r) = 31 540), each containing a single zinc atom. Cytidine deaminase from B. subtilis is a tetramer of identical subunits (M(r) = 14 800). After purification from an overexpressing strain, the enzyme from B. subtilis is found to contain a single atom of zinc per enzyme subunit by flame atomic absorption spectroscopy. Fluorescence titration indicates that each of the four subunits contains a binding site for the transition state analogue inhibitor 5-fluoro-3,4-dihydrouridine. A region of amino acid sequence homology, containing residues that are involved in zinc coordination in the enzyme from E. coli, strongly suggests that in the enzyme from B. subtilis, zinc is coordinated by the thiolate side chains of three cysteine residues (Cys-53, Cys-86, and Cys-89) [Song, B. H., and Neuhard, J. (1989) Mol. Gen. Genet. 216, 462-468]. This pattern of zinc coordination appears to be novel for a hydrolytic enzyme, and might be expected to reduce the reactivity of the active site substantially compared with that of the enzyme from E. coli (His-102, Cys-129, and Cys-132). Instead, the B. subtilis and E. coli enzymes are found to be similar in their activities, and also in their relative binding affinities for a series of structurally related inhibitors with binding affinities that span a range of 6 orders of magnitude. In addition, the apparent pK(a) value of the active site is shifted upward by less than 1 unit. Sequence alignments, together with model building, suggest one possible mechanism of compensation.

Mouse steroid sulfotransferases: substrate specificity and preliminary X-ray crystallographic analysis.

Three mouse cytosolic sulfotransferases were expressed in Escherichia coli cells in order to study their substrate specificities toward natural as well as synthetic steroid hormones. The Km and Vmax values confirmed the high substrate specificity of estrogen and hydroxysteroid sulfotransferases toward estradiol and dehydroepiandrosterone, respectively. In sharp contrast, the synthetic estrogen diethylstilbestrol was metabolized efficiently by both enzymes to its disulfate ester. These sulfotransferases display highly stereospecific sulfotransferase activity for sulfating only the trans-isomer of diethylstilbestrol. Crystals suitable for high-resolution structure determination of estrogen sulfotransferase were grown with polyethylene glycol. The crystals belong to the orthorhombic space group P2(1)2(1)2, and diffracted to 2.5 A.

Crystal structure of estrogen sulphotransferase.

The structure of estrogen sulphotransferase has been solved in the presence of inactive cofactor PAP and substrate 17 beta-estradiol. This structure reveals structural similarities between cytosolic sulphotransferases and nucleotide kinases.

The structure of the cytidine deaminase-product complex provides evidence for efficient proton transfer and ground-state destabilization.

Crystal structures of the cytidine deaminase-uridine product complex prepared either by cocrystallizing enzyme with uridine or by diffusing cytidine into ligand-free crystals show that the product binds as a 4-ketopyrimidine. They reveal four additional features of the catalytic process. (1) A water molecule bound to a site previously observed to bind the incoming 4-NH2 group represents the site for the leaving ammonia molecule. The conserved Pro 128 accommodates both moieties by orienting the carbonyl group of the previous residue. (2) The Glu 104 carboxylate group rotates from its hydrogen bond to the O4 hydroxyl group in transition-state analog complexes, forming a new hydrogen bond to the leaving group moiety. Thus, after stabilizing the hydroxyl group in the transition state, Glu 104 transfers a proton from that group to the leaving amino group, promoting enol-to-keto isomerization of the product. (3) Difference Fourier comparisons with transition-state complexes indicate that the pyrimidine ring rotates toward the zinc by approximately 10 degrees. The active site thus "pulls" the ring and 4-NH2 group in opposite directions during catalysis. To preserve coplanarity of the 4-keto group with the pyrimidine ring, the N1-C1' glycosidic bond bends by approximately 19 degrees out of the ring plane. This distortion may "spring-load" the product complex and promote dissociation. Failure to recognize a similar distortion could explain an earlier crystallographic interpretation of the adenosine deaminase-inosine complex [Wilson, D. K., & Quiocho, F. A. (1994) Nat. Struct. Biol. 1, 691-694]. (4) The Zn-Sgamma132 bond, which lengthens in transition-state complexes, shortens as the O4 atom returns to a state of lower negative charge in the planar product, consistent with our previous proposal that this bond buffers the zinc bond valence, compensating buildup of negative charge on the oxygen nucleophile during catalysis.

Chemical modifications of Bacillus subtilis tryptophanyl-tRNA synthetase.

A concerted conformational change in Bacillus subtilis tryptophanyl-tRNA synthetase (TrpRS) was evident from previous fluorescence on the quenching of the single Trp residue Trp-92 in the 4FTrp-AMP complexed enzyme. In this study, chemical modifications of the B. subtilis TrpRS were employed to further characterize this conformational change, with the single Trp residue serving as a marker for monitoring the change. Modifications of the enzyme by means of the Trp-specific agent N-bromosuccinimide (NBS) or 3-bromo-3-methyl-2-(2-nitrophenylmercapto)-3H-indole (BNPS-skatole) inactivated the enzyme in accord with the essential role of Trp-92, as identified previously by site-directed mutagenesis. ATP sensitized TrpRS toward inactivation by NBS and BNPS-skatole, which suggested a conformational change that resulted in greater accessibility of Trp-92 toward modifications. In contrast, the cognate tRNATrp substrate exerted a specific protective effect against inactivation by both of the reagents, indicating that the TrpRS-tRNATrp interaction reduces the accessibility of Trp-92 under our experimental conditions. By comparison, modification of sulfhydryl groups by means of iodoacetamide did not reduce TrpRS activity. Observations on Trp-specific modification and substrate protection effects are discussed in the context of the Bacillus stearothermophilus TrpRS crystal structure.

A local approximation to supersaturation affords a useful coordinate transformation for the study of crystal growth.

Two of the most important experimental variables in the search for appropriate crystallization conditions are the initial concentrations of macromolecule and crystallizing agent. Previously, it has been suggested that the coordinate transformation { [crystallizing agent], [macromolecule] } --> { [macromolecule] x [crystallizing agent], [macromolecule]} be used to sample crystal growth conditions. Here, it is shown that this transformation is a special case of a generally applicable transformation. The initial supersaturation can be represented locally by a rectangular hyperbola involving multiples of the product of macromolecule and crystallizing agent concentrations. The coordinate system for the solubility diagram, ([crystallizing agent] versus [macromolecule]), can thus be transformed analytically to an alternative coordinate system in which the independent variables are local approximations to the initial supersaturation and the reservoir of soluble macromolecule available to feed a growing seed. In the new coordinate system the 'nucleation zone' is 'orthogonalized', so it can be sampled efficiently on a rectangular grid, with greater assurance that experiments will give rise to crystals. Moreover, since these new coordinate directions segregate fundamental effects on nucleation from effects on growth, using them in experimental designs should improve data analysis for response surface experiments.

A concerted tryptophanyl-adenylate-dependent conformational change in Bacillus subtilis tryptophanyl-tRNA synthetase revealed by the fluorescence of Trp92.

A semi-conserved tryptophan residue of Bacillus subtilis tryptophanyl-tRNA synthetase (TrpRS) was previously asserted to be an essential residue and directly involved in tRNATrp binding and recognition. The crystal structure of the Bacillus stearothermophilus TrpRS tryptophanyl-5'-adenylate complex (Trp-AMP) shows that the corresponding Trp91 is buried and in the dimer interface, contrary to the expectations of the earlier assertation. Here we examine the role of this semi-conserved tryptophan residue using fluorescence spectroscopy. B. subtilis TrpRS has a single tryptophan residue, Trp92. 4-Fluorotryptophan (4FW) is used as a non-fluorescent substrate analog, allowing characterization of Trp92 fluorescence in the 4-fluorotryptophanyl-5'-adenylate (4FW-AMP) TrpRS complex. Complexation causes the Trp92 fluorescence to become quenched by 70%. Titrations, forming this complex under irreversible conditions, show that this quenching is essentially complete after half of the sites are filled. This indicates that a substrate-dependent mechanism exists for the inter-subunit communication of conformational changes. Trp92 fluorescence is not efficiently quenched by small solutes in either the apo- or complexed form. From this we conclude that this tryptophan residue is not solvent exposed and that binding of the Trp92 to tRNATrp is unlikely. Time-resolved fluorescence indicates conformational heterogeneity of B. subtilis Trp92 with the fluorescence decay being best described by three discrete exponential decay times. The decay-associated spectra (DAS) of the apo- and complexed-TrpRS show large variations of the concentration of individual fluorescence decay components. Based on recent correlations of these data with changes in the local secondary structure of the backbone containing the fluorescent tryptophan residue, we conclude that changes observed in Trp92 time-resolved fluorescence originate primarily from large perturbations of its local secondary structure. The quenching of Trp92 in the 4FW-AMP complex is best explained by the crystal structure conformation, in which the tryptophan residue is found in an alpha-helix. The amino acid residue cysteine is observed clearly within the quenching radius (3.6 angstroms) of the conserved tryptophan residue. These tryptophan and cysteine residues are neighbors, one helical turn apart. If this local alpha-helix was disrupted in the apo-TrpRS, this disruption would concomitantly relieve the putative cysteine quenching by separating the two residues. Hence we propose a substrate-dependent local helix-coil transition to explain both the observed time-resolved and steady-state fluorescence of Trp92. A mechanism can be further inferred for the inter-subunit communication involving the substrate ligand Asp132 and a small alpha-helix bridging the substrate tryptophan residue and the conserved tryptophan residue of the opposite subunit. This putative mechanism is also consistent with the observed pH dependence of TrpRS crystal growth and substrate binding. We observe that the mechanism of TrpRS has a dynamic component, and contend that conformational dynamics of aminoacyl-tRNA synthetases must be considered as part of the molecular basis for the recognition of cognate tRNA.

Incomplete factorial and response surface methods in experimental design: yield optimization of tRNA(Trp) from in vitro T7 RNA polymerase transcription.

We have studied the yield of Escherichia coli tRNA(Trp) obtained from in vitro T7 RNA polymerase transcription using incomplete factorial and response surface methods. Incomplete factorial experiments were first used to estimate the relative impact of six variables on the yield of tRNA(Trp). Fifteen trials were performed according to a balanced and randomized design. The correlation between observed yield and all experimental variables was identified by stepwise multiple linear regression analysis. The concentrations of T7 RNA polymerase, DNA template, NTP and MgCl2 proved to be significantly correlated with the yield of tRNA(Trp). We then optimized the yield with respect to each of these four variables simultaneously with a designed, response surface experiment based on the Hardin-Sloane minimum prediction variance algorithm. Twenty experiments were performed, in duplicate, to sample the quadratic surface relating the yield to the four significant variables. Coefficients of the quadratic function with all two-factor interactions were evaluated by stepwise regression using least squares, and significant coefficients were retained. Partial differentiation of the resulting quadratic model showed it to possess an optimum. Transcription performed at the corresponding conditions yielded 6-fold more tRNA(Trp) than the initial conditions, confirming the predictive value of the experimentally determined response surface.

Cytidine deaminase complexed to 3-deazacytidine: a "valence buffer" in zinc enzyme catalysis.

The cytidine deaminase substrate analog inhibitor 3-deazacytidine binds with its 4-amino group inserted into a site previously identified as a probable binding site for the leaving ammonia group. Binding to this site shifts the pyrimidine ring significantly further from the activated water molecule than the position it occupies in either of two complexes with compounds capable of hydrogen bonding at the 3-position of the ring [Xiang et al. (1995) Biochemistry 34, 4516-4523]. Difference Fourier maps between the deazacytidine, dihydrozebularine, and zebularine--hydrate inhibitor complexes suggest that the ring itself moves successively toward the activated water, leaving the amino group behind in this site as the substrate complex approaches the transition state. They also reveal systematic changes in a single zinc-sulfur bond distance. These correlate with chemical changes expected as the substrate approaches the tetrahedral transition state, in which the zinc-activated hydroxyl group develops maximal negative charge and forms a short hydrogen bond to the neighboring carboxylate group of Glu 104. Empirical bond valence relationships suggest that the Zn-S gamma 132 bond functions throughout the reaction as a "valence buffer" that accommodates changing negative charge on the hydroxyl group. Similar structural features in alcohol dehydrogenase suggest that analogous mechanisms may be a general feature of catalysis by zinc enzymes.

Representing stereochemical information in macromolecular electron-density distributions by multi-dimensional histograms.

Previous studies have demonstrated the value of ideal electron-density histograms as targets for the corresponding histograms of experimental electron-density maps. The electron-density histogram makes use of density values as independent objects, and no relationship between them is taken into account. Extension to include the relationships between neighboring density values leads naturally to a multi-dimensional histogram defined as the joint frequency of the density values and their higher order derivatives. We show here that the multi-dimensional histogram including additional dimensions composed of the gradient magnitude and Laplacian of the density is minimally dependent on molecular folding and packing, and captures substantially more stereochemical information than the conventional electron-density histogram. The gradient histogram appears to be much more sensitive to phase errors than the conventional electron-density histogram. Potential uses of the multi-dimensional histogram include improved targets for density modification and more reliable figures of merit for evaluating correct phases.

Transition-state selectivity for a single hydroxyl group during catalysis by cytidine deaminase.

Cytidine deaminase binds transition-state analog inhibitors approximately 10(7) times more tightly than corresponding 3,4-dihydro analogs containing a proton in place of the 4-hydroxyl group. X-ray crystal structures of complexes with the two matched inhibitors differ only near a "trapped" water molecule in the complex with the 3,4-dihydro analog, where contacts are substantially less favorable than those with the hydroxyl group of the transition-state analog. The hydrogen bond between the hydroxyl group and the Glu 104 carboxylate shortens in that complex, and may become a "low-barrier" hydrogen bond, since at the same time the bond between zinc and the Cys 132 thiolate ligand lengthens. These differences must therefore account for most of the differential binding affinity related to catalysis. Moreover, the trapped water molecule retains some of the binding energy stabilizing the hydroxyl group in the transition-state analog complex. To this extent, the ratio of binding affinities for the two compounds is smaller than the true contribution of the hydroxyl group, a conclusion with significant bearing on interpreting difference free energies derived from substituent effects arising from chemical modification and/or mutagenesis.