Bending and curvature calculations in B-DNA.

A simple program, BEND, has been written to calculate the magnitude of local bending and macroscopic curvature at each point along an arbitrary B-DNA sequence, using any desired bending model that specifies values of twist, roll and tilt as a function of sequence. The program has been used to evaluate six different DNA bending models in three categories. Two are bent non-A-tract models: (a) A new model based on the nucleosome positioning data of Satchwell et al 1986 (J. Mol. Biol. 191, 659-675), (b) The model of Calladine et al 1988 (J. Mol. Biol. 201, 127-137). Three are bent A-tract models: (c) The wedge model of Bolshoy et al 1991 (Proc. Natl. Acad. Sci. USA 88, 2312-2316), (d) The model of Cacchione et al 1989 (Biochem. 28, 8706-8713), (e) A reversed version of model (b). The last is a junction model: (f) The model of Koo & Crothers 1988 (Proc. Natl. Acad. Sci. USA 85, 1763-1767). Although they have widely different assumptions and values for twist, roll and tilt, all six models correctly predict experimental A-tract curvature as measured by gel retardation and cyclization kinetics, but only the new nucleosome positioning model is successful in predicting curvature in regions containing phased GGGCCC sequences. This model--showing local bending at mixed sequence DNA, strong bends at the sequence GGC, and straight, rigid A-tracts--is the only model consistent with both solution data from gel retardation and cyclization kinetics and structural data from x-ray crystallography.     

A theoretical model for the prediction of sequence-dependent nucleosome thermodynamic stability

A theoretical model for predicting nucleosome thermodynamic stability in terms of DNA sequence is advanced. The model is based on a statistical mechanical approach, which allows the calculation of the canonical ensemble free energy involved in the competitive nucleosome reconstitution. It is based on the hypothesis that nucleosome stability mainly depends on the bending and twisting elastic energy to transform the DNA intrinsic superstructure into the nucleosomal structure. The ensemble average free energy is calculated starting from the intrinsic curvature, obtained by integrating the dinucleotide step deviations from the canonical B-DNA and expressed in terms of a Fourier series, in the framework of first-order elasticity. The sequence-dependent DNA flexibility is evaluated from the differential double helix thermodynamic stability. A large number of free-energy experimental data, obtained in different laboratories by competitive nucleosome reconstitution assays, are successfully compared to the theoretical results. They support the hypothesis that the stacking energies are the major factor in DNA rigidity and could be a measure of DNA stiffness. A dual role of DNA intrinsic curvature and flexibility emerges in the determination of nucleosome stability. The difference between the experimental and theoretical (elastic) nucleosome-reconstitution free energy for the whole pool of investigated DNAs suggests a significant role for the curvature-dependent DNA hydration and counterion interactions, which appear to destabilize nucleosomes in highly curved DNAs. This model represents an attempt to clarify the main features of the nucleosome thermodynamic stability in terms of physical-chemical parameters and suggests that in molecular systems with a large degree of complexity, the average molecular properties dominate over the local features, as in a statistical ensemble.     

RFLP analysis of 1-aminocyclopropane-1-carboxylate synthase ACC2 and ACC4 genes from Polish cultivars of tomato

An important trait of tomato is the rate of fruit ripening, strongly dependent on ethylene production. The ripening-related ethylene synthesis in tomato is controlled mainly by 1-aminocyclopropane-1-carboxylate synthase LE-ACS2 and LE-ACS4 isoenzymes (Rottmann et al., 1991, J. Mol. Biol. 222: 937; Lincoln et al., 1993, J. Biol. Chem. 268: 19422; Barry et al., 2000, Plant Physiol. 123: 979). In spite of numerous reports on the LE-ACS2 and LE-ACS4 gene expression, only ones considered the genomic organisation each of these genes (Rottmann et al., 1991; Lincoln et al., 1993) reported one copy of each of these genes in tomato cv VF36. In this article we suggest that the genomic organisation of LE-ACS2 and LE-ACSS4 genes may depend on tomato cultivars and may differ from that described by the above authors. The results of Southern analyses of genomic DNAs from 17-day old seedlings (cultivars Jaga, Halicz, Betalux, New Yorker) imply that the genomic organisation of LE-ACS2 and LE-ACS4 genes in Polish cultivars differs from that reported for cv VF36.     

Evidence for the extent of insertion of the active site loop of intact alpha 1 proteinase inhibitor in beta-sheet A.

The extent of insertion of beta-strand s4A into sheet A in intact serpin alpha 1-proteinase inhibitor (alpha 1PI has been probed by peptide annealing experiments [Schulze et al. (1990) Eur. J. Biochem. 194, 51-56]. Twelve synthetic peptides of systematically varied length corresponding in sequence to the unprimed (N-terminal) side of the active site loop were complexed with alpha 1PI. The complexes were then characterized by circular dichroism spectroscopy and tested for inhibitory activity. Four peptides formed complexes which retained inhibitory activity, one of which was nearly as effective as the native protein. Comparison with the three dimensional structures of cleaved alpha 1PI [L?bermann et al. (1984) J. Mol. Biol. 177, 531-556] and plakalbumin [Wright et al. (1990) J. Mol. Biol. 213, 513-528] supports a model in which alpha 1PI requires the insertion of a single residue, Thr345, into sheet A for activity.     

Dual role of DNA intrinsic curvature and flexibility in determining nucleosome stability

A statistical mechanistic approach to evaluate the sequence-dependent thermodynamic stability of nucleosomes is proposed. The model is based on the calculation of the DNA intrinsic curvature, obtained by integrating the nucleotide step deviations from the canonical B-DNA structure, and on the evaluation of the first order elastic distortion energy to reach the nucleosomal superstructure. Literature data on the free energy of nucleosome formation as obtained by competitive nucleosome reconstitution of a significant pool of different DNA sequences were compared with the theoretical results, and a satisfactorily good correlation was found. A striking result of the comparison is the emergence of two opposite roles of the DNA intrinsic curvature and flexibility in determining nucleosome stability. Finally, the obtained results suggest that the curvature-dependent DNA hydration should play a relevant role in the sequence-dependent nucleosome stability.     

Nuclear Overhauser data and stereochemical considerations suggest that netropsin binds symmetrically within the minor groove of poly(dA).poly(dT), forming hydrogen bonds with both strands of the double helix

Sarma et al. (J. Biomol. Str. and Dynam. 2, 1085 (1985) have proposed, on the basis of nuclear magnetic resonance experiments on the complex of netropsin with poly(dA).poly(dT), that the drug molecule lies asymmetrically along the dA side of the minor groove and makes hydrogen bonds only with the dA strand. If the crystal structure analyses of B-DNA (Fratini et al., J. Biol. Chem. 257, 14686 (1982] and of its complex with netropsin (Kopka et al., J. Mol. Biol. 183, 553 (1985] are any guide, this off-center, wide-groove model is stereochemically unlikely. More to the point, the off-center model is unnecessary to explain the observed nmr data. All of the nuclear Overhauser and other observations are fully explained by the structure seen in the x-ray crystal analysis, in which netropsin sits squarely centered within the minor groove, making bifurcated hydrogen bonds with both strands.     

Unrestrained stochastic dynamics simulations of the UUCG tetraloop using an implicit solvation model.

Three unrestrained stochastic dynamics simulations have been carried out on the RNA hairpin GGAC[UUCG] GUCC, using the AMBER94 force field (Cornell et al., 1995. J. Am. Chem. Soc. 117:5179-5197) in MacroModel 5.5 (Mohamadi et al., 1990. J. Comp. Chem. 11:440-467) and either the GB/SA continuum solvation model (Still et al., 1990. J. Am. Chem. Soc. 112:6127-6129) or a linear distance-dependent dielectric (1/R) treatment. The linear distance-dependent treatment results in severe distortion of the nucleic acid structure, restriction of all hydroxyl dihedrals, and collapse of the counterion atmosphere over the course of a 5-ns simulation. An additional vacuum simulation without counterions shows somewhat improved behavior. In contrast, the two GB/SA simulations (1.149 and 3.060 ns in length) give average structures within 1.2 A of the initial NMR structure and in excellent agreement with results of an earlier explicit solvent simulation (Miller and Kollman, 1997. J. Mol. Biol. 270:436-450). In a 3-ns GB/SA simulation starting with the incorrect UUCG tetraloop structure (Cheong et al., 1990. Nature. 346:680-682), this loop conformation converts to the correct loop geometry (Allain and Varani, 1995. J. Mol. Biol. 250:333-353), suggesting enhanced sampling relative to the previous explicit solvent simulation. Thermodynamic effects of 2'-deoxyribose substitutions of loop nucleotides were experimentally determined and are found to correlate with the fraction of time the ribose 2'-OH is hydrogen bonded and the distribution of the hydroxyl dihedral is observed in the GB/SA simulations. The GB/SA simulations thus appear to faithfully represent structural features of the RNA without the computational expense of explicit solvent.     

A model for human cardiac troponin C and for modulation of its Ca2+ affinity by drugs

Calcium sensitizers are drugs which increase force development in striated muscle by sensitizing myofilaments to Ca2+. This can happen by increasing Ca2+ affinity of the regulatory domain of Ca2+ binding protein troponin C. High resolution crystal structures of two calcium binding proteins, calmodulin (Babu et al.: J. Mol. Biol. 203:191-204, 1988) and skeletal troponin C (Satyshur et al.: J. Biol. Chem. 263:1628-1647, 1988; Herzber et al.: J. Mol. Biol. 203:761-779, 1988), have recently been published. This makes it possible to model in detail the calcium-sensitizing action of drugs on troponin C. In this study a model of human cardiac troponin C in three-calcium state has been constructed. When calcium is bound to calcium site II of cardiac troponin C an open conformation of the protein results, which has a hydrophobic pocket surrounded by a few polar side chains. Complexation of three drugs, trifluoperazine, bepridil, and pimobendan, to the hydrophobic pocket is studied using energy minimization techniques. Two different binding modes are found, which differ in the location of a strong electrostatic interaction. In analogy with the crystal structure of skeletal troponin C it is hypothezed that in cardiac troponin C an interaction occurs between Gln-50 and Asp-88, which has a long-range effect on calcium binding. The binding modes of drugs, where a strong interaction with Asp-88 exists, can effectively prevent the interaction between Asp-88 and Gln-50 in the protein, and are proposed to be responsible for the calcium-sensitizing properties of the studied drugs.     

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.     

Use of pair potentials across protein interfaces in screening predicted docked complexes.

Empirical residue-residue pair potentials are used to screen possible complexes for protein-protein dockings. A correct docking is defined as a complex with not more than 2.5 A root-mean-square distance from the known experimental structure. The complexes were generated by "ftdock" (Gabb et al. J Mol Biol 1997;272:106-120) that ranks using shape complementarity. The complexes studied were 5 enzyme-inhibitors and 2 antibody-antigens, starting from the unbound crystallographic coordinates, with a further 2 antibody-antigens where the antibody was from the bound crystallographic complex. The pair potential functions tested were derived both from observed intramolecular pairings in a database of nonhomologous protein domains, and from observed intermolecular pairings across the interfaces in sets of nonhomologous heterodimers and homodimers. Out of various alternate strategies, we found the optimal method used a mole-fraction calculated random model from the intramolecular pairings. For all the systems, a correct docking was placed within the top 12% of the pair potential score ranked complexes. A combined strategy was developed that incorporated "multidock," a side-chain refinement algorithm (Jackson et al. J Mol Biol 1998;276:265-285). This placed a correct docking within the top 5 complexes for enzyme-inhibitor systems, and within the top 40 complexes for antibody-antigen systems.     

Differential scanning calorimetric studies of E. coli aspartate transcarbamylase. III. The denaturational thermodynamics of the holoenzyme with single-site mutations in the catalytic chain

Aspartate transcarbamylase (EC 2.1.3.2) from E. coli is a multimeric enzyme consisting of two catalytic subunits and three regulatory subunits whose activity is regulated by subunit interactions. Differential scanning calorimetric (DSC) scans of the wild-type enzyme consist of two peaks, each comprised of at least two components, corresponding to denaturation of the catalytic and regulatory subunits within the intact holoenzyme (Vickers et al., J. Biol. Chem. 253 (1978) 8493; Edge et al., Biochemistry 27 (1988) 8081). We have examined the effects of nine single-site mutations in the catalytic chains. Three of the mutations (Asp-100-Gly, Glu-86-Gln, and Arg-269-Gly) are at sites at the C1: C2 interface between c chains within the catalytic subunit. These mutations disrupt salt linkages present in both the T and R states of the molecule (Honzatko et al., J. Mol. Biol. 160 (1982) 219; Krause et al., J. Mol. Biol. 193 (1987) 527). The remainder (Lys-164-Ile, Tyr-165-Phe, Glu-239-Gln, Glu-239-Ala, Tyr-240-Phe and Asp-271-Ser) are at the C1: C4 interface between catalytic subunits and are involved in interactions which stabilize either the T or R state. DSC scans of all of the mutants except Asp-100-Gly and Arg-269-Gly consisted of two peaks. At intermediate concentrations, Asp-100-Gly and Arg-269-Gly had only a single peak near the Tm of the regulatory subunit transition in the holoenzyme, although their denaturational profiles were more complex at high and low protein concentrations. The catalytic subunits of Glu-86-Gln, Lys-164-Ile and Asp-271-Ser appear to be significantly destabilized relative to wild-type protein while Tyr-165-Phe and Tyr-240-Phe appear to be stabilized. Values of delta delta G degree cr, the difference between the subunit interaction energy of wild-type and mutant proteins, evaluated as suggested by Brandts et al. (Biochemistry 28 (1989) 8588) range from -3.7 kcal mol-1 for Glu-86-Gln to 2.4 kcal mol-1 for Tyr-165-Phe.     

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.     

A theoretical study of the stability of DNA binding with cis/trans platin

Both cis- and trans-platins are known to form intra- and interstrand cross-linking with DNA. Since the nature and strength of binding is different, it makes their efficacy as anti-tumour drug different. In the present communication, we report theoretical analysis by using an amended Zimm and Bragg theory, to explain the melting behaviour and heat capacity of DNA with and without platin binding. The sharpness of transition has been examined in terms of half width and sensitivity parameter (deltaH/sigma). The experimental measurements of Pilch et al (J Mol Biol 2000, 296, 803) and Ctirad and Brabec (J Biol Chem 2001, 276, 9655) have been used.     

Unique features of the structural model of 'hard' cuticle proteins: implications for chitin-protein interactions and cross-linking in cuticle

Cuticular proteins are one of the determinants of the physical properties of cuticle. A common consensus region (extended R&R Consensus) in these proteins binds to chitin, the other major component of cuticle. We previously predicted the preponderance of beta-pleated sheet in the consensus region and proposed its responsibility for the formation of helicoidal cuticle (Iconomidou et al., Insect Biochem. Mol. Biol. 29 (1999) 285). Subsequently, we verified experimentally the abundance of antiparallel beta-pleated sheet in the structure of cuticle proteins (Iconomidou et al., Insect Biochem. Mol. Biol. 31 (2001) 877). Homology modelling of soft (RR-1) cuticular proteins using bovine plasma retinol binding protein (RBP) as a template revealed an antiparallel beta-sheet half-barrel structure as the basic folding motif (Hamodrakas et al., Insect Biochem. Molec. Biol. 32 (2002) 1577). The RR-2 proteins characteristic of hard cuticle, have a far more conserved consensus and frequently more histidine residues. Extension of modelling to this class of consensus, in this work, reveals in detail several unique features of the proposed structural model to serve as a chitin binding structural motif, thus providing the basis for elucidating cuticle's overall architecture and chitin-protein interactions in cuticle.     

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).     

Nuclear Overhauser Data and Stereochemical Considerations Suggest that Netropsin Binds Symmetrically Within the Minor Groove of Poly(dA)·Poly(dT), Forming Hydrogen Bonds with Both Strands of the Double Helix

Abstract

Sarma et al. (J. Biomol. Str. and Dynam. 2, 1085 (1985) have proposed, on the basis of nuclear magnetic resonance experiments on the complex of netropsin with poly(dA)·poly(dT), that the drug molecule lies asymmetrically along the dA side of the minor groove and makes hydrogen bonds only with the dA strand. If the crystal structure analyses of B-DNA (Fratini et al., J. Biol. Chem. 257, 14686 (1982)) and of its complex with netropsin (Kopka et al., J. Mol. Biol. 183, 553 (1985)) are any guide, this off-center, wide-groove model is stereochemically unlikely. More to the point, the off-center model is unnecessary to explain the observed nmr data. All of the nuclear Overhauser and other observations are fully explained by the structure seen in the x-ray crystal analysis, in which netropsin sits squarely centered within the minor groove, making bifurcated hydrogen bonds with both strands.     

Comparative study of satellite sequences and phylogeny of five species from the genus Palorus (Insecta, Coleoptera).

Major satellite sequences are analysed in the three tenebrionid beetles Palorus cerylonoides, P. genalis, and P. ficicola, and compared with the ones from P. ratzeburgii and P. subdepressus reported elsewhere. All of them are A+T rich, pericentromerically located, and with lengths of about 150 bp, either in the form of monomers or formed by more complex repeating units. A preliminary phylogenetic analysis of Palorus species using the 3' end of the mitochondrial Cytochrome Oxidase I gene shows that the five Palorus species have been diverging for a considerable amount of evolutionary time, with the pair P. ratzeburgii and P. genalis being the most closely related. Only these two taxa showed some similarity between their respective high-copy-number satellite sequences, while other satellites are mutually unrelated and might have originated independently. However, all the satellites have in common tertiary structure induced by intrinsic DNA curvature, a characteristic which is conserved within the genus. Palorus major satellites were previously detected in the genomes of congeneric species as low-copy-number clusters (Mestrovi? et al., Mol. Biol. Evol. 15: 1062-1068. 1998). Given the divergences between the analysed species, the substitution rate deduced from high- and low-copy-number repeats is unexpectedly low. The presence of sequence-induced DNA curvature in all Palorus satellites and similar satellite DNAs in the species pair P. ratzeburgii and P. genalis suggest (i) that constraints are at the tertiary structure; and (ii) that the satellite DNA evolutionary turnover can be dependent on the history of the taxa under study, resulting in retention of similar satellites in related taxa.