Limits of NMR structure determination using variable target function calculations: ribonuclease T1, a case study

Improved NMR-data-compliant protein structure modeling captures context-dependent variations and expands the scope of functional inference

Nuclear magnetic resonance (NMR) spectroscopy can reveal conformational states of a protein in physiological conditions. However, sparsely available NMR data for a protein with large degrees of freedom can introduce structural artifacts in the built models. Currently used state-of-the-art methods deriving protein structure and conformation from NMR deploy molecular dynamics (MD) coupled with simulated annealing for building models. We provide an alternate graph-based modeling approach, where we first build substructures from NMR-derived distance-geometry constraints combined in one shot to form the core structure. The remaining molecule with inadequate data is modeled using a hybrid approach respecting the observed distance-geometry constraints. One-shot structure building is rarely undertaken for large and sparse data systems, but our data-driven bottom-up approach makes this uniquely feasible by suitable partitioning of the problem. A detailed comparison of select models with state-of-art methods reveals differences in the secondary structure regions wherein the correctness of our models is confirmed by NMR data. Benchmarking of 106 protein-folds covering 38-282 length structures shows minimal experimental-constraint violations while conforming to other structure quality parameters such as the proper folding, steric clash, and torsion angle violation based on Ramachandran plot criteria. Comparative MD studies using select protein models from a state-of-art method and ours under identical experimental parameters reveal distinct conformational dynamics that could be attributed to protein structure-function. Our work is thus useful in building enhanced NMR-evidence-based models that encapsulate the contextual secondary and tertiary structure variations present during the experimentation and expand the scope of functional inference.

Protein structure determination using a Riemannian approach

Protein NMR structure determination is one of the most extensively studied problems. Here, we adopt a novel method based on a matrix completion technique - the Riemannian approach - to rebuild the protein structure from the nuclear Overhauser effect distance restraints and the dihedral angle restraints. In comparison with the cyana method, the results generated via the Riemannian approach are more similar to the standard X-ray crystallographic structures as a result of the simple but powerful internal calculation processing function. In addition, our results demonstrate that the Riemannian approach has a comparable or even better performance than the cyana method on other structural assessment metrics, including the stereochemical quality and restraint violations. The Riemannian approach software is available at: https://github.com/xubiaopeng/Protein_Recon_MCRiemman.

Fabrication of RNA 3D Nanoprisms for Loading and Protection of Small RNAs and Model Drugs

Constructing containers with defined shape and size to load and protect therapeutics and subsequently control their release in the human body has long been a dream. The fabrication of 3D RNA prisms, characterized by atomic force microscopy, cryo-electron microscopy, dynamic light scattering, and polyacrylamide gel electrophoresis, is reported for the loading and protection of small molecules, proteins, small RNA molecules, and their controlled release.

CABS-flex predictions of protein flexibility compared with NMR ensembles

Motivation: Identification of flexible regions of protein structures is important for understanding of their biological functions. Recently, we have developed a fast approach for predicting protein structure fluctuations from a single protein model: the CABS-flex. CABS-flex was shown to be an efficient alternative to conventional all-atom molecular dynamics (MD). In this work, we evaluate CABS-flex and MD predictions by comparison with protein structural variations within NMR ensembles.

Results: Based on a benchmark set of 140 proteins, we show that the relative fluctuations of protein residues obtained from CABS-flex are well correlated to those of NMR ensembles. On average, this correlation is stronger than that between MD and NMR ensembles. In conclusion, CABS-flex is useful and complementary to MD in predicting protein regions that undergo conformational changes as well as the extent of such changes.

Availability and implementation: The CABS-flex is freely available to all users at http://biocomp.chem.uw.edu.pl/CABSflex.

Contact: sekmi@chem.uw.edu.pl

Supplementary information:Supplementary data are available at Bioinformatics online.

Extensive terminal and asymmetric processing of small RNAs from rRNAs, snoRNAs, snRNAs, and tRNAs

Deep sequencing studies frequently identify small RNA fragments of abundant RNAs. These fragments are thought to represent degradation products of their precursors. Using sequencing, computational analysis, and sensitive northern blot assays, we show that constitutively expressed non-coding RNAs such as tRNAs, snoRNAs, rRNAs and snRNAs preferentially produce small 5' and 3' end fragments. Similar to that of microRNA processing, these terminal fragments are generated in an asymmetric manner that predominantly favors either the 5' or 3' end. Terminal-specific and asymmetric processing of these small RNAs occurs in both mouse and human cells. In addition to the known processing of some 3' terminal tRNA-derived fragments (tRFs) by the RNase III endonuclease Dicer, we show that several RNase family members can produce tRFs, including Angiogenin that cleaves the TψC loop to generate 3' tRFs. The 3' terminal tRFs but not the 5' tRFs are highly complementary to human endogenous retroviral sequences in the genome. Despite their independence from Dicer processing, these tRFs associate with Ago2 and are capable of down regulating target genes by transcript cleavage in vitro. We suggest that endogenous 3' tRFs have a role in regulating the unwarranted expression of endogenous viruses through the RNA interference pathway.

Tryptophan conformations associated with partial unfolding in ribonuclease T1

The origin of the biexponential fluorescence decay of Trp in ribonuclease T1 under mildly destabilizing conditions, such as increased pH or temperature, or the presence of detergent, is still not understood. We have performed two extended replica-exchange molecular dynamics simulations to obtain a detailed representation of the native state at two protonation states corresponding to a high and low pH. At high pH, the appearance of partially unfolded states is evident. We found that this pH-induced destabilization originates from increased global repulsion as well as reduced local favorable electrostatic interactions and reduced H-bonding strength of His(27), His(40), and His(92). At high pH, alternative tryptophan rotamers appear and are linked to a distorted environment of the tryptophan, which also acts as a separate source of ground-state heterogeneity. The total population of these alternative conformations agrees well with the amplitude of the experimentally observed secondary fluorescence lifetime.

NMR and X-ray analysis of structural additivity in metal binding site-swapped hybrids of rubredoxin

Background

Chimeric hybrids derived from the rubredoxins of Pyrococcus furiosus (Pf) and Clostridium pasteurianum (Cp) provide a robust system for the characterization of protein conformational stability and dynamics in a differential mode. Interchange of the seven nonconserved residues of the metal binding site between the Pf and Cp rubredoxins yields a complementary pair of hybrids, for which the sum of the thermodynamic stabilities is equal to the sum for the parental proteins. Furthermore, the increase in amide hydrogen exchange rates for the hyperthermophile-derived metal binding site hybrid is faithfully mirrored by a corresponding decrease for the complementary hybrid that is derived from the less thermostable rubredoxin, indicating a degree of additivity in the conformational fluctuations that underlie these exchange reactions.

Results

Initial NMR studies indicated that the structures of the two complementary hybrids closely resemble "cut-and-paste" models derived from the parental Pf and Cp rubredoxins. This protein system offers a robust opportunity to characterize differences in solution structure, permitting the quantitative NMR chemical shift and NOE peak intensity data to be analyzed without recourse to the conventional conversion of experimental NOE peak intensities into distance restraints. The intensities for 1573 of the 1652 well-resolved NOE crosspeaks from the hybrid rubredoxins were statistically indistinguishable from the intensities of the corresponding parental crosspeaks, to within the baseplane noise level of these high sensitivity data sets. The differences in intensity for the remaining 79 NOE crosspeaks were directly ascribable to localized dynamical processes. Subsequent X-ray analysis of the metal binding site-swapped hybrids, to resolution limits of 0.79 Å and 1.04 Å, demonstrated that the backbone and sidechain heavy atoms in the NMR-derived structures lie within the range of structural variability exhibited among the individual molecules in the crystallographic asymmetric unit (~0.3 Å), indicating consistency with the "cut-and-paste" structuring of the hybrid rubredoxins in both crystal and solution.

Conclusion

Each of the significant energetic interactions in the metal binding site-swapped hybrids appears to exhibit a 1-to-1 correspondence with the interactions present in the corresponding parental rubredoxin structure, thus providing a structural basis for the observed additivity in conformational stability and dynamics. The congruence of these X-ray and NMR experimental data offers additional support for the interpretation that the conventional treatment of NOE distance restraints contributes substantially to the systematic differences that are commonly reported between NMR- and X-ray-derived protein structures.

Fungal ribotoxins: molecular dissection of a family of natural killers

RNase T1 is the best known representative of a large family of ribonucleolytic proteins secreted by fungi, mostly Aspergillus and Penicillium species. Ribotoxins stand out among them by their cytotoxic character. They exert their toxic action by first entering the cells and then cleaving a single phosphodiester bond located within a universally conserved sequence of the large rRNA gene, known as the sarcin-ricin loop. This cleavage leads to inhibition of protein biosynthesis, followed by cellular death by apoptosis. Although no protein receptor has been found for ribotoxins, they preferentially kill cells showing altered membrane permeability, such as those that are infected with virus or transformed. Many steps of the cytotoxic process have been elucidated at the molecular level by means of a variety of methodological approaches and the construction and purification of different mutant versions of these ribotoxins. Ribotoxins have been used for the construction of immunotoxins, because of their cytotoxicity. Besides this activity, Aspf1, a ribotoxin produced by Aspergillus fumigatus, has been shown to be one of the major allergens involved in allergic aspergillosis-related pathologies. Protein engineering and peptide synthesis have been used in order to understand the basis of these pathogenic mechanisms as well as to produce hypoallergenic proteins with potential diagnostic and immunotherapeutic applications.

Hydrogen Bonding Markedly Reduces the pK of Buried Carboxyl Groups in Proteins

The ionizable groups in proteins with the lowest pKs are the carboxyl groups of aspartic acid side-chains. One of the lowest, pK=0.6, is observed for Asp76 in ribonuclease T1. This low pK appeared to result from hydrogen bonds to a water molecule and to the side-chains of Asn9, Tyr11, and Thr91. The results here confirm this by showing that the pK of Asp76 increases to 1.7 in N9A, to 4.0 in Y11F, to 4.2 in T91V, to 4.4 in N9A+Y11F, to 4.9 in N9A+T91V, to 5.9 in Y11F+T91V, and to 6.4 in the triple mutant: N9A+Y11F+T91V. In ribonuclease Sa, the lowest pK=2.4 for Asp33. This pK increases to 3.9 in T56A, which removes the hydrogen bond to Asp33, and to 4.4 in T56V, which removes the hydrogen bond and replaces the -OH group with a -CH(3) group. It is clear that hydrogen bonds are able to markedly lower the pK values of carboxyl groups in proteins. These same hydrogen bonds make large contributions to the conformational stability of the proteins. At pH 7, the stability of D76A ribonuclease T1 is 3.8 kcal mol(-1) less than wild-type, and the stability of D33A ribonuclease Sa is 4.1 kcal mol(-1) less than wild-type. There is a good correlation between the changes in the pK values and the changes in stability. The results suggest that the pK values for these buried carboxyl groups would be greater than 8 in the absence of hydrogen bonds, and that the hydrogen bonds and other interactions of the carboxyl groups contribute over 8 kcal mol(-1) to the stability.

Conserved asparagine residue 54 of alpha-sarcin plays a role in protein stability and enzyme activity

Asparagine 54 of alpha-sarcin is a conserved residue within the proteins of the ribotoxin family of microbial ribonucleases. It is located in loop 2 of the protein, which lacks repetitive secondary structure elements but exhibits a well-defined conformation. Five mutant variants at this residue have been produced and characterized. The spectroscopic characterization of these proteins indicates that the overall conformation is not changed upon mutation. Activity and denaturation assays show that Asn-54 largely contributes to protein stability, and its presence is a requirement for the highly specific inhibitory activity of these ribotoxins on ribosomes.

NMR structure of the noncytotoxic alpha-sarcin mutant Delta(7-22): the importance of the native conformation of peripheral loops for activity

The deletion mutant Delta(7-22) of alpha-sarcin, unlike its wild-type protein counterpart, lacks the specific ability to degrade rRNA in intact ribosomes and exhibits an increased unspecific ribonuclease activity and decreased interaction with lipid vesicles. In trying to shed light on these differences, we report here on the three-dimensional structure of the Delta(7-22) alpha-sarcin mutant using NMR methods. We also evaluated its dynamic properties on the basis of theoretical models and measured its correlation time (6.2 nsec) by time-resolved fluorescence anisotropy. The global fold characteristic of ribotoxins is preserved in the mutant. The most significant differences with respect to the alpha-sarcin structure are concentrated in (1) loop 2, (2) loop 3, which adopts a new orientation, and (3) loop 5, which shows multiple conformations and an altered dynamics. The interactions between loop 5 and the N-terminal hairpin are lost in the mutant, producing increased solvent accessibility of the active-site residues. The degree of solvent exposure of the catalytic His 137 is similar to that shown by His 92 in RNase T1. Additionally, the calculated order parameters of residues belonging to loop 5 in the mutant correspond to an internal dynamic behavior more similar to RNase T1 than alpha-sarcin. On the other hand, changes in the relative orientation of loop 3 move the lysine-rich region 111-114, crucial for substrate recognition, away from the active site. All of the structural and dynamic data presented here reveal that the mutant is a hybrid of ribotoxins and noncytotoxic ribonucleases, consistent with its biological properties.

Thermal unfolding of ribonuclease T1 studied by multi-dimensional NMR spectroscopy

Thermal unfolding of ribonculease (RNase) T1 was studied by 1H nuclear Overhauser enhancement spectroscopy (NOESY) and 1H- 15N heteronuclear single-quantum coherence (HSQC) NMR spectroscopy at various temperatures. Native RNase T1 is a single-chain molecule of 104 amino acid residues, and has a single alpha-helix and two beta-sheets, A and B, which consist of two and five strands, respectively. Singular value decomposition analysis based on temperature-dependent HSQC spectra revealed that the thermal unfolding of RNase T1 can be described by a two-state transition model. The midpoint temperature and the change in enthalpy were determined as 54.0 degrees C and 696 kJ/mol, respectively, which are consistent with results obtained by other methods. To analyze the transition profile in more detail, we investigated local structural changes using temperature-dependent NOE intensities. The results indicate that the helical region starts to unfold at lower temperature than some beta-strands (B3, B4, and B5 in beta-sheet B). These beta-strands correspond to the hydrophobic cluster region, which had been expected to be a folding core. This was confirmed by structure calculations using the residual NOEs observed at 56 degrees C. Thus, the two-state transition of RNase T1 appears to involve locally different conformational changes.

Determination of the NMR structure of Gln25-ribonuclease T1

Ribonuclease (RNase) T1 is a guanyloribonuclease, having two isozymes in nature, Gln25- and Lys25-RNase T1. Between these two isozymes, there is no difference in catalytic activity and three-dimensional structure; however, Lys25-RNase T1 is slightly more stable than Gln25-RNase T1. Recently, it has been suggested that the existence of a salt bridge between Lys25 and Asp29/Glu31 in Lys25-RNase T1 contributes to the stability. To elucidate the effects of the replacement of Lys25 with a Gln on the conformation and microenvironments of RNase T1 in detail, the three-dimensional solution structure of Gln25-RNase T1 was determined by simulated-annealing calculations. As a result, the topology of the overall folding was shown to be very similar to that of the Lys25-isozyme except for some differences. In particular, there were two differences in the property of torsion angles of the two disulfide bonds and the conformations of the residues 11-13, 63-66, and 92-93. With regard to the residues 11-13, the lack of the above-mentioned salt bridge in Gln25-RNase T1 was thought to induce the conformational difference of this segment as compared with the Lys25-isozyme. Furthermore, it was proposed that the perturbation of this segment might transfer to the residues 92-93 via the two disulfide bonds.

Tautomeric state of alpha-sarcin histidines. Ndelta tautomers are a common feature in the active site of extracellular microbial ribonucleases

Extracellular fungal RNases, including ribotoxins such as alpha-sarcin, constitute a family of structurally related proteins represented by RNase T1. The tautomeric preferences of the alpha-sarcin imidazole side chains have been determined by nuclear magnetic resonance and electrostatic calculations. Histidine residues at the active site, H50 and H137, adopt the Ndelta tautomer, which is less common in short peptides, as has been found for RNase T1. Comparison with tautomers predicted from crystal structures of other ribonucleases suggests that two active site histidine residues with the Ndelta tautomer are a conserved feature of microbial ribonucleases and that this is related to their ribonucleolytic function.

Charge-charge interactions are the primary determinants of the pK values of the ionizable groups in Ribonuclease T1

Coulomb's law and a finite difference Poisson-Boltzmann based analysis are used to predict the pK values for 15 ionizable side chains (6 Asp, 6 Glu and 3 His) in ribonuclease T1. These predicted values are compared to the measured pK values to gain insight into the most important factors that influence the pK values of the ionizable groups in proteins. Charge-charge interactions are clearly the most important factor that determines the pK values of most ionizable groups in ribonuclease T1. However, pK values can be shifted by several pK units by the Born self energy associated with burying ionizable groups and by favorable intramolecular hydrogen bonding.

The structure and dipole moment of globular proteins in solution and crystalline states: use of NMR and X-ray databases for the numerical calculation of dipole moment

The large dipole moment of globular proteins has been well known because of the detailed studies using dielectric relaxation and electro-optical methods. The search for the origin of these dipolemoments, however, must be based on the detailed knowledge on protein structure with atomic resolutions. At present, we have two sources of information on the structure of protein molecules: (1) x-ray databases obtained in crystalline state; (2) NMR databases obtained in solution state. While x-ray databases consist of only one model, NMR databases, because of the fluctuation of the protein folding in solution, consist of a number of models, thus enabling the computation of dipole moment repeated for all these models. The aim of this work, using these databases, is the detailed investigation on the interdependence between the structure and dipole moment of protein molecules. The dipole moment of protein molecules has roughly two components: one dipole moment is due to surface charges and the other, core dipole moment, is due to polar groups such as N--H and C==O bonds. The computation of surface charge dipole moment consists of two steps: (A) calculation of the pK shifts of charged groups for electrostatic interactions and (B) calculation of the dipole moment using the pK corrected for electrostatic shifts. The dipole moments of several proteins were computed using both NMR and x-ray databases. The dipole moments of these two sets of calculations are, with a few exceptions, in good agreement with one another and also with measured dipole moments.

The highly refined solution structure of the cytotoxic ribonuclease alpha-sarcin reveals the structural requirements for substrate recognition and ribonucleolytic activity

alpha-Sarcin selectively cleaves a single phosphodiester bond in a universally conserved sequence of the major rRNA, that inactivates the ribosome. The elucidation of the three-dimensional solution structure of this 150 residue enzyme is a crucial step towards understanding alpha-sarcin's conformational stability, ribonucleolytic activity, and its exceptionally high level of specificity. Here, the solution structure has been determined on the basis of 2658 conformationally relevant distances restraints (including stereoespecific assignments) and 119 torsional angular restraints, by nuclear magnetic resonance spectroscopy methods. A total of 60 converged structures have been computed using the program DYANA. The 47 best DYANA structures, following restrained energy minimization by GROMOS, represent the solution structure of alpha-sarcin. The resulting average pairwise root-mean-square-deviation is 0.86 A for backbone atoms and 1.47 A for all heavy atoms. When the more variable regions are excluded from the analysis, the pairwise root-mean-square deviation drops to 0.50 A and 1.00 A, for backbone and heavy atoms, respectively. The alpha-sarcin structure is similar to that reported for restrictocin, although some differences are clearly evident, especially in the loop regions. The average rmsd between the structurally aligned backbones of the 47 final alpha-sarcin structures and the crystal structure of restrictocin is 1.46 A. On the basis of a docking model constructed with alpha-sarcin solution structure and the crystal structure of a 29-nt RNA containing the sarcin/ricin domain, the regions in the protein that could interact specifically with the substrate have been identified. The structural elements that account for the specificity of RNA recognition are located in two separate regions of the protein. One is composed by residues 51 to 55 and loop 5, and the other region, located more than 11 A away in the structure, is the positively charged segment formed by residues 110 to 114.

Role of histidine-50, glutamic acid-96, and histidine-137 in the ribonucleolytic mechanism of the ribotoxin alpha-sarcin

alpha-Sarcin is a ribotoxin secreted by the mold Aspergillus giganteus that degrades the ribosomal RNA by acting as a cyclizing ribonuclease. Three residues potentially involved in the mechanism of catalysis--histidine-50, glutamic acid-96, and histidine-137--were changed to glutamine. Three different single mutation variants (H50Q, E96Q, H137Q) as well as a double variant (H50/137Q) and a triple variant (H50/137Q/E96Q) were prepared and isolated to homogeneity. These variants were spectroscopically (circular dichroism, fluorescence emission, and proton nuclear magnetic resonance) characterized. According to these results, the three-dimensional structure of these variants of alpha-sarcin was preserved; only very minor local changes were detected. All the variants were inactive when assayed against either intact ribosomes or poly(A). The effect of pH on the ribonucleolytic activity of alpha-sarcin was evaluated against the ApA dinucleotide. This assay revealed that only the H50Q variant still retained its ability to cleave a phosphodiester bond, but it did so to a lesser extent than did wild-type alpha-sarcin. The results obtained are interpreted in terms of His137 and Glu96 as essential residues for the catalytic activity of alpha-sarcin (His137 as the general acid and Glu96 as the general base) and His50 stabilizing the transition state of the reaction catalyzed by alpha-sarcin.

The role of high-resolution structural studies in the development of commercial enzymes

Recent developments in both NMR and X-ray crystallography allow the analysis of commercial enzymes in unprecedented detail. The novel methods provide detailed insights into protein dynamics, establish the existence of special catalytic hydrogen bonds and define the ionization states at the enzyme active site. A more detailed understanding of how the changes in structure are related to altered function should facilitate the design of future commercial enzymes with improved performance for different environmental conditions.

Exploring the dynamic information content of a protein NMR structure: comparison of a molecular dynamics simulation with the NMR and X-ray structures of Escherichia coli ribonuclease HI

The multiconformer nature of solution nuclear magnetic resonance (NMR) structures of proteins results from the effects of intramolecular dynamics, spin diffusion and an uneven distribution of structural restraints throughout the molecule. A delineation of the former from the latter two contributions is attempted in this work for an ensemble of 15 NMR structures of the protein Escherichia coli ribonuclease HI (RNase HI). Exploration of the dynamic information content of the NMR ensemble is carried out through correlation with data from two crystal structures and a 1.7-ns molecular dynamics (MD) trajectory of RNase HI in explicit solvent. Assessment of the consistency of the crystal and mean MD structures with nuclear Overhauser effect (NOE) data showed that the NMR ensemble is overall more compatible with the high-resolution (1.48 A) crystal structure than with either the lower-resolution (2.05 A) crystal structure or the MD simulation. Furthermore, the NMR ensemble is found to span more conformational space than the MD simulation for both the backbone and the sidechains of RNase HI. Nonetheless, the backbone conformational variability of both the NMR ensemble and the simulation is especially consistent with NMR relaxation measurements of two loop regions that are putative sites of substrate recognition. Plausible side-chain dynamic information is extracted from the NMR ensemble on the basis of (i) rotamericity and syn-pentane character of variable torsion angles, (ii) comparison of the magnitude of atomic mean-square fluctuations (msf) with those deduced from crystallographic thermal factors, and (iii) comparison of torsion angle conformational behavior in the NMR ensemble and the simulation. Several heterogeneous torsion angles, while adopting non-rotameric/syn-pentane conformations in the NMR ensemble, exist in a unique conformation in the simulation and display low X-ray thermal factors. These torsions are identified as sites whose variability is likely to be an artifact of the NMR structure determination procedure. A number of other torsions show a close correspondence between the conformations sampled in the NMR and MD ensembles, as well as significant correlations among crystallographic thermal factors and atomic msf calculated from the NMR ensemble and the simulation. These results indicate that a significant amount of dynamic information is contained in the NMR ensemble. The relevance of the present findings for the biological function of RNase HI, protein recognition studies, and previous investigations of the motional content of protein NMR structures are discussed.

Characterization of pKa values and titration shifts in the cytotoxic ribonuclease alpha-sarcin by NMR. Relationship between electrostatic interactions, structure, and catalytic function

The electrostatic behavior of titrating groups in alpha-sarcin was investigated using 1H NMR spectroscopy. A total of 209 chemical shift titration curves corresponding to different protons in the molecule were determined over the pH range of 3.0-8.5. Nonlinear least-squares fits of the data to simple relationships derived from the Henderson-Hasselbalch equation led to the unambiguous determination of pKa values for all glutamic acid and histidine residues, as well as for the C-terminal carboxylate and most of the aspartic acids in the free enzyme. The ionization constants of catalytically relevant histidines, His50 and His137, and glutamic acid, Glu96, in the alpha-sarcin-2'-GMP complex were also determined. The pKa values of 15 ionizable groups (C-carboxylate, six aspartic acids, four glutamic acids, and four histidines) were found to be close to their normal values. On the other hand, a number of side chain groups, including those in the active center, showed pKa values far from their intrinsic values. Thus, the pKa values for active site residues His50, Glu96, and His137 were 7.7, 5.2, and 5.8 in the free enzyme and 7.6, approximately 4.8, and 6.8 in the alpha-sarcin-2'-GMP complex, respectively. The pKa values and the activity profile against ApA, as a function of pH, are in agreement with the proposed enzymatic mechanism (in common with RNase T1 and the family of the microbial ribonucleases), in which Glu96 and His137 act as a general base and general acid, respectively. In almost all microbial ribonucleases, a Phe-His interaction is present, which affects the pKa of one of the His residues at the active site (His137). The absence of this interaction in alpha-sarcin would explain the lower pKa value of this His residue, and provides an explanation for the decreased RNase activity of this protein as compared to those of other microbial ribonucleases.

Solution structure of Desulfovibrio vulgaris (Hildenborough) ferrocytochrome c3: structural basis for functional cooperativity

Desulfovibrio vulgaris cytochrome c3 is a 14 kDa tetrahaem cytochrome that plays a central role in energy transduction. The three-dimensional structure of the ferrocytochrome at pH 8.5 was solved through two-dimensional 1H-NMR. The structures were calculated using a large amount of experimental information, which includes upper and lower distance limits as well as dihedral angle restraints. The analysis allows for fast-flipping aromatic residues and flexibility in the haem plane. The structure was determined using 2289 upper and 2390 lower distance limits, 63 restricted ranges for the phi torsion angle, 88 stereospecific assignments out of the 118 stereopairs with non-degenerate chemical shifts (74.6%), and 115 out of the 184 nuclear Overhauser effects to fast-flipping aromatic residues (62.5%), which were pseudo-stereospecifically assigned to one or the other side of the ring. The calculated NMR structures are very well defined, with an average root-mean-square deviation value relative to the mean coordinates of 0.35 A for the backbone atoms and 0.70 A for all heavy-atoms. Comparison of the NMR structures of the ferrocytochrome at pH 8.5 with the available X-ray structure of the ferricytochrome at pH 5.5 reveals that the general fold of the molecule is very similar, but that there are some distinct differences. Calculation of ring current shifts for the residues with significantly different conformations confirms that the NMR structures represent better its solution structure in the reduced form. Some of the localised differences, such as a reorientation of Thr24, are thought to be state-dependent changes that involve alterations in hydrogen bond networks. An important rearrangement in the vicinity of the propionate groups of haem I and involving the covalent linkage of haem II suggests that this is the critical region for the functional cooperativities of this protein.