Conformon-driven biopolymer shape changes in cell modeling

Conceptual models of the atom preceded the mathematical model of the hydrogen atom in physics in the second decade of the 20th century. The computer modeling of the living cell in the 21st century may follow a similar course of development. A conceptual model of the cell called the Bhopalator was formulated in the mid-1980s, along with its twin theories known as the conformon theory of molecular machines and the cell language theory of biopolymer interactions [Ann. N.Y. Acad. Sci. 227 (1974) 211; BioSystems 44 (1997) 17; Ann. N.Y. Acad. Sci. 870 (1999a) 411; BioSystems 54 (2000) 107; Semiotica 138 (1-4) (2002a) 15; Fundamenta Informaticae 49 (2002b) 147]. The conformon theory accounts for the reversible actions of individual biopolymers coupled to irreversible chemical reactions, while the cell language theory provides a theoretical framework for understanding the complex networks of dynamic interactions among biopolymers in the cell. These two theories are reviewed and further elaborated for the benefit of both computational biologists and computer scientists who are interested in modeling the living cell and its functions. One of the critical components of the mechanisms of cell communication and cell computing has been postulated to be space- and time-organized teleonomic (i.e. goal-directed) shape changes of biopolymers that are driven by exergonic (free energy-releasing) chemical reactions. The generalized Franck-Condon principle is suggested to be essential in resolving the apparent paradox arising when one attempts to couple endergonic (free energy-requiring) biopolymer shape changes to the exergonic chemical reactions that are catalyzed by biopolymer shape changes themselves. Conformons, defined as sequence-specific mechanical strains of biopolymers first invoked three decades ago to account for energy coupling in mitochondria, have been identified as shape changers, the agents that cause shape changes in biopolymers. Given a set of space- and time-organized teleonomic shape changes of biopolymers driven by conformons, all of the functions of the cell can be accounted for in molecular terms-at least in principle. To convert a conceptual model of the cell into a computer model, it is necessary to represent the conceptual model in an algebraic language. To this end, we have begun to apply the process algebra of Milner [Communicating and Mobile Systems: The pi-calculus, Cambridge University Press, Cambridge, 1999] to develop what is here called the "shape algebra," capable of describing complex and mobile patterns of interactions among biomolecules leading to cell functions.     

Xenopus egg lysates repair heat-generated DNA nicks with an average patch size of 36 nucleotides.

Base excision repair (BER) is an essential DNA repair pathway since it processes spontaneous (endogenous) DNA damage such as abasic sites, oxidized and alkylated bases, as well as mismatches arising from deamination of cytosine and 5-methylcytosine. Some of these lesions are repaired by the exchange of a single deoxynucleotide [Dianov, G. et al. (1992) Mol. Cell. Biol. 12, 1605-1612; Wiebauer, K. and Jiricny, J. (1990) Proc. Natl. Acad. Sci. USA, 87, 5842-5845] or a few deoxynucleotides [Matsumoto, Y. et al. (1994) Mol. Cell. Biol., 14 6187-6197]. Here we report that DNA single strand breaks induced by hyperthermic conditions are repaired with an average patch size of approximately 36 nt in Xenopus laevis egg lysates.     

Tubulin folding cofactor D is a microtubule destabilizing protein

A rapid switch between growth and shrinkage at microtubule ends is fundamental for many cellular processes. The main structural components of microtubules, the alphabeta-tubulin heterodimers, are generated through a complex folding process where GTP hydrolysis [Fontalba et al. (1993) J. Cell Sci. 106, 627-632] and a series of molecular chaperones are required [Sternlicht et al. (1993) Proc. Natl. Acad. Sci. USA 90, 9422-9426; Campo et al. (1994) FEBS Lett. 353, 162-166; Lewis et al. (1996) J. Cell Biol. 132, 1-4; Lewis et al. (1997) Trends Cell Biol. 7, 479-484; Tian et al. (1997) J. Cell Biol. 138, 821-823]. Although the participation of the cofactor proteins along the tubulin folding route has been well established in vitro, there is also evidence that these protein cofactors might contribute to diverse microtubule processes in vivo [Schwahn et al. (1998) Nature Genet. 19, 327-332; Hirata et al. (1998) EMBO J. 17, 658-666; Fanarraga et al. (1999) Cell Motil. Cytoskel. 43, 243-254]. Microtubule dynamics, crucial during mitosis, cellular motility and intracellular transport processes, are known to be regulated by at least four known microtubule-destabilizing proteins. OP18/Stathmin and XKCM1 are microtubule catastrophe-inducing factors operating through different mechanisms [Waters and Salmon (1996) Curr. Biol. 6, 361-363; McNally (1999) Curr. Biol. 9, R274-R276]. Here we show that the tubulin folding cofactor D, although it does not co-polymerize with microtubules either in vivo or in vitro, modulates microtubule dynamics by sequestering beta-tubulin from GTP-bound alphabeta-heterodimers.     

Energy coupling to active transport in anaerobically grown mutants of Escherichia Coli K12.

1. Anaerobic uptake of proline requires either the presence of a coupled Mg2+-stimulated adenosine triphosphatase or anaerobic electron transport. 2. Anaerobic uptake of glutamine does not require anaerobic electron transport even in the absence of a coupled Mg+2-stimulated adenosine triphosphatase. 3. These results support previous suggestions [Berger (1973) Proc. Natl. Acad. Sci. U.S.A. 70, 1514--1518; Berger & Heppel (1974) J. Biol. Chem. 249, 7747-7755; Kobayashi, Kin & Anraku (1974) J. Biochem. (Tokyo) 76, 251-261] that two distinct mechanisms of energy coupling to active transport exist in Escherichia coli in that energization of anaerobic proline uptake requires the 'high-energy membrane state', whereas the energization of anaerobic glutamine uptake does not.     

Horizontal gene transfer in microbial genome evolution

Horizontal gene transfer is the collective name for processes that permit the exchange of DNA among organisms of different species. Only recently has it been recognized as a significant contribution to inter-organismal gene exchange. Traditionally, it was thought that microorganisms evolved clonally, passing genes from mother to daughter cells with little or no exchange of DNA among diverse species. Studies of microbial genomes, however, have shown that genomes contain genes that are closely related to a number of different prokaryotes, sometimes to phylogenetically very distantly related ones. (Doolittle et al., 1990, J. Mol. Evol. 31, 383-388; Karlin et al., 1997, J. Bacteriol. 179, 3899-3913; Karlin et al., 1998, Annu. Rev. Genet. 32, 185-225; Lawrence and Ochman, 1998, Proc. Natl. Acad. Sci. USA 95, 9413-9417; Rivera et al., 1998, Proc. Natl. Acad. Sci. USA 95, 6239-6244; Campbell, 2000, Theor. Popul. Biol. 57 71-77; Doolittle, 2000, Sci. Am. 282, 90-95; Ochman and Jones, 2000, Embo. J. 19, 6637-6643; Boucher et al. 2001, Curr. Opin., Microbiol. 4, 285-289; Wang et al., 2001, Mol. Biol. Evol. 18, 792-800). Whereas prokaryotic and eukaryotic evolution was once reconstructed from a single 16S ribosomal RNA (rRNA) gene, the analysis of complete genomes is beginning to yield a different picture of microbial evolution, one that is wrought with the lateral movement of genes across vast phylogenetic distances. (Lane et al., 1988, Methods Enzymol. 167, 138-144; Lake and Rivera, 1996, Proc. Natl. Acad. Sci. USA 91, 2880-2881; Lake et al., 1999, Science 283, 2027-2028).     

Drosophila chk2 plays an important role in a mitotic checkpoint in syncytial embryos

Drosophila chk2 (Dmchk2, also called Dmnk) plays a crucial role in the DNA damage response pathway mediating cell cycle arrest and apoptosis [Xu et al., FEBS Lett. 508 (2001) 394-398; Peters et al., Proc. Natl. Acad. Sci. USA 99 (2002) 11305-11310]. In this study, the role of Dmchk2 in early embryogenesis was investigated. In the absence of Dmchk2 function, abnormal nuclei accumulate in the cortex of the syncytial embryo. We show that the abnormal nuclei result from a failure of chromosome segregation probably due to damaged or incomplete replicated DNA. Importantly, this Dmchk2 phenotype is partially suppressed by reducing the gene dosage of polo or stg. As Polo-like kinase was shown to colocalize and coimmunoprecipitate with Chk2 [Tsvetkov et al., J. Biol. Chem. 278 (2003) 8468-8475] in mammals, these observations suggest that polo might be a key target of Dmchk2 in regulating mitotic entry in response to DNA damage or replication block.     

The complete amino-acid sequence of a trout-testis non-histone protein, H6, localized in a subset of nucleosomes and its similarity to calf-thymus non-histone proteins HMG-14 and HMG-17.

The complete amino acid sequence of a basic non-histone protein, H6, isolated from the chromatin of rainbow trout (Salmo gairdnerii) testis cells, has been determined. Protein H6, first described by D. T. Wigle and G. H. Dixon [J. Biol. Chem. 246, 5636--5644 (1971)] was extracted with 5% trichloracetic acid and purified by ion-exchange chromatography on carboxymethyl-cellulose (CM-52). Sequence analysis was performed by automatic Edman degradation of the amino terminus of the intact protein and a series of large fragments derived by cleavage with chymotrypsin, staphylococcal protease and with mild acid to cleave at aspartic acid residues. Protein H6 possesses 69 residues and shows considerable similarities to the 89-residue calf thymus HMG-17 protein previously sequenced [Walker, J. M., Hastings, J. R. B. & Johns, E. W. (1977) Eur. J. Biochem. 76, 461--468]. B. Levy W. and G. H. Dixon [Proc. Natl Acad. Sci. U.S.A. 74, 2810--2814 (1977)] have shown that H6 is selectively solubilized when trout testis nuclei (or chromatin) are digested with DNase I under conditions which preferentially hydrolyze that portion of DNA enriched in transcribed sequences [Levy, W. B. & Dixon, G. H. (1977) Nucleic Acids Res. 4, 883--898]. Recently H6 has been located as a stoichiometric component of a distinct subset of trout testis nucleosomes that are complexed with a core nucleosome comprising 140 base pairs of DNA and the inner histones H2A, H2B, H3 and H4 [Levy, W. B., Connor, W. & Dixon, G. H. (1979) J. Biol. Chem., in the press].     

Electrostatic switches that mediate the pH-dependent conformational change of "short" recombinant human pseudocathepsin D

Human cathepsin D (hCatD) is an aspartic peptidase with a low pH optimum. X-ray crystal structures have been solved for an active, low pH (pH 5.1) form (CatD(lo)) [Baldwin, E. T., Bhat, T. N., Gulnik, S., Hosur, M. V., Sowder, R. C., Cachau, R. E., Collins, J., Silva, A. M., and Erickson, J. W. (1993) Proc. Natl. Acad. Sci. U.S.A. 90, 6796-6800] and an inactive, high pH (pH 7.5) form (CatD(hi)) [Lee, A. Y., Gulnik, S. V., and Erickson, J. W. (1998) Nat. Struct. Biol. 5, 866-871]. It has been suggested that ionizable switches involving the carboxylate side chains of E5, E180, and D187 may mediate the reversible interconversion between CatD(hi) and CatD(lo) and that Y10 stabilizes CatD(hi) [Lee, A. Y., Gulnik, S. V., and Erickson, J. W. (1998) Nat. Struct. Biol. 5, 866-871]. To test these hypotheses, we generated single point mutants in "short" recombinant human pseudocathepsin D (srCatD), a model kinetically similar to hCatD [Beyer, B. M., and Dunn, B. M. (1996) J. Biol. Chem. 271, 15590-15596]. E180Q, Y10F, and D187N exhibit significantly higher kcat/Km values (2-, 3-, and 6-fold, respectively) at pH 3.7 and 4.75 compared to srCatD, indicating that these residues are important in stabilizing the CatD(hi). E5Q exhibits a 2-fold lower kcat/Km compared to srCatD at both pH values, indicating the importance of E5 in stabilizing the CatD(lo). Accordingly, full time-course "pH-jump" (pH 5.5-4.75) studies of substrate hydrolysis indicate that E180Q, D187N, and Y10F have shorter kinetic lag phases that represent the change from CatD(hi) to CatD(lo) compared to srCatD and E5Q. Intrinsic tryptophan fluorescence reveals that the variants have a native-like structure over the pH range of our assays. The results indicate that E180 and D187 participate as an electrostatic switch that initiates the conformational change of CatD(lo) to CatD(hi) and Y10 stabilizes CatD(hi) by hydrogen bonding to the catalytic Asp 33. E5 appears to play a less significant role as an ionic switch that stabilizes CatD(lo).     

A simple and efficient method for the separation and detection of small DNA fragments by electrophoresis in formamide containing agarose gels and Southern blotting to DBM-paper.

We show in this report that DNA fragments smaller than 300 bp are separated with high resolution by electrophoresis in concentrated (up to 7%) agarose gels containing 50% formamide. The separated DNA fragments can subsequently be quantitatively transferred to DBM-paper [Alwine, J.C. et al., Proc. Natl. Acad. Sci. USA (1977) 74, 5350-54] using the Southern technique [Southern, E.M., J. Mol. Biol. (1975) 98, 503-517], while preserving the sharpness of the original gel pattern. Since thin (0.2-0.4 mm) and thick (up to 5 mm) agarose slab gels can be easily handled in vertical or horizontal apparatus, this method should prove to be a very useful extention of the Southern technique, applicable to a variety of analytical and preparative purposes.     

Involvement of phenylalanine 272 of DNA polymerase beta in discriminating between correct and incorrect deoxynucleoside triphosphates

DNA polymerase beta is a small monomeric polymerase that participates in base excision repair and meiosis [Sobol, R., et al. (1996) Nature 379, 183-186; Plug, A., et al. (1997) Proc. Natl. Acad. Sci. U.S.A. 94, 1327-1331]. A DNA polymerase beta mutator mutant, F272L, was identified by an in vivo genetic screen [Washington, S., et al. (1997) Proc. Natl. Acad. Sci. U.S.A. 94, 1321-1326]. Residue 272 is located within the deoxynucleoside triphosphate (dNTP) binding pocket of DNA polymerase beta according to the known DNA polymerase beta crystal structures [Pelletier, H., et al. (1994) Science 264, 1891-1893; Sawaya, M., et al. (1997) Biochemistry 36, 11205-11215]. The F272L mutant produces errors at a frequency 10-fold higher than that of wild type in vivo and in the in vitro HSV-tk gap-filling assay. F272L shows an increase in the frequency of both base substitution mutations and frameshift mutations. Single-enzyme turnover studies of misincorporation by wild type and F272L DNA polymerase beta demonstrate that there is a 4-fold decrease in fidelity of the mutant as compared to that of the wild type enzyme for a G:A mismatch. The decreased fidelity is due primarily to decreased discrimination between the correct and incorrect dNTP during ground-state binding. These results suggest that the phenylalanine 272 residue is critical for maintaining fidelity during the binding of the dNTP.     

The DNA-binding protein Hdf1p (a putative Ku homologue) is required for maintaining normal telomere length in Saccharomyces cerevisiae.

In mammalian cells, the Ku autoantigen is an end- binding DNA protein required for the repair of DNA breaks [Troelstra, C. and Jaspers, N.G.J. (1994) Curr. Biol., 4, 1149- 1151]. A yeast gene (HDF1) encoding a putative homologue of the 70 kDa subunit of Ku has recently been identified [Feldmann, H. and Winnacker, E. L. (1993) J. Biol. Chem., 268, 12895- 12900]. We find that hdf1 mutant strains have substantially shorter telomeres than wild-type strains. We speculate that Hdf1p may bind the natural ends of the chromosome, in addition to binding to the ends of broken DNA molecules. Strains with both an hdf1 mutation and a mutation in TEL 1 (a gene related to the human ataxia telangiectasia gene) have extremely short telomeres and grow slowly.     

Selection of a cDNA clone which contains the complete coding sequence for the mature form of ornithine transcarbamylase from rat liver: expression of the cloned protein in Escherichia coli. Molecular cloning of rat ornithine transcarbamylase

A cDNA clone corresponding to the mature form of ornithine transcarbamylase (OTCase) was selected from a rat liver cDNA library constructed in bacteriophage lambda gt10. OTCase clones were selected using a synthetic DNA probe of 15 bases corresponding to the 3' end of the OTCase mRNA [Horwich, A. L., Kraus, J.P., Williams, K., Kalousek, F., K?nigsberg, W. & Rosenberg, L.E. (1983) Proc. Natl Acad. Sci. USA, 80, 4258-4262]. Putative OTCase clones were subcloned into the expression vector, pUC9, and the identity of inserts confirmed by colony immunoassay and by electrophoretic transfer of cloned proteins from sodium dodecyl sulphate/polyacrylamide gels to nitrocellulose filters followed by probing with monospecific anti-OTCase antibodies and 125I-labelled protein A. A clone corresponding to the full-length mature form of rat liver OTCase (plus 15 amino acids from Escherichia coli beta-galactosidase) was obtained and the identity of the clone was confirmed by comparison of the 5' sequence with a limited N-terminal amino acid sequence [Lusty, C., Jilka, R. L. & Nietsch, E. H. (1979) J. Biol. Chem. 254, 10030-10036]. A sequence discrepancy between the published sequence (Lusty et al.) and the sequence predicted from the cDNA structure is noted.     

Characterization of two forms of mouse salivary androgen-binding protein (ABP): implications for evolutionary relationships and ligand-binding function

Mouse salivary androgen-binding protein (ABP) is a member of the secretoglobin family produced in the submaxillary glands of house mice (Mus musculus). We report the cDNA sequences and amino acid sequences of the beta and gamma subunits of ABP from a mouse cDNA library, identifying the two subunits by their pIs and molecular weights. An anomalously high molecular weight of the alpha subunit is likely due to glycosylation at a single site. A phylogenetic comparison of the three subunits of ABP with the chains of other mammalian secretoglobins shows that ABP is most closely related to mouse lachrymal protein and to the major cat allergen Fel dI. An evaluation of the most conserved residues in ABP and the other secretoglobins, in light of structural data reported by others [Callebaut, I., Poupon, A., Bally, R., Demaret, J.-P., Housset, D., Delettre, J., Hossenlopp, P., and Mornon, J.-P. (2000) Ann. N.Y. Acad. Sci. 923, 90-112; Pattabiraman, N., Matthews, J., Ward, K., Mantile-Selvaggi, G., Miele, L., and Mukherjee, A. (2000) Ann. N.Y. Acad. Sci. 923, 113-127], allows us to draw conclusions about the critical residues important in ligand binding by the two different ABP dimers and to assess the importance of ligand binding in the function of the molecule. In addition to the cDNAs, which represent those of the musculus subspecies of Mus musculus, we also report the coding regions of the beta and gamma subunit cDNAs from two other mouse inbred strains which represent the other two subspecies: M. musculus domesticus and M. musculus castaneus. The high nonsynonymous/synonymous substitution rate ratios (K(a)/K(s)) for both the beta and gamma subunits suggest that these two proteins are evolving under strong directional selection, as has been reported for the alpha subunit [Hwang, J., Hofstetter, J., Bonhomme, F., and Karn, R. (1997) J. Hered. 88, 93-97; Karn, R., and Clements, M. (1999) Biochem. Genet. 37, 187-199].     

The yeast DNA repair proteins RAD1 and RAD7 share similar putative functional domains

Sequence information on eukaryotic DNA repair proteins provided so far only few clues concerning possible functional domains. Since the DNA repair process involves a strict sequential complex formation of several proteins [1988) FASEB J. 2, 2696-2701], we searched for special protein-protein interacting domains, which consist of tandemly repeated leucine rich motifs (LRM). Search algorithms, capable of detecting even largely divergent repeats by assessing their significance due to the tandem repetitivity, revealed that the yeast DNA repair proteins RAD1 and RAD7 contain 9 and 12 tandem LRM repeats, respectively. These results represent the first clues concerning specific domains in these proteins and assign them to the LRM superfamily, which includes such members as yeast adenylate cyclase, cell surface protein receptors and ribonuclease/angiogenin inhibitor, all exerting their function by specific protein-protein interactions involving LRM domains [( 1988) EMBO J. 7, 4151-4156; (1990) Proc. Natl. Acad. Sci. USA 87, 8711-8715; (1989) Science 245, 494-499; (1990) Mol. Cell. Biol. 10, 6436-6444; (1989) Proc. Natl. Acad. Sci. USA 86, 6773-6777].     

Phased psoralen cross-links do not bend the DNA double helix

Although the chemical reaction of psoralens with nucleic acids is well understood, the structure of psoralen-DNA cross-linked products is still not clear. Model building studies base on the crystal structure of the psoralen-thymine monoadduct suggest that each cross-link bends the DNA double helix by 46.5 degrees [Pearlman, D. A., Holbrook, S. R., Pirkle, D. H., & Kim, S.-H. (1985) Science (Washington, D.C.) 227, 1304-1308]. On the other hand, Sinden and Hagerman [Sinden, R. R., & Hagerman, P. J. (1984) Biochemistry 23, 6299-6303] find that, in solution, psoralen cross-linked DNA is not bent. Here we use gel electrophoresis to test the validity of the current models. We have synthesized a series of DNA fragments (21-24 base pairs in length), each containing one unique T-A site for 4'-(hydroxymethyl)-4,5',8-trimethylpsoralen (HMT) cross-linking. Because of an estimated 28 degrees unwinding of the helix by HMT [Wiesehahn, G., & Hearst, J. E. (1978) Proc. Natl. Acad. Sci. U.S.A. 75, 2703-2707], one expects that the 22-bp cross-linked fragment will be repeated nearly in phase with the average helical screw when multimerized. In that sequence ligation will maximally amplify any deformation to the double helix. We find that the ligated multimers of cross-linked DNA migrate close to the multimers of non-cross-linked DNA on polyacrylamide gels. Our observations place an upper limit of 10 degrees on DNA bending induced by psoralen cross-linking and indicate unwinding by about 1 bp, as well as stiffening of the double helix. These properties are not unexpected for classical intercalators.     

Sequence similarities among monkey ori-enriched (ors) fragments

Nucleotide sequences have been determined for eight ors (ori-enriched sequence) fragments isolated from monkey DNA by a method that was designed to enrich for origins of DNA replication [Kaufmann et al., Mol. Cell. Biol. 5 (1985) 721-727]. Evidence has been presented that some or possibly all of these sequences can serve, albeit inefficiently, as oris in vivo [Frappier and Zannis-Hadjopoulos, Proc. Natl. Acad. Sci. USA 84 (1987) 6668-6672]. Two of the fragments were found to contain the long terminal repeat-like elements of the 'O-family' of moderately repetitive sequences that are present in human DNA as a transposon-like element [Paulson et al., Nature 315 (1985) 359-361]. Extensive pair-wise comparisons of the sequences failed to detect any statistically significant common sequences, except for long asymmetrically distributed A + T-rich stretches. Nonetheless, when the ors fragments were examined for the presence of published consensus sequences, seven of eight were found to contain the control sequence described by Dierks et al. [Cell 32 (1983) 695-706], and the same seven of eight were found to contain both the scaffold attachment region T consensus [Gasser and Laemmli, Cell 46 (1986) 521-530] and the minimal Saccharomyces cerevisiae autonomously replicating sequence consensus [e.g., Palzkill and Newlon, Cell 53 (1988) 441-450].