Binding of phosphate and pyrophosphate ions at the active site of human angiogenin as revealed by X-ray crystallography

NMR Characterization of Angiogenin Variants and tRNAAla Products Impacting Aberrant Protein Oligomerization

Protein oligomerization is key to countless physiological processes, but also to abnormal amyloid conformations implicated in over 25 mortal human diseases. Human Angiogenin (h-ANG), a ribonuclease A family member, produces RNA fragments that regulate ribosome formation, the creation of new blood vessels and stress granule function. Too little h-ANG activity leads to abnormal protein oligomerization, resulting in Amyotrophic Lateral Sclerosis (ALS) or Parkinson's disease. While a score of disease linked h-ANG mutants has been studied by X-ray diffraction, some elude crystallization. There is also a debate regarding the structure that RNA fragments adopt after cleavage by h-ANG. Here, to better understand the beginning of the process that leads to aberrant protein oligomerization, the solution secondary structure and residue-level dynamics of WT h-ANG and two mutants i.e., H13A and R121C, are characterized by multidimensional heteronuclear NMR spectroscopy under near-physiological conditions. All three variants are found to adopt well folded and highly rigid structures in the solution, although the elements of secondary structure are somewhat shorter than those observed in crystallography studies. R121C alters the environment of nearby residues only. By contrast, the mutation H13A affects local residues as well as nearby active site residues K40 and H114. The conformation characterization by CD and 1D 1H NMR spectroscopies of tRNAAla before and after h-ANG cleavage reveals a retention of the duplex structure and little or no G-quadruplex formation.

Cryo-EM structures of the SARS-CoV-2 endoribonuclease Nsp15 reveal insight into nuclease specificity and dynamics

Nsp15, a uridine specific endoribonuclease conserved across coronaviruses, processes viral RNA to evade detection by host defense systems. Crystal structures of Nsp15 from different coronaviruses have shown a common hexameric assembly, yet how the enzyme recognizes and processes RNA remains poorly understood. Here we report a series of cryo-EM reconstructions of SARS-CoV-2 Nsp15, in both apo and UTP-bound states. The cryo-EM reconstructions, combined with biochemistry, mass spectrometry, and molecular dynamics, expose molecular details of how critical active site residues recognize uridine and facilitate catalysis of the phosphodiester bond. Mass spectrometry revealed the accumulation of cyclic phosphate cleavage products, while analysis of the apo and UTP-bound datasets revealed conformational dynamics not observed by crystal structures that are likely important to facilitate substrate recognition and regulate nuclease activity. Collectively, these findings advance understanding of how Nsp15 processes viral RNA and provide a structural framework for the development of new therapeutics.

How good are comparative models in the understanding of protein dynamics?

The 3D structure of a protein is essential to understand protein dynamics. If experimentally determined structure is unavailable, comparative models could be used to infer dynamics. However, the effectiveness of comparative models, compared to experimental structures, in inferring dynamics is not clear. To address this, we compared dynamics features of ~800 comparative models with their crystal structures using normal mode analysis. Average similarity in magnitude, direction, and correlation of residue motions is >0.8 (where value 1 is identical) indicating that the dynamics of models and crystal structures are highly similar. Accuracy of 3D structure and dynamics is significantly higher for models built on multiple and/or high sequence identity templates (>40%). Three-dimensional (3D) structure and residue fluctuations of models are closer to that of crystal structures than to templates (TM score 0.9 vs 0.7 and square inner product 0.92 vs 0.88). Furthermore, long-range molecular dynamics simulations on comparative models of RNase 1 and Angiogenin showed significant differences in the conformational sampling of conserved active-site residues that characterize differences in their activity levels. Similar analyses on two EGFR kinase variant models highlight the effect of mutations on the functional state-specific αC helix motions and these results corroborate with the previous experimental observations. Thus, our study adds confidence to the use of comparative models in understanding protein dynamics.

Evidence for Novel Action at the Cell-Binding Site of Human Angiogenin Revealed by Heteronuclear NMR Spectroscopy, in silico and in vivo Studies

A member of the ribonuclease A superfamily, human angiogenin (hAng) is a potent angiogenic factor. Heteronuclear NMR spectroscopy combined with induced-fit docking revealed a dual binding mode for the most antiangiogenic compound of a series of ribofuranosyl pyrimidine nucleosides that strongly inhibit hAng's angiogenic activity in vivo. While modeling suggests the potential for simultaneous binding of the inhibitors at the active and cell-binding sites, NMR studies indicate greater affinity for the cell-binding site than for the active site. Additionally, molecular dynamics simulations at 100 ns confirmed the stability of binding at the cell-binding site with the predicted protein-ligand interactions, in excellent agreement with the NMR data. This is the first time that a nucleoside inhibitor is reported to completely inhibit the angiogenic activity of hAng in vivo by exerting dual inhibitory activity on hAng, blocking both the entrance of hAng into the cell and its ribonucleolytic activity.

The ammonium sulfate inhibition of human angiogenin

In this study, we investigate the inhibition of human angiogenin by ammonium sulfate. The inhibitory potency of ammonium sulfate for human angiogenin (IC50 = 123.5 ± 14.9 mm) is comparable to that previously reported for RNase A (119.0 ± 6.5 mm) and RNase 2 (95.7 ± 9.3 mm). However, analysis of two X-ray crystal structures of human angiogenin in complex with sulfate anions (in acidic and basic pH environments, respectively) indicates an entirely distinct mechanism of inhibition. While ammonium sulfate inhibits the ribonucleolytic activity of RNase A and RNase 2 by binding to the active site of these enzymes, sulfate anions bind only to peripheral substrate anion-binding subsites of human angiogenin, and not to the active site.

New insights into the role of angiogenin in actin polymerization

Angiogenin is a potent stimulator of angiogenesis. It interacts with endothelial cells and induces a wide range of cellular responses initiating a process of blood vessel formation. One important target of angiogenin is endothelial cell-surface actin, and their interaction might be one of crucial steps in angiogenin-induced neovascularization. Recently, it was shown that angiogenin inhibits polymerization of G-actin and changes the physical properties of F-actin. These observations suggest that angiogenin may cause changes in the cell cytoskeleton. This chapter reviews the current state of the literature regarding angiogenin structure and function and discusses the relationship between the angiogenin and actin and possible functional roles of their interaction.

A new RNase sheds light on the RNase/angiogenin subfamily from zebrafish

Recently, extracellular RNases of the RNase A superfamily, with the characteristic CKxxNTF sequence signature, have been identified in fish. This has led to the recognition that these RNases are present in the whole vertebrate subphylum. In fact, they comprise the only enzyme family unique to vertebrates. Four RNases from zebrafish (Danio rerio) have been previously reported and have a very low RNase activity; some of these are endowed, like human angiogenin, with powerful angiogenic and bactericidal activities. In the present paper, we report the three-dimensional structure, the thermodynamic behaviour and the biological properties of a novel zebrafish RNase, ZF-RNase-5. The investigation of its structural and functional properties, extended to all other subfamily members, provides an inclusive description of the whole zebrafish RNase subfamily.

Influence of naturally-occurring 5'-pyrophosphate-linked substituents on the binding of adenylic inhibitors to ribonuclease a: an X-ray crystallographic study

Ribonuclease A is the archetype of a functionally diverse superfamily of vertebrate-specific ribonucleases. Inhibitors of its action have potential use in the elucidation of the in vivo roles of these enzymes and in the treatment of pathologies associated therewith. Derivatives of adenosine 5'-pyrophosphate are the most potent nucleotide-based inhibitors known. Here, we use X-ray crystallography to visualize the binding of four naturally-occurring derivatives that contain 5'-pyrophosphate-linked extensions. 5'-ATP binds with the adenine occupying the B(2) subsite in the manner of an RNA substrate but with the gamma-phosphate at the P(1) subsite. Diadenosine triphosphate (Ap(3)A) binds with the adenine in syn conformation, the beta-phosphate as the principal P(1) subsite ligand and without order beyond the gamma-phosphate. NADPH and NADP(+) bind with the adenine stacked against an alternative rotamer of His119, the 2'-phosphate at the P(1) subsite, and without order beyond the 5'-alpha-phosphate. We also present the structure of the complex formed with pyrophosphate ion. The structural data enable existing kinetic data on the binding of these compounds to a variety of ribonucleases to be rationalized and suggest that as the complexity of the 5'-linked extension increases, the need to avoid unfavorable contacts places limitations on the number of possible binding modes.

Structure of murine angiogenin: features of the substrate- and cell-binding regions and prospects for inhibitor-binding studies

Angiogenin is an unusual member of the pancreatic ribonuclease superfamily that induces blood-vessel formation and is a promising anticancer target. The three-dimensional structure of murine angiogenin (mAng) has been determined by X-ray crystallography. Two structures are presented: one is a complex with sulfate ions (1.5 Angstroms resolution) and the other a complex with phosphate ions (1.6 Angstroms resolution). Residues forming the putative B(1), P(1) and B(2) subsites occupy positions similar to their hAng counterparts and are likely to play similar roles. The anions occupy the P(1) subsite, sulfate binding conventionally and phosphate adopting two orientations, one of which is novel. The B(1) subsite is obstructed by Glu116 and Phe119, with the latter assuming a less invasive position than its hAng counterpart. Hydrophobic interactions between the C-terminal segment and the main body of the protein are more extensive than in hAng and may underly the lower enzymatic activity of the murine protein. Elsewhere, the structure of the H3-B2 loop supports the view that hAng Asn61 interacts directly with cell-surface molecules and does not merely stabilize adjacent regions of the hAng structure. mAng crystals appear to offer small-molecule inhibitors a clear route to the active site and may even withstand a reorientation of the C-terminal segment that provides access to the cryptic B(1) subsite. These features represent considerable advantages over crystalline hAng and bAng.

High-resolution crystal structures of ribonuclease A complexed with adenylic and uridylic nucleotide inhibitors. Implications for structure-based design of ribonucleolytic inhibitors

The crystal structures of bovine pancreatic ribonuclease A (RNase A) in complex with 3',5'-ADP, 2',5'-ADP, 5'-ADP, U-2'-p and U-3'-p have been determined at high resolution. The structures reveal that each inhibitor binds differently in the RNase A active site by anchoring a phosphate group in subsite P1. The most potent inhibitor of all five, 5'-ADP (Ki = 1.2 microM), adopts a syn conformation (in contrast to 3',5'-ADP and 2',5'-ADP, which adopt an anti), and it is the beta- rather than the alpha-phosphate group that binds to P1. 3',5'-ADP binds with the 5'-phosphate group in P1 and the adenosine in the B2 pocket. Two different binding modes are observed in the two RNase A molecules of the asymmetric unit for 2',5'-ADP. This inhibitor binds with either the 3' or the 5' phosphate groups in subsite P1, and in each case, the adenosine binds in two different positions within the B2 subsite. The two uridilyl inhibitors bind similarly with the uridine moiety in the B1 subsite but the placement of a different phosphate group in P1 (2' versus 3') has significant implications on their potency against RNase A. Comparative structural analysis of the RNase A, eosinophil-derived neurotoxin (EDN), eosinophil cationic protein (ECP), and human angiogenin (Ang) complexes with these and other phosphonucleotide inhibitors provides a wealth of information for structure-based design of inhibitors specific for each RNase. These inhibitors could be developed to therapeutic agents that could control the biological activities of EDN, ECP, and ANG, which play key roles in human pathologies.

Subtle functional collective motions in pancreatic-like ribonucleases: from ribonuclease A to angiogenin

The analysis of the dynamic behavior of enzymes is fundamental to structural biology. A direct relationship between protein flexibility and biological function has been shown for bovine pancreatic ribonuclease (RNase A) (Rasmussen et al., Nature 1992;357:423-424). More recently, crystallographic studies have shown that functional motions in RNase A involve the enzyme beta-sheet regions that move concertedly on substrate binding and release (Vitagliano et al., Proteins 2002;46:97-104). These motions have been shown to correspond to intrinsic dynamic properties of the native enzyme by molecular dynamics (MD) simulations. To unveil the occurrence of these collective motions in other members of pancreatic-like superfamily, we carried out MD simulations on human angiogenin (Ang). Essential dynamics (ED) analyses performed on the trajectories reveal that Ang exhibits collective motions similar to RNase A, despite the limited sequence identity (33%) of the two proteins. Furthermore, we show that these collective motions are also present in ensembles of experimentally determined structures of both Ang and RNase A. Finally, these subtle concerted beta-sheet motions were also observed for other two members of the pancreatic-like superfamily by comparing the ligand-bound and ligand-free structures of these enzymes. Taken together, these findings suggest that pancreatic-like ribonucleases share an evolutionary conserved dynamic behavior consisting of subtle beta-sheet motions, which are essential for substrate binding and release.

Inhibition of mammalian ribonucleases by endogenous adenosine dinucleotides

The most potent low molecular weight inhibitors of pancreatic RNase superfamily enzymes reported to date are synthetic derivatives of adenosine 5(')-pyrophosphate. Here we have investigated the effects of six natural nucleotides that also incorporate this moiety (NADP(+), NADPH, ATP, Ap(3)A, Ap(4)A, and Ap(5)A) on the activities of RNase A and two of its homologues, eosinophil-derived neurotoxin and angiogenin. With eosinophil-derived neurotoxin and angiogenin, Ap(5)A is comparable to the tightest binding inhibitors identified previously (K(i) values at pH 5.9 are 370 nM and 100 microM, respectively); it ranks among the strongest small antagonists of RNase A as well (K(i)=230 nM). The K(i) for NADPH with angiogenin is similar to that of Ap(5)A. These findings suggest that Ap(5)A and NADPH may serve as useful new leads for inhibitor design. Examination of inhibition under physiological conditions indicates that NADPH, ATP, and Ap(5)A may suppress intracellular RNase activity significantly in vivo.

Virtual screening to enrich hit lists from high-throughput screening: a case study on small-molecule inhibitors of angiogenin

"Hit lists" generated by high-throughput screening (HTS) typically contain a large percentage of false positives, making follow-up assays necessary to distinguish active from inactive substances. Here we present a method for improving the accuracy of HTS hit lists by computationally based virtual screening (VS) of the corresponding chemical libraries and selecting hits by HTS/VS consensus. This approach was applied in a case study on the target-enzyme angiogenin, a potent inducer of angiogenesis. In conjunction with HTS of the National Cancer Institute Diversity Set and ChemBridge DIVERSet E (approximately 18,000 compounds total), VS was performed with two flexible library docking/scoring methods, DockVision/Ludi and GOLD. Analysis of the results reveals that dramatic enrichment of the HTS hit rate can be achieved by selecting compounds in consensus with one or both of the VS functions. For example, HTS hits ranked in the top 2% by GOLD included 42% of the true hits, but only 8% of the false positives; this represents a sixfold enrichment over the HTS hit rate. Notably, the HTS/VS method was effective in selecting out inhibitors with midmicromolar dissociation constants typical of leads commonly obtained in primary screens.

A small-molecule inhibitor of the ribonucleolytic activity of human angiogenin that possesses antitumor activity

The results of previous preclinical and clinical studies have identified angiogenin (ANG) as a potentially important target for anticancer therapy. Here we report the design and implementation of a high-throughput screening assay to identify small molecules that bind to the ribonucleolytic active site of ANG, which is critically involved in the induction of angiogenesis by this protein. Screening of 18,310 compounds from the National Cancer Institute (NCI) Diversity Set and ChemBridge DIVERSet yielded 15 hits that inhibit the enzymatic activity of ANG with K(i) values <100 microM. One of these, NCI compound 65828 [8-amino-5-(4'-hydroxybiphenyl-4-ylazo)naphthalene-2-sulfonate; K(i) = 81 microM], was selected for more detailed studies. Minor changes in ANG or ligand structure markedly reduced potency, demonstrating that inhibition reflects active-site rather than nonspecific binding; these observations are consistent with a computationally generated model of the ANG.65828 complex. Local treatment with modest doses of 65828 significantly delayed the formation of s.c. tumors from two distinct human cancer cell types in athymic mice. ANG is the likely target involved because (i) a 65828 analogue with much lower potency against the enzymatic activity of ANG failed to exert any antitumor effect, (ii) tumors from 65828-treated mice had fewer interior blood vessels than those from control mice, and (iii) 65828 appears to have no direct effect on the tumor cells. Our findings provide considerable support for the targeting of the enzymatic active site of ANG as a strategy for developing new anticancer drugs.