High-resolution three-dimensional structure of ribonuclease A in solution by nuclear magnetic resonance spectroscopy

Conserved structural topologies in RNase A-like and trypsin-like serine proteases: a sequence-based folding analysis

BACKGROUND

Protein folding is a complex process in which amino acid sequences encode the information required for a polypeptide chain to fold into its functional three-dimensional (3D) structure. Many proteins share common substructures and recurring secondary structure elements that contribute to similar 3D folding patterns, even across different protein families. This study examines two distinct groups of proteins, the RNase A-like fold and the trypsin-like serine protease fold, classified by SCOPe. These proteins share only some substructures that contribute to their folding cores. Despite minimal sequence identity, they exhibit partial structural similarities in their 3D topologies. We used a sequence-based approach, including inter-residue average distance statistics and contact frequency prediction, to explore these folding characteristics. Structural observations guided further analyses of conserved hydrophobic residue packing, highlighting key folding units within each fold.

RESULTS

Our analysis predicted two compact regions within each protein group. Interactions between these regions form a partially shared topology. We identified conserved hydrophobic residues critical to these interactions, suggesting a common mechanism for establishing these structural features. Despite overall structural differences between the RNase A-like and trypsin-like folds, our findings emphasize the presence of a shared partial folding core.

CONCLUSIONS

The partially shared structural features in the RNase A-like and trypsin-like serine protease folds reflect a convergent folding mechanism. This mechanism underscores the evolutionary adaptation of protein folding, where distinct folds can still retain critical, conserved structural motifs. These findings highlight how proteins with overall different topologies can evolve to share key folding features, demonstrating the elegance and efficiency of protein evolution.

On the use of dioxane as reference for determination of the hydrodynamic radius by NMR spectroscopy

Measuring the compaction of a protein or complex is key to our understanding of the interactions within and between biomolecules. Experimentally, protein compaction is often probed either by estimating the radius of gyration (Rg) obtained from small-angle x-ray scattering (SAXS) experiments or the hydrodynamic radius (Rh) obtained, for example, by pulsed field gradient NMR (PFG NMR) spectroscopy. PFG NMR experiments generally report on the translational diffusion coefficient, which in turn can be used to estimate Rh using an internal standard to account for sample viscosity and uncertainty about the gradient strength. 1,4-Dioxane is one such commonly used internal standard, and the reference value of Rh is therefore important. We have revisited the basis for the commonly used reference value for the Rh of dioxane (2.12 Å) that is used to convert measured diffusion coefficients into a hydrodynamic radius. We followed the same approach that was used to establish the current reference value by measuring SAXS and PFG NMR data for a set of seven different proteins and using these as standards. Our analysis shows that the current Rh reference value for dioxane Rh is underestimated, and we instead suggest a new value of 2.27 ± 0.04 Å. Using this updated reference value results in a ∼7% increase in Rh values for proteins whose hydrodynamic radii have been measured by PFG NMR. These results are particularly important when the absolute value of Rh is of interest such as when determining or validating ensemble descriptions of intrinsically disordered proteins.

Benchmarking predictive methods for small-angle X-ray scattering from atomic coordinates of proteins using maximum likelihood consensus data

Stimulated by informal conversations at the XVII International Small Angle Scattering (SAS) conference (Traverse City, 2017), an international team of experts undertook a round-robin exercise to produce a large dataset from proteins under standard solution conditions. These data were used to generate consensus SAS profiles for xylose isomerase, urate oxidase, xylanase, lysozyme and ribonuclease A. Here, we apply a new protocol using maximum likelihood with a larger number of the contributed datasets to generate improved consensus profiles. We investigate the fits of these profiles to predicted profiles from atomic coordinates that incorporate different models to account for the contribution to the scattering of water molecules of hydration surrounding proteins in solution. Programs using an implicit, shell-type hydration layer generally optimize fits to experimental data with the aid of two parameters that adjust the volume of the bulk solvent excluded by the protein and the contrast of the hydration layer. For these models, we found the error-weighted residual differences between the model and the experiment generally reflected the subsidiary maxima and minima in the consensus profiles that are determined by the size of the protein plus the hydration layer. By comparison, all-atom solute and solvent molecular dynamics (MD) simulations are without the benefit of adjustable parameters and, nonetheless, they yielded at least equally good fits with residual differences that are less reflective of the structure in the consensus profile. Further, where MD simulations accounted for the precise solvent composition of the experiment, specifically the inclusion of ions, the modelled radius of gyration values were significantly closer to the experiment. The power of adjustable parameters to mask real differences between a model and the structure present in solution is demonstrated by the results for the conformationally dynamic ribonuclease A and calculations with pseudo-experimental data. This study shows that, while methods invoking an implicit hydration layer have the unequivocal advantage of speed, care is needed to understand the influence of the adjustable parameters. All-atom solute and solvent MD simulations are slower but are less susceptible to false positives, and can account for thermal fluctuations in atomic positions, and more accurately represent the water molecules of hydration that contribute to the scattering profile.

A round-robin approach provides a detailed assessment of biomolecular small-angle scattering data reproducibility and yields consensus curves for benchmarking

Through an expansive international effort that involved data collection on 12 small-angle X-ray scattering (SAXS) and four small-angle neutron scattering (SANS) instruments, 171 SAXS and 76 SANS measurements for five proteins (ribonuclease A, lysozyme, xylanase, urate oxidase and xylose isomerase) were acquired. From these data, the solvent-subtracted protein scattering profiles were shown to be reproducible, with the caveat that an additive constant adjustment was required to account for small errors in solvent subtraction. Further, the major features of the obtained consensus SAXS data over the q measurement range 0-1 Å-1 are consistent with theoretical prediction. The inherently lower statistical precision for SANS limited the reliably measured q-range to <0.5 Å-1, but within the limits of experimental uncertainties the major features of the consensus SANS data were also consistent with prediction for all five proteins measured in H2O and in D2O. Thus, a foundation set of consensus SAS profiles has been obtained for benchmarking scattering-profile prediction from atomic coordinates. Additionally, two sets of SAXS data measured at different facilities to q > 2.2 Å-1 showed good mutual agreement, affirming that this region has interpretable features for structural modelling. SAS measurements with inline size-exclusion chromatography (SEC) proved to be generally superior for eliminating sample heterogeneity, but with unavoidable sample dilution during column elution, while batch SAS data collected at higher concentrations and for longer times provided superior statistical precision. Careful merging of data measured using inline SEC and batch modes, or low- and high-concentration data from batch measurements, was successful in eliminating small amounts of aggregate or interparticle interference from the scattering while providing improved statistical precision overall for the benchmarking data set.

Long term storage of miRNA at room and elevated temperatures in a silica sol-gel matrix

Storage of biospecimens in their near native environment at room temperature can have a transformative global impact, however, this remains an arduous challenge to date due to the rapid degradation of biospecimens over time. Currently, most isolated biospecimens are refrigerated for short-term storage and frozen (-20 °C, -80 °C, liquid nitrogen) for long-term storage. Recent advances in room temperature storage of purified biomolecules utilize anhydrobiosis. However, a near aqueous storage solution that can preserve the biospecimen nearly "as is" has not yet been achieved by any current technology. Here, we demonstrate an aqueous silica sol-gel matrix for optimized storage of biospecimens. Our technique is facile, reproducible, and has previously demonstrated stabilization of DNA and proteins, within a few minutes using a standard benchtop microwave. Herein, we demonstrate complete integrity of miRNA 21, a highly sensitive molecule at 4, 25, and 40 °C over a period of ∼3 months. In contrast, the control samples completely degrade in less than 1 week. We attribute excellent stability to entrapment of miRNA within silica-gel matrices.

Atomic-level reconstruction of biomolecules by a rigid-fragment- and local-frame-based (RF-LF) strategy

Coarse-grained (CG) model has been a powerful tool in bridging the gap between theoretical studies and experimental phenomena in biological computing field. The reconstruction from a CG model to an atomic-detail structure is especially important in CG studies of biological systems. In this work, a rigid-fragment- and local-frame-based (RF-LF) backmapping method was proposed to achieve reverse mapping from CG models to atomic-level structures. The initial atomic-level structures were further refined to yield the final backmapping ones. With the popular Martini force field, the performance of the RF-LF method was extensively examined in the CG → AA (CG to AA) backmapping of protein/DNA/RNA systems. Besides, the RF-LF method was also extended to the backmapping of the TMFF model. Numerical results illustrate that the RF-LF backmapping method is generic and parameter-free and can provide a promising way to tackle atomic-level studies in CG models.

Native mass spectrometry beyond ammonium acetate: effects of nonvolatile salts on protein stability and structure

Structures and stabilities of proteins investigated with native mass spectrometry can be affected by nonvolatile salts, including Tris buffer, in solution.

Structural and functional changes in RNAse A originating from tyrosine and histidine cross-linking and oxidation induced by singlet oxygen and peroxyl radicals

Oxidation can be induced by multiple processes in biological samples, with proteins being important targets due to their high abundance and reactivity. Oxidant reactions with proteins are not comprehensively understood, but it is known that structural and functional changes may be a cause, or a consequence, of disease. The mechanisms of oxidation of the model protein RNAse A by singlet oxygen (¹O₂) were examined and compared to peroxyl radical (ROO^•) oxidation, both common biological oxidants. This protein is a prototypic member of the RNAse family that exhibits antiviral activity by cleaving single-stranded RNA. RNAse A lacks tryptophan and cysteine residues which are major oxidant targets, but contains multiple histidine, tyrosine and methionine residues; these were therefore hypothesized to be the major sites of damage. ¹O₂ and ROO^• induce different patterns and extents of damage; both induce cross-links and side-chain oxidation, and ¹O₂ exposure modulates enzymatic activity. Multiple products have been characterized including methionine sulfoxide and sulfone, alcohols, DOPA, 2-oxohistidine, histidine-derived ring-opened species and inter- and intra-molecular cross-links (di-tyrosine, histidine-lysine, histidine-arginine, tyrosine-lysine). In addition to methionine modification, which appears not to be causative to activity loss, singlet oxygen also induces alteration to specific histidine, tyrosine and proline residues, including modification and cross-linking of the active site histidine, His12. The high homology among the RNAse family suggests that similar modifications may occur in humans, and be associated with the increased risk of viral infections in people with diabetes, as markers for ¹O₂ have been found in early stages of this pathology.

Two-bead polarizable water models combined with a two-bead multipole force field (TMFF) for coarse-grained simulation of proteins

(a) Four water molecules contained in the polarizable CG water models in (b) two-bead polarizable water model 1 (TPW1) and (c) two-bead polarizable water model 2 (TPW2).

Tautomeric stabilities of 4-fluorohistidine shed new light on mechanistic experiments with labeled ribonuclease A

Ribonuclease A is the oldest model for studying enzymatic mechanisms, yet questions remain about proton transfer within the active site. Seminal work by Jackson et al. (Science, 1994) labeled Ribonuclease A with 4-fluorohistidine, concluding that active-site histidines act as general acids and bases. Calculations of 4-fluorohistidine indicate that the π-tautomer is predominant in all simulated environments (by ~17 kJ/mol), strongly suggesting that fluoro-labeled ribonuclease A functions with His119 in π-tautomer. The tautomeric form of His119 during proton transfer and tautomerism as a putative mechanistic step in wild-type RNase A remain open questions and should be considered in future mechanistic studies.

A phthalimidation protocol that follows protein defined parameters

This work outlines the first phthalimidation protocol suitable for protein labeling and performed in aqueous media at room temperature and neutral pH with no catalyst or co-reagent required. The methodology is suitable for a range of amines and its efficiency was determined with chemoselective and site-selective protein labeling.

Targeting tumor micro-environment for design and development of novel anti-angiogenic agents arresting tumor growth

Angiogenesis: a process of generation of new blood vessels has been proved to be necessary for sustained tumor growth and cancer progression. Inhibiting angiogenesis pathway has long been remained a significant hope for the development of novel, effective and target orientated antitumor agents arresting the tumor proliferation and metastasis. The process of neoangiogenesis as a biological process is regulated by several pro- and anti-angiogenic factors, especially vascular endothelial growth factor, fibroblast growth factor, epidermal growth factor, hypoxia inducible factor 1 and transforming growth factor. Every endothelial cell destined for vessel formation is equipped with receptors for these angiogenic peptides. Moreover, numerous other angiogenic cytokines such as platelet derived growth factor (PGDF), placenta growth factor (PGF), nerve growth factor (NGF), stem-cell factor (SCF), and interleukins-2, 4, 6 etc. These molecular players performs critical role in regulating the angiogenic switch. Couple of decade's research in molecular aspects of tumor biology has unraveled numerous structural and functional mysteries of these angiogenic peptides. In present article, a detailed update on the functional and structural peculiarities of the various angiogenic peptides is described focusing on structural opportunities made available that has potential to be used to modulate function of these angiogenic peptides in developing therapeutic agents targeting neoplastic angiogenesis. The data may be useful in the mainstream of developing novel anticancer agents targeting tumor angiogenesis. We also discuss major therapeutic agents that are currently used in angiogenesis associated therapies as well as those are subject of active research or are in clinical trials.

All-atom ensemble modeling to analyze small-angle x-ray scattering of glycosylated proteins

The flexible and heterogeneous nature of carbohydrate chains often renders glycoproteins refractory to traditional structure determination methods. Small-angle X-ray scattering (SAXS) can be a useful tool for obtaining structural information of these systems. All-atom modeling of glycoproteins with flexible glycan chains was applied to interpret the solution SAXS data for a set of glycoproteins. For simpler systems (single glycan, with a well-defined protein structure), all-atom modeling generates models in excellent agreement with the scattering pattern and reveals the approximate spatial occupancy of the glycan chain in solution. For more complex systems (several glycan chains, or unknown protein substructure), the approach can still provide insightful models, though the orientations of glycans become poorly determined. Ab initio shape reconstructions appear to capture the global morphology of glycoproteins but in most cases offer little information about glycan spatial occupancy. The all-atom modeling methodology is available as a web server at http://salilab.org/allosmod-foxs.

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.

Crystal structure of Onconase at 1.1 Å resolution--insights into substrate binding and collective motion

Onconase® (ONC) is an amphibian member of the pancreatic ribonuclease superfamily that is selectively toxic to tumor cells. It is a much less efficient enzyme than the archetypal ribonuclease A and, in an attempt to gain further insight, we report the first atomic resolution crystal structure of ONC, determined in complex with sulfate ions at 100 K. The electron density map is of a quality sufficient to reveal significant nonplanarity in several peptide bonds. The majority of active site residues are very well defined, with the exceptions being Lys31 from the catalytic triad and Lys33 from the B₁ subsite, which are relatively mobile but rigidify upon nucleotide binding. Cryocooling causes a compaction of the unit cell and the protein contained within. This is principally the result of an inward movement of one of the lobes of the enzyme (lobe 2), which also narrows the active site cleft. Binding a nucleotide in place of sulfate is associated with an approximately perpendicular movement of lobe 2 and has little further effect on the cleft width. Aspects of this deformation are present in the principal axes of anisotropy extracted from C^α atomic displacement parameters, indicating its intrinsic nature. The three lowest-frequency modes of ONC motion predicted by an anisotropic network model are compaction/expansion variations in which lobe 2 is the prime mover. Two of these have high similarity to the cryocooling response and imply that the essential ‘breathing’ motion of ribonuclease A is conserved in ONC. Instead, shifts in conformational equilibria may contribute to the reduced ribonucleolytic activity of ONC.

Comparative NMR analysis of an 80-residue G protein-coupled receptor fragment in two membrane mimetic environments

Fragments of integral membrane proteins have been used to study the physical chemical properties of regions of transporters and receptors. Ste2p(G31-T110) is an 80-residue polypeptide which contains a portion of the N-terminal domain, transmembrane domain 1 (TM1), intracellular loop 1, TM2 and part of extracellular loop 1 of the α-factor receptor (Ste2p) from Saccharomyces cerevisiae. The structure of this peptide was previously determined to form a helical hairpin in lyso-palmitoylphosphatidyl-glycerol micelles (LPPG) [1]. Herein, we perform a systematic comparison of the structure of this protein fragment in micelles and trifluoroethanol (TFE):water in order to understand whether spectra recorded in organic:aqueous medium can facilitate the structure determination in a micellar environment. Using uniformly labeled peptide and peptide selectively protonated on Ile, Val and Leu methyl groups in a perdeuterated background and a broad set of 3D NMR experiments we assigned 89% of the observable atoms. NOEs and chemical shift analysis were used to define the helical regions of the fragment. Together with constraints from paramagnetic spin labeling, NOEs were used to calculate a transiently folded helical hairpin structure for this peptide in TFE:water. Correlation of chemical shifts was insufficient to transfer assignments from TFE:water to LPPG spectra in the absence of further information.

The crystal structure of ribonuclease A in complex with thymidine-3'-monophosphate provides further insight into ligand binding

Thymidine-3'-monophosphate (3'-TMP) is a competitive inhibitor analogue of the 3'-CMP and 3'-UMP natural product inhibitors of bovine pancreatic ribonuclease A (RNase A). Isothermal titration calorimetry experiments show that 3'-TMP binds the enzyme with a dissociation constant (K(d)) of 15 microM making it one of the strongest binding members of the five natural bases found in nucleic acids (A, C, G, T, and U). To further investigate the molecular properties of this potent natural affinity, we have determined the crystal structure of bovine pancreatic RNase A in complex with 3'-TMP at 1.55 A resolution and we have performed NMR binding experiments with 3'-CMP and 3'-TMP. Our results show that binding of 3'-TMP is very similar to other natural and non-natural pyrimidine ligands, demonstrating that single nucleotide affinity is independent of the presence or absence of a 2'-hydroxyl on the ribose moiety of pyrimidines and suggesting that the pyrimidine binding subsite of RNase A is not a significant contributor of inhibitor discrimination. Accumulating evidence suggests that very subtle structural, chemical, and potentially motional variations contribute to ligand discrimination in this enzyme.

NMR structural determinants of eosinophil cationic protein binding to membrane and heparin mimetics

Eosinophil cationic protein (ECP) is a highly stable, cytotoxic ribonuclease with the ability to enter and disrupt membranes that participates in innate immune defense against parasites but also kills human cells. We have used NMR spectroscopy to characterize the binding of ECP to membrane and heparin mimetics at a residue level. We believe we have identified three Arg-rich surface loops and Trp(35) as crucial for membrane binding. Importantly, we have provided evidence that the interaction surface of ECP with heparin mimetics is extended with respect to that previously described (fragment 34-38). We believe we have identified new sites involved in the interaction for the first time, and shown that the N-terminal alpha-helix, the third loop, and the first and last beta-strands are key for heparin binding. We have also shown that a biologically active ECP N-terminal fragment comprising the first 45 residues (ECP1-45) retains the capacity to bind membrane and heparin mimetics, thus neither the ECP tertiary structure nor its high conformational stability are required for cytotoxicity.

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.

Mapping the stability clusters in bovine pancreatic ribonuclease A

In the present work, we have thermodynamically characterized the thermally induced unfolding of 20 variants of bovine pancreatic ribonuclease A (RNase A) to experimentally describe the residues and the regions that are critical for the stability of the enzyme. The achieved results, complemented with previous studies by our group, allowed us to define the significance of the two hydrophobic nuclei present in the RNase A structure, as well as the contribution of the participating residues within each nucleus, to the global enzyme stability. We propose a structural model for the major and the minor hydrophobic nuclei of RNase A. The major nucleus is composite and located in the cavity delimited by alpha-helices 1 and 3, and the beta-sheet that is formed by strands 2, 3, 5, and 6. It consists of a central tight packed part constituted by residues Phe8, Met13, Val54, Val57, Ile106, Val108, and Phe120. This central part is surrounded by a layer formed by residues Val63, Tyr73, Met79, Ile107, Val116, and Val118. The minor nucleus, although less complex, is also constituted by a tight packing that involves the side chains of residues Tyr25, Met29, Met30, Leu35, Phe46, and Tyr97, which fill the cavity that originates the beta-sheet formed by beta-strands 1, 4, and 5 together with alpha-helix2.

The (1)H, (13)C, (15)N resonance assignment, solution structure, and residue level stability of eosinophil cationic protein/RNase 3 determined by NMR spectroscopy

Eosinophil cationic protein (ECP)/human RNase 3, a member of the RNase A family, is a remarkably cytotoxic protein implicated in asthma and allergies. These activities are probably due to ECP's ability to interact with and disrupt membranes and depend on two Trp, 19 Arg, and possibly an extremely high conformational stability. Here, we have used NMR spectroscopy to assign essentially all (1)H, (15)N, and backbone (13)C resonances, to solve the 3D structure in aqueous solution and to quantify its residue-level stability. The NMR solution structure was determined on the basis of 2316 distance constraints and is well-defined (backbone RMSD = 0.81 A). The N-terminus and the loop composed of residues 114-123 are relatively well-ordered; in contrast, conformational diversity is observed for the loop segments 17-22, 65-68, and 92-95 and most exposed sidechains. The side chain NH groups of the two Trp and 19 Arg showed no significant protection against hydrogen/deuterium exchange. The most protected NH groups belong to the first and last two beta-strands, and curiously, the first alpha-helix. Analysis of their exchange rates reveals a strikingly high global stability of 11.8 kcal/mol. This value and other stability measurements are used to better quantify ECP's unfolding thermodynamics.

Multiple solvent crystal structures of ribonuclease A: an assessment of the method

The multiple solvent crystal structures (MSCS) method uses organic solvents to map the surfaces of proteins. It identifies binding sites and allows for a more thorough examination of protein plasticity and hydration than could be achieved by a single structure. The crystal structures of bovine pancreatic ribonuclease A (RNAse A) soaked in the following organic solvents are presented: 50% dioxane, 50% dimethylformamide, 70% dimethylsulfoxide, 70% 1,6-hexanediol, 70% isopropanol, 50% R,S,R-bisfuran alcohol, 70% t-butanol, 50% trifluoroethanol, or 1.0M trimethylamine-N-oxide. This set of structures is compared with four sets of crystal structures of RNAse A from the protein data bank (PDB) and with the solution NMR structure to assess the validity of previously untested assumptions associated with MSCS analysis. Plasticity from MSCS is the same as from PDB structures obtained in the same crystal form and deviates only at crystal contacts when compared to structures from a diverse set of crystal environments. Furthermore, there is a good correlation between plasticity as observed by MSCS and the dynamic regions seen by NMR. Conserved water binding sites are identified by MSCS to be those that are conserved in the sets of structures taken from the PDB. Comparison of the MSCS structures with inhibitor-bound crystal structures of RNAse A reveals that the organic solvent molecules identify key interactions made by inhibitor molecules, highlighting ligand binding hot-spots in the active site. The present work firmly establishes the relevance of information obtained by MSCS.

Refolding of ribonuclease A monitored by real-time photo-CIDNP NMR spectroscopy

Photo-CIDNP NMR spectroscopy is a powerful method for investigating the solvent accessibility of histidine, tyrosine and tryptophan residues in a protein. When coupled to real-time NMR, this technique allows changes in the environments of these residues to be used as a probe of protein folding. In this paper we describe experiments performed to monitor the refolding of ribonuclease A following dilution from a high concentration of chemical denaturant. These experiments provide a good example of the utility of this technique which provides information that is difficult to obtain by other biophysical methods. Real-time photo-CIDNP measurements yield residue-specific kinetic data pertaining to the folding reaction, interpreted in terms of current knowledge of the folding of bovine pancreatic ribonuclease A.

Disulfide-linked protein folding pathways

Determining the mechanism by which proteins attain their native structure is an important but difficult problem in basic biology. The study of protein folding is difficult because it involves the identification and characterization of folding intermediates that are only very transiently present. Disulfide bond formation is thermodynamically linked to protein folding. The availability of thiol trapping reagents and the relatively slow kinetics of disulfide bond formation have facilitated the isolation, purification, and characterization of disulfide-linked folding intermediates. As a result, the folding pathways of several disulfide-rich proteins are among the best known of any protein. This review discusses disulfide bond formation and its relationship to protein folding in vitro and in vivo.

The structural basis of the difference in sensitivity for PNGase F in the de-N-glycosylation of the native bovine pancreatic ribonucleases B and BS

In glycoanalysis protocols, N-glycans from glycoproteins are most frequently released with peptide- N (4)-( N-acetyl-beta-glucosaminyl)asparagine amidase F (PNGase F). As the enzyme is an amidase, it cleaves the NH-CO linkage between the Asn side chain and the Asn-bound GlcNAc residue. Usually, the enzyme has a low activity, or is not active at all, on native glycoproteins. A typical example is native bovine pancreatic ribonuclease B (RNase B) with oligomannose-type N-glycans at Asn-34. However, native RNase BS, generated by subtilisin digestion of native RNase B, which comprises amino acid residues 21-124 of RNase B, is sensitive to PNGase F digestion. The same holds for carboxymethylated RNase B (RNase B (cm)). In this study, NMR spectroscopy and molecular modeling have been used to explain the differences in PNGase F activity for native RNase B, native RNase BS, and RNase B (cm). NMR analysis combined with literature data clearly indicated that the N-glycan at Asn-34 is more mobile in RNase BS than in RNase B. MD simulations showed that the region around Asn-34 in RNase B is not very flexible, whereby the alpha-helix of the amino acid residues 1-20 has a stabilizing effect. In RNase BS, the alpha-helix formed by amino acid residues 23-32 is significantly more flexible. Using these data, the possibilities for complex formation of both RNase B and RNase BS with PNGase F were studied, and a model for the RNase BS-PNGase F complex is proposed.

Conformational changes below the Tm: molecular dynamics studies of the thermal pretransition of ribonuclease A

Recent work suggests that some native conformations of proteins can vary with temperature. To obtain an atomic-level description of this structural and conformational variation, we have performed all-atom, explicit-solvent molecular dynamics simulations of bovine pancreatic ribonuclease A (RNase A) up to its melting temperature (Tm approximately 337 K). RNase A has a thermal pretransition near 320 K [Stelea, S. D., Pancoska, P., Benight, A. S., and Keiderling, T. A. (2001) Protein Sci. 10, 970-978]. Our simulations identify a conformational change that coincides with this pretransition. Between 310 and 320 K, there is a small but significant decrease in the number of native contacts, beta-sheet hydrogen bonding, and deviation of backbone conformation from the starting structure, and an increase in the number of nonnative contacts. Native contacts are lost in beta-sheet regions and in alpha1, partially due to movement of alpha1 away from the beta-sheet core. At 330 and 340 K, a nonnative helical segment of residues 15-20 forms, corresponding to a helix observed in the N-terminal domain-swapped dimer [Liu, Y. S., Hart, P. J., Schulnegger, M. P., and Eisenberg, D. (1998) Proc. Natl. Acad. Sci. U.S.A. 95, 3437-3432]. The conformations observed at the higher temperatures possess nativelike topology and overall conformation, with many native contacts, but they have a disrupted active site. We propose that these conformations may represent the native state at elevated temperature, or the N' state. These simulations show that subtle, functionally important changes in protein conformation can occur below the Tm.

Destabilizing mutations alter the hydrogen exchange mechanism in ribonuclease A

The effect of strongly destabilizing mutations, I106A and V108G of Ribonuclease A (RNase A), on its structure and stability has been determined by NMR. The solution structures of these variants are essentially equivalent to RNase A. The exchange rates of the most protected amide protons in RNase A (35 degrees C), the I106A variant (35 degrees C), and the V108G variant (10 degrees C) yield stability values of 9.9, 6.0, and 6.8 kcal/mol, respectively, when analyzed assuming an EX2 exchange mechanism. Thus, the destabilization induced by these mutations is propagated throughout the protein. Simulation of RNase A hydrogen exchange indicates that the most protected protons in RNase A and the V108G variant exchange via the EX2 regime, whereas those of I106A exchange through a mixed EX1 + EX2 process. It is striking that a single point mutation can alter the overall exchange mechanism. Thus, destabilizing mutations joins high temperatures, high pH and the presence of denaturating agents as a factor that induces EX1 exchange in proteins. The calculations also indicate a shift from the EX2 to the EX1 mechanism for less protected groups within the same protein. This should be borne in mind when interpreting exchange data as a measure of local stability in less protected regions.

Oxygen accessibility to ribonuclease a: quantitative interpretation of nuclear spin relaxation induced by a freely diffusing paramagnet

The nuclear spin relaxation induced by a freely diffusing paramagnetic center provides a direct measure of intermolecular accessibility. A number of factors are involved in a quantitative interpretation of relaxation data including excluded volume effects, solvation differences, and the details of the electron spin relaxation in the paramagnetic center. In the case where the electron relaxation time is short compared with correlation times describing the electron-nuclear coupling, the nuclear spin relaxation rates may be related to the effective local concentration of the paramagnetic center at different locations about the solute of interest. The local concentrations may in turn be related to differences in the local free energies of interaction between the diffusing paramagnet and the cosolute. We demonstrate this approach for the case of ribonuclease A and deduce surface free energy differences for a large number of protein proton sites. We find that the oxygen accessibility is poorly represented by hard-sphere models such as computed solvent or steric accessibility. There is a distribution of local intermolecular interactions with a width of the order of RT that dominates the report of the intermolecular exploration of the protein by this simple solute.

Tracking local conformational changes of ribonuclease A using picosecond time-resolved fluorescence of the six tyrosine residues

The six tyrosine residues of ribonuclease A (RNase A) are used as individual intrinsic probes for tracking local conformational changes during unfolding. The fluorescence decays of RNase A are well described by sums of three exponentials with decay times (tau(1) = 1.7 ns, tau(2) = 180 ps, and tau(3) = 30 ps) and preexponential coefficients (A(1) = 1, A(2) = 1, and A(3) = 4) at pH 7, 25 degrees C. The decay times are controlled by photo-induced electron transfer from individual tyrosine residues to the nearest disulphide (-SS-), bridge, which is distance (R) dependent. We assign tau(1) to Tyr-76 (R = 12.8 A), tau(2) to Tyr-115 (R = 6.9 A), and tau(3) to Tyr-25, Tyr-73, Tyr-92, and Tyr-97 (all four at R = 5.5 +/- 0.3 A) at 23 degrees C. On the basis of this assignment, the results show that, upon thermal or chemical unfolding only Tyr-25, Tyr-92, and Tyr-76 undergo significant displacement from their nearest -SS- bridge. Despite reporting on different regions of the protein, the concordance between the transition temperatures, T(m), obtained from Tyr-76 (T(m) = 59.2 degrees C) and Tyr-25 and Tyr-92 (T(m) = 58.2 degrees C) suggests a single unfolding event in this temperature range that affects all these regions similarly.

Conformation of 3'CMP bound to RNase A using TrNOESY

The conditions for accurately determining distance constraints from TrNOESY data on a small ligand (3'CMP) bound to a small protein (RNase A, <14 kDa) are described. For small proteins, normal TrNOESY conditions of 10:1 ligand:protein or greater can lead to inaccurate structures for the ligand-bound conformation due to the contribution of the free ligand to the TrNOESY signals. By using two ligand:protein ratios (2:1 and 5:1), which give the same distance constraints, a conformation of 3'CMP bound to RNase A was determined (glycosidic torsion angle, chi=-166 degrees ; pseudorotational phase angle, 0 degrees < or = P < or =36 degrees ). Ligand-protein NOESY cross peaks were also observed and used to dock 3'CMP into the binding pocket of the apo-protein (7rsa). After energy minimization, the conformation of the 3'CMP:RNase A complex was similar to the X-ray structure (1 rpf) except that a C3'-endo conformation for the ribose ring (rather than C2'-exo conformation) was found in the TrNOESY structure.

Investigation of spectral diffusion in ribonuclease by photolabeling of intrinsic aromatic amino acids

Spectral diffusion dynamics in ribonuclease A was observed via the broadening of photochemical holes burned into the absorption spectrum of intrinsic tyrosine residues. Unlike previous results based on hole burning of chromophores in the pockets of heme proteins, where spectral diffusion develops according to a power law in time, the dynamics in ribonuclease follow a logarithmic law. The results suggest that the experiment preferentially labels the tyrosines located on the surface of the protein where the two-level system dynamics of the glass host matrix exert a strong influence.

Formation of methionine sulfoxide during glycoxidation and lipoxidation of ribonuclease A

Chemical modification of proteins by reactive oxygen species affects protein structure, function and turnover during aging and chronic disease. Some of this damage is direct, for example by oxidation of amino acids in protein by peroxide or other reactive oxygen species, but autoxidation of ambient carbohydrates and lipids amplifies both the oxidative and chemical damage to protein and leads to formation of advanced glycoxidation and lipoxidation end-products (AGE/ALEs). In previous work, we have observed the oxidation of methionine during glycoxidation and lipoxidation reactions, and in the present work we set out to determine if methionine sulfoxide (MetSO) in protein was a more sensitive indicator of glycoxidative and lipoxidative damage than AGE/ALEs. We also investigated the sites of methionine oxidation in a model protein, ribonuclease A (RNase), in order to determine whether analysis of the site specificity of methionine oxidation in proteins could be used to indicate the source of the oxidative damage, i.e. carbohydrate or lipid. We describe here the development of an LC/MS/MS for quantification of methionine oxidation at specific sites in RNase during glycoxidation or lipoxidation by glucose or arachidonate, respectively. Glycoxidized and lipoxidized RNase were analyzed by tryptic digestion, followed by reversed phase HPLC and mass spectrometric analysis to quantify methionine and methionine sulfoxide containing peptides. We observed that: (1) compared to AGE/ALEs, methionine sulfoxide was a more sensitive biomarker of glycoxidative or lipoxidative damage to proteins; (2) regardless of oxidizable substrate, the relative rate of oxidation of methionine residues in RNase was Met29>Met30>Met13, with Met79 being resistant to oxidation; and (3) arachidonate produced a significantly greater yield of MetSO, compared to glucose. The methods developed here should be useful for assessing a protein's overall exposure to oxidative stress from a variety of sources in vivo.

Increase of RNase a N-terminus polarity or C-terminus apolarity changes the two domains' propensity to swap and form the two dimeric conformers of the protein

Do the polarities of the N-terminus or the apolarity of the C-terminus of bovine RNase A influence the relative yields of its two 3D domain-swapped dimeric conformers, the N-dimer and C-dimer? We have addressed this question by substituting Ala-4 or Ala-5 with serine (A4S and A5S mutants) or Ser-123 with alanine (S123A mutant) through site-directed mutagenesis. Both the polarity of the N-terminus and the apolarity of the C-terminus of RNase A were, therefore, increased. CD spectra revealed no significant differences between the secondary structures of the mutants and native RNase A. According to thermal denaturation analyses, the A4S and A5S mutants are less stable, and the S123A mutant is more stable than wild type RNase A. By subjecting the mutants under mild or drastic denaturing conditions, side-by-side with native and recombinant RNase A, to a thermally induced oligomerization procedure, the following results were obtained. (i) The N-terminal mutants showed a higher propensity, with respect to the native protein, to form N-dimers under mild unfolding conditions. (ii) The C-terminal mutant showed a higher propensity to form the C-dimer under severely unfolding conditions. These results are discussed in light of the relative stabilities of the various RNase A species under different environmental conditions, and we conclude that the hydrophilic or hydrophobic character of the RNase N-terminus or C-terminus can be an important variable governing the oligomerization of RNase A and possibly other proteins through the 3D domain-swapping mechanism.

Correlation of folding kinetics with the number and isomerization states of prolines in three homologous proteins of the RNase family

Several studies attribute the slower phases in protein folding to prolyl isomerizations, and several others do not. A correlation exists between the number of prolines in a protein and the complexity of the mechanism with which it folds. In this study, we have demonstrated a direct correlation between the number of cis-prolyl bonds in a native protein and the complexity with which it folds via slower phases by studying the folding of three structurally homologous proteins of the ribonuclease family, namely RNase A, onconase and angiogenin, which differ in the number and isomerization states of their proline residues.

Pressure-jump-induced kinetics reveals a hydration dependent folding/unfolding mechanism of ribonuclease A

Pressure-jump (p-jump)-induced relaxation kinetics was used to explore the energy landscape of protein folding/unfolding of Y115W, a fluorescent variant of ribonuclease A. Pressure-jumps of 40 MPa amplitude (5 ms dead-time) were conducted both to higher (unfolding) and to lower (folding) pressure, in the range from 100 to 500 MPa, between 30 and 50 degrees C. Significant deviations from the expected symmetrical protein relaxation kinetics were observed. Whereas downward p-jumps resulted always in single exponential kinetics, the kinetics induced by upward p-jumps were biphasic in the low pressure range and monophasic at higher pressures. The relative amplitude of the slow phase decreased as a function of both pressure and temperature. At 50 degrees C, only the fast phase remained. These results can be interpreted within the framework of a two-dimensional energy surface containing a pressure- and temperature-dependent barrier between two unfolded states differing in the isomeric state of the Asn-113-Pro-114 bond. Analysis of the activation volume of the fast kinetic phase revealed a temperature-dependent shift of the unfolding transition state to a larger volume. The observed compensation of this effect by glycerol offers an explanation for its protein stabilizing effect.

Chemically accurate protein structures: validation of protein NMR structures by comparison of measured and predicted pKa values

A new method is presented for evaluating the quality of protein structures obtained by NMR. This method exploits the dependence between measurable chemical properties of a protein, namely pKa values of acidic residues, and protein structure. The accurate and fast empirical computational method employed by the PROPKA program ( http://www.propka.chem.uiowa.edu) allows the user to test the ability of a given structure to reproduce known pKa values, which in turn can be used as a criterion for the selection of more accurate structures. We demonstrate the feasibility of this novel idea for a series of proteins for which both NMR and X-ray structures, as well as pKa values of all ionizable residues, have been determined. For the 17 NMR ensembles used in this study, this criterion is shown effective in the elimination of a large number of NMR structure ensemble members.

Enzyme dynamics along the reaction coordinate: critical role of a conserved residue

Conformational flexibility of the enzyme architecture is essential for biological function. These structural transitions often encompass significant portions of the enzyme molecule. Here, we present a detailed study of functionally relevant RNase A dynamics in the wild type and a D121A mutant form by NMR spin-relaxation techniques. In the wild-type enzyme, the dynamic properties are largely conserved in the apo, enzyme-substrate, and enzyme-product complexes. In comparison, mutation of aspartic acid 121 to alanine disrupts the timing of active-site dynamics, the product-release step, and global conformational changes, indicating that D121 plays a significant role in coordinating the dynamic events in RNase A. In addition, this mutation results in 90% loss of catalytic activity despite the absence of direct participation of D121 in the chemical reaction or in interactions with the substrate. These data suggest that one role of this conserved residue is to facilitate important millisecond protein dynamics.