Genetic selection for critical residues in ribonucleases

Effects of Pathogenic Mutants of the Neuroprotective RNase 5-Angiogenin in Amyotrophic Lateral Sclerosis (ALS)

Amyotrophic Lateral Sclerosis (ALS) is a fatal neurodegenerative disease that affects the motoneurons. More than 40 genes are related with ALS, and amyloidogenic proteins like SOD1 and/or TDP-43 mutants are directly involved in the onset of ALS through the formation of polymorphic fibrillogenic aggregates. However, efficacious therapeutic approaches are still lacking. Notably, heterozygous missense mutations affecting the gene coding for RNase 5, an enzyme also called angiogenin (ANG), were found to favor ALS onset. This is also true for the less-studied but angiogenic RNase 4. This review reports the substrate targets and illustrates the neuroprotective role of native ANG in the neo-vascularization of motoneurons. Then, it discusses the molecular determinants of many pathogenic ANG mutants, which almost always cause loss of function related to ALS, resulting in failures in angiogenesis and motoneuron protection. In addition, ANG mutations are sometimes combined with variants of other factors, thereby potentiating ALS effects. However, the activity of the native ANG enzyme should be finely balanced, and not excessive, to avoid possible harmful effects. Considering the interplay of these angiogenic RNases in many cellular processes, this review aims to stimulate further investigations to better elucidate the consequences of mutations in ANG and/or RNase 4 genes, in order to achieve early diagnosis and, possibly, successful therapies against ALS.

An ancient evolutionary connection between Ribonuclease A and EndoU families

The ribonuclease A family of proteins is well studied from the biochemical and biophysical points of view, but its evolutionary origins are obscure, as no sequences homologous to this family have been reported outside of vertebrates. Recently, the spatial structure of the ribonuclease domain from a bacterial polymorphic toxin was shown to be closely similar to the structure of vertebrate ribonuclease A. The absence of sequence similarity between the two structures prompted a speculation of convergent evolution of bacterial and vertebrate ribonuclease A-like enzymes. We show that bacterial and homologous archaeal polymorphic toxin ribonucleases with a known or predicted ribonuclease A-like fold are distant homologs of the ribonucleases from the EndoU family, found in all domains of cellular life and in viruses. We also detected a homolog of vertebrate ribonucleases A in the transcriptome assembly of the sea urchin Mesocentrotus franciscanus These observations argue for the common ancestry of prokaryotic ribonuclease A-like and ubiquitous EndoU-like ribonucleases, and suggest a better-grounded scenario for the origin of animal ribonucleases A, which could have emerged in the deuterostome lineage, either by an extensive modification of a copy of an EndoU gene, or, more likely, by a horizontal acquisition of a prokaryotic immunity-mediating ribonuclease gene.

Human angiogenin is a potent cytotoxin in the absence of ribonuclease inhibitor

Angiogenin (ANG) is a secretory ribonuclease that promotes the proliferation of endothelial cells, leading to angiogenesis. This function relies on its ribonucleolytic activity, which is low for simple RNA substrates. Upon entry into the cytosol, ANG is sequestered by the ribonuclease inhibitor protein (RNH1). We find that ANG is a potent cytotoxin for RNH1-knockout HeLa cells, belying its inefficiency as a nonspecific catalyst. The toxicity does, however, rely on the ribonucleolytic activity of ANG and a cytosolic localization, which lead to the accumulation of particular tRNA fragments (tRFs), such as tRF-5 Gly-GCC. These up-regulated tRFs are highly cytotoxic at physiological concentrations. Although ANG is well-known for its promotion of cell growth, our results reveal that ANG can also cause cell death.

Conserved amino acid networks modulate discrete functional properties in an enzyme superfamily

In this work, we applied the sequence-based statistical coupling analysis approach to characterize conserved amino acid networks important for biochemical function in the pancreatic-type ribonuclease (ptRNase) superfamily. This superfamily-wide analysis indicates a decomposition of the RNase tertiary structure into spatially distributed yet physically connected networks of co-evolving amino acids, termed sectors. Comparison of this statistics-based description with new NMR experiments data shows that discrete amino acid networks, termed sectors, control the tuning of distinct functional properties in different enzyme homologs. Further, experimental characterization of evolutionarily distant sequences reveals that sequence variation at sector positions can distinguish homologs with a conserved dynamic pattern and optimal catalytic activity from those with altered dynamics and diminished catalytic activities. Taken together, these results provide important insights into the mechanistic design of the ptRNase superfamily, and presents a structural basis for evolutionary tuning of function in functionally diverse enzyme homologs.

Molecular basis for the autonomous promotion of cell proliferation by angiogenin

Canonical growth factors act indirectly via receptor-mediated signal transduction pathways. Here, we report on an autonomous pathway in which a growth factor is internalized, has its localization regulated by phosphorylation, and ultimately uses intrinsic catalytic activity to effect epigenetic change. Angiogenin (ANG), a secreted vertebrate ribonuclease, is known to promote cell proliferation, leading to neovascularization as well as neuroprotection in mammals. Upon entering cells, ANG encounters the cytosolic ribonuclease inhibitor protein, which binds with femtomolar affinity. We find that protein kinase C and cyclin-dependent kinase phosphorylate ANG, enabling ANG to evade its inhibitor and enter the nucleus. After migrating to the nucleolus, ANG cleaves promoter-associated RNA, which prevents the recruitment of the nucleolar remodeling complex to the ribosomal DNA promoter. The ensuing derepression of rDNA transcription promotes cell proliferation. The biochemical basis for this unprecedented mechanism of signal transduction suggests new modalities for the treatment of cancers and neurological disorders.

Perturbation of the Conformational Dynamics of an Active-Site Loop Alters Enzyme Activity

The role of internal dynamics in enzyme function is highly debated. Specifically, how small changes in structure far away from the reaction site alter protein dynamics and overall enzyme mechanisms is of wide interest in protein engineering. Using RNase A as a model, we demonstrate that elimination of a single methyl group located >10 Å away from the reaction site significantly alters conformational integrity and binding properties of the enzyme. This A109G mutation does not perturb structure or thermodynamic stability, both in the apo and ligand-bound states. However, significant enhancement in conformational dynamics was observed for the bound variant, as probed over nano- to millisecond timescales, resulting in major ligand repositioning. These results illustrate the large effects caused by small changes in structure on long-range conformational dynamics and ligand specificities within proteins, further supporting the importance of preserving wild-type dynamics in enzyme systems that rely on flexibility for function.

Network of long-range concerted chemical shift displacements upon ligand binding to human angiogenin

Molecular recognition models of both induced fit and conformational selection rely on coupled networks of flexible residues and/or structural rearrangements to promote protein function. While the atomic details of these motional events still remain elusive, members of the pancreatic ribonuclease superfamily were previously shown to depend on subtle conformational heterogeneity for optimal catalytic function. Human angiogenin, a structural homologue of bovine pancreatic RNase A, induces blood vessel formation and relies on a weak yet functionally mandatory ribonucleolytic activity to promote neovascularization. Here, we use the NMR chemical shift projection analysis (CHESPA) to clarify the mechanism of ligand binding in human angiogenin, further providing information on long-range intramolecular residue networks potentially involved in the function of this enzyme. We identify two main clusters of residue networks displaying correlated linear chemical shift trajectories upon binding of substrate fragments to the purine- and pyrimidine-specific subsites of the catalytic cleft. A large correlated residue network clusters in the region corresponding to the V1 domain, a site generally associated with the angiogenic response and structural stability of the enzyme. Another correlated network (residues 40-42) negatively affects the catalytic activity but also increases the angiogenic activity. (15) N-CPMG relaxation dispersion experiments could not reveal the existence of millisecond timescale conformational exchange in this enzyme, a lack of flexibility supported by the very low-binding affinities and catalytic activity of angiogenin. Altogether, the current report potentially highlights the existence of long-range dynamic reorganization of the structure upon distinct subsite binding events in human angiogenin.

Altered endoribonuclease activity of apurinic/apyrimidinic endonuclease 1 variants identified in the human population

Apurinic/apyrimidinic endonuclease 1 (APE1) is the major mammalian enzyme in the DNA base excision repair pathway and cleaves the DNA phosphodiester backbone immediately 5′ to abasic sites. APE1 also has 3′-5′ DNA exonuclease and 3′ DNA phosphodiesterase activities, and regulates transcription factor DNA binding through its redox regulatory function. The human APE1 has recently been shown to endonucleolytically cleave single-stranded regions of RNA. Towards understanding the biological significance of the endoribonuclease activity of APE1, we examined eight different amino acid substitution variants of APE1 previously identified in the human population. Our study shows that six APE1 variants, D148E, Q51H, I64V, G241R, R237A, and G306A, exhibit a 76–85% reduction in endoribonuclease activity against a specific coding region of the c-myc RNA, yet fully retain the ability to cleave apurinic/apyrimidinic DNA. We found that two APE1 variants, L104R and E126D, exhibit a unique RNase inhibitor-resistant endoribonuclease activity, where the proteins cleave c-myc RNA 3′ of specific single-stranded guanosine residues. Expression of L104R and E126D APE1 variants in bacterial Origami cells leads to a 60–80% reduction in colony formation and a 1.5-fold increase in cell doubling time, whereas the other variants, which exhibit diminished endoribonuclease activity, had no effect. These data indicate that two human APE1 variants exhibit a unique endoribonuclease activity, which correlates with their ability to induce cytotoxicity or slow down growth in bacterial cells and supports the notion of their biological functionality.

Identification of a novel missense mutation in angiogenin in a Chinese amyotrophic lateral sclerosis cohort

Abstract Angiogenin (ANG) gene mutations have been identified in both familial and sporadic amyotrophic lateral sclerosis (ALS) patients from multiple European and North American populations. However, no ANG mutation has yet been reported in Asian ALS populations. Here, we screened for ANG mutations in a Chinese ALS cohort. The entire coding region of the ANG gene was sequenced in 10 familial ALS pedigrees, 202 sporadic ALS patients, and 151 healthy controls. All patients were negative for SOD1, FUS, and TARDBP mutations. We identified a novel missense mutation, c.379G > A (p.V103I), in one sporadic ALS patient, but not in the controls. No mutations were found in the familial ALS patients. A novel missense variant, c.323A > G (p.H84R), was detected in one healthy individual. We identified the presence of the known single nucleotide polymorphism, rs11701 (T/G), in both ALS cases and controls. However, no significant association of the G allele with ALS susceptibility was demonstrated. In conclusion, ANG mutations accounted for 0.5% of our SOD1-, FUS-, TARDBP- mutation-negative ALS cohort. Our findings highlight that the genetic background of ALS differs between different populations, and suggest that ANG mutation may be involved in the aetiology of ALS in the Han Chinese population.

Rational design and evaluation of mammalian ribonuclease cytotoxins

Mammalian pancreatic-type ribonucleases (ptRNases) comprise an enzyme family that is remarkably well suited for therapeutic exploitation. ptRNases are robust and prodigious catalysts of RNA cleavage that can naturally access the cytosol. Instilling cytotoxic activity requires endowing them with the ability to evade a cytosolic inhibitor protein while retaining other key attributes. These efforts have informed our understanding of ptRNase-based cytotoxins, as well as the action of protein-based drugs with cytosolic targets. Here, we address the most pressing problems encountered in the design of cytotoxic ptRNases, along with potential solutions. In addition, we describe assays that can be used to evaluate a successful design in vitro, in cellulo, and in vivo. The emerging information validates the continuing development of ptRNases as chemotherapeutic agents.

Ribonuclease A suggests how proteins self-chaperone against amyloid fiber formation

Genomic analyses have identified segments with high fiber-forming propensity in many proteins not known to form amyloid. Proteins are often protected from entering the amyloid state by molecular chaperones that permit them to fold in isolation from identical molecules; but, how do proteins self-chaperone their folding in the absence of chaperones? Here, we explore this question with the stable protein ribonuclease A (RNase A). We previously identified fiber-forming segments of amyloid-related proteins and demonstrated that insertion of these segments into the C-terminal hinge loop of nonfiber-forming RNase A can convert RNase A into the amyloid state through three-dimensional domain-swapping, where the inserted fiber-forming segments interact to create a steric zipper spine. In this study, we convert RNase A into amyloid-like fibers by increasing the loop length and hence conformational freedom of an endogenous fiber-forming segment, SSTSAASS, in the N-terminal hinge loop. This is accomplished by sandwiching SSTSAASS between inserted Gly residues. With these inserts, SSTSAASS is now able to form the steric zipper spine, allowing RNase A to form amyloid-like fibers. We show that these fibers contain RNase A molecules retaining their enzymatic activity and therefore native-like structure. Thus, RNase A appears to prevent fiber formation by limiting the conformational freedom of this fiber-forming segment from entering a steric zipper. Our observations suggest that proteins have evolved to self-chaperone by using similar protective mechanisms.

The relationship between relative solvent accessibility and evolutionary rate in protein evolution

Recent work with Saccharomyces cerevisiae shows a linear relationship between the evolutionary rate of sites and the relative solvent accessibility (RSA) of the corresponding residues in the folded protein. Here, we aim to develop a mathematical model that can reproduce this linear relationship. We first demonstrate that two models that both seem reasonable choices (a simple model in which selection strength correlates with RSA and a more complex model based on RSA-dependent amino acid distributions) fail to reproduce the observed relationship. We then develop a model on the basis of observed site-specific amino acid distributions and show that this model behaves appropriately. We conclude that evolutionary rates are directly linked to the distribution of amino acids at individual sites. Because of this link, any future insight into the biophysical mechanisms that determine amino acid distributions will improve our understanding of evolutionary rates.

Alteration of hydrogen bonding in the vicinity of histidine 48 disrupts millisecond motions in RNase A

The motion of amino acid residues on the millisecond (ms) time scale is involved in the tight regulation of catalytic function in numerous enzyme systems. Using a combination of mutational, enzymological, and relaxation-compensated (15)N Carr-Purcell-Meiboom-Gill (CPMG) methods, we have previously established the conformational significance of the distant His48 residue and the neighboring loop 1 in RNase A function. These studies suggested that RNase A relies on an intricate network of hydrogen bonding interactions involved in propagating functionally relevant, long-range ms motions to the catalytic site of the enzyme. To further investigate the dynamic importance of this H-bonding network, this study focuses on the individual replacement of Thr17 and Thr82 with alanine, effectively altering the key H-bonding interactions that connect loop 1 and His48 to the rest of the protein. (15)N CPMG dispersion studies, nuclear magnetic resonance (NMR) chemical shift analysis, and NMR line shape analysis of point mutants T17A and T82A demonstrate that the evolutionarily conserved single H-bond linking His48 to Thr82 is essential for propagating ms motions from His48 to the active site of RNase A on the time scale of catalytic turnover, whereas the T17A mutation increases the off rate and conformational exchange motions in loop 1. Accumulating evidence from our mutational studies indicates that residues experiencing conformational exchange in RNase A can be grouped into two separate clusters displaying distinct dynamical features, which appear to be independently affected by mutation. Overall, this study illuminates how tightly controlled and finely tuned ms motions are in RNase A, suggesting that designed modulation of protein motions may be possible.

In planta mutagenesis determines the functional regions of the wheat puroindoline proteins

In planta analysis of protein function in a crop plant could lead to improvements in understanding protein structure/function relationships as well as selective agronomic or end product quality improvements. The requirements for successful in planta analysis are a high mutation rate, an efficient screening method, and a trait with high heritability. Two ideal targets for functional analysis are the Puroindoline a and Puroindoline b (Pina and Pinb, respectively) genes, which together compose the wheat (Triticum aestivum L.) Ha locus that controls grain texture and many wheat end-use properties. Puroindolines (PINs) together impart soft texture, and mutations in either PIN result in hard seed texture. Studies of the PINs' mode of action are limited by low allelic variation. To create new Pin alleles and identify critical function-determining regions, Pin point mutations were created in planta via EMS treatment of a soft wheat. Grain hardness of 46 unique PIN missense alleles was then measured using segregating F(2):F(3) populations. The impact of individual missense alleles upon PIN function, as measured by grain hardness, ranged from neutral (74%) to intermediate to function abolishing. The percentage of function-abolishing mutations among mutations occurring in both PINA and PINB was higher for PINB, indicating that PINB is more critical to overall Ha function. This is contrary to expectations in that PINB is not as well conserved as PINA. All function-abolishing mutations resulted from structure-disrupting mutations or from missense mutations occurring near the Tryptophan-rich region. This study demonstrates the feasibility of in planta functional analysis of wheat proteins and that the Tryptophan-rich region is the most important region of both PINA and PINB.

A model for protein sequence evolution based on selective pressure for protein stability: application to hemoglobins

Negative selection against protein instability is a central influence on evolution of proteins. Protein stability is maintained over evolution despite changes in underlying sequences. An empirical all-site stability-based model of evolution was developed to focus on the selection of residues arising from their contributions to protein stability. In this model, site rates could vary. A structure-based method was used to predict stationary frequencies of hemoglobin residues based on their propensity to promote protein stability at a site. Sites with destabilizing residues were shown to change more rapidly in hemoglobins than sites with stabilizing residues. For diverse proteins the results were consistent with stability-based selection. Maximum likelihood studies with hemoglobins supported the stability-based model over simple Poisson-based methods. These observations are consistent with suggestions that purifying selection to maintain protein structural stability plays a dominant role in protein evolution.

Short protein segments can drive a non-fibrillizing protein into the amyloid state

Protein fibrils termed amyloid-like are associated with numerous degenerative diseases as well as some normal cellular functions. Specific short segments of amyloid-forming proteins have been shown to form fibrils themselves. However, it has not been shown in general that these segments are capable of driving a protein from its native structure into the amyloid state. We applied the 3D profile method to identify fibril-forming segments within the amyloid-forming human proteins tau, alpha-synuclein, PrP prion and amyloid-beta. Ten segments, six to eight residues in length, were chosen and inserted into the C-terminal hinge loop of the highly constrained enzyme RNase A, and tested for fibril growth and Congo red birefringence. We find that all 10 unique inserts cause RNase A to form amyloid-like fibrils which display characteristic yellow to apple-green Congo red birefringence when observed with cross polarizers. These six to eight residue inserts can fibrillize RNase A and are sufficient for amyloid fibril spine formation.

Ribonuclease inhibitor regulates neovascularization by human angiogenin

Human angiogenin (ANG) is a homologue of bovine pancreatic ribonuclease (RNase A) that induces neovascularization. ANG is the only human angiogenic factor that possesses ribonucleolytic activity. To stimulate blood vessel growth, ANG must be transported to the nucleus and must retain its catalytic activity. Like other mammalian homologues of RNase A, ANG forms a femtomolar complex with the cytosolic ribonuclease inhibitor protein (RI). To determine whether RI affects ANG-induced angiogenesis, we created G85R/G86R ANG, which possesses 10(6)-fold lower affinity for RI but retains wild-type ribonucleolytic activity. The neovascularization of rabbit corneas by G85R/G86R ANG was more pronounced and more rapid than by wild-type ANG. These findings provide the first direct evidence that RI serves to regulate the biological activity of ANG in vivo.

Relative tolerance of an enzymatic molten globule and its thermostable counterpart to point mutation

Enzyme structures reflect the complex interplay between the free energy of unfolding (DeltaG) and catalytic efficiency. Consequently, the effects of point mutations on structure, stability, and function are difficult to predict. It has been proposed that the mutational robustness of homologous enzymes correlates with a higher initial DeltaG. To examine this issue, we compared the tolerance of a natural thermostable chorismate mutase and an engineered molten globular variant to targeted mutation. These mutases possess similar sequence, structure, and catalytic efficiency but dramatically different DeltaG values. We find that analogous point mutations can have widely divergent effects on catalytic activity in these scaffolds. In a set of five rationally designed single-amino acid changes, the thermostable scaffold suffers activity losses ranging from 50-fold smaller, for an aspartate-to-glycine substitution at the active site, to 2-fold greater, for a phenylalanine-to-tryptophan substitution in the hydrophobic core, versus that of the molten globular scaffold. However, biophysical characterization indicates that the variations in catalytic efficiency are not caused by losses of either secondary structural integrity or thermodynamic stability. Rather, the activity differences between variant pairs are very much context-dependent and likely stem from subtle changes in the fine structure of the active site. Thus, in many cases, it may be more productive to focus on changes in local conformation than on global stability when attempting to understand and predict how enzymes respond to point mutations.

Detection of SNPs in fish DNA: application of the fluorogenic ribonuclease protection (FRIP) assay for the authentication of food contents

The fluorogenic ribonuclease protection (FRIP) assay was used to detect single nucleotide polymorphisms (SNPs) in commercially produced fish products. By using fluorescence resonance energy transfer (FRET) between fluorophore and quencher labeled probes, the species-specific cleavage of sample RNA was detected by measuring the fluorescence intensity during the FRIP assay. We were able to discriminate raw and thermally processed eel and tuna species using the FRIP-based SNP detection method. Furthermore, the intensity of fluorescence was correlated with the mutant/wild-type ratio. These results suggest that the FRIP assay is a useful method for the in situ confirmation of labels of fishery foods during food production.

Genetic selection for peptide inhibitors of angiogenin

The improper regulation of angiogenesis is implicit in a variety of diseases, including cancer. Angiogenin is unique among angiogenic factors in having ribonucleolytic activity. Inhibitors of this activity could serve as chemotherapeutics. The ribonucleolytic activity of angiogenin is toxic to the Origami strain of Escherichia coli. Herein, this cytotoxicity was used to identify inhibitors from a random nonapeptide library tethered to the C-terminus of human angiogenin. The selected sequences fell into three classes: (i) extremely hydrophobic, (ii) putative protease (ClpXP) substrates and (iii) slightly anionic. Two peptides corresponding to sequences in the last class were synthesized chemically and found to inhibit the ribonucleolytic activity of human angiogenin in vitro with micromolar values of Ki. Both peptides also inhibit bovine pancreatic ribonuclease, a homolog of angiogenin, though one exhibits selectivity for angiogenin. The affinity and selectivity of these peptides are comparable with the best known inhibitors of angiogenin. Moreover, the strategy used to identify them is general and could be applied to other cytotoxins.

Intraspecies regulation of ribonucleolytic activity

The evolutionary rate of proteins involved in obligate protein-protein interactions is slower and the degree of coevolution higher than that for nonobligate protein-protein interactions. The coevolution of the proteins involved in certain nonobligate interactions is, however, essential to cell survival. To gain insight into the coevolution of one such nonobligate protein pair, the cytosolic ribonuclease inhibitor (RI) proteins and secretory pancreatic-type ribonucleases from cow (Bos taurus) and human (Homo sapiens) were produced in Escherichia coli and purified, and their physicochemical properties were analyzed. The two intraspecies complexes were found to be extremely tight (bovine Kd = 0.69 fM; human Kd = 0.34 fM). Human RI binds to its cognate ribonuclease (RNase 1) with 100-fold greater affinity than to the bovine homologue (RNase A). In contrast, bovine RI binds to RNase 1 and RNase A with nearly equal affinity. This broader specificity is consistent with there being more pancreatic-type ribonucleases in cows (20) than humans (13). Human RI (32 cysteine residues) also has 4-fold less resistance to oxidation by hydrogen peroxide than does bovine RI (29 cysteine residues). This decreased oxidative stability of human RI, which is caused largely by Cys74, implies a larger role for human RI as an antioxidant. The conformational and oxidative stabilities of both RIs increase upon complex formation with ribonucleases. Thus, RI has evolved to maintain its inhibition of invading ribonucleases, even when confronted with extreme environmental stress. That role appears to take precedence over its role in mediating oxidative damage.

Structural and functional implications of positive selection at the primate angiogenin gene

BACKGROUND

Angiogenesis, the formation of new blood vessels, is a primordial process in development and its dysregulation has a central role in the pathogenesis of many diseases. Angiogenin (ANG), a peculiar member of the RNase A superfamily, is a potent inducer of angiogenesis involved in many different types of cancer, amyotrophic lateral sclerosis and also with a possible role in the innate immune defense. The evolutionary path of this family has been a highly dynamic one, where positive selection has played a strong role. In this work we used a combined gene and protein level approach to determine the main sites under diversifying selection on the primate ANG gene and analyze its structural and functional implications.

RESULTS

We obtained evidence for positive selection in the primate ANG gene. Site specific analysis pointed out 15 sites under positive selection, most of which also exhibited drastic changes in amino acid properties. The mapping of these sites in the ANG 3D-structure described five clusters, four of which were located in functional regions: two in the active site region, one in the nucleolar location signal and one in the cell-binding site. Eight of the 15 sites under selection in the primate ANG gene were highly or moderately conserved in the RNase A family, suggesting a directed event and not a simple consequence of local structural or functional permissiveness. Moreover, 11 sites were exposed to the surface of the protein indicating that they may influence the interactions performed by ANG.

CONCLUSION

Using a maximum likelihood gene level analysis we identified 15 sites under positive selection in the primate ANG genes, that were further corroborated through a protein level analysis of radical changes in amino acid properties. These sites mapped onto the main functional regions of the ANG protein. The fact that evidence for positive selection is present in all ANG regions required for angiogenesis may be a good indication that angiogenesis is the process under selection. However, other possibilities to be considered arise from the possible involvement of ANG in innate immunity and the potential influence or co-evolution with its interacting proteins and ligands.

Back to the future: Ribonuclease A

Pancreatic ribonuclease A (EC 3.1.27.5, RNase) is, perhaps, the best-studied enzyme of the 20th century. It was isolated by René Dubos, crystallized by Moses Kunitz, sequenced by Stanford Moore and William Stein, and synthesized in the laboratory of Bruce Merrifield, all at the Rockefeller Institute/University. It has proven to be an excellent model system for many different types of experiments, both as an enzyme and as a well-characterized protein for biophysical studies. Of major significance was the demonstration by Chris Anfinsen at NIH that the primary sequence of RNase encoded the three-dimensional structure of the enzyme. Many other prominent protein chemists/enzymologists have utilized RNase as a dominant theme in their research. In this review, the history of RNase and its offspring, RNase S (S-protein/S-peptide), will be considered, especially the work in the Merrifield group, as a preface to preliminary data and proposed experiments addressing topics of current interest. These include entropy-enthalpy compensation, entropy of ligand binding, the impact of protein modification on thermal stability, and the role of protein dynamics in enzyme action. In continuing to use RNase as a prototypical enzyme, we stand on the shoulders of the giants of protein chemistry to survey the future.

Genetic selection reveals the role of a buried, conserved polar residue

The burial of nonpolar surface area is known to enhance markedly the conformational stability of proteins. The contribution from the burial of polar surface area is less clear. Here, we report on the tolerance to substitution of Ser75 of bovine pancreatic ribonuclease (RNase A), a residue that has the unusual attributes of being buried, conserved, and polar. To identify variants that retain biological function, we used a genetic selection based on the intrinsic cytotoxicity of ribonucleolytic activity. Cell growth at 30 degrees C, 37 degrees C, and 44 degrees C correlated with residue size, indicating that the primary attribute of Ser75 is its small size. The side-chain hydroxyl group of Ser75 forms a hydrogen bond with a main-chain nitrogen. The conformational stability of the S75A variant, which lacks this hydrogen bond, was diminished by DeltaDeltaG = 2.5 kcal/mol. Threonine, which can reinstate this hydrogen bond, provided a catalytically active RNase A variant at higher temperatures than did some smaller residues (including aspartate), indicating that a secondary attribute of Ser75 is the ability of its uncharged side chain to accept a hydrogen bond. These results provide insight on the imperatives for the conservation of a buried polar residue.