Theoretical study of the binding of aliphatic diamines to the minor groove of a B-DNA (dA-dT)11 oligomer

Design and synthesis of a peptide derivative of ametantrone targeting the major groove of the d(GGCGCC)2palindromic sequence

We report the synthesis of a peptide derivative of antitumor anthraquinones, designed to target GC-rich palindromic sequences. It has micromolar activities on three cancer cell lines and is fifty times less toxic than mitoxantrone on a healthy line.

Anisotropic, Polarizable Molecular Mechanics Studies of Inter- and Intramolecular Interactions and Ligand-Macromolecule Complexes. A Bottom-Up Strategy

We present an overview of the SIBFA polarizable molecular mechanics procedure, which is formulated and calibrated on the basis of quantum chemistry (QC). It embodies nonclassical effects such as electrostatic penetration, exchange-polarization, and charge transfer. We address the issues of anisotropy, nonadditivity, and transferability by performing parallel QC computations on multimolecular complexes. These encompass multiply H-bonded complexes and polycoordinated complexes of divalent cations. Recent applications to the docking of inhibitors to Zn-metalloproteins are presented next, namely metallo-beta-lactamase, phosphomannoisomerase, and the nucleocapsid of the HIV-1 retrovirus. Finally, toward third-generation intermolecular potentials based on density fitting, we present the development of a novel methodology, the Gaussian electrostatic model (GEM), which relies on ab initio-derived fragment electron densities to compute the components of the total interaction energy. As GEM offers the possibility of a continuous electrostatic model going from distributed multipoles to densities, it allows an inclusion of short-range quantum effects in the molecular mechanics energies. The perspectives of an integrated SIBFA/GEM/QM procedure are discussed.

DNA compaction by divalent cations: structural specificity revealed by the potentiality of designed quaternary diammonium salts

DNA interaction with quaternary diammonium dications, R(CH(3))(2)N(+)(CH(2))(n)N(+)(CH(3))(2)R, having various intercharge distances, lengths, and branching, and the chemical nature of the hydrophobic substituents were investigated by fluorescent microscopy and circular dichroism (CD) spectroscopy to reveal their structural specificity for binding to DNA. The conformational behavior of DNA was found to be highly sensitive to the structure of the dications with separated charges. The distance between two ammonium groups greatly influences the compaction activity of the dications. To explain this situation, we proposed a model that demonstrates that the charge density of the dication and the geometric fit between DNA phosphates and the ammonium groups in the dications play an important role in providing efficient DNA collapse. Elongation of the alkyl substituents (R) in the diammonium salts from ethyl to hexyl did not generate any significant alterations in the compaction activities, whereas the branching of substituents caused a drastic decrease in their compaction ability. Based on the results of CD spectroscopy, it was found that the ability of the dications to provoke a DNA transition from the B-form to A-form was also specific: it depended on their intercharge distances and was independent of the length of alkyl substituents.

Theoretical design of a bis-orthopepetide derivative of a tetracationic porphyrin targeted toward a six-base pair sequence of DNA

Tetra-(4-N-methylpyridyl)-porphyrin, (T4MPyP), is a tetracationic porphyrin that binds to G-C sequences of DNA by means of an intercalative mode. In order to extend its selective sequence recognition capacity for bases beyond the intercalation site, and in the major groove, we have undertaken the theoretical design of bis-ortho peptide derivatives of T4MPyP. In these, two ortho N-methylpiperidinium nitrogens are linked to a cationic residue, L-Lys, L-Orn, or L-Arg. The binding energetics of these novel compounds were compared for six distinct double-stranded palindromic hexanucleotide sequences. Four distinct modes of binding were compared: a) major, b) minor groove binding of both peptidic arms; c) a straddling mode in which each arm is in a different groove of DNA; d) exclusive binding of the arms to the sugar-phosphate backbone. For our most promising compound, that with Lys side-chains, a distinctive energetical advantage was computed in favor of an all-major groove binding to sequence d(CCCGGG)2. The corresponding complex is separated by an energy gap of 12 kcal/mol, with respect to the second-best sequence bound in the major groove, d(GGCGCC)2, and of 20 kcal/mol with respect to minor groove binding to sequence d(TACGTA)2. The results obtained with such a prototypic compound indicate that it is fully possible to design sequence selective (> 6 base-pairs) photosensitizers as peptide derivatives of T4MPyP and prompt the engineering of further, more complex analogs thereof.

DNA structure and perturbation by drug binding

Images

The flexibility of the nucleic acids: (III). The interaction of an aliphatic diamine, putrescine, with flexible B-DNA

A theoretical modelling of the interaction of putrescine (H3+N-(CH2)4-(+NH3) with DNA is carried out, introducing two new features which make the simulation of this interaction considerably more realistic. Firstly, the DNA to which putrescine is bound is fully flexible and thus able to respond to the distorting influence of the ligand. Secondly, the effect of changing the ratio of DNA base pairs per bound ligand is explicitly modelled. In this way, we have been able to confirm the experimentally known preference of putrescine binding with AT base pairs in B-DNA, but we also show, through the new features introduced, that the nature of the binding site of the ligand and the resulting impact on DNA conformation is strongly modified by the ligand binding density.

A theoretical investigation on the sequence selective binding of adriamycin to double-stranded polynucleotides

Theoretical computations are performed on the structural and energetical factors involved in the sequence selective binding of adriamycin (ADM) to five self-complementary double-stranded hexanucleotides. Among the two regularly alternating hexanucleotides d (TATATA)2 and d (CGCGCG)2, a stronger binding is predicted for the former. The strongest complex is computed, however, for the mixed hexanucleotide d (CGTACG)2, containing the intercalation site between two CG base pairs and an adjacent TA base pair. The overall sequence preference is the result of an intricate interplay of sequence preferences of the constituents in particular of daunosamine and the 9-OH substituent. Altogether, the selective base pair recognition by adriamycin cannot be defined in terms of the two base pairs implicated in the intercalation site alone but must be expressed in terms of a triplet of base pairs.

A theoretical investigation on the sequence selective binding of daunomycin to double-stranded polynucleotides

Theoretical computations are performed on the structural and energetical factors involved in the sequence selective binding of daunomycin (DNM) to six representative self-complementary double-stranded hexanucleotides: d(CGTACG)2,d(CGATCG)2,d(CITACI)2, d(TATATA)2, d(CGCGCG)2 and d(TACGTA)2. The conformational angles of the hexanucleotides are fixed in values found in the representative crystal structure of the d(CGTACG)2-DNM complex. The intermolecular DNM-hexanucleotide interaction energies and the conformational energy changes of DNM upon binding are computed and optimized in the framework of the SIBFA procedure, which uses empirical formulas based on ab initio SCF computations. Among the two regularly alternating hexanucleotides, d(TATATA)2 and d(CGCGCG)2, a stronger binding is predicted for the former, in agreement with experimental results obtained with poly(dA-dT).poly(dA-dT) and poly(dG-dC).poly(dG-dC). Altogether, however, among the six investigated sequences, the strongest complexes are computed for the mixed hexanucleotides d(CGATCG)2 and d(CGTACG)2, containing the intercalation site between two CG base pairs and an adjacent TA base pair. This situation may be related to the increased affinity of DNM for GC rich DNA's and to the situation in the crystal structure of the DNM-d(CGTACG)2 complex. Analysis of the intrinsic base sequence preferences expressed by the individual constituents of DNM, namely the daunosamine side chain, the chromophore ring and its two 9-hydroxyl and 9-acetoxy substituents, reveals that the overall sequence preference found is the result of a rather intricate interplay of intrinsic sequence preferences, in particular at the level of daunosamine and the 9-hydroxyl substituent.(ABSTRACT TRUNCATED AT 250 WORDS)