A theoretical study of the sequence specificity in binding of lexitropsins to B-DNA

A theoretical study is presented on the binding to B-DNA of a series of lexitropsins, these ligands being netropsin derivatives in which one or both of the pyrrole rings have been replaced by imidazoles. The best complexes have been located by energy minimisation taking into account nucleic acid flexibility, ligand flexibility, explicit, mobile counterions and solvent dielectric effects. Calculations have been performed for two homopolymeric DNA receptor sequences, AT base sequence, which only decreases in the imidazole derivatives. These results emphasize the decisive role of the molecular electrostatic potential of the nucleic acid in determining the sequence selectivity of these ligands, as opposed to the postulated role of adenine C2 - pyrrole beta hydrogen contacts.     

A Theoretical Study of the Sequence Specificity in Binding of Lexitropsins to B-DNA

Abstract

A theoretical study is presented on the binding to B-DNA of a series of lexitropsins, these ligands being netropsin derivatives in which one or both of the pyrrole rings have been replaced by imidazoles. The best complexes have been located by energy minimisation taking into account nucleic acid flexibility, ligand flexibility, explicit, mobile counterions and solvent dielectric effects. Calculations have been performed for two homopolymeric DNA receptor sequences, AT and GC. All the compounds studied exhibit an overall binding preference for the AT base sequence, which only decreases in the imidazole derivatives. These results emphasize the decisive role of the molecular electrostatic potential of the nucleic acid in determining the sequence selectivity of these ligands, as opposed to the postulated role of adenine C2 - pyrrole β hydrogen contacts.     

Anthramycin binding to deoxyribonucleic acid-mitomycin C complexes. Evidence for drug-induced deoxyribonucleic acid conformational change and cooperativity in mitomycin C binding

Anthramycin and mitomycin C (MC) are two DNA reactive drugs, which bind covalently to GC pairs producing different effects on DNA: anthramycin stiffening and MC distorsion. This paper describes experiments in which we have used anthramycin as a probe to sense quantitatively the effects on DNA of MC binding. Saturation binding experiments show that both anthramycin and MC partially inhibit the binding of the other drug to DNA (maximum inhibition by MC and anthramycin, 22.4% and 19.7%, respectively) but by a mechanism other than direct site exclusion. This suggests that MC binds in the major groove of DNA, since anthramycin is known to bind in the minor groove. An abrupt reduction in the binding of anthramycin to DNA-MC complexes occurs between MC binding ratios of 0.030 and 0.035, which parallels and probably results from sudden intensification of a MC-induced DNA conformational change occurring between these binding ratios. Dialysis measurements indicate that anthramycin is very possibly binding at sites distant from MC sites and suggest a clustering of closely bound MC chromophores resulting from possible cooperative binding. S1 nuclease digest experiments demonstrate an initial enhancement of nuclease activity in DNA-MC complexes, the magnitude of which correlates well with the reduction of anthramycin binding, relative to the degree of MC binding. The enhanced nuclease activity in these complexes indicates regions of exposed DNA or helix base distortion which is related to or is the result of conformational change.     

Shuffling of amino acid sequence: an important control in synthetic peptide studies of nucleic acid-binding domains. Binding properties of fragments of a conserved eukaryotic RNA binding motif.

We have used synthetic peptides to study a conserved RNA binding motif in yeast poly(A)-binding protein. Two peptides, 45 and 44 amino acids in length, corresponding to amino and carboxyl halves of a 90-amino acid RNA-binding domain in the protein were synthesized. While the amino-terminal peptide had no significant affinity for nucleic acids, the carboxyl-terminal peptide-bound nucleic acids with similar characteristics to that for the entire 577 residue yeast poly(A)-binding protein. In 100 mM NaCl, the latter peptide retained over 50% of the intrinsic binding free energy of the protein, as well as, similar RNA versus DNA binding specificity. However, shuffling of the sequence of this 44 residue peptide had surprisingly little effect on its nucleic acid binding properties suggesting the overriding importance of amino acid composition as opposed to primary sequence. Deletion studies on the 44 residue peptide with the "correct" sequence succeeded in identifying amino acids important for conferring RNA specificity and for increasing our understanding of the molecular basis for nucleic acid binding by synthetic peptides. The shuffled peptide study, however, clearly indicates that considerable caution must be exercised before extrapolating results of structure/function studies on synthetic peptide analogues to the parent protein.     

Comparison of properties of the in vitro and cellular anthramycin-DNA adducts and characterization of the reaction of anthramycin with chromatin

The reaction of anthramycin with DNA has been examined to determine the chemical identity of the adduct which forms in a living cell and to observe the effects of the nucleosome structure of chromatin on drug binding. The chemical identity of the cellular adduct was probed by comparing various properties of the cellular adduct to properties of the known, in vitro adduct. The effect of the histones on anthramycin binding was investigated by time-course binding reactions. Results indicate that the properties of the cellular anthramycin-DNA adduct are similar to the in vitro adduct. The histone proteins associated with DNA in chromatin were found to decrease both the reaction kinetics and the final levels of anthramycin binding. Anthramycin reacts appreciably with nucleosome core DNA, but appears to exhibit a preference for linker DNA.     

Proteins that shape DNA

During the last five years, considerable accumulation of data on nucleic acids metabolism leads to the discovery of a number of proteins designed to change the conformation of DNA and to "shape" it. Experimental results emphasize the importance of the conformation and the flexibility of DNA itself in such interactions. The mutual recognition of nucleic acids by proteins may be or not dependent on the nnucleotide sequence and in most cases is accompanied by conformational changes in the proteins involved. Among these are proteins that bind in stoichiometric amounts to DNA, proteins that promote the separation of the two strands in a duplex, and finally proteins that change the topology of DNA.     

Physical studies of tyrosine and tryptophan residues in mammalian A1 heterogeneous nuclear ribonucleoprotein. Support for a segmented structure.

The mammalian heterogeneous ribonucleoprotein (hnRNP) A1 and its constituent N-terminal domain, termed UP1, have been studied by steady-state and dynamic fluorimetry, as well as phosphorescence and optically detected magnetic resonance (ODMR) spectroscopy at cryogenic temperatures. The results of these diverse techniques coincide in assigning the site of the single tryptophan residue of A1, located in the UP1 domain, to a partially solvent-exposed site distal to the protein's nucleic acid binding surface. In contrast, tyrosine fluorescence is significantly perturbed when either protein associates with single-stranded polynucleotides. Tyr to Trp energy transfer at the singlet level is found for both UP1 and A1 proteins. Single-stranded polynucleotide binding induces a quenching of their intrinsic fluorescence emission, which can be attributed to a significant reduction (greater than 50%) of the Tyr contribution, while Trp emission is only quenched by approximately 15%. Tyrosine quenching effects of similar magnitude are seen upon polynucleotide binding by either UP1 (1 Trp, 4 Tyr) or A1 (1 Trp, 12 Tyr), strongly suggesting that Tyr residues in both the N-terminal and C-terminal domain of A1 are involved in the binding process. Tyr phosphorescence emission was strongly quenched in the complexes of UP1 with various polynucleotides, and was attributed to triplet state energy transfer to nucleic acid bases located in the close vicinity of the fluorophore. These results are consistent with stacking of the tyrosine residues with the nucleic acid bases. While the UP1 Tyr phosphorescence lifetime is drastically shortened in the polynucleotide complex, no change of phosphorescence emission maximum, phosphorescence decay lifetime or ODMR transition frequencies were observed for the single Trp residue. The results of dynamic anisotropy measurements of the Trp fluorescence have been interpreted as indicative of significant internal flexibility in both UP1 and A1, suggesting a flexible linkage connecting the two sub-domains in UP1. Theoretical calculations based on amino acid sequence for chain flexibility and other secondary structural parameters are consistent with this observation, and suggest that flexible linkages between sub-domains may exist in other RNA binding proteins. While the dynamic anisotropy data are consistent with simultaneous binding of both the C-terminal and the N-terminal domains to the nucleic acid lattice, no evidence for simultaneous binding of both UP1 sub-domains was found.     

Theoretical study of the sequence specificity in the covalent binding of the antitumor drug CC-1065 to DNA

A theoretical modelling is presented of the covalent adducts of the antitumor agent CC-1065 with B-DNA. The optimal complexes are obtained by energy minimisation, taking into account full structure flexibility, including the flexible rings of the ligand and DNA. The binding preference of CC-1065 with respect to base sequence is studied. The results obtained elucidate the origin of the preference for two AT base pairs on the 5'side of the modified adenine. The modifications of the DNA structure upon ligand covalent binding are discussed.     

Recognition of the DNA helix stabilizing anthramycin-N2 guanine adduct by UVRABC nuclease

The binding of the anti-tumor antibiotic anthramycin to a defined linear DNA fragment was investigated using both exonuclease III and lambda exonuclease. We show that most of the guanine residues are reactive toward anthramycin; however, several guanine residues showed preferential reactivity for the drug. Using purified UVRA, UVRB and UVRC proteins we present evidence that these three proteins in concert are able to recognize and produce specific strand cleavage flanking anthramycin-DNA adducts. The cleavage of anthramycin adducts by UVRABC nuclease is specific and results in strand breaks at five or six bases 5' and three or four bases 3'-flanking an adduct. At some guanine residues single incisions were observed only on one side of the adduct. The 5' strand breaks observed often occurred as doublet bands on sequencing gels, indicating plasticity in the site of 5' cleavage whereas the 3' cleavage did not show this effect. When DNA fragments modified with elevated levels of anthramycin were used as substrates the activity of the UVRABC nuclease toward the anthramycin adducts decreased. Possible mechanisms for the recognition and specific cleavage of the helix-stabilizing anthramycin DNA adduct and other helix destabilizing lesions by the UVRABC nuclease are discussed.     

Pyrrolo(1,4)benzodiazepine antitumor antibiotics. In vitro interaction of anthramycin, sibiromycin and tomaymycin with DNA using specifically radiolabelled molecules.

Anthramycin, tomaymycin and sibiromycin are pyrrolo(1,4)benzodiazepine antitumor antibiotics. These compounds react with DNA and other guanine-containing polydeoxynucleotides to form covalently bound antibiotic - polydeoxynucleotide complexes. Experiments utilizing radiolabelled antibiotics have led to the following conclusions: 1. Sibiromycin reacts much faster than either anthramycin or tomaymycin with DNA. 2. At saturation binding the final antibiotic to base ratios for sibiromycin, anthramycin and tomaymycin are 1 : 8.8,1: 12.9, and 1 : 18.2, respectively. 3. No reaction with RNA or protein occurs with the pyrrolo(1,4)benzodiazepine antibiotics. 4. Sibiromycin effectively competes for the same DNA binding sites as anthramycin and tomaymycin; however, there is only partial overlap for the same binding sites between anthramycin and tomaymycin. 5. Whereas all three pyrrolo(1,4)benzodiazepine antibiotic-DNA complexes are relatively stable to alkaline conditions, their stability under acidic conditions increases in the order tomaymycin, anthramycin and sibiromycin. 6. No loss of non-exchangeable hydrogens in either the pyrrol ring or the side chains of these antibiotics occurs upon formation of their complexes with DNA. 7. Unchanged antibiotic has been demonstrated to be released upon acid treatment of the anthramycin-DNA and tomaymycin-DNA complexes. 8. A Schiff base linkage between the antibiotics and DNA has been eliminated. The comparative reactivity of the three antibiotics towards DNA and the stability of their DNA complexes is discussed in relation to their structures. A working hypothesis for the formation of the antibiotic-DNA covalent complexes is proposed based upon the available information.     

Peptide nucleic acids: deoxyribonucleic acid mimics with a peptide backbone

This tutorial briefly describes a new class of synthetic biopolymer, which is referred to as peptide nucleic acid (PNA). In PNA, individual nucleobases are linked to an achiral neutral peptide backbone. PNA exhibits the hybridization characteristic (e.g., Watson—Crick duplex formation) of DNA. The achiral peptide backbone provides similar interbase distances as natural DNA, and adequate flexibility to permit base pair interactions with complementary RNA or DNA strands. Several potential applications of PNA oligomers in biotechnology are suggested. These include the use of PNAs as a probe for specific recognition of a DNA or RNA sequence selective, purification of nucleic acids via designed high affinity binding to PNA, screening for DNA mutations, and as possible therapeutic agents.     

ATP binding modulates the nucleic acid affinity of hepatitis C virus helicase

The helicase of hepatitis C virus (HCV) unwinds nucleic acid using the energy of ATP hydrolysis. The ATPase cycle is believed to induce protein conformational changes to drive helicase translocation along the length of the nucleic acid. We have investigated the energetics of nucleic acid binding by HCV helicase to understand how the nucleotide ligation state of the helicase dictates the conformation of its nucleic acid binding site. Because most of the nucleotide ligation states of the helicase are transient due to rapid ATP hydrolysis, several compounds were analyzed to find an efficient unhydrolyzable ATP analog. We found that the beta-gamma methylene/amine analogs of ATP, ATPgammaS, or [AlF4]ADP were not effective in inhibiting the ATPase activity of HCV helicase. On the other hand, [BeF3]ADP was found to be a potent inhibitor of the ATPase activity, and it binds tightly to HCV helicase with a 1:1 stoichiometry. Equilibrium binding studies showed that HCV helicase binds single-stranded nucleic acid with a high affinity in the absence of ATP or in the presence of ADP. Upon binding to the ATP analog, a 100-fold reduction in affinity for ssDNA was observed. The reduction in affinity was also observed in duplex DNA with 3' single-stranded tail and in RNA but not in duplex DNA. The results of this study indicate that the nucleic acid binding site of HCV helicase is allosterically modulated by the ATPase reaction. The binding energy of ATP is used to bring HCV helicase out of a tightly bound state to facilitate translocation, whereas ATP hydrolysis and product release steps promote tight rebinding of the helicase to the nucleic acid. On the basis of these results we propose a Brownian motor model for unidirectional translocation of HCV helicase along the nucleic acid length.     

Minor groove deformability of DNA: a molecular dynamics free energy simulation study

The conformational deformability of nucleic acids can influence their function and recognition by proteins. A class of DNA binding proteins including the TATA box binding protein binds to the DNA minor groove, resulting in an opening of the minor groove and DNA bending toward the major groove. Explicit solvent molecular dynamics simulations in combination with the umbrella sampling approach have been performed to investigate the molecular mechanism of DNA minor groove deformations and the indirect energetic contribution to protein binding. As a reaction coordinate, the distance between backbone segments on opposite strands was used. The resulting deformed structures showed close agreement with experimental DNA structures in complex with minor groove-binding proteins. The calculated free energy of minor groove deformation was approximately 4-6 kcal mol(-1) in the case of a central TATATA sequence. A smaller equilibrium minor groove width and more restricted minor groove mobility was found for the central AAATTT and also a significantly ( approximately 2 times) larger free energy change for opening the minor groove. The helical parameter analysis of trajectories indicates that an easier partial unstacking of a central TA versus AT basepair step is a likely reason for the larger groove flexibility of the central TATATA case.     

Dissociation of minor groove binders from DNA: insights from metadynamics simulations

We have used metadynamics to investigate the mechanism of noncovalent dissociation from DNA by two representatives of alkylating and noncovalent minor groove (MG) binders. The compounds are anthramycin in its anhydrous form (IMI) and distamycin A (DST), which differ in mode of binding, size, flexibility and net charge. This choice enables to evaluate the influence of such factors on the mechanism of dissociation. Dissociation of IMI requires an activation free energy of approximately 12 kcal/mol and occurs via local widening of the MG and loss of contacts between the drug and one DNA strand, along with the insertion of waters in between. The detachment of DST occurs at a larger free energy cost, approximately 16.5 or approximately 18 kcal/mol depending on the binding mode. These values compare well with that of 16.6 kcal/mol extracted from stopped-flow experiments. In contrast to IMI, an intermediate is found in which the ligand is anchored to the DNA through its amidinium tail. From this conformation, binding and unbinding occur almost at the same rate. Comparison between DST and with kinetic models for the dissociation of Hoechst 33258 from DNA uncovers common characteristics across different classes of noncovalent MG ligands.     

Static and Dynamic Conformational Properties of AT sequences in B-DNA

Abstract

A theoretical study of the optimal conformations of nucleic acid oligomers containing tracts of AT base pairs is presented. The oligomers are studied in isolation and complexed with netropsin, a minor groove binding ligand. The flexibility of the oligomers and of their complexes is calculated by adiabatic mapping with respect to the total winding angle. The results of this study show that in uncomplexed oligomers the dinucleotide junctions AA, AT and TA have very different structural parameters and different responses to winding stress. The TA junction is clearly the most flexible and is the principal site for accommodating the imposed overwinding. Complexation by netropsin leads to two important effects: firstly, the three junctions adopt more uniform structures, the largest changes again being observed for TA, secondly, the differences in flexibility as a function of sequence are strongly attenuated.     

Static and dynamic conformational properties of AT sequences in B-DNA.

A theoretical study of the optimal conformations of nucleic acid oligomers containing tracts of AT base pairs is presented. The oligomers are studied in isolation and complexed with netropsin, a minor groove binding ligand. The flexibility of the oligomers and of their complexes is calculated by adiabatic mapping with respect to the total winding angle. The results of this study show that in uncomplexed oligomers the dinucleotide junctions AA, AT and TA have very different structural parameters and different responses to winding stress. The TA junction is clearly the most flexible and is the principal site for accommodating the imposed overwinding. Complexation by netropsin leads to two important effects: firstly, the three junctions adopt more uniform structures, the largest changes again being observed for TA, secondly, the differences in flexibility as a function of sequence are strongly attenuated.     

DNA curvature and phosphate neutralization: an important aspect of specific protein binding

A theoretical study is presented of the influence of salt bridges between protein cationic side chains and DNA phosphates on DNA conformation and flexibility. Two DNA sequences are studied containing respectively the HNF3 and CAP binding sites. The effect of salt bridges is modelled by the neutralisation of net phosphate charges for the groups involved in such interactions in the complex. Energy optimised conformations are obtained by molecular mechanics calculations using the JUMNA program. Base sequence dependence is studied by moving the phosphate neutralisation pattern along the sequence, while normal mode analysis is used to evaluate DNA flexibility. The results show that phosphate neutralisation has a strong influence on DNA conformation. For the HNF3 binding sequence, the free oligomer is bent in direction very different from that observed in the protein complex. Phosphate neutralisation changes this direction by 45 degrees to within only 4 degrees of the direction in the complex, without changing the bending angle. For the CAP binding sequence, the free oligomer is already intrinsically curved in the direction favoured by the protein, but phosphate neutralisation increases the bending angle. For both oligomers studied these effects are strongly sequence dependent. It is also shown that oligomer flexibility cannot be explained by a simple superposition of the properties of successive dinucleotide steps. Important long range coupling effects are observed. However, for both sequence studied, phosphate neutralisation however leads to a reduction in oligomer flexibility.     

Solution structure of the single-stranded DNA binding protein of the filamentous Pseudomonas phage Pf3: similarity to other proteins binding to single-stranded nucleic acids.

The three-dimensional structure of the homodimeric single-stranded DNA binding protein encoded by the filamentous Pseudomonas bacteriophage Pf3 has been determined using heteronuclear multidimensional NMR techniques and restrained molecular dynamics. NMR experiments and structure calculations have been performed on a mutant protein (Phe36 --> His) that was successfully designed to reduce the tendency of the protein to aggregate. The protein monomer is composed of a five-stranded antiparallel beta-sheet from which two beta-hairpins and a large loop protrude. The structure is compared with the single-stranded DNA binding protein encoded by the filamentous Escherichia coli phage Ff, a protein with a similar biological function and DNA binding properties, yet quite different amino acid sequence, and with the major cold shock protein of Escherichia coli, a single-stranded DNA binding protein with an entirely different sequence, biological function and binding characteristics. The amino acid sequence of the latter is highly homologous to the nucleic acid binding domain (i.e. the cold shock domain) of proteins belonging to the Y-box family. Despite their differences in amino acid sequence and function, the folds of the three proteins are remarkably similar, suggesting that this is a preferred folding pattern shared by many single-stranded DNA binding proteins.