 Research article
 Open Access
 Published:
Rational polynomial representation of ribonucleotide reductase activity
BMC Biochemistry volume 6, Article number: 8 (2005)
Abstract
Background
Recent data suggest that ribonucleotide reductase (RNR) exists not only as a heterodimer R1_{2}R2_{2} of R1_{2} and R2_{2} homodimers, but also as tetramers R1_{4}R2_{4} and hexamers R1_{6}R2_{6}. Recent data also suggest that ATP binds the R1 subunit at a previously undescribed hexamerization site, in addition to its binding to previously described dimerization and tetramerization sites. Thus, the current view is that R1 has four NDP substrate binding possibilities, four dimerization site binding possibilities (dATP, ATP, dGTP, or dTTP), two tetramerization site binding possibilities (dATP or ATP), and one hexamerization site binding possibility (ATP), in addition to possibilities of unbound site states. This large number of internal R1 states implies an even larger number of quaternary states. A mathematical model of RNR activity which explicitly represents the states of R1 currently exists, but it is complicated in several ways: (1) it includes up to sixfold nested sums; (2) it uses different mathematical structures under different substratemodulator conditions; and (3) it requires root solutions of high order polynomials to determine R1 proportions in mono, di, tetra and hexamer states and thus RNR activity as a function of modulator and total R1 concentrations.
Results
We present four (one for each NDP) rational polynomial models of RNR activity as a function of substrate and reaction rate modifier concentrations. The new models avoid the complications of the earlier model without compromising curve fits to recent data.
Conclusion
Compared to the earlier model of recent data, the new rational polynomial models are simpler, adequately fitting, and likely better suited for biochemical network simulations.
Background
Ribonucleotide reductase (RNR) is a key component of de novo deoxynucleotide (dNTP) metabolism and an important target of cancer therapies [1]. This enzyme, which reduces ribonucleoside diphosphates into corresponding deoxyribonucleoside diphosphates, is exquisitely controlled to properly balance dNTP fluxes in the face of changing scheduled (S phase) and unscheduled (DNA damage/repair) dNTP synthesis demands [2].
Recent data [3–6] suggest that ribonucleotide reductase (RNR) exists not only as a heterodimer R1_{2}R2_{2} of R1_{2} and R2_{2} homodimers [2], but also as a R1_{4}R2_{4} tetramer and as a R1_{6}R2_{6} hexamer, where hexamer formation is driven by ATP binding to a previously undescribed hexamerization site. Thus, in addition to its four substrate binding possibilities in ADP, GDP, CDP, or UDP, and four dimerization/specificity site binding possibilities in dATP, ATP, dGTP, or dTTP, the current view (Figure 1) is that R1 has two tetramerization/inhibitory site binding possibilities in dATP or ATP, and one hexamerization/activation site binding possibility in ATP, in addition to possibilities of unbound site states. The resulting large number of possible R1 states implies an even larger number of quaternary states, and this leads to a complicated mathematical model of RNR activity [3–6]. This model, although useful for explaining RNR activity data, is not useful for biochemical network simulations because: a) it is unwieldy (including up to sixfold nested sums), b) it uses different mathematical structures under different substratemodulator conditions, and c) it requires root solutions of high order polynomials to determine R1 proportions in mono, di, tetra and hexamer quaternary states, and thus RNR activity, as a function of modulator and total R1 concentrations. Simpler mathematical reaction rate models of RNR are needed if deoxynucleotide metabolism [7] is to be represented using Systems Biology Markup Language (SBML) [8–10], a standard which requires single algebraic expression reaction rate laws in some applications [11, 12]. Based on recent data from Cooperman's group [4–6], such expressions are provided here for RNR.
Results
Reaction activities are viewed here as weighted sums of enzyme state specific activities multiplied by probabilities of enzymes being in specific states. For example, a MichaelisMenten reaction rate law is viewed as
where the probability that the enzyme is in a loaded/reactive state (with activity k_{cat}) is P(EA) and the probability that the enzyme is in an empty/unreactive state (with no activity) is P(E).
The RNR models presented here are based on the following four enzyme state probability assumptions:
1. The probability that a particular R1 subunit is bound to NDP is assumed to be
2. The probability that L ∈ {ATP, dATP, dTTP, dGTP} is also bound to the dimerization/specificity site, conditional on NDP binding, is assumed to be
where, if L = dTTP and NDP = GDP for example,
is the probability that dTTP is bound to the dimerization/specificity site given that GDP is bound to the substrate site; the probability of an empty dimerization/specificity site is thus
3. The probability that the tetramerization site is either empty, occupied by ATP, or occupied by dATP, is assumed to be, respectively,
4. Finally, the probability that the hexamerization site is occupied by ATP is assumed to be
The subscripts a (activation), i (inactivation), and s (specificity) on the binding constants correspond to h (hexamerization), t (tetramerization), and d (dimerization) subscripts on P, respectively. That parameter values differ depending upon which NDP substrate is bound to the active site (see Table 1) is indicated by the conditional probabilities.
Previous work [3–6] has shown that the dimer and hexamer states are active, that the tetramer state is slightly active for ADP and GDP and is otherwise inactive, that dimer state activity for CDP and UDP exists when ATP or dATP is bound to the dimerization/specificity site, and that an empty dimerization/specificity site still permits the formation of some dimer with CDP reductase activity, see Table 5 of [5]. Thus, based on the enzyme state probabilities given above, for k_{cat} implicitly defined through
we propose the following expressions:
ADP reduction
GDP reduction
CDP reduction
UDP reduction
In these equations, for ADP and GDP, the first factor is the probability that the dimer site is occupied, and for CDP and UDP, the first factor is the expectation of k_{cat} conditional on R1 being in a dimer state (i.e. having an empty tetramerization site). In the ADP and GDP models, the second factor is the conditional expectation of k_{cat} given that the dimerization site is occupied: the first term of this second factor has in its numerator the statement that k_{cat} = k_{2} if the tetramerization site is empty, or k_{cat} = k_{4} if it is occupied by either dATP or ATP, and the second term states that k_{cat} = k_{6} if the hexamerization site is occupied by ATP. For the CDP and UDP models, the first term of the second factor is the probability of an empty tetramerization site (the event that the corresponding first factor was conditioned on), and the second term states that if the ATP concentration is high enough that the hexamerization/second term dominates the tetramerization/first term whilst the first factor approaches k_{2A}, k_{6} is the overall k_{cat}. This rationale served as our model selection guide. Importantly, the models fitted recent data [3–6] very well, see Figure 2 and Table 1.
Discussion
In general, when an integrated system is engineered from component subsystems, the behavior of the overall system depends on component inputoutput specifications more so than the details of component implementations. By analogy, when enzymological data are applied to biochemical network modeling, rather than the elucidation of reaction mechanisms, it can be expected that the reaction surfaces themselves (i.e. the enzyme's inputoutput characteristics) determine network behavior more so than the details of how such surfaces are represented. Thus, for applications to systems biology, large confidence intervals (CI) in the model parameter estimates of Table 1 (not shown) are not a problem because only goodnessoffit (Fig. 2) really matters; this claim assumes an operating range within the data range, since similarly fitting models often veer apart when used in extrapolations. If reaction mechanism inferences were instead being sought, the large CI in the model parameter estimates would have been a problem, e.g. the squared terms in the model suggest cooperative binding, but this choice provides only slightly better curve fits compared to linear terms (not shown), so cooperative binding cannot be inferred from this model.
In the RNR model presented here, the proportion of R1 units existing in monomer, dimer, tetramer, or hexamer states, and thus the RNR activity per unit enzyme, depends on site binding occupancies but does not depend on the total R1 concentration. In the more complicated previous model [3–6], higher total R1 concentrations favor higher order quaternary states. The degree to which this is so is illustrated by plots of predicted GDP reductase activity as a function of ATP concentration at various R1 concentrations (Figure 3). Consistent with the formation of higher order quaternary R1 states, these plots contract to the left as the total R1 concentrations increase from 1 μM to 100 μM. In future work, the model given here will be altered to capture such trends without losing its simplicity; the total R1 concentration will enter such a model not only as a linear modulator of the reaction surface amplitude (i.e. E_{0} in Eq. 10), but also as a modifier of reaction surface shape parameters, e.g. K_{aATP} will be replaced by a decreasing function of R1.
Conclusion
We identified a rational polynomial model of RNR activity that has single algebraic expressions for each reductase reaction rate law. The expressions provide reasonably good fits (Fig. 2) to recent data [3–6]. Compared to previous reaction rate expressions for this data [3–6], the new expressions are simpler and thus better suited for biochemical network simulations, particularly those constrained to use enzyme reaction rate laws defined as single algebraic expressions [11, 12].
Methods
The parameter estimates shown in Table 1 were obtained through a trialanderror iterative process of nonlinear least squares curve fitting under various, convergence enabling, parameter fixations (i.e. profile searches). In the end, the curve fits were those of Figure 2 with corresponding parameter estimates in Table 1; large 95% confidence intervals (not shown) allowed rounding of the parameter estimates to somewhat arbitrary choices. Nonlinear least squares parameter estimations were performed using the optimization method of Nelder and Mead [13] and the statistical computing environment R [14]. All parameters were estimated as e^{c} to assure positive values. For additional details, R scripts are available with the data as supplementary material [15].
Abbreviations
 RNR:

= ribonucleotide reductase
 dNTP:

= deoxynucleoside tripshospate
 dNDP:

= deoxynucleoside dipshospate
 a :

= activation
 i :

= inactivation
 s :

= specificity
 h :

= hexamerization
 t :

= tetramerization
 d :

= dimerization
 CI:

= confidence interval
 SBML:

= Systems Biology Markup Language.
References
 1.
Tsimberidou AM, Alvarado Y, Giles FJ: Evolving role of ribonucleoside reductase inhibitors in hematologic malignancies. Expert Rev Anticancer Ther. 2002, 2 (4): 437448. 10.1586/14737140.2.4.437.
 2.
Kolberg M, Strand KR, Graff P, Andersson KK: Structure, function, and mechanism of ribonucleotide reductases. Biochim Biophys Acta. 2004, 1699 (12): 134.
 3.
Cooperman BS, Kashlan OB: A comprehensive model for the allosteric regulation of Class Ia ribonucleotide reductases. Adv Enzyme Regul. 2003, 43: 167182. 10.1016/S00652571(02)000353.
 4.
Kashlan OB, Cooperman BS: Comprehensive model for allosteric regulation of mammalian ribonucleotide reductase: refinements and consequences. Biochemistry. 2003, 42 (6): 16961706. 10.1021/bi020634d.
 5.
Kashlan OB, Scott CP, Lear JD, Cooperman BS: A comprehensive model for the allosteric regulation of mammalian ribonucleotide reductase. Functional consequences of ATP and dATPinduced oligomerization of the large subunit. Biochemistry. 2002, 41 (2): 462474. 10.1021/bi011653a.
 6.
Scott CP, Kashlan OB, Lear JD, Cooperman BS: A quantitative model for allosteric control of purine reduction by murine ribonucleotide reductase. Biochemistry. 2001, 40 (6): 16511661. 10.1021/bi002335z.
 7.
Jackson RC: The Theoretical Foundations of Cancer Chemotherapy Introduced by Computer Models. 1992, New York , Academic Press, 447
 8.
Hucka M, Finney A, Sauro HM, Bolouri H, Doyle JC, Kitano H, Arkin AP, Bornstein BJ, Bray D, CornishBowden A, Cuellar AA, Dronov S, Gilles ED, Ginkel M, Gor V, Goryanin II, Hedley WJ, Hodgman TC, Hofmeyr JH, Hunter PJ, Juty NS, Kasberger JL, Kremling A, Kummer U, Le Novere N, Loew LM, Lucio D, Mendes P, Minch E, Mjolsness ED, Nakayama Y, Nelson MR, Nielsen PF, Sakurada T, Schaff JC, Shapiro BE, Shimizu TS, Spence HD, Stelling J, Takahashi K, Tomita M, Wagner J, Wang J: The systems biology markup language (SBML): a medium for representation and exchange of biochemical network models. Bioinformatics. 2003, 19 (4): 524531. 10.1093/bioinformatics/btg015.
 9.
Finney A, Hucka M: Systems biology markup language: Level 2 and beyond. Biochem Soc Trans. 2003, 31 (Pt 6): 14721473.
 10.
Systems Biology Markup Language. [http://sbml.org]
 11.
Radivoyevitch T: A twoway interface between limited Systems Biology Markup Language and R. BMC Bioinformatics. 2004, 5 (1): 19010.1186/147121055190.
 12.
Radivoyevitch T: SBMLR. [http://www.bioconductor.org/repository/devel/package/html/SBMLR.html]
 13.
Nelder JA, Mead R: A simplex algorithm for function minimization. Computer Journal. 1965, 7: 308313.
 14.
Ihaka R, Gentleman R: R:a language for data analysis and graphics. Journal of Computational and graphical statistics. 1996, 5: 299314.
 15.
Radivoyevitch T, Kashlan OB, Cooperman BS: Rational Polynomial Representation of Ribonucleotide Reductase Activity. [http://epbiradivot.cwru.edu/cooperman/]
Acknowledgements
We thank Dr. Charles Scott for sharing his data. TR was supported by the Biostatistics Core Facility of the Comprehensive Cancer Center of Case Western Reserve University and University Hospitals of Cleveland (P30 CA43703), the American Cancer Society (IRG9102209), the National Cancer Institute's Integrative Cancer Biology Program (P20 CA11296301) and NIH grant 1K25 CA10479101A1. BSC was supported by NIH grant CA 58567. OBK was a recipient of a postdoctoral fellowship award from the PennsylvaniaDelaware Affiliate of the AHA, and was supported by grants DK061296 and DK066883 from the NIH.
Author information
Affiliations
Corresponding author
Additional information
Authors' contributions
TR performed the curve fits to the data and explored various model choices.
OBK and BSC provided the original model, its simulations (Figure 3) and the data.
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.
Rights and permissions
About this article
Cite this article
Radivoyevitch, T., Kashlan, O.B. & Cooperman, B.S. Rational polynomial representation of ribonucleotide reductase activity. BMC Biochem 6, 8 (2005). https://doi.org/10.1186/1471209168
Received:
Accepted:
Published:
Keywords
 Ribonucleotide Reductase
 System Biology Markup Language
 Dimerization Site
 Ribonucleotide Reductase Activity
 Quaternary State