SPINE: an integrated tracking database and data mining approach for identifying feasible targets in high-throughput structural proteomics.

High-throughput structural proteomics is expected to generate considerable amounts of data on the progress of structure determination for many proteins. For each protein this includes information about cloning, expression, purification, biophysical characterization and structure determination via NMR spectroscopy or X-ray crystallography. It will be essential to develop specifications and ontologies for standardizing this information to make it amenable to retrospective analysis. To this end we created the SPINE database and analysis system for the Northeast Structural Genomics Consortium. SPINE, which is available at bioinfo.mbb.yale.edu/nesg or nesg.org, is specifically designed to enable distributed scientific collaboration via the Internet. It was designed not just as an information repository but as an active vehicle to standardize proteomics data in a form that would enable systematic data mining. The system features an intuitive user interface for interactive retrieval and modification of expression construct data, query forms designed to track global project progress and external links to many other resources. Currently the database contains experimental data on 985 constructs, of which 740 are drawn from Methanobacterium thermoautotrophicum, 123 from Saccharomyces cerevisiae, 93 from Caenorhabditis elegans and the remainder from other organisms. We developed a comprehensive set of data mining features for each protein, including several related to experimental progress (e.g. expression level, solubility and crystallization) and 42 based on the underlying protein sequence (e.g. amino acid composition, secondary structure and occurrence of low complexity regions). We demonstrate in detail the application of a particular machine learning approach, decision trees, to the tasks of predicting a protein's solubility and propensity to crystallize based on sequence features. We are able to extract a number of key rules from our trees, in particular that soluble proteins tend to have significantly more acidic residues and fewer hydrophobic stretches than insoluble ones. One of the characteristics of proteomics data sets, currently and in the foreseeable future, is their intermediate size ( approximately 500-5000 data points). This creates a number of issues in relation to error estimation. Initially we estimate the overall error in our trees based on standard cross-validation. However, this leaves out a significant fraction of the data in model construction and does not give error estimates on individual rules. Therefore, we present alternative methods to estimate the error in particular rules.

Structure of Thermotoga maritima stationary phase survival protein SurE: a novel acid phosphatase.

BACKGROUND: The rpoS, nlpD, pcm, and surE genes are among many whose expression is induced during the stationary phase of bacterial growth. rpoS codes for the stationary-phase RNA polymerase sigma subunit, and nlpD codes for a lipoprotein. The pcm gene product repairs damaged proteins by converting the atypical isoaspartyl residues back to L-aspartyls. The physiological and biochemical functions of surE are unknown, but its importance in stress is supported by the duplication of the surE gene in E. coli subjected to high-temperature growth. The pcm and surE genes are highly conserved in bacteria, archaea, and plants. RESULTS: The structure of SurE from Thermotoga maritima was determined at 2.0 A. The SurE monomer is composed of two domains; a conserved N-terminal domain, a Rossman fold, and a C-terminal oligomerization domain, a new fold. Monomers form a dimer that assembles into a tetramer. Biochemical analysis suggests that SurE is an acid phosphatase, with an optimum pH of 5.5-6.2. The active site was identified in the N-terminal domain through analysis of conserved residues. Structure-based site-directed point mutations abolished phosphatase activity. T. maritima SurE intra- and intersubunit salt bridges were identified that may explain the SurE thermostability. CONCLUSIONS: The structure of SurE provided information about the protein's fold, oligomeric state, and active site. The protein possessed magnesium-dependent acid phosphatase activity, but the physiologically relevant substrate(s) remains to be identified. The importance of three of the assigned active site residues in catalysis was confirmed by site-directed mutagenesis.

Chemical Biology Approaches for Characterization of Epigenetic Regulators.

Chemical biology approaches are a powerful means to functionally characterize epigenetic regulators such as histone modifying enzymes. We outline experimental protocols and best practices for the cellular characterization and use of "chemical probes" that selectively inhibit protein methyltransferases, many of which methylate histones to regulate heritable gene expression patterns. We describe biomarker assays to validate the probes in specific cellular systems, and provide guidelines for their use in functional characterization of methyltransferases including detailed protocols, examples, and controls. Together these techniques enable precision manipulation of cellular epigenomes and the exploration of the therapeutic potential of epigenetic targets in human disease.

New insights into p53 function from structural studies.

Recent structural analysis of p53 has greatly enhanced our understanding of the biochemical activities of this protein by presenting us with a detailed picture of the chemical groups in the protein that are involved in protein stability, conformation and functional interactions. The current structures form the basis for the design of potential therapeutics which could, for example, revert a DNA-binding mutant back to a DNA-binding competent conformation. The structure of the tet domain forms the basis for designing an active therapeutic p53 with an oligomerization domain which would not cross react with a DNA-binding mutant p53. However, as useful as these structures have been in providing insight into the structure/function relationship for p53, a complete understanding of this protein awaits more detailed information on the full-length protein. In this respect, one of the most useful roles for future structural studies will be to help identify the nature of the conformational transition between latent and active p53, and how it can be modulated.

Characterization of the oligomerization defects of two p53 mutants found in families with Li-Fraumeni and Li-Fraumeni-like syndrome.

Recently two germline mutations in the oligomerization domain of p53 have been identified in patients with Li-Fraumeni and Li-Fraumeni-like Syndromes. We have used biophysical and biochemical methods to characterize these two mutants in order to better understand their functional defects and the role of the p53 oligomerization domain (residues 325-355) in oncogenesis. We find that residues 310-360 of the L344P mutant are monomeric, apparently unfolded and cannot interact with wild-type (WT) p53. The full length L344P protein is unable to bind sequence specifically to DNA and is therefore an inactive, but not a dominant negative mutant. R337C, on the other hand, can form dimers and tetramers, can hetero-oligomerize with WTp53 and can bind to a p53 consensus element. However, the thermal stability of R337C is much lower than that of WTp53 and at physiological temperatures more than half of this mutant is less than tetrameric. Thus, the R337C mutant retains some functional activity yet leads to a predisposition to cancer, suggesting that even partial inactivation of p53 oligomerization is sufficient for accelerated tumour progression.

Structure and function in the p53 family.

The recent discovery of several p53 homologs has uncovered a p53 superfamily of transcription factors that can trigger cell cycle arrest and apoptosis. The challenge now is to understand the similarities and differences between family members especially in terms of their regulation and potential for physical or genetic interactions with one another. This review summarizes recent progress in understanding the structure-function relationship within the p53 family. The new family members, p63 and p73, have an additional conserved domain at their C-termini which may have a regulatory function. The structure of this domain (a SAM domain) suggests that it is a protein-protein interaction module that may be involved in developmental processes. The oligomerization domains of p53 family members, while conserved in sequence and three-dimensional structure do not interact appreciably with other family members, but do mediate interactions between the multiple splice variants from an individual gene.

Refined solution structures of the Escherichia coli trp holo- and aporepressor.

The solution structures of the trp-repressor from Escherichia coli in both the liganded (holo-) and unliganded (apo-) form, have been refined by restrained molecular dynamics with simulated annealing using the program XPLOR and additional experimental constraints. The ensemble of refined holorepressor structures have a root-mean-square deviation (r.m.s.d.) of 0.8 A relative to the average structure for the backbone of the dimer core (helices A, B, C, A', B', C') and 2.5 A for the helix-turn-helix DNA-binding domain (helices D and E). The corresponding values for the aporepressor are 0.9 A for the backbone of the ABC-dimer core and 3.2 A for the DE helix-turn-helix. The r.m.s.d. of the average structures from the corresponding crystal structures are 2.3 A for the holorepressor ABC core and 4.2 A for its DE region; 2.3 A for the aporepressor core and 5.5 A for its DE region. The relative disorder of the DNA-binding domain is reflected in a number of experimental parameters including substantially more rapid backbone proton exchange rates, exchange-limited relaxation times and crystallographic B-factors. The stabilizing effect of the L-Trp ligand is evident in these measurements, as it is in the higher precision of the holorepressor structure.

Structure and functionality of a designed p53 dimer.

P53 is a homotetrameric tumor suppressor protein involved in transcriptional control of genes that regulate cell proliferation and death. In order to probe the role that oligomerization plays in this capacity, we have previously designed and characterized a series of p53 proteins with altered oligomeric states through hydrophilc substitution of residues Met340 or Leu344 in the normally tetrameric oligomerization domain. Although such mutations have little effect on the overall secondary structural content of the oligomerization domain, both solubility and the resistance to thermal denaturation are substantially reduced relative to that of the wild-type domain. Here, we report the design and characterization of a double-mutant p53 with alterations of residues at positions Met340 and Leu344. The double-mutations Met340Glu/Leu344Lys and Met340Gln/Leu344Arg resulted in distinct dimeric forms of the protein. Furthermore, we have verified by NMR structure determination that the double-mutant Met340Gln/Leu344Arg is essentially a "half-tetramer". Analysis of the in vivo activities of full-length p53 oligomeric mutants reveals that while cell-cycle arrest requires tetrameric p53, transcriptional transactivation activity of monomers and dimers retain roughly background and half of the wild-type activity, respectively.

The solution structure of bacteriophage lambda protein W, a small morphogenetic protein possessing a novel fold.

Protein W (gpW) from bacteriophage lambda is required for the stabilization of DNA within the phage head and for attachment of tails onto the head during morphogenesis. Although comprised of only 68 residues, it likely interacts with at least two other proteins in the mature phage and with DNA. Thus, gpW is an intriguing subject for detailed structural studies. We have determined its solution structure using NMR spectroscopy and have found it to possesses a novel fold consisting of two alpha-helices and a single two-stranded beta-sheet arranged around a well-packed hydrophobic core. The 14 C-terminal residues of gpW, which are essential for function, are unstructured in solution.

Cooperative interaction between the DNA-binding domains of PU.1 and IRF4.

The two lymphoid-specific transcription factors PU.1 and IRF4 form a cooperative ternary complex at immunoglobulin enhancer elements such as the lambdaB and kappaE3' sites. We report here that the synergy of this interaction can be reconstituted in part with the DNA-binding domains of the two proteins. The minimal DNA binding-domain of IRF4 was mapped to residues 20 to 137, corresponding to the conserved DNA-binding region of other interferon regulatory factors (IRFs). This domain can bind weakly to a synthetic murine lambdaB element, while IRF4 constructs that contain residues 1 to 19 require the presence of PU.1 for DNA-binding at similar concentrations. Fluorescence polarization of fluorescein-labelled DNA was used to show that the presence of residues 1 to 19 decreases the binding affinity of IRF4 N-terminal constructs from two- to fivefold. However, all constructs bound better to the lambdaB element in the presence of the DNA-binding domain of PU.1. This cooperative interaction was not dependent on phosphorylation of the PEST domain of PU.1, but was dependent on the proper spacing of the binding sites for PU.1 and IRF4. These data suggest that at least part of the cooperative interaction between full-length PU.1 and IRF4 involves the DNA-binding domains of the two proteins. NMR spectroscopy of 15N-labelled PU.1 and IRF4 constructs indicates that the PEST domain of PU.1 and residues 1 to 19 of IRF4 may be unstructured in the isolated proteins.

Structure-based functional classification of hypothetical protein MTH538 from Methanobacterium thermoautotrophicum.

The structure of MTH538, a previously uncharacterized hypothetical protein from Methanobacterium thermoautotrophicum, has been determined by NMR spectroscopy. MTH538 is one of numerous structural genomics targets selected in a genome-wide survey of uncharacterized sequences from this organism. MTH538 is a so-called singleton, a sequence not closely related to any other (known) sequences. The structure of MTH538 closely resembles the known structures of receiver domains from two component response regulator systems, such as CheY, and is similar to the structures of flavodoxins and GTP-binding proteins. Tests on MTH538 for characteristic activities of CheY and flavodoxin were negative. MTH538 did not become phosphorylated in the presence of acetyl phosphate and Mg(2+), although it appeared to bind Mg(2+). MTH538 also did not bind flavin mononucleotide (FMN) or coenzyme F(420). Nevertheless, sequence and structure parallels between MTH538/CheY and two families of ATPase/phosphatase proteins suggest that MTH538 may have a role in a phosphorylation-independent two-component response regulator system.

Solution structure and dynamics of yeast elongin C in complex with a von Hippel-Lindau peptide.

Elongin is a transcription elongation factor that stimulates the rate of elongation by suppressing transient pausing by RNA polymerase II at many sites along the DNA. It is heterotrimeric in mammals, consisting of elongins A, B and C subunits, and bears overall similarity to a class of E3 ubiquitin ligases known as SCF (Skp1-Cdc53 (cullin)-F-box) complexes. A subcomplex of elongins B and C is a target for negative regulation by the von Hippel-Lindau (VHL) tumor-suppressor protein. Elongin C from Saccharomyces cerevisiae, Elc1, exhibits high sequence similarity to mammalian elongin C. Using NMR spectroscopy we have determined the three-dimensional structure of Elc1 in complex with a human VHL peptide, VHL(157-171), representing the major Elc1 binding site. The bound VHL peptide is entirely helical. Elc1 utilizes two C-terminal helices and an intervening loop to form a binding groove that fits VHL(157-171). Chemical shift perturbation and dynamics analyses reveal that a global conformational change accompanies Elc1/VHL(157-171) complex formation. Moreover, the disappearance of conformational exchange phenomena on the microsecond to millisecond time scale within Elc1 upon VHL peptide binding suggests a role for slow internal motions in ligand recognition.

Study of a noncovalent trp repressor: DNA operator complex by electrospray ionization time-of-flight mass spectrometry.

Electrospray ionization time-of-flight mass spectrometry (ESI-TOF MS) has been used to study noncovalent interactions between the trp apo-repressor (TrpR), its co-repressor tryptophan and its specific operator DNA. In 5 mM ammonium acetate, TrpR was detected as a partially unfolded monomer. In the presence of a 21-base-pair DNA possessing the two symmetrically arranged CTAG consensus sequences required for specific TrpR binding, a homodimer-dsDNA complex with a 1:1 stoichiometry was observed. Co-repressor was not needed for the complex to form under our experimental conditions. Collision induced dissociation (CID-MS) revealed that this complex was very stable in the gas phase since dissociation was achieved only at energies that also broke covalent bonds. We saw no evidence for the presence of the six water molecules that mediate the interaction between the protein and the DNA in the crystal structure. To check the binding specificity of the TrpR for its target DNA, a competitive experiment was undertaken: the protein was mixed with an equimolar amount of three different DNAs in which the two CTAG sequences were separated by 2, 4, and 6 bp, respectively. Only the DNA with the correct consensus spacing of 4 bp was able to form stable interactions with TrpR. This experiment demonstrates the potential of ESI-MS to test the sequence-specificity of protein-DNA complexes. The interactions between the TrpR-DNA complex and 5-methyl-, L- and D-tryptophan were also investigated. Two molecules of 5-methyl- or L-tryptophan were bound with high affinity to the TrpR-DNA complex. On the other hand, D-tryptophan appeared to bind to the complex with poor specificity and poor affinity.

Structural analysis and comparison of the C-terminal transport signal domains of hemolysin A and leukotoxin A.

NMR spectroscopy was used to study the structure of the C-terminal signal sequences of the bacterial toxins, hemolysin A(HlyA) and leukotoxin A (LktA). The two signals share little sequence homology; however, both can direct toxin transport with equal efficiency. We report here that in a membrane mimetic environment both peptides form two short non-interacting alpha-helices separated by a short loop. This higher order structure may be a common feature of C-terminal signals and may be required for interaction with the membrane associated transporter complex.

Optimization of resolution and sensitivity of 4D NOESY using multi-dimensional decomposition.

Highly resolved multi-dimensional NOE data are essential for rapid and accurate determination of spatial protein structures such as in structural genomics projects. Four-dimensional spectra contain almost no spectral overlap inherently present in lower dimensionality spectra and are highly amenable to application of automated routines for spectral resonance location and assignment. However, a high resolution 4D data set using conventional uniform sampling usually requires unacceptably long measurement time. Recently we have reported that the use of non-uniform sampling and multi-dimensional decomposition (MDD) can remedy this problem. Here we validate accuracy and robustness of the method, and demonstrate its usefulness for fully protonated protein samples. The method was applied to 11 kDa protein PA1123 from structural genomics pipeline. A systematic evaluation of spectral reconstructions obtained using 15-100% subsets of the complete reference 4D 1H-13C-13C-1H NOESY spectrum has been performed. With the experimental time saving of up to six times, the resolution and the sensitivity per unit time is shown to be similar to that of the fully recorded spectrum. For the 30% data subset we demonstrate that the intensities in the reconstructed and reference 4D spectra correspond with a correlation coefficient of 0.997 in the full range of spectral amplitudes. Intensities of the strong, middle and weak cross-peaks correlate with coefficients 0.9997, 0.9965, and 0.83. The method does not produce false peaks. 2% of weak peaks lost in the 30% reconstruction is in line with theoretically expected noise increase for the shorter measurement time. Together with good accuracy in the relative line-widths these translate to reliable distance constrains derived from sparsely sampled, high resolution 4D NOESY data.

DNA binding specificity studies of four ETS proteins support an indirect read-out mechanism of protein-DNA recognition.

Members of the ETS family of transcription factors are involved in several developmental and physiological processes, and, when overexpressed or misexpressed, can contribute to a variety of cancers. Each family member has a conserved DNA-binding domain that recognizes DNA sequences containing a G-G-A trinucleotide. Discrimination between potential ETS-binding sites appears to be governed by both the nucleotides flanking the G-G-A sequence and protein-protein interactions. We have used an adaptation of the "length-encoded multiplex" approach (Desjarlais, J. R., and Berg, J. M. (1994) Proc. Natl. Acad. Sci. U. S. A. 91, 11099-11103) to define DNA binding specificities for four ETS proteins: Fli-1, SAP-1, PU.1, and TEL. Our results support a model in which cooperative effects among neighboring bases flanking the central G-G-A site contribute to the formation of stable ETS/DNA complexes. These results are consistent with a mechanism for specific DNA binding that is partially governed by an indirect read-out of the DNA sequence, in which a sequence-specific DNA conformation is sensed or induced.

Secretion and circular dichroism analysis of the C-terminal signal peptides of HlyA and LktA.

The secretion of the 107 kDa hemolysin A (HlyA) from Escherichia coli is mediated by membrane proteins hemolysin B (HlyB) and hemolysin D (HlyD). The signal for transport has been mapped to the C-terminal 60 amino acids of the HylA molecule. We have shown previously that the C-terminal 70 amino acids of leukotoxin (LktA) from Pasteurella hemolytica can substitute functionally for the HlyA signal sequence. This 70 amino acid peptide contains little primary sequence similarity to the HlyA signal sequence, and we have hypothesized that these signal sequences assume a similar higher-order structure which is recognized by the HlyB/D transporter. In the present study, we have expressed and purified small peptides containing the C-terminal 61 amino acids of HlyA and the C-terminal 70 amino acids of LktA. We show that these signal peptides are sufficient for secretion from E. coli in a HlyB/D dependent manner. Circular dichroism analyses show that both molecules exhibit common biophysical properties. In aqueous solution, they appear to be mainly unstructured, but in a membrane mimetic environment they assume a helical secondary structure. The conformational change observed for both peptides going from an aqueous to a membrane mimetic environment may be an important feature of these signal sequences necessary for their recognition and transport.

Solution structure of a conserved C-terminal domain of p73 with structural homology to the SAM domain.

p73 and p63 are two recently cloned genes with homology to the tumor suppressor p53, whose protein product is a key transcriptional regulator of genes involved in cell cycle arrest and apoptosis. While all three proteins share conserved transcriptional activation, DNA-binding and oligomerization domains, p73 and p63 have an additional conserved C-terminal region. We have determined the three-dimensional solution structure of this conserved C-terminal domain of human p73. The structure reveals a small five-helix bundle with striking similarity to the SAM (sterile alpha motif) domains of two ephrin receptor tyrosine kinases. The SAM domain is a putative protein-protein interaction domain found in a variety of cytoplasmic signaling proteins and has been shown to form both homo- and hetero-oligomers. However, the SAM-like C-terminal domains of p73 and p63 are monomeric and do not interact with one another, suggesting that this domain may interact with additional, as yet uncharacterized proteins in a signaling and/or regulatory role.

New perceptions of transcription factor properties from NMR.

The complementarity of NMR and X-ray crystallography for biomacromolecular studies has been particularly evident in analysis of transcription factor structures and interactions. While X-ray crystallography can be used to tackle relatively complicated structural problems including multicomponent (three and higher) complexes, NMR studies have provided new insights into the nature of protein-DNA and protein-protein interactions that would be difficult to obtain by other biophysical methods. We describe herein some of the novel and important information recently derived from NMR studies of transcription factors.