Interaction of unfused tricyclic aromatic cations with DNA: a new class of intercalators

Unfused tricyclic aromatic ring systems 1-6 with one or two cationic side chains have been synthesized and their interactions with DNA and synthetic polymers probed with a variety of techniques. Molecular mechanics calculations indicate that the torsional angle between ring planes in the minimum energy conformation of the tricyclic molecules can range from 0 degree to as high as 50 degrees depending on the type of rings and substituents. Viscometric titrations with linear and supercoiled DNA, linear dichroism, and NMR studies indicated that all compounds with torsional angles of approximately 20 degrees or less bind to DNA by intercalation. The more highly twisted intercalators caused significant perturbation of DNA structure. Unfused intercalators with twist angles of approximately 20 degrees have reduced binding constants, suggesting that they could not form an optimum interaction with the DNA base pairs. Unfused intercalators with twist less than 20 degrees formed strong complexes with DNA. The structures of these unfused intercalators are more analogous to typical groove-binding molecules, and an analysis of their interaction with DNA provides a better understanding of the subtle differences between intercalation and groove-binding modes for aromatic cations. The results indicate that intercalation and groove-binding modes should be viewed as two potential wells on a continuous energy surface. The results also suggest design strategies for intercalators that can optimally complement DNA base pair propeller twist or that can induce bends in DNA at the intercalation site.     

A molecular mechanical model to predict the helix twist angles of B-DNA

We present here a model for the prediction of helix twist angles in B-DNA, a model composed of a collection of torsional springs. Statistically averaged conformational energy calculations show that, for a specified basepair step, the basepair-basepair conformational energy is quadratically dependent on the helix twist angle, so the calculations provide the spring parameters for the basepair-basepair interactions. Torsional springs can also be used to model the effects of the backbone on the helix twist, and the parameters for those springs are derived by fitting the model to experimental data. The model predicts a macroscopic torsional stiffness and a longitudinal compressibility (Young's modulus) which are both in good agreement with experiment. One biological consequence of the model is examined, the sequence specificity of the Eco RI restriction endonuclease, and it is shown that the discriminatory power of the enzyme receives a substantial contribution from the energetic cost of torsional deformations of the DNA when wrong sequences are forced into the enzyme binding site.     

Structure of DNA sequence d-TGATCA by two-dimensional nuclear magnetic resonance spectroscopy and restrained molecular dynamics

The 5' d-TpG 3' element is a part of DNA sequences involved in regulation of gene expression and is also a site for intercalation of several anticancer drugs. Solution conformation of DNA duplex d-TGATCA containing this element has been investigated by two-dimensional NMR spectroscopy. Using a total of 12 torsional angles and 121 distance constraints, structural refinement has been carried out by restrained molecular dynamics (rMDs) in vacuum up to 100 ps. The structure is characterized by a large positive roll at TpG/CpA base pair step and large negative propeller twist for AT and TA base pairs. The backbone torsional angle, gamma(O5'-C5'-C4'-C3'), of T1 residue adopts a trans-conformation which is corroborated by short intra nucleotide T1H6-T1H5' (3.7A) distance in nuclear overhauser effect spectroscopy (NOESY) spectra while the backbone torsional angle, beta(P-O5'-C5'-C4'), exists in trans as well as gauche state for T1 and C5 residues. There is evidence of significant flexibility of the sugar-phosphate backbone with rapid inter-conversion between two different conformers at TpG/CpA base pair step. The base sequence dependent variations and local structural heterogeneity have important implications in specific recognition of DNA by ligands.     

Torsional constant of 27-mer DNA oligomers of different sequences.

We have studied the torsional elastic constant (alpha) of short DNA (27mer) oligomers of various sequence by fluorescence polarization anysotropy (FPA) measurements. The lowest alpha values were found in samples with sequence rich in AA dinucleotides or containing the alternating d(A-T) x d(A-T) motif. The torsional rigidity of our DNA samples was compared to that calculated according to the current values of twist angle fluctuations derived for ten dinucleotide steps by recent analyses of DNA crystal structure database. The values of torsional rigidity derived from crystals are higher than our experimental ones, obtained by FPA analysis, suggesting that packing force in crystals may notably hinder the dinucleotide twist angle fluctuations that occur in solution. This behaviour is more evident for samples containing AA, TA and AT steps. In all the samples there is about a twofold change of the alpha value in the 10-40 degrees C range. An activation enthalpy (Delta H (#)) of about 17.4 kJ mol(-1), on average, was obtained for the temperature dependence of eight of the ten samples studied. A correlation with the stacking energy is discussed.     

Structure of DNA hexamer sequence d-CGATCG by two-dimensional nuclear magnetic resonance spectroscopy and restrained molecular dynamics

Solution conformation of self-complementary DNA duplex d-CGATCG, containing 5' d-CpG 3' site for intercalation of anticancer drug, daunomycin and adriamycin, has been investigated by nuclear magnetic resonance (NMR) spectroscopy. Complete resonance assignments of all the protons (except some H5'/H5" protons) have been obtained following standard procedures based on double quantum filtered correlation spectroscopy (dQF COSY) and two-dimensional nuclear Overhauser effect (NOE) spectra. Analysis of sums of coupling constants in one-dimensional NMR spectra, cross peak patterns in dQF COSY spectra and inter proton distances shows that the DNA sequence assumes a conformation close to the B-DNA family. The deoxyribose sugar conformation is in dynamic equilibrium with predominantly S-type conformer and a minor N-type conformer with N<-->S equilibrium varying with temperature. At 325 K, the mole fraction of the N-conformer increases for some of the residues by approximately 9%. Using a total of 10 spin-spin coupling constants and 112 NOE intensities, structural refinement has been carried out using Restrained Molecular Dynamics (rMD) with different starting structures, potential functions and rMD protocols. It is observed that pseudorotation phase angle of deoxyribose sugar for A3 and T4 residues is approximately 180 degrees and approximately 120 degrees, respectively while all other residues are close to C2'endo-conformation. A large propeller twist (approximately -18 degrees) and smallest twist angle (approximately 31 degrees) at A3pT4 step, in the middle of the sequence, a wider (12 A) and shallower (3.0 A) major groove with glycosidic bond rotation as high anti at both the ends of hexanucleotide are observed. The structure shows base-sequence dependent variations and hence strong local structural heterogeneity, which may have implications in ligand binding.     

A zinc-binding region in Vif binds Cul5 and determines cullin selection

Human immunodeficiency virus-1 (HIV-1) Vif overcomes the anti-viral activity of APOBEC3G by targeting it for ubiquitination via a Cullin 5-ElonginB-ElonginC (Cul5-EloBC) E3 ligase. Vif associates with Cul5-EloBC through a BC-box motif that binds EloC, but the mechanism by which Vif selectively recruits Cul5 is poorly understood. Here we report that a region of Vif (residues 100-142) upstream of the BC-box binds selectively to Cul5 in the absence of EloC. This region contains a zinc coordination site HX5CX17-18CX3-5H (HCCH), with His/Cys residues at positions 108, 114, 133, and 139 coordinating one zinc ion. The HCCH zinc coordination site, which is conserved among primate lentivirus Vif proteins, does not correspond to any known class of zinc-binding motif. Mutations of His/Cys residues in the HCCH motif impair zinc coordination, Cul5 binding, and APOBEC3G degradation. Mutations of conserved hydrophobic residues (Ile-120, Ala-123, and Leu-124) located between the two Cys residues in the HCCH motif disrupt binding of the zinc-coordinating region to Cul5 and inhibit APOBEC3G degradation. The Vif binding site maps to the first cullin repeat in the N terminus of Cul5. These data suggest that the zinc-binding region in Vif is a novel cullin interaction domain that mediates selective binding to Cul5. We propose that the HCCH zinc-binding motif facilitates Vif-Cul5 binding by playing a structural role in positioning hydrophobic residues for direct contact with Cul5.     

Solution structure of an RNA fragment with the P7/P9.0 region and the 3'-terminal guanosine of the tetrahymena group I intron

In the second step of the two consecutive transesterifications of the self-splicing reaction of the group I intron, the conserved guanosine at the 3' terminus of the intron (omegaG) binds to the guanosine-binding site (GBS) in the intron. In the present study, we designed a 22-nt model RNA (GBS/omegaG) including the GBS and omegaG from the Tetrahymena group I intron, and determined the solution structure by NMR methods. In this structure, omegaG is recognized by the formation of a base triple with the G264 x C311 base pair, and this recognition is stabilized by the stacking interaction between omegaG and C262. The bulged structure at A263 causes a large helical twist angle (40 +/- 80) between the G264 x C311 and C262 x G312 base pairs. We named this type of binding pocket with a bulge and a large twist, formed on the major groove, a "Bulge-and-Twist" (BT) pocket. With another twist angle between the C262 x G312 and G413 x C313 base pairs (45 +/- 100), the axis of GBS/omegaG is kinked at the GBS region. This kinked axis superimposes well on that of the corresponding region in the structure model built on a 5.0 A resolution electron density map (Golden et al., Science, 1998, 282:345-358). This compact structure of the GBS is also consistent with previous biochemical studies on group I introns. The BT pockets are also found in the arginine-binding site of the HIV-TAR RNA, and within the 16S rRNA and the 23S rRNA.     

F-actin retains a memory of angular order

Modifications can be made to F-actin that do not interfere with the binding of myosin but inhibit force generation, suggesting that actin's internal dynamics are important for muscle contraction. Observations from electron microscopy and x-ray diffraction have shown that subunits in F-actin have a relatively fixed axial rise but a variable twist. One possible explanation for this is that the actin subunits randomly exist in different discrete states of "twist, " with a significant energy barrier separating these states. This would result in very slow torsional transitions. Paracrystals impose increased order on F-actin filaments by reducing the variability in twist. By looking at filaments that have recently been dissociated from paracrystals, we find that F-actin retains a "memory" of its previous environment that persists for many seconds. This would be consistent with slow torsional transitions between discrete states of twist.     

Editing and diagonal peak suppression in three-dimensional HCCH protein NMR correlation experiments

A novel three-dimensional (3D) HCCH NMR experiment is introduced. It involves ¹³C-¹³C COSY or TOCSY coherence transfer plus two independent editing steps according to the number of protons attached to the individual carbons before and after the ¹³C-¹³C homonuclear mixing. This double editing leads to simplification of HCCH protein side chain spectra that otherwise are prone to spectral overlap. Another interesting feature is amino acid selectivity, i.e. that the presence of certain correlations in a doubly edited HCCH subspectrum gives a clue as to assignment to a particular subgroup of amino acids or segments thereof. Finally, the selection of two different multiplicities in the two editing steps leads to diagonal peak suppression in the ¹H-¹H (3D spectrum recorded with two ¹H and one ¹³C dimension) or the ¹³C-¹³C (3D spectrum recorded with one ¹H and two ¹³C dimensions) two-dimensional projection. The new experiment is demonstrated using a ¹³C,¹⁵N-labeled protein sample, chymotrypsin inhibitor 2, at 500 MHz.     

1H NMR study of self-association and restricted internal rotation of the C8-substituted deoxyguanosine 5'-monophosphate adduct of the carcinogen 2-(acetylamino)fluorene

The high-field 1H NMR spectra of a nucleotide-carcinogen adduct formed from 2-(acetylamino)fluorene (8-(N-fluoren-2-ylacetamido)-2'-deoxyguanosine 5'-monophosphate) have been examined in aqueous solution as a function of concentration at high and low temperatures. An anomalous concentration dependence of NMR spectra was observed at concentration levels over 1 mM. These spectral characteristics have been analyzed in terms of changes in self-association and in the interconversions between torsional diastereomers associated with the central nitrogen. Association constants have been computed. Stacking interactions, which involve both the fluorene and guanine rings, are strong, cooperative and highly temperature-dependent. Deacetylation alters the mode of stacking. Several effects of solvent and aggregation on the conformation at the central nitrogen are discussed.     

Structural analysis and dynamics of retinal chromophore in dark and meta I states of rhodopsin from 2H NMR of aligned membranes

Rhodopsin is a prototype for G protein-coupled receptors (GPCRs) that are implicated in many biological responses in humans. A site-directed (2)H NMR approach was used for structural analysis of retinal within its binding cavity in the dark and pre-activated meta I states. Retinal was labeled with (2)H at the C5, C9, or C13 methyl groups by total synthesis, and was used to regenerate the opsin apoprotein. Solid-state (2)H NMR spectra were acquired for aligned membranes in the low-temperature lipid gel phase versus the tilt angle to the magnetic field. Data reduction assumed a static uniaxial distribution, and gave the retinylidene methyl bond orientations plus the alignment disorder (mosaic spread). The dark-state (2)H NMR structure of 11-cis-retinal shows torsional twisting of the polyene chain and the beta-ionone ring. The ligand undergoes restricted motion, as evinced by order parameters of approximately 0.9 for the spinning C-C(2)H(3) groups, with off-axial fluctuations of approximately 15 degrees . Retinal is accommodated within the rhodopsin binding pocket with a negative pre-twist about the C11=C12 double bond that explains its rapid photochemistry and the trajectory of 11-cis to trans isomerization. In the cryo-trapped meta I state, the (2)H NMR structure shows a reduction of the polyene strain, while torsional twisting of the beta-ionone ring is maintained. Distortion of the retinal conformation is interpreted through substituent control of receptor activation. Steric hindrance between trans retinal and Trp265 can trigger formation of the subsequent activated meta II state. Our results are pertinent to quantum and molecular mechanics simulations of ligands bound to GPCRs, and illustrate how (2)H NMR can be applied to study their biological mechanisms of action.     

Assignments of 31P NMR resonances in oligodeoxyribonucleotides: origin of sequence-specific variations in the deoxyribose phosphate backbone conformation and the 31P chemical shifts of double-helical nucleic acids

It is now possible to unambiguously assign all 31P resonances in the 31P NMR spectra of oligonucleotides by either two-dimensional NMR techniques or site-specific 17O labeling of the phosphoryl groups. Assignment of 31P signals in tetradecamer duplexes, (dTGTGAGCGCTCACA)2, (dTAT-GAGCGCTCATA)2, (dTCTGAGCGCTCAGA)2, and (dTGTGTGCGCACACA)2, and the dodecamer duplex d(CGTGAATTCGCG)2 containing one base-pair mismatch, combined with additional assignments in the literature, has allowed an analysis of the origin of the sequence-specific variation in 31P chemical shifts of DNA. The 31P chemical shifts of duplex B-DNA phosphates correlate reasonably well with some aspects of the Dickerson/Calladine sum function for variation in the helical twist of the oligonucleotides. Correlations between experimentally measured P-O and C-O torsional angles and results from molecular mechanics energy minimization calculations show that these results are consistent with the hypothesis that sequence-specific variations in 31P chemical shifts are attributable to sequence-specific changes in the deoxyribose phosphate backbone. The major structural variation responsible for these 31P shift perturbations appears to be P-O and C-O backbone torsional angles which respond to changes in the local helical structure. Furthermore, 31P chemical shifts and JH3'-P coupling constants both indicate that these backbone torsional angle variations are more permissive at the ends of the double helix than in the middle. Thus 31P NMR spectroscopy and molecular mechanics energy minimization calculations appear to be able to support sequence-specific structural variations along the backbone of the DNA in solution.     

Solid-state NMR determination of sugar ring pucker in (13)C-labeled 2'-deoxynucleosides

The H3'-C3'-C4'-H4' torsional angles of two microcrystalline 2'-deoxynucleosides, thymidine and 2'-deoxycytidine.HCl, doubly (13)C-labeled at the C3' and C4' positions of the sugar ring, have been measured by solid-state magic-angle-spinning nuclear magnetic resonance (NMR). A double-quantum heteronuclear local field experiment with frequency-switched Lee-Goldberg homonuclear decoupling was used. The H3'-C3'-C4'-H4' torsional angles were obtained by comparing the experimental curves with numerical simulations, including the two (13)C nuclei, the directly bonded (1)H nuclei, and five remote protons. The H3'-C3'-C4'-H4' angles were converted into sugar pucker angles and compared with crystallographic data. The delta torsional angles determined by solid-state NMR and x-ray crystallography agree within experimental error. Evidence is also obtained that the proton positions may be unreliable in the x-ray structures. This work confirms that double-quantum solid-state NMR is a feasible tool for studying sugar pucker conformations in macromolecular complexes that are unsuitable for solution NMR or crystallography.     

The helical repeat of DNA at high temperature.

The increasing number of studies on thermophilic organisms addressed the question of DNA double helix parameters at high temperature. The present study shows that the helix rotation angle per base pair omega of an unconstrained DNA decreases linearly upon temperature increase, up to the premelting range. In the ionic conditions tested, this rule extends to temperatures up to 85 degrees C, which is a common growth temperature for many hyperthermophilic organisms. In addition, the torsional constant K of DNA decreases with temperature, indicating that the energy required to modify the DNA twist is lower at high temperature. These findings have several implications for people working on the structure and enzymology of DNA at high temperature.     

Backbone dynamics of membrane proteins in lipid bilayers: the effect of two-dimensional array formation as revealed by site-directed solid-state 13C NMR studies on [3-13C]Ala- and [1-13C]Val-labeled bacteriorhodopsin

We have recorded site-directed solid-state 13C NMR spectra of [3-13C]Ala- and [1-13C]Val-labeled bacteriorhodopsin (bR) as a typical membrane protein in lipid bilayers, to examine the effect of formation of two-dimensional (2D) lattice or array of the proteins toward backbone dynamics, to search the optimum condition to be able to record full 13C NMR signals from whole area of proteins. Well-resolved 13C NMR signals were recorded for monomeric [3-13C]Ala-bR in egg phosphatidylcholine (PC) bilayer at ambient temperature, although several 13C NMR signals from the loops and transmembrane alpha-helices were still suppressed. This is because monomeric bR reconstituted into egg PC, dimyristoylphosphatidylcholine (DMPC) or dipalmytoylphosphatidylcholine (DPPC) bilayers undergoes conformational fluctuations with frequency in the order of 10(4)-10(5) Hz at ambient temperature, which is interfered with frequency of magic angle spinning or proton decoupling. It turned out, however, that the 13C NMR signals of purple membrane (PM) were almost fully recovered in gel phase lipids of DMPC or DPPC bilayers at around 0 degrees C. This finding is interpreted in terms of aggregation of bR in DMPC or DPPC bilayers to 2D hexagonal array in the presence of endogenous lipids at low temperature, resulting in favorable backbone dynamics for 13C NMR observation. It is therefore concluded that [3-13C]Ala-bR reconstituted in egg PC, DMPC or DPPC bilayers at ambient temperature, or [3-13C]Ala- and [1-13C]Val-bR at low temperature gave rise to well-resolved 13C NMR signals, although they are not always completely the same as those of 2D hexagonal lattice from PM.     

Sensitivity of NMR internucleotide distances to B-DNA conformation: underlying mechanics.

Nuclear magnetic resonance (NMR) spectroscopy, combining correlated spectroscopy (COSY) coupling constant measurements with nuclear Overhauser effect spectroscopy (NOESY) interatomic distances, should make it possible to determine an averaged solution structure for DNA oligomers. However, even if such data could be obtained with high accuracy, it is not clear which structural parameters of DNA would be determined. Here, the relationships between measurable internucleotide distances and helical parameters are systematically studied through molecular modelling. Investigations are carried out using four representative sequences, (ACGT)n, (TCGA)n, (AGCT)n and (TGCA)n, composed of repeated tetranucleotides belonging to oligomers previously studied by NMR. Correlations between interatomic distances become evident and strong connections between distances and inter-base helical parameters are observed. Results imply that twist, roll, shift and slide values can be accurately determined from NMR data. Sequence independent mechanical coupling which link backbone and sugar conformations to helical twist are also described.     

Regio-selective detection of dynamic structure of transmembrane alpha-helices as revealed from (13)C NMR spectra of [3-13C]Ala-labeled bacteriorhodopsin in the presence of Mn2+ ion.

13C Nuclear magnetic resonance (NMR) spectra of [3-(13)C]Ala-labeled bacteriorhodopsin (bR) were edited to give rise to regio-selective signals from hydrophobic transmembrane alpha-helices by using NMR relaxation reagent, Mn(2+) ion. As a result of selective suppression of (13)C NMR signals from the surfaces in the presence of Mn(2+) ions, several (13)C NMR signals of Ala residues in the transmembrane alpha-helices were identified on the basis of site-directed mutagenesis without overlaps from (13)C NMR signals of residues located near the bilayer surfaces. The upper bound of the interatomic distances between (13)C nucleus in bR and Mn(2+) ions bound to the hydrophilic surface to cause suppressed peaks by the presence of Mn(2+) ion was estimated as 8.7 A to result in the signal broadening to 100 Hz and consistent with the data based on experimental finding. The Ala C(beta) (13)C NMR peaks corresponding to Ala-51, Ala-53, Ala-81, Ala-84, and Ala-215 located around the extracellular half of the proton channel and Ala-184 located at the kink in the helix F were successfully identified on the basis of (13)C NMR spectra of bR in the presence of Mn(2+) ion and site-directed replacement of Ala by Gly or Val. Utilizing these peaks as probes to observe local structure in the transmembrane alpha-helices, dynamic conformation of the extracellular half of bR at ambient temperature was examined, and the local structures of Ala-215 and 184 were compared with those elucidated at low temperature. Conformational changes in the transmembrane alpha-helices induced in D85N and E204Q and its long-range transmission from the proton release site to the site around the Schiff base in E204Q were also examined.     

Temperature dependence of dielectric relaxations in alpha-elastin coacervate: evidence for a peptide librational mode

The dielectric permittivity of alpha-elastin coacervate is reported over the frequency range of 1 MHz to 1000 MHz and the temperature dependence from 6.8 degrees C to 70 degrees C is also reported. A temperature-dependent simple Debye-type relaxation is observed with a correlation time of 8 nsec (40 degrees C) which is similar to that of the polypentapeptide of elastin (i.e. 7 nsec at 40 degrees C) where the band has been assigned to a peptide librational mode. By analogy this allows for the first assignment of a peptide librational mode in a naturally occurring polypeptide or protein. The strong spectrally localized band indicates a regularity of structure. The low temperature dependence of the correlation time, giving a 1.7 kcal/mole enthalpy of activation, is consistent with torsional motions associated with a peptide librational mode.     

NMR spectroscopy of RNA duplexes containing pseudouridine in supercooled water

We have performed NMR experiments in supercooled water in order to decrease the temperature-dependent exchange of protons in RNA duplexes. NMR spectra of aqueous samples of RNA in bundles of narrow capillaries that were acquired at temperatures as low as -18 degrees C reveal resonances of exchangeable protons not seen at higher temperatures. In particular, we detected the imino protons of terminal base pairs and the imino proton of a non-base-paired pseudouridine in a duplex representing the eukaryotic pre-mRNA branch site helix. Analysis of the temperature dependence of chemical shift changes (thermal coefficients) for imino protons corroborated hydrogen bonding patterns observed in the NMR-derived structural model of the branch site helix. The ability to observe non-base-paired imino protons of RNA is of significant value in structure determination of RNA motifs containing loop and bulge regions.     

Conformation and rigidity of DNA microcircles containing waf1 response element for p53 regulatory protein

The tumor-suppressor activity of p53 is closely related to its DNA-binding properties. It binds a number of DNA response-elements and it is likely that these share a common structural feature. Here, we present a new, general method to determine the absolute twist of flexible DNA promoter sequences based on direct imaging of the topology of microcircles containing the sequences. We have used magnetically driven dynamic force microscopy ("MacMode" AFM) to observe, in solution, the conformation of 168 base-pair DNA microcircles, each containing four equally spaced copies of the waf1/cip1/p21 p53 response-element. Analysis of the images showed that the microcircles are markedly puckered with a small excess of negatively writhed molecules. The average measured values of writhe are 0.109+/-0.013 (for 60 positively writhed molecules) and -0.098+/-0.011 (for 65 negatively writhed molecules). These values lead directly to a difference in linking number for the positively and negatively writhed molecules prior to ligation, from which we derive a twist mismatch of 178 degrees (overtwist). This is 44.5 degrees for each 42-mer precursor containing a single waf1/cip1/p21 p53 response-element, in good agreement with the range of values deduced by indirect biochemical techniques. The two values of writhe may also be used to determine the ratio of the bending (B) to twisting (C) rigidity, yielding B/C=0.23. This is about one-third of the value for long, random-sequence DNA, suggesting that the waf1/cip1/p21 p53 response-element is extremely flexible, a result that is also consistent with indirect biochemical experiments. These results support the idea, proposed by us earlier, that torsional stress may play a role in the regulation of p53 binding through modulation of twist at the binding site.