Evolution and Structural Theory: The Frontier Between Chemistry and Biology. In: Dugas H, editor. Bioorganic Chemistry Frontiers. Berlin: Springer Berlin Heidelberg; 1990. pp. 1-70. [DOI: 10.1007/978-3-642-75256-8_1]

Artificially Expanded Genetic Information Systems for New Aptamer Technologies

Directed evolution was first applied to diverse libraries of DNA and RNA molecules a quarter century ago in the hope of gaining technology that would allow the creation of receptors, ligands, and catalysts on demand. Despite isolated successes, the outputs of this technology have been somewhat disappointing, perhaps because the four building blocks of standard DNA and RNA have too little functionality to have versatile binding properties, and offer too little information density to fold unambiguously. This review covers the recent literature that seeks to create an improved platform to support laboratory Darwinism, one based on an artificially expanded genetic information system (AEGIS) that adds independently replicating nucleotide “letters” to the evolving “alphabet”.

Molecular paleoscience: systems biology from the past

Experimental paleomolecular biology, paleobiochemistry, and paleogenetics are closely related emerging fields that infer the sequences of ancient genes and proteins from now-extinct organisms, and then resurrect them for study in the laboratory. The goal of paleogenetics is to use information from natural history to solve the conundrum of modern genomics: How can we understand deeply the function of biomolecular structures uncovered and described by modern chemical biology? Reviewed here are the first 20 cases where biomolecular resurrections have been achieved. These show how paleogenetics can lead to an understanding of the function of biomolecules, analyze changing function, and put meaning to genomic sequences, all in ways that are not possible with traditional molecular biological studies.

Functional inferences from reconstructed evolutionary biology involving rectified databases--an evolutionarily grounded approach to functional genomics

If bioinformatics tools are constructed to reproduce the natural, evolutionary history of the biosphere, they offer powerful approaches to some of the most difficult tasks in genomics, including the organization and retrieval of sequence data, the updating of massive genomic databases, the detection of database error, the assignment of introns, the prediction of protein conformation from protein sequences, the detection of distant homologs, the assignment of function to open reading frames, the identification of biochemical pathways from genomic data, and the construction of a comprehensive model correlating the history of biomolecules with the history of planet Earth.

The B12-dependent ribonucleotide reductase from the archaebacterium Thermoplasma acidophila: an evolutionary solution to the ribonucleotide reductase conundrum

A coenzyme B12-dependent ribonucleotide reductase was purified from the archaebacterium Thermoplasma acidophila and partially sequenced. Using probes derived from the sequence, the corresponding gene was cloned, completely sequenced, and expressed in Escherichia coli. The deduced amino acid sequence shows that the catalytic domain of the B12-dependent enzyme from T. acidophila, some 400 amino acids, is related by common ancestry to the diferric tyrosine radical iron(III)-dependent ribonucleotide reductase from E. coli, yeast, mammalian viruses, and man. The critical cysteine residues in the catalytic domain that participate in the thiyl radical-dependent reaction have been conserved even though the cofactor that generates the radical is not. Evolutionary bridges created by the T. acidophila sequence and that of a B12-dependent reductase from Mycobacterium tuberculosis establish homology between the Fe-dependent enzymes and the catalytic domain of the Lactobacillus leichmannii B12-dependent enzyme as well. These bridges are confirmed by a predicted secondary structure for the Lactobacillus enzyme. Sequence similarities show that the N-terminal domain of the T. acidophila ribonucleotide reductase is also homologous to the anaerobic ribonucleotide reductase from E. coli, which uses neither B12 nor Fe cofactors. A predicted secondary structure of the N-terminal domain suggests that it is predominantly helical, as is the domain in the aerobic E. coli enzyme depending on Fe, extending the homologous family of proteins to include anaerobic ribonucleotide reductases, B12 ribonucleotide reductases, and Fe-dependent aerobic ribonucleotide reductases. A model for the evolution of the ribonucleotide reductase family is presented; in this model, the thiyl radical-based reaction mechanism is conserved, but the cofactor is chosen to best adapt the host organism to its environment. This analysis illustrates how secondary structure predictions can assist evolutionary analyses, each important in "post-genomic" biochemistry.

Pseudogenes in ribonuclease evolution: a source of new biomacromolecular function?

Bovine seminal ribonuclease (RNase) diverged from pancreatic RNase after a gene duplication ca. 35 million years ago. Members of the seminal RNase gene family evidently remained as unexpressed pseudogene for much of its evolutionary history. Between 5 and 10 million years ago, however, after the divergence of kudu but before the divergence of ox, evidence suggests that the pseudogene was repaired and expressed. Intriguingly, detailed analysis of the sequences suggests that the repair may have involved gene conversion, transfer of information from the pancreatic gene to the RNase pseudogene. Further, the ratio of non-silent to silent substitutions suggests that the pancreatic RNases are divergently evolving under functional constraints, the seminal RNase pseudogenes are diverging under no functional constraints, while the genes expressed in the seminal plasma are evolving extremely rapidly in their amino acid sequences, as if to fulfil a new physiological role.

A predicted consensus structure for the protein kinase C2 homology (C2H) domain, the repeating unit of synaptotagmin

A secondary structure has been predicted for the protein kinase C2 regulatory domain found in homologous form in synaptotagmin, some phospholipases, and some GTP activated proteins. The proposed structure is built from seven consecutive beta strands followed by a terminal alpha helix. Considerations of overall surface exposure of individual secondary structural elements suggest that these are packed into a 2-sheet beta sandwich structure, with one of only three of the many possible folds being preferred.

A consensus prediction of the secondary structure for the 6-phospho-beta-D-galactosidase superfamily

Two separate unrefined models for the secondary structure of two subfamilies of the 6-phospho-beta-D-galactosidase superfamily were independently constructed by examining patterns of variation and conservation within homologous protein sequences, assigning surface, interior, parsing, and active site residues to positions in the alignment, and identifying periodicities in these. A consensus model for the secondary structure of the entire superfamily was then built. The prediction tests the limits of an unrefined prediction made using this approach in a large protein with substantial functional and sequence divergence within the family. The protein belongs to the (alpha-beta class), with the core beta strands aligned parallel. The supersecondary structural elements that are readily identified in this model is a parallel beta sheet built by strands C, D, and E, with helices 2 and 3 connecting strands (C+D) and (D+E), respectively, and an analogous beta-alpha unit (strand G and helix 7) toward the end of the sequence. The resemblance of the supersecondary model to the tertiary structure formed by 8-fold alpha-beta barrel proteins is almost certainly not coincidental.

Predicting the conformation of proteins from sequences. Progress and future progress

Recent progress in structure prediction has allowed bona fide predictions, those made and published before an experimental structure is determined, to be remarkably accurate. The most successful methods rely on an analysis of patterns of conservation and variation within homologous protein sequences, extract tertiary structural information before secondary structure is predicted, and avoid 'three state per residue scores' as a tool for evaluating a prediction, focusing instead on efforts to understand why a prediction is successful when it is successful, and why it fails when it fails.

A secondary structure model of the integrin alpha subunit N-terminal domain based on analysis of multiple alignments

The integrins are alpha/beta heterodimeric proteins which mediate cell-matrix and cell-cell interactions. Current data indicate that the N-terminal moiety of the alpha subunit is involved in ligand binding. This region of the receptor is made up of a seven-fold repeated sequence of unknown structure which contains EF-hand-like putative divalent cation-binding sites. Recent studies have shown that multiple sequence alignments can be analysed to yield secondary structure predictions. Therefore, to obtain a model structure for the integrin alpha subunit N-terminal domain repeat, a large alignment of the seven repeats from sixteen integrin sequences was generated. Two methods of analysis were used: First, Chou and Fasman and Garnier, Osguthorpe and Robson predictions were carried out for individual sequences and the consensus predictions derived. Consensus hydrophobicity and chain flexibility data were also used to provide additional data. Second, sites of conservation and variation were analysed by a computer program STAMA (STructure After Multiple Alignment) to yield a secondary structure prediction. The two analyses gave essentially the same predicted structure: undefined region, loop, alpha-helix, beta-strand, divalent cation-binding loop, beta-strand, putative turn, loop, beta-strand. This is the first model structure to be presented for an integrin domain. Its implications for integrin function are discussed.

Predicting the conformation of proteins from sequences. Progress and future progress

A new paradigm for predicting the secondary and tertiary structure of functional proteins from sequence data has emerged from detailed models of how natural selection, conservation, and neutral drift, the three fundamental factors in molecular evolution, leave their mark upon protein sequences. Structural information is extracted from a set of aligned homologous sequences via an analysis of patterns of conservation and variation between proteins with quantitatively defined evolutionary relationships. Tertiary structural information is obtained prior to the assignment of secondary structure, where it plays an important role. Throughout, structural predictions are made with the active involvement of a biochemist whose expertise and insight is critical both for making the prediction and in analyzing its successful and unsuccessful parts. Secondary structure predictions are evaluated based on their ability to sustain an effort to model tertiary structure. Several predictions made using the new paradigm can now be compared with those made under the classical paradigm, including a neural network. The results obtained from the new paradigm are clearly superior to those obtained with the classical paradigm, at least within the protein families that were examined.

Predicting the conformation of proteins. Man versus machine

Two types of approaches for predicting the conformation of proteins from sequence data have lately received attention: 'black box' tools that generate fully automated predictions of secondary structure from a set of homologous protein sequences, and methods involving the expertise of a human biochemist who is assisted, but not replaced, by computer tools. A friendly controversy has emerged as to which approach offers a brighter future. In fact, both are necessary. Nevertheless, a snapshot of the controversy at this instant offers much insight into the structure prediction problem itself.

The nitrogenase MoFe protein. A secondary structure prediction

Surface residues, interior residues, and parsing residues, together with a secondary structure derived from these, are predicted for the MoFe nitrogenase protein in advance of a crystal structure of the protein, scheduled shortly to appear in Nature. By publishing this prediction, we test our method for predicting the conformation of proteins from patterns in the divergent evolution of homologous protein sequences in a way that places the method 'at risk'.

Patterns of divergence in homologous proteins as indicators of secondary and tertiary structure: a prediction of the structure of the catalytic domain of protein kinases

The secondary structure and elements of tertiary structure have been predicted for the catalytic domain of protein kinases using a method that extracts structural information from the patterns of conservation and variation in an alignment of homologous proteins. The central features of this structural prediction are: (a) the catalytic domains of protein kinases do not incorporate a Rossmann fold; (b) the core of the structure is founded on beta sheets built from pairs of bent antiparallel beta strands; (c) five helices, including an especially long helix (alignment positions 129-152) that lie on the outside of the folded core. These proteins are important in many aspects of metabolic regulation.