Structural genomics: opportunities and challenges

A Strategy for High-Throughput Assay Development Using Leads Derived from Nuclear Magnetic Resonance-Based Screening

A strategy is described for the development of high-throughput screening assays against targets of unknown function that involves the use of nuclear magnetic resonance (NMR) spectroscopy. Using this approach, molecules that bind to the protein target are identified from an NMR-based screen of a library of substrates, cofactors, and other compounds that are known to bind to many proteins and enzymes. Once a ligand has been discovered, a fluorescent or radiolabeled analog of the ligand is synthesized that can be used in a high-throughput screen. The approach is illustrated in the development of a high-throughput screening assay against HI-0033, a conserved protein from Haemophilus influenzae whose function is currently unknown. Adenosine was found to bind to HI-0033 by NMR, and fluorescent analogs were rapidly identified that bound to HI-0033 in the submicromolar range. Using these fluorescent compounds, a fluorescence polarization assay was developed that is suitable for high-throughput screening and obtaining detailed structure-activity relationships for lead optimization.

Utilisation of adsorption and desorption for simultaneously improving protein crystallisation success rate and crystal quality

High-quality protein crystals of suitable size are an important prerequisite for applying X-ray crystallography to determine the 3-dimensional structure of proteins. However, it is often difficult to obtain protein crystals of appropriate size and quality because nucleation and growth processes can be unsuccessful. Here, we show that by adsorbing proteins onto porous polystyrene-divinylbenzene microspheres (SDB) floating on the surface of the crystallisation solution, a localised high supersaturation region at the surface of the microspheres and a low supersaturation region below the microspheres can coexist in a single solution. The crystals will easily nucleate in the region of high supersaturation, but when they grow to a certain size, they will sediment to the region of low supersaturation and continue to grow. In this way, the probability of crystallisation and crystal quality can be simultaneously increased in a single solution without changing other crystallisation parameters.

Applications of the second virial coefficient: protein crystallization and solubility

This article begins by highlighting some of the ground-based studies emanating from NASA's Microgravity Protein Crystal Growth (PCG) program. This is followed by a more detailed discussion of the history of and the progress made in one of the NASA-funded PCG investigations involving the use of measured second virial coefficients (B values) as a diagnostic indicator of solution conditions conducive to protein crystallization. A second application of measured B values involves the determination of solution conditions that improve or maximize the solubility of aqueous and membrane proteins. These two important applications have led to several technological improvements that simplify the experimental expertise required, enable the measurement of membrane proteins and improve the diagnostic capability and measurement throughput.

What an rRNA secondary structure tells about phylogeny of fungi in Ascomycota with emphasis on evolution of major types of ascus

Background

RNA secondary structure is highly conserved throughout evolution. The higher order structure is fundamental in establishing important structure-function relationships. Nucleotide sequences from ribosomal RNA (rRNA) genes have made a great contribution to our understanding of Ascomycota phylogeny. However, filling the gaps between molecular phylogeny and morphological assumptions based on ascus dehiscence modes and type of fruitbodies at the higher level classification of the phylum remains an unfulfilled task faced by mycologists.

Methodology/Principal Findings

We selected some major groups of Ascomycota to view their phylogenetic relationships based on analyses of rRNA secondary structure. Using rRNA secondary structural information, here, we converted nucleotide sequences into the structure ones over a 20-symbol code. Our structural analyses together with ancestral character state reconstruction produced reasonable phylogenetic position for the class Geoglossomycetes as opposed to the classic nucleotide analyses. Judging from the secondary structure analyses with consideration of mode of ascus dehiscence and the ability of forming fruitbodies, we draw a clear picture of a possible evolutionary route for fungal asci and some major groups of fungi in Ascomycota. The secondary structure trees show a more reasonable phylogenetic position for the class Geoglossomycetes.

Conclusions

Our results illustrate that asci lacking of any dehiscence mechanism represent the most primitive type. Passing through the operculate and Orbilia-type asci, bitunicate asci occurred. The evolution came to the most advanced inoperculate type. The ascus-producing fungi might be derived from groups lacking of the capacity to form fruitbodies, and then evolved multiple times. The apothecial type of fruitbodies represents the ancestral state, and the ostiolar type is advanced. The class Geoglossomycetes is closely related to Leotiomycetes and Sordariomycetes having a similar ascus type other than it was originally placed based on nucleotide sequence analyses.

Cell-free protein production system with the E. coli crude extract for determination of protein folds

Escherichia coli cell extract-based coupled transcription-translation cell-free system has been developed for large-scale production of protein samples for both X-ray crystallography (selenomethionine substitution) and NMR (stable-isotope labeling). For both cases, higher labeling/substitution efficiency can be achieved compared with the production using cell-based expression system. In addition, as the system is easily adapted to automated and/or high-throughput procedures, it is an especially suitable protein expression method for structural genomics and proteomics project. In this chapter, the procedure for large-scale protein production for structure determination using our E. coli cell-free system is presented.

Fishing with (Proto)Net-a principled approach to protein target selection

Structural genomics strives to represent the entire protein space. The first step towards achieving this goal is by rationally selecting proteins whose structures have not been determined, but that represent an as yet unknown structural superfamily or fold. Once such a structure is solved, it can be used as a template for modelling homologous proteins. This will aid in unveiling the structural diversity of the protein space. Currently, no reliable method for accurate 3D structural prediction is available when a sequence or a structure homologue is not available. Here we present a systematic methodology for selecting target proteins whose structure is likely to adopt a new, as yet unknown superfamily or fold. Our method takes advantage of a global classification of the sequence space as presented by ProtoNet-3D, which is a hierarchical agglomerative clustering of the proteins of interest (the proteins in Swiss-Prot) along with all solved structures (taken from the PDB). By navigating in the scaffold of ProtoNet-3D, we yield a prioritized list of proteins that are not yet structurally solved, along with the probability of each of the proteins belonging to a new superfamily or fold. The sorted list has been self-validated against real structural data that was not available when the predictions were made. The practical application of using our computational–statistical method to determine novel superfamilies for structural genomics projects is also discussed.

Identification of genes encoding tRNA modification enzymes by comparative genomics

As the molecular adapters between codons and amino acids, transfer-RNAs are pivotal molecules of the genetic code. The coding properties of a tRNA molecule do not reside only in its primary sequence. Posttranscriptional nucleoside modifications, particularly in the anticodon loop, can modify cognate codon recognition, affect aminoacylation properties, or stabilize the codon-anticodon wobble base pairing to prevent ribosomal frameshifting. Despite a wealth of biophysical and structural knowledge of the tRNA modifications themselves, their pathways of biosynthesis had been until recently only partially characterized. This discrepancy was mainly due to the lack of obvious phenotypes for tRNA modification-deficient strains and to the difficulty of the biochemical assays used to detect tRNA modifications. However, the availability of hundreds of whole-genome sequences has allowed the identification of many of these missing tRNA-modification genes. This chapter reviews the methods that were used to identify these genes with a special emphasis on the comparative genomic approaches. Methods that link gene and function but do not rely on sequence homology will be detailed, with examples taken from the tRNA modification field.

Improving expression and solubility of rice proteins produced as fusion proteins in Escherichia coli

For proteins of higher eukaryotes, such as plants, which have large genomes, recombinant protein expression and purification are often difficult. Expression levels tend to be low and the expressed proteins tend to misfold and aggregate. We tested seven different expression vectors in Escherichia coli for rapid subcloning of rice genes and for protein expression and solubility levels. Each expressed gene product has an N-terminal fusion protein and/or tag, and an engineered protease site upstream of the mature rice protein. Several different fusion proteins/tags and protease sites were tested. We found that the fusion proteins and the protease sites have significant and varying effects on expression and solubility levels. The expression vector with the most favorable characteristics is pDEST-trx. The vector, which is a modified version of the commercially available expression vector, pET-32a, contains an N-terminal thioredoxin fusion protein and a hexahistidine tag, and is adapted to the Gateway expression system. However, addition of an engineered protease site could drastically change the expression and solubility properties. We selected 135 genes corresponding to potentially interesting rice proteins, transferred the genes from cDNAs to expression vectors, and engineered in suitable protease sites N-terminal to the mature proteins. Of 135 genes, 131 (97.0%) could be expressed and 72 (53.3%) were soluble when the fusion proteins/tags were present. Thirty-eight mature-length rice proteins and domains (28.1%) are suitable for NMR solution structure studies and/or X-ray crystallography. Our expression systems are useful for the production of soluble plant proteins in E. coli to be used for structural genomics studies.

A limited universe of membrane protein families and folds

One of the goals of structural genomics is to obtain a structural representative of almost every fold in nature. A recent estimate suggests that 70%-80% of soluble protein domains identified in the first 1000 genome sequences should be covered by about 25,000 structures-a reasonably achievable goal. As no current estimates exist for the number of membrane protein families, however, it is not possible to know whether family coverage is a realistic goal for membrane proteins. Here we find that virtually all polytopic helical membrane protein families are present in the already known sequences so we can make an estimate of the total number of families. We find that only approximately 700 polytopic membrane protein families account for 80% of structured residues and approximately 1700 cover 90% of structured residues. While apparently a finite and reachable goal, we estimate that it will likely take more than three decades to obtain the structures needed for 90% residue coverage, if current trends continue.

Target selection and deselection at the Berkeley Structural Genomics Center

At the Berkeley Structural Genomics Center (BSGC), our goal is to obtain a near-complete structural complement of proteins in the minimal organisms Mycoplasma genitalium and M. pneumoniae, two closely related pathogens. Current targets for structure determination have been selected in six major stages, starting with those predicted to be most tractable to high throughput study and likely to yield new structural information. We report on the process used to select these proteins, as well as our target deselection procedure. Target deselection reduces experimental effort by eliminating targets similar to those recently solved by the structural biology community or other centers. We measure the impact of the 69 structures solved at the BSGC as of July 2004 on structure prediction coverage of the M. pneumoniae and M. genitalium proteomes. The number of Mycoplasma proteins for which the fold could first be reliably assigned based on structures solved at the BSGC (24 M. pneumoniae and 21 M. genitalium) is approximately 25% of the total resulting from work at all structural genomics centers and the worldwide structural biology community (94 M. pneumoniae and 86 M. genitalium) during the same period. As the number of structures contributed by the BSGC during that period is less than 1% of the total worldwide output, the benefits of a focused target selection strategy are apparent. If the structures of all current targets were solved, the percentage of M. pneumoniae proteins for which folds could be reliably assigned would increase from approximately 57% (391 of 687) at present to around 80% (550 of 687), and the percentage of the proteome that could be accurately modeled would increase from around 37% (254 of 687) to about 64% (438 of 687). In M. genitalium, the percentage of the proteome that could be structurally annotated based on structures of our remaining targets would rise from 72% (348 of 486) to around 76% (371 of 486), with the percentage of accurately modeled proteins would rise from 50% (243 of 486) to 58% (283 of 486). Sequences and data on experimental progress on our targets are available in the public databases TargetDB and PEPCdb.

Improvements to the throughput of recombinant protein expression in the baculovirus/insect cell system

Recombinant baculoviruses have proved to be a very useful means to express many proteins over the last 20 years. Since their introduction, there have been a number of significant improvements that have simplified and speeded up the construction of baculoviruses. One of the most commonly used methods relies upon recombination with the baculovirus genome maintained in Escherichia coli. In this paper, we report the conversion of nearly all the steps in this process including the expression testing and purification to a multi-well plate format. This enables a significant increase in the number of constructs that can be processed in a shorter period of time and an order of magnitude increase in the number of expression conditions that can be analysed. A key step in our process is that the transfection is done in suspension rather than adherent cells, which gives a much higher virus titre than in the standard methods.

Kinomics: characterizing the therapeutically validated kinase space

The annotation and visualization of medicinally relevant kinase space revealed that kinase inhibitors in the clinic are, on average, of higher molecular weight and more lipophilic than all other clinically investigated drugs. Tyrosine kinases from the vascular endothelial growth factor and epidermal growth factor receptor families are the most pursued targets. Furthermore, oncological indications account for 75% of all kinase-related clinical interest. In addition, analysis of the similarity between kinase targets with respect to sequence, selectivity and structure has revealed that kinases with > or =60% sequence identity are most likely to be inhibited by the same classes of compounds and have similar ATP-binding sites. The identification of this threshold, together with the widely accepted representation of the sequence-based kinase space, is expanding our understanding of the clinical and structural space of the kinome.

Solution structure of At3g04780.1-des15, an Arabidopsis thaliana ortholog of the C-terminal domain of human thioredoxin-like protein

The structure of At3g04780.1-des15, an Arabidopsis thaliana ortholog of the C-terminal domain of human thioredoxin-like protein, was determined by NMR spectroscopy. The structure is dominated by a beta-barrel sandwich. A two-stranded anti-parallel beta-sheet, which seals off one end of the beta-barrel, is flanked by two flexible loops rich in acidic amino acids. Although this fold often provides a ligand binding site, the structure did not reveal an appreciable cavity inside the beta-barrel. The three-dimensional structure of At3g04780.1-des15 provides an entry point for understanding its functional role and those of its mammalian homologs.

From Penicillin to the Ribosome: Revolutions in the Determination and Use of Molecular Structure in Chemistry and Biology

A revolution in structural analysis is in progress in the biological sciences that parallels a similar revolution that took place in chemistry 40–50 years ago. This has major implications for chemistry, offering exciting opportunities at the interface between chemistry and biology. The advances are driven by the value of structural information in biology, for understanding biological function, and for applications in structure-based drug design and structural genomics. Two directions are apparent: towards technically challenging biological structures and assemblies, typified by the potassium channel and the ribosome; and towards high-throughput structure determination of many, smaller, proteins, as in structural genomics. In this review, the advances in molecular biology and in structure determination by X-ray crystallography that make these developments possible are discussed, together with appropriate examples.

High-throughput screening of structural proteomics targets using NMR

We applied a high-throughput strategy for the screening of targets for structural proteomics of Xanthomonas axonopodis pv citri. This strategy is based on the rapid (1)H-(15)N HSQC NMR analysis of bacterial lysates containing selectively (15)N-labelled heterologous proteins. Our analysis permitted us to classify the 19 soluble candidates in terms of 'foldedness', that is, the extent to which they present a well-folded solution structure, as reflected by the quality of their NMR spectra. This classification allowed us to define a priority list to be used as a guide to select protein candidates for further structural studies.

Role for NMR in structural genomics

The 2nd EMSL Workshop on Structural Genomics was held on 28th and 29th July 2000 at the Environmental Molecular Sciences Laboratory at the Department of Energy's Pacific Northwest National Laboratory in Richland, WA. The workshop focused on four topics: 1. The role for NMR in structural and functional genomics; 2. The technical challenges NMR faces for structural and functional genomics; 3. The potential need for a national NMR center for structural and functional genomics in the United States; and 4. Organization of the NMR community. This report summarizes the workshop proceedings and conclusions reached regarding the role of NMR in the emerging fields of structural and functional genomics.

Specific radiation damage can be used to solve macromolecular crystal structures

The use of third generation synchrotron sources has led to renewed concern about the effect of ionizing radiation on crystalline biological samples. In general, the problem is seen as one to be avoided. However, in this paper, it is shown that, far from being a hindrance to successful structure determination, radiation damage provides an opportunity for phasing macromolecular structures. This is successfully demonstrated for both a protein and an oligonucleotide, by way of which complete models were built automatically. The possibility that, through the exploitation of radiation damage, the phase problem could become less of a barrier to macromolecular crystal structure determination is discussed.

Protein expression systems for structural genomics and proteomics

One of the key steps of structural genomics and proteomics is high-throughput expression of many target proteins. Gene cloning, especially by ligation-independent cloning techniques, and recombinant protein expression using microbial hosts such as Escherichia coli and the yeast Pichia pastoris are well optimized and further robotized. Cell-free protein synthesis systems have been developed for large-scale production of protein samples for NMR (stable-isotope labeling) and X-ray crystallography (selenomethionine substitution). Protein folding is still a major bottleneck in protein expression. Cell-based and cell-free methods for screening of suitable samples for structure determination have been developed for achieving a high success rate.

Structuring the universe of proteins

High-throughput sequencing of human genomes and those of important model organisms (mouse, Drosophila melanogaster, Caenorhabditis elegans, fungi, archaea) and bacterial pathogens has laid the foundation for another "big science" initiative in biology. Together, X-ray crystallographers, nuclear magnetic resonance (NMR) spectroscopists, and computational biologists are pursuing high-throughput structural studies aimed at developing a comprehensive three-dimensional view of the protein structure universe. The new science of structural genomics promises more than 10,000 experimental protein structures and millions of calculated homology models of related proteins. The evolutionary underpinnings and technological challenges of automating target selection, protein expression and purification, sample preparation, NMR and X-ray data measurement/analysis, homology modeling, and structure/function annotation are discussed in detail. An informative case study from one of the structural genomics centers funded by the National Institutes of Health and the National Institute of General Medical Sciences (NIH/NIGMS) demonstrates how this experimental/computational pipeline will reveal important links between form and function in biology and provide new insights into evolution and human health and disease.

MAMMOTH (matching molecular models obtained from theory): an automated method for model comparison

Advances in structural genomics and protein structure prediction require the design of automatic, fast, objective, and well benchmarked methods capable of comparing and assessing the similarity of low-resolution three-dimensional structures, via experimental or theoretical approaches. Here, a new method for sequence-independent structural alignment is presented that allows comparison of an experimental protein structure with an arbitrary low-resolution protein tertiary model. The heuristic algorithm is given and then used to show that it can describe random structural alignments of proteins with different folds with good accuracy by an extreme value distribution. From this observation, a structural similarity score between two proteins or two different conformations of the same protein is derived from the likelihood of obtaining a given structural alignment by chance. The performance of the derived score is then compared with well established, consensus manual-based scores and data sets. We found that the new approach correlates better than other tools with the gold standard provided by a human evaluator. Timings indicate that the algorithm is fast enough for routine use with large databases of protein models. Overall, our results indicate that the new program (MAMMOTH) will be a good tool for protein structure comparisons in structural genomics applications. MAMMOTH is available from our web site at http://physbio.mssm.edu/~ortizg/.

Tracing the evolution of RNA structure in ribosomes

The elucidation of ribosomal structure has shown that the function of ribosomes is fundamentally confined to dynamic interactions established between the RNA components of the ribosomal ensemble. These findings now enable a detailed analysis of the evolution of ribosomal RNA (rRNA) structure. The origin and diversification of rRNA was studied here using phylogenetic tools directly at the structural level. A rooted universal tree was reconstructed from the combined secondary structures of large (LSU) and small (SSU) subunit rRNA using cladistic methods and considerations in statistical mechanics. The evolution of the complete repertoire of structural ribosomal characters was formally traced lineage-by-lineage in the tree, showing a tendency towards molecular simplification and a homogeneous reduction of ribosomal structural change with time. Character tracing revealed patterns of evolution in inter-subunit bridge contacts and tRNA-binding sites that were consistent with the proposed coupling of tRNA translocation and subunit movement. These patterns support the concerted evolution of tRNA-binding sites in the two subunits and the ancestral nature and common origin of certain structural ribosomal features, such as the peptidyl (P) site, the functional relay of the penultimate stem helix of SSU rRNA, and other structures participating in ribosomal dynamics. Overall results provide a rare insight into the evolution of ribosomal structure.

Alternative splicing: combinatorial output from the genome

Alternative splicing has emerged as a mechanism that can account for a large proportion of the disparity between the modest number of genes in the human genome and the much higher complexity of the expressed proteome. At least a third, and probably the majority, of human genes are alternatively spliced, and some genes can generate thousands of protein isoforms by complex alternative splicing events. Analysis of the transcriptome will therefore require the development of massively parallel technologies that are able to encompass the complexity arising from alternative splicing.

Deriving folds of macromolecular complexes through electron cryomicroscopy and bioinformatics approaches

Intermediate-resolution (7-9A) structures of large macromolecular complexes can be obtained by electron cryomicroscopy. This structural information, combined with bioinformatics data for the individual protein components or domains, can lead to a fold model for the entire complex. Such approaches have been demonstrated with the 6.8 A structure of the rice dwarf virus to derive models for the major capsid shell proteins.

High-throughput crystallization and structure determination in drug discovery

High-throughput protein X-ray crystallography offers an unprecedented opportunity to facilitate drug discovery. The key bottlenecks in the path from target gene to three-dimensional protein structure determination are defined. Special emphasis is placed on the concept that drug discovery projects are typically directed at a key protein target whose structure must be solved within a reasonable time frame to have an impact on the drug discovery process. The time-sensitive nature of structural data has placed growing pressure on the need to automate all aspects of protein crystallography, from gene identification to model building and refinement. Current technological innovations and strategies for automation are discussed with respect to the bottleneck they are intended to eliminate.

A tour of structural genomics

Structural genomics projects aim to provide an experimental or computational three-dimensional model structure for all of the tractable macromolecules that are encoded by complete genomes. To this end, pilot centres worldwide are now exploring the feasibility of large-scale structure determination. Their experimental structures and computational models are expected to yield insight into the molecular function and mechanism of thousands of proteins. The pervasiveness of this information is likely to change the use of structure in molecular biology and biochemistry.