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.     

Induction of DNA bending by bifunctional intercalating agents of the 7H-pyridocarbazole family

The structures and binding energetics of selected complexes formed between the deoxynucleotides d(CpGpGpCpG).d(CpGpCpCpG), d(CpGpApTpCpG)2, d(GpCpGpCpCpG).d(CpGpGpCpGpC), and d(CpGpCpCpCpG)2 with the DNA bifunctional intercalating agent ditercalinium and three of its rigid linking chain derivatives have been investigated theoretically by means of a molecular mechanics approach that takes into account nucleic acid flexibility, ligand flexibility and solvent dielectric effects (R. Lavery, in: Unusual DNA structures, eds S. Harvey and R. Wells (Pergamon, New York, 1988) p. 189; R. Lavery, in: DNA bending and curvature, eds W.K. Olson et al. (Adenine Press, New York, 1988) p. 191). The piperidinium chains of the bis-intercalating ligands are always located in the major groove of DNA. For the energy-minimized complexes the ligand proceeds to bind following preferentially the 5'-pyrimidine-purine-3' alternating sequence, thus dictating the number of internal exclusion sites. The complexes with three exclusion sites will present (i) a bending of the structure towards the major groove, and (ii) a non-ideal distribution of unwinding angles; complexes with less than three exclusion sites will remain essentially linear. The absence of a bend does not preclude other types of local deformations of the base-pairs such as inclination, buckle and tip. The proposed structures of the d(CpGpApTpCpG)2 complexes are in agreement with NMR structural results. The possible relevance of these findings to a previously proposed mode of interaction for ditercalinium-like DNA ligands is discussed.     

A theoretical investigation of the sequence specificity in the binding of the antitumor drug anthramycin to DNA

A theoretical study is presented concerning DNA-anthramycin adducts. By explicit energy minimisations using a semi-empirical energy formula and an advanced algorithm the structural properties and the energetics of this system are analysed. The results obtained demonstrate that the formation of a covalently bound adduct in which anthramycin is attached to the N2 site of a guanine within a DNA fragment is accompanied by a considerable change in the nucleic acid conformation as confirmed by recent experimental evidence. With the use of the "SIR" methodology for treating DNA flexibility the general features of this change are characterised. The sequence specificity of anthramycin binding is investigated and the important role of sequence dependent nucleic acid flexibility brought to light. This theoretical treatment thus provides new elements for the interpretation of the origins of ligand binding specificities.     

HIV-1跨膜蛋白gp41与N-取代吡咯衍生物的结合模式研究

采用分子对接,分子动力学(MD)模拟和分子力学/泊松-波尔兹曼溶剂可有面积方法与分子力学/广义伯恩溶剂可及面积方法(MM-PBSA/MM-GBSA),预测两种N-取代吡咯衍生物与HIV-1 跨膜蛋白gp41疏水口袋的结合模式与作用机理．分子对接采用多种受体构象,并从结果中选取几种可能的结合模式进行MD 模拟,然后通过MM-PBSA计算结合能的方法识别最优的结合模式. MM-PBSA计算结果表明,范德华相互作用是结合的主要驱动力,而极性相互作用决定了配体在结合过程中的取向．进一步的结合能分解显示,配体的羧基与gp41残基Arg579的静电相互作用对结合有重要贡献．上述工作为进一步优化N-取代吡咯衍生物类的HIV-1融合抑制剂建立了良好的理论基础．     

Synthesis, photophysical properties, and nucleic acid binding of phenanthridinium derivatives based on ethidium

A series of substituted phenanthridine derivatives has been synthesized by converting the amines at the 3- and 8-positions of ethidium bromide into guanidine, pyrrole, urea, and various substituted ureas. The resulting derivatives exhibit unique spectral properties that change upon binding nucleic acids. The compounds were analyzed for their ability to inhibit the HIV-1 Rev-Rev Response Element (RRE) interaction, as well as for their affinity to calf thymus DNA. One derivative (3,8-bis-urea-ethylenediamine-5-ethyl-6-phenylphenanthridinium trifuroracetate) has an enhanced affinity and specificity for HIV-1 RRE as compared to ethidium bromide. These results indicate that the nucleic acid affinity and specificity of an intercalating agent can be tuned by synthetic modification of its exocyclic amines.     

Nucleic acid controlled catalysts of carboxylic esters hydrolysis

Two Cu(2+)-binding ligands were covalently attached to termini of short DNAs. The optimal compound of this type forms a catalytically inert complex with Cu(2+). In the presence of a complementary nucleic acid the complex is decomposed forming products, which may catalyze hydrolysis of carboxylic acid esters. We have demonstrated that this process can be applied for sequence specific detection of nucleic acids.     

Bile acid interactions with rabbit ileal lipid binding protein and an engineered helixless variant reveal novel ligand binding properties of a versatile beta-clam shell protein scaffold

The intracellular ileal lipid binding proteins (ILBPs) are involved in the transport and enterohepatic circulation of bile acids. ILBPs from different species show high sequence and structural homology and have been shown to bind multiple bile acid ligands with differing degrees of selectivity and positive co-operativity. Human ILBP binds bile acid derivatives in a well-characterised 2:1 ligand:protein complex, however, we show that the highly homologous rabbit ILBP (82% sequence identity) with seven conservative substitutions preferentially binds multiple conjugated deoxycholate ligands in a novel 3:1 binding mode essentially within the same beta-clam shell structure. We have extended these studies to investigate the role of the alpha-helical capping motif (residues 9-35) in controlling the dimensions of the binding cavity and ligand uptake. Substituting the alpha-helical motif (residues 9-35) with a short Gly-Gly-Ser-Gly linker dramatically affects the protein stability such that under physiological conditions the mutant (Deltaalpha-ILBP) is highly disordered. However, we show that the inability of the mutant to adopt a stable three-dimensional structure under these conditions is no barrier to binding ligands with near-native affinity. These structural modifications not only demonstrate the possibility of strong coupling between ligand binding and protein folding, but result in changes in bile acid selectivity and binding stoichiometry, which we characterise in detail using isothermal calorimetry and mass spectrometry.     

Allosteric aptamers: targeted reversibly attenuated probes

Aptamers are unique nucleic acids with regulatory potentials that differ markedly from those of proteins. A significant feature of aptamers not possessed by proteins is their ability to participate in at least two different types of three-dimensional structure: a single-stranded folded structure that makes multiple contacts with the aptamer target and a double-helical structure with a complementary nucleic acid sequence. We have made use of this structural flexibility to develop an aptamer-based biosensor (a targeted reversibly attenuated probe, TRAP) in which hybridization of a cis-complementary regulatory nucleic acid (attenuator) controls the ability of the aptamer to bind to its target molecule. The central portion of the TRAP, between the aptamer and the attenuator, is complementary to a target nucleic acid, such as an mRNA, which is referred to as a regulatory nucleic acid (regNA) because it regulates the activity of the aptamer in the TRAP by hybridization with the central (intervening) sequence. The studies reported here of the ATP-DNA TRAP suggest that, as well as inhibiting the aptamer, the attenuator also acts as a structural guide, much like a chaperone, to promote proper folding of the TRAP such that it can be fully activated by the regDNA. We also show that activation of the aptamer in the TRAP by the complementary nucleic acid at physiological temperatures is sensitive to single-base mismatches. Aptamers that can be regulated by a specific nucleic sequence such as in an mRNA have potential for many in vivo applications including regulating a particular enzyme or signal transduction pathway or imaging gene expression in vivo.     

Stabilization of DNA Triple Helices Using Conjugates of Oligonucleotides and Synthetic Ligands

The possibility is discussed of stabilizing a DNA triple helix by covalent conjugation to the third strand (through its terminal phosphate) of ligands that have affinity to double and triple helices. Two types of stabilizers are considered: minor groove binders based on oligopyrroles, and triplex-specific intercalators. As a target, a synthetic 29-mer duplex containing a natural polypurine sequence of the human immunodeficiency provirus was employed. The stabilization with minor groove binders requires several conditions to be respected: a sufficiently long linker capable of reaching the minor groove from the major groove, a specific double-stranded structure of the oligopyrrole fragment, and its in-phase fitness to the target sequence. The best stabilizers of a triplex were novel conjugates in which two parallel molecules containing six pyrrole units each are linked to the same 5"-phosphate of a 16-mer triplex-forming oligonucleotide. The stabilizing properties of these derivatives were comparable to those of benzoindoloquinoline (BIQ) intercalators attached to the terminal phosphate of triple-helix forming oligonucleotides.     

Electron Correlated Ab Initio Study of Amino Group Flexibility for Improvement of Molecular Mechanics Simulations on Nucleic Acid Conformations and Interactions

High level ab initio studies demonstrate substantial conformational flexibility of amino groups of nucleic acid bases. This flexibility is important for biological functions of DNA. Existing force field models of molecular mechanics do not describe this phenomenon due to a lack of quantitative experimental data necessary for an adjustment of empirical parameters. We have performed extensive calculations of nucleic acid bases at the MP2/6-31G(d,p) level of ab initio theory for broad set of amino group configurations. Two-dimensional maps of energy and geometrical characteristics as functions of two amino hydrogen torsions have been constructed. We approximate the maps by polynomial expressions, which can be used in molecular mechanics calculations. Detailed considerations of these maps enable us to propose a method for determination of numerical coefficients in the developed formulae using restricted sets of points obtained via higher-level calculations.     

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.     

Theoretical Study of the Sequence Selectivity of Isolexins,Isohelical DNA Groove Binding Ligands. Proposal for the GC Minor Groove Specific Compounds

Abstract

Atheoretical study is presented of complex formation between DNA fragments of different base sequences and isolexins, “isohelical base reading polymers”, formed of heteroaromatic pentagonal rings joined by appropriate linkers. Extensive computations are performed for the isolexin composed of the furan-pyrrole-furan sequence. They involve charged ligands with propioamidinium groups at both ends as well eis neutral molecules with terminal methyl, carbonyl and amino groups. Two different groups (C=O and NH) are used as linkers between the base reading moieties. The role of these elements on the binding preference of the ligands has been examined. The results show that the mere possibility of formation of hydrogen bonds between a ligand and the nucleic acid bases is not sufficient to ensure its binding specificity which is determined largely by the interplay of electrostatic factors. Thus the dicationic isolexins uniformly prefer AT sequences. For the neutral isolexins the nature of the groups forming the linkers is a major factor in defining the specificity, although these groups do not participate directly in the interaction with DNA The C=O linkers favour binding to AT sequence while the N-H linkers permit preferential binding to the GAG sequence. Finally, for the first time in theoretical computations, a ligand is proposed which should bind preferentially to the minor groove of GC sequences: this ligand is a neutral isolexin composed of three furan rings linked by two N-H groups. This ligand is considered as an improvable prototype. Altogether the results presented open the path for the designing of minor groove ligands specific for any desirable DNA base sequence.     

Stabilization of DNA triple helix using conjugates of oligonucleotides and synthetic ligands

Possibility of stabilization of DNA triple helix is discussed using a covalent conjugation to the third strand (through its terminal phosphate) of ligands that have affinity to double and triple helices. Two types of stabilizers are considered: minor groove binders based on oligopyrroles and triplex-specific interacalators. As a target, a synthetic 29-mer duplex containing a natural polypurinic sequence of the human immunodeficiency provirus was employed. The stabilization with minor groove binders requires several conditions to be respected: a sufficiently long linker capable of reaching out the minor groove from the major one, a specific double-stranded structure of the oligopyrrole fragment and its in-phase fitness to the target sequence. The best stabilizers of a triplex turned out to be novel conjugates in which two parallel molecules containing six pyrrole units each are linked to the same 5'-phosphate of a 16-mer triplex-forming oligonucleotide. The stabilizing properties of these derivatives were comparable with those of benzoindoloquinoline (BIQ) intercalators attached to the terminal phosphate of triple-helix forming oligonucleotides.     

By-Product Formation in Repetitive PCR Amplification of DNA Libraries during SELEX

The selection of nucleic acid aptamers is an increasingly important approach to generate specific ligands binding to virtually any molecule of choice. However, selection-inherent amplification procedures are prone to artificial by-product formation that prohibits the enrichment of target-recognizing aptamers. Little is known about the formation of such by-products when employing nucleic acid libraries as templates. We report on the formation of two different forms of by-products, named ladder- and non-ladder-type observed during repetitive amplification in the course of in vitro selection experiments. Based on sequence information and the amplification behaviour of defined enriched nucleic acid molecules we suppose a molecular mechanism through which these amplification by-products are built. Better understanding of these mechanisms might help to find solutions minimizing by-product formation and improving the success rate of aptamer selection.     

多组分反应在吡咯衍生物合成中的应用进展

多组分反应(Multicomponent Coupling Reactions,MCRs)是将三种或三种以上的相对简单的原料加入"一锅煮"反应中,不经过中间体分离而直接得到结构复杂的分子的方法。具有操作简单、原料易得、环境友好等优点。多组分反应发展至今,已经在包括药物合成的多个领域中起到了重要作用。吡咯衍生物作为一类在自然界中广泛分布、有着重要生物医药价值和材料科学价值的重要化合物,其合成方法一直是有机合成研究中的热点之一。近年来用于合成吡咯衍生物的多组分反应也不断地被报道。多组分反应能够很好地解决传统吡咯合成路线中取代基引入的问题并能提高反应收率,为吡咯衍生物的合成与研究带来便利。本文从合成方法上概述了近年来多组分反应在吡咯合成中的应用,为合成已有吡咯衍生物和探寻全新结构的吡咯衍生物提供合成方法上的新思路。     

RNase activity of a DNA minor groove binder with a minimalist catalytic motif from RNase A

Imidazole and compounds containing imidazole residues have been shown to cleave RNA in an RNase A-mimicking manner. Di-imidazole lexitropsin is a compound which is derived from the polyamide drugs distamycin and netropsin essentially by the replacement of two pyrrole heterocycles with N-methyl-imidazole residues. This enables it to bind to the minor groove of B-DNA in a sequence-specific manner. We demonstrate here that this lexitropsin derivative has RNA cleavage activity, as tested on model RNAs. Optimal cleavage conditions and cleavage specificity resemble those known from other imidazole conjugates and are thus consistent with an RNase A type cleavage mechanism. The optimum concentration of the compound for cleavage is similar to previously investigated imidazole-based RNase mimics. As a whole new class of chemical compounds capable of interacting with nucleic acids through extensive hydrogen bonding, these imidazole containing compounds constitute promising scaffolds and ligands, for the construction of novel RNase mimics with high affinity.     

Minimal RNA Aptamer Sequences That Can Inhibit or Alleviate Noncompetitive Inhibition of the Muscle-Type Nicotinic Acetylcholine Receptor

Combinatorially synthesized nucleotide polymers have been used during the last decade to find ligands that bind to specific sites on biological molecules, including membrane-bound proteins such as the nicotinic acetylcholine receptors (nAChRs). The neurotransmitter receptors belong to a group of four structurally related proteins that regulate signal transmission between ~10¹¹ neurons of the mammalian nervous system. The nAChRs are inhibited by compounds such as the anticonvulsant MK-801 [(+)-dizocilpine] and abused drugs such as cocaine. Based on predictions arising from the mechanism of receptor inhibition by MK-801 and cocaine, we developed two classes of RNA aptamers: class I members, which inhibit the nAChR, and class II members, which alleviate inhibition of the receptor by MK-801 and cocaine. The systematic evolution of ligands by the exponential enrichment (SELEX) method was used to obtain these compounds. Here, we report that we have truncated RNA aptamers in each class to determine the minimal nucleic acid sequence that retains the characteristic function for which the aptamer was originally selected. We demonstrate that a truncated class I aptamer containing a sequence of seven nucleotides inhibits the nAChR and that a truncated class II aptamer containing a sequence of only four nucleotides can alleviate MK-801 inhibition.