Links to Exploration step
Le document en format XML
<record><TEI><teiHeader><fileDesc><titleStmt><title xml:lang="en">Automated mass action model space generation and analysis methods for two-reactant combinatorially complex equilibriums: An analysis of ATP-induced ribonucleotide reductase R1 hexamerization data</title>
<author><name sortKey="Radivoyevitch, Tomas" sort="Radivoyevitch, Tomas" uniqKey="Radivoyevitch T" first="Tomas" last="Radivoyevitch">Tomas Radivoyevitch</name>
<affiliation><nlm:aff id="I1">Department of Epidemiology and Biostatistics, Case Western Reserve University, Cleveland, Ohio 44106, USA</nlm:aff>
</affiliation>
</author>
</titleStmt>
<publicationStmt><idno type="wicri:source">PMC</idno>
<idno type="pmid">20003203</idno>
<idno type="pmc">2799446</idno>
<idno type="url">http://www.ncbi.nlm.nih.gov/pmc/articles/PMC2799446</idno>
<idno type="RBID">PMC:2799446</idno>
<idno type="doi">10.1186/1745-6150-4-50</idno>
<date when="2009">2009</date>
<idno type="wicri:Area/Pmc/Corpus">000B28</idno>
<idno type="wicri:explorRef" wicri:stream="Pmc" wicri:step="Corpus" wicri:corpus="PMC">000B28</idno>
</publicationStmt>
<sourceDesc><biblStruct><analytic><title xml:lang="en" level="a" type="main">Automated mass action model space generation and analysis methods for two-reactant combinatorially complex equilibriums: An analysis of ATP-induced ribonucleotide reductase R1 hexamerization data</title>
<author><name sortKey="Radivoyevitch, Tomas" sort="Radivoyevitch, Tomas" uniqKey="Radivoyevitch T" first="Tomas" last="Radivoyevitch">Tomas Radivoyevitch</name>
<affiliation><nlm:aff id="I1">Department of Epidemiology and Biostatistics, Case Western Reserve University, Cleveland, Ohio 44106, USA</nlm:aff>
</affiliation>
</author>
</analytic>
<series><title level="j">Biology Direct</title>
<idno type="eISSN">1745-6150</idno>
<imprint><date when="2009">2009</date>
</imprint>
</series>
</biblStruct>
</sourceDesc>
</fileDesc>
<profileDesc><textClass></textClass>
</profileDesc>
</teiHeader>
<front><div type="abstract" xml:lang="en"><sec><title>Background</title>
<p>Ribonucleotide reductase is the main control point of dNTP production. It has two subunits, R1, and R2 or p53R2. R1 has 5 possible catalytic site states (empty or filled with 1 of 4 NDPs), 5 possible <italic>s</italic>
-site states (empty or filled with ATP, dATP, dTTP or dGTP), 3 possible <italic>a</italic>
-site states (empty or filled with ATP or dATP), perhaps two possible <italic>h</italic>
-site states (empty or filled with ATP), and all of this is folded into an R1 monomer-dimer-tetramer-hexamer equilibrium where R1 j-mers can be bound by variable numbers of R2 or p53R2 dimers. Trillions of RNR complexes are possible as a result. The problem is to determine which are needed in models to explain available data. This problem is intractable for 10 reactants, but it can be solved for 2 and is here for R1 and ATP.</p>
</sec>
<sec><title>Results</title>
<p>Thousands of ATP-induced R1 hexamerization models with up to three (<italic>s</italic>
, <italic>a </italic>
and <italic>h</italic>
) ATP binding sites per R1 subunit were automatically generated via hypotheses that complete dissociation constants are infinite and/or that binary dissociation constants are equal. To limit the model space size, it was assumed that <italic>s</italic>
-sites are always filled in oligomers and never filled in monomers, and to interpret model terms it was assumed that <italic>a</italic>
-sites fill before <italic>h</italic>
-sites. The models were fitted to published dynamic light scattering data. As the lowest Akaike Information Criterion (AIC) of the 3-parameter models was greater than the lowest of the 2-parameter models, only models with up to 3 parameters were fitted. Models with sums of squared errors less than twice the minimum were then partitioned into two groups: those that contained no occupied <italic>h</italic>
-site terms (508 models) and those that contained at least one (1580 models). Normalized AIC densities of these two groups of models differed significantly in favor of models that did not include an <italic>h</italic>
-site term (Kolmogorov-Smirnov p < 1 × 10<sup>-15</sup>
); consistent with this, 28 of the top 30 models (ranked by AICs) did not include an <italic>h</italic>
-site term and 28/30 > 508/2088 with p < 2 × 10<sup>-15</sup>
. Finally, 99 of the 2088 models did not have any terms with ATP/R1 ratios >1.5, but of the top 30, there were 14 such models (14/30 > 99/2088 with p < 3 × 10<sup>-16</sup>
), i.e. the existence of R1 hexamers with >3 <italic>a</italic>
-sites occupied by ATP is also not supported by this dataset.</p>
</sec>
<sec><title>Conclusion</title>
<p>The analysis presented suggests that three <italic>a</italic>
-sites may not be occupied by ATP in R1 hexamers under the conditions of the data analyzed. If <italic>a</italic>
-sites fill before <italic>h</italic>
-sites, this implies that the dataset analyzed can be explained without the existence of an <italic>h</italic>
-site.</p>
</sec>
<sec><title>Reviewers</title>
<p>This article was reviewed by Ossama Kashlan <italic>(nominated by Philip Hahnfeldt)</italic>
, Bin Hu <italic>(nominated by William Hlavacek) </italic>
and Rainer Sachs.</p>
</sec>
</div>
</front>
<back><div1 type="bibliography"><listBibl><biblStruct><analytic><author><name sortKey="Radivoyevitch, T" uniqKey="Radivoyevitch T">T Radivoyevitch</name>
</author>
<author><name sortKey="Loparo, Ka" uniqKey="Loparo K">KA Loparo</name>
</author>
<author><name sortKey="Jackson, Rc" uniqKey="Jackson R">RC Jackson</name>
</author>
<author><name sortKey="Sedwick, Wd" uniqKey="Sedwick W">WD Sedwick</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Manegold, C" uniqKey="Manegold C">C Manegold</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Kroep, Jr" uniqKey="Kroep J">JR Kroep</name>
</author>
<author><name sortKey="Van Moorsel, Cj" uniqKey="Van Moorsel C">CJ Van Moorsel</name>
</author>
<author><name sortKey="Veerman, G" uniqKey="Veerman G">G Veerman</name>
</author>
<author><name sortKey="Voorn, Da" uniqKey="Voorn D">DA Voorn</name>
</author>
<author><name sortKey="Schultz, Rm" uniqKey="Schultz R">RM Schultz</name>
</author>
<author><name sortKey="Worzalla, Jf" uniqKey="Worzalla J">JF Worzalla</name>
</author>
<author><name sortKey="Tanzer, Lr" uniqKey="Tanzer L">LR Tanzer</name>
</author>
<author><name sortKey="Merriman, Rl" uniqKey="Merriman R">RL Merriman</name>
</author>
<author><name sortKey="Pinedo, Hm" uniqKey="Pinedo H">HM Pinedo</name>
</author>
<author><name sortKey="Peters, Gj" uniqKey="Peters G">GJ Peters</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Shewach, Ds" uniqKey="Shewach D">DS Shewach</name>
</author>
<author><name sortKey="Reynolds, Kk" uniqKey="Reynolds K">KK Reynolds</name>
</author>
<author><name sortKey="Hertel, L" uniqKey="Hertel L">L Hertel</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Wang, J" uniqKey="Wang J">J Wang</name>
</author>
<author><name sortKey="Lohman, Gj" uniqKey="Lohman G">GJ Lohman</name>
</author>
<author><name sortKey="Stubbe, J" uniqKey="Stubbe J">J Stubbe</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Ferrandina, G" uniqKey="Ferrandina G">G Ferrandina</name>
</author>
<author><name sortKey="Mey, V" uniqKey="Mey V">V Mey</name>
</author>
<author><name sortKey="Nannizzi, S" uniqKey="Nannizzi S">S Nannizzi</name>
</author>
<author><name sortKey="Ricciardi, S" uniqKey="Ricciardi S">S Ricciardi</name>
</author>
<author><name sortKey="Petrillo, M" uniqKey="Petrillo M">M Petrillo</name>
</author>
<author><name sortKey="Ferlini, C" uniqKey="Ferlini C">C Ferlini</name>
</author>
<author><name sortKey="Danesi, R" uniqKey="Danesi R">R Danesi</name>
</author>
<author><name sortKey="Scambia, G" uniqKey="Scambia G">G Scambia</name>
</author>
<author><name sortKey="Del Tacca, M" uniqKey="Del Tacca M">M Del Tacca</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Ohtaka, K" uniqKey="Ohtaka K">K Ohtaka</name>
</author>
<author><name sortKey="Kohya, N" uniqKey="Kohya N">N Kohya</name>
</author>
<author><name sortKey="Sato, K" uniqKey="Sato K">K Sato</name>
</author>
<author><name sortKey="Kitajima, Y" uniqKey="Kitajima Y">Y Kitajima</name>
</author>
<author><name sortKey="Ide, T" uniqKey="Ide T">T Ide</name>
</author>
<author><name sortKey="Mitsuno, M" uniqKey="Mitsuno M">M Mitsuno</name>
</author>
<author><name sortKey="Miyazaki, K" uniqKey="Miyazaki K">K Miyazaki</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Flanagan, Sa" uniqKey="Flanagan S">SA Flanagan</name>
</author>
<author><name sortKey="Robinson, Bw" uniqKey="Robinson B">BW Robinson</name>
</author>
<author><name sortKey="Krokosky, Cm" uniqKey="Krokosky C">CM Krokosky</name>
</author>
<author><name sortKey="Shewach, Ds" uniqKey="Shewach D">DS Shewach</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Muller, A" uniqKey="Muller A">A Muller</name>
</author>
<author><name sortKey="Fishel, R" uniqKey="Fishel R">R Fishel</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Leung, Wk" uniqKey="Leung W">WK Leung</name>
</author>
<author><name sortKey="Kim, Jj" uniqKey="Kim J">JJ Kim</name>
</author>
<author><name sortKey="Kim, Jg" uniqKey="Kim J">JG Kim</name>
</author>
<author><name sortKey="Graham, Dy" uniqKey="Graham D">DY Graham</name>
</author>
<author><name sortKey="Sepulveda, Ar" uniqKey="Sepulveda A">AR Sepulveda</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Nakata, B" uniqKey="Nakata B">B Nakata</name>
</author>
<author><name sortKey="Wang, Yq" uniqKey="Wang Y">YQ Wang</name>
</author>
<author><name sortKey="Yashiro, M" uniqKey="Yashiro M">M Yashiro</name>
</author>
<author><name sortKey="Nishioka, N" uniqKey="Nishioka N">N Nishioka</name>
</author>
<author><name sortKey="Tanaka, H" uniqKey="Tanaka H">H Tanaka</name>
</author>
<author><name sortKey="Ohira, M" uniqKey="Ohira M">M Ohira</name>
</author>
<author><name sortKey="Ishikawa, T" uniqKey="Ishikawa T">T Ishikawa</name>
</author>
<author><name sortKey="Nishino, H" uniqKey="Nishino H">H Nishino</name>
</author>
<author><name sortKey="Hirakawa, K" uniqKey="Hirakawa K">K Hirakawa</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Hartmann, A" uniqKey="Hartmann A">A Hartmann</name>
</author>
<author><name sortKey="Zanardo, L" uniqKey="Zanardo L">L Zanardo</name>
</author>
<author><name sortKey="Bocker Edmonston, T" uniqKey="Bocker Edmonston T">T Bocker-Edmonston</name>
</author>
<author><name sortKey="Blaszyk, H" uniqKey="Blaszyk H">H Blaszyk</name>
</author>
<author><name sortKey="Dietmaier, W" uniqKey="Dietmaier W">W Dietmaier</name>
</author>
<author><name sortKey="Stoehr, R" uniqKey="Stoehr R">R Stoehr</name>
</author>
<author><name sortKey="Cheville, Jc" uniqKey="Cheville J">JC Cheville</name>
</author>
<author><name sortKey="Junker, K" uniqKey="Junker K">K Junker</name>
</author>
<author><name sortKey="Wieland, W" uniqKey="Wieland W">W Wieland</name>
</author>
<author><name sortKey="Knuechel R Rueschoff, J" uniqKey="Knuechel R Rueschoff J">J Knuechel R Rueschoff</name>
</author>
<author><name sortKey="Hofstaedter, F" uniqKey="Hofstaedter F">F Hofstaedter</name>
</author>
<author><name sortKey="Fishel, R" uniqKey="Fishel R">R Fishel</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Ichikawa, Y" uniqKey="Ichikawa Y">Y Ichikawa</name>
</author>
<author><name sortKey="Lemon, Sj" uniqKey="Lemon S">SJ Lemon</name>
</author>
<author><name sortKey="Wang, S" uniqKey="Wang S">S Wang</name>
</author>
<author><name sortKey="Franklin, B" uniqKey="Franklin B">B Franklin</name>
</author>
<author><name sortKey="Watson, P" uniqKey="Watson P">P Watson</name>
</author>
<author><name sortKey="Knezetic, Ja" uniqKey="Knezetic J">JA Knezetic</name>
</author>
<author><name sortKey="Bewtra, C" uniqKey="Bewtra C">C Bewtra</name>
</author>
<author><name sortKey="Lynch, Ht" uniqKey="Lynch H">HT Lynch</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Simpkins, Sb" uniqKey="Simpkins S">SB Simpkins</name>
</author>
<author><name sortKey="Bocker, T" uniqKey="Bocker T">T Bocker</name>
</author>
<author><name sortKey="Swisher, Em" uniqKey="Swisher E">EM Swisher</name>
</author>
<author><name sortKey="Mutch, Dg" uniqKey="Mutch D">DG Mutch</name>
</author>
<author><name sortKey="Gersell, Dj" uniqKey="Gersell D">DJ Gersell</name>
</author>
<author><name sortKey="Kovatich, Aj" uniqKey="Kovatich A">AJ Kovatich</name>
</author>
<author><name sortKey="Palazzo, Jp" uniqKey="Palazzo J">JP Palazzo</name>
</author>
<author><name sortKey="Fishel, R" uniqKey="Fishel R">R Fishel</name>
</author>
<author><name sortKey="Goodfellow, Pj" uniqKey="Goodfellow P">PJ Goodfellow</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Leung, Sy" uniqKey="Leung S">SY Leung</name>
</author>
<author><name sortKey="Chan, Tl" uniqKey="Chan T">TL Chan</name>
</author>
<author><name sortKey="Chung, Lp" uniqKey="Chung L">LP Chung</name>
</author>
<author><name sortKey="Chan, As" uniqKey="Chan A">AS Chan</name>
</author>
<author><name sortKey="Fan, Yw" uniqKey="Fan Y">YW Fan</name>
</author>
<author><name sortKey="Hung, Kn" uniqKey="Hung K">KN Hung</name>
</author>
<author><name sortKey="Kwong, Wk" uniqKey="Kwong W">WK Kwong</name>
</author>
<author><name sortKey="Ho, Jw" uniqKey="Ho J">JW Ho</name>
</author>
<author><name sortKey="Yuen, St" uniqKey="Yuen S">ST Yuen</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Thelander, L" uniqKey="Thelander L">L Thelander</name>
</author>
<author><name sortKey="Reichard, P" uniqKey="Reichard P">P Reichard</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Tanaka, H" uniqKey="Tanaka H">H Tanaka</name>
</author>
<author><name sortKey="Arakawa, H" uniqKey="Arakawa H">H Arakawa</name>
</author>
<author><name sortKey="Yamaguchi, T" uniqKey="Yamaguchi T">T Yamaguchi</name>
</author>
<author><name sortKey="Shiraishi, K" uniqKey="Shiraishi K">K Shiraishi</name>
</author>
<author><name sortKey="Fukuda, S" uniqKey="Fukuda S">S Fukuda</name>
</author>
<author><name sortKey="Matsui, K" uniqKey="Matsui K">K Matsui</name>
</author>
<author><name sortKey="Takei, Y" uniqKey="Takei Y">Y Takei</name>
</author>
<author><name sortKey="Nakamura, Y" uniqKey="Nakamura Y">Y Nakamura</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Nakano, K" uniqKey="Nakano K">K Nakano</name>
</author>
<author><name sortKey="Balint, E" uniqKey="Balint E">E Balint</name>
</author>
<author><name sortKey="Ashcroft, M" uniqKey="Ashcroft M">M Ashcroft</name>
</author>
<author><name sortKey="Vousden, Kh" uniqKey="Vousden K">KH Vousden</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Feng, Ij" uniqKey="Feng I">IJ Feng</name>
</author>
<author><name sortKey="Radivoyevitch, T" uniqKey="Radivoyevitch T">T Radivoyevitch</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Kashlan, Ob" uniqKey="Kashlan O">OB Kashlan</name>
</author>
<author><name sortKey="Scott, Cp" uniqKey="Scott C">CP Scott</name>
</author>
<author><name sortKey="Lear, Jd" uniqKey="Lear J">JD Lear</name>
</author>
<author><name sortKey="Cooperman, Bs" uniqKey="Cooperman B">BS Cooperman</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Ihaka, R" uniqKey="Ihaka R">R Ihaka</name>
</author>
<author><name sortKey="Gentleman, R" uniqKey="Gentleman R">R Gentleman</name>
</author>
</analytic>
</biblStruct>
<biblStruct></biblStruct>
<biblStruct><analytic><author><name sortKey="Akaike, H" uniqKey="Akaike H">H Akaike</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Burnham, Kp" uniqKey="Burnham K">KP Burnham</name>
</author>
<author><name sortKey="Anderson, Dr" uniqKey="Anderson D">DR Anderson</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Radivoyevitch, T" uniqKey="Radivoyevitch T">T Radivoyevitch</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Xu, H" uniqKey="Xu H">H Xu</name>
</author>
<author><name sortKey="Faber, C" uniqKey="Faber C">C Faber</name>
</author>
<author><name sortKey="Uchiki, T" uniqKey="Uchiki T">T Uchiki</name>
</author>
<author><name sortKey="Fairman, Jw" uniqKey="Fairman J">JW Fairman</name>
</author>
<author><name sortKey="Racca, J" uniqKey="Racca J">J Racca</name>
</author>
<author><name sortKey="Dealwis, C" uniqKey="Dealwis C">C Dealwis</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Rofougaran, R" uniqKey="Rofougaran R">R Rofougaran</name>
</author>
<author><name sortKey="Vodnala, M" uniqKey="Vodnala M">M Vodnala</name>
</author>
<author><name sortKey="Hofer, A" uniqKey="Hofer A">A Hofer</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Roach, Ls" uniqKey="Roach L">LS Roach</name>
</author>
<author><name sortKey="Song, H" uniqKey="Song H">H Song</name>
</author>
<author><name sortKey="Ismagilov, Rf" uniqKey="Ismagilov R">RF Ismagilov</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Song, H" uniqKey="Song H">H Song</name>
</author>
<author><name sortKey="Ismagilov, Rf" uniqKey="Ismagilov R">RF Ismagilov</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Huebner, A" uniqKey="Huebner A">A Huebner</name>
</author>
<author><name sortKey="Olguin, Lf" uniqKey="Olguin L">LF Olguin</name>
</author>
<author><name sortKey="Bratton, D" uniqKey="Bratton D">D Bratton</name>
</author>
<author><name sortKey="Whyte, G" uniqKey="Whyte G">G Whyte</name>
</author>
<author><name sortKey="Huck, Wt" uniqKey="Huck W">WT Huck</name>
</author>
<author><name sortKey="De Mello, Aj" uniqKey="De Mello A">AJ de Mello</name>
</author>
<author><name sortKey="Edel, Jb" uniqKey="Edel J">JB Edel</name>
</author>
<author><name sortKey="Abell, C" uniqKey="Abell C">C Abell</name>
</author>
<author><name sortKey="Hollfelder, F" uniqKey="Hollfelder F">F Hollfelder</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Huebner, A" uniqKey="Huebner A">A Huebner</name>
</author>
<author><name sortKey="Bratton, D" uniqKey="Bratton D">D Bratton</name>
</author>
<author><name sortKey="Whyte, G" uniqKey="Whyte G">G Whyte</name>
</author>
<author><name sortKey="Yang, M" uniqKey="Yang M">M Yang</name>
</author>
<author><name sortKey="Demello, Aj" uniqKey="Demello A">AJ Demello</name>
</author>
<author><name sortKey="Abell, C" uniqKey="Abell C">C Abell</name>
</author>
<author><name sortKey="Hollfelder, F" uniqKey="Hollfelder F">F Hollfelder</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Atkinson, Ac" uniqKey="Atkinson A">AC Atkinson</name>
</author>
<author><name sortKey="Donev, An" uniqKey="Donev A">AN Donev</name>
</author>
<author><name sortKey="Tobias, Rd" uniqKey="Tobias R">RD Tobias</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Hucka, M" uniqKey="Hucka M">M Hucka</name>
</author>
<author><name sortKey="Finney, A" uniqKey="Finney A">A Finney</name>
</author>
<author><name sortKey="Sauro, Hm" uniqKey="Sauro H">HM Sauro</name>
</author>
<author><name sortKey="Bolouri, H" uniqKey="Bolouri H">H Bolouri</name>
</author>
<author><name sortKey="Doyle, Jc" uniqKey="Doyle J">JC Doyle</name>
</author>
<author><name sortKey="Kitano, H" uniqKey="Kitano H">H Kitano</name>
</author>
<author><name sortKey="Arkin, Ap" uniqKey="Arkin A">AP Arkin</name>
</author>
<author><name sortKey="Bornstein, Bj" uniqKey="Bornstein B">BJ Bornstein</name>
</author>
<author><name sortKey="Bray, D" uniqKey="Bray D">D Bray</name>
</author>
<author><name sortKey="Cornish Bowden, A" uniqKey="Cornish Bowden A">A Cornish-Bowden</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Le Novere, N" uniqKey="Le Novere N">N Le Novere</name>
</author>
<author><name sortKey="Bornstein, B" uniqKey="Bornstein B">B Bornstein</name>
</author>
<author><name sortKey="Broicher, A" uniqKey="Broicher A">A Broicher</name>
</author>
<author><name sortKey="Courtot, M" uniqKey="Courtot M">M Courtot</name>
</author>
<author><name sortKey="Donizelli, M" uniqKey="Donizelli M">M Donizelli</name>
</author>
<author><name sortKey="Dharuri, H" uniqKey="Dharuri H">H Dharuri</name>
</author>
<author><name sortKey="Li, L" uniqKey="Li L">L Li</name>
</author>
<author><name sortKey="Sauro, H" uniqKey="Sauro H">H Sauro</name>
</author>
<author><name sortKey="Schilstra, M" uniqKey="Schilstra M">M Schilstra</name>
</author>
<author><name sortKey="Shapiro B Snoep, Jl" uniqKey="Shapiro B Snoep J">JL Shapiro B Snoep</name>
</author>
<author><name sortKey="Hucka, M" uniqKey="Hucka M">M Hucka</name>
</author>
</analytic>
</biblStruct>
<biblStruct></biblStruct>
</listBibl>
</div1>
</back>
</TEI>
<pmc article-type="research-article"><pmc-dir>properties open_access</pmc-dir>
<front><journal-meta><journal-id journal-id-type="nlm-ta">Biol Direct</journal-id>
<journal-title-group><journal-title>Biology Direct</journal-title>
</journal-title-group>
<issn pub-type="epub">1745-6150</issn>
<publisher><publisher-name>BioMed Central</publisher-name>
</publisher>
</journal-meta>
<article-meta><article-id pub-id-type="pmid">20003203</article-id>
<article-id pub-id-type="pmc">2799446</article-id>
<article-id pub-id-type="publisher-id">1745-6150-4-50</article-id>
<article-id pub-id-type="doi">10.1186/1745-6150-4-50</article-id>
<article-categories><subj-group subj-group-type="heading"><subject>Research</subject>
</subj-group>
</article-categories>
<title-group><article-title>Automated mass action model space generation and analysis methods for two-reactant combinatorially complex equilibriums: An analysis of ATP-induced ribonucleotide reductase R1 hexamerization data</article-title>
</title-group>
<contrib-group><contrib contrib-type="author" corresp="yes" id="A1"><name><surname>Radivoyevitch</surname>
<given-names>Tomas</given-names>
</name>
<xref ref-type="aff" rid="I1">1</xref>
<email>txr24@case.edu</email>
</contrib>
</contrib-group>
<aff id="I1"><label>1</label>
Department of Epidemiology and Biostatistics, Case Western Reserve University, Cleveland, Ohio 44106, USA</aff>
<pub-date pub-type="collection"><year>2009</year>
</pub-date>
<pub-date pub-type="epub"><day>9</day>
<month>12</month>
<year>2009</year>
</pub-date>
<volume>4</volume>
<fpage>50</fpage>
<lpage>50</lpage>
<history><date date-type="received"><day>2</day>
<month>12</month>
<year>2009</year>
</date>
<date date-type="accepted"><day>9</day>
<month>12</month>
<year>2009</year>
</date>
</history>
<permissions><copyright-statement>Copyright ©2009 Radivoyevitch; licensee BioMed Central Ltd.</copyright-statement>
<copyright-year>2009</copyright-year>
<copyright-holder>Radivoyevitch; licensee BioMed Central Ltd.</copyright-holder>
<license license-type="open-access" xlink:href="http://creativecommons.org/licenses/by/2.0"><license-p>This is an Open Access article distributed under the terms of the Creative Commons Attribution License (<ext-link ext-link-type="uri" xlink:href="http://creativecommons.org/licenses/by/2.0">http://creativecommons.org/licenses/by/2.0</ext-link>
), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.</license-p>
</license>
</permissions>
<self-uri xlink:href="http://www.biology-direct.com/content/4/1/50"></self-uri>
<abstract><sec><title>Background</title>
<p>Ribonucleotide reductase is the main control point of dNTP production. It has two subunits, R1, and R2 or p53R2. R1 has 5 possible catalytic site states (empty or filled with 1 of 4 NDPs), 5 possible <italic>s</italic>
-site states (empty or filled with ATP, dATP, dTTP or dGTP), 3 possible <italic>a</italic>
-site states (empty or filled with ATP or dATP), perhaps two possible <italic>h</italic>
-site states (empty or filled with ATP), and all of this is folded into an R1 monomer-dimer-tetramer-hexamer equilibrium where R1 j-mers can be bound by variable numbers of R2 or p53R2 dimers. Trillions of RNR complexes are possible as a result. The problem is to determine which are needed in models to explain available data. This problem is intractable for 10 reactants, but it can be solved for 2 and is here for R1 and ATP.</p>
</sec>
<sec><title>Results</title>
<p>Thousands of ATP-induced R1 hexamerization models with up to three (<italic>s</italic>
, <italic>a </italic>
and <italic>h</italic>
) ATP binding sites per R1 subunit were automatically generated via hypotheses that complete dissociation constants are infinite and/or that binary dissociation constants are equal. To limit the model space size, it was assumed that <italic>s</italic>
-sites are always filled in oligomers and never filled in monomers, and to interpret model terms it was assumed that <italic>a</italic>
-sites fill before <italic>h</italic>
-sites. The models were fitted to published dynamic light scattering data. As the lowest Akaike Information Criterion (AIC) of the 3-parameter models was greater than the lowest of the 2-parameter models, only models with up to 3 parameters were fitted. Models with sums of squared errors less than twice the minimum were then partitioned into two groups: those that contained no occupied <italic>h</italic>
-site terms (508 models) and those that contained at least one (1580 models). Normalized AIC densities of these two groups of models differed significantly in favor of models that did not include an <italic>h</italic>
-site term (Kolmogorov-Smirnov p < 1 × 10<sup>-15</sup>
); consistent with this, 28 of the top 30 models (ranked by AICs) did not include an <italic>h</italic>
-site term and 28/30 > 508/2088 with p < 2 × 10<sup>-15</sup>
. Finally, 99 of the 2088 models did not have any terms with ATP/R1 ratios >1.5, but of the top 30, there were 14 such models (14/30 > 99/2088 with p < 3 × 10<sup>-16</sup>
), i.e. the existence of R1 hexamers with >3 <italic>a</italic>
-sites occupied by ATP is also not supported by this dataset.</p>
</sec>
<sec><title>Conclusion</title>
<p>The analysis presented suggests that three <italic>a</italic>
-sites may not be occupied by ATP in R1 hexamers under the conditions of the data analyzed. If <italic>a</italic>
-sites fill before <italic>h</italic>
-sites, this implies that the dataset analyzed can be explained without the existence of an <italic>h</italic>
-site.</p>
</sec>
<sec><title>Reviewers</title>
<p>This article was reviewed by Ossama Kashlan <italic>(nominated by Philip Hahnfeldt)</italic>
, Bin Hu <italic>(nominated by William Hlavacek) </italic>
and Rainer Sachs.</p>
</sec>
</abstract>
</article-meta>
</front>
<body><sec><title>Background</title>
<sec><title>Introduction</title>
<p>The dNTP supply system produces dNTPs at rates that match those demanded by DNA replication and repair. With respect to ribose ring moieties, it is comprised of both a <italic>de novo </italic>
system whose substrates are ribonucleoside diphosphates (NDPs) and a salvage system whose substrates are deoxynucleosides (dNs). The <italic>de novo </italic>
system includes ribonucleotide reductase (RNR), deoxycytidylate deaminase (DCTD), and thymidylate synthetase (TS), and the salvage system includes deoxycytidine kinase (dCK), thymidine kinase 1 (TK1), deoxyguanosine kinase (dGK) and thymidine kinase 2 (TK2), see Fig. <xref ref-type="fig" rid="F1">1</xref>
. The dNTP supply system is important because many anticancer agents target or traverse it (e.g. gemcitabine, hydroxyurea, triapine, 5-FU) or damage DNA directly (e.g. ionizing radiation, alkylating agents, oxaliplatin) and thus place demands on it for replacement dNTPs. In the future, mathematical models of cancer relevant systems will be needed to optimize multi-agent anticancer dose timings [<xref ref-type="bibr" rid="B1">1</xref>
]. For example, gemcitabine (dFdC, diflourodeoxycytdine) [<xref ref-type="bibr" rid="B2">2</xref>
] absorption is rate limited by dCK [<xref ref-type="bibr" rid="B3">3</xref>
,<xref ref-type="bibr" rid="B4">4</xref>
], dFdC targets RNR [<xref ref-type="bibr" rid="B5">5</xref>
], dFdC resistance is associated with RNR over expression [<xref ref-type="bibr" rid="B6">6</xref>
,<xref ref-type="bibr" rid="B7">7</xref>
], and differential ionizing radiation (IR) sensitivity that dFdC imparts onto mismatch repair (MMR) defective cells may be due to mismatches caused by dNTP pool imbalances caused by RNR inhibition, rather than differential dFdC incorporation into DNA [<xref ref-type="bibr" rid="B8">8</xref>
], so mathematical models of dNTP supply will be needed to optimize dFdC-IR therapies of MMR defective cancers; MMR defective cancers are significant as they comprise ~10% of colorectal [<xref ref-type="bibr" rid="B9">9</xref>
], gastric [<xref ref-type="bibr" rid="B10">10</xref>
], pancreatic [<xref ref-type="bibr" rid="B11">11</xref>
], urinary [<xref ref-type="bibr" rid="B12">12</xref>
], gynecologic [<xref ref-type="bibr" rid="B13">13</xref>
,<xref ref-type="bibr" rid="B14">14</xref>
] and glioma [<xref ref-type="bibr" rid="B15">15</xref>
] cancers.</p>
<fig id="F1" position="float"><label>Figure 1</label>
<caption><p><bold>The dNTP supply system</bold>
. Thick lines are fluxes, thin solid lines are activations, thin dashed lines are inhibitions. Key enzymes are described in the text. An RNR <italic>s</italic>
-site mediated large positive feedback loop ATP → dCTP → dUMP → dTTP → dGTP → dATP terminates when dATP binds the R1 <italic>a</italic>
-site to inhibit all 4 RNR reductions. Models of enzymes of this system will eventually be useful in anticancer drug dose time course optimizations [<xref ref-type="bibr" rid="B1">1</xref>
].</p>
</caption>
<graphic xlink:href="1745-6150-4-50-1"></graphic>
</fig>
<p>The dNTP supply system is ideal for cancer systems biology research because, among cancer relevant processes, it is perhaps the best understood. This is important because, intuitively, the more understanding a mathematical model captures, the more likely it is to be more useful than a conceptual model. Thus, the dNTP supply system is well poised to be successfully controlled better with mathematical modeling than without, and because of this, this system could become a standard of success in systems biology; the basis of this argument is prior success in the use of mathematical models to improve the control of well understood systems such as power plants and airplanes.</p>
<p>RNR (NDP → dNDP) [<xref ref-type="bibr" rid="B16">16</xref>
] has two subunits, R1, and R2 or p53R2 [<xref ref-type="bibr" rid="B17">17</xref>
,<xref ref-type="bibr" rid="B18">18</xref>
]. On short time scales of seconds to minutes, RNR is controlled through two R1 regulatory sites, a selectivity (<italic>s</italic>
-) site that is somewhat analogous to a radio tuning control knob, and an activity (<italic>a</italic>
-) site that can be thought of as a volume control knob. Complicated positive and negative dNTP-mediated feedback loops (Fig. <xref ref-type="fig" rid="F1">1</xref>
) impinge upon these two sites to implement a sophisticated solution to a challenging dNTP pool balance regulation problem; if pool balance regulation performance varies across individuals and MMR performance also varies, individuals compromised in both systems may be predisposed to cancer [<xref ref-type="bibr" rid="B19">19</xref>
]. RNR functional complexity is mirrored by the combinatorial complexity of its R1 subunit: its catalytic site can be empty or filled with 1 of 4 NDP substrates, its <italic>s</italic>
-site can be empty or filled with ATP, dATP, dTTP or dGTP, its <italic>a</italic>
-site can be empty or bound by ATP for activation or dATP for inactivation, it may have an <italic>h</italic>
-site that can be empty or bound by ATP [<xref ref-type="bibr" rid="B20">20</xref>
], and all of this is folded by an R1 monomer-dimer-tetramer-hexamer equilibrium where R1 j-mers may also be bound by variable numbers of R2 (or p53R2) dimers. As a result, trillions of R1 complexes are possible if R1, R2, UDP, CDP, GDP, ADP, ATP, dATP, dTTP and dGTP are all present (in this case ~10<sup>2 </sup>
R1 monomers implies ~10<sup>12 </sup>
R1 hexamers) and the problem then is to determine which are needed in models to explain the data at hand. To appreciate the magnitude of the problem, if 10<sup>12 </sup>
complexes are possible, the number of possible complete dissociation constant models is 2 raised to the 10<sup>12 </sup>
(i.e. 1 followed by ~300 billion zeros), since each complex, and its corresponding complete dissociation constant K, can either be in the model (estimated) or out (set to infinity if the model hypothesizes that the concentration of the complex is approximately zero across all of the experimental conditions of the dataset). This huge number of models is even greater if, in addition to hypotheses that complete dissociation constants are infinite, hypotheses that binary dissociation constants equal each other are also considered. Though this problem is intractable for 10-reactants, 2-reactant solutions are feasible and may yield insights needed to enable 3-reactant solutions, and so on.</p>
<p>The R [<xref ref-type="bibr" rid="B21">21</xref>
] package Combinatorially Complex Equilibrium Model Selection (<monospace>ccems</monospace>
) [<xref ref-type="bibr" rid="B22">22</xref>
] is used here to automatically generate thousands of ATP induced R1 hexamerization models partitioned into two classes: those that include model terms of complexes with ATP occupied <italic>h</italic>
-sites (i.e. models that would support the existence of <italic>h</italic>
-sites if selected) and those that do not. Comparisons of the Akaike Information Criterion (AIC) [<xref ref-type="bibr" rid="B23">23</xref>
,<xref ref-type="bibr" rid="B24">24</xref>
] of these two classes of models were then used to assess the extent of <italic>h</italic>
-site evidence strictly in the ATP-induced R1 hexamerization dynamic light scattering (DLS) data in figure <xref ref-type="fig" rid="F1">1</xref>
of reference [<xref ref-type="bibr" rid="B20">20</xref>
]. No evidence was found. As discussed in Kashlan's review below, evidence for an <italic>h</italic>
-site may, however, exist under different experimental conditions.</p>
</sec>
<sec><title>A Simple Example</title>
<p>To introduce concepts of model space generation, standard models of competitive and non-competitive inhibition are derived below as instances of models in two systematically defined model spaces, one of spur graphs which focus on <italic>complete </italic>
dissociation constants and hypotheses that they are infinite, and one of grid graphs which focus on <italic>binary </italic>
dissociation constants and hypotheses that they are equal.</p>
<sec><title>Spur Graphs</title>
<p>Consider the enzyme-substrate-inhibitor (ESI) models of Fig. <xref ref-type="fig" rid="F2">2</xref>
. The full spur graph at the top of this figure is represented by the following total concentration constraint (TCC) system of coupled free concentration polynomials:<disp-formula id="bmcM1"><label>(1)</label>
<graphic xlink:href="1745-6150-4-50-i1.gif"></graphic>
</disp-formula>
</p>
<fig id="F2" position="float"><label>Figure 2</label>
<caption><p><bold>ESI model graphs</bold>
. The full spur graph at the top generates the seven models/graphs below it via hypotheses taken one at a time, two at a time, etc, that dissociation constants are infinite. The C-shaped grid graph is a data-fitting equivalent of the full spur graph. It generates the non-competitive inhibition model to its right where parallel edge binary K's are equal.</p>
</caption>
<graphic xlink:href="1745-6150-4-50-2"></graphic>
</fig>
<p>where [E<sub>T</sub>
], [S<sub>T</sub>
] and [I<sub>T</sub>
] are the known total concentrations, or system inputs, [E], [S] and [I] are model dependent free concentrations, which can be thought of as predicted latent system state variables, the dissociation constants <italic>K </italic>
are <italic>complete </italic>
dissociation constants, and implicit in these equations is the equilibrium concept of mass action, e.g. [ESI] = [E][S][I]/K<sub>ESI</sub>
; in spur graphs every complex is reached from free E in a single step and thus every term in the corresponding equations has only one <italic>K </italic>
parameter in its denominator, see [<xref ref-type="bibr" rid="B25">25</xref>
]. The model in Eq. (1) is called a full model because it is fully parameterized to the extent that no constraints have yet been placed on its parameters. Model spaces are then generated from full models by applying constraints to them. For example, the other spur graphs in Fig. <xref ref-type="fig" rid="F2">2</xref>
are obtained from the full spur graph by applying <italic>K </italic>
= ∞ constraints to remove nodes in the graph, and thus columns of <italic>K </italic>
terms in Eq. (1), one at a time, two at a time, and three at a time. Each resulting graph/model is then a hypothesis that some <italic>K</italic>
's are approximately infinitely large, or that concentrations of the deleted complexes are so small, throughout all of the experimental conditions of the dataset, that they can be approximated as zero, e.g. the competitive inhibition model, where substrate and inhibitor bind the same site and thus cannot bind simultaneously, hypothesizes that [ESI] = 0, or equivalently, that <italic>K</italic>
<sub>ESI </sub>
= ∞, and its equations are<disp-formula><graphic xlink:href="1745-6150-4-50-i2.gif"></graphic>
</disp-formula>
</p>
<p>Note that this model corresponds to the full spur graph less the ESI node/edge and that it can also be viewed as the pair (<italic>K</italic>
<sub>ESI </sub>
= ∞, Eq. 1).</p>
</sec>
<sec><title>Grid Graphs</title>
<p>Figure <xref ref-type="fig" rid="F2">2</xref>
also shows a C-shaped full grid graph with system equations<disp-formula id="bmcM2"><label>(2)</label>
<graphic xlink:href="1745-6150-4-50-i3.gif"></graphic>
</disp-formula>
</p>
<p>where underscores in subscripts indicate specific binary reactions. In grid graphs, because the dissociation constants are binary, equation terms that represent complexes of <italic>n </italic>
reactants have <italic>n </italic>
- 1 <italic>K </italic>
parameters in their denominators.</p>
<p>There is a one-to-one mapping between the <italic>K</italic>
's in grid graph system (2) and those in spur graph system (1), namely, <italic>K</italic>
<sub>ES </sub>
= <italic>K</italic>
<sub>E_S</sub>
, <italic>K</italic>
<sub>EI </sub>
= <italic>K</italic>
<sub>E_I</sub>
, and <italic>K</italic>
<sub>ESI </sub>
= <italic>K</italic>
<sub>E_I</sub>
<italic>K</italic>
<sub>EI_S </sub>
or <italic>K</italic>
<sub>EI_S </sub>
= <italic>K</italic>
<sub>ESI</sub>
/<italic>K</italic>
<sub>EI</sub>
. It follows then that (1) and (2) are data-fitting equivalents. The added value of full grid graph systems such as system (2) is their ability to spawn new hypotheses that cannot be generated by corresponding full spur graph systems such as (1). The additional hypotheses are binary <italic>K </italic>
equality hypotheses. In the example here there is only one such hypothesis/model, the non-competitive inhibition model <italic>K</italic>
<sub>E_S </sub>
= <italic>K</italic>
<sub>EI_S </sub>
where inhibitor binding has no detectable effect on substrate binding.</p>
</sec>
<sec><title>Solutions</title>
<p>System equations are solved as the steady state of a parent system of ordinary differential equations (ODEs) [<xref ref-type="bibr" rid="B25">25</xref>
], e.g. system (1) is solved by simulation of<disp-formula id="bmcM3"><label>(3)</label>
<graphic xlink:href="1745-6150-4-50-i4.gif"></graphic>
</disp-formula>
</p>
<p>to large τ where τ has nothing to do with real time and the state trajectory is merely algorithmic and thus not a biophysical path. Free concentration solutions then map to complex concentrations via mass action laws and these are then mapped to expected measurements, e.g. see Eq. (6) below and Eqs. (16-18) in [<xref ref-type="bibr" rid="B25">25</xref>
].</p>
</sec>
</sec>
</sec>
<sec><title>Results</title>
<sec><title>Limitations</title>
<p>The methods presented here are currently limited to biochemical systems where one central hub protein mediates all of the interactions and total concentrations of the reactants are approximately known exactly. It is assumed that the latter condition is adequately met in analyses of data derived from systems that were reconstituted from purified reactants. If the hub protein has more than one binding site for the same ligand, as R1 does for ATP, to interpret model terms, a specified sequence of site filling must be assumed. This assumption, made due to lack of a better option, may not hold. Automated model space generation is currently limited to two-reactant systems.</p>
</sec>
<sec><title>ATP-induced R1 Hexamerization Models</title>
<sec><title>Full model</title>
<p>To generate a space of ATP induced R1 hexamerization models, the first step is to pick a full model and the second step is to apply <italic>K </italic>
hypotheses to it [<xref ref-type="bibr" rid="B25">25</xref>
]. Full models that include <italic>s</italic>
, <italic>a </italic>
and <italic>h </italic>
ATP binding sites generate two classes of models: those that include at least one occupied <italic>h</italic>
-site complex and that thus support the existence of an <italic>h</italic>
-site, if selected, and those that do not, i.e. those that allege that all occupied <italic>h</italic>
-site complex concentrations are approximately zero and that thus support claims that the <italic>h</italic>
-site is not needed to explain the data, if selected. The full model below generates both of these model types.</p>
<p>Full models place an upper bound on the complexity and size of the model space and should thus be no more complicated than needed to answer the question of interest, e.g. models with a fourth ATP binding site should not be considered as there is no evidence that such a site exists. Lower bounds can also be placed on the simplicity of complexes, and this too can reduce the size of the model space. Thus, based on the crystal structure of yeast R1 dimers that shows that the <italic>s</italic>
-site is created at the R1 dimer interface [<xref ref-type="bibr" rid="B26">26</xref>
], and based on dTTP induced R1 dimerization being well represented by free reactants forming (R1)<sub>2</sub>
(dTTP)<sub>2 </sub>
directly [<xref ref-type="bibr" rid="B25">25</xref>
] (dTTP binds only at the <italic>s</italic>
-site), it is reasonable to assume that R1 oligomer <italic>s</italic>
-sites are always fully occupied (i.e. that oligomers cannot form without full <italic>s</italic>
-site occupancy) and that R1 monomer <italic>s </italic>
sites are always unoccupied (i.e. that the <italic>s</italic>
-site does not exist in R1 monomers). With these restrictions, denoting ATP and R1 by X and R respectively, and using<disp-formula id="bmcM4"><label>(4)</label>
<graphic xlink:href="1745-6150-4-50-i5.gif"></graphic>
</disp-formula>
</p>
<p>the full spur graph system equations are:<disp-formula id="bmcM5"><label>(5)</label>
<graphic xlink:href="1745-6150-4-50-i6.gif"></graphic>
</disp-formula>
</p>
<p>where, in each equation, first sum limits assume <italic>s</italic>
-sites cannot be bound in monomers and other sum limits assume <italic>s</italic>
-sites must be bound in oligomers; here X = ATP is used because a = dATP and A = ADP are being reserved for subsequent RNR models and R<sup>j</sup>
X<sup>i </sup>
is used instead of R<sub>j</sub>
X<sub>i </sub>
to stress connections to polynomials.</p>
</sec>
<sec><title>Interpretations</title>
<p>The hub protein monomer complex RX can be interpreted as X bound to either the <italic>a</italic>
- or <italic>h</italic>
-site. Because the <italic>a</italic>
-site is known to exist, <italic>a</italic>
-site binding will be assumed. RXX is then a monomer with both the <italic>a</italic>
- and <italic>h</italic>
-sites occupied. For oligomers, in addition to all of the <italic>s</italic>
-sites being prefilled, it is assumed that: R<sup>2</sup>
X<sup>2 </sup>
through R<sup>2</sup>
X<sup>4</sup>
, R<sup>4</sup>
X<sup>4 </sup>
through R<sup>4</sup>
X<sup>8</sup>
, and R<sup>6</sup>
X<sup>6 </sup>
through R<sup>6</sup>
X<sup>12</sup>
, have zero to full <italic>a</italic>
-site occupancies and no <italic>h</italic>
-site occupancies, and that R<sup>2</sup>
X<sup>5 </sup>
and R<sup>2</sup>
X<sup>6</sup>
, R<sup>4</sup>
X<sup>9 </sup>
through R<sup>4</sup>
X<sup>12</sup>
, and R<sup>6</sup>
X<sup>13 </sup>
through R<sup>6</sup>
X<sup>18</sup>
, have partial to full <italic>h</italic>
-site occupancies in addition to completely filled <italic>a</italic>
-sites (and <italic>s</italic>
-sites). Model inferences will be based on these interpretations.</p>
</sec>
<sec><title>Output linkage</title>
<p>The fitted output measurements are mass-weighted average protein masses<disp-formula id="bmcM6"><label>(6)</label>
<graphic xlink:href="1745-6150-4-50-i7.gif"></graphic>
</disp-formula>
</p>
<p>where [<italic>R</italic>
<sub>jT</sub>
] is the total <italic>j</italic>
-mer concentration (i.e. the sum of the concentrations of all <italic>j</italic>
-mer hub protein complexes),<disp-formula><graphic xlink:href="1745-6150-4-50-i8.gif"></graphic>
</disp-formula>
</p>
<p>is the total hub protein concentration, M<sub>1 </sub>
is the mass of a monomer (90 kDa for the R1 subunit of RNR), <italic>ε </italic>
is noise with constant variance and zero mean, and the <italic>j</italic>
<sup>2 </sup>
in the numerator includes one factor of <italic>j </italic>
because the mass of a <italic>j</italic>
-mer is <italic>j </italic>
times that of a monomer, and another factor of <italic>j </italic>
because light scattering is proportional to mass; ligand masses are treated as negligible relative to protein masses.</p>
</sec>
<sec><title>Implementation</title>
<p>Nonlinear least squares was used to fit the models. The fitted models were then rank ordered by their AICs [<xref ref-type="bibr" rid="B24">24</xref>
]. Because the number of data points <italic>N </italic>
is small at 15, the small sample size corrected version of the AIC was used: AIC = 2*<italic>P </italic>
+ 2*<italic>P</italic>
(<italic>P</italic>
+1)/(<italic>N</italic>
-<italic>P</italic>
-1) + <italic>N</italic>
*log (2π) + <italic>N</italic>
*log (SSE/<italic>N</italic>
) + <italic>N </italic>
where <italic>P </italic>
is the number of estimated parameters (including the variance) and SSE is the sum of squared errors [<xref ref-type="bibr" rid="B24">24</xref>
]. In parameter optimizations (i.e. SSE minimizations, see Methods) the initial complete dissociation constant values were 100 μM raised to the sum of the powers of the numerators in Eq. (4) minus one, i.e. <italic>j </italic>
+ <italic>i </italic>
- 1. This was critical, as it increased the number of models that converged from roughly 10% (when 1 μM was used uniformly) to nearly 100%. Models were fitted in parallel in a load balanced manner using R [<xref ref-type="bibr" rid="B21">21</xref>
]; the R package <monospace>ccems</monospace>
uses the R package <monospace>snow</monospace>
(small network of workstations) to accomplish this. R scripts that were used to produce the results in this and the accompanying paper are available as examples in the <monospace>papers</monospace>
directory of <monospace>ccems</monospace>
.</p>
</sec>
<sec><title>Spur graphs</title>
<p>The number of complexes represented in Eq. (5) is 2 + 5 + 9 + 13 = 29 and this implies that the number of spur models is 2<sup>29 </sup>
= ~500 million. Relevant here, however, is the number of 1-, 2- and 3-parameter models. There are 29 single-edge models, 406 (29 choose 2) two-edge models, and 3654 (29 choose 3) three-edge models. The lowest AIC of the 3-parameter models (144.4) was higher than the lowest AIC of the 2-parameter models (142.7), so 4-parameter models were not fitted; as parameter numbers increase AICs typically first decrease as SSEs decrease, but then rise and continue to rise due to over-parameterization.</p>
</sec>
<sec><title>Grid graphs</title>
<p>Binary dissociation constants that are alleged equal to one another must be defined on a per site basis in terms of <italic>k</italic>
<sub>off</sub>
/<italic>k</italic>
<sub>on</sub>
. For example, using the equilibrium property that net fluxes between any two complexes must vanish, for <italic>a</italic>
-sites in tetramers, the on and off fluxes of the 1<sup>st </sup>
ligand yield<disp-formula><graphic xlink:href="1745-6150-4-50-i9.gif"></graphic>
</disp-formula>
</p>
<p>where factors of 1 and 4 arise because there is 1 occupied site for the dissociation reaction and 4 unoccupied sites for the association reaction, respectively. Similarly, for second, third and fourth ligand bindings to a tetramer <italic>a</italic>
-site, the per site binary dissociation constants <italic>K </italic>
are:<disp-formula><graphic xlink:href="1745-6150-4-50-i10.gif"></graphic>
</disp-formula>
</p>
<p>These binary <italic>K </italic>
can be hypothesized to equal each other and same-site per-site binary <italic>K </italic>
of other oligomers such as<disp-formula><graphic xlink:href="1745-6150-4-50-i11.gif"></graphic>
</disp-formula>
</p>
<p>Similar arguments apply to <italic>h</italic>
-sites, with <italic>X</italic>
<sup><italic>i</italic>
+<italic>j </italic>
</sup>
replacing <italic>X</italic>
<sup><italic>i </italic>
</sup>
in <italic>j</italic>
-mers.</p>
<p>To introduce a few additional concepts in their simplest forms, the next two paragraphs consider dTTP induced R1 dimerization. That not all <italic>K </italic>
equality hypotheses take the form <italic>K </italic>
= <italic>K</italic>
' is seen in the penultimate column of the n-shaped graph pairs in the top half of Fig. <xref ref-type="fig" rid="F3">3</xref>
where K<sup>2</sup>
<sub>R_t </sub>
= [R]<sup>2</sup>
[t]<sup>2</sup>
/[Rt]<sup>2 </sup>
= K<sub>RR_t</sub>
K<sub>RRt_t </sub>
= [RR][t]<sup>2</sup>
/[RRtt]; here t denotes dTTP. This graph pair can, however, be restated as the equivalent <italic>K </italic>
= <italic>K</italic>
' E-shaped graph pair shown in the corresponding positions below it. Though earlier work focused on E-shaped graphs [<xref ref-type="bibr" rid="B25">25</xref>
], based on the paragraph that follows, it suffices to consider only n-shaped graphs and natural extensions thereof (e.g. see Fig. <xref ref-type="fig" rid="F4">4</xref>
).</p>
<fig id="F3" position="float"><label>Figure 3</label>
<caption><p><bold>dTTP-induced R1 dimerization K equality models</bold>
. The n-shaped graphs are equal to their corresponding E-shaped graphs below them. The rightmost three columns are very unlikely (see text). R = R1, t = dTTP and edges marked = or -- are alleged equal.</p>
</caption>
<graphic xlink:href="1745-6150-4-50-3"></graphic>
</fig>
<fig id="F4" position="float"><label>Figure 4</label>
<caption><p><bold><italic>K </italic>
equality RX model space</bold>
. The graphs shown are the same for both <italic>a</italic>
- and <italic>h</italic>
-site models. The top two rows have independent threads and the bottom two rows have at least two threads that have equal binary <italic>K </italic>
values indicated by =, --, or x. Bridge edges in the horizontal bars of each graph (i.e. the curtain rods) are spur edges from the hub, rather than binary K. Left to right, threads on curtain rods correspond to monomers, dimers, tetramers and hexamers.</p>
</caption>
<graphic xlink:href="1745-6150-4-50-4"></graphic>
</fig>
<p>The first two columns of <italic>K </italic>
equality models/hypotheses in Fig. <xref ref-type="fig" rid="F3">3</xref>
are plausible, as a protein could be so rigid that a ligand binding site is unchanged with respect to binding affinity regardless of other bindings (first column) and this could be true within <italic>j</italic>
-mers but not between them (second column). The third column hypotheses are less likely, however, as they claim that binding of R (to R), which is massive relative to ligand and thus more likely to cause alterations upon binding, causes no change in ligand site affinities, yet, binding of the first small ligand to the dimer alters the dimer structure enough to change ligand binding at the second site. Continuing, fourth column hypotheses are even less likely, as they claim that the first ligand binds dimer differently than monomer, yet, after it binds, by chance, the second ligand binding energy exactly equals the amount needed for the product of these <italic>K</italic>
's to equal the square of the monomer ligand <italic>K</italic>
. Equivalent E-shaped graphs (same column) support this claim of unlikelihood, as they claim that although dimerization energies are different between R + R and Rt + R, they somehow return to the R + R value when both reactants are Rt. Finally, in the fifth column, it would be remarkable if ligand binding to free dimer differs from ligand binding to free monomer, yet somehow, binding of the first ligand returns the unoccupied dimer subunit to a state indistinguishable from that of the free monomer (with respect to its ligand binding affinity). The third and fifth columns can be interpreted in terms of rigid asymmetric dimers where one subunit holds its monomer shape and the other has a different shape with either tighter (fifth column) or weaker (third column) ligand binding. From this perspective, it is very unlikely that all of the dimerization induced shape changes (deformation energies) would fall strictly onto one of two identical subunits. Thus, the grid model space used here will only include <italic>K </italic>
equality hypotheses that are analogous to the first 2 columns in Fig. <xref ref-type="fig" rid="F3">3</xref>
, and it suffices to consider n-shaped graphs.</p>
<p>Within a site type, binary <italic>K</italic>
's of <italic>j</italic>
-mers can be depicted as threads hanging from a curtain rod as shown in Fig. <xref ref-type="fig" rid="F4">4</xref>
. In the accompanying paper, binary <italic>K </italic>
equalities in contiguous stretches within threads are considered. Here, each thread is homogeneous in its binary <italic>K </italic>
value (i.e. only full thread length contiguous stretches are considered) and <italic>K </italic>
equality models are instead generated by considering thread <italic>K </italic>
values as independent of other threads (top 2 rows in Fig. <xref ref-type="fig" rid="F4">4</xref>
), infinite (graphs with missing threads), or equal to those of other threads of the same site type (bottom 3 rows in Fig. <xref ref-type="fig" rid="F4">4</xref>
) within contiguous stretches of threads, the idea being that if one protein oligomerization step alters a ligand <italic>K</italic>
, it is unlikely that an additional step would return it to one of its previous values.</p>
<p>Since R binds R to form R<sup>2 </sup>
via one protein surface, and since it is likely that R<sup>2 </sup>
binds R<sup>2 </sup>
and R<sup>4 </sup>
using two different protein surfaces (or patches thereof), no hypotheses of <italic>K </italic>
equivalence will be considered between the saturated <italic>s</italic>
-site complexes R<sup>2</sup>
X<sup>2</sup>
, R<sup>4</sup>
X<sup>4 </sup>
and R<sup>6</sup>
X<sup>6 </sup>
(i.e. complexes in the curtain rods in Fig. <xref ref-type="fig" rid="F4">4</xref>
). Thus, all of the <italic>K </italic>
equality hypotheses explored will be with respect to ligand binding site constants in threads.</p>
<p>The binary <italic>K </italic>
equality model space of interest here is the product of a space of <italic>a</italic>
-site models and a completely analogous space of <italic>h</italic>
-site models. The 32 models shown in Fig. <xref ref-type="fig" rid="F4">4</xref>
thus imply a <italic>K </italic>
equality space of 1024 models. If thread head nodes within curtain rods are allowed to remain in models where all other nodes in the same thread have infinite <italic>K</italic>
, the number of models increases: models missing hexamer threads split into two models (there are 8 of these in Fig. <xref ref-type="fig" rid="F4">4</xref>
) and models missing both tetramer and hexamer threads (there are 3 of these in Fig. <xref ref-type="fig" rid="F4">4</xref>
) split into four models. The total number of grid models is then (32 + 8 + 9)<sup>2 </sup>
= 49<sup>2 </sup>
= 2401.</p>
</sec>
<sec><title>Models that contain hexamer terms</title>
<p>Since external data [<xref ref-type="bibr" rid="B27">27</xref>
] confirms ATP induced R1 hexamerization, the model space was reduced to only models that contain at least one hexamer term. This reduced the number of grid models with 2 and 3 parameters to 2 and 15 (from 7 and 36) and the number of spur models with 1, 2 or 3 parameters to 13, 286 and 3094 (from 29, 406 and 3654), i.e. the number of models is now 17 + 3393 = 3410.</p>
</sec>
<sec><title>Competitive models</title>
<p>Of 3410 ATP-induced R1 hexamerization models automatically fitted to the DLS data [<xref ref-type="bibr" rid="B20">20</xref>
], four failed to converge (these all involved X<sup>17 </sup>
and X<sup>18</sup>
), 966 (Fig. <xref ref-type="fig" rid="F5">5B</xref>
) converged but had singular Hessians (12 and 954 of these were 2- and 3-parameter models and all had infinite upper confidence limits in all parameters) and of the remaining 2440 fits (Fig. <xref ref-type="fig" rid="F5">5A</xref>
), 1200 had no infinite upper bounds, 680 had 1 of 3, 91 had 1 of 2, and 369 had 2 of 3. Models with no infinite upper bounds comprised 64%, 6% and <1% of the lowest to highest AIC clusters shown in Fig. <xref ref-type="fig" rid="F5">5A</xref>
. To purge the space of problematic models (e.g. those that were either incorrect or sensitive to initial parameter values), the space was reduced to models with SSEs that were less than twice the minimum SSE of the fitted models. AICs of the resulting 2088 models are shown in Figs. <xref ref-type="fig" rid="F5">5C</xref>
and <xref ref-type="fig" rid="F5">5D</xref>
.</p>
<fig id="F5" position="float"><label>Figure 5</label>
<caption><p><bold>AIC model densities versus AIC model rank</bold>
. <bold>A) </bold>
and <bold>C) </bold>
show models that had non-singular Hessians (matrices of objective function second derivatives) at their optimums; <bold>B) </bold>
and <bold>D) </bold>
are models with singular (determinant = 0) Hessians, i.e. models that converged onto likelihood surfaces that were flat in one direction. <bold>C) </bold>
and <bold>D) </bold>
show 1613 and 475 models with SSEs < twice the minimum SSE.</p>
</caption>
<graphic xlink:href="1745-6150-4-50-5"></graphic>
</fig>
</sec>
<sec><title>h-site existence</title>
<p>Fig. <xref ref-type="fig" rid="F6">6A</xref>
shows normalized AIC densities of the models in Figs. <xref ref-type="fig" rid="F5">5C</xref>
and <xref ref-type="fig" rid="F5">5D</xref>
partitioned into two groups, those that do not represent any complexes that have occupied <italic>h</italic>
-sites (508 models) and those that do (1580 models). That these densities differ is apparent by inspection, a Kolmogorov-Smirnov p < 10<sup>-16 </sup>
and by 28 of the top 30 models not including an <italic>h</italic>
-site term, i.e. 28/30 > 508/2088 with p < 1 × 10<sup>-15</sup>
. The densities were then decomposed into models with <3 parameters (6B), 3 parameters (6C) and singular Hessians (6D). In each case the same conclusion held, the DLS data did not support the existence of an <italic>h</italic>
-site.</p>
<fig id="F6" position="float"><label>Figure 6</label>
<caption><p><bold>Normalized AIC densities of models with SSEs less than twice the minimum</bold>
. <bold>A) </bold>
The models of Figs. 5C and 5D. There are 508 (1580) models without (with) an occupied <italic>h</italic>
-site. <bold>B) </bold>
The 1- and 2-parameter models of A). <bold>C) </bold>
The 3-parameter models of A). <bold>D) </bold>
The models in A) that have singular Hessians (i.e. the 475 models in Fig. 5D). In all cases a difference in <italic>h</italic>
-site hypothesis densities is supported by a two-sample Kolmogorov-Smirnov test, P < 10<sup>-15 </sup>
(A, C, D) and P < 10<sup>-5 </sup>
(B).</p>
</caption>
<graphic xlink:href="1745-6150-4-50-6"></graphic>
</fig>
<p>Fig. <xref ref-type="fig" rid="F7">7</xref>
shows the fits of only the 1-parameter models. The model R<sup>6</sup>
X<sup>8 </sup>
fits the data the best and is immediately flanked by R<sup>6</sup>
X<sup>7 </sup>
(less steep) and R<sup>6</sup>
X<sup>9 </sup>
(steeper) which also fit reasonably well, but beyond these, the fits become noticeably poorer; since these models all have the same number of parameters, AICs in the legend reflect SSEs and the gap between the 3<sup>rd </sup>
and 4<sup>th </sup>
model thus reflects a lack of fit. That R<sup>6</sup>
X<sup>6 </sup>
is a poor fit supports the assumption that <italic>s</italic>
-sites are pre-filled, and that R<sup>6</sup>
X<sup>10 </sup>
is a poor fit supports a hypothesis that only 1 to 3 hexamer <italic>a</italic>
-sites are filled. That R<sup>6</sup>
X<sup>11 </sup>
and higher models provide worse and worse fits with increasing numbers of bound ATP supports the Fig. <xref ref-type="fig" rid="F6">6</xref>
conclusion that this dataset does not support the existence of an <italic>h</italic>
-site.</p>
<fig id="F7" position="float"><label>Figure 7</label>
<caption><p><bold>Fits of the 1-parameter models</bold>
. The legend in the plot indicates the model order ranked by AICs (values shown). The top 5 models are indicated by thicker lines. Beyond 5 or more <italic>a</italic>
-sites occupied by ATP and with increasing numbers of <italic>h</italic>
-sites occupied, the fits become worse and worse.</p>
</caption>
<graphic xlink:href="1745-6150-4-50-7"></graphic>
</fig>
<p>Fits of the top 3 models are shown in Fig. <xref ref-type="fig" rid="F8">8</xref>
and their parameter estimates are given in Table <xref ref-type="table" rid="T1">1</xref>
. The second of these is a 2-parameter model that differs from the other two in that it uses its second parameter to avoid its obligation to reach an average mass of 540 kDa in the limit of large [ATP]; when this model was extrapolated to [ATP] = 1 M (Fig. <xref ref-type="fig" rid="F9">9</xref>
) the predicted average mass in this limit was 180 kDa. It was conjectured then that complexes with the highest ratios of ATP bound per R1 completely dominate the distribution in the limit of high ligand concentrations and that only in cases of maximum ratio ties does a balance result (i.e. that the system's objective in this limit is to partition as much ATP as possible away from its free form and into a bound form). Representative models of both a balance and of hexamer dominance support this conjecture (Fig. <xref ref-type="fig" rid="F9">9</xref>
). If it is asserted then that the ATP per R1 ratio cannot decrease as higher R1 oligomers are formed, the space of 2088 models reduces to 1420 models, but the calculations of Fig. <xref ref-type="fig" rid="F7">7</xref>
still yield the same conclusion, i.e. no <italic>h</italic>
-site existence (Fig. <xref ref-type="fig" rid="F10">10</xref>
). If the ratio is forced to strictly increase with oligomerization, the number of models is 1287, and again, the same conclusion holds (plots not shown).</p>
<table-wrap id="T1" position="float"><label>Table 1</label>
<caption><p>The top 3 models (lowest AIC) of the RX model space.</p>
</caption>
<table frame="hsides" rules="groups"><thead><tr><th align="left">Model</th>
<th align="left">Parameter</th>
<th align="left">Initial Value</th>
<th align="left">Optimal Value</th>
<th align="left">Confidence Interval</th>
</tr>
</thead>
<tbody><tr><td align="left"><bold>1 IIIIIIIIIIIJIIIIIIIIIIIIIIIII</bold>
</td>
<td align="left">R6X8</td>
<td align="left">100.000^13</td>
<td align="left">63.101^13</td>
<td align="left">(59.878^13, 66.175^13)</td>
</tr>
<tr><td colspan="5"><hr></hr>
</td>
</tr>
<tr><td></td>
<td align="left">SSE</td>
<td align="left">840487.830</td>
<td align="left">7726.693</td>
<td></td>
</tr>
<tr><td colspan="5"><hr></hr>
</td>
</tr>
<tr><td></td>
<td align="left">AIC</td>
<td align="left">211.573</td>
<td align="left">141.234</td>
<td></td>
</tr>
<tr><td colspan="5"><hr></hr>
</td>
</tr>
<tr><td></td>
<td align="left">CPU</td>
<td align="left">0.000</td>
<td align="left">3.465</td>
<td align="left">fit succeeded</td>
</tr>
<tr><td colspan="5"><hr></hr>
</td>
</tr>
<tr><td align="left"><bold>2 IIIIIJIIIIIJIIIIIIIIIIIIIIIII</bold>
</td>
<td align="left">R2X4</td>
<td align="left">100.000^5</td>
<td align="left">432.997^5</td>
<td align="left">(311.064^5, 601.845^5)</td>
</tr>
<tr><td colspan="5"><hr></hr>
</td>
</tr>
<tr><td></td>
<td align="left">R6X8</td>
<td align="left">100.000^13</td>
<td align="left">62.796^13</td>
<td align="left">(59.878^13, 66.175^13)</td>
</tr>
<tr><td colspan="5"><hr></hr>
</td>
</tr>
<tr><td></td>
<td align="left">SSE</td>
<td align="left">778605.476</td>
<td align="left">6873.411</td>
<td></td>
</tr>
<tr><td colspan="5"><hr></hr>
</td>
</tr>
<tr><td></td>
<td align="left">AIC</td>
<td align="left">213.608</td>
<td align="left">142.660</td>
<td></td>
</tr>
<tr><td colspan="5"><hr></hr>
</td>
</tr>
<tr><td></td>
<td align="left">CPU</td>
<td align="left">0.000</td>
<td align="left">5.302</td>
<td align="left">fit succeeded</td>
</tr>
<tr><td colspan="5"><hr></hr>
</td>
</tr>
<tr><td align="left"><bold>3 IIIIIIIIIIIIJIIIIIIIIIIIIIIII</bold>
</td>
<td align="left">R6X9</td>
<td align="left">100.000^14</td>
<td align="left">70.367^14</td>
<td align="left">(66.686^14, 74.228^14)</td>
</tr>
<tr><td colspan="5"><hr></hr>
</td>
</tr>
<tr><td></td>
<td align="left">SSE</td>
<td align="left">841284.881</td>
<td align="left">8624.411</td>
<td></td>
</tr>
<tr><td colspan="5"><hr></hr>
</td>
</tr>
<tr><td></td>
<td align="left">AIC</td>
<td align="left">211.588</td>
<td align="left">142.883</td>
<td></td>
</tr>
<tr><td colspan="5"><hr></hr>
</td>
</tr>
<tr><td></td>
<td align="left">CPU</td>
<td align="left">0.000</td>
<td align="left">2.903</td>
<td align="left">fit succeeded</td>
</tr>
</tbody>
</table>
</table-wrap>
<fig id="F8" position="float"><label>Figure 8</label>
<caption><p><bold>Top 3 fits to the data of Kashlan <italic>et al</italic>
</bold>
. [<xref ref-type="bibr" rid="B20">20</xref>
]. <bold>A) </bold>
Fits of the top 3 models. The second model uses its dimer term to capture a slight downturn in average mass at high [ATP] (see Fig. 9), consistent with its geometric mean binding constant being greater than that of the hexamer term in Table 1. Not shown in this plot is the point (0 μM, 90 kDa) which all models fit perfectly if M<sub>1 </sub>
= 90 is fixed (as it is here) rather than estimated. <bold>B) </bold>
Residuals of the top model R6X8. The reduced variance and positive mean of the first 4 residuals may be due to bias arising from prior knowledge that the monomer is 90 kDa and thus too prior knowledge that the average mass must increase from 90 kDa. Non-weighted least squares gives less weight to these 4 points which, coincidentally, is advantageous given these suspicions. PlotDigitizer (Methods) was used to obtain the data from a pdf of the original paper, and as this step involves human intervention, it too may have introduced some bias and random error.</p>
</caption>
<graphic xlink:href="1745-6150-4-50-8"></graphic>
</fig>
<fig id="F9" position="float"><label>Figure 9</label>
<caption><p><bold>Model limits at large ligand concentrations</bold>
. In the second model in Fig. 8 and Table 1 (solid line here) the dimer term R2X4 causes a below expectation peak (510 instead of 540) at high [ATP]. In the limit of very high [ATP], this model is dominated by the R2X4 term (average mass approaches 180 kDa) because this term partitions more X into a bound state with R than the hexamer R6X8, i.e. 4/2 = 2 > 8/6 = 1.3. These ratios are both 1.5 in the model R2X3.R6X9 (dashed line) which has a balanced population in this limit (with a limiting average mass of 523 kDa). These ratios are 1.5 and 2 in the model R2X3.R6X12 (dotted) which yields pure hexamers (average mass = 540 kDa) in the limit of infinite ATP.</p>
</caption>
<graphic xlink:href="1745-6150-4-50-9"></graphic>
</fig>
<fig id="F10" position="float"><label>Figure 10</label>
<caption><p><bold>Normalized model number densities of AICs of models with SSEs less than twice the minimum SSE and monotonic non-decreasing ratios of ATP per R1 as oligomer sizes increase</bold>
. <bold>A) </bold>
The complete set of such models. <bold>B) </bold>
The 1- and 2-parameter models. <bold>C) </bold>
The 3-parameter models. <bold>D) </bold>
Models with singular Hessians. In all cases a difference in <italic>h</italic>
-site hypothesis densities is supported by a two-sample Kolmogorov-Smirnov test, P < 10<sup>-15 </sup>
(A, C, D) and P < 2 × 10<sup>-5 </sup>
(B). Compare to Fig. 6.</p>
</caption>
<graphic xlink:href="1745-6150-4-50-10"></graphic>
</fig>
</sec>
<sec><title>Unoccupied a-sites</title>
<p>The space of 2088 models contains 99 models that do not have any terms with ratios of ATP bound per R1 > 1.5 (i.e. models consistent with ≤ 50% <italic>a</italic>
-site filling in oligomers). Of the top 30 models, however, there were 14 such models, which is significant, 14/30 > 99/2088 with p < 3 × 10<sup>-16</sup>
. In the reduced spaces of 1420 and 1287 models (of the previous paragraph), the proportions are 15/30 > 75/1420 and 13/30 > 54/1287, which also yield p < 3 × 10<sup>-16</sup>
. Thus, this dataset does not support the existence of hexamers with >3 <italic>a</italic>
-sites occupied by ATP. It should be noted that this statement implies a lack of <italic>h</italic>
-site evidence if <italic>a</italic>
-sites fill before <italic>h</italic>
-sites.</p>
<p>In the spaces of 1420 and 1287 models, the top 2 models are R<sup>6</sup>
X<sup>8 </sup>
and R<sup>6</sup>
X<sup>9</sup>
. If the 13 single edge spur models are fitted with M<sub>1 </sub>
of Eq. (6) estimated, a probability <italic>p </italic>
that R1 is capable of oligomerizing estimated (in Eq. 5 <italic>p </italic>
would then multiply [R<sub>T</sub>
]), or both, the top 3 models are R<sup>6</sup>
X<sup>9</sup>
, R<sup>6</sup>
X<sup>8 </sup>
and R<sup>6</sup>
X<sup>10 </sup>
(see Fig. <xref ref-type="fig" rid="F11">11</xref>
), i.e. the main claim is still supported but subsequent studies may find that a 4<sup>th </sup>
hexamer <italic>a</italic>
-site can also be filled by ATP.</p>
<fig id="F11" position="float"><label>Figure 11</label>
<caption><p><bold>The 13 single-edge spur models of Fig. 7(A) with <italic>p </italic>
(B), M<sub>1 </sub>
(C) or both (D) estimated</bold>
. The plots show that R<sup>6</sup>
X<sup>9 </sup>
should perhaps be trusted more than R<sup>6</sup>
X<sup>8</sup>
. AICs in the B-D legends suggest that R<sup>6</sup>
X<sup>10 </sup>
(which allows a 4<sup>th </sup>
filled a-site) is more likely than R<sup>6</sup>
X<sup>7</sup>
.</p>
</caption>
<graphic xlink:href="1745-6150-4-50-11"></graphic>
</fig>
</sec>
</sec>
</sec>
<sec><title>Discussion</title>
<p>No terms higher than R<sup>6</sup>
X<sup>9 </sup>
were needed to explain the ATP induced R1 hexamerization data found in figure <xref ref-type="fig" rid="F1">1</xref>
of reference [<xref ref-type="bibr" rid="B20">20</xref>
]. If <italic>s</italic>
-sites fill before <italic>a</italic>
-sites, this implies that ~1/2 of the hexamer <italic>a</italic>
-sites are not bound by ATP under the experimental conditions of this dataset. If <italic>h</italic>
-sites fill after <italic>a</italic>
-sites, this also implies that <italic>h</italic>
-sites need not exist to explain this dataset. Since the <italic>s</italic>
-site is at the dimer interface in yeast [<xref ref-type="bibr" rid="B26">26</xref>
], and since it is reasonable that hexamers form as trimers of dimers, it is likely that <italic>s</italic>
-sites do fill first.</p>
<p>If it is true that <italic>s</italic>
-sites are always bound in oligomers and never bound in monomers, dNTP access to hexamer <italic>s</italic>
-sites, as is needed for RNR control, implies that either the monomer-dimer-tetramer-hexamer equilibrium is rapid enough that changes in ligand bound to the <italic>s</italic>
-site can occur via the monomer-dimer equilibrium, or perhaps the hexamer stabilizes internal dimers enough that hexamer <italic>s</italic>
-sites can vacate without hexamer decomposition. The latter case would complicate the analysis as the term R<sup>6</sup>
X<sup>9 </sup>
for example might then describe more than 3 filled <italic>a</italic>
-sites.</p>
<p>Regarding <italic>a</italic>
- before <italic>h</italic>
-site filling, since <italic>a</italic>
-sites are known to exist [<xref ref-type="bibr" rid="B16">16</xref>
] and <italic>h</italic>
-sites are in question, this is a reasonable default. The alternative, to assume <italic>h</italic>
-site existence and instead challenge <italic>a</italic>
-site existence, is much less reasonable.</p>
<p>The most important short time constant (i.e. allosteric) feedback control of RNR is via dATP. This statement is based on ATP being too broadly used in cells for its level to be manipulated to control dNTP supply, and dTTP and dGTP being only selectivity controllers while when dATP controls selectivity, it also closes a large positive feedback loop that threads through dTTP and dGTP in series, see Fig <xref ref-type="fig" rid="F1">1</xref>
). This may help dNTP pools fill uniformly and rapidly at the onset of S-phase. Once the dNTP pools are filled, dATP also has the responsibility of shutting off its <italic>s</italic>
-site mediated positive feedback loop through <italic>a</italic>
-site mediated inhibition of all four reaction rates (note that this argues in favor of dATP binding the <italic>s</italic>
-site more tightly than the <italic>a</italic>
-site). This picture suggests that prompt <italic>a</italic>
-site mediated inhibition in response to changes in [dATP] is important since without a rapid response, the <italic>s</italic>
-site mediated dATP positive feedback loop may cause [dATP] overshoot. Given that [dATP] << [ATP] implies that dATP collisions with <italic>a</italic>
-sites are much rarer than ATP <italic>a</italic>
-site collisions, when they do occur, it would help the circuit respond rapidly if at least half of the time the site was empty and thus ready to be filled. This leads to interesting speculations: R1 hexamers may have two types of <italic>a</italic>
-sites, one for ATP and one for dATP, and beyond ligand differences, with this view (Fig. <xref ref-type="fig" rid="F12">12</xref>
) dATP inhibition versus ATP activation could in part be due to differences in binding pockets of the two types of <italic>a</italic>
-sites. Indeed, this may be the reason that R1 hexamers exist.</p>
<fig id="F12" position="float"><label>Figure 12</label>
<caption><p><bold>R1 hexamer model</bold>
. R1 hexamer formation could result in the creation of two types of <italic>a</italic>
-sites.</p>
</caption>
<graphic xlink:href="1745-6150-4-50-12"></graphic>
</fig>
<p>The approach used selects model terms (and thus parameters) based on how needed they are to explain the data analyzed. If, in solution, hexamers rarely have more than 3 ATPs bound to their <italic>a</italic>
-sites, no parameters are allocated to complexes with higher numbers of bound ATP. The analysis presented does not claim that ≥ 3 <italic>a</italic>
-sites will remain unoccupied under R1 crystallization conditions that may differ greatly from those used to generate the data analyzed here.</p>
<p>For the data analyzed, [R<sub>T</sub>
] is 7 μM (i.e. yielding up to 21 μM of ATP binding sites if <italic>h</italic>
-sites exist) and the minimum [X<sub>T</sub>
] is 46 μM, so the approximation [X] = [X<sub>T</sub>
] would not have been valid for this lowest [ATP] data point. The value of such approximations is less with ≥ 2 oligomerization states than with one (which has an analytic solution, see accompanying paper), as the univariate polynomial that results still requires a numerical solution (e.g. [X] = [X<sub>T</sub>
] in Eq. (5) yields a 6<sup>th </sup>
order polynomial in free [R]), but univariate polynomials are solved much faster than multivariate polynomials, e.g. using ODEs as in Eq. (3), so the computational savings are worthwhile if the approximation is tolerable. For the data used here, as [X] = [X<sub>T</sub>
] is valid for most of the data points, this approximation caused deviations of only 1% in the parameter estimates of Table <xref ref-type="table" rid="T1">1</xref>
but it gave a 30-fold increase in computation speed.</p>
<p>Since there were 13, 286 and 3094 <italic>K </italic>
infinity spur graphs with 1, 2 or 3 parameters, compared to 0, 2 and 15 <italic>K </italic>
equality grid graphs, and since models with few parameters have an AIC advantage when dataset sizes are modest, it is not surprising that with 15 data points, the top models were all spur models. In the future, as automation affords richer datasets, grid graphs may become more competitive. Thus, though the grid graph enumeration efforts expended in this paper did not pay immediate dividends, they may in the future.</p>
<p>Contiguous stretches of equal binary <italic>K </italic>
parameters within threads were not explored because binary <italic>K </italic>
models were already non-competitive due to over-parameterization, and because additional ATP ligands on a <italic>j</italic>
-mer would not have changed DLS masses detectably, so <italic>K </italic>
cooperativity within threads would not have been detectable.</p>
<p>Microfluidic chip technology [<xref ref-type="bibr" rid="B28">28</xref>
-<xref ref-type="bibr" rid="B31">31</xref>
] will eventually enable 5-dimensional RNR studies where [ATP] and [NDP]s are fixed to <italic>in vivo </italic>
levels and [R1], [R2], [dATP], [dTTP] and [dGTP] vary across ranges centered about physiological operating points. If <monospace>ccems</monospace>
can automatically analyze new RNR data as it arrives, it could find uses in sequential experimental designs [<xref ref-type="bibr" rid="B32">32</xref>
] where the chip conditions of the next measurements are determined in real time to implement efficient 5-D sampling strategies.</p>
<p>As protein expression and purification core facilities become more common, reconstituted network analyses where alleged protein-protein interactions are mathematically characterized for applications in systems biology [<xref ref-type="bibr" rid="B33">33</xref>
,<xref ref-type="bibr" rid="B34">34</xref>
] will eventually also become more common. It is anticipated here that many of these interactions will be combinatorially complex and that <monospace>ccems</monospace>
will then find broader uses.</p>
</sec>
<sec><title>Conclusion</title>
<p>No terms higher than R<sup>6</sup>
X<sup>9 </sup>
were needed to explain the ATP induced R1 hexamerization data found in figure <xref ref-type="fig" rid="F1">1</xref>
of reference [<xref ref-type="bibr" rid="B20">20</xref>
]. This suggests that under the experimental conditions of this dataset, ~1/2 of the hexamer <italic>a</italic>
-sites are not bound by ATP, and that if <italic>a</italic>
-sites fill before <italic>h</italic>
-sites, that <italic>h</italic>
-sites need not exist to explain this dataset.</p>
<p>The R package <monospace>ccems</monospace>
currently solves 2-reactant problems where total reactant concentrations are known and manipulated, free reactant concentrations are determined by a system of mass action-based total concentration constraint polynomials, expected measurements are determined by model predicted complex concentrations, and the number of models is large due to combinatorial complexity. This is a generic <italic>in vitro </italic>
synthetic biochemical system problem statement, so <monospace>ccems</monospace>
could have a broad impact.</p>
</sec>
<sec sec-type="methods"><title>Methods</title>
<p>Data were digitized by plotDigitizer [<xref ref-type="bibr" rid="B35">35</xref>
] and analyzed using <monospace>ccems</monospace>
[<xref ref-type="bibr" rid="B22">22</xref>
]. Hessians of SSEs obtained using <monospace>optim</monospace>
were divided by 2, inverted, multiplied by SSE/(<italic>N </italic>
- <italic>P</italic>
), and the square roots of the main diagonal were then multiplied by 1.96 to form 95% Wald CI. Parameters were estimated in exponentiated forms to constrain them to positive values.</p>
</sec>
<sec><title>Competing interests</title>
<p>The author declares that he has no competing interests.</p>
</sec>
<sec><title>Authors' contributions</title>
<p>TR performed the work and wrote the paper.</p>
</sec>
<sec><title>Reviewer's comments</title>
<sec><title>Reviewer's report 1</title>
<p>Ossama B. Kashlan, University of Pittsburgh (nominated by Philip Hahnfeldt, Tufts)</p>
<p>As you've shown, you don't need to invoke the <italic>h</italic>
-site to fit our figure <xref ref-type="fig" rid="F1">1</xref>
data. But we did need it for the other data in the paper, e.g. the global fit of DLS and activity data in figure <xref ref-type="fig" rid="F5">5</xref>
of our paper (with dTTP saturating the <italic>s</italic>
-site). Since we wanted to use a single model for all the data, we therefore used an <italic>h</italic>
-site to fit our figure <xref ref-type="fig" rid="F1">1</xref>
data. You should include a discussion of the potential of the <italic>h</italic>
-site-less models to account for our figure <xref ref-type="fig" rid="F5">5</xref>
DLS data (as you noted, modeling activity data greatly increases the model space).</p>
</sec>
<sec><title>Radivoyevitch's Response</title>
<p>Let us denote by R the dimer complex (R1)<sub>2</sub>
(GDP)<sub>2</sub>
(dTTP)<sub>2</sub>
. The full model is then<disp-formula id="bmcMR1"><label>(R1)</label>
<graphic xlink:href="1745-6150-4-50-i12.gif"></graphic>
</disp-formula>
</p>
<p>and in this case, the top 3 models are R<sup>2</sup>
X<sup>7 </sup>
+ R<sup>3</sup>
X<sup>12</sup>
, R<sup>2</sup>
X<sup>5 </sup>
+ R<sup>3</sup>
X<sup>9 </sup>
and R<sup>2</sup>
X<sup>6 </sup>
+ R<sup>3</sup>
X<sup>10 </sup>
+ R<sup>3</sup>
X<sup>12</sup>
, i.e. models with <italic>h</italic>
-site terms. Thus, you are indeed correct that there appears to be evidence for an <italic>h</italic>
-site in your figure <xref ref-type="fig" rid="F5">5</xref>
(of ref [<xref ref-type="bibr" rid="B20">20</xref>
]) data. Fits of these models are shown in the new Fig. <xref ref-type="fig" rid="F13">13</xref>
. If I remove the three datapoints with [ATP] = 5 mM, 7 mM and 10 mM, i.e. values higher than those in your figure <xref ref-type="fig" rid="F1">1</xref>
data, similar results are obtained. Thus, perhaps the <italic>h</italic>
-site is created by the presence of bound GDP substrate, or bound dTTP, or both, but yes, different inferences are drawn from your figure <xref ref-type="fig" rid="F1">1</xref>
and <xref ref-type="fig" rid="F5">5</xref>
datasets when they are analyzed individually. It should perhaps be noted that all of the tetramer terms above also have occupied <italic>h</italic>
-sites, and that the model that you used did not, i.e. the methods presented may have value in analyses of your figure <xref ref-type="fig" rid="F5">5</xref>
data as well. Further, since no 1-parameter models were among the top 100, they likely all yielded poor fits.</p>
<fig id="F13" position="float"><label>Figure 13</label>
<caption><p><bold>Fits of the top 3 models with [dTTP] and [GDP] at saturating levels</bold>
. In response to reviewer Kashlan's comment, the top 3 models of the data in figure 5 of his paper are shown. All models yield m = 180 at [ATP] = 0; m = 181 was measured.</p>
</caption>
<graphic xlink:href="1745-6150-4-50-13"></graphic>
</fig>
</sec>
<sec><title>Kashlan's Response</title>
<p>As you concluded from the DLS data with saturating dTTP, an <italic>h</italic>
-site is necessary to fit this data. However, I find your interpretation that the presence of GDP and/or dTTP possibly creates the <italic>h</italic>
-site to be unlikely. More likely is that ligands bound to one site have heterotropic binding effects on the other sites. Under this logic, the underlying binding K's for ATP in the absence of GDP and dTTP are such that the DLS data under these conditions lack the complexity to require the invocation of an <italic>h</italic>
-site. However, in the presence of GDP and dTTP, the underlying K's are such that this complexity is unmasked. I also have a few additional comments:</p>
<p>1. Regarding the assumption that <italic>s</italic>
- fills before <italic>a</italic>
- fills before <italic>h</italic>
-site, and that oligomers always have certain sites filled. This assumption ignores a few important observations and possibilities. First is the ability of R1 to dimerize in the absence of a filled <italic>s</italic>
-site, e.g. we observed that CDP reduction occurs (at a low rate) in the absence of an <italic>s</italic>
-site ligand. Second, as you pointed out, is the ability to switch <italic>s</italic>
-site ligands while the <italic>a</italic>
- (and/or <italic>h</italic>
-) site(s) are occupied. This is important, because you base your conclusion of a lack of <italic>h</italic>
-site evidence on the fact that not enough ATP are bound to fill the <italic>a</italic>
-sites. But unless <italic>h</italic>
-site binding is dramatically weaker than <italic>a</italic>
-site binding, binding the first few <italic>h</italic>
-sites may be more favorable than filling the last few <italic>a</italic>
-sites.</p>
<p>2. The minimal models presented as 'best', and from which physical conclusions are drawn, should be able to account for both (all) datasets. Having a different model framework for each given set of ligands adds a new level of complexity. Can you combine the AIC scores for the fits to both sets of data to find the best model(s), based at least on these data?</p>
<p>3. The conclusion on p. 3 should be edited to reflect the above comments and your response to my previous comment.</p>
<p>4. The background on p. 3 should be edited to read, "ATP binds to both of these sites and there is some evidence, based on RNR *mass and* activity versus [ATP] data,."</p>
</sec>
<sec><title>Radivoyevitch's Response</title>
<p>Two comments: 1) if we let site creation include increasing the affinity so that an infinite K is now finite, we may be saying the same things; and 2) you have 9 data points in your figure <xref ref-type="fig" rid="F5">5</xref>
and 16 in your figure <xref ref-type="fig" rid="F1">1</xref>
, and it would have been better to have these sample numbers reversed if indeed the conditions of your figure <xref ref-type="fig" rid="F5">5</xref>
yield more complex data. My inclination is to trust a 1-parameter model fitted to 15 points more than a 2-parameter model fitted to 9. Regarding your other two remarks:</p>
<p>1. By stating model assumptions I did not mean to say that I thought they must be true (no model is ever correct). What I meant to say is that my inferences are all conditional on their truth. The situation is such that unless such an assumption is made about the binding order, the meaning of a polynomial term that exists in a model is ambiguous. This is a big weakness, but I do not see any way around it. I now state this weakness at the end of the Limitations Section.</p>
<p>2. Note that we really have two different hub proteins for these two datasets, one with R as R1 monomer and the other with R as dTTP and GDP saturated R1 dimer. All complete K parameters downstream of these two hubs would thus be independent and there would thus be no way to pool the parameter estimates (beyond the error estimate). The joint model would thus be the sum of the terms.</p>
<p>3. The conclusion in the abstract has been softened to reflect conditionality on the experimental conditions of your figure <xref ref-type="fig" rid="F1">1</xref>
dataset.</p>
<p>4. Agreed, but this sentence no longer exists.</p>
</sec>
<sec><title>Reviewer's report 2</title>
<p>Bin Hu, Los Alamos National Laboratory (nominated by William Hlavacek, LANL)</p>
<p>In this work, the author did a theoretical study on the equilibrium of ATP-induced hexamerization of the R1 subunit of ribonucleotide reductase (RNR), where combinatorial complexity is observed. Statistical hypotheses with assumptions were used to generate an array of models of the hexamerization process. Statistical comparison of the model structures suggests that a-sites may not be occupied by ATP in R1 hexamers at physiological ATP concentrations and the <italic>h</italic>
-site may not exist. Results from this work suggest that the work in reference [<xref ref-type="bibr" rid="B20">20</xref>
] did not consider the possibility of a-sites not being occupied, which allows the authors in [<xref ref-type="bibr" rid="B20">20</xref>
] to suggest that the <italic>h</italic>
-site exists. A final judgment about whether the <italic>h</italic>
-site exists; however, can only be reached through additional experimental work. Although the statistical analysis approach used in this work is interesting, the author may wish to add a discussion comparing results from this work to those of Ref [<xref ref-type="bibr" rid="B20">20</xref>
]. Minor comments are as follows:</p>
<p>1. An introduction to RNR, its function and regulation is needed. Currently there is only one sentence in the beginning of the Introduction about RNR.</p>
<p>2. It is not mentioned in the paper how the R1 subunit forms RNR with R2 units and what kind of multimers the R1 subunit can form and their biological importance. References, especially for the crystallography data, are needed.</p>
<p>3. In the Results section, the meaning of "complete dissociation constants" is not clear.</p>
<p>4. The R package ccems was first introduced without reference.</p>
<p>5. The author submitted two papers to this journal simultaneously. Instead of using "the accompanying paper," the author or editorial office may want to change it to some other description that may help readers to find out which paper the author is referencing.</p>
<p>6. I cannot tell whether the assumption "s-sites are always filled in oligomers and never filled in monomers" is acceptable in this study.</p>
<p>7. It would be interesting if the author compared the parameters used this work with those used in [<xref ref-type="bibr" rid="B20">20</xref>
].</p>
</sec>
<sec><title>Radivoyevitch's Responses</title>
<p>The model of ATP induced R1 hexamerization previously proposed by Kashlan <italic>et al</italic>
. [<xref ref-type="bibr" rid="B20">20</xref>
] assumed: a) that the binary dissociation constants <italic>K </italic>
of the ATP binding sites <italic>s</italic>
, <italic>a </italic>
and <italic>h </italic>
are the same within oligomers (within site types); b) that these <italic>K</italic>
<sub>s</sub>
, <italic>K</italic>
<sub>a </sub>
and <italic>K</italic>
<sub>h </sub>
are infinite in structures smaller than dimers, tetramers and hexamers, respectively (note that all three are thus infinite in monomers); c) that finite <italic>K </italic>
are equal wherever it is plausible that they might be, i.e. that <italic>K</italic>
<sub>a </sub>
in tetramers equals <italic>K</italic>
<sub>a </sub>
in hexamers and that <italic>K</italic>
<sub>s </sub>
is the same across dimers, tetramers and hexamers; d) that the dissociation constants for R1 binding to itself (<italic>K</italic>
<sub>R_R</sub>
), R1 dimers binding to themselves (<italic>K</italic>
<sub>R2_R2</sub>
), and R1 tetramers binding to R1 dimers (<italic>K</italic>
<sub>R4_R2</sub>
), are independently adjustable; and e) that R1 tetramers can isomerize with an isomerization constant <italic>K</italic>
<sub>is</sub>
. Assumptions a) to c) constrain the model and d) to e) broaden it. When Kashlan <italic>et al</italic>
. fitted their model to their DLS data, <italic>K</italic>
<sub>R_R</sub>
, <italic>K</italic>
<sub>R2_R2</sub>
, <italic>K</italic>
<sub>R4_R2 </sub>
and <italic>K</italic>
<sub>is </sub>
were treated as being independent of R1 ligands, and consistent with these assumptions, the data in their figures <xref ref-type="fig" rid="F1">1</xref>
and 3-9 were fitted to single values of these constants such that the fits in these figures did not appear too poor). With respect to their figure <xref ref-type="fig" rid="F1">1</xref>
data, however, the first five residuals of their fit were negative and thus correlated, and although the residuals were small and thus subtle, the fit was thus poor. This paper focuses on their figure <xref ref-type="fig" rid="F1">1</xref>
data alone.</p>
<p>Regarding proof that an h-site does or does not exist, I agree that binding studies must be performed to see how many ATPs actually bind to R1, but such studies may be difficult, as evidenced by the fact that they have not yet been performed. My hunch is that experimental challenges are associated with weak ATP binding and thus the high [ATP] needed to achieve binding, which may make changes in free [ATP] due to free ATP losses to ATP bound to R1 difficult to detect. Regarding comparisons of results, since the single model that they fitted to their DLS data was also based on their RNR activity data, and since activity data is much more complicated to analyze than mass data because conjectures about different activity parameters being zero or equal to each other greatly expands the model space further, and since my software is not yet ready for more than one type of oligomer in activity data analyses (in the accompanying paper the enzyme TK1 was strictly tetrameric), our results cannot be compared. Regarding your points:</p>
<p>1. Background material regarding dNTP supply and RNR have now been added to the Introduction. However, since I do not model RNR activity data, this work has limited relevance to dNTP supply metabolism. My focus is thus on the methods developed. Indeed, it may be best to view R1 as merely some protein that has either 2 or 3 binding sites on it for a ligand that induces its hexamerization.</p>
<p>2. A) R2 is irrelevant here since this paper does not delve into activity data and since R2 was not present in the experiment that yielded the DLS average mass data analyzed. B) The first eukaryotic (yeast) R1 structure showed a dimer and this was referenced [<xref ref-type="bibr" rid="B26">26</xref>
]. Though a dATP induced R1 tetramer was observed in Ref [<xref ref-type="bibr" rid="B20">20</xref>
], it was not observed in [<xref ref-type="bibr" rid="B27">27</xref>
], and no lab has observed it directly using the more relevant ligand (for this paper) ATP. Thus, tetramers could perhaps have been left out of the model space, but there is strong support for R1 monomers (e.g. the low [ATP] DLS data in Fig. <xref ref-type="fig" rid="F7">7</xref>
), dimers (the structure in [<xref ref-type="bibr" rid="B26">26</xref>
]) and hexamers (e.g. the high [ATP] DLS data in Fig. <xref ref-type="fig" rid="F7">7</xref>
).</p>
<p>3. By complete dissociation constants I mean those where the Gibbs Free energies are with respect to all reactants being completely separated from one another by infinite distances. In contrast, by binary dissociation constants I mean situations where only one reactant (or perhaps a subcomplex) is separated out at infinity while all of the other reactants remain bound together.</p>
<p>4. The link to my ccems page is now referenced earlier.</p>
<p>5. They should end up back-to-back and I hope readers will read, and know of, both.</p>
<p>6. The conclusions remain the same if I drop <italic>h</italic>
-site terms and blow up the model space by introducing s-site terms, i.e. there is some support for the assumption besides dTTP induced R1dimerization results in [<xref ref-type="bibr" rid="B25">25</xref>
] and the structure in [<xref ref-type="bibr" rid="B26">26</xref>
].</p>
<p>7. The model used in [<xref ref-type="bibr" rid="B20">20</xref>
] has 7 parameters. Their binary K values for ATP binding were 100 μM for dimers and tetramers and 1.1 mM for hexamers. For R1 oligomerization and isomerization their values were <italic>K</italic>
<sub>R_R </sub>
= 170 μM, <italic>K</italic>
<sub>R2_R2 </sub>
= 2-5 mM, <italic>K</italic>
<sub>R4_R2 </sub>
= 2-6 mM and <italic>K</italic>
<sub>is </sub>
= 10-40 where ranges depend on different tetrameric activity assumptions. Without the same parameters in my models, comparisons are difficult. The best model in Table <xref ref-type="table" rid="T1">1</xref>
is R6X8 and it has a geometric mean binary binding constant of 63 μM. Since all of the binary K values of [<xref ref-type="bibr" rid="B20">20</xref>
] are ≥ 100 μM, their geometric mean must also be ≥ 100 μM. Indeed, aassuming ATP fills 6 s-sites and 2 a-sites in R6X8, one obtains a geometric mean binary K of 190 μM = ((100)<sup>8</sup>
(170)<sup>3</sup>
(3000)<sup>1</sup>
(3000)<sup>1</sup>
)<sup>1/13</sup>
, i.e. there is a difference of a factor of 3 in parameter estimates between analyses.</p>
</sec>
<sec><title>Reviewer's report 3</title>
<sec><title>Rainer K. Sachs, UC Berkeley</title>
<p>In general, is there some systematic rationale and/or underlying reasoning on what criteria to use to distinguish hypotheses? Do you estimate that almost any criterion would give the same final answers qualitatively? How did you decide on 30 models in your comparison of model proportions? What is your main motivation to study RNR? Can your software handle more than two reactants? Can you provide a code use example?</p>
</sec>
</sec>
<sec><title>Radivoyevitch's Responses</title>
<p>The AIC was used (without consideration of alternatives such as the Bayesian Information Criterion, BIC) only because it is the most popular. The idea was to pick a criterion to present my main contribution, which is in model space generation rather than model selection. Though there may be reasons to switch to a different criterion that I have yet to learn of, in the interim, the AIC is my default. I suspect that the conclusions made here would be robust to such changes.</p>
<p>The choice of 30 models involved data snooping as you may have guessed, i.e. 30 looked like a good breaking point for a claim that the data does not support an <italic>h</italic>
-site. Thus, the difference in proportions p-value that I reported may be overly significant. Nevertheless, the Kolmogorov-Smironov Test is with respect to entire distributions in Fig. <xref ref-type="fig" rid="F6">6</xref>
, so the conclusion that this particular dataset does not demand the existence of an <italic>h</italic>
-site is robust. There is no indication that this conclusion will hold under different experimental conditions, however, e.g. see Kashlan's review above.</p>
<p>Three paragraphs were added to the Introduction to motivate RNR research. The automated model space generation capabilities of <monospace>ccems</monospace>
are indeed limited to two reactants and this is now stated in the title and elsewhere. If R and <monospace>ccems</monospace>
are installed, the R command <monospace>load(ccems)</monospace>
followed by <monospace>?ccems</monospace>
yields help which includes the code example in Fig. <xref ref-type="fig" rid="F14">14</xref>
which can be pasted into the R command line to create the model space used here.</p>
<fig id="F14" position="float"><label>Figure 14</label>
<caption><p><bold>ccems code example</bold>
. These codes generate Table 1 and the RX model space used in this paper.</p>
</caption>
<graphic xlink:href="1745-6150-4-50-14"></graphic>
</fig>
</sec>
</sec>
</body>
<back><sec><title>Acknowledgements</title>
<p>The project described was supported by Award Number K25CA104791 from the National Cancer Institute. The content is solely the responsibility of the author and does not necessarily represent the official views of the National Cancer Institute or the National Institutes of Health.</p>
</sec>
<ref-list><ref id="B1"><mixed-citation publication-type="journal"><name><surname>Radivoyevitch</surname>
<given-names>T</given-names>
</name>
<name><surname>Loparo</surname>
<given-names>KA</given-names>
</name>
<name><surname>Jackson</surname>
<given-names>RC</given-names>
</name>
<name><surname>Sedwick</surname>
<given-names>WD</given-names>
</name>
<article-title>On systems and control approaches to therapeutic gain</article-title>
<source>BMC Cancer</source>
<year>2006</year>
<volume>6</volume>
<issue>1</issue>
<fpage>104</fpage>
<pub-id pub-id-type="doi">10.1186/1471-2407-6-104</pub-id>
<pub-id pub-id-type="pmid">16638124</pub-id>
</mixed-citation>
</ref>
<ref id="B2"><mixed-citation publication-type="journal"><name><surname>Manegold</surname>
<given-names>C</given-names>
</name>
<article-title>Gemcitabine (Gemzar) in non-small cell lung cancer</article-title>
<source>Expert Rev Anticancer Ther</source>
<year>2004</year>
<volume>4</volume>
<issue>3</issue>
<fpage>345</fpage>
<lpage>360</lpage>
<pub-id pub-id-type="doi">10.1586/14737140.4.3.345</pub-id>
<pub-id pub-id-type="pmid">15161434</pub-id>
</mixed-citation>
</ref>
<ref id="B3"><mixed-citation publication-type="journal"><name><surname>Kroep</surname>
<given-names>JR</given-names>
</name>
<name><surname>Van Moorsel</surname>
<given-names>CJ</given-names>
</name>
<name><surname>Veerman</surname>
<given-names>G</given-names>
</name>
<name><surname>Voorn</surname>
<given-names>DA</given-names>
</name>
<name><surname>Schultz</surname>
<given-names>RM</given-names>
</name>
<name><surname>Worzalla</surname>
<given-names>JF</given-names>
</name>
<name><surname>Tanzer</surname>
<given-names>LR</given-names>
</name>
<name><surname>Merriman</surname>
<given-names>RL</given-names>
</name>
<name><surname>Pinedo</surname>
<given-names>HM</given-names>
</name>
<name><surname>Peters</surname>
<given-names>GJ</given-names>
</name>
<article-title>Role of deoxycytidine kinase (dCK), thymidine kinase 2 (TK2), and deoxycytidine deaminase (dCDA) in the antitumor activity of gemcitabine (dFdC)</article-title>
<source>Adv Exp Med Biol</source>
<year>1998</year>
<volume>431</volume>
<fpage>657</fpage>
<lpage>660</lpage>
<pub-id pub-id-type="pmid">9598147</pub-id>
</mixed-citation>
</ref>
<ref id="B4"><mixed-citation publication-type="journal"><name><surname>Shewach</surname>
<given-names>DS</given-names>
</name>
<name><surname>Reynolds</surname>
<given-names>KK</given-names>
</name>
<name><surname>Hertel</surname>
<given-names>L</given-names>
</name>
<article-title>Nucleotide specificity of human deoxycytidine kinase</article-title>
<source>Mol Pharmacol</source>
<year>1992</year>
<volume>42</volume>
<issue>3</issue>
<fpage>518</fpage>
<lpage>524</lpage>
<pub-id pub-id-type="pmid">1406603</pub-id>
</mixed-citation>
</ref>
<ref id="B5"><mixed-citation publication-type="journal"><name><surname>Wang</surname>
<given-names>J</given-names>
</name>
<name><surname>Lohman</surname>
<given-names>GJ</given-names>
</name>
<name><surname>Stubbe</surname>
<given-names>J</given-names>
</name>
<article-title>Enhanced subunit interactions with gemcitabine-5'-diphosphate inhibit ribonucleotide reductases</article-title>
<source>Proc Natl Acad Sci USA</source>
<year>2007</year>
<volume>104</volume>
<issue>36</issue>
<fpage>14324</fpage>
<lpage>14329</lpage>
<pub-id pub-id-type="doi">10.1073/pnas.0706803104</pub-id>
<pub-id pub-id-type="pmid">17726094</pub-id>
</mixed-citation>
</ref>
<ref id="B6"><mixed-citation publication-type="other"><name><surname>Ferrandina</surname>
<given-names>G</given-names>
</name>
<name><surname>Mey</surname>
<given-names>V</given-names>
</name>
<name><surname>Nannizzi</surname>
<given-names>S</given-names>
</name>
<name><surname>Ricciardi</surname>
<given-names>S</given-names>
</name>
<name><surname>Petrillo</surname>
<given-names>M</given-names>
</name>
<name><surname>Ferlini</surname>
<given-names>C</given-names>
</name>
<name><surname>Danesi</surname>
<given-names>R</given-names>
</name>
<name><surname>Scambia</surname>
<given-names>G</given-names>
</name>
<name><surname>Del Tacca</surname>
<given-names>M</given-names>
</name>
<article-title>Expression of nucleoside transporters, deoxycitidine kinase, ribonucleotide reductase regulatory subunits, and gemcitabine catabolic enzymes in primary ovarian cancer</article-title>
<source>Cancer Chemother Pharmacol</source>
<year>2009</year>
<pub-id pub-id-type="pmid">19639316</pub-id>
</mixed-citation>
</ref>
<ref id="B7"><mixed-citation publication-type="journal"><name><surname>Ohtaka</surname>
<given-names>K</given-names>
</name>
<name><surname>Kohya</surname>
<given-names>N</given-names>
</name>
<name><surname>Sato</surname>
<given-names>K</given-names>
</name>
<name><surname>Kitajima</surname>
<given-names>Y</given-names>
</name>
<name><surname>Ide</surname>
<given-names>T</given-names>
</name>
<name><surname>Mitsuno</surname>
<given-names>M</given-names>
</name>
<name><surname>Miyazaki</surname>
<given-names>K</given-names>
</name>
<article-title>Ribonucleotide reductase subunit M1 is a possible chemoresistance marker to gemcitabine in biliary tract carcinoma</article-title>
<source>Oncol Rep</source>
<year>2008</year>
<volume>20</volume>
<issue>2</issue>
<fpage>279</fpage>
<lpage>286</lpage>
<pub-id pub-id-type="pmid">18636187</pub-id>
</mixed-citation>
</ref>
<ref id="B8"><mixed-citation publication-type="journal"><name><surname>Flanagan</surname>
<given-names>SA</given-names>
</name>
<name><surname>Robinson</surname>
<given-names>BW</given-names>
</name>
<name><surname>Krokosky</surname>
<given-names>CM</given-names>
</name>
<name><surname>Shewach</surname>
<given-names>DS</given-names>
</name>
<article-title>Mismatched nucleotides as the lesions responsible for radiosensitization with gemcitabine: a new paradigm for antimetabolite radiosensitizers</article-title>
<source>Mol Cancer Ther</source>
<year>2007</year>
<volume>6</volume>
<issue>6</issue>
<fpage>1858</fpage>
<lpage>1868</lpage>
<pub-id pub-id-type="doi">10.1158/1535-7163.MCT-07-0068</pub-id>
<pub-id pub-id-type="pmid">17575114</pub-id>
</mixed-citation>
</ref>
<ref id="B9"><mixed-citation publication-type="journal"><name><surname>Muller</surname>
<given-names>A</given-names>
</name>
<name><surname>Fishel</surname>
<given-names>R</given-names>
</name>
<article-title>Mismatch repair and the hereditary non-polyposis colorectal cancer syndrome (HNPCC)</article-title>
<source>Cancer Invest</source>
<year>2002</year>
<volume>20</volume>
<issue>1</issue>
<fpage>102</fpage>
<lpage>109</lpage>
<pub-id pub-id-type="doi">10.1081/CNV-120000371</pub-id>
<pub-id pub-id-type="pmid">11852992</pub-id>
</mixed-citation>
</ref>
<ref id="B10"><mixed-citation publication-type="journal"><name><surname>Leung</surname>
<given-names>WK</given-names>
</name>
<name><surname>Kim</surname>
<given-names>JJ</given-names>
</name>
<name><surname>Kim</surname>
<given-names>JG</given-names>
</name>
<name><surname>Graham</surname>
<given-names>DY</given-names>
</name>
<name><surname>Sepulveda</surname>
<given-names>AR</given-names>
</name>
<article-title>Microsatellite instability in gastric intestinal metaplasia in patients with and without gastric cancer</article-title>
<source>Am J Pathol</source>
<year>2000</year>
<volume>156</volume>
<issue>2</issue>
<fpage>537</fpage>
<lpage>543</lpage>
<pub-id pub-id-type="pmid">10666383</pub-id>
</mixed-citation>
</ref>
<ref id="B11"><mixed-citation publication-type="journal"><name><surname>Nakata</surname>
<given-names>B</given-names>
</name>
<name><surname>Wang</surname>
<given-names>YQ</given-names>
</name>
<name><surname>Yashiro</surname>
<given-names>M</given-names>
</name>
<name><surname>Nishioka</surname>
<given-names>N</given-names>
</name>
<name><surname>Tanaka</surname>
<given-names>H</given-names>
</name>
<name><surname>Ohira</surname>
<given-names>M</given-names>
</name>
<name><surname>Ishikawa</surname>
<given-names>T</given-names>
</name>
<name><surname>Nishino</surname>
<given-names>H</given-names>
</name>
<name><surname>Hirakawa</surname>
<given-names>K</given-names>
</name>
<article-title>Prognostic value of microsatellite instability in resectable pancreatic cancer</article-title>
<source>Clin Cancer Res</source>
<year>2002</year>
<volume>8</volume>
<issue>8</issue>
<fpage>2536</fpage>
<lpage>2540</lpage>
<pub-id pub-id-type="pmid">12171881</pub-id>
</mixed-citation>
</ref>
<ref id="B12"><mixed-citation publication-type="journal"><name><surname>Hartmann</surname>
<given-names>A</given-names>
</name>
<name><surname>Zanardo</surname>
<given-names>L</given-names>
</name>
<name><surname>Bocker-Edmonston</surname>
<given-names>T</given-names>
</name>
<name><surname>Blaszyk</surname>
<given-names>H</given-names>
</name>
<name><surname>Dietmaier</surname>
<given-names>W</given-names>
</name>
<name><surname>Stoehr</surname>
<given-names>R</given-names>
</name>
<name><surname>Cheville</surname>
<given-names>JC</given-names>
</name>
<name><surname>Junker</surname>
<given-names>K</given-names>
</name>
<name><surname>Wieland</surname>
<given-names>W</given-names>
</name>
<name><surname>Knuechel R Rueschoff</surname>
<given-names>J</given-names>
</name>
<name><surname>Hofstaedter</surname>
<given-names>F</given-names>
</name>
<name><surname>Fishel</surname>
<given-names>R</given-names>
</name>
<article-title>Frequent microsatellite instability in sporadic tumors of the upper urinary tract</article-title>
<source>Cancer Res</source>
<year>2002</year>
<volume>62</volume>
<issue>23</issue>
<fpage>6796</fpage>
<lpage>6802</lpage>
<pub-id pub-id-type="pmid">12460887</pub-id>
</mixed-citation>
</ref>
<ref id="B13"><mixed-citation publication-type="journal"><name><surname>Ichikawa</surname>
<given-names>Y</given-names>
</name>
<name><surname>Lemon</surname>
<given-names>SJ</given-names>
</name>
<name><surname>Wang</surname>
<given-names>S</given-names>
</name>
<name><surname>Franklin</surname>
<given-names>B</given-names>
</name>
<name><surname>Watson</surname>
<given-names>P</given-names>
</name>
<name><surname>Knezetic</surname>
<given-names>JA</given-names>
</name>
<name><surname>Bewtra</surname>
<given-names>C</given-names>
</name>
<name><surname>Lynch</surname>
<given-names>HT</given-names>
</name>
<article-title>Microsatellite instability and expression of MLH1 and MSH2 in normal and malignant endometrial and ovarian epithelium in hereditary nonpolyposis colorectal cancer family members</article-title>
<source>Cancer Genet Cytogenet</source>
<year>1999</year>
<volume>112</volume>
<issue>1</issue>
<fpage>2</fpage>
<lpage>8</lpage>
<pub-id pub-id-type="doi">10.1016/S0165-4608(98)00252-0</pub-id>
<pub-id pub-id-type="pmid">10432927</pub-id>
</mixed-citation>
</ref>
<ref id="B14"><mixed-citation publication-type="journal"><name><surname>Simpkins</surname>
<given-names>SB</given-names>
</name>
<name><surname>Bocker</surname>
<given-names>T</given-names>
</name>
<name><surname>Swisher</surname>
<given-names>EM</given-names>
</name>
<name><surname>Mutch</surname>
<given-names>DG</given-names>
</name>
<name><surname>Gersell</surname>
<given-names>DJ</given-names>
</name>
<name><surname>Kovatich</surname>
<given-names>AJ</given-names>
</name>
<name><surname>Palazzo</surname>
<given-names>JP</given-names>
</name>
<name><surname>Fishel</surname>
<given-names>R</given-names>
</name>
<name><surname>Goodfellow</surname>
<given-names>PJ</given-names>
</name>
<article-title>MLH1 promoter methylation and gene silencing is the primary cause of microsatellite instability in sporadic endometrial cancers</article-title>
<source>Hum Mol Genet</source>
<year>1999</year>
<volume>8</volume>
<issue>4</issue>
<fpage>661</fpage>
<lpage>666</lpage>
<pub-id pub-id-type="doi">10.1093/hmg/8.4.661</pub-id>
<pub-id pub-id-type="pmid">10072435</pub-id>
</mixed-citation>
</ref>
<ref id="B15"><mixed-citation publication-type="journal"><name><surname>Leung</surname>
<given-names>SY</given-names>
</name>
<name><surname>Chan</surname>
<given-names>TL</given-names>
</name>
<name><surname>Chung</surname>
<given-names>LP</given-names>
</name>
<name><surname>Chan</surname>
<given-names>AS</given-names>
</name>
<name><surname>Fan</surname>
<given-names>YW</given-names>
</name>
<name><surname>Hung</surname>
<given-names>KN</given-names>
</name>
<name><surname>Kwong</surname>
<given-names>WK</given-names>
</name>
<name><surname>Ho</surname>
<given-names>JW</given-names>
</name>
<name><surname>Yuen</surname>
<given-names>ST</given-names>
</name>
<article-title>Microsatellite instability and mutation of DNA mismatch repair genes in gliomas</article-title>
<source>Am J Pathol</source>
<year>1998</year>
<volume>153</volume>
<issue>4</issue>
<fpage>1181</fpage>
<lpage>1188</lpage>
<pub-id pub-id-type="pmid">9777949</pub-id>
</mixed-citation>
</ref>
<ref id="B16"><mixed-citation publication-type="journal"><name><surname>Thelander</surname>
<given-names>L</given-names>
</name>
<name><surname>Reichard</surname>
<given-names>P</given-names>
</name>
<article-title>Reduction of ribonucleotides</article-title>
<source>Annu Rev Biochem</source>
<year>1979</year>
<volume>48</volume>
<fpage>133</fpage>
<lpage>158</lpage>
<pub-id pub-id-type="doi">10.1146/annurev.bi.48.070179.001025</pub-id>
<pub-id pub-id-type="pmid">382982</pub-id>
</mixed-citation>
</ref>
<ref id="B17"><mixed-citation publication-type="journal"><name><surname>Tanaka</surname>
<given-names>H</given-names>
</name>
<name><surname>Arakawa</surname>
<given-names>H</given-names>
</name>
<name><surname>Yamaguchi</surname>
<given-names>T</given-names>
</name>
<name><surname>Shiraishi</surname>
<given-names>K</given-names>
</name>
<name><surname>Fukuda</surname>
<given-names>S</given-names>
</name>
<name><surname>Matsui</surname>
<given-names>K</given-names>
</name>
<name><surname>Takei</surname>
<given-names>Y</given-names>
</name>
<name><surname>Nakamura</surname>
<given-names>Y</given-names>
</name>
<article-title>A ribonucleotide reductase gene involved in a p53-dependent cell-cycle checkpoint for DNA damage</article-title>
<source>Nature</source>
<year>2000</year>
<volume>404</volume>
<issue>6773</issue>
<fpage>42</fpage>
<lpage>49</lpage>
<pub-id pub-id-type="doi">10.1038/35003506</pub-id>
<pub-id pub-id-type="pmid">10716435</pub-id>
</mixed-citation>
</ref>
<ref id="B18"><mixed-citation publication-type="journal"><name><surname>Nakano</surname>
<given-names>K</given-names>
</name>
<name><surname>Balint</surname>
<given-names>E</given-names>
</name>
<name><surname>Ashcroft</surname>
<given-names>M</given-names>
</name>
<name><surname>Vousden</surname>
<given-names>KH</given-names>
</name>
<article-title>A ribonucleotide reductase gene is a transcriptional target of p53 and p73</article-title>
<source>Oncogene</source>
<year>2000</year>
<volume>19</volume>
<issue>37</issue>
<fpage>4283</fpage>
<lpage>4289</lpage>
<pub-id pub-id-type="doi">10.1038/sj.onc.1203774</pub-id>
<pub-id pub-id-type="pmid">10980602</pub-id>
</mixed-citation>
</ref>
<ref id="B19"><mixed-citation publication-type="other"><name><surname>Feng</surname>
<given-names>IJ</given-names>
</name>
<name><surname>Radivoyevitch</surname>
<given-names>T</given-names>
</name>
<article-title>SNP-SNP interactions between the dNTP supply and mismatch DNA repair systems in breast cancer</article-title>
<source>Ohio Collaborative Conference on Bioinformatics</source>
<year>2009</year>
</mixed-citation>
</ref>
<ref id="B20"><mixed-citation publication-type="journal"><name><surname>Kashlan</surname>
<given-names>OB</given-names>
</name>
<name><surname>Scott</surname>
<given-names>CP</given-names>
</name>
<name><surname>Lear</surname>
<given-names>JD</given-names>
</name>
<name><surname>Cooperman</surname>
<given-names>BS</given-names>
</name>
<article-title>A comprehensive model for the allosteric regulation of mammalian ribonucleotide reductase. Functional consequences of ATP- and dATP-induced oligomerization of the large subunit</article-title>
<source>Biochemistry</source>
<year>2002</year>
<volume>41</volume>
<issue>2</issue>
<fpage>462</fpage>
<lpage>474</lpage>
<pub-id pub-id-type="doi">10.1021/bi011653a</pub-id>
<pub-id pub-id-type="pmid">11781084</pub-id>
</mixed-citation>
</ref>
<ref id="B21"><mixed-citation publication-type="journal"><name><surname>Ihaka</surname>
<given-names>R</given-names>
</name>
<name><surname>Gentleman</surname>
<given-names>R</given-names>
</name>
<article-title>R:a language for data analysis and graphics</article-title>
<source>Journal of Computational and graphical statistics</source>
<year>1996</year>
<volume>5</volume>
<fpage>299</fpage>
<lpage>314</lpage>
<pub-id pub-id-type="doi">10.2307/1390807</pub-id>
</mixed-citation>
</ref>
<ref id="B22"><mixed-citation publication-type="other"><article-title>Combinatorially Complex Equilibrium Model Selection</article-title>
<ext-link ext-link-type="uri" xlink:href="http://epbi-radivot.cwru.edu/ccems/overview.html">http://epbi-radivot.cwru.edu/ccems/overview.html</ext-link>
</mixed-citation>
</ref>
<ref id="B23"><mixed-citation publication-type="journal"><name><surname>Akaike</surname>
<given-names>H</given-names>
</name>
<article-title>A new look at the statistical model identification</article-title>
<source>IEEE Transactions on Automatic Control</source>
<year>1974</year>
<volume>19</volume>
<fpage>716</fpage>
<lpage>723</lpage>
<pub-id pub-id-type="doi">10.1109/TAC.1974.1100705</pub-id>
</mixed-citation>
</ref>
<ref id="B24"><mixed-citation publication-type="book"><name><surname>Burnham</surname>
<given-names>KP</given-names>
</name>
<name><surname>Anderson</surname>
<given-names>DR</given-names>
</name>
<source>Model Selection and Multimodel Inference: A Practical-Theoretic Approach</source>
<year>2002</year>
<publisher-name>Springer-Verlag</publisher-name>
</mixed-citation>
</ref>
<ref id="B25"><mixed-citation publication-type="journal"><name><surname>Radivoyevitch</surname>
<given-names>T</given-names>
</name>
<article-title>Equilibrium model selection: dTTP induced R1 dimerization</article-title>
<source>BMC Syst Biol</source>
<year>2008</year>
<volume>2</volume>
<issue>1</issue>
<fpage>15</fpage>
<pub-id pub-id-type="doi">10.1186/1752-0509-2-15</pub-id>
<pub-id pub-id-type="pmid">18248678</pub-id>
</mixed-citation>
</ref>
<ref id="B26"><mixed-citation publication-type="journal"><name><surname>Xu</surname>
<given-names>H</given-names>
</name>
<name><surname>Faber</surname>
<given-names>C</given-names>
</name>
<name><surname>Uchiki</surname>
<given-names>T</given-names>
</name>
<name><surname>Fairman</surname>
<given-names>JW</given-names>
</name>
<name><surname>Racca</surname>
<given-names>J</given-names>
</name>
<name><surname>Dealwis</surname>
<given-names>C</given-names>
</name>
<article-title>Structures of eukaryotic ribonucleotide reductase I provide insights into dNTP regulation</article-title>
<source>Proc Natl Acad Sci USA</source>
<year>2006</year>
<volume>103</volume>
<issue>11</issue>
<fpage>4022</fpage>
<lpage>4027</lpage>
<pub-id pub-id-type="doi">10.1073/pnas.0600443103</pub-id>
<pub-id pub-id-type="pmid">16537479</pub-id>
</mixed-citation>
</ref>
<ref id="B27"><mixed-citation publication-type="journal"><name><surname>Rofougaran</surname>
<given-names>R</given-names>
</name>
<name><surname>Vodnala</surname>
<given-names>M</given-names>
</name>
<name><surname>Hofer</surname>
<given-names>A</given-names>
</name>
<article-title>Enzymatically active mammalian ribonucleotide reductase exists primarily as an alpha6beta2 octamer</article-title>
<source>J Biol Chem</source>
<year>2006</year>
<volume>281</volume>
<issue>38</issue>
<fpage>27705</fpage>
<lpage>27711</lpage>
<pub-id pub-id-type="doi">10.1074/jbc.M605573200</pub-id>
<pub-id pub-id-type="pmid">16861739</pub-id>
</mixed-citation>
</ref>
<ref id="B28"><mixed-citation publication-type="journal"><name><surname>Roach</surname>
<given-names>LS</given-names>
</name>
<name><surname>Song</surname>
<given-names>H</given-names>
</name>
<name><surname>Ismagilov</surname>
<given-names>RF</given-names>
</name>
<article-title>Controlling nonspecific protein adsorption in a plug-based microfluidic system by controlling interfacial chemistry using fluorous-phase surfactants</article-title>
<source>Anal Chem</source>
<year>2005</year>
<volume>77</volume>
<issue>3</issue>
<fpage>785</fpage>
<lpage>796</lpage>
<pub-id pub-id-type="doi">10.1021/ac049061w</pub-id>
<pub-id pub-id-type="pmid">15679345</pub-id>
</mixed-citation>
</ref>
<ref id="B29"><mixed-citation publication-type="journal"><name><surname>Song</surname>
<given-names>H</given-names>
</name>
<name><surname>Ismagilov</surname>
<given-names>RF</given-names>
</name>
<article-title>Millisecond kinetics on a microfluidic chip using nanoliters of reagents</article-title>
<source>J Am Chem Soc</source>
<year>2003</year>
<volume>125</volume>
<issue>47</issue>
<fpage>14613</fpage>
<lpage>14619</lpage>
<pub-id pub-id-type="doi">10.1021/ja0354566</pub-id>
<pub-id pub-id-type="pmid">14624612</pub-id>
</mixed-citation>
</ref>
<ref id="B30"><mixed-citation publication-type="journal"><name><surname>Huebner</surname>
<given-names>A</given-names>
</name>
<name><surname>Olguin</surname>
<given-names>LF</given-names>
</name>
<name><surname>Bratton</surname>
<given-names>D</given-names>
</name>
<name><surname>Whyte</surname>
<given-names>G</given-names>
</name>
<name><surname>Huck</surname>
<given-names>WT</given-names>
</name>
<name><surname>de Mello</surname>
<given-names>AJ</given-names>
</name>
<name><surname>Edel</surname>
<given-names>JB</given-names>
</name>
<name><surname>Abell</surname>
<given-names>C</given-names>
</name>
<name><surname>Hollfelder</surname>
<given-names>F</given-names>
</name>
<article-title>Development of quantitative cell-based enzyme assays in microdroplets</article-title>
<source>Anal Chem</source>
<year>2008</year>
<volume>80</volume>
<issue>10</issue>
<fpage>3890</fpage>
<lpage>3896</lpage>
<pub-id pub-id-type="doi">10.1021/ac800338z</pub-id>
<pub-id pub-id-type="pmid">18399662</pub-id>
</mixed-citation>
</ref>
<ref id="B31"><mixed-citation publication-type="journal"><name><surname>Huebner</surname>
<given-names>A</given-names>
</name>
<name><surname>Bratton</surname>
<given-names>D</given-names>
</name>
<name><surname>Whyte</surname>
<given-names>G</given-names>
</name>
<name><surname>Yang</surname>
<given-names>M</given-names>
</name>
<name><surname>Demello</surname>
<given-names>AJ</given-names>
</name>
<name><surname>Abell</surname>
<given-names>C</given-names>
</name>
<name><surname>Hollfelder</surname>
<given-names>F</given-names>
</name>
<article-title>Static microdroplet arrays: a microfluidic device for droplet trapping, incubation and release for enzymatic and cell-based assays</article-title>
<source>Lab Chip</source>
<year>2009</year>
<volume>9</volume>
<issue>5</issue>
<fpage>692</fpage>
<lpage>698</lpage>
<pub-id pub-id-type="doi">10.1039/b813709a</pub-id>
<pub-id pub-id-type="pmid">19224019</pub-id>
</mixed-citation>
</ref>
<ref id="B32"><mixed-citation publication-type="book"><name><surname>Atkinson</surname>
<given-names>AC</given-names>
</name>
<name><surname>Donev</surname>
<given-names>AN</given-names>
</name>
<name><surname>Tobias</surname>
<given-names>RD</given-names>
</name>
<source>Optimum Experimental Designs, with SAS</source>
<year>2007</year>
<publisher-name>Oxford: Oxford University Press</publisher-name>
</mixed-citation>
</ref>
<ref id="B33"><mixed-citation publication-type="journal"><name><surname>Hucka</surname>
<given-names>M</given-names>
</name>
<name><surname>Finney</surname>
<given-names>A</given-names>
</name>
<name><surname>Sauro</surname>
<given-names>HM</given-names>
</name>
<name><surname>Bolouri</surname>
<given-names>H</given-names>
</name>
<name><surname>Doyle</surname>
<given-names>JC</given-names>
</name>
<name><surname>Kitano</surname>
<given-names>H</given-names>
</name>
<name><surname>Arkin</surname>
<given-names>AP</given-names>
</name>
<name><surname>Bornstein</surname>
<given-names>BJ</given-names>
</name>
<name><surname>Bray</surname>
<given-names>D</given-names>
</name>
<name><surname>Cornish-Bowden</surname>
<given-names>A</given-names>
</name>
<article-title>The systems biology markup language (SBML): a medium for representation and exchange of biochemical network models</article-title>
<source>Bioinformatics</source>
<year>2003</year>
<volume>19</volume>
<issue>4</issue>
<fpage>524</fpage>
<lpage>531</lpage>
<pub-id pub-id-type="doi">10.1093/bioinformatics/btg015</pub-id>
<pub-id pub-id-type="pmid">12611808</pub-id>
</mixed-citation>
</ref>
<ref id="B34"><mixed-citation publication-type="other"><name><surname>Le Novere</surname>
<given-names>N</given-names>
</name>
<name><surname>Bornstein</surname>
<given-names>B</given-names>
</name>
<name><surname>Broicher</surname>
<given-names>A</given-names>
</name>
<name><surname>Courtot</surname>
<given-names>M</given-names>
</name>
<name><surname>Donizelli</surname>
<given-names>M</given-names>
</name>
<name><surname>Dharuri</surname>
<given-names>H</given-names>
</name>
<name><surname>Li</surname>
<given-names>L</given-names>
</name>
<name><surname>Sauro</surname>
<given-names>H</given-names>
</name>
<name><surname>Schilstra</surname>
<given-names>M</given-names>
</name>
<name><surname>Shapiro B Snoep</surname>
<given-names>JL</given-names>
</name>
<name><surname>Hucka</surname>
<given-names>M</given-names>
</name>
<article-title>BioModels Database: a free, centralized database of curated, published, quantitative kinetic models of biochemical and cellular systems</article-title>
<source>Nucleic Acids Res</source>
<year>2006</year>
<issue>34 Database</issue>
<fpage>D689</fpage>
<lpage>691</lpage>
<pub-id pub-id-type="doi">10.1093/nar/gkj092</pub-id>
<pub-id pub-id-type="pmid">16381960</pub-id>
</mixed-citation>
</ref>
<ref id="B35"><mixed-citation publication-type="other"><article-title>Plot Digitizer</article-title>
<ext-link ext-link-type="uri" xlink:href="http://plotdigitizer.sourceforge.net/">http://plotdigitizer.sourceforge.net/</ext-link>
</mixed-citation>
</ref>
</ref-list>
</back>
</pmc>
</record>
Pour manipuler ce document sous Unix (Dilib)
EXPLOR_STEP=$WICRI_ROOT/Sante/explor/MersV1/Data/Pmc/Corpus
HfdSelect -h $EXPLOR_STEP/biblio.hfd -nk 000B28 | SxmlIndent | more
Ou
HfdSelect -h $EXPLOR_AREA/Data/Pmc/Corpus/biblio.hfd -nk 000B28 | SxmlIndent | more
Pour mettre un lien sur cette page dans le réseau Wicri
{{Explor lien |wiki= Sante |area= MersV1 |flux= Pmc |étape= Corpus |type= RBID |clé= |texte= }}
![]() | This area was generated with Dilib version V0.6.33. | ![]() |