Integrated stoichiometric, thermodynamic and kinetic modelling of steady state metabolism.

The quantitative analysis of biochemical reactions and metabolites is at frontier of biological sciences. The recent availability of high-throughput technology data sets in biology has paved the way for new modelling approaches at various levels of complexity including the metabolome of a cell or an organism. Understanding the metabolism of a single cell and multi-cell organism will provide the knowledge for the rational design of growth conditions to produce commercially valuable reagents in biotechnology. Here, we demonstrate how equations representing steady state mass conservation, energy conservation, the second law of thermodynamics, and reversible enzyme kinetics can be formulated as a single system of linear equalities and inequalities, in addition to linear equalities on exponential variables. Even though the feasible set is non-convex, the reformulation is exact and amenable to large-scale numerical analysis, a prerequisite for computationally feasible genome scale modelling. Integrating flux, concentration and kinetic variables in a unified constraint-based formulation is aimed at increasing the quantitative predictive capacity of flux balance analysis. Incorporation of experimental and theoretical bounds on thermodynamic and kinetic variables ensures that the predicted steady state fluxes are both thermodynamically and biochemically feasible. The resulting in silico predictions are tested against fluxomic data for central metabolism in Escherichia coli and compare favourably with in silico prediction by flux balance analysis.

Eukaryotic DNA polymerases, a growing family.

In eukaryotic cells, DNA polymerases are required to maintain the integrity of the genome during processes, such as DNA replication, various DNA repair events, translesion DNA synthesis, DNA recombination, and also in regulatory events, such as cell cycle control and DNA damage checkpoint function. In the last two years, the number of known DNA polymerases has increased to at least nine (called alpha, beta, gamma, delta, epsilon, zeta, eta, t and iota), and yeast Saccharomyces cerevisiae contains REV1 deoxycytidyl transferase.

Human Cdc45 is a proliferation-associated antigen.

Cell division cycle protein 45 (Cdc45) plays a critical role in DNA replication to ensure that chromosomal DNA is replicated only once per cell cycle. We analysed the expression of human Cdc45 in proliferating and nonproliferating cells. Our findings show that Cdc45 protein is absent from long-term quiescent, terminally differentiated and senescent human cells, although it is present throughout the cell cycle of proliferating cells. Moreover, Cdc45 is much less abundant than the minichromosome maintenance (Mcm) proteins in human cells, supporting the concept that origin binding of Cdc45 is rate limiting for replication initiation. We also show that the Cdc45 protein level is consistently higher in human cancer-derived cells compared with primary human cells. Consequently, tumour tissue is preferentially stained using Cdc45-specific antibodies. Thus, Cdc45 expression is tightly associated with proliferating cell populations and Cdc45 seems to be a promising candidate for a novel proliferation marker in cancer cell biology.

Human DNA polymerase alpha gene: sequences controlling expression in cycling and serum-stimulated cells.

We have investigated the DNA polymerase alpha promoter sequence requirements for the expression of a heterologous gene in actively cycling cells and following serum addition to serum-deprived cells. An 11.4-kb genomic clone that spans the 5' end of this gene and includes 1.62 kb of sequence upstream from the translation start site was isolated. The transcription start site was mapped at 46 +/- 1 nucleotides upstream from the translation start site. The upstream sequence is GC rich and lacks a TATA sequence but has a CCAAT sequence on the opposite strand. Analysis of a set of deletion constructs in transient transfection assays demonstrated that efficient expression of the reporter in cycling cells requires 248 bp of sequence upstream from the cap site. Clustered within these 248 nucleotides are sequences similar to consensus sequences for Sp1-, Ap1-, Ap2-, and E2F-binding sites. The CCAAT sequence and the potential E2F- and Ap1-binding sites are shown to be protected from DNase I digestion by partially purified nuclear proteins. The DNA polymerase alpha promoter can confer upon the reporter an appropriate, late response to serum addition. No single sequence element could be shown to confer serum inducibility. Rather, multiple sequence elements appear to mediate the full serum response.

Species-specific replication of simian virus 40 DNA in vitro requires the p180 subunit of human DNA polymerase alpha-primase.

Human cell extracts efficiently support replication of simian virus 40 (SV40) DNA in vitro, while mouse cell extracts do not. Since human DNA polymerase alpha-primase is the major species-specific factor, we set out to determine the subunit(s) of DNA polymerase alpha-primase required for this species specificity. Recombinant human, mouse, and hybrid human-mouse DNA polymerase alpha-primase complexes were expressed with baculovirus vectors and purified. All of the recombinant DNA polymerase alpha-primases showed enzymatic activity and efficiently synthesized the complementary strand on an M13 single-stranded DNA template. The human DNA polymerase alpha-primase (four subunits [HHHH]) and the hybrid DNA polymerase alpha-primase HHMM (two human subunits and two mouse subunits), containing human p180 and p68 and mouse primase, initiated SV40 DNA replication in a purified system. The human and the HHMM complex efficiently replicated SV40 DNA in mouse extracts from which DNA polymerase alpha-primase was deleted, while MMMM and the MMHH complex did not. To determine whether the human p180 or p68 subunit was required for SV40 DNA replication, hybrid complexes containing only one human subunit, p180 or p68, together with three mouse subunits (HMMM and MHMM) or three human subunits and one mouse subunit (MHHH and HMHH) were tested for SV40 DNA replication activity. The hybrid complexes HMMM and HMHH synthesized oligoribonucleotides in the SV40 initiation assay with purified proteins and replicated SV40 DNA in depleted mouse extracts. In contrast, the hybrid complexes containing mouse p180 were inactive in both assays. We conclude that the human p180 subunit determines host-specific replication of SV40 DNA in vitro.

Two immunologically distinct human DNA polymerase alpha-primase subpopulations are involved in cellular DNA replication.

Metabolic labeling of primate cells revealed the existence of phosphorylated and hypophosphorylated DNA polymerase alpha-primase (Pol-Prim) populations that are distinguishable by monoclonal antibodies. Cell cycle studies showed that the hypophosphorylated form was found in a complex with PP2A and cyclin E-Cdk2 in G1, whereas the phosphorylated enzyme was associated with a cyclin A kinase in S and G2. Modification of Pol-Prim by PP2A and Cdks regulated the interaction with the simian virus 40 origin-binding protein large T antigen and thus initiation of DNA replication. Confocal microscopy demonstrated nuclear colocalization of hypophosphorylated Pol-Prim with MCM2 in S phase nuclei, but its presence preceded 5-bromo-2'-deoxyuridine (BrdU) incorporation. The phosphorylated replicase exclusively colocalized with the BrdU signal, but not with MCM2. Immunoprecipitation experiments proved that only hypophosphorylated Pol-Prim associated with MCM2. The data indicate that the hypophosphorylated enzyme initiates DNA replication at origins, and the phosphorylated form synthesizes the primers for the lagging strand of the replication fork.

The mouse DNA polymerase alpha-primase subunit p48 mediates species-specific replication of polyomavirus DNA in vitro.

Mouse cell extracts support vigorous replication of polyomavirus (Py) DNA in vitro, while human cell extracts do not. However, the addition of purified mouse DNA polymerase alpha-primase to human cell extracts renders them permissive for Py DNA replication, suggesting that mouse polymerase alpha-primase determines the species specificity of Py DNA replication. We set out to identify the subunit of mouse polymerase alpha-primase that mediates this species specificity. To this end, we cloned and expressed cDNAs encoding all four subunits of mouse and human polymerase alpha-primase. Purified recombinant mouse polymerase alpha-primase and a hybrid DNA polymerase alpha-primase complex composed of human subunits p180 and p68 and mouse subunits p58 and p48 supported Py DNA replication in human cell extracts depleted of polymerase alpha-primase, suggesting that the primase heterodimer or one of its subunits controls host specificity. To determine whether both mouse primase subunits were required, recombinant hybrid polymerase alpha-primases containing only one mouse primase subunit, p48 or p58, together with three human subunits, were assayed for Py replication activity. Only the hybrid containing mouse p48 efficiently replicated Py DNA in depleted human cell extracts. Moreover, in a purified initiation assay containing Py T antigen, replication protein A (RP-A) and topoisomerase I, only the hybrid polymerase alpha-primase containing the mouse p48 subunit initiated primer synthesis on Py origin DNA. Together, these results indicate that the p48 subunit is primarily responsible for the species specificity of Py DNA replication in vitro. Specific physical association of Py T antigen with purified recombinant DNA polymerase alpha-primase, mouse DNA primase heterodimer, and mouse p48 suggested that direct interactions between Py T antigen and primase could play a role in species-specific initiation of Py replication.

Cell cycle-dependent regulation of human DNA polymerase alpha-primase activity by phosphorylation.

DNA polymerase alpha-primase is known to be phosphorylated in human and yeast cells in a cell cycle-dependent manner on the p180 and p68 subunits. Here we show that phosphorylation of purified human DNA polymerase alpha-primase by purified cyclin A/cdk2 in vitro reduced its ability to initiate simian virus 40 (SV40) DNA replication in vitro, while phosphorylation by cyclin E/cdk2 stimulated its initiation activity. Tryptic phosphopeptide mapping revealed a family of p68 peptides that was modified well by cyclin A/cdk2 and poorly by cyclin E/cdk2. The p180 phosphopeptides were identical with both kinases. By mass spectrometry, the p68 peptide family was identified as residues 141 to 160. Cyclin A/cdk2- and cyclin A/cdc2-modified p68 also displayed a phosphorylation-dependent shift to slower electrophoretic mobility. Mutation of the four putative phosphorylation sites within p68 peptide residues 141 to 160 prevented its phosphorylation by cyclin A/cdk2 and the inhibition of replication activity. Phosphopeptide maps of the p68 subunit of DNA polymerase alpha-primase from human cells, synchronized and labeled in G1/S and in G2, revealed a cyclin E/cdk2-like pattern in G1/S and a cyclin A/cdk2-like pattern in G2. The slower-electrophoretic-mobility form of p68 was absent in human cells in G1/S and appeared as the cells entered G2/M. Consistent with this, the ability of DNA polymerase alpha-primase isolated from synchronized human cells to initiate SV40 replication was maximal in G1/S, decreased as the cells completed S phase, and reached a minimum in G2/M. These results suggest that the replication activity of DNA polymerase alpha-primase in human cells is regulated by phosphorylation in a cell cycle-dependent manner.

RPA subunit arrangement near the 3'-end of the primer is modulated by the length of the template strand and cooperative protein interactions.

To analyze the interaction of human replication protein A (RPA) and its subunits with the DNA template-primer junction in the DNA replication fork, we designed several template-primer systems differing in the size of the single-stranded template tail (4, 9, 13, 14, 19 and 31 nt). Base substituted photoreactive dNTP analogs-5-[ N -(2-nitro-5-azidobenzoyl)- trans -3-amino-propenyl-1]-2'-deoxyuridine-5'-triphosphate (NAB-4-dUTP) and 5-[ N -[ N -(2-nitro-5-azidobenzoyl)glycyl]- trans -3-aminopropenyl-1]-2'-deoxyuridine-5'-triphosphate (NAB-7-dUTP)-were used as substrates for elongation of radiolabeled primer-template by DNA polymerases in the presence or absence of RPA. Subsequent UV crosslinking showed that the pattern of p32 and p70 RPA subunit labeling, and consequently their interaction with the template-primer junction, is strongly dependent on the template extension length at a particular RPA concentration, as well as on the ratio of RPA to template concentration. Our results suggest a model of changes in the RPA configuration modulating by the length of the template extension in the course of nascent DNA synthesis.

A human topoisomerase I cleavage complex is recognized by an additional human topisomerase I molecule in vitro.

Several recent studies have shown that human topoisomerase I (htopoI) can recognize various DNA lesions and thereby form a covalent topoisomerase I-DNA complex, which is known to be detrimental to cells. We have investigated whether htopoI recognizes another htopoI that is covalently trapped on a DNA substrate. For this purpose we created an artificial DNA substrate containing a specific topoisomerase I binding sequence, where the enzyme was trapped in the covalently bound form. We demonstrate that, in vitro, free htopoI stimulates the formation of an additional cleavage complex immediately upstream of the covalently bound topoisomerase I. The predominant distance between the two cleavage sites is 13 nt. In addition we find that these two enzymes may form direct protein-protein contacts and we propose that these may be mediated through the formation of a dimer by domain swapping involving the C-terminal and the core domains. Finally, we discuss the possibility that the double cleavage reaction may be the initial step for the removal of the recognized cleavage complex.

Phosphorylation of DNA polymerase alpha-primase by cyclin A-dependent kinases regulates initiation of DNA replication in vitro.

DNA polymerase alpha-primase is the only known eukaryotic enzyme that can start DNA replication de novo. In this study, we investigated the regulation of DNA replication by phosphorylation of DNA polymerase alpha-primase. The p180 and the p68 subunits of DNA polymerase alpha-primase were phosphorylated using Cyclin A-, B- and E- dependent kinases. This phosphorylation did not influence its DNA polymerase activity on activated DNA, but slightly stimulated primase activity using poly(dT) single-stranded DNA (ssDNA) without changing the product length of primers. In contrast, site-specific initiation of replication on plasmid DNA containing the SV40 origin is affected: Cyclin A-Cdk2 and Cyclin A-Cdc2 inhibited initiation of SV40 DNA replication in vitro, Cyclin B-Cdc2 had no effect and Cyclin E-Cdk2 stimulated the initiation reaction. DNA polymerase alpha-primase that was pre-phosphorylated by Cyclin A-Cdk2 was completely unable to initiate the SV40 DNA replication in vitro; Cyclin B-Cdc2-phosphorylated enzyme was moderately inhibited, while Cyclin E-Cdk2-treated DNA polymerase alpha-primase remained fully active in the initiation reaction.

Surface plasmon resonance measurements reveal stable complex formation between p53 and DNA polymerase alpha.

Surface plasmon resonance measurements were used for detecting and quantifying protein-protein interactions between the tumor suppressor protein p53, the SV40 large T antigen (T-ag), the cellular DNA polymerase alpha-primase complex (pol-prim), and the cellular single-strand DNA binding protein RPA. Highly purified p53 protein bound to immobilized T-ag with an apparent binding constant of 2 x 10(8) M(-1). Binding of p53 to RPA was in the same order of magnitude with a binding constant of 4 x 10(8) M(-1), when RPA was coupled to the sensor chip via its smallest subunit, and 1 x 10(8) M(-1), when RPA was coupled via its p70 subunit. Furthermore, p53 bound human DNA polymerase alpha-primase complex (pol-prim) with a K(A) value of 1 x 10(10) m(-1). Both the p68 subunit and the p180 subunit of pol-prim could interact with p53 displaying binding constants of 2 x 10(10) m1(-1) and 5 X 10(9) M(-1), respectively. Complex formation was also observed with a p180/p68 heterodimer, and again with a binding constant similar. Hence, there was no synergistic effect when p53 bound to higher order complexes of pol-prim. A truncated form of p53, consisting of amino acids 1-320, bound pol-prim by four orders of magnitude less efficiently. Therefore, an intact C-terminus of p53 seems to be important for efficient binding to pol-prim. It was also tried to measure complex formation between p53, pol-prim, and T-ag. However there was no evidence for the existence of a ternary complex consisting of T-ag, pol-prim, and p53.

DNA polymerase alpha-DNA primase from human lymphoblasts.

The DNA polymerase alpha-DNA primase complex from the human lymphoblast line HSC93 has been enriched to near homogeneity by using an immunoaffinity purification protocol which was developed earlier for the purification of the calf thymus enzyme (Nasheuer, H.-P. and Grosse, F. (1987) Biochemistry 26, 8458-8466). Immunoaffinity purified polymerase-primase from human cells consisted of four subunits displaying molecular weights of 195,000 and 180,000 for the DNA synthesizing alpha-subunit, of 68,000 for the beta-subunit, and of 55,000 and 48,000 for the primase-carrying gamma- and delta-subunit, respectively. The isoelectric pH values for the individual subunits were estimated from non-equilibrium pH gradients to be between 5.9 and 5.7 for the alpha-subunit, at 5.5 for the beta-subunit, and at 7.5 and 8.0 for the gamma- and delta-subunit, respectively. The purified polymerase-primase converted single-stranded phi X174 DNA into the double-stranded form in a primase-initiated reaction. During this process, 3-10 RNA primers were formed. RNA primers were about 11 nucleotides long. Elongation of existing RNA primers by the human polymerase-primase was semi-processive; following primer binding the DNA polymerase continuously incorporated 20 to 50 nucleotides, then it dissociated from the template DNA.

Control of complex formation of DNA polymerase alpha-primase and cell-free DNA replication by the C-terminal amino acids of the largest subunit p180.

DNA polymerase alpha-primase is a heterotetrameric complex essential for simian vacuolating virus 40 (SV40) DNA replication. We show that the C-terminal 67 amino acid residues of the human p180 subunit are essential for SV40 DNA replication as they are required for binding of the p68 subunit and play a role in the interaction with the primase subunits, p48 and p58. Furthermore, we demonstrate that exchanging these residues to those of mouse origin can only partially rescue the SV40 DNA replication activity of DNA polymerase alpha-primase.

DNA polymerase alpha-DNA primase from human placenta. Immunoaffinity purification and preliminary characterization.

Highly purified DNA polymerase alpha-DNA primase from normal human tissue (human placenta) has been prepared by immunoaffinity purification on immobilized anti-human DNA polymerase alpha monoclonal antibody SJK 287-38. According to data from SDS electrophoresis this preparation consists of subunits of 180, 160, 145, 140 kDa (a cluster of DNA-polymerizing subunits), 73 kDa (function unknown) and 59, 52 kDa (corresponding to primase). Three active enzyme forms of 270, 460 and 575 kDa have been revealed using native electrophoresis followed by detection of DNA polymerase activity.

Alternative conformations of human replication protein A are detected by crosslinks with primers carrying a photoreactive group at the 3'-end.

To analyze the influence of single-stranded template extension of DNA duplex on the conformation of human replication protein A (RPA) bound to DNA we have designed two template-primer systems differing by the size of the single-stranded template tail (9 and 19 nucleotides (nt)). Base-substituted photoreactive dUTP analogs were used as substrates for elongation of radiolabeled template-primer by DNA polymerase beta in the absence or in the presence of RPA. Following UV-crosslinking it was demonstrated that the pattern of RPA subunit labeling and consequently RPA arrangement near the 3'-end of the primer is strongly dependent upon the length of the template extension.

Interaction of the p70 subunit of RPA with a DNA template directs p32 to the 3'-end of nascent DNA.

Human replication protein A is a heterotrimeric protein involved in various processes of DNA metabolism. To understand the contribution of replication protein A individual subunits to DNA binding, we have expressed them separately as soluble maltose binding protein fusion proteins. Using a DNA construct that had a photoreactive group incorporated at the 3'-end of the primer strand, we show that the p70 subunit on its own is efficiently cross-linked to the primer at physiological concentrations. In contrast, crosslinking of the p32 subunit required two orders of magnitude higher protein concentrations. In no case was the p14 subunit labelled above background. p70 seems to be the predominant subunit to bind single-stranded DNA and this interaction positions the p32 subunit to the 3'-end of the primer.

Quantitative assignment of reaction directionality in constraint-based models of metabolism: application to Escherichia coli.

Constraint-based modeling is an approach for quantitative prediction of net reaction flux in genome-scale biochemical networks. In vivo, the second law of thermodynamics requires that net macroscopic flux be forward, when the transformed reaction Gibbs energy is negative. We calculate the latter by using (i) group contribution estimates of metabolite species Gibbs energy, combined with (ii) experimentally measured equilibrium constants. In an application to a genome-scale stoichiometric model of Escherichia coli metabolism, iAF1260, we demonstrate that quantitative prediction of reaction directionality is increased in scope and accuracy by integration of both data sources, transformed appropriately to in vivo pH, temperature and ionic strength. Comparison of quantitative versus qualitative assignment of reaction directionality in iAF1260, assuming an accommodating reactant concentration range of 0.02-20mM, revealed that quantitative assignment leads to a low false positive, but high false negative, prediction of effectively irreversible reactions. The latter is partly due to the uncertainty associated with group contribution estimates. We also uncovered evidence that the high intracellular concentration of glutamate in E. coli may be essential to direct otherwise thermodynamically unfavorable essential reactions, such as the leucine transaminase reaction, in an anabolic direction.

Immunoaffinity-purified DNA polymerase alpha displays novel properties.

The purification and characterization of a novel and more intact form of the DNA polymerase alpha-primase complex from calf thymus are described. The polymerase-primase was enriched 10,000-fold to apparent homogeneity by chromatography on phosphocellulose, heparin-Sepharose, and an immobilized anti-human DNA polymerase alpha monoclonal antibody [SJK287-38; Tanaka, S., Hu, S., Wang, T. S.-F., & Korn, D. (1982) J. Biol. Chem. 257, 8386-8390]. A quantitative elution from the antibody column was achieved by shifting the pH from neutrality to between 12.5 and 13. From 1 kg of calf thymus, the procedure yields 1-2 mg of polymerase-primase with a specific activity of 30,000-40,000 units/mg for the polymerase and 15,000-20,000 units/mg for the primase. The complex sediments at 9 S through a sucrose gradient and exhibits a Stokes radius of 6.0 nm, yielding a native molecular mass of 335,000. Denaturing gel electrophoresis of the complex gives bands of Mrs 180,000, 155,000, 148,000, 73,000, 59,000, and 48,000 with a relative abundance of the two smallest subunits. Primase activity was partially resolved from the complex by centrifugation through sucrose gradients. The primer-forming activity was found to be associated with the Mr 59,000 and 48,000 polypeptides. In contrast to conventional preparations, the immunopurified polymerase displays several features which show it is the most intact form of the enzyme known to date. The deoxynucleoside triphosphate Km values are all within the range of 0.6-0.9 microM. The Km for binding to a single RNA primer on M13 DNA is 3.5 nM; the Ki for nonspecific binding to unprimed DNA is 70 microM (nucleotide).(ABSTRACT TRUNCATED AT 250 WORDS)

DNA polymerase alpha-primase from calf thymus. Determination of the polypeptide responsible for primase activity.

Immunoaffinity-purified DNA polymerase alpha-primase complex from calf thymus consists of subunits with molecular weights of 148,000-180,000, 73,000, 59,000, and 48,000 (Nasheuer, H.-P., and Grosse, F. (1987) Biochemistry 26, 8458-8466). Primase activity was separated from the immobilized complex by washing extensively with 2 M KCl or, alternatively, by shifting to pH 11.5 in the presence of 1 M KCl. From both elution procedures, the primase activity was found to be associated with the polypeptides with molecular weights of 59,000 and 48,000. The specific activity, using either elution procedure, was 30,000 units/mg. Both polypeptides sedimented together at 5.7 S upon zonal centrifugation on a sucrose gradient. Primase activity was found in the flow-through fraction after DEAE-cellulose chromatography of the free primase. Analysis of this fraction by sodium dodecyl sulfate gel electrophoresis revealed only one band with a Mr of 48,000. Polyclonal antibodies were raised against the Mr 59,000 and 48,000 polypeptides. The anti-Mr 59,000 antibody affected the primase activity only marginally, whereas the anti-Mr 48,000 antibody inhibited the primase activity nearly completely. UV cross-linking of the DNA polymerase alpha-primase complex with alpha-32P-labeled GTP revealed a binding site at the Mr 48,000 polypeptide, but none at the other subunits of the complex. Taken together, these results suggest that the Mr 48,000 polypeptide bears the active site of the DNA primase activity. The Mr 59,000 polypeptide stabilizes the primase activity.