High throughput crystallography of TB drug targets.

Tuberculosis (TB) infects one-third of the world population. Despite 50 years of available drug treatments, TB continues to increase at a significant rate. The failure to control TB stems in part from the expense of delivering treatment to infected individuals and from complex treatment regimens. Incomplete treatment has fueled the emergence of multi-drug resistant (MDR) strains of Mycobacterium tuberculosis (Mtb). Reducing non-compliance by reducing the duration of chemotherapy will have a great impact on TB control. The development of new drugs that either kill persisting organisms, inhibit bacilli from entering the persistent phase, or convert the persistent bacilli into actively growing cells susceptible to our current drugs will have a positive effect. We are taking a multidisciplinary approach that will identify and characterize new drug targets that are essential for persistent Mtb. Targets are exposed to a battery of analyses including microarray experiments, bioinformatics, and genetic techniques to prioritize potential drug targets from Mtb for structural analysis. Our core structural genomics pipeline works with the individual laboratories to produce diffraction quality crystals of targeted proteins, and structural analysis will be completed by the individual laboratories. We also have capabilities for functional analysis and the virtual ligand screening to identify novel inhibitors for target validation. Our overarching goals are to increase the knowledge of Mtb pathogenesis using the TB research community to drive structural genomics, particularly related to persistence, develop a central repository for TB research reagents, and discover chemical inhibitors of drug targets for future development of lead compounds.

SAD phasing by combination of direct methods with the SOLVE/RESOLVE procedure.

In the initial stage of SAD phasing, the essential point is to break the intrinsic phase ambiguity. The presence of two kinds of phase information enables the discrimination of phase doublets from SAD data prior to density modification. One is from the heavy atoms (anomalous scatterers), while the other is from the direct-methods phase relationships. The former can be expressed by the Sim distribution, while the latter can be expressed by the Cochran distribution. Typically, only the Sim distribution has been used to yield initial phases for subsequent density modification. However, it has been demonstrated that using direct-methods phases based on the product of the Sim and Cochran distributions can lead to improved initial phases. In this paper, the direct-methods phasing procedure OASIS has been improved and combined with the SOLVE/RESOLVE procedure. Experimental SAD data from three known proteins with expected Bijvoet ratios / in the range 1.4-7.0% were used as test cases. In all cases, the phases obtained using the program RESOLVE beginning with initial phases based on experimental phases plus Sim and direct-methods information were more accurate than those based on experimental plus Sim phase information alone.

The TB structural genomics consortium: a resource for Mycobacterium tuberculosis biology.

The TB Structural Genomics Consortium is an organization devoted to encouraging, coordinating, and facilitating the determination and analysis of structures of proteins from Mycobacterium tuberculosis. The Consortium members hope to work together with other M. tuberculosis researchers to identify M. tuberculosis proteins for which structural information could provide important biological information, to analyze and interpret structures of M. tuberculosis proteins, and to work collaboratively to test ideas about M. tuberculosis protein function that are suggested by structure or related to structural information. This review describes the TB Structural Genomics Consortium and some of the proteins for which the Consortium is in the progress of determining three-dimensional structures.

Maximum-likelihood density modification using pattern recognition of structural motifs.

The likelihood-based approach to density modification [Terwilliger (2000), Acta Cryst. D56, 965-972] is extended to include the recognition of patterns of electron density. Once a region of electron density in a map is recognized as corresponding to a known structural element, the likelihood of the map is reformulated to include a term that reflects how closely the map agrees with the expected density for that structural element. This likelihood is combined with other aspects of the likelihood of the map, including the presence of a flat solvent region and the electron-density distribution in the protein region. This likelihood-based pattern-recognition approach was tested using the recognition of helical segments in a largely helical protein. The pattern-recognition method yields a substantial phase improvement over both conventional and likelihood-based solvent-flattening and histogram-matching methods. The method can potentially be used to recognize any common structural motif and incorporate prior knowledge about that motif into density modification.

Map-likelihood phasing.

The recently developed technique of maximum-likelihood density modification [Terwilliger (2000), Acta Cryst. D56, 965-972] allows a calculation of phase probabilities based on the likelihood of the electron-density map to be carried out separately from the calculation of any prior phase probabilities. Here, it is shown that phase-probability distributions calculated from the map-likelihood function alone can be highly accurate and that they show minimal bias towards the phases used to initiate the calculation. Map-likelihood phase probabilities depend upon expected characteristics of the electron-density map, such as a defined solvent region and expected electron-density distributions within the solvent region and the region occupied by a macromolecule. In the simplest case, map-likelihood phase-probability distributions are largely based on the flatness of the solvent region. Though map-likelihood phases can be calculated without prior phase information, they are greatly enhanced by high-quality starting phases. This leads to the technique of prime-and-switch phasing for removing model bias. In prime-and-switch phasing, biased phases such as those from a model are used to prime or initiate map-likelihood phasing, then final phases are obtained from map-likelihood phasing alone. Map-likelihood phasing can be applied in cases with solvent content as low as 30%. Potential applications of map-likelihood phasing include unbiased phase calculation from molecular-replacement models, iterative model building, unbiased electron-density maps for cases where 2F(o) - F(c) or sigma(A)-weighted maps would currently be used, structure validation and ab initio phase determination from solvent masks, non-crystallographic symmetry or other knowledge about expected electron density.

Solution structure of Pyrobaculum aerophilum DsrC, an archaeal homologue of the gamma subunit of dissimilatory sulfite reductase.

The solution structure of DsrC, an archaeal homologue of the gamma subunit of dissimilatory sulfite reductase, has been determined by NMR spectroscopy. This 12.7-kDa protein from the hyperthermophilic archaeon Pyrobaculum aerophilum adopts a novel fold consisting of an orthogonal helical bundle with a beta hairpin along one side. A portion of the structure resembles the helix-turn-helix DNA-binding motif common in transcriptional regulator proteins. The protein contains two disulfide bonds but remains folded following reduction of the disulfides. DsrC proteins from organisms other than Pyrobaculum species do not contain these disulfide bonds. A conserved cysteine next to the C-terminus, which is not involved in the disulfide bonds, is located on a seven-residue C-terminal arm that is not part of the globular protein and is likely to dynamically sample more than one conformation.

Ff gene 5 protein has a high binding affinity for single-stranded phosphorothioate DNA.

The gene 5 protein (g5p) of Ff bacteriophages is a well-studied model ssDNA-binding protein that binds cooperatively to the Ff ssDNA genome and single-stranded polynucleotides. Its affinity, K omega (the intrinsic binding constant times a cooperativity factor), can differ by several orders of magnitude for ssDNAs of different nearest-neighbor base compositions [Mou, T. C., Gray, C. W., and Gray, D. M. (1999) Biophys. J. 76, 1537-1551]. We found that the DNA backbone can also dramatically affect the binding affinity. The K omega for binding phosphorothioate-modified S-d(A)(36) was >300-fold higher than for binding unmodified P-d(A)(36) at 0.2 M NaCl. CD titrations showed that g5p bound phosphorothioate-modified oligomers with the same stoichiometry as unmodified oligomers. The CD spectrum of S-d(A)(36) underwent the same qualitative change upon protein binding as did the spectrum of unmodified DNA, and the phosphorothioate-modified DNA appeared to bind in the normal g5p binding site. Oligomers of d(A)(36) with different proportions of phosphorothioate nucleotides had binding affinities and CD perturbations intermediate to those of the fully modified and unmodified sequences. The influence of phosphorothioation on binding affinity was nearly proportional to the extent of the modification, with a small nearest-neighbor dependence. These and other results using d(ACC)(12) oligomers and mutant proteins indicated that the increased binding affinity of g5p for phosphorothioate DNA was not a polyelectrolyte effect and probably was not an effect due to the altered nucleic acid structure, but was more likely a general effect of the properties of the sulfur in the context of the phosphorothioate group.

Structural genomics in North America.

Structural genomics in North America has moved remarkably quickly from ideas to pilot projects. Just three years ago, the field was only a concept, independently being discussed by its many inventors. Now it is already a well-organized, increasingly-funded, consortium-based effort to determine protein structures on a large scale.

Maximum-likelihood density modification.

A likelihood-based approach to density modification is developed that can be applied to a wide variety of cases where some information about the electron density at various points in the unit cell is available. The key to the approach consists of developing likelihood functions that represent the probability that a particular value of electron density is consistent with prior expectations for the electron density at that point in the unit cell. These likelihood functions are then combined with likelihood functions based on experimental observations and with others containing any prior knowledge about structure factors to form a combined likelihood function for each structure factor. A simple and general approach to maximizing the combined likelihood function is developed. It is found that this likelihood-based approach yields greater phase improvement in model and real test cases than either conventional solvent flattening and histogram matching or a recent reciprocal-space solvent-flattening procedure [Terwilliger (1999), Acta Cryst. D55, 1863-1871].

Haloalkane dehalogenases: structure of a Rhodococcus enzyme.

The hydrolytic haloalkane dehalogenases are promising bioremediation and biocatalytic agents. Two general classes of dehalogenases have been reported from Xanthobacter and Rhodococcus. While these enzymes share 30% amino acid sequence identity, they have significantly different substrate specificities and halide-binding properties. We report the 1.5 A resolution crystal structure of the Rhodococcus dehalogenase at pH 5.5, pH 7.0, and pH 5.5 in the presence of NaI. The Rhodococcus and Xanthobacter enzymes have significant structural homology in the alpha/beta hydrolase core, but differ considerably in the cap domain. Consistent with its broad specificity for primary, secondary, and cyclic haloalkanes, the Rhodococcus enzyme has a substantially larger active site cavity. Significantly, the Rhodococcus dehalogenase has a different catalytic triad topology than the Xanthobacter enzyme. In the Xanthobacter dehalogenase, the third carboxylate functionality in the triad is provided by D260, which is positioned on the loop between beta7 and the penultimate helix. The carboxylate functionality in the Rhodococcus catalytic triad is donated from E141. A model of the enzyme cocrystallized with sodium iodide shows two iodide binding sites; one that defines the normal substrate and product-binding site and a second within the active site region. In the substrate and product complexes, the halogen binds to the Xanthobacter enzyme via hydrogen bonds with the N(eta)H of both W125 and W175. The Rhodococcusenzyme does not have a tryptophan analogous to W175. Instead, bound halide is stabilized with hydrogen bonds to the N(eta)H of W118 and to N(delta)H of N52. It appears that when cocrystallized with NaI the Rhodococcus enzyme has a rare stable S-I covalent bond to S(gamma) of C187.

In vivo screening of haloalkane dehalogenase mutants.

Haloalkane dehalogenase (Dh1A) from Xanthobacter autotrophicus GJ10 catalyzes the dehalogenation of short chain primary alkyl halides. Due to the high Km and low turnover, wild type Dh1A is not optimal for applications in bioremediation. We have developed an in vivo screen, based on a colorimetric pH indicator, to identify Dh1A mutant with improved catalytic activity. After screening 50,000 colonies, we identified a Dh1A mutant with a lower pH optimum. Sequence analysis of the mutant revealed a single substitution, alanine 149 to threonine, which is located close to the active site of Dh1A. Replacement of alanine 149 via site-directed mutagenesis with threonine, serine or cysteine retained the mutant phenotype. Other substitutions at position 149 show little or no activity.

Reciprocal-space solvent flattening.

Solvent flattening is a powerful tool for improving crystallographic phases for macromolecular structures obtained at moderate resolution, but uncertainties in the optimal weighting of experimental phases and modified phases make it difficult to extract all the phase information possible. Solvent flattening is essentially an iterative method for maximizing a likelihood function which consists of (i) experimental phase information and (ii) information on the likelihood of various arrangements of electron density in a map, but the likelihood function is generally not explicitly defined. In this work, a procedure is described for reciprocal-space maximization of a likelihood function based on experimental phases and characteristics of the electron-density map. The procedure can readily be applied to phase improvement based on solvent flattening and can potentially incorporate information on a wide variety of other characteristics of the electron-density map.

Evaluation of macromolecular electron-density map quality using the correlation of local r.m.s. density.

It has recently been shown that the standard deviation of local r.m. s. electron density is a good indicator of the presence of distinct regions of solvent and protein in macromolecular electron-density maps [Terwilliger & Berendzen (1999). Acta Cryst. D55, 501-505]. Here, it is demonstrated that a complementary measure, the correlation of local r.m.s. density in adjacent regions on the unit cell, is also a good measure of the presence of distinct solvent and protein regions. The correlation of local r.m.s. density is essentially a measure of how contiguous the solvent (and protein) regions are in the electron-density map. This statistic can be calculated in real space or in reciprocal space and has potential uses in evaluation of heavy-atom solutions in the MIR and MAD methods as well as for evaluation of trial phase sets in ab initio phasing procedures.

Rapid protein-folding assay using green fluorescent protein.

Formation of the chromophore of green fluorescent protein (GFP) depends on the correct folding of the protein. We constructed a "folding reporter" vector, in which a test protein is expressed as an N-terminal fusion with GFP. Using a test panel of 20 proteins, we demonstrated that the fluorescence of Escherichia coli cells expressing such GFP fusions is related to the productive folding of the upstream protein domains expressed alone. We used this fluorescent indicator of protein folding to evolve proteins that are normally prone to aggregation during expression in E. coli into closely related proteins that fold robustly and are fully soluble and functional. This approach to improving protein folding does not require functional assays for the protein of interest and provides a simple route to improving protein folding and expression by directed evolution.

Sigma2R, a reciprocal-space measure of the quality of macromolecular electron-density maps.

It has previously been shown that the presence of distinct regions of solvent and protein in macromolecular crystals leads to a high value of the standard deviation of local r.m.s. electron density and that this can in turn be used as a reliable measure of the quality of macromolecular electron-density maps [Terwilliger & Berendzen (1999a). Acta Cryst. D55, 501-505]. Here, it is demonstrated that a similar measure, sigmaR2, the variance of the local roughness of the electron density, can be calculated in reciprocal space. The formulation is suitable for rapid evaluation of macromolecular crystallographic phases, for phase improvement and for ab initio phasing procedures.

Automated MAD and MIR structure solution.

Obtaining an electron-density map from X-ray diffraction data can be difficult and time-consuming even after the data have been collected, largely because MIR and MAD structure determinations currently require many subjective evaluations of the qualities of trial heavy-atom partial structures before a correct heavy-atom solution is obtained. A set of criteria for evaluating the quality of heavy-atom partial solutions in macromolecular crystallography have been developed. These have allowed the conversion of the crystal structure-solution process into an optimization problem and have allowed its automation. The SOLVE software has been used to solve MAD data sets with as many as 52 selenium sites in the asymmetric unit. The automated structure-solution process developed is a major step towards the fully automated structure-determination, model-building and refinement procedure which is needed for genomic scale structure determinations.

Discrimination of solvent from protein regions in native Fouriers as a means of evaluating heavy-atom solutions in the MIR and MAD methods.

An automated examination of the native Fourier is tested as a means of evaluation of a heavy-atom solution in MAD and MIR methods for macromolecular crystallography. It is found that the presence of distinct regions of high and low density variation in electron-density maps is a good indicator of the correctness of a heavy-atom solution in the MIR and MAD methods. The method can be used to evaluate heavy-atom solutions during MAD and MIR structure solutions and to determine the handedness of the structure if anomalous data have been measured.

Exploring structure space. A protein structure initiative.

The genome projects are changing biology by providing the genetic blueprints of entire organisms. The blueprints are tantalizing but we cannot deduce everything we need to know from them, including the structures and detailed functions of proteins. In this paper we describe an approach for obtaining structural information about proteins on a genomic scale. We describe how structural and functional information might eventually be put together to form a basis for describing life at many levels. We then describe how structural information fits into this picture and classes of proteins for which structural information would be useful in a genomic context. We conclude with a proposal for an initiative to determine protein structures on a very large scale.

A nondenaturing purification scheme for the DNA-binding domain of poly(ADP-ribose) polymerase, a structure-specific DNA-binding protein.

Poly(ADP-ribose) polymerase (PARP) is thought to be involved in DNA repair given its ability to recognize and bind to DNA strand breaks. During apoptosis, PARP is proteolytically cleaved into two stable fragments, the N-terminal 25-kDa DNA-binding domain (DBD) and the 85-kDa fragment containing the automodification and catalytic domains. To understand the DNA-binding properties of PARP, we expressed a recombinant hexahistidine tagged protein (His-DBD) in Escherichia coli. We modified expression to facilitate protein folding by including zinc and reducing the induction temperature. Properly folded, the DNA-binding domain of PARP binds to single- and double-stranded DNA in a structure-specific manner. To eliminate contamination with bacterial DNA that occurred during the purification process, a purification procedure was developed to produce DNA-free protein. In addition, to remove the hexahistidine tag from the recombinant protein, thrombin cleavage was carried out while the recombinant protein was bound to a DNA column. This procedure stabilized the recombinant protein and resulted in nearly 100% cleavage with no appreciable loss to unwanted proteolytic degradation. This nondenaturing purification scheme results in high-quality, native PARP-DBD for use in structural and biochemical studies.

Class-directed structure determination: foundation for a protein structure initiative.

The recent sequencing of many complete genomes, combined with the development of methods that allow rapid structure determination for many proteins, has changed the way in which protein structure determinations can be approached. One-by-one determinations of individual protein structures will soon be augmented by class-directed structure analyses in which a group of proteins is targeted and structures of representative members are determined and used to represent the entire group. Such a shift in approach would be the foundation for a broad protein structure initiative targeting classes of proteins important for biotechnology and for a fundamental understanding of protein function.