An Interactive Dot Matrix System for Locating Potentially Significant Features in Nucleic Acid Molecules

Abstract

An interactive computer system using a dot matrix approach has been developed and used to determine potentially significant features due to distortions in the B-DNA helix as a result of variations of purine and pyrimidine patterns. Sequences were compared using matrices which were generated using the Calladine-Dickerson rules (C.R. Calladine, J. Mol. Biol. 161, 343–352, 1982 and R.E. Dickerson, J. Mol. Biol. 766, 419–441, 1983). Having control over various parameters to enhance different aspects of the visual appearance of these matrices was helpful in discovering patterns that were not known a priori. Specifically, it was found that a pattern of alternating doublets of purines and pyrimidines appear to exist in regulatory regions. This event is shown to be beyond probabilistic expectation.     

Transgenic analysis of a hybrid poplar wound-inducible promoter reveals developmental patterns of expression similar to that of storage protein genes.

The wound-inducible win3 multigene family from hybrid poplars (Populus trichocarpa x Populus deltoides) encodes proteins with structural similarities with Kunitz-type protease inhibitors (H.D. Bradshaw Jr., J.B. Hollick, T.J. Parsons, H.R.G. Clarke, M.P. Gordon [1990] Plant Mol Biol 14: 51-59), and at least one member, win3.12, is transcribed de novo in the injured and uninjured leaves of wounded trees (J.B. Hollick, M.P. Gordon [1993] Plant Mol Biol 22: 561-572). A previous study demonstrated that 1352 bp of 5' flanking DNA from the win3.12 gene confers local wound-regulated expression of the beta-glucuronidase gene in transgenic tobacco (Nicotiana tabacum cv Xanthi n.c.) (J.B. Hollick, M.P. Gordon [1993] Plant Mol Biol 22: 561-572). We extend this transgenic analysis here by examining the developmental regulation and systemic wound induction conferred by the same transgene construct in tobacco. Biochemical and histochemical surveys of beta-glucuronidase activity are described for four, independent transgenic lines. The observed spatial and temporal expression patterns coincide with dormant storage tissues and with previously described expression patterns for both seed and vegetative storage protein genes. Developmental northern blot analysis of win3 RNA levels in poplar seeds confirms that proper temporal expression of the reporter gene is maintained during tobacco seed maturation. These results demonstrate that a putative Kunitz-type protease inhibitor can be wound inducible in addition to being expressed in developing seeds.     

A comparison of muscle thin filament models obtained from electron microscopy reconstructions and low-angle X-ray fibre diagrams from non-overlap muscle

The regulation of striated muscle contraction involves changes in the interactions of troponin and tropomyosin with actin thin filaments. In resting muscle, myosin-binding sites on actin are thought to be blocked by the coiled-coil protein tropomyosin. During muscle activation, Ca2+ binding to troponin alters the tropomyosin position on actin, resulting in cyclic actin-myosin interactions that accompany muscle contraction. Evidence for this steric regulation by troponin-tropomyosin comes from X-ray data [Haselgrove, J.C., 1972. X-ray evidence for a conformational change in the actin-containing filaments of verterbrate striated muscle. Cold Spring Habor Symp. Quant. Biol. 37, 341-352; Huxley, H.E., 1972. Structural changes in actin and myosin-containing filaments during contraction. Cold Spring Habor Symp. Quant. Biol. 37, 361-376; Parry, D.A., Squire, J.M., 1973. Structural role of tropomyosin in muscle regulation: analysis of the X-ray diffraction patterns from relaxed and contracting muscles. J. Mol. Biol. 75, 33-55] and electron microscope (EM) data [Spudich, J.A., Huxley, H.E., Finch, J., 1972. Regulation of skeletal muscle contraction. II. Structural studies of the interaction of the tropomyosin-troponin complex with actin. J. Mol. Biol. 72, 619-632; O'Brien, E.J., Gillis, J.M., Couch, J., 1975. Symmetry and molecular arrangement in paracrystals of reconstituted muscle thin filaments. J. Mol. Biol. 99, 461-475; Lehman, W., Craig, R., Vibert, P., 1994. Ca2+-induced tropomyosin movement in Limulus thin filaments revealed by three-dimensional reconstruction. Nature 368, 65-67] each with its own particular strengths and limitations. Here we bring together some of the latest information from EM analysis of single thin filaments from Pirani et al. [Pirani, A., Xu, C., Hatch, V., Craig, R., Tobacman, L.S., Lehman, W. (2005). Single particle analysis of relaxed and activated muscle thin filaments. J. Mol. Biol. 346, 761-772], with synchrotron X-ray data from non-overlapped muscle fibres to refine the models of the striated muscle thin filament. This was done by incorporating current atomic-resolution structures of actin, tropomyosin, troponin and myosin subfragment-1. Fitting these atomic coordinates to EM reconstructions, we present atomic models of the thin filament that are entirely consistent with a steric regulatory mechanism. Furthermore, fitting the atomic models against diffraction data from skinned muscle fibres, stretched to non-overlap to preclude crossbridge binding, produced very similar results, including a large Ca2+-induced shift in tropomyosin azimuthal location but little change in the actin structure or apparent alteration in troponin position.     

Structural transition of alpha 1-antitrypsin by a peptide sequentially similar to beta-strand s4A

Crystal structure studies have shown that cleaved and intact serpins differ essentially in the topology of beta-sheet A. This is five-stranded in the intact molecules and six-stranded after cleavage by insertion of strand s4A whose C-terminus has become free [L?bermann, H., Tokuoka, R., Deisenhofer, J. & Huber, R. (1984) J. Mol. Biol. 177, 531-556; Wright, T. H., Qian, H. X. & Huber, R. (1990) J. Mol. Biol. 213, 513-528]. The structural transition is accompanied by changes in spectral properties and an increase in thermal stability. We show here that an N alpha-acetyl-tetradecapeptide with the amino acid sequence of strand s4A, residues 345-358 of human alpha 1-antitrypsin, associates with intact alpha 1-antitrypsin and forms a stoichiometric complex with properties very similar to cleaved alpha 1-antitrypsin. Complex generation has the characteristics of a folding process.     

Hydrogen exchange of individual amide protons in the E. coli lac repressor DNA-binding domain: a nuclear magnetic resonance study

Proton exchange in lac repressor headpiece was studied by COSY and 2D NOE spectroscopy. The exchange rates of amide protons, stabilized by the hydrogen bonds of the three alpha-helices of the headpiece, could be determined quantitatively. The exchange rates in these helices showed repetitive patterns of about three to four residues. A correlation with the position of the amide proton in the interior or the exterior of the alpha-helix of the protein was found. The exchange data strongly support the validity of the three-dimensional structure, as determined recently (Kaptein, R. et al., J. Mol. Biol. 182, 179-182 (1985)).     

Analysis of the high- and low-spin Soret bands of horse-heart metmyoglobin complexes

Interest in the thermodynamics of the iron-binding site in hemoproteins has increased in recent years due to refinements in x-ray crystallographic studies of hemoproteins [see Deathage, J. F., Lee, R. S., Anderson, C. M. & Moffat, K. (1976) J. Mol. Biol. 104 , 687–706; Heidner, E. J., Ladner, R. C. & Perutz, M. F. (1976) J. Mol. Biol. 104 , 707–722; Deathage, J. F., Lee, R. S. & Moffat, K. (1976) J. Mol. Biol. 104 , 723–728; Ladner, R. C., Heidner, E. J. & Perutz, M. F. (1976) J. Mol. Biol. 114 , 385–414; Fermi, G. & Perutz, M. F. (1977) J. Mol. Biol. 114 , 421–431; Takano, T. (1977) J. Mol. Biol. 110 , 537–568 and 569–589], the synthesis and x-ray analysis of model heme compounds [see Scheidt, W. R. (1977) Acc. Chem. Res. 10 , 339–345; Kastner, M. E., Scheidt, W. R., Mashino, T. & Reed, C. A. (1978) J. Am. Chem. Soc. 100 , 666–667; Mashiko, T., Kastner, M. E., Spartalian, K., Scheidt, W. R. & Reed, C. A. (1978) J. Am. Chem. Soc. 100 , 6354–6362; Hill, H. A. O., Skite, P. P., Buchler, J. W., Luchr, H., Tonn, M., Gregson, A. K. & Pellizer, G. (1979) Chem. Commun. 4 , 151–152; and Scheidt, W. R., Cohen, I. A. & Kastner, M. E. (1979) Biochemistry 18 , 3546–3556], and the numerous data on heme–protein interactions that account for the differences observed in ligand binding between the various species of animals. Numerous probes have been used and provide information about the structure and thermodynamics of the binding site, but no single probe can provide the complete picture [see Iizuka, T. & Yonetani, T. (1970) Adv. Biophys. 1 , 157–182; Smith, D. W. & Williams, R. J. P. (1970) Struct. Bond. 7 , 1–45; and Spiro, T. G. (1975) Biochim. Biophys. Acta 416 , 169–189].     

Structural instability of IncP-1 plasmids in Pseudomonas aeruginosa PAT involves interaction with plasmid pVS1

The structural instability exhibited by IncP-1 plasmids in Pseudomonas aeruginosa strain PAT was shown to be Rec+ dependent and involved interaction with the resident plasmid pVS1. Structural instability resulted from deletion of plasmid deoxyribonucleic acid at a frequency of ca. 10(-2)/cell per generation. Deletants could be stabilized by transduction into P. aeruginosa strain PAO, but in strain PAT deletants had only a transient existence, as continued deletion led eventually to the loss of the entire plasmid. The patterns of markers lost in PAT were used to demonstrate a marker order for R68 similar to that published elsewhere for RP4 (Barth and Grinter, J. Mol. Biol. 113:455-474, 1977), except that only one Tra region was found. R68 also exhibited Rec+-dependent structural instability in PAO(pVS1) derivatives but, unlike the case in PAT, instability was not accompanied by chromosome mobilization. We isolated deletants of pVS1 which were unable to promote structural instability.     

Folding of homologous proteins: conservation of the folding mechanism of the alpha subunit of tryptophan synthase from Escherichia coli, Salmonella typhimurium, and five interspecies hybrids

The equilibrium and kinetic properties for the urea-induced unfolding of the alpha subunit of tryptophan synthase from Escherichia coli, Salmonella typhimurium, and five interspecies hybrids were compared to determine the role of protein folding in evolution. The parent proteins differ at 40 positions in the sequence of 268 amino acids, and the hybrids differ by up to 15 amino acids from the Escherichia coli alpha subunit. The results show that all the proteins follow the same folding mechanism and are consistent with a previously proposed hypothesis [Hollecker, M., & Creighton, T. E. (1983) J. Mol. Biol. 168, 409; Krebs, H., Schmid, F. X., & Jaenicke, R. (1983) J. Mol. Biol. 169, 619] that the folding mechanisms are conserved in homologous proteins. Analysis of the kinetic data suggests that the 15 positions at which the parent proteins differ in the amino folding unit, residues 1-188, do not play a role in a rate-limiting step in folding that has been previously identified as the association of the amino and carboxyl folding units [Beasty, A. M., Hurle, M. R., Manz, J. T., Stackhouse, T. S., Onuffer, J. J., & Matthews, C. R. (1986) Biochemistry 25, 2965]. One or more of the 25 positions at which the parent proteins differ in the carboxyl folding unit, residues 189-268, do appear to play a role in this same rate-limiting step.     

The complete amino-acid sequence of the bilin-binding protein from Pieris brassicae and its similarity to a family of serum transport proteins like the retinol-binding proteins

The amino-acid sequence from the bilin binding protein (BBP) of the butterfly Pieris brassicae has been determined. The apoprotein with a length of 173 amino-acid residues has a molecular mass of 19,676 Da. The sequence analysis was performed by automated Edman degradation of the intact apoprotein and of fragments as large as possible generated from different digestions. The 3-dimensional structure of BBP, determined by Huber et al. (Huber, R., Schneider, M., Epp, O., Mayr, I., Messerschmidt, A., Pflugrath, J. & Kayser, H. (1987) J. Mol. Biol. 195, 423-434 and Huber, R., Schneider, M., Mayr, I., Müller, R., Deutzmann, R., Suter, F., Zuber, H., Falk, H. & Kayser, H. (1987) J. Mol. Biol. 198, 499-513) down to 2-A resolution, exhibits a similar conformation to the human retinol binding protein. Sawyer (Sawyer, L. (1987) Nature (London) 327, 659) demonstrated that proteins from a wide variety of sources can be gathered into a "superfamily". Computer searches of data banks yielded in a new member of this superfamily, namely human alpha 1-acid glycoprotein. One of the functions of the listed proteins is to bind and transport small hydrophobic molecules in serum.     

The TGN38 glycoprotein contains two non-overlapping signals that mediate localization to the trans-Golgi network

The membrane-spanning and cytoplasmic domains of CD4 and CD8 were replaced by those of TGN38. After transient expression in HeLa cells, the location of the hybrid proteins was determined using immunofluorescence and quantitative immuno-electron microscopy, FACS analysis and metabolic labeling. The membrane-spanning domain was found to contain a signal that localized hybrid proteins to the TGN. This was in addition to the signal previously identified in the cytoplasmic domain (Bos, K., C. Wraight, and K. Stanley. 1993. EMBO (Eur. Mol. Biol. Organ) J. 12:2219-2228. Humphrey, J. S., P. J. Peters, L. C. Yuan, and J. S. Bonifacino. 1993. J. Cell Biol. 120:1123-1135. Wong, S. H., and W. Hong. 1993. J. Biol. Chem. 268:22853-22862). The different properties of these two signals suggest that each operates by a different mechanism.     

Mechanism of inhibition of HIV-1 integrase by G-tetrad-forming oligonucleotides in Vitro

The G-tetrad-forming oligonucleotides and have been identified as potent inhibitors of human immunodeficiency virus type 1 integrase (HIV-1 IN) activity (Rando, R. F., Ojwang, J., Elbaggari, A., Reyes, G. R., Tinder, R., McGrath, M. S., and Hogan, M. E. (1995) J. Biol. Chem. 270, 1754-1760; Mazumder, A., Neamati, N., Ojwang, J. O., Sunder, S., Rando, R. F., and Pommier, Y. (1996) Biochemistry 35, 13762-13771; Jing, N., and Hogan, M. E. (1998) J. Biol. Chem. 273, 34992-34999). To understand the inhibition of HIV-1 IN activity by the G-quartet inhibitors, we have designed the oligonucleotides and, composed of three and four G-quartets with stem lengths of 19 and 24 A, respectively. The fact that increasing the G-quartet stem length from 15 to 24 A kept inhibition of HIV-1 IN activity unchanged suggests that the binding interaction occurs between a GTGT loop domain of the G-quartet inhibitors and a catalytic site of HIV-1 IN, referred to as a face-to-face interaction. Docking the NMR structure of (Jing and Hogan (1998)) into the x-ray structure of the core domain of HIV-1 IN, HIV-1 IN-(51-209) (Maignan, S., Guilloteau, J.-P. , Qing, Z.-L., Clement-Mella, C., and Mikol, V. (1998) J. Mol. Biol. 282, 359-368), was performed using the GRAMM program. The statistical distributions of hydrogen bonding between HIV-1 IN and were obtained from the analyses of 1000 random docking structures. The docking results show a high probability of interaction between the GTGT loop residues of the G-quartet inhibitors and the catalytic site of HIV-1 IN, in agreement with the experimental observation.     

Comparison of analyses of DNA curvature

Popular programs for characterizing DNA structure include Curves 5.1 (Lavery, R. and Sklenar, H., J. Biomol. Struct. Dyn. 6, 63-91, 1988; Lavery, R. and Sklenar, H., J. Biomol. Struct. Dyn. 6, 655-67, 1989) and Freehelix98 (Dickerson, R. E., Nucleic Acids Res. 26, 1906-1926, 1998), along with the more recent 3DNA (X. J. Lu, Z. Shakked and W. K. Olson., J. Mol. Biol. 300, 819-840 (2000). Given input of structural coordinates, all of these programs return values of the local helical parameters, such as roll, tilt, twist, etc. The first two programs also provide characterization of global curvature. Madbend (Strahs, D. and Schlick, T., J. Mol. Biol. 301, 643-663, 2000), a program that computes global curvature from local roll, tilt, and twist parameters, can be applied to the output of all three structural programs. We have compared the curvature predicted by the three programs with and without the use of Madbend. Global bend magnitudes and directions as well as values of helical kinks were calculated for four high-resolution DNA structures and four model DNA helices. Global curvature determined by Curves 5.1 without Madbend was found to differ from values obtained using Freehelix98 with or without Madbend or 3DNA and Curves 5.1 with Madbend. Using model helices, this difference was attributed the fact that Curves 5.1 is the only program sensitive to changes in axial displacement, such as shift and slide. Madbend produced robust values of bend magnitude and direction, and displayed little sensitivity to axis displacement or the source of local helical parameters. Madbend also appears to be the method of choice for bending comparisons of high-resolution structures with results from cyclization kinetics, a method that measures DNA curvature as a vectorial sum of local roll and tilt angles.     

Quaternary structure of rice nonsymbiotic hemoglobin

Plant nonsymbiotic hemoglobins are hexacoordinate heme proteins found in all plants. Although expression is linked with hypoxic environmental conditions (Taylor, E. R., Nie, X. Z., Alexander, W. M., and Hill, R. D. (1994) Plant Mol. Biol. 24, 853-862), no discrete physiological function has yet been attributed to this family of proteins. The crystal structure of a nonsymbiotic hemoglobin from rice has recently been determined. The crystalline protein is homodimeric and hexacoordinate with two histidine side chains coordinating the heme iron atom. Despite the fact that the amino acids responsible for the subunit interface are relatively conserved among the nonsymbiotic hemoglobins, previous work suggests that this group of proteins might display variability in quaternary structure (Duff, S. M. G., Wittenberg, J. B., and Hill, R. D. (1997) J. Biol. Chem. 272, 16746-16752; Arredondo-Peter, R., Hargrove, M. S., Sarath, G., Moran, J. F., Lohrman, J., Olson, J. S., and Klucas, R. V. (1997) Plant Physiol. 115, 1259-1266). Analytical ultracentrifugation and size exclusion high pressure liquid chromatography were used to investigate the quaternary structure of rice nonsymbiotic hemoglobin at various states of ligation and oxidation. Additionally, site-directed mutagenesis was used to test the role of several interface amino acids in dimer formation and ligand binding. Results were analyzed in light of possible physiological functions and indicate that the plant nonsymbiotic hemoglobins are not oxygen transport proteins but more closely resemble known oxygen sensors.     

Actin's propensity for dynamic filament patterning

Actin, through its various forms of assembly, provides the basic framework for cell motility, cell shape and intracellular organization in all eukaryotic cells. Many other cellular processes, for example endocytosis and cytokinesis, are also associated with dynamic changes of the actin cytoskeleton. Important prerequisites for actin's functional diversity are its intrinsic ability to rapidly assemble and disassemble filaments and its spatially and temporally well-controlled supramolecular organization. A large number of proteins that interact with actin, collectively referred to as actin-binding proteins (ABPs), carefully orchestrate different scenarios. Since its isolation in 1994 [Machesky, L.M. et al. (1994) J. Cell Biol. 127, 107-115], the Arp2/3 complex containing the actin-related proteins Arp2 and Arp3 has evolved to be one of the main players in the assembly and maintenance of many actin-based structures in the cell (for review see [Borths, E.L. and Welch, M.D. (2002) Structure 10, 131-135; May, R.C. (2001) Cell Mol. Life Sci. 58, 1607-1626; Pollard, T.D. et al. (2000) Rev. Biophys. Biomol. Struct. 29, 545-576; Welch, M.D. (1999) Trends Cell Biol. 11, 423-427]). In particular, when it comes to the assembly of the intricate branched actin network at the leading edge of lamellipodia, the Arp2/3 complex seems to have received all the attention in recent years. In parallel, but not so much in the spotlight, several reports showed that actin on its own can assume different conformations [Bubb, M.R. et al. (2002) J. Biol. Chem. 277, 20999-21006; Schoenenberger, C.-A. et al. (1999) Microsc. Res. Tech. 47, 38-50; Steinmetz, M.O. et al. (1998) J. Mol. Biol. 278, 793-811; Steinmetz, M.O. et al. (1997) J. Cell Biol. 138, 559-574; Millonig, R., Salvo, H. and Aebi, U. (1988) J. Cell Biol. 106, 785-796] through which it drives its supramolecular patterning, and which ultimately generate its functional diversity.     

Solvation effects on the sequence variability of DNA double helical conformations

The role of solvation on the sequence dependent conformational variabilities in DNA has been studied by calculating hydration free energies from solvent accessible surface areas for several base steps, as a function of various helical parameters, roll, twist and propeller twist. The results of roll calculations suggest opposite trends for AA and GG steps, with the former tending to have a compressed minor groove and the latter a compressed major groove. These trends are consistent with the experimental findings on sequence preferences and the nature of anisotropic bending of DNA observed in nucleosomes (Drew, H.R. and Travers, A.A., J. Mol. Biol. 186, 773-790 (1985); Satchwell, S.C., Drew, H.R. and Travers, A.A., J. Mol. Biol. 191, 659-675 (1986)) and CAP-DNA interactions (Gartenberg, M.R. and Crothers, D.M., Nature 333, 824-829, (1988)). Solvation energy profiles also indicate preferences for the base pairs in GG and AA steps to adopt low and high propeller twists, respectively. Such agreements may either reflect a coincidence of solvation effects with other energy terms or a dominance of solvent effects. The results are discussed in the context of the crystallographic observations of structural tendencies.     

NMR study of the exchange rates of allosterically responsive labile protons in the heme pockets of hemoglobin A.

1H NMR spectroscopy has been used to measure the proximal histidyl labile ring proton (NH) rates of exchange with bulk solvent in the individual subunits of hemoglobin (Hb) A. These protons displayed a substantial decrease in their exchange rates in comparison with related monomeric proteins and exhibited sensitivity to the quarternary state. With the beta subunit NH, the exchange behaviour was similar to an allosterically responsive subset of protons, which have been identified using 1H-3H methods (Englander, J.J., R. Rogero, and S. W. Englander, 1983, J. Mol. Biol. 169:325-344). Assuming similar exchange mechanisms for the two subunits, the NMR data suggested a more flexible alpha than beta subunit in Hb A.