Molecular modelling study of changes induced by netropsin binding to nucleosome core particles.

It is well known that certain sequence-dependent modulators in structure appear to determine the rotational positioning of DNA on the nucleosome core particle. That preference is rather weak and could be modified by some ligands as netropsin, a minor-groove binding antibiotic. We have undertaken a molecular modelling approach to calculate the relative energy of interaction between a DNA molecule and the protein core particle. The histones particle is considered as a distribution of positive charges on the protein surface that interacts with the DNA molecule. The molecular electrostatic potentials for the DNA, simulated as a discontinuous cylinder, were calculated using the values for all the base pairs. Computing these parameters, we calculated the relative energy of interaction and the more stable rotational setting of DNA. The binding of four molecules of netropsin to this model showed that a new minimum of energy is obtained when the DNA turns toward the protein surface by about 180 degrees, so a new energetically favoured structure appears where netropsin binding sites are located facing toward the histones surface. The effect of netropsin could be explained in terms of an induced change in the phasing of DNA on the core particle. The induced rotation is considered to optimize non-bonded contacts between the netropsin molecules and the DNA backbone.     

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.     

Theoretical studies of the selective binding to DNA of two non-intercalating ligands: netropsin and SN 18071.

A theoretical study of the binding to DNA of netropsin and a bisquaternary ammonium heterocycle, SN 18071, is undertaken with an energy minimizing program based on empirical potential functions. The positioning of the ligand is achieved by force and torque calculations and its internal flexibility is taken into account. The binding preference of both drugs studied for the AT minor groove of B-DNA is shown to depend on both the electrostatic potential generated by the base sequence and the quality of the steric fit of the ligand in the groove. Ligand-DNA hydrogen bonds are shown to aid binding, but not to be essential in establishing binding preferences.     

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.     

Theoretical exploration of netropsin binding to tRNA(Phe)

Theoretical exploration of the possible interaction of netropsin with tRNAPhe indicates that binding should occur preferentially with the major groove of the T psi C stem of the macromolecule, specifically with the bases G51, U52, G53 and phosphates 52, 53, 61 and 62. This agrees with the recent crystallographic result of Rubin and Sundaralingam. It is demonstrated that the difference with respect to netropsin binding with B-DNA, where it occurs specifically in the minor groove of AT sequences, is due to the differences in the distribution of the electrostatic molecular potential generated by these different types of DNA: this potential is sequence dependent in B-DNA (located in the minor groove of AT sequences and the major groove of GC sequences), while it is sequence independent and always located in the major groove in A-RNA. The result demonstrates the major role of electrostatics in determining the location of the binding site.     

Effect of ionic strength and cationic DNA affinity binders on the DNA sequence selective alkylation of guanine N7-positions by nitrogen mustards

Large variations in alkylation intensities exist among guanines in a DNA sequence following treatment with chemotherapeutic alkylating agents such as nitrogen mustards, and the substituent attached to the reactive group can impose a distinct sequence preference for reaction. In order to understand further the structural and electrostatic factors which determine the sequence selectivity of alkylation reactions, the effect of increased ionic strength, the intercalator ethidium bromide, AT-specific minor groove binders distamycin A and netropsin, and the polyamine spermine on guanine N7-alkylation by L-phenylalanine mustard (L-Pam), uracil mustard (UM), and quinacrine mustard (QM) was investigated with a modification of the guanine-specific chemical cleavage technique for DNA sequencing. For L-Pam and UM, increased ionic strength and the cationic DNA affinity binders dose dependently inhibited the alkylation. QM alkylation was less inhibited by salt (100 mM NaCl), ethidium (10 microM), and spermine (10 microM). Distamycin A and netropsin (100 microM) gave an enhancement of overall QM alkylation. More interestingly, the pattern of guanine N7-alkylation was qualitatively altered by ethidium bromide, distamycin A, and netropsin. The result differed with both the nitrogen mustard (L-Pam less than UM less than QM) and the cationic agent used. The effect, which resulted in both enhancement and suppression of alkylation sites, was most striking in the case of netropsin and distamycin A, which differed from each other. DNA footprinting indicated that selective binding to AT sequences in the minor groove of DNA can have long-range effects on the alkylation pattern of DNA in the major groove.     

Theoretical Exploration of Netropsin Binding to tRNAPhe

Abstract

Theoretical exploration of the possible interaction of netropsin with tRNA^Phe indicates that binding should occur preferentially with the major groove of the TψC stem of the macromolecule, specifically with the bases G51, U52, G53 and phosphates 52, 53, 61 and 62. This agrees with the recent crystallographic result of Rubin and Sundaralingam. It is demonstrated that the difference with respect to netropsin binding with B-DNA, where it occurs specifically in the minor groove of AT sequences, is due to the differences in the distribution of the electrostatic molecular potential generated by these different types of DNA: this potential is sequence dependent in B-DNA (located in the minor groove of AT sequences and the major groove of GC sequences), while it is sequence independent and always located in the major groove in A-RNA. The result demonstrates the major role of electrostatics in determining the location of the binding site.     

Selective Binding of Synthetic Polypeptides to DNA of Varying Composition and Sequence: Effect of Minor Groove Binding Drugs

Abstract

Repetitive basic polypeptides containing lysine or arginine as every third amino acid were shown to cause DNA condensation at physiological salt concentration connected with selective DNA binding with respect to DNA composition and sequence. This selectivity is very similar to that existing in the case of histone H1 and other basic proteins and does not depend on polypeptide chain conformation. The effect of the minor groove binding drugs netropsin and distamycin was tested to elucidate the origin of the binding selectivity. The results suggest that the binding preferences are due to the variations in the conformation in various types of B-DNA that depend on DNA composition and sequence. The most important factor affecting the selectivity is probably the value of the negative electrostatic potential in the minor groove.     

Netropsin, distamycin A, bis-netropsins as selective inhibitors of restrictases and DNAse I

The simultaneous analysis of DNAase I "footprinting" data and restriction endonucleases inhibition data was performed on the same DNA end-labelled fragment. The inhibition induced by netropsin, a number of bis-netropsins and distamycin A was investigated. These experiments led us to the following conclusions. The restriction endonucleases inhibition by the ligands is caused by the ligand molecules binding in the close vicinity to the restriction endonuclease recognition sequence. The zone of +/- 4 bp from the center of the restriction endonuclease recognition sequence can be defined as the zone of the influence of the bounded ligand on the restriction endonuclease. But in this case the intersection of recognition sequence and the binding site occupied by a single ligand molecule is not sufficient for the inhibition to occur. Restriction endonuclease cutting sites protected by netropsin can be predicted basing upon known nucleotide sequence specificity of netropsin. Netropsin and bis-netropsins show different nucleotide sequence specificity. This fact can be used for selective inhibition of restriction endonucleases.     

Thermodynamics of interaction of a fluorescent DNA oligomer with the anti-tumour drug netropsin.

Fluorescence spectroscopy was used to study the interaction between the minor-groove-binding drug netropsin and the self-complementary oligonucleotide d(CTGAnPTTCAG)2 containing the fluorescent base analogue 2-aminopurine (nP). The binding of netropsin to this oligonucleotide causes strong quenching of the 2-aminopurine fluorescence, observed by steady-state as well as time-resolved spectroscopy. From fluorescence titrations, binding isotherms were recorded and evaluated. The parameters showed one netropsin binding site/oligonucleotide duplex and an association constant of about 10(5) M-1 at 25 degrees C, 3-4 orders of magnitude weaker than for an exclusive adenine/thymine host sequence. From the temperature dependence of the association constant the thermodynamic parameters were obtained as delta G = -29 kJ/mol, delta H = -12 kJ/mol and delta S = +55 J.mol-1.K-1 at 25 degrees C. These parameters resemble those of the interaction of poly[(dG-dC).(dG-dC)] with netropsin, indicating a mainly entropy-driven reaction. The amino group of 2-aminopurine, like that of guanine, resides in the minor groove of DNA. Therefore the relatively weak binding of netropsin to d(CTGAnPTTCAG)2 is probably related to partial blockage of the tight fit of netropsin into the preferred minor groove of an exclusive adenine/thymine host sequence.     

Molecular modelling of the interaction of carbocyclic analogues of netropsin and distamycin with d(CGCGAATTCGCG)2

A molecular mechanics and molecular dynamics approach was used to examine the structure of complexes formed between the d(CGCGAATTCGCG)2 duplex and netropsin, distamycin, and four carbocyclic analogues of netropsin and distamycin (1-4). The resulting structures of the ligand-DNA model complexes and their energetics were examined. It is predicted that the compounds 1-4 should have a decreased affinity for the minor groove of AT-rich regions in comparison to netropsin and distamycin. From the energetic analysis it appears that van der Waals and electrostatic interactions are more important than specific hydrogen bonds in stabilizing the ligand-duplex complexes. We predict that compounds 1 and 2 are effectively isohelical with the DNA minor groove. The superior DNA-binding afforded by 1 and 2 in comparison to 3 and 4 results from their more effective penetration into the minor groove and smaller perturbation of molecular structure upon complex formation.     

Identification and separation of components of calf thymus DNA using a CsCl-netropsin density gradient

Calf thymus DNA containing satellite components of various densities was used as a model to study the effect of netropsin on the density of DNA in a CsCl gradient. The binding of netropsin resulted in a decrease in density which depended upon the quantity of netropsin added and on the average composition of the DNA. Differences in density of DNA components were higher in CsCl - netropsin gradients than in simple CsCl gradients. By use of netropsin a main band and four satellite bands could be differentiated in calf thymus DNA. Satellite DNA's were isolated using preparative CsCl - netropsin gradient centrifugation and were characterised by density and homogeneity in native and in reassociated state. Two of the satellite components, with densities of 1.722 and 1.714 g/cm³, are probably of homogenous sequence, the other two components of densities 1.709 and 1.705 g/cm³ appear to be heterogeneous.     

Netropsin . dG-dG-dA-dA-dT-dT-dC-dC complex. Antibiotic binding at adenine . thymine base pairs in the minor groove of the self-complementary octanucleotide duplex.

The structure of the netropsin . dG-dG-dA-dA-dT-dT-dC-dC complex (one antibiotic molecule/self-complementary octanucleodide duplex) and its dynamics as a function of temperature have been monitored by the nuclear magnetic resonances of the Watson-Crick protons, the nonexchangeable base and sugar protons and the backbone phosphates. The antibiotic forms a complex with the nucleic acid duplex at the dA . dT-containing tetranucleotide segment dA-dA-dT-dT, with slow migration amongst potential binding sites at low temperature. The downfield shifts in the exchangeable protons of netropsin on complex formation demonstrate the contributions of hydrogen-bonding interactions between the antibiotic and the nucleic acid to the stability of the complex. Complex formation results in changes in the glycosidic torsion angles of both thymidine residues and one deoxyadenosine residue as monitored by chemical shift changes in the thymine C-6 and adenine C-8 protons. The close proximity of the pyrrole rings of the antibiotic and the base-pair edges in the minor groove is manifested in the downfield shifts (0.3--0.5 ppm) of the pyrrole C-3 protons of netropsin and one adenine C-2 proton and one thymine N-3 base-pair proton on complex formation. The internucleotide phosphates of the octanucleotide undergo 31P chemical shift changes on addition of netropsin and these may reflect, in part, contributions from electrostatic interactions between the charged ends of the antibiotic and the backbone phosphates of the nucleic acid.     

Variation of DNA sequence specificity of DNA-oligopeptide binding ligands related to netropsin: imidazole-containing lexitropsins

CD binding studies of nonintercalative oligopeptides related to netropsin, named lexitropsins, have been carried out with synthetic duplex DNAs and natural DNA. While netropsin possesses a high dA.dT sequence specificity, these ligands show a progressive lowering of the ability to bind to dA.dT basepairs in DNA and a dramatic reduction of the sequence specificity seen at high salt concentration due to a replacement of pyrrole moieties by imidazoles. This variation in DNA sequence specificity of lexitropsins is mirrored in corresponding large differences in the template inactivation of poly(dA-dT).poly(dA-dT) in the RNA polymerase reaction by these drugs. The presence of imidazole permits binding of the oligopeptide to dG.dC pairs, which is most effective for the triimidazole peptide. Results at increasing salt concentration reveal, however, that a tight binding to pure dG.dC sequences does not occur. A proper sequence containing dG.dC and dA.dT pairs is supposed to be required for a higher specificity. The CD data accord well with previously reported melting studies and are in favor of recent theoretical results suggesting that the diminished AT preference may be due to an increase in the complexation energy with the dG.dC pairs.     

Specific protection of DNA by distamycin A, netropsin and bis-netropsins against the action of DNAse I

Interaction of netropsin, distamycin A and a number of bis-netropsins with DNA fragments of definite nucleotide sequence was studied by footprinting technique. The nuclease protection experiments were made at fixed DNA concentration and varying ligand concentrations. The affinity of ligand for a DNA site was estimated from measurements of ligand concentration that causes 50% protection of the DNA site. Distribution pattern of the protected and unprotected regions along the DNA fragment was compared with the theoretically expected arrangement of the ligand along the same DNA. The comparison led us to the following conclusions: 1. Footprinting experiments show that at high levels of binding the arrangement of netropsin molecules along the DNA corresponds closely to the distribution pattern expected from theoretical calculations based on the known geometry of netropsin--DNA complex. However, the observed differences in the affinity of netropsin for various DNA sequences is markedly greater than that expected from theoretical calculations. 2. Netropsin exhibits a greater selectivity of binding than that expected for a ligand with three specific reaction centers associated with the antibiotic amide groups. It binds preferentially to DNA regions containing four or more successive AT pairs. Among 13 putative binding sites for netropsin with four or more successive AT pairs there are 11 strong binding sites and two weaker sites which are occupied at 2 D/P less than or equal to 1/9 and 2 D/P = 1/4, respectively. 3. The extent of specificity manifested by distamycin A is comparable to that shown by netropsin although the molecule of distamycin A contains four rather than three amide groups. At high levels of binding distamycin A occupies the same binding sites on DNA as netropsin does. 4. The binding specificity of bis-netropsins is greater than that of netropsin. Bis-netropsins can bind to DNA in such a way that the two netropsin-like fragments are implicated in specific interaction with DNA base pairs. However, the apparent affinity of bis-netropsins estimated from footprinting experiments is comparable with that of netropsin for the same DNA region. 5. At high levels of binding bis-netropsins and distamycin A (but not netropsin) can occupy any potential site on DNA irrespectively of the DNA sequence. 6. Complex formation with netropsin increases sensitivity to DNase I at certain DNA sites along with the protection effect observed at neighboring sites.     

Binding of an antitumor drug to DNA, Netropsin and C-G-C-G-A-A-T-T-BrC-G-C-G

The antitumor antibiotic netropsin has been co-crystallized with a double-helical B-DNA dodecanucleotide of sequence: C-G-C-G-A-A-T-T-BrC-G-C-G, and the structure of the complex has been solved by X-ray diffraction at a resolution of 2.2 A. The structure has been refined independently by Jack-Levitt and Hendrickson-Konnert least-squares methods, leading to a final residual error of 0.257 by the Jack-Levitt approach (0.211 for two-sigma data) or 0.248 by the Hendrickson-Konnert approach, with no significant difference between refined structures. The netropsin molecule displaces the spine of hydration and fits snugly within the minor groove in the A-A-T-T center. It widens the groove slightly and bends the helix axis back by 8 degrees, but neither unwinds nor elongates the double helix. The drug molecule is held in place by amide NH hydrogen bonds that bridge adenine N-3 and thymine O-2 atoms, exactly as with the spine of hydration. The requirement of A X T base-pairs in the binding site arises because the N-2 amino group of guanine would demand impermissibly close contacts with netropsin. It is proposed that substitution of imidazole for pyrrole in netropsin should create a family of "lexitropsins" capable of reading G X C-containing base sequences.     

Normal mode calculation of a netropsin–DNA complex: Effect of structural deformation on vibrational spectrum

Calculations of the normal mode spectrum of a netropsin-DNA complex, an isolated DNA helix of the same sequence, and a free netropsin molecule were carried out and compared with observations. We find that the frequency of modes of the complex that are primarily vibrations localized to the helix are relatively unchanged from the frequencies found for similar modes of the isolated helix. On the other hand, the frequencies of those modes of the complex that are found to be primarily in the netropsin are found to change significantly in the complex from that of the free netropsin. Further analysis indicates that the changes in the “netropsin modes” are primarily due to the conformational deformation induced in the netropsin by formation of the complex. All the frequencies and shifts in frequency are in good agreement with recent Raman measurements. © 1995 John Wiley & Sons, Inc.     

Free Energy Calculations Predict Sequence Specificity in DNA-Drug Complexes

Abstract

Molecular dynamics simulations of DNA-netropsin complexes in water were performed using the thermodynamic cycle-perturbation method to calculate the free energy difference between complexes with an adenine-containing binding site and corresponding complexes where adenines are replaced by 2,6-diaminopurines (dap). The calculations predict a free energy difference of 3.7±0.9 kcal/mol (at 300K) in favour of netropsin binding to an (AATT)₂ DNA sequence compared to a (dapdapTT)₂ sequence.     

DNA conformational effects on the interaction of netropsin with A-tract sequences

The influence of cosolutes and DNA sequence on the interaction of netropsin with three duplexes has been studied by isothermal titration calorimetry. In buffer, netropsin forms two complexes with a net stoichiometry of 1:1 in the minor groove of the oligonucleotide (GCGCGAATTCGCGC)2. One complex has a weaker affinity and is more enthalpically favored relative to the other one, consistent with previous studies [Freyer, M. W., et al. (2006) Biophys. Chem. 126, 186-196]. With the cosolutes betaine and 2-methyl-2,4-pentanediol, the enthalpy and heat capacity changes indicate that the complex with weaker affinity is disfavored relative to the complex with higher affinity. With (CGCGCAATTGCGCG)2, netropsin has one binding mode in buffer, and complex formation is not influenced by the cosolutes. The similarities of the enthalpy and heat capacity changes suggest that netropsin interacts similarly with these two oligonucleotides in the presence of cosolutes. The oligonucleotide (GCGCAAATTTGCGC)2 also forms two complexes with netropsin, and the complex with weaker affinity is again disfavored by the cosolutes. Thus, the interaction of netropsin with these A/T binding sites is influenced both by the bases adjacent to the binding site and by cosolutes. We suggest that these two factors influence the conformation of the minor-groove binding site of DNA.