r in. ...,., .33.... . 3.... 4 2. .4 . :. .3. {3.13: a . . “ ....—.. 2 ES 1 2.; .. r... . u‘ .1 r’ I: ' 4mm ,- w nu: 3.22.3.1). llnlmunll‘l‘lfilllilllillnliliumlml 3 1293 01789 3144 LIBRARY Michigan State University This is to certify that the thesis entitled 1H NMR SPECTROSCOPIC STUDIES OF DINUCLEAR TRANSITION METAL CARBOXYLATE ADDUCTS OF DNA OLIGONUCLEOTIDES presented by Elizabeth Ursula Lozada Carrasco has been accepted towards fulfillment of the requirements for Masters Chemistry degree in fi/fl %4fi24 £42,“ / Major professor Date 7/07 j/QJ 0-7639 MS U is an Aflirmative Action/Equal Opportunity Institution PLACE IN RETURN BOX to remove this checkout from your record. TO AVOID FINES return on or before date due. MAY BE RECALLED with earlier due date if requested. DATE DUE DATE DUE DATE DUE 1M clam-9659.14 1H NMR SPECTROSCOPIC STUDIES OF DINUCLEAR TRANSITION METAL CARBOXYLATE ADDUCTS OF DNA OLIGONUCLEOTIDES By Elizabeth Ursula Lozada Carrasco A THESIS Submitted to Michigan State University in partial fulfillment of the requirements for the degree of MASTER OF SCIENCE Department of Chemistry 1998 ABSTRACT 1H NMR SPECTROSCOPIC STUDIES OF DINUCLEAR TRANSITION METAL CARBOXYLATE ADDUCTS OF DNA OLIGONUCLEOTIDES By Elizabeth Ursula Lozada Carrasco Cisplatin, cis-Pt(NH3)2Clz, is in clinical use as an effective treatment for cancer. Its ability to bind covalently at the N7 position of guanine bases and form intrastrand cross-links with DNA is thought to be responsible for its biological activity. Similar metal-DNA interactions involving other transition metal compounds have been implicated, but none have been elucidated to date. Studies in our laboratories revealed an unprecedented bridging mode for dimetal carboxylate compounds of 9-ethylguanine and 9- ethyladenine involving the N7 and O6/N6 atoms. These results led us to investigate the type of interactions that occur with oligonucleotides. 1H NMR spectroscopy was used to determine the solution structure of the duplex d(5 ’-CCTCTGGTCTCC-3 ’) - d(5 ’-GGAGACCAGAGG-3 ’) before and after reaction with [Rh2(OzCCH3)2(CH3CN)6][BF4]2. Differences in chemical shifts and intensities of the cross-peaks provided information about the perturbations of the oligonucleotide structure caused by covalent interactions with the metal at the GG sites of d(5’-CCTCTGGTCTCC-3’). TABLE OF CONTENTS . LIST OF TABLES V.“ V CHAPTER 1 Mechanism of Action of Cisplatin ...... 4 Dirhodium Carboxylate Compounds as Potential Anticancer Agents 14 NMR Analysis of Cisplatin-DNA Adducts 18 DJ CHAPTER 2 APPLICATION OF 1H NMR SPECTROSCOPY TO THE STRUCTURAL ELUCIDATION OF METAL-DNA ADDUCTS 21 Solvent Suppression Techniques .._.._..._ w 23 Two Dimensional 1H NMR Spectroscopy..-.._._,.__.-.-..-- 29 CHAPTER 3 EXPERIMENTAL SECTION 39 DNA Purification 40 NMR Sample Preparation .-__._ 45 1H NMR Analysis of a DNA Dodecanucleotide-.......-.-.............,.-.....--.-...._.-. 46 iii Reaction of the DNA Dodecanucleotide with a Dirhodium Acetate Complex 70 1H NMR Analysis of a RhZ-DNA Adduct 74 CHAPTER 4 MOLECULAR MODELING STUDIES 88 CONCLUSIONS 98 LIST OF RERERENCES 102 iv --.'!'} Table 1 . Table 2 . Table 3 . Table 4 . LIST OF TABLES Chemical Shifts for Exchangeable and Non-Exchangeable Protons in DNA DNA 12-mer : Relaxation and Calculated Correlation Times for AH8 Protons Base-to-Base and Base-to-Hl ’ Intra and Inter-residue Distances ha-DNA Adduct : Relaxation and Calculated Correlation Times for AH8 Protons 47 59 64 76 LIST OF FIGURES Figure l . Schematic of the Cisplatin Molecule 3 Figure 2 . Nucleotide Structure Showing Labeling Scheme .................................... 5 Figure 3 . Schematics Depicting The Basic Structure of DNA ............................ 7 Figure 4 . Cisplatin-DNA Inter— and Intrastrand Adducts ........................................... 9 Figure 5 . Biologically Relevant Pt Compounds 13 Figure 6 . Comparison of Cisplatin and Dimetal Carboxylate Binding Sites 15 Figure 7 . X-Ray Structures of Dirhodium Compounds with Guanine Bases ~ 17 Figure 8 . DNA Octamer Sequence Studied by 1H NMR Spectroscopy 19 Figure 9 . Sugar Conformations 19 Figure 10 . NMR Principle : Magnetization Vector 22 Figure 11 . NMR Pulse Experiment 23 Figure 12 . Binomial 1-1 Pulse Sequence 25 Figure 13 . Solvent Magnetization Vector for 1-1 Sequence .......................... 25 Figure 14 . Signal Intensity as a Function of Frequency Offset ..................... 26 vi Figure 15 . Figure 16. Figure 17. Figure 18. Figure 19 . Figure 20 . Figure 21 Figure 22 . Figure 23 . Figure 24 . Figure 25 . Figure 26 . Figure 27 . Figure 28 . Figure 29 . Figure 30 . Gradient Pulse Effect on Magnetization 27 vii WATERGATE Pulse Sequence 27 Spin State Observable and Non Observable Transitions ' 3O NOESY Pulse Sequence 32 T1 Measurements : Inversion Recovery Experiment 33 TOCSY Pulse Sequence 36 . 2Q COSY Pulse Sequence 38 HPLC Chromatogram of DNA Strands 42 DNA 12-mer 1H NMR Spectrum in 90% H20 at 25 °C 48 Imino Region of DNA Duplex at Different Temperatures 50 Base-to-Hl ’ NOE Connections 52 DNA Duplex WATERGATE NOESY Spectrum at 25 °C 54 2D Imino Region of DNA Duplex 56 DNA 12-mer 1H NMR Spectrum in 99% D20 at 25 °C 57 DNA 12-mer NOESY Spectrum in 99% D20 at 25 °C 60 DNA 12-mer 1H NMR Spectrum in 99% D20 at 30 0C 62 Figure 31 . DNA 12-mer NOESY Spectrum Figure 32 Figure 33 Figure 34 . Figure 35 . Figure 36 . Figure 37 . Figure 38 . Figure 39 Figure 40 Figure 41 . Figure 42 . Figure 43 Figure 44 . in 99% D20 at 30 °C 63 . DNA Duplex TOCSY Spectrum in 99% D20 at 30 °C 65 . DNA Duplex 2Q COSY Spectrum in 99% D20 at 25 °C 66 Base-to-Hl ’ NOESY Walk 68 Base-to-H3’ NOESY Walk 69 1H NMR Spectrum of [Rh2(02CCH3)2(CH3CN)6][BF412 in D20 71 HPLC Chromatogram of ha-DNA Adduct ...................................... 73 Imino Region of ha-DNA Adduct at Different Temperatures 75 . ha-DNA Adduct 1H NMR Spectrum in 99% D20 at 30 °C 77 . ha-DNA Adduct NOESY Spectrum in 99% D20 at 30 °C 80 ha-DNA Adduct 1H NMR Spectrum in 99% D20 at 20 °C 81 ha-DNA Adduct NOESY Spectrum in 99% D20 at 20 °C 82 . Base-to-Hl ’ Region of RhZ-DNA Adduct 83 Imino Region of ha-DNA Adduct in Buffer Solution 87 viii Figure 45 . DNA Backbone Showing Torsion Angles ............................................ 91 Figure 46 . DNA Structure Showing NOE-Distance Distance Restraints 94 Figure 47 . Superimposed Structures of d(CCTCTGGTCTCC)-d(GGAGACCAGAGG) 96 using Simulated Annealing ix LIST OF ABBREVIATIONS cis-DDP cis-Diamminedichloroplatinum (II) RNA Ribonucleic Acid DNA Deoxyribonucleic Acid A Adenine C Cytosine G Guanine T Thymine p Phosphate Group HMG High Mobility Group SRY Sex-determining Region Y protein LEF-l Lymphoid Enhancer Binding Factor 1 tsHMG Murine Testis Specific High Mobility Group trans-DDP trans-Diamminedichloroplatinum TBP TATA Binding Protein NMR Nuclear Magnetic Resonance JM2 1 6 Bis-Acetato-Ammine-Dichlorocyclohexylamine- platinum (IV) NOE Nuclear Overhauser Effects FID Free Induction Decay M0 Magnetization Vector My Transverse Magnetization Vector Bo Magnetic Field B1 External Magnetic Field WATERGATE Water Suppression by Gradient Tailored Excitation NOESY Nuclear Overhauser Effect Spectroscopy W1 First Order Transition Wo Zero Order Transition W2 Second Order Transition S Spin T1 Longitudinal Relaxation Time T2 Transverse Relaxation Time Tc Correlation Time rm Mixing Time t1 Delay Time t2 Acquisition Time xi F1 First Dimension F 2 Second Dimension FID Free Induction Decay LP Linear Prediction TOCSY Total Correlation Spectroscopy ms milliseconds ns nanoseconds ps picoseconds 2Q COSY Two Quantum COSY 1H Proton d Deoxy HPLC High Performance Liquid Chromatography nm nanometers min minutes mL milliliters g grams cm centimeters M Molar mM millimolar kDa kiloDaltons xii PPm SA RMSD MD Revolutions per minute parts per million Simulated Annealing Root Mean Square Deviation Molecular Dynamics xiii INTRODUCTION Inorganic compounds have been used extensively as medicinal agents for centuries. Presently, metal compounds are in clinical use for the treatment of a number of diseases such as arthritis, hypertension, bacterial infections and cancer. Regarding the latter, a great deal of research has been conducted in the last 30 years due to the discovery that cisplatin and other related complexes exhibit remarkable efficacy in the chemotherapeutic treatment against several types of cancer. Although the success of cisplatin as a chemotherapeutic agent has been well established and, indeed, it is currently in widespread use, its toxic side effects represent a considerable health risk for patients being treated with this drug. This disadvantage has provoked the search for new compounds with the same or higher activities but with lower toxicities than cisplatin. The use of other transition metals as a possible alternative to platinum has been proposed in a number of studies over the years. In this dissertation, a structural study of the principal binding target of cisplatin, namely, the intrastrand GpG cross-link is being investigated for dirhodium compounds. The main objective of this research is the elucidation of the mode of action of dimetal antitumor active compounds by understanding the structural distortions of DNA after metal binding has occurred. It is hoped that a better understanding of the DNA perturbations at the molecular level will lead to the discovery of new compounds that will satisfy the conditions necessary for inhibiting DNA replication in cancer cells. Chapter 1 1. BACKGROUND The history of the anticancer properties of the compound cisplatin, PtC12(NH3)2 or cis-DDP (Figure 1), dates back to 1965 when Professor Barnett Rosenberg and his research group observed filamentous growth of E.coli in experiments designed to study the effects of an electric field on the properties of these bacteria.1 The bacterial cells, which normally divide very rapidly, grew up to 300 times their size without undergoing division. This observation was eventually recognized to be due to the presence of platinum-ammine complexes formed by in situ electrolysis of the Pt electrodes used in the experiment. One of these complexes is cis- diamminedichloroplatinum(II) or cisplatin which has been shown to be effective in treating testicular, ovarian, lung, head and neck, bladder and cervical cancers. H 3 N $531 ’0... Pt . \\‘ / \ H 3 N Cl Figure 1. Schematic of the Cisplatin Molecule A. Mechanism of Action for Cisplatin After many decades of investigations, researchers have compiled a large database of information about how cisplatin is metabolized. The chloride ligands of cisplatin are stable with respect to substitution at high chloride concentrations such as those found in the extracellular matrix, but when the compound diffuses into the cell, the lower chloride concentration causes loss of chloride ligands and replacement with water.2 The resulting mono and bis-aqua complexes are reactive with nucleophilic sites of cellular macromolecules. The effect of this inorganic compound on the synthesis of RNA, DNA, and proteins has been evaluated, and it was found that at high concentrations and long incubation periods, the synthesis of these three biomolecules is inhibited.3 DNA synthesis is preferentially inhibited at all concentrations, thus DNA is thought to be the most likely target of cisplatin with respect to its chemotherapeutic properties. In order to understand the suppression of DNA synthesis by cisplatin, it is important to consider the nature of this macromolecule. DNA is a polynucleotide chain composed of nucleotides that are phosphate esters of pentose with nitrogenous bases linked to them (Figure 2). The standard nomenclature for nucleic acids establishes numbering from 1 to 9 for the *IIZ'LAK' “or. -. "'1' . base and numbers from 1’ to 5’ for the sugar moieties. The hydrogen atoms retain the same label of the carbon to which they are attached. In the case of where there is more than one hydrogen atom on the same carbon atom, ie, CX’ (X=H), the second hydrogen atom is labeled as CX”. The base is attached at the 1’ position of the sugar residue, which, in this case, is a deoxyribose since it has only one hydroxyl group located in the 3’ position. The bases are derivatives of purines and pyrimidines, and are responsible for the specificity of the sequence of DNA. This sequence defines the structure and coding for further transcription of DNA into RNA and final translation into proteins. The purine bases are guanine (G) and adenine (A) while the pyrimidine bases are cytosine (C) and thymine (T). The polynucleotide chain presents an orientation from the 5’ end to the 3’ end which are the positions where the phosphate group binds to the deoxyribose. N1 Figure 2. Nucleotide Structure Showing Labeling Scheme The most common form of DNA consists of two strands wound about a common axis in an antiparallel and helical fashion forming a right handed 5’ to 3’ helix. The strands are held together by hydrogen bonds formed between the bases of each strand, thereby forming base pairs. The most stable base pairs are those formed between guanine (G) and cytosine (C) and between adenine (A) and thymine (T). These bases are called the Watson- Crick base pairs. The double stranded DNA is then formed with the phosphate backbone at the periphery (hydrophilic region) in contact with the solvent, keeping the bases in the core (hydrophobic region) of the helix. Because the base-pairs are not exactly in the center of the helical axis, the DNA helix possesses two grooves with different widths and depths; these are the major and minor grooves (Figure 3). Double helical DNA can be classified into three types, namely, A, B, and Z-DNA, which exhibit fundamental differences in base stacking along the strand and in the conformation of the sugar-phosphate backbone. Among these three types, B-DNA is considered the native form because its X-ray pattern resembles that of DNA in intact sperm heads.4 .7-4——-—__m—_ a ‘_—\__‘ .. - v- IL ,H N 0 N ONIH N'H / H / Bar H—< I NAS- H:l‘< N/l N/LO Guanine Adenine Thymine Cytosine Watson-Crick Minor base pairs groove Figure 3. Schematics Depicting The Basic Structure of DNA The biological action of cisplatin is linked to the formation of stable Pt-DNA adducts. Major products of Pt-DNA binding involve the (N7,N7)- didcntate cisplatin cross-link between two intrastrand adjacent guanines (GpG) or adenine-guanine (ApG) (Figure 4a,b). Minor products are intrastrand cross-links between two guanine residues separated by one nucleotide residue and interstrand cross-lian between two guanine residues (Figure 4c,d).5 These cross-links have been implicated for the biological activity of cisplatin because they produce distortions in the conformational structure of DNA.6 Distortions observed in an X-ray structure of a DNA oligonucleotide of twelve base pairs bound to cisplatin was published in 1996 by S. J. Lippard and co-workers. In this structure, the DNA oligonucleotide exhibits a bent conformation with an angle of 35-40".7 Unwinding of the duplex was also observed. It is presumed that the distortions caused by cisplatin have deleterious effects on replication and transcription, thereby causing mutation.8 It is hypothesized that these lesions are recognized by certain proteins whose binding prevents the lesions from being repaired, a situation that ultimately leads to cell death. Such recognition is specific and requires a “kinked” DNA conformation.9 Figure 4. Cisplatin-DNA Inter- and Intrastrand Adducts 10 Some of the binding proteins have been identified as the high mobility group (HMG) domain proteins. Recently, it has been discovered that proteins responsible for mediating the biological processes of recombination, replication and transcription require the presence of DNA in an unwound and bent form and contain one or more HMG domains.10 HMG domains are a-helical structures that contain a high proportion of proline, aromatic and basic amino acids. Some of these HMG domain proteins, such as SRY and LEF-l, bind to specific DNA sequences via their hydrophobic residues. These domains behave as intercalators, causing unstacking of several base pairs, which leads to unwinding and bending of the DNA.11 Similar structural distortions have been observed when cisplatin binds to a DNA duplex.12 One of the latest reports on this subject, published in 1997, established the high binding affinity and specificity of the murine testis-specific high mobility group protein, tsHMG, to cisplatin- modified DNA.” It has been shown that the trans isomer of cisplatin, namely trans- DDP, enters the cells and binds to DNA like cisplatin.l4 Experiments performed in vitro have shown that, initially, trans-DDP binds more rapidly than cisplatin, but, with time, the quantity of bound trans-DDP diminishes 5 while the quantity of bound cisplatin continues to rise.1 These results 11 suggest that preferential repair of trans-DDP adducts occurs over those of cisplatin, although the details of how trans-DDP is removed are still not well understood. Interestingly, the HMG domain proteins studies also helped to explain the inactivity of the trans isomer. Now it is believed that the inactivity of trans-DDP is due to the specific binding of the HMG-domain proteins to cisplatin modified DNA. Although trans-DDP is capable of binding to only one strand in duplex DNA, it does so by forming a 1,3-intrastrand cross-link of guanine bases, trans-[Pt(NH3)2{ d(G*prG*)}]16, which leads to different structural alterations that are thought to be repaired more efficiently than the distortions produced by cisplatin. More recently, experiments using in vitro transcription challenge competition assays demonstrated the capability of the promoter recognition factor involved in transcription, the TATA binding protein (TBP), to preferentially bind cisplatin damaged DNA in favor of normal binding sites. The consequence of this is the prevention (or at least reduction) of transcription.17 Despite its high anticancer activity, the use of cisplatin as a drug is mediated due to its highly toxic side effects. Patients who receive the drug as part of their therapeutic treatment are likely to suffer from renal toxicity, 12 hematologic deficiencies, hearing loss, nausea, neurotoxicity and cardiac abnormalities. In response to these problems, other platinum compounds with reduced toxicity such as carboplatin have been developed (Figure 5a). Results indicate, however, that carboplatin may be less effective against testicular cancer with respect to long term cures.18 Presently, other compounds such as JM216 19 (Figure 5b) and oxaliplatin20 (Figure 5c) are in various stages of clinical testing. The former has been shown to exhibit lower toxicities than cisplatin and carboplatin. The latter appears to be promising for colon cancer treatment, but its use might be limited due to its high levels of neurotoxicity. In spite of the promise of metal-based drugs, the battle against cancer continues, with studies being focused on both improved platinum compounds and on non-platinum metal compounds. New compounds may exhibit reduced toxicities, improved activities against cisplatin-resistant tumors or efficacy against cancers that are currently untreatable. 13 Cl In, \\\ N ” ...... ‘\ 2 OIO> Pt/ ""2 CI NH3 0 O \ NH2 JM216 Oxaliplatin O (b) (c) 2+ 0 I’m", Pt,.\\\NH 3 NH 3/01.. Pt..\\\\Cl NHa’ ‘NH2(CH2).H2N/ Na, (d) 2+ NH3I/,,..Pt...\\\Cl Cllm..,Pt,..\\NH3 NH,’ ‘NH2(CH2),,H2N/ ‘NH3 (8) Figure 5. Biologically Relevant Platinum Compounds . d—--.—l-—r‘l—‘ — - 1 ‘ 14 B. Rhodium Carboxylates as Potential Anticancer Agents In the search for new chemotherapeutic agents, studies were initiated with antitumor active compounds that may be capable of binding to bases in a manner similar to cisplatin. Dinuclear platinum and rhodium compounds fall into this category. Because of their wider “bite angle” for binding as compared to cis-DDP, there is a greater possibility that dinuclear compounds can form both intrastrand and interstrand cross-links.21 Indeed, dinuclear platinum complexes structurally related to cisplatin such as the cis- and trans- isomers of [{PtCl(NH3)2}2H2N(CH2),,NH2]2+ (Figures 5d and Se), have been tested for anticancer activity in mice and were shown to be effective against cisplatin resistant tumors. 22 The DNA binding mechanism of these compounds may be different from that of cisplatin, since it was found that in a d(GCGC) sequence, the preferred binding sites are interstrand over the favorite intrastrand -GG- cross-link of cisplatin.23 Some of these dinuclear complexes promote a change in the conformation of DNA from the B type to Z-DNA.24 Like cisplatin, they unwind DNA but they do not bend DNA in the same way. This may play a role in their efficacy towards cisplatin- resistant cells. 15 Among the group of non-platinum compounds that may exhibit binding modes similar to cisplatin are dimetal carboxylates. Instead of one square planar metal unit, they are composed of two approximate square planar units joined by a metal-metal bond (Figure 6). Over the years, it has been shown in our laboratories and others that the dimetal tetracarboxylate compounds react with incoming bases by displacing two cis RCOz' ligands, which means that they could react in a similar manner to cisplatin. /><’\O H3N/ NH3 0 O O \ .Pt’ — L Cl Cl 0 ’ 0 90° 900 Figure 6. Comparison of Cisplatin and Dimetal Carboxylate Binding Sites During the 1970’s, studies of dirhodium carboxylate compounds revealed their potential as chemotherapeutic agents against certain tumors including Ehrlich ascites and L1210 tumor cells in mice. The strong inhibitory effect of the compounds on DNA but not on RNA synthesis in viva was also demonstrated.25 In these experiments, the anticancer activity of the compounds Rh2(OzCCR)4L2 (R = CH3, C3H7, C4H9) was explored and it was found that the propionate and butyrate compounds exhibited higher 16 potencies than the acetate complex. The level of toxicity however, was invariant. Among the dinuclear rhodium carboxylates, the compound Rh;(DTolF)2(OzCCF3)2(HzO)2 (DTolF = N,N’-di-pftolylfonnamidinate) was also tested for antitumor activity with promising results against Yoshida sarcoma and T8 sarcoma of Guerin being obtained.26 This compound exhibits reduced toxicity as compared to Rh2(OzCCH3)4 and cisplatin. The early literature regarding the DNA binding studies of dirhodium carboxylates reported that these compounds were unreactive towards poly G DNA sequences but reactive towards poly A sequences and single-stranded DNA. This conclusion was based on the absence of color change during the reaction with guanine bases, while in the reactions with adenine bases a dramatic color change from blue-green to pink was observed. Reactions performed with 9-ethylguanine in our laboratories revealed that gaunine does indeed bind and in an unprecedented manner wherein two guanine bases are bound in a bridging mode. X-ray structures were determined and the compounds were found to contain cis- guanines with the N7 and 06 atoms coordinated in a “head-to-head” bridging mode (Figure 7a) or a “head-to-tail” bridging mode (Figure 7b).27 The reactions were performed with 9-ethyladenine as well and it was found that N7,N6 bridges formed. These results support the conclusion that dirhodium tetraacetate has l7 0 V - .' 1" . ll 3" "’ :’ H8 06’ ‘ . sf. 7C8 H89 07 ' Rh2 "' os 01 1'3- 1004 ‘a . 2 ,’__ 0 &4 (b) Rh2(02CCH3)2(9-EIG)2(MCOH)2.ZMCOH NHz F N HN f .s\\\\\\ T M‘— T u (y Head-to-Head N/ e . HN /) \\\\\\\ ‘3 (“—07 V Head-to-Tail T T. Figure 7. X-Ray Structures of Dirhodium Compounds with Guanine Bases 18 equatorial sites available for binding as well as axial sites. C. N MR Analysis of Cisplatin-DNA Adducts The aforementioned studies of dimetal compounds, which included single cystal X-ray structures of model purine reaction products, were the starting point for our studies of the DNA binding of dimetal carboxylate compounds. We began our current investigation from the perspective of what is known about cisplatin in order to make comparisons and contrasts to the behavior of dirhodium carboxylates under the same conditions. Determination of the structural distortions of cisplatin-DNA adducts at the molecular level is important for the eventual understanding of cisplatin’s activity. In addition to X-ray crystallography, NMR spectroscopy is an extremely valuable tool for probing DNA-adducts of cisplatin. The full mechanism of action of cisplatin as an anticancer agent may eventually be understood on the basis of both solution and solid-state structures of DNA/cisplatin/protein interacions. Initial analyses of the distortions created by platinum complexes bound to DNA were conducted with di- and tn'nucleotides. The small size of these adducts renders them excellent candidates for study by NMR l9 spectroscopy, molecular mechanics and X-ray crystallography.28 Distortions upon platination of DNA duplexes of ten, eleven and twelve base pairs were also investigated by 1D and 2D NMR spectroscopy.29 In 1995, a more detailed 2D NMR analysis of an intrastrand cis-DDP-GG cross-link structure of an octameric DNA duplex was performed.30 The DNA sequence selected for these spectroscopic studies was d(CCTGGTCC) d(GGACCAGG), represented in Figure 8. 1234 s 678 5’CCTG*G*TCC 3’ 3’GGACCAGG5’ 1615 14 13 12 11 1o 9 G*G* = Platinated site Figure 8. DNA Octamer Sequence Studied by 1H NMR Spectroscopy Some of the results of this work are the following: (1) The conversion of the sugar conformation of the 5’ platinated guanosine from C2’ endo to C3’ endo was confirmed (Figure 9), although the remainder of the sugar residues of the double helix were in the original B type DNA conformation (C2’ endo). cs' 2' N C5' 3' 0 w 4' 3' 1' 4' 2. 1' C2' endo; S type C3' endo; N type Figure 9. Sugar Conforrnations 20 (2) The connectivities of both strands were generally uninterrupted with some variations in the intensity of the NOE cross-peaks from protons near the lesion site (Pt-G4pG5). (3) The chemical shifts of G4*H8 and G5*H8 resonances were downfield in comparison with the other G resonances, indicating that both bases are platinated at their respective N7 positions. (4) The G4 imino proton was not observed, presumably due to its rapid exchange with water. The NOE restrained, refined structure of the modified duplex showed that the helix is kinked approximately 58° and unwound by —21°. It should be noted that the structural features of this DNA sequence were obtained by using 1H NMR spectroscopic data only, underscoring the importance of this analytical technique for these types of structural studies. Chapter 2 1. APPLICATION OF NMR SPECTROSCOPY TO STRUCTURAL ELUCIDATION OF METAL-DNA ADDUCTS The application of NMR spectroscopy to the study of biological macromolecules involves the analysis of NMR spectra in two or more dimensions. In these experiments, it is possible to observe through-bond couplings for nuclei separated by two, three or sometimes four bonds. Also, one can observe through-space connections between resonances of nuclei that are close to each other in space (55 A), and exchange processes between resonances of the same nucleus in two different chemical environments. The structural information is principally obtained from measurements of the Nuclear Overhauser Effect (NOE), which provides constraints on internuclear distances, especially for small oligonucleotides. In addition to the determination of complete three-dimensional structures for macromolecules in solution, NMR provides valuable information on local structure, conformational dynamics and structural aspects of interactions with small molecules. 21 22 The principles of NMR spectroscopy are based on the behavior of the nuclear spins when they are placed in a static magnetic field (Bo). Each nucleus precesses around the field axis and by adding up all the 2 components of the nuclear magnetic moments, a macroscopic magnetization vector (Mo) along the field direction is obtained (Figure 10). In an NMR experiment, transitions are induced between the energy levels by irradiating the nuclei with a external magnetic field (Bl) originated from an electromagnetic wave of the suitable radiofrequency. Mo= Magnetization a vector I AB = y(h/21t)Bo m = +1/2 [3 m = - 1/2 Figure 10. NMR Principle : Magnetization Vector In the pulse experiment, a radio frequency pulse of duration 1,, (Figure 11a) is obtained by applying the magnetic field (B1) along a specific 23 direction. This can be observed in Figure 11b where the x and y axes are rotating at the same frequency of B1. The magnetization vector M0 is tipped away from the z axis by an angle that is dependent on the duration of the pulse and the strength of the magnetic field applied. The resultant vector (My) is called the transverse magnetization vector. This vector is important because it represents the NMR signal since the detector contains the receiver coil aligned along the y direction. (a) * —-> "p <—- My=Mo Bl = external magnetic field applied Figure 11. Pulse Experiment A. Solvent Suppression Techniques For a thorough understanding of NMR processes in nucleic acids, it is important to differentiate between two types of protons in these biomolecules. This classification stems from the ability or inability of these protons to exchange with protons of the solvent. The exchangeable types 24 include the protons A (NH2)6, G (NH2)2 and G(NH)1, T(NH)3, and C(NH2)4 with the remainder of the protons present in the molecule being non-exchangeable. In order to detect the exchangeable protons, the sample must be dissolved in a protonated, rather than deuterated, solvent. Unfortunately, this presents the problem of dealing with a very high concentration of water compared with the sample concentration, diminishing tremendously the dynamic range of the instrument and the intensity of the resonances of interest. For this type of analysis either selective excitation or solvent suppression NMR techniques must be performed in order to overcome this problem. Among the selective excitation methods, the binomial method is commonly used. The method uses a selective excitation pulse sequence which produces a non-excited region around the solvent chemical shift while causing significant excitation in the region of interest. The magnetization process can be explained by looking at one of the simplest pulse sequences, the 1-1 sequence (Figure 12). In this sequence, the excitation pulses are placed at one end of the spectrum, separated by N Hz, where N represents the distance from the solvent resonance frequency to the region to be excited. All the pulses have the same phase and are separated by a delay time (2N'l ) in seconds. 25 1c/4x n/4 , 1 tD Acquire tD = 2N sec Figure 12. Binomial 1-1 Pulse Sequence By observing the magnetization process for the solvent signal using the vector diagram shown in Figure 13, one can observe that after the two pulses, the magnetization vector returns to its initial position along the z axis, therefore yielding no net excitation of this signal. 112 ‘1 3.?!L/<%\}P,%\flf />\ r ‘4 \4/ \q/ ’\.. \ ., T _. .. _ L F 2 .1 .1 T .1... ‘ 0. mm a 03.x? a 6.58% E22 E. 58$ 5flo 4:0 3:0 2:0 1 HI Fl (ppm) Figure 42. ha-DNA Adduct NOESY Spectrum at 20 °C l . 0 F2 (Ppm) Hl 83 ’43:}: H8H 0H g. g -“l i‘k’ -‘ a 2'1! Inge” 5 - 3‘: - ”)3! 9:: -6 f q ”A rq‘,. ._ _ P‘ : * AL¢§LP 5 he?“ 5.29"} .A‘ Q " G24l~l8,Hl’ “ >~ Figure 43. Base-to-Hl ’ Region of ha-DNA Adduct 84 frequency before returning to their original chemical shifi, they will be out of phase with the spins that did not exchange. The accumulation of these phase errors over time will lead to destructive interference in the FID and loss of the signal. As a result, the FID decays rapidly and the Fourier transformation produces a broad line. This is known as exchange broadening. The greater the chemical exchange rate, the more frequently the spin visits another site with different precession frequency, and the more rapidly will be the loss of phase coherence and therefore the greater the line broadening will be. Several attempts to overcome the severe instability of the ha-DNA adduct were performed. One of the parameters that was modified to remedy this situation is temperature. Due to the fact that it was possible to obtain a NOESY walk of the H1’-to-base region of one of the strands when the temperature was lowered from 30 °C to 20 °C, it was decided to lower the temperature even further, down to 5 °C. Focusing on the same H1 ’-to-base region, it was observed that the breadth of the lines due to the decrease in the correlation time of the molecule at this temperature negatively affected the cross-peaks of the complimentary strand while causing only a small improvement in the resolution of resonances corresponding to the metallated strand. 85 A second factor that was modified is the ionic strength of the solution. By increasing the salt concentration of the sample, the negative repulsion of the phosphate backbone of the oligonucleotide is diminished. The objective is to minimize the instability of the metallated DNA as much as possible. By decreasing the electrostatic repulsion, it is possible to stabilize the adduct such that a 2D spectrum with more defined cross-peaks may be obtained. A buffer solution of 1M NazHPO4/NaH2PO4 (pH = 6.4) and 1M NaCl was added to the D20 sample to obtain a final salt concentration of 50 mM NazHPO4/NaH2PO4, 150 mM NaCl and 0.1 mM NaC2H6As02 at pH = 6.5 in a final volume of 0.250 mL. It is important to note that by lowering the pH of the solution, chemical exchange may be diminished as well. A 2D NOESY spectrum in at 5 °C was measured under these new buffer conditions. After data processing, the H1 ’-to-base region of the spectrum showed the cross-peaks of both strands mixed together in the same ill-defined region. We take this to mean that the degree of stability of the sample was dramatically changed to the point of causing aggregation which would lead to poor resolution of the 2D spectrum. By adding the same buffer solution to a sample in 90% H20, the change in stability of the sample was confirmed by observing the imino region at different temperatures (Figure 44). The breadth of the resonances increased and, as a consequence, 86 the signal-to-noise ratio diminished; however, the stability of the metallated duplex increased. This is clearly observed in the spectra. Even at 40 °C, some of the imino proton resonances could still be observed. From these aforementioned experiments, it was determined that 20 °C is the best temperature for measuring the 2D NOESY spectrum. The NOESY spectrum at 20 °C, however, revealed that the molecule is quite unstable especially in the region of the metallated strand. This instability is exacerbated by the decrease in signal-to-noise ratio, which we were able to partially overcome by acquiring twice the number of scans. Unfortunatly this only increased the signal—to-noise ratio by 40 percent. For this reason, more assignments of the cross-peaks of the spectrum are not possible to make at this time. It is evident that further manipulation of the sample is required if we are to obtain a better resolved spectrum with more dispersed cross-peaks in both dimensions. The attempts that have been made to date lead us to conclude that this will not be a straightforward task, and that additional factors such as concentration and use of buffer solutions will have to be taken into account. 5°C fl 1 r I g i . ' {l f. i. ;' “ a. [ . r} l g / «l M \ f l . lwn F \l-A l _ . ~wkm‘ufiw 1* ”M “flail-w “‘40 "‘1‘; 14.5 14.0 13.5 13.0 12.5 12.0 ppm 15 °C , .7) ‘ ,1. .': ,"V V l ,1 l g 'I A ,I ' / "‘1 M ‘ I l \ f ‘. N” l r \ V)- . ‘ WWL‘Jx‘V W QK'IVI N‘WMI M”. 19"“; —'\ {k 1 14.5 14.0 13.5 13.0 12.5 12.0 ppm 20 °C ' .' H .. ,1 W \J V r» l \ I l f l j , If . I“. \\ ("V V l f" V)‘ NM‘MWW _ MW" 'MJM‘WYW‘Y 14.5 14.0 13.5 13.0 12.5 12.0 ppm 25 °C .1 :1 W! IV" R) lb!) l l ,1 \ (l ) 1 1V! 0 \‘n Wim‘i’fi V '..~~f". _,.___“Tv~"f¢ ."W 14.5 14.0 13.5 13.0 12.5 12.0 ppm 87 30 °C 7 1" WMM'J/ )VA‘NM/ A- . . 14.5 14.0 13.5 13.0 40 °C ll ‘- W" ‘J ¢.‘)V~‘A~1FVJI‘.N Y‘ k’kfl‘fhfi'fil‘fi ‘A‘VWFJV’WM 14.5 14.0 13.5 13.0125 12.0 ppm 45 °C 1.1" H) V II) I it 1' v x . l x. I l i ‘N r N LN M.-wt~.~i‘-m~rw‘»"*f*f .‘Ww-N, , . , _ .Wwwfi 14.5 14.0 13.5 13.0 12.5 12.0 ppm 49 °C 14‘ J» l H A ‘ A I If “‘4 JY .- W‘s"A‘i‘.’--i-’v~xt‘r’t7in:it";. .“‘. MW: “5*. , ... WMMH 14.5 14.0 13.5 13.0 12.5 12.0 ppm Figure 44. Imino Region of ha-DNA Adduct in Buffer Solution MM 12.5 12.0 ppm Chapter 4 1. MOLECULAR MODELING STUDIES The ultimate goal of 1H NMR spectroscopic measurements on large biomolecules is the modeling of structural features that will satisfy the NMR data as well as chemical information. The structure of the DNA double helix varies according to its sequence and physical environment. Intemuclear distances, torsion angles, bond lengths and bond angles are characteristics of the molecule that provide important restraints for the determination of its structure. Structural determination involves taking the restraints together with the covalent structure and inputting into them a computational algorithm that searches the conformational space. This method provides a series of possible or feasible structures that satisfy the assigned structural restraints. Determination of DNA structures by 1H NMR spectroscopy is an alternative to X-ray crystallography, and also provides valuable complementary information about solution structure. The main difference between these two techniques is that in X-ray crystallography typically only one structure with a minimum energy in the solid state is obtained. 88 89 However, with NMR spectroscopy, according to the quality and number of restraints that can be extracted from the multidimensional experiments, a set of similar structures of varying energies are obtained as expected for molecules in the solution state. New techniques and computational programs are continually being developed in order to establish more accurate restraints and, consequently, to narrow the set of possible structures for the molecules. As stated earlier, the intensity of the NOE cross-peaks is inversely related to the distance between two nuclei. Therefore, an accurate measurement of the volume of each of the cross-peaks will give rise to an array of inter-proton distance restraints. Interestingly, the two protons that give rise to a cross-peak are not alone in the molecule.» They belong to a coupled network. However, assuming that each cross-peak is due to two isolated nuclear spins is a good approximation; this is referred to as the two spin or Isolated Spin Pair Approximation (ISPA) represented by the following equation:37 r ab: r..r(a ref/a ab)1/6 where rah is the interproton distance between nuclei a and b that needs to be determined, aab is the corresponding NOE cross-peak intensity and ad and ab are the known interproton distance and cross-peak intensity respectively. 90 The reference distances can be obtained from fixed covalent distances within the molecule under investigation. These distances can be related to the cross-peak intensities in the NOE spectrum. However, it is important to note that the mixing time (rm) of the NOE experiment plays a critical role in the accuracy of the NOE measurements that can be obtained. A limitation of ISPA is the omission of factors such as spin diffusion that may give rise to cross-peaks of greater intensities than those expected, which would lead to an inaccurate assessment of distances. Structural information can be obtained by using torsion angle restraints and NOE distance restraints. To obtain data of the former, correlation spectroscopy (COSY) in any of its variations must be performed. From these experiments, measurements of coupling constants provide torsion angle restraints in the sugar deoxyribose rings 38 which consequently help to define the conformation of the DNA sugar phosphate backbone (Figure 45). 91 Figure 45. DNA Backbone Showing Torsion Angles Regarding NOE distance restraints, simulated annealing (SA), which is a form of restrained molecular dynamics, relates the experimental structural restraints with energy considerations. In this approach, the potential energy is calculated for a set of atomic coordinates using a force field whose limits are determined by the values established for the restraint bounds and standard covalent bond parameters (lengths, angles). The farther the NOE distance restraints are from these bounds, the higher the potential energy of the molecule and therefore, the greater the degree of instability of the structure. 92 In simulated annealing, an initial structure is selected. This structure is allowed to search conformational space. After that, the set of restraints, (torsion and/or distance) is applied. These restraints will be responsible for the final structure obtained. By repeating the (process several times, the structures produced at the end, although not identical, should be very similar with a small atomic root mean square deviation (RMSD) between them (<1A). The different structures are averaged and a final structure is reported. F urtherrnore, several structures can be considered as initial structures, e.g., A and B-type DNA, with the expectation that similar structures will be obtained if the restraint set is sufficient. For the dodecanucleotide under study, the NOESY data obtained from the spectrum measured at 30 °C gives rise to assignments of 718 NOE cross- peaks from both sides of the diagonal. These assignments were introduced in a table under FELIX 97.0 (Molecular Simulations Inc.) which generated the volumes table that represents the intensity of each cross-peak. The reference distance and cross-peak intensity were obtained from the C19 residue, specifically, Cl9:H5-C19zH6 covalent bond distance (2.464 A) and C l9:H2’-Cl9:H2” distance (1.785 A). From these tables, the restraints file was obtained with a total of 320 restraints calculated. Some of the restraints 93 were ignored due to overlapping problems that would leads to erroneous distance information. The restraints file was exported to the structural calculation program InsightII (Molecular Simulation Inc.). Within InsightII, the simulated annealing module was selected for the refinement calculations. Together with the NOE distance restraints obtained from the NOE spectrum in D20, some additional restraints regarding the imino-imino cross-peaks from the NOESY spectrum in HZO were included. To enforce proper base pairing, hydrogen bonding restraints were also used for all Watson-Crick base pairs in the sequence. The starting structure selected was the sequence in standard B form conformation. The set of restraints was displayed on this molecule and this result is presented in Figure 46. The simulated annealing protocol contains a Molecular Dynamics schedule (MD_Schedule) edit command that provides a flexible and easy- to-use interface for designing new protocols according to the requirements of the molecule under study. For the DNA oligonucleotide, the restrained molecular dynamics protocol begins with a minimization of the B-form DNA duplex with gradual application of the restraints and the force constants by 10 percent increments. The molecule was allowed to adopt randomized velocities from the beginning of the protocol. The simulated 94 wuflabmvm + Ohio—Cum fiafibmom 853me was oozSflQflOZ wagonm BBQ—Em <75 .3. 0.53m 33.38% 8:8me 28 ooqfimaéoz Baogm <75 N,“ V C u—A‘. b ‘7»-..3 ‘0.!‘- I 9......“ a A A... "‘2‘; u A... "‘1‘\. '4“le ...4" r'.’~'1.| 95 annealing Quartic system applies the non-bond repulsion energy values gently. This first step in the protocol is followed by two dynamic processes starting with a steepest descent minimization algorithm. The molecular structure is heated from O to 300 K for 1000 ps under fully restrained conditions. The second dynamics process involves a conjugate gradient minimization algorithm that uses the cartesian coordinates of the atoms as the space variables, keeping the temperature at 300 K. These conditions led to sufficiently randomized starting geometries that, after a final conjugate gradient minimization, gave rise to a structure that lies well within the standard B-form DNA family. Some structural distortions are observed as were expected, especially at the ends of the molecule. The final structure obtained had a total energy of approximately 2202 Kcal. Energy increments due to NOE distance violations were lower than 200 Kcal comprising less than 10% of the overall energy of the molecular structure. Figure 47 shows superimposed structures resulting from ten simulations, each starting from a different trajectory. Convergence to an atomic RMSD of 0.87 A was achieved using this protocol. The ten structures overlapped with fairly good reproducibility. Most of the differences among the structures are observed near the ends of the duplex. This is due to the lower number of restraints in these regions of the 96