A tentative model of the intercalative binding of the neocarzinostatin chromophore to double-stranded tetranucleotides.

Theoretical computations are performed of the intercalative binding of the neocarzinostatin chromophore (NCS) with the double-stranded oligonucleotides d(CGCG)2, d(GCGC)2, d(TATA)2 and d(ATAT)2. Minor groove binding is preferred over major groove binding. It is found that the long axis of the stacked naphtoate ring lies approximately parallel to the long axis of the base pairs of the intercalation site. The galactosamine ammonium group interacts with specific sites of the groove (O2/N3 of bases 2 and O1' of sugar S3), whereas the dodecadyine ring system wraps around the groove towards the backbone. An overall AT versus GC preference is derived. Intercalation in a central purine-(3', 5')-pyrimidine sequence appears to be preferred over that in a central pyrimidine-(3', 5')-purine sequence.     

Modelling basic features of specificity in the binding of a dicationic steroid diamine to double-stranded oligonucleotides.

An investigation of the intrinsically preferred binding modes of a steroid diamine, dipyrandium, to the double-stranded hexanucleotides d(TATATA)2, d(ATATAT)2, and d(CGCGCG)2 is carried out by the energy minimization procedure JUMNA. Several alternative binding modes are compared: groove binding in which the conformation of the oligonucleotide remains close to that of B-DNA, intercalation between base-pairs and interaction with variously kinked structures in which base pairs of dinucleoside steps open towards the groove in which the binding occurs. The favored binding configuration occurs at the d(TpA) step of the AT kinked nucleotides in which the kink opens the base pairs towards the minor groove. Thus, for the d(T1A2T3A4T5A6)2 sequences the preferred complexation involves the kink at the T3A4 step facing the cyclohexane rings A, B, and C of the ligand. For the d(A1T2A3T4A5T6)2 sequence, the kink occurs at the T2A3 step facing the cationic pyrrolidine ring linked to ring A. The binding of dipyrandium to d(CGCGCG)2 is found to be considerably less favourable than for either of the two (AT) sequences.     

Drug recognition of DNA. Proposal for GC minor groove specific ligands: vinylexins

In a previous publication in this journal we have proposed an isolexin-like prototype of a GC minor groove specific ligand. The present paper is devoted to refinements of this prototype (increase in specificity and in DNA binding energy). It is shown that only a very limited improvement can be obtained by increasing the proton accepting capabilities of the heteroaromatic ring systems of the prototype, although these rings interact directly with the proton donating NH2 group of guanine. On the other hand a significant increase both in GC specificity and in DNA binding energy is obtained by replacing the NH linkers of the isolexin by C = C double bonds (yielding what we term "vinylexins"). Specificity is still largely conserved and the DNA binding energy is significantly increased in monocationic vinylexins, which should thus be efficient GC minor groove specific ligands. The outstanding importance for the GC specificity of the C = C linkers is evidenced by the disappearance of this specificity when these linkers are replaced by peptide bonds (peptilexins). On the other hand vinylexins with proton donating heteroaromatic rings are, as expected, AT specific. The vinylexin family may thus represent universal minor groove binding agents susceptible to bind to any given base pair sequence of DNA, following the positioning of their proton donor and proton acceptor rings. This study confirms the insufficiency of purely geometrical and/or hydrogen bonding considerations for the correct estimation of GC versus AT specificity of groove binding ligands. These can only be accounted for by taking into consideration the overall electronic properties of the interacting species and explicitly calculating the energies of complex formation including all the relevant contributions.     

The effect of the amino-acid side chains on the energy profiles for ion transport in the gramicidin A channel

Computations on the energy profiles for Na+ in the gramicidin A (GA) channel have been extended by introducing the effect, previously neglected, of the amino acid side chains of GA, fixed in their most stable conformations. The calculations have been performed in two approximations: 1) with the ethanolamine tail fixed in its most stable conformation, 2) with the tail allowed to optimize its conformation upon the progression of the ion. In both approximations the overall shape of the energy profile is very similar to that obtained in the absence of the side chains. One observes, however, a general lowering of the profile upon the adjunction of the side chains. The analysis of the factors responsible for this energy lowering indicates that it is due essentially to the electrostatic and polarisation components of the interaction which interplay differently, however, in the different parts of the channel. A particular role is attributed in this respect to the tryptophan residues of GA. The role of the 4 tryptophans present, Trp 15, 13, 11 and 9, is individualized by stripping of one of them at a time. The strongest effect on the energy deepening is due to Trp 13 and is particularly prominent in the entrance zone at 14.5A from the center of the channel. The result indicates the possibility of investigating theoretically the effect on the energy profiles of the substitution of the "natural" side chain by others.     

Theoretical exploration of netropsin binding to tRNA(Phe)

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

The dependence of the surface electrostatic potential of B-DNA on environmental factors

The electrostatic potential of B-DNA is calculated on its surface envelope for two homopolymeric base pair sequences using models representing the effects of both counterion binding and of aqueous solution. The influence of these two factors on the resulting potentials is established and the significance of calculations which omit such effects is discussed.     

A theoretical study of anthracene and phenanthrene derivatives acting as A-T specific intercalators.

Theoretical computations are performed on the comparative A-T versus G-C binding selectivities of two DNA intercalating molecules recently synthesized by Wilson et al. These are derivatives of phenanthrene and anthracene with side chains containing an hydroxy group bound to its C alpha carbon and a cationic amino group bound to its C beta carbon. We have optimized the binding energies of these phenanthrene and anthracene derivatives (1 and 2, respectively) to the double-stranded tetramers d(ATAT)2 and d(GCGC)2, the intercalation occurring in the central pyrimidine (3'-5') purine sequence. The sum of the intercalator-oligonucleotide intermolecular interaction energy plus the conformational energy variation of the intercalator upon binding were computed by the SIBFA procedures, which use empirical formulas based on ab initio SCF computations. Both compounds are found to bind more favourably to the AT sequence than to the GC one. Moreover, the affinity of 1 for the AT oligomer is computed to be larger than that of 2, whereas conversely that of 2 is larger than that of 1 for the GC oligomer. The AT versus GC binding selectivity of 1 is significantly larger than that of 2. These results are in excellent agreement with the experimental findings of Wilson et al. However, contrary to the suggestion of these authors the alpha-hydroxy group of the side chain of the intercalators does not seem to play a decisive role in determining the A-T specificity.     

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.     

Molecular orbital calculations on the conformation of polypeptides and proteins. V. Conformational energy maps and stereochemical rotational states of aliphatic residues

The quantum-mechanical calculations by the PCILO method on the conformation of amino acid residues of proteins have been extended to the valyl, leucyl, and isoleucyl residues. In distinction to the earlier “empirical” computations, the quantum-mechanical results indicate very similar energy contours for the stable conformations of the three residues. Their general outline is also similar to that of the alanyl residue, although reduced by about 25%. Contrary to the “empirical” computations, the present results predict that the region corresponding to the α-helix should be one of great stability for the three residues and in particular for the valyl residue. The quantum-mechanical results are in excellent agreement with the experimental conformations of the aliphatic residues in lysozyme and myoglobin. Their prediction as to the ready availability of the valyl residue in the α-helical conformation agrees moreover with Ptitsyn's statistical evaluation of the participation of this residue in the inner turns of the helical regions in six globular proteins. The maximum conformational space allowed for the aliphatic residues is somewhat smaller than that allowed for the aromatic ones, while the minimum conformational space (region of stability common to all the residues) is similar in both groups.