Understanding selectivity of hard and soft metal cations within biological systems using the subvalence concept. I. Application to blood coagulation: direct cation-protein electronic effects vs. indirect interactions through water networks.

Following a previous study by de Courcy et al. ((2009) Interdiscip. Sci. Comput. Life Sci. 1, 55-60), we demonstrate in this contribution, using quantum chemistry, that metal cations exhibit a specific topological signature in the electron localization of their density interacting with ligands according to its "soft" or "hard" character. Introducing the concept of metal cation subvalence, we show that a metal cation can split its outer-shell density (the so-called subvalent domains or basins) according to it capability to form a partly covalent bond involving charge transfer. Such behaviour is investigated by means of several quantum chemical interpretative methods encompasing the topological analysis of the Electron Localization Function (ELF) and Bader's Quantum Theory of Atoms in Molecules (QTAIM) and two energy decomposition analyses (EDA), namely the Restricted Variational Space (RVS) and Constrained Space Orbital Variations (CSOV) approaches. Further rationalization is performed by computing ELF and QTAIM local properties such as electrostatic distributed moments and local chemical descriptors such as condensed Fukui Functions and dual descriptors. These reactivity indexes are computed within the ELF topological analysis in addition to QTAIM offering access to non atomic reactivity local index, for example on lone pairs. We apply this "subvalence" concept to study the cation selectivity in enzymes involved in blood coagulation (GLA domains of three coagulation factors). We show that the calcium ions are clearly able to form partially covalent charge transfer networks between the subdomain of the metal ion and the carboxylate oxygen lone pairs whereas magnesium does not have such ability. Our analysis also explains the different role of two groups (high affinity and low affinity cation binding sites) present in GLA domains. If the presence of Ca(II) is mandatory in the central "high affinity" region to conserve a proper folding and a charge transfer network, external sites are better stabilised by Mg(II), rather than Ca(II), in agreement with experiment. The central role of discrete water molecules is also discussed in order to understand the stabilities of the observed X-rays structures of the Gla domain. Indeed, the presence of explicit water molecules generating indirect cation-protein interactions through water networks is shown to be able to reverse the observed electronic selectivity occuring when cations directly interact with the Gla domain without the need of water.

Importance of lone pair interactions/redistribution in hard and soft ligands within the active site of alcohol dehydrogenase Zn-metalloenzyme: Insights from electron localization function.

As a continuation of our previous work (de Courcy et al., 2008. J. Chem. Theo. Comput. 4 1659), lone pair-cation interactions were quantum-mechanically studied within the active site of the alcohol dehydrogenase Zn(II)-metalloenzyme by means of the topological analysis of the Electron Localization Function (ELF) and the Reduced Variational Space (RVS) energy decomposition analysis. Ligands lone pairs in direct interaction with the metal were shown to control the physical nature of the interaction as it appears to be dominated by polarization when the number of interacting lone pairs increases. Furthermore, we observed a peculiar behaviour of the cysteinate S(-) lone pairs which can redistribute and merge, thereby reducing their number to accommodate the zinc cation which also exhibits a consequent plasticity of its density outer shells which can delocalize towards ligands. Such observations should allow a deeper understanding of the usual softness/hardness concept of ions and ligands.

Design of peptoid analogue dimers and measure of their affinity for Grb2 SH3 domains.

This paper describes the design of the highest affinity ligands for Grb2 SH3 domains reported so far. These compounds were designed by combining N-alkyl amino acid incorporation in a proline-rich sequence with subsequent dimerization of the peptoid sequence based on structural data and molecular modeling. Optimization of the linker size is discussed, and the N-alkyl amino acid incorporation into both monomeric halves is reported. Because the affinity for Grb2 of the optimized compounds was too high to be measured using the fluorescent modifications that they induce on the Grb2 emission spectrum, a competition assay was developed. In this test, Grb2 is pulled down from a cellular extract by the initial VPPPVPPRRR peptide bound to Sepharose beads. In the presence of competitors, the test quantifies the amount of Grb2 displaced from the beads. It has enabled us to determine a K(i) value in the 10(-10) M range for the highest affinity Grb2 peptoid analogue dimer.

Small peptides containing phosphotyrosine and adjacent alphaMe-phosphotyrosine or its mimetics as highly potent inhibitors of Grb2 SH2 domain.

A series of small peptides with the sequence mAZ-pTyr-Xaa-Asn-NH(2), where Xaa denotes alpha-methylphosphotyrosine or its carboxylic mimetics, were synthesized as inhibitors of the Grb2 SH2 domain. Peptide 3 with (alpha-Me)pTyr as Xaa has the highest affinity for Grb2 (K(d) = 3 +/- 1 nM) and exhibits to date the best inhibitory activity (IC(50) = 11 +/- 1 nM) to displace PSpYVNVQN-Grb2 interaction in an ELISA test. The lower affinities of peptides with (alpha-Me)Tyr, (alpha-Me)Phe(4-CO(2)H), or (alpha-Me)Phe(4-CH(2)CO(2)H) as Xaa demonstrate the importance of a double charged phosphate group at the pY+1 position. Molecular modeling showed additional hydrogen bond interactions provided by the (alpha-Me)pTyr residue with the Grb2 SH2 domain. These results thus show that the effect of hydrophobic pY+1 residues, initially put forth to increased the binding affinities, can be further enhanced by a (-Me)pTyr residue which has both hydrophobic and hydrophilic properties.

A mimic of HIV-1 nucleocapsid protein impairs reverse transcription and displays antiviral activity.

Combined inhibition of HIV-1 reverse transcriptase and protease has significantly improved the treatment of AIDS. Nevertheless, resistance to these drugs occurs rapidly because of viral mutations, emphasizing the importance of identifying novel retroviral targets to develop new drug combinations. The critical role played by the nucleocapsid protein NCp7 of HIV-1 at different steps of the retrovirus life cycle makes it an attractive target for the development of new antiviral agents. NCp7 contains two highly conserved zinc fingers and is characterized by a three-dimensional structure that cannot be modified without a complete loss of infectivity of mutated viruses. Based on these structural data, we report that RB 2121, a cyclic peptide designed to mimic several essential biological determinants of NCp7, displays antiviral activity by inhibiting HIV-1 replication in CEM-4 cells infected by HIV-1. In vitro, RB 2121 does not interfere with HIV-1 cell entry and viral enzymes but is able to inhibit the annealing activities of NCp7 by recognizing nucleic acids. Analysis of proviral DNA synthesis by means of PCR has shown that RB 2121 acts at an early step of the retrovirus life cycle by inducing a dose-dependent reduction in transcribed DNA levels through inhibition of NCp7-reverse transcriptase interaction. Because of its original mechanism of action, RB 2121 provides an interesting lead for the rational development of new anti-HIV-1 agents that could be associated advantageously with enzyme inhibitors to counteract rapid virus mutations and resistance problems observed in tritherapies.

Major versus minor groove DNA binding of a bisarginylporphyrin hybrid molecule: a molecular mechanics investigation.

On the basis of theoretical computations, we have recently synthesised [Perrée-Fauvet, M. and Gresh, N., Tetrahedron Lett., 36 (1995) 4227] a bisarginyl conjugate of a tricationic porphyrin (BAP), designed to target, in the major groove of DNA, the d(GGC GCC)2 sequence which is part of the primary binding site of the HIV-1 retrovirus site [Wain-Hobson, S. et al., Cell, 40 (1985) 9]. In the theoretical model, the chromophore intercalates at the central d(CpG)2 step and each of the arginyl arms targets O6/N7 belonging to guanine bases flanking the intercalation site. Recent IR and UV-visible spectroscopic studies have confirmed the essential features of these theoretical predictions [Mohammadi, S. et al., Biochemistry, 37 (1998) 6165]. In the present study, we compare the energies of competing intercalation modes of BAP to several double-stranded oligonucleotides, according to whether one, two or three N-methylpyridinium rings project into the major groove. Correspondingly, three minor groove binding modes were considered, the arginyl arms now targeting N3, O2 sites belonging to the purine or pyrimidine bases flanking the intercalation site. This investigation has shown that: (i) in both the major and minor grooves, the best-bound complexes have the three N-methylpyridinium rings in the groove opposite to that of the phenyl group bearing the arginyl arms; (ii) major groove binding is preferred over minor groove binding by a significant energy (29 kcal/mol); and (iii) the best-bound sequence in the major groove is d(GGC GCC)2 with two successive guanines upstream from the intercalation. On the other hand, due to the flexibility of the arginyl arms, other GC-rich sequences have close binding energies, two of them being less stable than it by less than 8 kcal/mol. These results serve as the basis for the design of derivatives of BAP with enhanced sequence selectivities in the major groove.

Joint molecular modeling and spectroscopic studies of DNA complexes of a bis(arginyl) conjugate of a tricationic porphyrin designed to target the major groove.

To target selectively the major groove of double-stranded B DNA, we have designed and synthesized a bis(arginyl) conjugate of a tricationic porphyrin (BAP). Its binding energies with a series of double-stranded dodecanucleotides, having in common a central d(CpG)2 intercalation site were compared. The theoretical results indicated a significant energy preference favoring major groove over minor groove binding and a preferential binding to a sequence encompassing the palindrome GGCGCC encountered in the Primary Binding Site of the HIV-1 retrovirus. Spectroscopic studies were carried out on the complexes of BAP with poly(dG-dC) and poly(dA-dT) and a series of oligonucleotide duplexes having either a GGCGCC, CCCGGG, or TACGTA sequence. The results of UV-visible and circular dichroism spectroscopies indicated that intercalation of the porphyrin takes place in poly(dG-dC) and all the oligonucleotides. Thermal denaturation studies showed that BAP increased significantly the melting temperature of the oligonucleotides having the GGCGCC sequence, whereas it produced only a negligible stabilization of sequences having CCCGGG or TACGTA in place of GGCGCC. This indicates a preferential binding of BAP to GGCGCC, fully consistent with the theoretical predictions. IR spectroscopy on d(GGCGCC)2 indicated that the guanine absorption bands, C6=O6 and N7-C8-H, were shifted by the binding of BAP, indicative of the interactions of the arginine arms in the major groove. Thus, the de novo designed compound BAP constitutes one of the very rare intercalators which, similar to the antitumor drugs mitoxantrone and ditercalinium, binds DNA in the major groove rather than in the minor groove.

Modeling of inhibitor-metalloenzyme interactions and selectivity using molecular mechanics grounded in quantum chemistry.

We investigated the binding properties of the metalloprotease inhibitors hydroxamate, methanethiolate, and methylphosphoramidate to a model coordination site occurring in several Zn2+ metalloproteases, including thermolysin. This was carried out using both the SIBFA (sum of interactions between fragments ab initio-computed) molecular mechanics and the SCF/MP2 procedures for the purpose of evaluating SIBFA as a metalloenzyme modeling tool. The energy-minimized structures were closely similar to the X-ray crystallographic structures of related thermolysin-inhibitor complexes. We found that selectivity between alternative geometries and between inhibitors usually stemmed from multiple interaction components included in SIBFA. The binding strength sequence is hydroxamate > methanethiolate > or = methylphosphoramidate from multiple interaction components included in SIBFA. The trends in interaction energy components, rankings, and preferences for mono- or bidentate binding were consistent in both computational procedures. We also compared the Zn2+ vs. Mg2+ selectivities in several other polycoordinated sites having various "hard" and "soft" qualities. This included a hexahydrate, a model representing Mg2+/Ca2+ binding sites, a chlorophyll-like structure, and a zinc finger model. The latter three favor Zn2+ over Mg2+ by a greater degree than the hydrated state, but the selectivity varies widely according to the ligand "softness." SIBFA was able to match the ab initio binding energies by < 2%, with the SIBFA terms representing dispersion and charge-transfer contributing the most to Zn2+/Mg2+ selectivity. These results showed this procedure to be a very capable modeling tool for metalloenzyme problems, in this case giving valuable information about details and limitations of "hard" and "soft" selectivity trends.

Can a polyproline II helical motif be used in the context of sequence-selective major groove recognition of B-DNA? A molecular modelling investigation.

Proline-rich peptides are known to adopt preferentially the extended polyproline II (PPII) helical conformation, which is involved in several protein-protein recognition events. By resorting to molecular modelling techniques, we wished to investigate the extent to which PPII helices could be used for the formation of isochelical peptide-DNA complexes leading to the selective recognition of the major groove of B-DNA. For that purpose, we have grafted to a cationic intercalator, 9-amino-acridine, an oligopeptide having the sequence: Pro- Arg-Pro-Pro-Arg-Pro-Pro-Arg-Pro-Pro-Asp-Pro-Pro. Each residue in the sequence was set in the D configuration, to prevent enzymatic hydrolysis, and each Arg residue was designed to target O6/N7 of a guanine base following the intercalation site. The Asp residue was designed to target a cytosine base, whilst simultaneously forming a bidentate complex with the Arg three residues upstream. Energy-minimization, using the JUMNA procedure, led to the following conclusions : 1) major groove binding is favoured over minor groove or exclusive binding to the phosphates by large energy differences, of over 50 and 90 kcal/mole, respectively: 2) the two best bound sequences are those having three successive guanine bases on the same DNA strand, immediately adjacent to the intercalation site. Sequence d(CGGGC G), encountered in the Primer Binding Site of the HIV retrovirus, thus ranks amongst the best-bound sequences; 3) replacement of an individual guanine amongst the three ones upstream of the intercalation site, by an adenine base, weakens by > 6 kcal/mole the binding energetics; 4) the conformational rigidity of the DNA-bound PPII helix should enable for a modulation of the base sequence selectivity, by appropriate replacements of the Arg and Asp residues. Thus sequence CGGCAAG, also encountered in the HIV genome, could be targeted by an oligopeptide having the sequence Pro-Arg-Pro-Pro-Asp-Pro-Pro-Asn-Pro-Pro-Asn-Pro-Pro-Arg-Ala.

Theoretical design, chemical synthesis and footprinting analysis of a novel peptide derivative of the intercalator 7-H pyridocarbazole targeted towards the major groove of DNA.

In order to target the major groove of DNA, we have designed novel peptide derivatives of 7-H pyridocarbazole, which is the chromophoric ring of ditercalinium, a potent antitumor bisintercalator. We will present here the results obtained with a compound that has a D-Asn tethered to the pyridinium nitrogen of the ring by a protonated beta-alanyl-ethyl chain. We have investigated two alternative means of intercalation of the chromophore: first, into the (pyrpur) sequences, d(CpG)2 and d(CpA).d(TpG); second, into the (pur-pyr) sequences, d(GpC)2 and d(GpT).d(ApC). For the first intercalative mode, the best bound triplet sequences are d(ACG).d(CGT) and d(ACA)d(TGT), namely with an adenine immediately upstream from the intercalation site. In these complexes, the chromophore has its concave side in the major groove, its long axis nearly colinear with the mean long axis of the two base pairs of the intercalation site, and a bidentate H-bonded configuration occurs which involves the C = O and NH groups of the D-Asn side chain and HN6 and N7 (resp.) of the adenine base upstream. One alkylammonium proton is H-bonded to N7 of the guanine of the intercalation site, on the strand opposite to the one bearing the adenine. In the second intercalative mode, the chromophore's concave site now faces one DNA strand, and both alkylammonium protons are involved in H-bonds with N7 and O6 of the 3' guanine on the same strand. The peptide's complexes with sequences having A, G, or C upstream of this guanine were computed to be energetically competitive with those with the best (pyr-pur) triplets. This provides a rare example of energetically favourable drug intercalation in-between (pur-pyr) sequences as compared to the standard (pyr-pur) ones. The synthesis of this compound was performed, and a series of footprinting experiments undertaken on a total of approximately 300 nucleotides. These experiments were consistent with the inferences from the theoretical computations.

Structure and energetics in the complexes of a double-stranded B-DNA dodecamer with netropsin derivatives of a tricationic water-soluble porphyrin: a theoretical investigation.

The structural and energetical characteristics of the complexes formed between two auto-complementary DNA dodecamers, d(CGCGAATTCGCG)2, and d(GCGCAATTGCGC)2, and two novel netropsin (I) and glycine-netropsin (II) conjugates of a tricationic water-soluble porphyrin are investigated in detail by means of theoretical computations. This study was prompted by the successful chemical synthesis of II, which was recently reported (Anneheim-Herbelin, G., Perrée-Fauvet, M., Gaudemer, A., Hélissey, P., Giorgi-Renault, S. and Gresh, N., Tetrahedron Lett. 34, 7263 (1993)). The results indicate that: a) Intercalative binding of II does not entail significant distortions of the DNA backbone, and the Net moiety can bind tightly to the core of the minor groove. b) Intercalative binding of I is computed to energetically weaker than that of II. This is a consequence of the reduced length of the oligopeptide arm, such that the terminal propionamidinium group interacts less favorably with the fourth A-T base-pair than is the case with II. c) Nonintercalative binding of II produces considerable conformational distortions of the DNA. These results in a break of the DNA axis in between the fourth and the fifth base-pairs, namely, at the level where the long axis of the chromophore and of the oligopeptide intersect.

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.

A supermolecule study of the effect of hydration on the conformational behaviour of leucine-enkephalin.

A theoretical conformational study was performed on leu-enkephalin in its zwitterionic form, both in vacuo and in the presence of a number, n, of up to 13 water molecules saturating its first hydration shell. The intramolecular energy of enkephalin as well as the intermolecular enkephalin-water and water-water interaction energies were computed with the SIBFA procedure (Sum of Interactions Between Fragments Ab initio computed), which uses additive ab initio multipole systematics and analytical formulas grounded on ab initio SCF computations. Energy minimizations were performed with a polyvalent minimizer, Merlin, with which three distinct derivative and three distinct nonderivative minimizers can be activated in a sequential fashion. Eight different candidate conformations of enkephalin were used as starting points. These conformations are either those found in distinct X-ray structures, or those proposed on the basis of theoretical computations by other authors. In the absence of hydration, they converged towards distinct folded energy-minima, the best four ones being separated by an energy gap of 8.7 kcal/mol. In marked contrast, with up to n = 13, the energetical separation between the six best conformers narrowed down to congruent to 4 kcal/mol. They can be characterized by: (a) either a direct or a water-mediated ammonium-carboxylate interaction; b) either a close proximity (as in morphine) or a large separation between the aromatic rings of Tyr and Phe (intercenter separations of congruent to 4.5 A and congruent to 10.5 A, respectively), with each of the four mutual combinations of (a) and (b) being represented.

Theoretical design of a bistetrapeptide derivative of mitoxantrone targeted towards the double-stranded hexanucleotide sequence d(GGCGCC)2.

The hexanucleotide d(GGCGCC)2 is encountered in recurrent fashion within transcriptional activating sequences in retroviruses and protooncogenes. Our first theoretical design of novel oligopeptide derivatives of mitoxantrone, MTX (1), had enabled us to predict derivatives depsiGly-Lys(L) and depsiGly-Gly-Orn(D) to preferentially target the tetrameric core d(CCGG)2. Owing to the crucial importance of hexamer d(GGCGCC)2, we have attempted to extend the realm of our approach by now targeting this specific hexanucleotide. For that purpose, we undertook the design of further oligopeptide derivatives of MTX, in which each arm was identically amidated (rather than esterified as in (1)) by tri- or tetrapeptides of varying sequences and individual residue configurations. The binding affinities of these derivatives to the palindromic sequences d(GGCGCC)2, d(CGCGCG)2, d(GCCGGC)2 and d(CCCGGG)2, were compared by energy-minimization. We report here the results obtained with the most promising derivative, having the sequence Arg(L)-Gly-Val(L)-Glu(L), and displaying a considerable energy preference for d(GGCGCC)2 over the other candidate hexameric sites (referred to as I). In the corresponding complexes, the two arms are in two mutually antiparallel directions in the major groove, and adopt a beta-sheet like arrangement stabilized by two H-bonds involving the carbonyl and amide groups of the Gly residues. Each Arg side chain on a given arm chelates O6 and N7 atoms of G1, G2/G1', G2' with its imino and cis amino hydrogen, and is simultaneously bound through two amino hydrogens in a bidentate interaction with the Glu residue.

Theoretical variation of the H alpha chemical shift in acetyl-glycyl-N-methylamide and oligoglycines with molecular conformation and environment.

The sum of the magnetic anisotropy and polarization contributions to the magnetic shielding constants of the alpha protons is calculated as a function of the torsion angles about the NC alpha (phi) and C alpha C' (psi) bonds of the dipeptide. The results show that the polarization or electric field effect is several fold larger than the magnetic anisotropy contribution. The calculated variations are large enough to account for the spread of the values measured for these protons in peptides and proteins. The results obtained for polyglycine alpha helices and antiparallel beta sheets are discussed in relation with molecular conformation and environmental effects on the one hand and experimental data on the other.

Local interactions in peptides. 1H-1H, 13C-H coupling constants and calculations for the conformational analysis of N-acetyl-N'-methylamides of aliphatic amino acids.

We report the results of a joint NMR and theoretical investigation devoted to the conformational properties of N-acetyl-N'-methylamides of aliphatic amino acids with side chains of increasing bulkiness: Gly, Ala, Leu, Ile, and tert.Leu. In this series, determination of the coupling constants 3JHNC alpha H together with the coupling constants 3JC'NC alpha H (thanks to specific carbon-13 labeling of the N-acetyl carbonyl group) led to the derivation of alternative A, B, and C parameters in a Karplus-type relation expressing the dependence of 3JC'NC alpha H upon the phi dihedral angle. The value of the latter is found to increase regularly following the increase of the side-chain bulkiness. The theoretical conformational analysis is performed by applying the SIBFA procedure, which uses empirical formulas based on ab initio SCF computations. The conformational energy maps illustrate the progressive distortion of the backbone conformation incurred in the series Gly to tert.Leu. Theoretical values computed for 3JHNC alpha H and 3JC'NC alpha H are found to be in a good quantitative agreement with the experimental ones.

A joint 2D NMR and theoretical investigation of Ca2+ binding loops III and IV of calmodulin.

A 2D NMR NOESY spectrum of integral CaM in water(148 residues) reveals a series of downfield-shifted crosspeaks stemming from the NH protons of the Ca2(+)-binding loops III and IV. Their attribution, with the help of already assigned proton resonances of isolated tryptic fragments, was complemented by means of energy-minimizations on the Ca2+ complexes of loops III and IV. From these calculations, a set of two alternative, related conformations was obtained for each loop. The first type of conformation provides a coordination pattern for Ca2+ that is similar to that found in loop EF of parvalbumin. The computed interproton distances in both loops are fully compatible with the inferences from the sets of NOESY cross-peaks. Evidence is also provided for interloop interactions.

Theoretical design of novel, 4 base pair selective derivatives of mitoxantrone.

Mitoxantrone (MTX) is a recently synthesized antitumor intercalative molecule, currently in use in chemotherapy. Previous theoretical computations showed that the base pair selectivity of MTX is limited to the sole two base-pair sequence making up the intercalation site. In order to further extend the recognition site, we undertook, by means of theoretical computations, the design of novel MTX derivatives, in which the terminal hydroxyl group of each side chain is esterified with oligopeptides. We compare in the present study the binding affinities of two derivatives, depsiGly-Lys(D) and depsiGly-Gly-Orn(L), for the palindromic sequences d(CCCGGG)2, d(GCCGGC)2, d(GGCGCC)2, and d(CGCGCG)2. Major groove binding of the oligopeptide arms was shown to be significantly more favourable than either minor groove binding, or binding to the sole phosphate groups. With the two arms adopting two antiparallel directions, two distinct arrangements were investigated in the major groove: (a) the two oligopeptides are brought closer together by means of two hydrogen bonds involving the backbone of their second residue in a beta-sheet like arrangement; (b) the two arms are remote from each other so as to reduce their mutual electrostatic repulsion. Whatever the disposition, the optimal binding configurations were invariably found to be those in which the cationic side chains of the terminal residues chelate N7/O6 of two successive guanines, whenever present on a given strand. A distinct energetical preference for arrangement (a) was obtained with the depsiGly-Gly-Orn(L) derivative. Replacement of the central Gly residue by a Cys one, as in the sequence depsiGly-Cys-Orn(L), was proposed subsequently, so as to further stabilize such a beta-sheet arrangement by means of a disulfide bridge between the two Cys residues. The two investigated compounds were shown to preferentially bind sequences d(CCCGGG)2 and d(GCCGGC)2, with a tetrameric core CCGG rather than sequences d(GGCGCC)2 and d(CGCGCG)2, with a tetrameric core GCGC.

Modelling of the binding specificity in the interactions of cationic porphyrins with DNA.

A theoretical investigation is performed of the complexes of a tetracationic porphyrin, tetra-(4-N-methylpyridyl)-porphyrin, (T4MPyP), with the hexanucleotides d(CGCGCG)2 and d(TATATA)2, considering the possibility of both the intercalative and the groove binding interactions. These computations demonstrate that T4MPyP manifests a significant preference for intercalation in its complex with d(CGCGCG)2 but for non intercalative binding in the minor groove in its complex with d(TATATA)2. Such a dual binding behaviour of T4MPyP as a function of the sequence to which it is attached is fully consistent with available experimental data. It demonstrates that intercalation and groove binding may be viewed as two potential wells on a continuous energy surface. In agreement with experiment, the computations indicate that in the here considered case the deepest well is associated with intercalation.

A theoretical investigation of the base sequence preferences of monointercalating polymethylene carboxamide derivatives 9-aminoacridine.

Theoretical computations are performed of the comparative binding affinities of five polymethylene carboxamide derivatives of 9-aminoacridine to a series of double-stranded hexanucleotides. The purpose of this investigation is to ascertain whether minor groove recognition of a guanine base adjacent to the intercalation site can occur, and be preferentially stabilized, for a given length of the polymethylene side chain, encompassing from n = 2 up to n = 6 methylene groups. For that purpose, several representative sequences were investigated, in which intercalation of the 9-aminoacridine chromophore occurred at a central d(CpG) or d(TpA) step. Investigated were the self-complementary sequences d(CGCGCG)2, d(GCCGGC)2, d(TATATA)2 and d(ATTAAT)2, as well as the 'mixed' sequences d(ACTAAT) .d(ATTAGT) and d(TGTATA). d(TATACA). For n = 3 up to n = 6, such a recognition was enabled only when the guanine base was located downstream of the intercalation site, i.e. with steps d(CGG) and d(TAG). It occurred by means of a bidentate interaction involving, on the one hand, H(N2) and N3 of the base, and, on the other hand, the carbonyl oxygen and the cis amino hydrogen of the terminal formamide moiety of the ligand. Because of the flexibility of the side chain, however, alternative binding modes were also found to occur competitively, involving backbone-only interactions of the side chain. On the basis of the present computations, upon binding to the sequence d(GCCGGC)2, an optimal value of n = 5 could be derived, with the corresponding acridine derivative eliciting both a significant prevalence of the bidentate over backbone only binding mode, and the most favourable energy balance within the investigated series. This privileged value of n = 5 is fully consistent with the experimental results of Markovits et al. and Gaugain et al. The very flexibility of the side chain, however, hampered any preferential recognition of a triplet sequence with a downstream guanine, such as d(CGG) or d(TAG), to be elicited over sequences such as d(TAA), d(TAT) or d(TAC).