Structure-based assignment of the biochemical function of a hypothetical protein: A test case of structural genomics

Many small bacterial, archaebacterial, and eukaryotic genomes have been sequenced, and the larger eukaryotic genomes are predicted to be completely sequenced within the next decade. In all genomes sequenced to date, a large portion of these organisms’ predicted protein coding regions encode polypeptides of unknown biochemical, biophysical, and/or cellular functions. Three-dimensional structures of these proteins may suggest biochemical or biophysical functions. Here we report the crystal structure of one such protein, MJ0577, from a hyperthermophile, Methanococcus jannaschii, at 1.7-Å resolution. The structure contains a bound ATP, suggesting MJ0577 is an ATPase or an ATP-mediated molecular switch, which we confirm by biochemical experiments. Furthermore, the structure reveals different ATP binding motifs that are shared among many homologous hypothetical proteins in this family. This result indicates that structure-based assignment of molecular function is a viable approach for the large-scale biochemical assignment of proteins and for discovering new motifs, a basic premise of structural genomics.

Structure-based conformational preferences of amino acids

Proteins can be very tolerant to amino acid substitution, even within their core. Understanding the factors responsible for this behavior is of critical importance for protein engineering and design. Mutations in proteins have been quantified in terms of the changes in stability they induce. For example, guest residues in specific secondary structures have been used as probes of conformational preferences of amino acids, yielding propensity scales. Predicting these amino acid propensities would be a good test of any new potential energy functions used to mimic protein stability. We have recently developed a protein design procedure that optimizes whole sequences for a given target conformation based on the knowledge of the template backbone and on a semiempirical potential energy function. This energy function is purely physical, including steric interactions based on a Lennard-Jones potential, electrostatics based on a Coulomb potential, and hydrophobicity in the form of an environment free energy based on accessible surface area and interatomic contact areas. Sequences designed by this procedure for 10 different proteins were analyzed to extract conformational preferences for amino acids. The resulting structure-based propensity scales show significant agreements with experimental propensity scale values, both for α-helices and β-sheets. These results indicate that amino acid conformational preferences are a natural consequence of the potential energy we use. This confirms the accuracy of our potential and indicates that such preferences should not be added as a design criterion.

Recent improvements in prediction of protein structure by global optimization of a potential energy function

Recent improvements of a hierarchical ab initio or de novo approach for predicting both α and β structures of proteins are described. The united-residue energy function used in this procedure includes multibody interactions from a cumulant expansion of the free energy of polypeptide chains, with their relative weights determined by Z-score optimization. The critical initial stage of the hierarchical procedure involves a search of conformational space by the conformational space annealing (CSA) method, followed by optimization of an all-atom model. The procedure was assessed in a recent blind test of protein structure prediction (CASP4). The resulting lowest-energy structures of the target proteins (ranging in size from 70 to 244 residues) agreed with the experimental structures in many respects. The entire experimental structure of a cyclic α-helical protein of 70 residues was predicted to within 4.3 Å α-carbon (Cα) rms deviation (rmsd) whereas, for other α-helical proteins, fragments of roughly 60 residues were predicted to within 6.0 Å Cα rmsd. Whereas β structures can now be predicted with the new procedure, the success rate for α/β- and β-proteins is lower than that for α-proteins at present. For the β portions of α/β structures, the Cα rmsd's are less than 6.0 Å for contiguous fragments of 30–40 residues; for one target, three fragments (of length 10, 23, and 28 residues, respectively) formed a compact part of the tertiary structure with a Cα rmsd less than 6.0 Å. Overall, these results constitute an important step toward the ab initio prediction of protein structure solely from the amino acid sequence.

A test of the "jigsaw puzzle" model for protein folding by multiple methionine substitutions within the core of T4 lysozyme.

To test whether the structure of a protein is determined in a manner akin to the assembly of a jigsaw puzzle, up to 10 adjacent residues within the core of T4 lysozyme were replaced by methionine. Such variants are active and fold cooperatively with progressively reduced stability. The structure of a seven-methionine variant has been shown, crystallographically, to be similar to wild type and to maintain a well ordered core. The interaction between the core residues is, therefore, not strictly comparable with the precise spatial complementarity of the pieces of a jigsaw puzzle. Rather, a certain amount of give and take in forming the core structure is permitted. A simplified hydrophobic core sequence, imposed without genetic selection or computer-based design, is sufficient to retain native properties in a globular protein.

Chromatin structure and gene expression.

It is now well understood that chromatin structure is perturbed in the neighborhood of expressed genes. This is most obvious in the neighborhood of promoters and enhancers, where hypersensitivity to nucleases marks sites that no longer carry canonical nucleosomes, and to which transcription factors bind. To study the relationship between transcription factor binding and the generation of these hypersensitive regions, we mutated individual cis-acting regulatory elements within the enhancer that lies between the chicken beta- and epsilon-globin genes. Constructions carrying the mutant enhancer were introduced by stable transformation into an avian erythroid cell line. We observed that weakening the enhancer resulted in creation of two classes of site: those still completely accessible to nuclease attack and those that were completely blocked. This all-or-none behavior suggests a mechanism by which chromatin structure can act to sharpen the response of developmental systems to changing concentrations of regulatory factors. Another problem raised by chromatin structure concerns the establishment of boundaries between active and inactive chromatin domains. We have identified a DNA element at the 5' end of the chicken beta-globin locus, near such a boundary, that has the properties of an insulator; in test constructions, it blocks the action of an enhancer on a promoter when it is placed between them. We describe the properties and partial dissection of this sequence. A third problem is posed by the continued presence of nucleosomes on transcribed genes, which might prevent the passage of RNA polymerase. We show, however, that a prokaryotic polymerase can transcribe through a histone octamer on a simple chromatin template. The analysis of this process reveals that an octamer is capable of transferring from a position in front of the polymerase to one behind, without ever losing its attachment to the DNA.

Behavior predicts genes structure in a wild primate group.

The predictability of genetic structure from social structure and differential mating success was tested in wild baboons. Baboon populations are subdivided into cohesive social groups that include multiple adults of both sexes. As in many mammals, males are the dispersing sex. Social structure and behavior successfully predicted molecular genetic measures of relatedness and variance in reproductive success. In the first quantitative test of the priority-of-access model among wild primates, the reproductive priority of dominant males was confirmed by molecular genetic analysis. However, the resultant high short-term variance in reproductive success did not translate into equally high long-term variance because male dominance status was unstable. An important consequence of high but unstable short-term variance is that age cohorts will tend to be paternal sibships and social groups will be genetically substructured by age.

Antisense globin RNA in mouse erythroid tissues: structure, origin, and possible function.

The aim of the experiments described in this paper was to test for the presence of antisense globin RNA in mouse erythroid tissues and, if found, to characterize these molecules. The present study made use of a multistep procedure in which a molecular tag is attached to cellular RNA by ligation with a defined ribooligonucleotide. The act of ligation preserves the termini of RNA molecules, which become the junctions between cellular RNAs and the ligated ribooligonucleotide. It also unambiguously preserves the identity of cellular RNA as a sense or antisense molecule through all subsequent manipulations. Using this approach, we identified and characterized antisense beta-globin RNA in erythroid spleen cells and reticulocytes from anemic mice. We show in this paper that the antisense globin RNA is fully complementary to spliced globin mRNA, indicative of the template/transcript relationship. It terminates at the 5' end with a uridylate stretch, reflecting the presence of poly(A) at the 3' end of the sense globin mRNA. With respect to the structure of their 3' termini, antisense globin RNA can be divided into three categories: full-size molecules corresponding precisely to globin mRNA, truncated molecules lacking predominantly 14 3'-terminal nucleotides, and extended antisense RNA containing 17 additional 3'-terminal nucleotides. The full-size antisense globin RNA contains two 14-nt-long complementary sequences within its 3'-terminal segment corresponding to the 5'-untranslated region of globin mRNA. This, together with the nature of the predominant truncation, suggests a mechanism by which antisense RNA might give rise to new sense-strand globin mRNA.

Maintenance of pre-mRNA secondary structure by epistatic selection.

Linkage disequilibrium between polymorphisms in a natural population may result from various evolutionary forces, including random genetic drift due to sampling of gametes during reproduction, restricted migration between subpopulations in a subdivided population, or epistatic selection. In this report, we present evidence that the majority of significant linkage disequilibria observed in introns of the alcohol dehydrogenase locus (Adh) of Drosophila pseudoobscura are due to epistatic selection maintaining secondary structure of precursor mRNA (pre-mRNA). Based on phylogenetic-comparative analysis and a likelihood approach, we propose secondary structure models of Adh pre-mRNA for the regions of the adult intron and intron 2 where clustering of linkage disequilibria has been observed. Furthermore, we applied the likelihood ratio test to the phylogenetically predicted secondary structure in intron 1. In contrast to the other two structures, polymorphisms associated with the more conserved stem-loop structure of intron 1 are in low frequency, and linkage disequilibria have not been observed. These findings are qualitatively consistent with a model of compensatory fitness interactions. This model assumes that mutations disrupting pairing in a secondary structural element are individually deleterious if they destabilize a functionally important structure; a second "compensatory" mutation, however, may restabilize the structure and restore fitness.

Repacking protein cores with backbone freedom: structure prediction for coiled coils.

Progress in homology modeling and protein design has generated considerable interest in methods for predicting side-chain packing in the hydrophobic cores of proteins. Present techniques are not practically useful, however, because they are unable to model protein main-chain flexibility. Parameterization of backbone motions may represent a general and efficient method to incorporate backbone relaxation into such fixed main-chain models. To test this notion, we introduce a method for treating explicitly the backbone motions of alpha-helical bundles based on an algebraic parameterization proposed by Francis Crick in 1953 [Crick, F. H. C. (1953) Acta Crystallogr. 6, 685-689]. Given only the core amino acid sequence, a simple calculation can rapidly reproduce the crystallographic main-chain and core side-chain structures of three coiled coils (one dimer, one trimer, and one tetramer) to within 0.6-A root-mean-square deviations. The speed of the predictive method [approximately 3 min per rotamer choice on a Silicon Graphics (Mountain View, CA) 4D/35 computer] permits it to be used as a design tool.