Structure-based functional inference in structural genomics

The dramatically increasing number of new protein sequences arising from genomics 4 proteomics requires the need for methods to rapidly and reliably infer the molecular and cellular functions of these proteins. One such approach, structural genomics, aims to delineate the total repertoire of protein folds in nature, thereby providing three-dimensional folding patterns for all proteins and to infer molecular functions of the proteins based on the combined information of structures and sequences. The goal of obtaining protein structures on a genomic scale has motivated the development of high throughput technologies and protocols for macromolecular structure determination that have begun to produce structures at a greater rate than previously possible. These new structures have revealed many unexpected functional inferences and evolutionary relationships that were hidden at the sequence level. Here, we present samples of structures determined at Berkeley Structural Genomics Center and collaborators laboratories to illustrate how structural information provides and complements sequence information to deduce the functional inferences of proteins with unknown molecular functions.Two of the major premises of structural genomics are to discover a complete repertoire of protein folds in nature and to find molecular functions of the proteins whose functions are not predicted from sequence comparison alone. To achieve these objectives on a genomic scale, new methods, protocols, and technologies need to be developed by multi-institutional collaborations worldwide. As part of this effort, the Protein Structure Initiative has been launched in the United States (PSI; www.nigms.nih.gov/funding/psi.html). Although infrastructure building and technology development are still the main focus of structural genomics programs [1–6], a considerable number of protein structures have already been produced, some of them coming directly out of semi-automated structure determination pipelines [6–10]. The Berkeley Structural Genomics Center (BSGC) has focused on the proteins of Mycoplasma or their homologues from other organisms as its structural genomics targets because of the minimal genome size of the Mycoplasmas as well as their relevance to human and animal pathogenicity (http://www.strgen.org). Here we present several protein examples encompassing a spectrum of functional inferences obtainable from their three-dimensional structures in five situations, where the inferences are new and testable, and are not predictable from protein sequence information alone.     

Structural Genomics of Minimal Organisms and Protein Fold Space

The initial aim of the Berkeley Structural Genomics Center is to obtain a near-complete structural complement of two minimal organisms, closely related pathogens Mycoplasma genitalium and M. pneumoniae. The former has fewer than 500 genes and the latter fewer than 700 genes. To achieve this goal, the current protein targets have been selected starting with those predicted to be most tractable and likely to yield new structural and functional information. During the past 3 years, the semi-automated structural genomics pipeline has been set up from cloning, expression, purification, and ultimately to structural determination. The results from the pipeline substantially increased the coverage of the protein fold space of M. pneumoniae and M. genitalium. Furthermore, about 1/2 of the structures of ‘unique’ protein sequences revealed new and novel folds, and over 2/3 of the structures of previously annotated ‘hypothetical proteins’ inferred their molecular functions.     

Structures of proteins of biomedical interest from the Center for Eukaryotic Structural Genomics

The Center for Eukaryotic Structural Genomics (CESG) produces and solves the structures of proteins from eukaryotes. We have developed and operate a pipeline to both solve structures and to test new methodologies. Both NMR and X-ray crystallography methods are used for structure solution. CESG chooses targets based on sequence dissimilarity to known structures, medical relevance, and nominations from members of the scientific community. Many times proteins qualify in more than one of these categories. Here we review some of the structures that have connections to human health and disease.     

Structural genomics of highly conserved microbial genes of unknown function in search of new antibacterial targets

With more than 100 antibacterial drugs at our disposal in the 1980s, the problem of bacterial infection was considered solved. Today, however, most hospital infections are insensitive to several classes of antibacterial drugs, and deadly strains of Staphylococcus aureus resistant to vancomycin – the last resort antibiotic – have recently begin to appear. Other life-threatening microbes, such as Enterococcus faecalis and Mycobacterium tuberculosis are already able to resist every available antibiotic. There is thus an urgent, and continuous need for new, preferably large-spectrum, antibacterial molecules, ideally targeting new biochemical pathways. Here we report on the progress of our structural genomics program aiming at the discovery of new antibacterial gene targets among evolutionary conserved genes of uncharacterized function. A series of bioinformatic and comparative genomics analyses were used to identify a set of 221 candidate genes common to Gram-positive and Gram-negative bacteria. These genes were split between two laboratories. They are now submitted to a systematic 3-D structure determination protocol including cloning, protein expression and purification, crystallization, X-ray diffraction, structure interpretation, and function prediction. We describe here our strategies for the 111 genes processed in our laboratory. Bioinformatics is used at most stages of the production process and out of 111 genes processed – and 17 months into the project – 108 have been successfully cloned, 103 have exhibited detectable expression, 84 have led to the production of soluble protein, 46 have been purified, 12 have led to usable crystals, and 7 structures have been determined.     

Life in the fast lane for protein crystallization and X-ray crystallography

The common goal for structural genomic centers and consortiums is to decipher as quickly as possible the three-dimensional structures for a multitude of recombinant proteins derived from known genomic sequences. Since X-ray crystallography is the foremost method to acquire atomic resolution for macromolecules, the limiting step is obtaining protein crystals that can be useful of structure determination. High-throughput methods have been developed in recent years to clone, express, purify, crystallize and determine the three-dimensional structure of a protein gene product rapidly using automated devices, commercialized kits and consolidated protocols. However, the average number of protein structures obtained for most structural genomic groups has been very low compared to the total number of proteins purified. As more entire genomic sequences are obtained for different organisms from the three kingdoms of life, only the proteins that can be crystallized and whose structures can be obtained easily are studied. Consequently, an astonishing number of genomic proteins remain unexamined. In the era of high-throughput processes, traditional methods in molecular biology, protein chemistry and crystallization are eclipsed by automation and pipeline practices. The necessity for high-rate production of protein crystals and structures has prevented the usage of more intellectual strategies and creative approaches in experimental executions. Fundamental principles and personal experiences in protein chemistry and crystallization are minimally exploited only to obtain “low-hanging fruit” protein structures. We review the practical aspects of today's high-throughput manipulations and discuss the challenges in fast pace protein crystallization and tools for crystallography. Structural genomic pipelines can be improved with information gained from low-throughput tactics that may help us reach the higher-bearing fruits. Examples of recent developments in this area are reported from the efforts of the Southeast Collaboratory for Structural Genomics (SECSG).     

Crystal structures of a phosphotransacetylase from <Emphasis Type="Italic">Bacillus subtilis</Emphasis> and its complex with acetyl phosphate

Phosphotransacetylase (Pta) [EC 2.3.1.8] plays a major role in acetate metabolism by catalyzing the reversible transfer of the acetyl group between coenzyme A (CoA) and orthophosphate: CH(3)COSCoA+HPO(4)(2-)<-->CH(3)COOPO(3)(2-) +CoASH. In this study, we report the crystal structures of Pta from Bacillus subtilis at 2.75 A resolution and its complex with acetyl phosphate, one of its substrates, at 2.85 A resolution. In addition, the Pta activity of the enzyme has been assayed. The enzyme folds into an alpha/beta architecture with two domains separated by a prominent cleft, very similar to two other known Pta structures. The enzyme-acetyl phosphate complex structure reveals a few potential substrate binding sites. Two of them are located in the middle of the interdomain cleft: each one is surrounded by a region of strictly and highly conserved residues. High structural similarities are found with 4-hydroxythreonine-4-phosphate dehydrogenase (PdxA), and isocitrate and isopropylmalate dehydrogenases, all of which utilize NADP+ as their cofactor, which binds in the interdomain cleft. Their substrate binding sites are close to the acetyl phosphate binding sites of Pta in the cleft as well. These results suggest that the CoA is likely to bind to the interdomain cleft of Pta in a similar way as NADP+ binds to the other three enzymes.     

2007 Annual progress report synopsis of the Center for Structures of Membrane Proteins

A synopsis of the 2007 annual progress report for the Center for Structures of Membrane Proteins, a specialized center of the Protein Structure Initiative.     

Structural reorganization of molecular sheets derived from cellulose II by molecular dynamics simulations

We investigated structural reorganization of two different kinds of molecular sheets derived from the cellulose II crystal using molecular dynamics (MD) simulations, in order to identify the initial structure of the cellulose crystal in the course of its regeneration process from solution. After a one-nanosecond simulation, the molecular sheet formed by van der Waals forces along the () crystal plane did not change its structure in an aqueous environment, while the other one formed by hydrogen bonds along the (1 1 0) crystal plane changed into a van der Waals-associated molecular sheet, such as the former. The two structures that were calculated showed substantial similarities such as the high occupancy of intramolecular hydrogen bonds between O3^H and O5 of over 0.75, few intermolecular hydrogen bonds, and the high occurrence of hydrogen bonding with water. The convergence of the two structures into one denotes that the van der Waals-associated molecular sheet can be the initial structure of the cellulose crystal formed in solution. The main chain conformations were almost the same as those in the cellulose II crystal except for a −16° shift of φ (dihedral angle of O5-C1-O1-C4) and the gauche-gauche conformation of the hydroxymethyl side group appears probably due to its hydrogen bonding with water. These results suggest that the van der Waals-associated molecular sheet becomes stable in an aqueous environment with its hydrophobic inside and hydrophilic periphery. Contrary to this, a benzene environment preferred a hydrogen-bonded molecular sheet, which is expected to be the initial structure formed in benzene.     

Structural and mechanistic insight from high resolution structures of archaeal rhodopsins

Lipidic cubic phase-grown crystals yielded high resolution structures of a number of archaeal retinal proteins, the molecular mechanisms of which are being revealed as structures of photocycle intermediates become available. The structural basis for bacteriorhodopsin's mechanism of proton pumping is discussed, revealing a well-synchronized sequence of molecular events. Comparison with the high resolution structures of the halide pump halorhodopsin, as well as with the receptor sensory rhodopsin II, illustrates how small and localized structural changes result in functional divergence. Fundamental principles of energy transduction and sensory reception in the archaeal rhodopsins, which may have relevance to other systems, are discussed.     

Structural analysis of a novel class of R-M controller proteins: C.Csp231I from Citrobacter sp. RFL231

Controller proteins play a key role in the temporal regulation of gene expression in bacterial restriction-modification (R-M) systems and are important mediators of horizontal gene transfer. They form the basis of a highly cooperative, concentration-dependent genetic switch involved in both activation and repression of R-M genes. Here we present biophysical, biochemical, and high-resolution structural analysis of a novel class of controller proteins, exemplified by C.Csp231I. In contrast to all previously solved C-protein structures, each protein subunit has two extra helices at the C-terminus, which play a large part in maintaining the dimer interface. The DNA binding site of the protein is also novel, having largely AAAA tracts between the palindromic recognition half-sites, suggesting tight bending of the DNA. The protein structure shows an unusual positively charged surface that could form the basis for wrapping the DNA completely around the C-protein dimer.     

Structural genomics on membrane proteins: comparison of more than 100 GPCRs in 3 expression systems

Production of recombinant receptors has been one of the major bottlenecks in structural biology on G protein-coupled receptors (GPCRs). The MePNet (Membrane Protein Network) was established to overexpress a large number of GPCRs in three major expression systems, based on Escherichia coli, Pichia pastoris and Semliki Forest virus (SFV) vectors. Evaluation by immunodetection demonstrated that 50% of a total of 103 GPCRs were expressed in bacterial inclusion bodies, 94% in yeast cell membranes and 95% in SFV-infected mammalian cells. The expression levels varied from low to high and the various GPCR families and subtypes were analyzed for their expressability in each expression system. More than 60% of the GPCRs were expressed at milligram levels or higher in one or several systems, compatible to structural biology applications. Functional activity was determined by binding assays in yeast and mammalian cells and the correlation between immunodetection and binding activity was analyzed.     

Evaluation of models for the evolution of protein sequences and functions under structural constraint

In the field of evolutionary structural genomics, methods are needed to evaluate why genomes evolved to contain the fold distributions that are observed. In order to study the effects of population dynamics in the evolved genomes we need fast and accurate evolutionary models which can analyze the effects of selection, drift and fixation of a protein sequence in a population that are grounded by physical parameters governing the folding and binding properties of the sequence. In this study, various knowledge-based, force field, and statistical methods for protein folding have been evaluated with four different folds: SH2 domains, SH3 domains, Globin-like, and Flavodoxin-like, to evaluate the speed and accuracy of the energy functions. Similarly, knowledge-based and force field methods have been used to predict ligand binding specificity in SH2 domain. To demonstrate the applicability of these methods, the dynamics of evolution of new binding capabilities by an SH2 domain is demonstrated.     

Structural diversity in the six-fold redundant set of acyl-CoA carboxyltransferases in Mycobacterium tuberculosis

Mycobacterium tuberculosis contains multiple versions of the accA and accD genes that encode the - and β-subunits of at least three distinct multi-functional acyl-CoA carboxylase complexes. Because of its proposed involvement in pathogenic M. tuberculosis survival, the high-resolution crystal structure of the β-subunit gene accD5 product has been determined and reveals a hexameric 356 kDa complex. Analysis of the active site properties of AccD5 and homology models of the other five M. tuberculosis AccD homologues reveals unexpected differences in their surface composition, providing a molecular rational key for a sorting mechanism governing correct acyl-CoA carboxylase holo complex assembly in M. tuberculosis.     

Challenges at the frontiers of structural biology

Knowledge of the three-dimensional structures of proteins is the key to unlocking the full potential of genomic information. There are two distinct directions along which cutting-edge research in structural biology is currently moving towards this goal. On the one hand, tightly focused long-term research in individual laboratories is leading to the determination of the structures of macromolecular assemblies of ever-increasing size and complexity. On the other hand, large consortia of structural biologists, inspired by the pace of genome sequencing, are developing strategies to determine new protein structures rapidly, so that it will soon be possible to predict reasonably accurate structures for most protein domains. We anticipate that a small number of complex systems, studied in depth, will provide insights across the field of biology with the aid of genome-based comparative structural analysis.     

The impact of extremophiles on structural genomics (and vice versa)

The advent of the complete genome sequences of various organisms in the mid-1990s raised the issue of how one could determine the function of hypothetical proteins. While insight might be obtained from a 3D structure, the chances of being able to predict such a structure is limited for the deduced amino acid sequence of any uncharacterized gene. A template for modeling is required, but there was only a low probability of finding a protein closely-related in sequence with an available structure. Thus, in the late 1990s, an international effort known as structural genomics (SG) was initiated, its primary goal to “fill sequence-structure space” by determining the 3D structures of representatives of all known protein families. This was to be achieved mainly by X-ray crystallography and it was estimated that at least 5,000 new structures would be required. While the proteins (genes) for SG have subsequently been derived from hundreds of different organisms, extremophiles and particularly thermophiles have been specifically targeted due to the increased stability and ease of handling of their proteins, relative to those from mesophiles. This review summarizes the significant impact that extremophiles and proteins derived from them have had on SG projects worldwide. To what extent SG has influenced the field of extremophile research is also discussed.     

Structural stability of myoglobin and glycomyoglobin: a comparative molecular dynamics simulation study

Glycoproteins are formed as the result of enzymatic glycosylation or chemical glycation in the body, and produced in vitro in industrial processes. The covalently attached carbohydrate molecule(s) confer new properties to the protein, including modified stability. In the present study, the structural stability of a glycoprotein form of myoglobin, bearing a glucose unit in the N-terminus, has been compared with its native form by the use of molecular dynamics simulation. Both structures were subjected to temperatures of 300 and 500 K in an aqueous environment for 10 ns. Changes in secondary structures and RMSD were then assessed. An overall higher stability was detected for glycomyoglobin, for which the most stable segments/residues were highlighted and compared with the native form. The simple addition of a covalently bound glucose is suggested to exert its stabilizing effect via increased contacts with surrounding water molecules, as well as a different pattern of interactions with neighbor residues.

Electronic supplementary material

The online version of this article (doi:10.1007/s10867-015-9383-2) contains supplementary material, which is available to authorized users.     

Inaugural structure from the DUF3349 superfamily of proteins, Mycobacterium tuberculosis Rv0543c

The first structure for a member of the DUF3349 (PF11829) family of proteins, Rv0543c from Mycobacterium tuberculosis, has been determined using NMR-based methods and some of its biophysical properties characterized. Rv0543c is a 100 residue, 11.3 kDa protein that both size exclusion chromatography and NMR spectroscopy show to be a monomer in solution. The structure of the protein consists of a bundle of five α-helices, α1 (M1 – Y16), α2 (P21 – C33), α3 (S37 – G52), α4 (G58 – H65) and α5 (S72 – G87), held together by a largely conserved group of hydrophobic amino acid side chains. Heteronuclear steady-state {¹H}–¹⁵N NOE, T₁, and T₂ values are similar through-out the sequence indicating that the backbones of the five helices are in a single motional regime. The thermal stability of Rv0543c, characterized by circular dichroism spectroscopy, indicates that Rv0543c irreversibly unfolds upon heating with an estimated melting temperature of 62.5 °C. While the biological function of Rv0543c is still unknown, the presence of DUF3349 proteins predominately in Mycobacterium and Rhodococcus bacterial species suggests that Rv0543 may have a biological function unique to these bacteria, and consequently, may prove to be an attractive drug target to combat tuberculosis.