Mars Extant Life: What's Next? Conference Report.

On November 5-8, 2019, the "Mars Extant Life: What's Next?" conference was convened in Carlsbad, New Mexico. The conference gathered a community of actively publishing experts in disciplines related to habitability and astrobiology. Primary conclusions are as follows: A significant subset of conference attendees concluded that there is a realistic possibility that Mars hosts indigenous microbial life. A powerful theme that permeated the conference is that the key to the search for martian extant life lies in identifying and exploring refugia ("oases"), where conditions are either permanently or episodically significantly more hospitable than average. Based on our existing knowledge of Mars, conference participants highlighted four potential martian refugium (not listed in priority order): Caves, Deep Subsurface, Ices, and Salts. The conference group did not attempt to reach a consensus prioritization of these candidate environments, but instead felt that a defensible prioritization would require a future competitive process. Within the context of these candidate environments, we identified a variety of geological search strategies that could narrow the search space. Additionally, we summarized a number of measurement techniques that could be used to detect evidence of extant life (if present). Again, it was not within the scope of the conference to prioritize these measurement techniques-that is best left for the competitive process. We specifically note that the number and sensitivity of detection methods that could be implemented if samples were returned to Earth greatly exceed the methodologies that could be used at Mars. Finally, important lessons to guide extant life search processes can be derived both from experiments carried out in terrestrial laboratories and analog field sites and from theoretical modeling.

Exhaustive matching of the entire protein sequence database.

The entire protein sequence database has been exhaustively matched. Definitive mutation matrices and models for scoring gaps were obtained from the matching and used to organize the sequence database as sets of evolutionarily connected components. The methods developed are general and can be used to manage sequence data generated by major genome sequencing projects. The alignments made possible by the exhaustive matching are the starting point for successful de novo prediction of the folded structures of proteins, for reconstructing sequences of ancient proteins and metabolisms in ancient organisms, and for obtaining new perspectives in structural biochemistry.

Total synthesis and cloning of a gene coding for the ribonuclease S protein.

A gene for ribonuclease S protein, has been chemically synthesized and cloned. The gene is designed to have 25 specific restriction endonuclease sites spaced at short intervals, permitting its structure to be rapidly modified. This flexibility facilitates tests of hypotheses relating the primary structure of the enzyme to its physical and catalytic behavior.

The return of pancreatic ribonucleases.

A decade after losing favor as an 'uninteresting' digestive enzyme, pancreatic ribonuclease has been found to be homologous to a series of extracellular proteins that may influence tumor cell growth, neurological development and biological differentiation. One surprising outcome of these discoveries has been the confirmation of the hypothesis that extracellular 'communicator RNA' is a messenger important in cell growth and differentiation. The only question is: why wasn't this recognized earlier?

Hydroxymethanesulfonate from Volcanic Sulfur Dioxide: A "Mineral" Reservoir for Formaldehyde and Other Simple Carbohydrates in Prebiotic Chemistry.

While formaldehyde (HCHO) was likely generated in Earth's prebiotic atmosphere by ultraviolet light, electrical discharge, and/or volcano-created lightning, HCHO could not have accumulated in substantial amounts in prebiotic environments, including those needed for prebiotic processes that generate nucleosidic carbohydrates. HCHO at high concentrations in alkaline solutions self-reacts in the Cannizzaro reaction to give methanol and formate, neither having prebiotic value. Here, we explore the possibility that volcanic sulfur dioxide (SO2) might have generated a reservoir for Hadean HCHO by a reversible reaction with HCHO to give hydroxymethanesulfonate (HMS). We show that salts of HMS are stable as solids at 90°C and do not react with themselves in solution, even at high (>8 M) concentrations. This makes them effective stores of HCHO, since the reverse reaction slowly delivers HCHO back into an environment where it can participate in prebiotically useful reactions. Specifically, we show that in alkaline borate solutions, HCHO derived from HMS allows formation of borate-stabilized carbohydrates as effectively as free HCHO, without losing material to Cannizzaro products. Further, we show that SO2 can perform similar roles for glycolaldehyde and glyceraldehyde, two intrinsically unstable carbohydrates that are needed by various models as precursors for RNA building blocks. Zircons from the Hadean show that the Hadean mantle likely provided volcanic SO2 at rates at least as great as the rates of atmospheric HCHO generation, making the formation of Hadean HMS essentially unavoidable. Thus, hydroxymethylsulfonate adducts of formaldehyde, glycolaldehyde, and glyceraldehyde, including the less soluble barium, strontium, and calcium salts, are likely candidates for prebiotically useful organic minerals on early Earth.

Evolution, language and analogy in functional genomics.

Almost a century ago, Wittgenstein pointed out that theory in science is intricately connected to language. This connection is not a frequent topic in the genomics literature. But a case can be made that functional genomics is today hindered by the paradoxes that Wittgenstein identified. If this is true, until these paradoxes are recognized and addressed, functional genomics will continue to be limited in its ability to extrapolate information from genomic sequences.

Predicted secondary structure for the Src homology 3 domain.

A de novo secondary structure prediction has been prepared for Src homology domain 3, in advance of any crystallographic information concerning any member of this interesting protein family. The prediction can be compared with a crystal structure that will be published in Nature on October 29, 1992. The prediction is based on analysis of a multiple alignment of homologous proteins. The patterns of variation and conservation of amino acids across the alignment allow the determination of surface and internal positions, which then allow the assignment of secondary structure. The prediction is quite different both in method and, in this case, result from predictions based on propensities (e.g. Garnier-Osgurthorpe-Robson) of particular amino acids to appear in particular types of secondary structure.

Bona fide prediction of aspects of protein conformation. Assigning interior and surface residues from patterns of variation and conservation in homologous protein sequences.

Heuristics have been developed for analyzing patterns of conservation and variation within a set of aligned homologous protein sequences for the purpose of assigning amino acids whose side-chains lie on the surface and inside the folded structure of a protein. These were used in several recent bona fide predictions of the secondary structure of proteins from sequence data, made and published before crystallographic information became available. Heuristics based on concurrent hydrophilic variation identify positions that lie on the surface. Heuristics based on concurrent hydrophobic conservation and variation identify positions lying in the interior. These heuristics are described here in detail and their performance evaluated when applied to seven protein families with known three-dimensional structures. The performance of individual heuristics is shown to depend on the nature of the multiple alignment within the protein family, and a strategy is presented for obtaining surface and interior assignments useful for predicting secondary structure.

Empirical and structural models for insertions and deletions in the divergent evolution of proteins.

The exhaustive matching of the protein sequence database makes possible a broadly based study of insertions and deletions (indels) during divergent evolution. In this study, the probability of a gap in an alignment of a pair of homologous protein sequences was found to increase with the evolutionary distance measured in PAM units (number of accepted point mutations per 100 amino acid residues). A relationship between the average number of amino acid residues between indels and evolutionary distance suggests that a unit 30 to 40 amino acid residues in length remains, on average, undisrupted by indels during divergent evolution. Further, the probability of a gap was found to be inversely proportional to gap length raised to the 1.7 power. This empirical law fits closely over the entire range of gap lengths examined. Gap length distribution is largely independent of evolutionary distance. These results rule out the widely used linear gap penalty as a satisfactory formula for scoring gaps when constructing alignments. Further, the observed gap length distribution can be explained by a simple model of selective pressures governing the acceptance of indels during divergent evolution. Finally, this model provides theoretical support for using indels as part of "parsing algorithms", important in the de novo prediction of the folded structure of proteins from the sequence data.

Detecting compensatory covariation signals in protein evolution using reconstructed ancestral sequences.

When protein sequences divergently evolve under functional constraints, some individual amino acid replacements that reverse the charge (e.g. Lys to Asp) may be compensated by a replacement at a second position that reverses the charge in the opposite direction (e.g. Glu to Arg). When these side-chains are near in space (proximal), such double replacements might be driven by natural selection, if either is selectively disadvantageous, but both together restore fully the ability of the protein to contribute to fitness (are together "neutral"). Accordingly, many have sought to identify pairs of positions in a protein sequence that suffer compensatory replacements, often as a way to identify positions near in space in the folded structure. A "charge compensatory signal" might manifest itself in two ways. First, proximal charge compensatory replacements may occur more frequently than predicted from the product of the probabilities of individual positions suffering charge reversing replacements independently. Conversely, charge compensatory pairs of changes may be observed to occur more frequently in proximal pairs of sites than in the average pair. Normally, charge compensatory covariation is detected by comparing the sequences of extant proteins at the "leaves" of phylogenetic trees. We show here that the charge compensatory signal is more evident when it is sought by examining individual branches in the tree between reconstructed ancestral sequences at nodes in the tree. Here, we find that the signal is especially strong when the positions pairs are in a single secondary structural unit (e.g. alpha helix or beta strand) that brings the side-chains suffering charge compensatory covariation near in space, and may be useful in secondary structure prediction. Also, "node-node" and "node-leaf" compensatory covariation may be useful to identify the better of two equally parsimonious trees, in a way that is independent of the mathematical formalism used to construct the tree itself. Further, compensatory covariation may provide a signal that indicates whether an episode of sequence evolution contains more or less divergence in functional behavior. Compensatory covariation analysis on reconstructed evolutionary trees may become a valuable tool to analyze genome sequences, and use these analyses to extract biomedically useful information from proteome databases.

Expanding the genetic lexicon: incorporating non-standard amino acids into proteins by ribosome-based synthesis.

Only 20 amino acids are normally incorporated into proteins synthesized in living cells, and this has limited the structural range of proteins that can be prepared. New methods that allow the incorporation of amino acids that are not normally encoded by natural genes are being developed: these include reassigning functions within the existing genetic code, and expanding the genetic code by constructing additional, non-natural codons. Used in conjunction with recent major advances in understanding protein structure-function relationships, these approaches should extend the range of de novo protein designs that are possible.

Extracellular 'communicator RNA'.

Evidence is presented that RNA molecules exist that act as extracellular messages important in the development of higher organisms.

An improved system for expressing pancreatic ribonuclease in Escherichia coli.

An improved method for expressing and purifying bovine pancreatic ribonuclease from a synthetic gene using the lambda promoter controlled by a temperature-sensitive repressor is described. The procedure involves isolation in the presence of a refolding buffer containing oxidized and reduced glutathione, under conditions where RNase can refold, but where proteases presumably do not. Yields are approx. 2 mg purified protein per 1 ferment.

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.

Interferon-gamma activates the cleavage of double-stranded RNA by bovine seminal ribonuclease.

Bovine seminal ribonuclease (BS-RNase), a dimeric homologue of RNase A, cleaves both single- and double-stranded RNA and inhibits the growth of tumor cells. Its catalytic activity against double-stranded RNA, either homopolymeric ([3H]polyA/polyU) or mixed sequence, is enhanced by bovine or human recombinant interferon-gamma (IFN-gamma). Activation is seen with as little as 4-10 interferon units per assay. Enhancing the degradation of double-stranded RNA, an intermediate in the growth cycle of many viruses, could contribute to IFN-gamma's ability to control cell growth and induce an antiviral state.

Site-directed mutagenesis of bovine pancreatic ribonuclease: lysine-41 and aspartate-121.

Chemical modification studies suggest that two residues of bovine pancreatic ribonuclease A (RNase A), Lys-41 and Asp-121, are important for catalysis. Three mutants of RNase A have been prepared, two point mutants with Lys-41 altered to Arg-41 and Asp-121 altered to Glu-121, and a double mutant where both residues are altered. The Lys-41 Arg mutant has ca. 2% the catalytic activity (kcat/Km) of the native protein, while the Asp-121Glu mutant has ca. 17% the catalytic activity of the native protein. The double mutant has catalytic activity comparable to the Lys-41Arg mutant.

The ribonuclease from an extinct bovid ruminant.

The sequence of the ribonuclease from the ancestor of swamp buffalo, river buffalo, and ox, corresponding approximately to Pachyportax latidens, an extinct ruminant known from the fossil record, has been reconstructed using the rule of 'maximum parsimony'. This protein and two sequences that may have been intermediates in the evolution of modern ribonuclease have been constructed in the laboratory by site-directed mutagenesis, and their properties examined.

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'.

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.

Synthesis and biodistribution of a short nonionic oligonucleotide analogue in mouse with a potential to mimic peptides.

A nonionic RNA analogue of the sequence r(USO2GSO2ASO2C) has been synthesized where each bridging phosphate diester is replaced by a dimethylene sulfone unit (rSNA). The rSNA was synthesized in solution from 3',5'-bishomo-beta-ribonucleoside derivatives as building blocks. Full experimental procedures are provided, and the product and all synthetic intermediates are fully characterized. The tetramer is nonionic but highly dipolar due to multiple hydrogen bonding opportunities. It is freely soluble in water only at higher pH's, permitting it to be radiolabeled by exchange of the acidic protons alpha to the sulfones with tritiated water. The tritiated molecule was administered intravenously into the tail vein (2.6 mg/kg) of mice, and its distribution was monitored over 48 h. The rSNA was widely distributed in the biological tissues, including the brain, and excreted in both the feces and the urine. The accumulation of radioactivity was significantly higher in liver and kidney than in other tissues. Radiolabel was recovered from the urine, analyzed by HPLC, and shown to be intact oligonucleotide sulfone. This is the first bioavailability study on a short nonionic oligonucleotide analogue, a class of molecules with potential biomedical applications.