Quantifying force-dependent and zero-force DNA intercalation by single-molecule stretching

We used single DNA molecule stretching to investigate DNA intercalation by ethidium and three ruthenium complexes. By measuring ligand-induced DNA elongation at different ligand concentrations, we determined the binding constant and site size as a function of force. Both quantities depend strongly on force and, in the limit of zero force, converge to the known bulk solution values, when available. This approach allowed us to distinguish the intercalative mode of ligand binding from other binding modes and allowed characterization of intercalation with binding constants ranging over almost six orders of magnitude, including ligands that do not intercalate under experimentally accessible solution conditions. As ligand concentration increased, the DNA stretching curves saturated at the maximum amount of ligand intercalation. The results showed that the applied force partially relieves normal intercalation constraints. We also characterized the flexibility of intercalator-saturated dsDNA for the first time.     

The thermodynamics of polyamide-DNA recognition: hairpin polyamide binding in the minor groove of duplex DNA

Crescent-shaped synthetic ligands containing aromatic amino acids have been designed for specific recognition of predetermined DNA sequences in the minor groove of DNA. Simple rules have been developed that relate the side-by-side pairings of Imidazole (Im) and Pyrrole (Py) amino acids to their predicted target DNA sequences. We report here thermodynamic characterization of the DNA-binding properties of the six-ring hairpin polyamide, ImImPy-gamma-PyPyPy-beta-Dp (where gamma = gamma-aminobutyric acid, beta = beta-alanine, and Dp = dimethylaminopropylamide). Our data reveal that, at 20 degrees C, this ligand binds with a relatively modest 1.8-fold preference for the designated match site, 5'-TGGTA-3', over the single base pair mismatch site, 5'-TGTTA-3'. By contrast, we find that the ligand exhibits a 102-fold greater affinity for its designated match site relative to the double base pair mismatch site, 5'-TATTA-3'. These results demonstrate that the energetic cost of binding to a double mismatch site is not necessarily equal to twice the energetic cost of binding to a single mismatch site. Our calorimetrically measured binding enthalpies and calculated entropy data at 20 degrees C reveal the ligand sequence specificity to be enthalpic in origin. We have compared the DNA-binding properties of ImImPy-gamma-PyPyPy-beta-Dp with the hairpin polyamide, ImPyPy-gamma-PyPyPy-beta-Dp (an Im --> Py "mutant"). Our data reveal that both ligands exhibit high affinities for their designated match sites, consistent with the Dervan pairing rules. Our data also reveal that, relative to their corresponding single mismatch sites, ImImPy-gamma-PyPyPy-beta-Dp is less selective than ImPyPy-gamma-PyPyPy-beta-Dp for its designated match site. This result suggests, at least in this case, that enhanced binding affinity can be accompanied by some loss in sequence specificity. Such systematic comparative studies allow us to begin to establish the thermodynamic database required for the rational design of synthetic polyamides with predictable DNA-binding affinities and specificities.     

Enthalpy/entropy compensation: influence of DNA flanking sequence on the binding of 7-amino actinomycin D to its primary binding site in short DNA duplexes

The effect of the context of the flanking sequence on ligand binding to DNA oligonucleotides that contain consensus binding sites was investigated for the binding of the intercalator 7-amino actinomycin D. Seven self-complementary DNA oligomers each containing a centrally located primary binding site, 5'-A-G-C-T-3', flanked on either side by the sequences (AT)(n) or (AA)(n) (with n = 2, 3, 4) and AA(AT)(2), were studied. For different flanking sequences, (AA)(n)-series or (AT)(n)-series, differential fluorescence enhancements of the ligand due to binding were observed. Thermodynamic studies indicated that the flanking sequences not only affected DNA stability and secondary structure but also modulated ligand binding to the primary binding site. The magnitude of the ligand binding affinity to the primary site was inversely related to the sequence dependent stability. The enthalpy of ligand binding was directly measured by isothermal titration calorimetry, and this made it possible to parse the binding free energy into its energetic and entropic terms. Our results reveal a pronounced enthalpy-entropy compensation for 7-amino actinomycin D binding to this family of oligonucleotides and suggest that the DNA sequences flanking the primary binding site can strongly influence ligand recognition of specific sites on target DNA molecules.     

Simultaneous recognition of nucleobase and sites of DNA damage: effect of tethered cation on the binding affinity

BACKGROUND: The 3,5-diamino-N-(3-aminopropyl)-6-chloropyrazine-2-carboxamide (DCPC-NH(2)) has been synthesized and characterized by Mass and (1)H NMR. The selective binding of the ligand to thymine (T) target base is investigated by the melting temperature (T(m)) and fluorescence measurements. METHODS: Thermal denaturation study of DNA duplex containing T target base revealed the DeltaT(m) of 5.1 degrees C, while least influence was observed for other target bases. The fluorescence of the ligand DCPC-NH(2) is quenched only upon adding the DNA containing T target base. RESULTS: The binding constant for the interaction of the ligand to T target base containing DNA duplex was determined to be 4.7 (+/-0.3)x10(6) M(-1). The tethered cation in the ligand is found to enhance the binding constant. The ligand binds to both a target nucleotide and an AP site on the complimentary strand for the target strand in a DNA duplex. GENERAL SIGNIFICANCE: Interestingly, the electronic behavior of the ligand depends on the bases flanking the AP site. Its fluorescence is quenched with guanine flanking bases, while it is enhanced with DNA duplex containing T bases flanking an AP site. Finally, the binding modes were visualized by molecular modeling.     

Distamycin inhibition of topoisomerase I-DNA interaction: a mechanistic analysis.

Inhibition of eukaryotic DNA topoisomerase I by the minor groove binding ligand, distamycin A, was investigated. Low concentrations of the ligand selectively prevented catalytic action at a high affinity topoisomerase I binding sequence. A restriction enzyme protection assay indicated that the catalytic cycle was blocked at the binding step. Distamycin binding sites on DNA were localized by hydroxyl radical footprinting. A strongly preferred site mapped to a homopolymeric (dA).(dT)-tract partially included in the essential topoisomerase I binding region. Mutational elimination of the stable helix curvature associated with this ligand binding site demonstrated that (i) the intrinsic bend was unessential for efficient binding of topoisomerase I, and (ii) distamycin inhibition did not occur by deformation of a stable band. Alternative modes of inhibition are discussed.     

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.     

Mixed mode of ligand-DNA binding results in S-shaped binding curves

S-shaped binding curves often characterize interactions of ligands with nucleic acid molecules as analyzed by different physico-chemical and biophysical techniques. S-shaped experimental binding curves are usually interpreted as indicative of the positive cooperative interactions between the bound ligand molecules. This paper demonstrates that S-shaped binding curves may occur as a result of the "mixed mode" of DNA binding by the same ligand molecule. Mixed mode of the ligand-DNA binding can occur, for example, due to 1) isomerization or dimerization of the ligands in solution or on the DNA lattice, 2) their ability to intercalate the DNA and to bind it within the minor groove in different orientations. DNA-ligand complexes are characterized by the length of the ligand binding site on the DNA lattice (so-called "multiple-contact" model). We show here that if two or more complexes with different lengths of the ligand binding sites could be produced by the same ligand, the dependence of the concentration of the complex with the shorter length of binding site on the total concentration of ligand should be S-shaped. Our theoretical model is confirmed by comparison of the calculated and experimental CD binding curves for bis-netropsin binding to poly(dA-dT) poly(dA-dT). Bis-netropsin forms two types of DNA complexes due to its ability to interact with the DNA as monomers and trimers. Experimental S-shaped bis-netropsin-DNA binding curve is shown to be in good correlation with those calculated on the basis of our theoretical model. The present work provides new insight into the analysis of ligand-DNA binding curves.     

Footprinting of Chlorella virus DNA ligase bound at a nick in duplex DNA.

The 298-amino acid ATP-dependent DNA ligase of Chlorella virus PBCV-1 is the smallest eukaryotic DNA ligase known. The enzyme has intrinsic specificity for binding to nicked duplex DNA. To delineate the ligase-DNA interface, we have footprinted the enzyme binding site on DNA and the DNA binding site on ligase. The size of the exonuclease III footprint of ligase bound a single nick in duplex DNA is 19-21 nucleotides. The footprint is asymmetric, extending 8-9 nucleotides on the 3'-OH side of the nick and 11-12 nucleotides on the 5'-phosphate side. The 5'-phosphate moiety is essential for the binding of Chlorella virus ligase to nicked DNA. Here we show that the 3'-OH moiety is not required for nick recognition. The Chlorella virus ligase binds to a nicked ligand containing 2',3'-dideoxy and 5'-phosphate termini, but cannot catalyze adenylation of the 5'-end. Hence, the 3'-OH is important for step 2 chemistry even though it is not itself chemically transformed during DNA-adenylate formation. A 2'-OH cannot substitute for the essential 3'-OH in adenylation at a nick or even in strand closure at a preadenylated nick. The protein side of the ligase-DNA interface was probed by limited proteolysis of ligase with trypsin and chymotrypsin in the presence and absence of nicked DNA. Protease accessible sites are clustered within a short segment from amino acids 210-225 located distal to conserved motif V. The ligase is protected from proteolysis by nicked DNA. Protease cleavage of the native enzyme prior to DNA addition results in loss of DNA binding. These results suggest a bipartite domain structure in which the interdomain segment either comprises part of the DNA binding site or undergoes a conformational change upon DNA binding. The domain structure of Chlorella virus ligase inferred from the solution experiments is consistent with the structure of T7 DNA ligase determined by x-ray crystallography.     

Binding isotherms and the interaction between proflavine and a DNA of high G · C content

Binding isotherms corresponding to several situations of ligand binding to a linear polymer are calculated, including cases of cooperativity or anticooperativity between the bound ligand states, multiple binding modes that are competitive or non competitive, and possible exclusion of an arbitrary number of adjacent sites upon occupancy of a site by a single ligand. The sequence generating function method of Lifson and Bradley is used, requiring the assumption that no end effects are involved. The case of strong binding of the dye proflavine to a DNA of high G. C content, that of M. lysodeikticus, is considered in detail, and a single model capable of reconciling the available kinetic and equilibrium data on this system, involving two competing binding modes, is discussed.     

Cooperation effects in binding of large ligands to DNA. II. Contact interactions between adsorbed ligands

Cooperative effects arising upon binding of biologically active ligands to DNA are considered. Equations are derived which enable one to describe the binding of two different ligands to DNA. We also consider the case when ligand can form two type of DNA complexes. The cooperative binding of the ligand in the vicinity of saturation level of binding can be described with a good accuracy by equation derived for the non-cooperative adsorption of the same ligand with some effective binding constant Keff. It is shown that cooperative effects arising upon binding of proteins and other ligands to DNA can be divided into two groups depending on the symmetry of interactions between the bound ligand molecules. In particular, if such interactions favor the formation of dimeric ligand species on the DNA, Keff approximately a1/2, where a is the ligand-ligand interaction constant. If cooperative interactions favor the formation of aggregates of unrestricted size, then Keff approximately aL+Y, where L is the size of the binding site for the ligand on DNA.     

Mixed Mode of Ligand-DNA Binding Results in S-Shaped Binding Curves

Abstract

S-shaped binding curves often characterize interactions of ligands with nucleic acid molecules as analyzed by different physicochemical and biophysical techniques. S-shaped experimental binding curves are usually interpreted as indicative of the positive cooperative interactions between the bound ligand molecules. This paper demonstrates that S-shaped binding curves may occur as a result of the “mixed mode” of DNA binding by the same ligand molecule. Mixed mode of the ligand-DNA binding can occur, for example, due to 1) isomerization or dimerization of the ligands in solution or on the DNA lattice, 2) their ability to intercalate the DNA and to bind it within the minor groove in different orientations. DNA- ligand complexes are characterized by the length of the ligand binding site on the DNA lattice (so-called “multiple-contact” model). We show here that if two or more complexes with different lengths of the ligand binding sites could be produced by the same ligand, the dependence of the concentration of the complex with the shorter length of binding site on the total concentration of ligand should be S-shaped. Our theoretical model is confirmed by comparison of the calculated and experimental CD binding curves for bis-netropsin binding to poly(dA-dT) poly(dA-dT). Bis-netropsin forms two types of DNA complexes due to its ability to interact with the DNA as monomers and trimers. Experimental S-shaped bis-netropsin-DNA binding curve is shown to be in good correlation with those calculated on the basis of our theoretical model. The present work provides new insight into the analysis of ligand-DNA binding curves.     

Holo-sterol carrier protein-2. (13)C NMR investigation of cholesterol and fatty acid binding sites

Although sterol carrier protein-2 (SCP-2) stimulates sterol transfer in vitro, almost nothing is known regarding the identity of the putative cholesterol binding site. Furthermore, the interrelationship(s) between this SCP-2 ligand binding site and the recently reported SCP-2 long chain fatty acid (LCFA) and long chain fatty acyl-CoA (LCFA-CoA) binding site(s) remains to be established. In the present work, two SCP-2 ligand binding sites were identified. First, both [4-(13)C]cholesterol and 22-(N-(7-nitrobenz-2-oxa-1, 3-diazol-4-yl)amino)-23,24-bisnor-5-cholen-3beta-ol (NBD-cholesterol) binding assays were consistent with a single cholesterol binding site in SCP-2. This ligand binding site had high affinity for NBD-cholesterol, K(d) = 4.15 +/- 0.71 nM. (13)C NMR-labeled ligand competition studies demonstrated that the SCP-2 high affinity cholesterol binding site also bound LCFA or LCFA-CoA. However, only the LCFA-CoA was able to effectively displace the SCP-2-bound [4-(13)C]cholesterol. Thus, the ligand affinities at this SCP-2 binding site were in the relative order cholesterol = LCFA-CoA > LCFA. Second, (13)C NMR studies demonstrated the presence of another ligand binding site on SCP-2 that bound either LCFA or LCFA-CoA but not cholesterol. Photon correlation spectroscopy was consistent with SCP-2 being monomeric in both liganded and unliganded states. In summary, both (13)C NMR and fluorescence techniques demonstrated for the first time that SCP-2 had a single high affinity binding site that bound cholesterol, LCFA, or LCFA-CoA. Furthermore, results with (13)C NMR supported the presence of a second SCP-2 ligand binding site that bound either LCFA or LCFA-CoA but not cholesterol. These data contribute to our understanding of a role for SCP-2 in both cellular cholesterol and LCFA metabolism.     

Footprinting of echinomycin and actinomycin D on DNA molecules asymmetrically substituted with inosine and/or 2,6-diaminopurine.

In order to clarify the role of the purine 2-amino group in the recognition of DNA by small molecules we have examined the binding of actinomycin D and echinomycin to artificial DNA molecules asymmetrically substituted with inosine and/or 2,6-diaminopurine (DAP) in one of the complementary strands. These DNAs, prepared by a method based upon PCR, present various potential sites for antibiotic binding, including several containing only a single purine 2-amino group in different configurations. The results show unambiguously that the presence of two 2-amino groups is mandatory for binding of actinomycin D to double-stranded DNA. In the case of echinomycin only one purine 2-amino group is required for remarkably strong binding to the asymmetric TpDAP.TpA dinucleotide step, but the CpDAP.TpI step (which also contains only a single purine-2 amino group) does not afford a binding site. Evidently, removing a 2-amino group (G-->I substitution) is dominant over adding one (A-->DAP substitution). No sequences containing just a single guanine residue are acceptable. The possibility is raised that replacing guanosine with inosine may do more than remove a group endowed with hydrogen bonding capability and interfere with ligand binding in other ways. The new methodology developed to construct asymmetrically substituted DNA substrates for this work provides a novel strategy that should be generally applicable for studying ligand-DNA interactions, beyond the specific interest in drug binding to DNA, and may help to elucidate how proteins and oligonucleotides recognize their target sites.     

Determinants of the position of a Flp-induced DNA bend.

The Flp site-specific recombinase from Saccharomyces cerevisiae induces DNA bending upon interaction with the Flp recognition target (FRT) site. The minimal FRT site is comprised of two inverted binding elements which flank a central core region. Binding of a single monomer of Flp to DNA induces a DNA bend of 60 degrees. The position of this bend differed depending on whether the substrate contained a single binding element or a two-element FRT site. In the present work we tested and disproved a model in which a single Flp monomer interacts with both symmetry elements of a single FRT site. Likewise, we showed that a model in which a Flp monomer dissociates from a singly occupied FRT site and reassociates with the unbound element of another singly occupied FRT site during electrophoresis, does not account for the apparent shift in the position of the bend centre. It seems that the movement of a Flp monomer between the a and b elements of one FRT site during electrophoresis accounts for this anomaly. The position of the DNA bend resulting from the association of a Flp monomer with the FRT site is also influenced by the DNA sequences flanking the site. We conclude that attempts to measure the bend centre of a complex of one Flp molecule bound to a DNA containing two binding elements give misleading results. The position of the bend is more accurately measured in the presence of a single binding element.     

Theoretical analysis of the footprinting experiment

In the footprinting experiment, an end-radiolabeled DNA restriction fragment is subjected to digest by an endonuclease in the presence and absence of a ligand which alters the endonuclease cleavage rate at sites of ligand-DNA contact. The location of these sites, and the strength of the ligand binding, are then deduced from the measured concentrations of the different oligonucleotides produced by the digest. We analyze the experiment in terms of coupled kinetic equations which take into account the cutting rates of endonuclease for sites with ligand present and absent, and the rates of binding and dissociation of the ligand to a site. As long as the ligand concentration remains essentially constant (which occurs, for example, if digest is terminated early enough to assure that all fragments result from single cuts by the endonuclease), the oligonucleotide concentrations reflect only the ligand binding equilibrium constant (ratio of rate constants) and the cutting rates in the presence and absence of ligand. We also show how the measured oligonucleotide concentrations (from, e.g. an autoradiogram) can be used to deduce the ligand equilibrium binding constants for the various sites on the polymer.     

Interaction of berenil with the tyrT DNA sequence studied by footprinting and molecular modelling. Implications for the design of sequence-specific DNA recognition agents.

We have developed a technique of partially-restrained molecular mechanics enthalpy minimisation which enables the sequence-dependence of the DNA binding of a non-intercalating ligand to be studied for arbitrary sequences of considerable length (greater than = 60 base-pairs). The technique has been applied to analyse the binding of berenil to the minor groove of a 60 base-pair sequence derived from the tyrT promoter; the results are compared with those obtained by DNAse I and hydroxyl radical footprinting on the same sequence. The calculated and experimentally observed patterns of binding are in good agreement. Analysis of the modelling data highlights the importance of DNA flexibility in ligand binding. Further, the electrostatic component of the interaction tends to favour binding to AT-rich regions, whilst the van der Waals interaction energy term favours GC-rich ones. The results also suggest that an important contribution to the observed preference for binding in AT-rich regions arises from lower DNA perturbation energies and is not accompanied by reduced DNA structural perturbations in such sequences. It is therefore concluded that those modes of DNA distortion favourable to binding are probably more flexible in AT-rich regions. The structure of the modelled DNA sequence has also been analysed in terms of helical parameters. For the DNA energy-minimised in the absence of berenil, certain helical parameters show marked sequence-dependence. For example, purine-pyrimidine (R-Y) base pairs show a consistent positive buckle whereas this feature is consistently negative for Y-R pairs. Further, CG steps show lower than average values of slide while GC steps show lower than average values of rise. Similar analysis of the modelling data from the calculations including berenil highlights the importance of DNA flexibility in ligand binding. We observe that the binding of berenil induces characteristic responses in different helical parameters for the base-pairs around the binding site. For example, buckle and tilt tend to become more negative to the 5'-side of the binding site and more positive to the 3'-side, while the base steps at either side of the centre of the site show increased twist and decreased roll.     

Thermodynamic compensation upon binding to exosite 1 and the active site of thrombin

Several lines of experimental evidence including amide exchange and NMR suggest that ligands binding to thrombin cause reduced backbone dynamics. Binding of the covalent inhibitor dPhe-Pro-Arg chloromethyl ketone to the active site serine, as well as noncovalent binding of a fragment of the regulatory protein, thrombomodulin, to exosite 1 on the back side of the thrombin molecule both cause reduced dynamics. However, the reduced dynamics do not appear to be accompanied by significant conformational changes. In addition, binding of ligands to the active site does not change the affinity of thrombomodulin fragments binding to exosite 1; however, the thermodynamic coupling between exosite 1 and the active site has not been fully explored. We present isothermal titration calorimetry experiments that probe changes in enthalpy and entropy upon formation of binary ligand complexes. The approach relies on stringent thrombin preparation methods and on the use of dansyl-l-arginine-(3-methyl-1,5-pantanediyl)amide and a DNA aptamer as ligands with ideal thermodynamic signatures for binding to the active site and to exosite 1. Using this approach, the binding thermodynamic signatures of each ligand alone as well as the binding signatures of each ligand when the other binding site was occupied were measured. Different exosite 1 ligands with widely varied thermodynamic signatures cause a similar reduction in ΔH and a concomitantly lower entropy cost upon DAPA binding at the active site. The results suggest a general phenomenon of enthalpy-entropy compensation consistent with reduction of dynamics/increased folding of thrombin upon ligand binding to either the active site or exosite 1.     

Double threading through DNA: NMR structural study of a bis-naphthalene macrocycle bound to a thymine-thymine mismatch

The macrocyclic bis-naphthalene macrocycle (2,7-BisNP), belonging to the cyclobisintercalator family of DNA ligands, recognizes T-T mismatch sites in duplex DNA with high affinity and selectivity, as evidenced by thermal denaturation experiments and NMR titrations. The binding of this macrocycle to an 11-mer DNA oligonucleotide containing a T-T mismatch was studied using NMR spectroscopy and NMR-restrained molecular modeling. The ligand forms a single type of complex with the DNA, in which one of the naphthalene rings of the ligand occupies the place of one of the mismatched thymines, which is flipped out of the duplex. The second naphthalene unit of the ligand intercalates at the A-T base pair flanking the mismatch site, leading to encapsulation of its thymine residue via double stacking. The polyammonium linking chains of the macrocycle are located in the minor and the major grooves of the oligonucleotide and participate in the stabilization of the complex by formation of hydrogen bonds with the encapsulated thymine base and the mismatched thymine remaining inside the helix. The study highlights the uniqueness of this cyclobisintercalation binding mode and its importance for recognition of DNA lesion sites by small molecules.     

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.