Contribution of the active site histidine residues of ribonuclease A to nucleic acid binding

Folding Kinetics and Volume Variation of the β-Hairpin Peptide Chignolin upon Ultrafast pH-Jumps

In-depth characterization of fundamental folding steps of small model peptides is crucial for a better understanding of the folding mechanisms of more complex biomacromolecules. We have previously reported on the folding/unfolding kinetics of a model α-helix. Here, we study folding transitions in chignolin (GYDPETGTWG), a short β-hairpin peptide previously used as a model to study conformational changes in β-sheet proteins. Although previously suggested, until now, the role of the Tyr2-Trp9 interaction in the folding mechanism of chignolin was not clear. In the present work, pH-dependent conformational changes of chignolin were characterized by circular dichroism (CD), nuclear magnetic resonance (NMR), ultrafast pH-jump coupled with time-resolved photoacoustic calorimetry (TR-PAC), and molecular dynamics (MD) simulations. Taken together, our results present a comprehensive view of chignolin's folding kinetics upon local pH changes and the role of the Tyr2-Trp9 interaction in the folding process. CD data show that chignolin's β-hairpin formation displays a pH-dependent skew bell-shaped curve, with a maximum close to pH 6, and a large decrease in β-sheet content at alkaline pH. The β-hairpin structure is mainly stabilized by aromatic interactions between Tyr2 and Trp9 and CH-π interactions between Tyr2 and Pro4. Unfolding of chignolin at high pH demonstrates that protonation of Tyr2 is essential for the stability of the β-hairpin. Refolding studies were triggered by laser-induced pH-jumps and detected by TR-PAC. The refolding of chignolin from high pH, mainly due to the protonation of Tyr2, is characterized by a volume expansion (10.4 mL mol-1), independent of peptide concentration, in the microsecond time range (lifetime of 1.15 μs). At high pH, the presence of the deprotonated hydroxyl (tyrosinate) hinders the formation of the aromatic interaction between Tyr2 and Trp9 resulting in a more disorganized and dynamic tridimensional structure of the peptide. This was also confirmed by comparing MD simulations of chignolin under conditions mimicking neutral and high pH.

Sugar osmolyte inhibits and attenuates the fibrillogenesis in RNase A: An in vitro and in silico characterizations

Mechanisms of protein aggregation are of immense interest in therapeutic biology and neurodegenerative medicine. Biochemical processes within the living cell occur in a highly crowded environment. The phenomenon of macromolecular crowding affects the diffusional and conformational dynamics of proteins and modulates their folding. Macromolecular crowding is reported to cause protein aggregation in some cases, so it is a cause of concern as it leads to a plethora of neurodegenerative disorders and systemic amyloidosis. To divulge the mechanism of aggregation, it is imperative to study aggregation in well-characterized model proteins in the presence of macromolecular crowder. One such protein is ribonuclease A (RNase A), which deciphers neurotoxic function in humans; therefore we decided to explore the amyloid fibrillogenesis of this thermodynamically stable protein. To elucidate the impact of crowder, dextran-70 and its monomer glucose on the aggregation profile of RNase-A various techniques such as Absorbance, Fluorescence, Fourier Transforms Infrared, Dynamic Light Scattering and circular Dichroism spectroscopies along with imaging techniques like Atomic Force Microscopy and Transmission Electron Microscopy were employed. Thermal aggregation and fibrillation were further promoted by dextran-70 while glucose counteracted the effect of the crowding agent in a concentration-dependent manner. This study shows that glucose provides stability to the protein and prevents fibrillation. Intending to combat aggregation, which is the hallmark of numerous late-onset neurological disorders and systemic amyloidosis, this investigation unveils that naturally occurring osmolytes or other co-solutes can be further exploited in novel drug design strategies.

Impact of Peptide Sequence on Functional siRNA Delivery and Gene Knockdown with Cyclic Amphipathic Peptide Delivery Agents

Short-interfering RNA (siRNA) oligonucleotide therapeutics that modify gene expression by accessing RNA-interference (RNAi) pathways have great promise for the treatment of a range of disorders; however, their application in clinical settings has been limited by significant challenges in cellular delivery. Herein, we report a structure-function study using a series of modified cyclic amphipathic cell-penetrating peptides (CAPs) to determine the impact of peptide sequence on (1) siRNA-binding efficiency, (2) cellular delivery and knockdown efficiency, and (3) the endocytic uptake mechanism. Nine cyclic peptides of the general sequence Ac-C[XZ]4CG-NH2 in which X residues are hydrophobic/aromatic (Phe, Tyr, Trp, or Leu) and Z residues are charged/hydrophilic (Arg, Lys, Ser, or Glu) are assessed along with one acyclic peptide, Ac-(WR)4G-NH2. Cyclization is enforced by intramolecular disulfide bond formation between the flanking Cys residues. Binding analyses indicate that strong cationic character and the presence of aromatic residues that are competent to participate in CH-π interactions lead to CAP sequences that most effectively interact with siRNA. CAP-siRNA binding increases in the following order as a function of CAP hydrophobic/aromatic content: His < Phe < Tyr < Trp. Both cationic charge and disulfide-constrained cyclization of CAPs improve uptake of siRNA in vitro. Net neutral CAPs and an acyclic peptide demonstrate less-efficient siRNA translocation compared to the cyclic, cationic CAPs tested. All CAPs tested facilitated efficient siRNA target gene knockdown of at least 50% (as effective as a lipofectamine control), with the best CAPs enabling >80% knockdown. Significantly, gene knockdown efficiency does not strongly correlate with CAP-siRNA internalization efficiency but moderately correlates with CAP-siRNA-binding affinity. Finally, utilization of small-molecule inhibitors and targeted knockdown of essential endocytic pathway proteins indicate that most CAP-siRNA nanoparticles facilitate siRNA delivery through clathrin- and caveolin-mediated endocytosis. These results provide insight into the design principles for CAPs to facilitate siRNA delivery and the mechanisms by which these peptides translocate siRNA into cells. These studies also demonstrate the nature of the relationships between peptide-siRNA binding, cellular delivery of siRNA cargo, and functional gene knockdown. Strong correlations between these properties are not always observed, which illustrates the complexity in the design of optimal next-generation materials for oligonucleotide delivery.

A method for quantifying how the activity of an enzyme is affected by the net charge of its nearest crowded neighbor

The electrostatic effects of protein crowding have not been systematically explored. Rather, protein crowding is generally studied with co‐solvents or crowders that are electrostatically neutral, with no methods to measure how the net charge (Z) of a crowder affects protein function. For example, can the activity of an enzyme be affected electrostatically by the net charge of its neighbor in crowded milieu? This paper reports a method for crowding proteins of different net charge to an enzyme via semi‐random chemical crosslinking. As a proof of concept, RNase A was crowded (at distances ≤ the Debye length) via crosslinking to different heme proteins with Z = +8.50 ± 0.04, Z = +6.39 ± 0.12, or Z = −10.30 ± 1.32. Crosslinking did not disrupt the structure of proteins, according to amide H/D exchange, and did not inhibit RNase A activity. For RNase A, we found that the electrostatic environment of each crowded neighbor had significant effects on rates of RNA hydrolysis. Crowding with cationic cytochrome c led to increases in activity, while crowding with anionic “supercharged” cytochrome c or myoglobin diminished activity. Surprisingly, electrostatic crowding effects were amplified at high ionic strength (I = 0.201 M) and attenuated at low ionic strength (I = 0.011 M). This salt dependence might be caused by a unique set of electric double layers at the dimer interspace (maximum distance of 8 Å, which cannot accommodate four layers). This new method of crowding via crosslinking can be used to search for electrostatic effects in protein crowding.

Modifications of Ribonucleases in Order to Enhance Cytotoxicity in Anticancer Therapy

Ribonucleases (RNases) are a superfamily of enzymes that have been extensively studied since the 1960s. For a long time, this group of secretory enzymes was studied as an important model for protein chemistry such as folding, stability and enzymatic catalysis. Since it was discovered that RNases displayed cytotoxic activity against several types of malignant cells, recent investigation has focused mainly on the biological functions and medical applications of engineered RNases. In this review, we describe structures, functions and mechanisms of antitumor activity of RNases. They operate at the crossroads of transcription and translation, preferentially degrading tRNA. As a result, this inhibits protein synthesis, induces apoptosis and causes death of cancer cells. This effect can be enhanced thousands of times when RNases are conjugated with monoclonal antibodies. Such combinations, called immunoRNases, have demonstrated selective antitumor activity against cancer cells both in vitro and in animal models. This review summarizes the current status of engineered RNases and immunoRNases as promising novel therapeutic agents for different types of cancer. Also, we describe our experimental results from published or previously unpublished research and compare with other scientific information.

Human ribonuclease 1 serves as a secretory ligand of ephrin A4 receptor and induces breast tumor initiation

Human ribonuclease 1 (hRNase 1) is critical to extracellular RNA clearance and innate immunity to achieve homeostasis and host defense; however, whether it plays a role in cancer remains elusive. Here, we demonstrate that hRNase 1, independently of its ribonucleolytic activity, enriches the stem-like cell population and enhances the tumor-initiating ability of breast cancer cells. Specifically, secretory hRNase 1 binds to and activates the tyrosine kinase receptor ephrin A4 (EphA4) signaling to promote breast tumor initiation in an autocrine/paracrine manner, which is distinct from the classical EphA4-ephrin juxtacrine signaling through contact-dependent cell-cell communication. In addition, analysis of human breast tumor tissue microarrays reveals a positive correlation between hRNase 1, EphA4 activation, and stem cell marker CD133. Notably, high hRNase 1 level in plasma samples is positively associated with EphA4 activation in tumor tissues from breast cancer patients, highlighting the pathological relevance of the hRNase 1-EphA4 axis in breast cancer. The discovery of hRNase 1 as a secretory ligand of EphA4 that enhances breast cancer stemness suggests a potential treatment strategy by inactivating the hRNase 1-EphA4 axis.

Highly Selective Stable Isotope Labeling of Histidine Residues by Using a Novel Precursor in E. coli-Based Overexpression Systems

The importance of NMR spectroscopy in unraveling the structural and dynamic properties of proteins is ever-expanding owing to progress in experimental techniques, hardware development, and novel labeling approaches. Multiple sophisticated methods of aliphatic residue labeling can be found in the literature, whereas the selective incorporation of NMR active isotopes into other amino acids still holds the potential for improvement. In order to close this methodological gap, we present a novel metabolic precursor for cell-based protein overexpression to assemble ¹³ C/² H isotope patterns in the peptide backbone, as well as in side chain positions of a mechanistically distinguished histidine residue.

Exploring the mechanisms of action of human secretory RNase 3 and RNase 7 against Candida albicans

Human antimicrobial RNases, which belong to the vertebrate RNase A superfamily and are secreted upon infection, display a wide spectrum of antipathogen activities. In this work, we examined the antifungal activity of the eosinophil RNase 3 and the skin-derived RNase 7, two proteins expressed by innate cell types that are directly involved in the host defense against fungal infection. Candida albicans has been selected as a suitable working model for testing RNase activities toward a eukaryotic pathogen. We explored the distinct levels of action of both RNases on yeast by combining cell viability and membrane model assays together with protein labeling and confocal microscopy. Site-directed mutagenesis was applied to ablate either the protein active site or the key anchoring region for cell binding. This is the first integrated study that highlights the RNases' dual mechanism of action. Along with an overall membrane-destabilization process, the RNases could internalize and target cellular RNA. The data support the contribution of the enzymatic activity for the antipathogen action of both antimicrobial proteins, which can be envisaged as suitable templates for the development of novel antifungal drugs. We suggest that both human RNases work as multitasking antimicrobial proteins that provide a first line immune barrier.

The first crystal structure of human RNase 6 reveals a novel substrate-binding and cleavage site arrangement

We describe the first human RNase 6 crystal structure in complex with sulfate anions. Kinetic analysis, site-directed mutagenesis and molecular dynamics simulations identified novel substrate recognition and cleavage sites.

Human RNase 6 is a cationic secreted protein that belongs to the RNase A superfamily. Its expression is induced in neutrophils and monocytes upon bacterial infection, suggesting a role in host defence. We present here the crystal structure of RNase 6 obtained at 1.72 Å (1 Å=0.1 nm) resolution, which is the first report for the protein 3D structure and thereby setting the basis for functional studies. The structure shows an overall kidney-shaped globular fold shared with the other known family members. Three sulfate anions bound to RNase 6 were found, interacting with residues at the main active site (His¹⁵, His¹²² and Gln¹⁴) and cationic surface-exposed residues (His³⁶, His³⁹, Arg⁶⁶ and His⁶⁷). Kinetic characterization, together with prediction of protein–nucleotide complexes by molecular dynamics, was applied to analyse the RNase 6 substrate nitrogenous base and phosphate selectivity. Our results reveal that, although RNase 6 is a moderate catalyst in comparison with the pancreatic RNase type, its structure includes lineage-specific features that facilitate its activity towards polymeric nucleotide substrates. In particular, enzyme interactions at the substrate 5′ end can provide an endonuclease-type cleavage pattern. Interestingly, the RNase 6 crystal structure revealed a novel secondary active site conformed by the His³⁶–His³⁹ dyad that facilitates the polynucleotide substrate catalysis.

Integration of kinetic isotope effect analyses to elucidate ribonuclease mechanism

The well-studied mechanism of ribonuclease A is believed to involve concerted general acid-base catalysis by two histidine residues, His12 and His119. The basic features of this mechanism are often cited to explain rate enhancement by both protein and RNA enzymes that catalyze RNA 2'-O-transphosphorylation. Recent kinetic isotope effect analyses and computational studies are providing a more chemically detailed description of the mechanism of RNase A and the rate limiting transition state. Overall, the results support an asynchronous mechanism for both solution and ribonuclease catalyzed reactions in which breakdown of a transient dianoinic phosphorane intermediate by 5'OP bond cleavage is rate limiting. Relative to non-enzymatic reactions catalyzed by specific base, a smaller KIE on the 5'O leaving group and a less negative βLG are observed for RNase A catalysis. Quantum mechanical calculations consistent with these data support a model in which electrostatic and H-bonding interactions with the non-bridging oxygens and proton transfer from His119 render departure of the 5'O less advanced and stabilize charge buildup in the transition state. Both experiment and computation indicate advanced 2'OP bond formation in the rate limiting transition state. However, this feature makes it difficult to resolve the chemical steps involved in 2'O activation. Thus, modeling the transition state for RNase A catalysis underscores those elements of its chemical mechanism that are well resolved, as well as highlighting those where ambiguity remains. This article is part of a Special Issue entitled: Enzyme Transition States from Theory and Experiment.

Interpretation of pH-activity profiles for acid-base catalysis from molecular simulations

The measurement of reaction rate as a function of pH provides essential information about mechanism. These rates are sensitive to the pK(a) values of amino acids directly involved in catalysis that are often shifted by the enzyme active site environment. Experimentally observed pH-rate profiles are usually interpreted using simple kinetic models that allow estimation of "apparent pK(a)" values of presumed general acid and base catalysts. One of the underlying assumptions in these models is that the protonation states are uncorrelated. In this work, we introduce the use of constant pH molecular dynamics simulations in explicit solvent (CpHMD) with replica exchange in the pH-dimension (pH-REMD) as a tool to aid in the interpretation of pH-activity data of enzymes and to test the validity of different kinetic models. We apply the methods to RNase A, a prototype acid-base catalyst, to predict the macroscopic and microscopic pK(a) values, as well as the shape of the pH-rate profile. Results for apo and cCMP-bound RNase A agree well with available experimental data and suggest that deprotonation of the general acid and protonation of the general base are not strongly coupled in transphosphorylation and hydrolysis steps. Stronger coupling, however, is predicted for the Lys41 and His119 protonation states in apo RNase A, leading to the requirement for a microscopic kinetic model. This type of analysis may be important for other catalytic systems where the active forms of the implicated general acid and base are oppositely charged and more highly correlated. These results suggest a new way for CpHMD/pH-REMD simulations to bridge the gap with experiments to provide a molecular-level interpretation of pH-activity data in studies of enzyme mechanisms.

Evaluation of human pancreatic RNase as effector molecule in a therapeutic antibody platform

Human antibody-ribonuclease (RNase) fusion proteins, referred to as immunoRNases, have been proposed as an alternative to heterologous immunotoxins, without their immunogenicity and unspecific toxicity issues. In this study, we investigated if human pancreatic RNase will be suitable as effector component in a therapeutic antibody development platform. We generated several fusion proteins consisting of tumor-specific human immunoglobulins (IgGs) and human pancreatic RNase. Transient mammalian cell production was efficient and IgG-RNases were purified to homogeneity. Antigen binding was comparable to the parental antibodies and RNase catalytic activity was retained even in the presence of 50-fold molar excess of human cytosolic RNase inhibitor (RI). Serum stability, cell binding and internalization of IgG-RNases were comparable to the parental IgGs. Despite these promising properties, none of the IgG-RNases revealed significant inhibition of tumor cell growth in vitro even when targeting different antigens putatively employing different endocytotic pathways. The introduction of different linkers containing endosomal protease cleavage sites into the IgG-RNase did not enhance cytotoxicity. Similarly, RI evasive human pancreatic RNase variants mediated only small inhibiting effects on tumor cell growth at high concentrations, potentially reflecting inefficient cytosolic translocation. Taken together, human pancreatic RNase and variants did not prove to be generally suitable as effector component for a therapeutic antibody drug development platform.

Bovine pancreatic ribonuclease: fifty years of the first enzymatic reaction mechanism

Fifty years ago, the group of Tony Mathias and Bob Rabin at University College London deduced the first mechanism for catalysis by an enzyme, ribonuclease [Findlay, D., Herries, D. G., Mathias, A. P., Rabin, B. R., and Ross, C. A. (1961) Nature 190, 781-784]. Here, we celebrate this historic accomplishment by surveying knowledge of enzymology and protein science at that time, facts that led to the formulation of the mechanism, criticisms and alternative mechanisms, data that supported the proposed mechanism, and some of the refinements that have since provided a more precise picture of catalysis of RNA cleavage by ribonucleases. The Mathias and Rabin mechanism has appeared in numerous textbooks, monographs, and reviews and continues to have a profound impact on biochemistry.

Reduced protein adsorption by osmolytes

Osmolytes are substances that affect osmosis and are used by cells to adapt to environmental stress. Here, we report a neutron reflectivity study on the influence of some osmolytes on protein adsorption at solid-liquid interfaces. Bovine ribonuclease A (RNase) and bovine insulin were used as model proteins adsorbing at a hydrophilic silica and at a hydrophobic polystyrene surface. From the neutron reflectivity data, the adsorbed protein layers were characterized in terms of layer thickness, protein packing density, and adsorbed protein mass in the absence and presence of urea, trehalose, sucrose, and glycerol. All data point to the clear effect of these nonionic cosolvents on the degree of protein adsorption. For example, 1 M sucrose leads to a reduction of the adsorbed amount of RNase by 39% on a silica surface and by 71% on a polystyrene surface. Trehalose was found to exhibit activity similar to that of sucrose. The changes in adsorbed protein mass can be attributed to a decreased packing density of the proteins in the adsorbed layers. Moreover, we investigated insulin adsorption at a hydrophobic surface in the absence and presence of glycerol. The degree of insulin adsorption is decreased by even 80% in the presence of 4 M of glycerol. The results of this study demonstrate that nonionic cosolvents can be used to tune and control nonspecific protein adsorption at aqueous-solid interfaces, which might be relevant for biomedical applications.

On the information expressed in enzyme structure: more lessons from ribonuclease A

Brownian computations were directed at Ribonuclease A (RNase A) and variants in folded states so as to quantify information of the statistical type at the atom/covalent bond level. This advanced the research reported in this journal last year on the information properties of enzyme primary structure. Brownian computation data are illustrated for a sixteen-member library. The results identify signature traits that distinguish the folded wild type (WT) molecule from variants. The distinctions are explainable in terms of correlated information and dispersion energy. The Brownian tools used for this study can be directed at other protein families (e.g., kinases, isomerases, etc.) in rapid screening, QSAR, and design applications.

Functional and structural analyses of N-acylsulfonamide-linked dinucleoside inhibitors of RNase A

Molecular probes are useful for both studying and controlling the functions of enzymes and other proteins. The most useful probes have high affinity for their target, along with small size and resistance to degradation. Here, we report on new surrogates for nucleic acids that fulfill these criteria. Isosteres in which phosphoryl [R–O–P(O₂⁻)–O–R′] groups are replaced with N-acylsulfonamidyl [R–C(O)–N⁻–S(O₂)–R′] or sulfonimidyl [R–S(O₂)–N⁻–S(O₂)–R′] groups increase the number of nonbridging oxygens from two (phosphoryl) to three (N-acylsulfonamidyl) or four (sulfonimidyl). Six such isosteres were found to be more potent inhibitors of catalysis by bovine pancreatic RNase A than are parent compounds containing phosphoryl groups. The atomic structures of two RNase A·N-acylsulfonamide complexes were determined at high resolution by X-ray crystallography. The N-acylsulfonamidyl groups were observed to form more hydrogen bonds with active site residues than did the phosphoryl groups in analogous complexes. These data encourage the further development and use of N-acylsulfonamides and sulfonimides as antagonists of nucleic acid-binding proteins.

Ribonuclease S redux

The S-peptide and S-protein components of bovine pancreatic ribonuclease form a noncovalent complex with restored ribonucleolytic activity. Although this archetypal protein-fragment complementation system has been the object of historic work in protein chemistry, intrinsic limitations compromise its utility. Modern methods are shown to overcome those limitations and enable new applications.

Molecular dynamics simulation of bovine pancreatic ribonuclease A-CpA and transition state-like complexes

The mechanisms of enzymes are intimately connected with their overall structure and dynamics in solution. Experimentally, it is considerably challenging to provide detailed atomic level information about the conformational events that occur at different stages along the chemical reaction path. Here, theoretical tools may offer new potential insights that complement those obtained from experiments that may not yield an unambiguous mechanistic interpretation. In this study, we apply molecular dynamics simulations of bovine pancreatic ribonuclease A, an archetype ribonuclease, to study the conformational dynamics, structural relaxation, and differential solvation that occur at discrete stages of the transesterification and cleavage reaction. Simulations were performed with explicit solvation with rigorous electrostatics and utilize recently developed molecular mechanical force field parameters for transphosphorylation and hydrolysis transition state analogues. Herein, we present results for the enzyme complexed with the dinucleotide substrate cytidilyl-3',5'-adenosine (CpA) in the reactant, and transphosphorylation and hydrolysis transition states. A detailed analysis of active site structures and hydrogen-bond patterns is presented and compared. The integrity of the overall backbone structure is preserved in the simulations and supports a mechanism whereby His12 stabilizes accumulating negative charge at the transition states through hydrogen-bond donation to the nonbridge oxygens. Lys41 is shown to be highly versatile along the reaction coordinate and can aid in the stabilization of the dianionic transition state, while being poised to act as a general acid catalyst in the hydrolysis step.

Contributions of the histidine side chain and the N-terminal alpha-amino group to the binding thermodynamics of oligopeptides to nucleic acids as a function of pH

Interactions of histidine with nucleic acid phosphates and histidine pK(a) shifts make important contributions to many protein-nucleic acid binding processes. To characterize these phenomena in simplified systems, we quantified binding of a histidine-containing model peptide HWKK ((+)NH(3)-His-Trp-Lys-Lys-NH(2)) and its lysine analogue KWKK ((+)NH(3)-Lys-Trp-Lys-Lys-NH(2)) to a single-stranded RNA model, polyuridylate (polyU), by changes in tryptophan fluorescence as a function of salt concentration and pH. For both HWKK and KWKK, equilibrium binding constants, K(obs), and magnitudes of log-log salt derivatives, SK(obs) identical with (partial differential logK(obs)/partial differential log[Na(+)]), decreased with increasing pH in the manner expected for a titration curve model in which deprotonation of the histidine and alpha-amino groups weakens binding and reduces its salt-dependence. Fully protonated HWKK and KWKK exhibit the same K(obs) and SK(obs) within uncertainty, and these SK(obs) values are consistent with limiting-law polyelectrolyte theory for +4 cationic oligopeptides binding to single-stranded nucleic acids. The pH-dependence of HWKK binding to polyU provides no evidence for pK(a) shifts nor any requirement for histidine protonation, in stark contrast to the thermodynamics of coupled protonation often seen for these cationic residues in the context of native protein structure where histidine protonation satisfies specific interactions (e.g., salt-bridge formation) within highly complementary binding interfaces. The absence of pK(a) shifts in our studies indicates that additional Coulombic interactions across the nonspecific-binding interface between RNA and protonated histidine or the alpha-amino group are not sufficient to promote proton uptake for these oligopeptides. We present our findings in the context of hydration models for specific vs nonspecific nucleic acid binding.

X-ray crystallographic studies of RNase A variants engineered at the most destabilizing positions of the main hydrophobic core: further insight into protein stability

To investigate the structural origin of decreased pressure and temperature stability, the crystal structure of bovine pancreatic ribonuclease A variants V47A, V54A, V57A, I81A, I106A, and V108A was solved at 1.4-2.0 A resolution and compared with the structure of wild-type protein. The introduced mutations had only minor influence on the global structure of ribonuclease A. The structural changes had individual character that depends on the localization of mutated residue, however, they seemed to expand from mutation site to the rest of the structure. Several different parameters have been evaluated to find correlation with decrease of free energy of unfolding DeltaDeltaG(T), and the most significant correlation was found for main cavity volume change. Analysis of the difference distance matrices revealed that the ribonuclease A molecule is organized into five relatively rigid subdomains with individual response to mutation. This behavior consistent with results of unfolding experiments is an intrinsic feature of ribonuclease A that might be surviving remnants of folding intermediates and reflects the dynamic nature of the molecule.

On the information expressed in enzyme primary structure: lessons from Ribonuclease A

The information expressed in an enzyme's primary structure is investigated. Brownian computations are directed at Ribonuclease A (RNase A) so as to quantify the information at the atom/covalent bond level. The information content and distribution are crucial because the primary structure underpins the molecule's chemical functions. Brownian computation data are illustrated for the native protein, mutants, and sequence isomers. The results identify signature features of the active protein on new information grounds. The same tools offer rapid screening of proteins and polypeptides whereby several examples are illustrated.

Polymeric brushes as functional templates for immobilizing ribonuclease A: study of binding kinetics and activity

The ability to immobilize proteins with high binding capacities on surfaces while maintaining their activity is critical for protein microarrays and other biotechnological applications. We employed poly(acrylic acid) (PAA) brushes as templates to immobilize ribonuclease A (RNase A), which is commonly used to remove RNA from plasmid DNA preparations. The brushes are grown by surface-anchored atom-transfer radical polymerization (ATRP) initiators. RNase A was immobilized by both covalent esterification and a high binding capacity metal-ion complexation method to PAA brushes. The polymer brushes immobilized 30 times more enzyme compared to self-assembled monolayers. As the thickness of the brush increases, the surface density of the RNase A increases monotonically. The immobilization was investigated by ellipsometry, X-ray photoelectron spectroscopy (XPS), thermogravimetric analysis (TGA), and near-edge X-ray absorption fine structure spectroscopy (NEXAFS). The activity of the immobilized RNase A was determined using UV absorbance. As much as 11.0 microg/cm(2) of RNase A was bound to PAA brushes by metal-ion complexation compared to 5.8 microg/cm(2) by covalent immobilization which is 30 and 16 times the estimated mass bound in a monolayer. The calculated diffusion coefficient D was 0.63 x 10(-14) cm(2)/s for metal-ion complexation and 0.71 x 10(-14) cm(2)/s for covalent immobilization. Similar values of D indicate that the binding kinetics is similar, but the thermodynamic equilibrium coverage varies with the binding chemistry. Immobilization kinetics and thermodynamics were characterized by ellipsometry for both methods. A maximum relative activity of 0.70-0.80 was reached between five and nine monolayers of the immobilized enzyme. However, the relative activity for covalent immobilization was greater than that of metal-ion complexation. Covalent esterification resulted in similar temperature dependence as free enzyme, whereas metal-ion complexation showed no temperature dependence indicating a significant change in conformation.

Structural basis for catalysis by onconase

Onconase (ONC) is a homolog of bovine pancreatic ribonuclease (RNase A) from the frog Rana pipiens. ONC displays antitumoral activity and is in advanced clinical trials for the treatment of cancer. Here, we report the first atomic structures of ONC-nucleic acid complexes: a T89N/E91A ONC-5'-AMP complex at 1.65 A resolution and a wild-type ONC-d(AUGA) complex at 1.90 A resolution. The latter structure and site-directed mutagenesis were used to reveal the atomic basis for substrate recognition and turnover by ONC. The residues in ONC that are proximal to the scissile phosphodiester bond (His10, Lys31, and His97) and uracil nucleobase (Thr35, Asp67, and Phe98) are conserved from RNase A and serve to generate a similar bell-shaped pH versus k(cat)/K(M) profile for RNA cleavage. Glu91 of ONC forms two hydrogen bonds with the guanine nucleobase in d(AUGA), and Thr89 is in close proximity to that nucleobase. Installing a neutral or cationic residue at position 91 or an asparagine residue at position 89 virtually eliminated the 10(2)-fold guanine:adenine preference of ONC. A variant that combined such substitutions, T89N/E91A ONC, actually preferred adenine over guanine. In contrast, installing an arginine residue at position 91 increased the guanine preference and afforded an ONC variant with the highest known k(cat)/K(M) value. These data indicate that ONC discriminates between guanine and adenine by using Coulombic interactions and a network of hydrogen bonds. The structure of the ONC-d(AUGA) complex was also used to probe other aspects of catalysis. For example, the T5R substitution, designed to create a favorable Coulombic interaction between ONC and a phosphoryl group in RNA, increased ribonucleolytic activity by twofold. No variant, however, was more toxic to human cancer cells than wild-type ONC. Together, these findings provide a cynosure for understanding catalysis of RNA cleavage in a system of high medicinal relevance.

Recognition of ribonuclease A by 3'-5'-pyrophosphate-linked dinucleotide inhibitors: a molecular dynamics/continuum electrostatics analysis

The proteins of the pancreatic ribonuclease A (RNase A) family catalyze the cleavage of the RNA polymer chain. The development of RNase inhibitors is of significant interest, as some of these compounds may have a therapeutic effect in pathological conditions associated with these proteins. The most potent low molecular weight inhibitor of RNase reported to date is the compound 5'-phospho-2'-deoxyuridine-3-pyrophosphate (P-->5)-adenosine-3-phosphate (pdUppA-3'-p). The 3',5'-pyrophosphate group of this compound increases its affinity and introduces structural features which seem to be unique in pyrophosphate-containing ligands bound to RNase A, such as the adoption of a syn conformation by the adenosine base at RNase subsite B(2) and the placement of the 5'-beta-phosphate of the adenylate (instead of the alpha-phosphate) at subsite P(1) where the phosphodiester bond cleavage occurs. In this work, we study by multi-ns molecular dynamics simulations the structural properties of RNase A complexes with the ligand pdUppA-3'-p and the related weaker inhibitor dUppA, which lacks the 3' and 5' terminal phosphate groups of pdUppA-3'-p. The simulations show that the adenylate 5'-beta-phosphate binding position and the adenosine syn orientation constitute robust structural features in both complexes, stabilized by persistent interactions with specific active-site residues of subsites P(1) and B(2). The simulation structures are used in conjunction with a continuum-electrostatics (Poisson-Boltzmann) model, to evaluate the relative binding affinity of the two complexes. The computed relative affinity of pdUppA-3'-p varies between -7.9 kcal/mol and -2.8 kcal/mol for a range of protein/ligand dielectric constants (epsilon(p)) 2-20, in good agreement with the experimental value (-3.6 kcal/mol); the agreement becomes exact with epsilon(p) = 8. The success of the continuum-electrostatics model suggests that the differences in affinity of the two ligands originate mainly from electrostatic interactions. A residue decomposition of the electrostatic free energies shows that the terminal phosphate groups of pdUppA-3'-p make increased interactions with residues Lys(7) and Lys(66) of the more remote sites P(2) and P(0), and His(119) of site P(1).

Substrate interactions with human ferrochelatase

Ferrochelatase, the terminal enzyme in heme biosynthesis, catalyzes the insertion of ferrous iron into protoporphyrin IX to form protoheme IX. Human ferrochelatase is a homodimeric, inner mitochondrial membrane-associated enzyme that possesses an essential [2Fe-2S] cluster. In this work, we report the crystal structure of human ferrochelatase with the substrate protoporphyrin IX bound as well as a higher resolution structure of the R115L variant without bound substrate. The data presented reveal that the porphyrin substrate is bound deep within an enclosed pocket. When compared with the location of N-methylmesoporphyrin in the Bacillus subtilis ferrochelatase, the porphyrin is rotated by approximately 100 degrees and is buried an additional 4.5 A deeper within the active site. The propionate groups of the substrate do not protrude into solvent and are bound in a manner similar to what has been observed in uroporphyrinogen decarboxylase. Furthermore, in the substrate-bound form, the jaws of the active site mouth are closed so that the porphyrin substrate is completely engulfed in the pocket. These data provide insights that will aid in the determination of the mechanism for ferrochelatase.

Genetic selection for critical residues in ribonucleases

Homologous mammalian proteins were subjected to an exhaustive search for residues that are critical to their structure/function. Error-prone polymerase chain reactions were used to generate random mutations in the genes of bovine pancreatic ribonuclease (RNase A) and human angiogenin, and a genetic selection based on the intrinsic cytotoxicity of ribonucleolytic activity was used to isolate inactive variants. Twenty-three of the 124 residues in RNase A were found to be intolerant to substitution with at least one particular amino acid. Twenty-nine of the 123 residues in angiogenin were likewise intolerant. In both RNase A and angiogenin, only six residues appeared to be wholly intolerant to substitution: two histidine residues involved in general acid/base catalysis and four cysteine residues that form two disulfide bonds. With few exceptions, the remaining critical residues were buried in the hydrophobic core of the proteins. Most of these residues were found to tolerate only conservative substitutions. The importance of a particular residue as revealed by this genetic selection correlated with its sequence conservation, though several non-conserved residues were found to be critical for protein structure/function. Despite voluminous research on RNase A, the importance of many residues identified herein was unknown, and those can now serve as targets for future work. Moreover, a comparison of the critical residues in RNase A and human angiogenin, which share only 35% amino acid sequence identity, provides a unique perspective on the molecular evolution of the RNase A superfamily, as well as an impetus for applying this methodology to other ribonucleases.

Kinetic characterization and ligand binding studies of His351 mutants of Bacillus anthracis adenylate cyclase

Edema factor is a calmodulin dependent adenylyl cyclase secreted as one of the primary exotoxins by Bacillus anthracis. A histidine residue at position 351 located in its active site has been implicated in catalysis but direct evidence of its functional role is still lacking. In the present study, we introduced mutations in full-length edema factor (EF) to generate alanine (H351A), asparagine (H351N), and phenylalanine (H351F) variants. Spectral analysis of these variants displayed no gross structural deformities. Kinetic characterization showed that the adenylyl cyclase activity of H351N and H351F mutants decreased 34- and 40-fold, respectively, whereas H351A mutant completely lost activity. K(m) and K(i) values for ATP, pH activity profiles, and calmodulin activation curves of asparagine and phenylalanine mutants were not altered markedly. This kinetic data corroborated our ligand binding studies. Apparent K(d) values for calmodulin and ATP binding were found to be similar for wild-type EF and these active site variants. The effective substitution of H351 by asparagine and phenylalanine, albeit at a greatly reduced K(cat), without perturbing the ATP binding highlights the importance of this residue in transition-state stabilization. This was also evident from the positive free energy difference calculated for these mutants. However, equilibrium dialysis experiments revealed noticeable increase in ATP binding constant of H351A mutant, suggesting an additional role of H351 in precise substrate binding in the catalytic pocket. This is the first comprehensive study that describes the kinetic and ligand binding properties of H351 mutants and validates the importance of this residue in EF catalysis.

Release factors eRF1 and RF2: a universal mechanism controls the large conformational changes

Class I release factors 1 and 2 (RF1 and RF2) terminate protein synthesis by recognizing stop codons on the mRNA via their conserved amino acid motifs (NIKS in eRF1 and SPF in RF2) and by the conserved tripeptide (GGQ) interactions with the ribosomal peptidyltransferase center. Crystal structures of eRF1 and RF2 do not fit their ribosomal binding pocket (approximately 73 angstroms). Cryoelectron microscopy indicates large conformational changes in the ribosome-bound RF2. Here, we investigate the conformational dynamics of the eRF1 and RF2 using molecular dynamics simulation, structural alignment, and electrostatic analysis of domain interactions. We show that relaxed eRF1 has a shape remarkably similar to the ribosome-bound RF2 observed by cryoelectron microscopy. The similarity between the two release factors is as good as between elongation factor G and elongation factor Tu-guanosine-5'(beta,gamma-imido)triphosphate-tRNA. Further, the conformational transitions and dynamics of eRF1 and RF2 between the free and ribosome-bound states are most likely controlled by protonation of conserved histidines. For eRF1, the distance between the NIKS and GGQ motifs shrinks from 97.5 angstroms in the crystal to 70-80 angstroms. For RF2, the separation between SPF and GGQ elongates from 32 angstroms in the crystal to 50 angstroms. Coulombic interaction strongly favors the open conformation of eRF1; however, solvation and histidine protonation modulate the domain interactions, making the closed conformation of eRF1 more accessible. Thus, RF1 and RF2 function like molecular machines, most likely fueled by histidine protonation. The unified conformational control and the shapes of eRF1 and RF2 support the proposition that the termination of protein synthesis involves similar mechanisms across species.

The unswapped chain of bovine seminal ribonuclease: Crystal structure of the free and liganded monomeric derivative

Bovine seminal ribonuclease, a homodimeric enzyme joined covalently by two interchain disulphide bonds, is an equilibrium mixture of two conformational isomers, MxM and M=M. The major form, MxM, whose crystal structure has been previously determined at 1.9 A resolution, presents the swapping of the N-terminal segments (residues 1-15) and composite active sites formed by residues of different chains. The three-dimensional domain swapping does not occur in the M=M form. The different fold of each N-terminal tail is directed by the hinge loop (residue 16-22) connecting the swapping domain to the body of the protein. Reduction and alkylation of interchain disulphide bridges produce a monomeric derivative and a noncovalent swapped dimer, which are both active. The free and nucleotide-bound forms of the monomer have been crystallized at an alkaline pH and refined at 1.45 and 1.65 A resolution, respectively. In both cases, the N-terminal fragment is folded on the main body of the protein to produce an intact active site and a chain architecture very similar to that of bovine pancreatic ribonuclease. In this new fold of the seminal chain, the hinge loop is disordered. Despite the difference between the tertiary structure of the monomer and that of the chains in the MxM form, the active sites of the two enzymes are virtually indistinguishable. Furthermore, the structure of the liganded enzyme represents the first example of a ribonuclease complex studied at an alkaline pH and provides new information on the binding of a nucleotide when the catalytic histidines are deprotonated.

Adjacent cysteine residues as a redox switch

Oxidation of adjacent cysteine residues into a cystine forms a strained eight-membered ring. This motif was tested as the basis for an enzyme with an artificial redox switch. Adjacent cysteine residues were introduced into two different structural contexts in ribonuclease A (RNase A) by site-directed mutagenesis to produce the A5C/A6C and S15C/S16C variants. Ala5 and Ala6 are located in an alpha-helix, whereas Ser15 and Ser16 are located in a surface loop. Only A5C/A6C RNase A had the desired property. The catalytic activity of this variant decreases by 70% upon oxidation. The new disulfide bond also decreases the conformational stability of the A5C/A6C variant. Reduction with dithiothreitol restores full enzymatic activity. Thus, the insertion of adjacent cysteine residues in a proper context can be used to modulate enzymatic activity.