Aptaligner: automated software for aligning pseudorandom DNA X-aptamers from next-generation sequencing data

Next-generation sequencing results from bead-based aptamer libraries have demonstrated that traditional DNA/RNA alignment software is insufficient. This is particularly true for X-aptamers containing specialty bases (W, X, Y, Z, ...) that are identified by special encoding. Thus, we sought an automated program that uses the inherent design scheme of bead-based X-aptamers to create a hypothetical reference library and Markov modeling techniques to provide improved alignments. Aptaligner provides this feature as well as length error and noise level cutoff features, is parallelized to run on multiple central processing units (cores), and sorts sequences from a single chip into projects and subprojects.

X-aptamers: a bead-based selection method for random incorporation of druglike moieties onto next-generation aptamers for enhanced binding

By combining pseudorandom bead-based aptamer libraries with conjugation chemistry, we have created next-generation aptamers, X-aptamers (XAs). Several X-ligands can be added in a directed or random fashion to the aptamers to further enhance their binding affinities for the target proteins. Here we describe the addition of a drug (N-acetyl-2,3-dehydro-2-deoxyneuraminic acid), demonstrated to bind to CD44-HABD, to a complete monothioate backbone-substituted aptamer to increase its binding affinity for the target protein by up to 23-fold, while increasing the drug's level of binding 1-million fold.

NMR structure of human thymosin alpha-1

800 MHz NMR structure of the 28-residue peptide thymosin alpha-1 in 40% TFE/60% water (v/v) has been determined. Restrained molecular dynamic simulations with an explicit solvent box containing 40% TFE/60% TIP3P water (v/v) were used, in order to get the 3D model of the NMR structure. We found that the peptide adopts a structured conformation having two stable regions: an alpha-helix region from residues 14 to 26 and two double β-turns in the N-terminal twelve residues which form a distorted helical structure.

A pyrazolato-bridged dinuclear platinum(II) complex induces only minor distortions upon DNA-binding

The cytotoxic, pyrazolato-bridged dinuclear platinum(II) complex [(cis-{Pt(NH3)2})2(mu-OH)(mu-pz)]2+ (pz=pyrazolate) has been found to cross-link two adjacent guanines of a double-stranded DNA decamer without destabilizing the duplex and without changing the directionality of the helix axis. A 1H NMR study of the oligonucleotide d(CTCTG*G*TCTC)-d(GAGACCAGAG), cross-linked at the two G* guanines by [(cis-{Pt(NH3)2})2(mu-pz)]3+, and molecular dynamics simulations of the explicitly solvated duplex were performed to characterize the structural details of the adduct. The dinuclear platinum cross-link unwinds the helix by approximately 15 degrees , that is, to a similar extent as the widely used antitumor drug cisplatin, but, in contrast to the latter, induces no significant bend in the helix axis. The Watson-Crick base-pairing remains intact, and the melting temperature of the duplex is unaffected by the cross-link. The helical twist is considerably reduced between the two platinated bases, as becomes manifest in an unusually short sequential H1'-H1' distance. This unwinding also affects the sugar ring of the guanosine in the 3'-position to the cross-link, which presents an N<-->S equilibrium. This is the first cytotoxic platinum complex that has been successfully designed by envisioning the structural consequences of its binding to DNA.

Fixing the conformations of diamineplatinum(II)-GpG chelates: NMR and CD signatures of individual rotamers

The bulky, asymmetric analog of the antitumor drug cisplatin, [PtCl(2)(tmen)] (tmen = N,N,N'-trimethylethylenediamine), was used to produce crosslinks with the dinucleotide d(GpG), modeling the most frequent lesions that cisplatin and its analogs cause to DNA. The ligand tmen was chosen because it is expected to constrain the guanine cis to the NMe(2) group in the adduct [Pt(tmen){d(GpG)}](+) to an orientation perpendicular to the coordination plane and to stabilize the other guanine in an oblique orientation, thus maintaining a head-to-head geometry typical of cisplatin-d(GpG) crosslinks within single- and double-stranded DNA. Of the four possible combinations of tmen chirality (R or S symmetry of the coordinated NHMe group) and crosslink direction (5'-G bound cis to the secondary or the tertiary amino group of tmen), two isomers were preponderantly formed, [Pt(R-tmen){d(GpG)}](+) with 5'-G bound cis to NMe(2) and [Pt(S-tmen){d(GpG)}](+) with 5'-G bound cis to NHMe. The former was shown to have a right-handed R2 orientation of guanines similar to that found in duplex DNA, whereas the latter had a left-handed L1 orientation that modeled cisplatin-d(GpG) adducts within single-stranded DNA. The R2 rotamer was found to be in an equilibrium (as observed using EXSY spectroscopy) with a minor fraction (< or =4%) of a Delta-HT rotamer related to R2 by rotation of the 3'-G about the Pt-N7 bond. The major rotamers R2 and L1 were isolated using reverse-phase HPLC, and their NMR and CD signatures were compared to those of the corresponding rotamers of the less hindered adduct [Pt(dmen)(GpG)](+) (dmen = N,N-dimethylethylenediamine). From this and other comparisons with previously reported platinum dinucleotide complexes, and from molecular modeling, it could be concluded that both steric repulsion between guanine and substituents of the cis amino group and N-H...O6 hydrogen bonding are significant effects favoring the oblique orientation of one guanine base typical of the HH rotamers of [Pt(diamine){d(GpG)}](+) and [Pt(diamine)(GpG)](+) complexes.

Cross-links of quadruplex structures from human telomeric DNA by dinuclear platinum complexes show the flexibility of both structures

The folding of AG(3)(T(2)AG(3))(3) was investigated in the presence of Na(+) or K(+) ions, by using the dinuclear platinum complexes [{trans-PtCl(NH(3))(2)}(2)H(2)N(CH(2))(n)NH(2)]Cl(2) (n = 2 or 6). AG(3)(T(2)AG(3))(3) has been previously found to adopt two different quadruplex structures: the antiparallel one in a solution containing Na(+) and the parallel one in a K(+)-containing crystal. The two structures are strikingly distinct and are not expected to form the same platinum cross-links. Therefore, characterization of the cross-links formed with platinum complexes in solution allowed the predominant conformation(s) to be identified. The bases coordinating the platinum atoms were identified by chemical and 3'-exonuclease digestions. The observed cross-links showed that the parallel structure exists in solution whatever the cation and confirmed the existence of the antiparallel structure in the presence of both cations as previously reported from cross-linking experiments of AG(3)(T(2)AG(3))(3) by mononuclear platinum complexes. Furthermore, the major platinum cross-links were unexpectedly formed between two guanines belonging to the same G-quartet. Their formation was rationalized using molecular dynamics simulations in implicit solvent of the two quadruplex structures. It was shown that they were flexible, allowing some guanines to leave reversibly the top G-quartet and thus rendering their N(7) atom accessible to platinum complexes. Our results also suggest that the human telomere sequence could be a target for such platinum complexes.

Platinum cross-linking of adenines and guanines on the quadruplex structures of the AG3(T2AG3)3 and (T2AG3)4 human telomere sequences in Na+ and K+ solutions

The quadruplex structures of the human telomere sequences AG3(T2AG3)3 I and (T2AG3)4 II were investigated in the presence of Na+ and K+ ions, through the cross-linking of adenines and guanines by the cis- and trans-[Pt(NH3)2(H2O)2](NO3)2 complexes 1 and 2. The bases involved in chelation of the cis- and trans-Pt(NH3)2 moieties were identified by chemical and 3'-exonuclease digestions of the products isolated after denaturing gel electrophoresis. These are the four adenines of each sequence and four out of the 12 guanines. Two largely different structures have been reported for I: A from NMR data in Na+ solution and B from X-ray data of a K+-containing crystal. Structure A alone agrees with our conclusions about the formation of the A1-G10, A13-G22, A1-A13 platinum chelates at the top of the quadruplex and A7-A19, G4-A19 and A7-G20 at the bottom, whether the Na+ or K+ ion is present. At variance with a recent proposal that structures A and B could be the major species in Na+ and K+ solutions, respectively, our results suggest that structure A exists predominantly in the presence of both ions. They also suggest that covalent platinum cross-linking of a human telomere sequence could be used to inhibit telomerase.

Recognition of DNA interstrand cross-link of antitumor cisplatin by HMGB1 protein

Several proteins that specifically bind to DNA modified by cisplatin, including those containing HMG-domains, mediate antitumor activity of this drug. Oligodeoxyribonucleotide duplexes containing a single, site-specific interstrand cross-link of cisplatin were probed for recognition by the rat chromosomal protein HMGB1 and its domains A and B using the electrophoretic mobility-shift assay. It has been found that the full-length HMGB1 protein and its domain B to which the lysine-rich region (seven amino acid residues) of the A/B linker is attached at the N-terminus (the domain HMGB1b7) specifically recognize DNA interstrand cross-linked by cisplatin. The affinity of these proteins to the interstrand cross-link of cisplatin is not very different from that to the major 1,2-GG intrastrand cross-link of this drug. In contrast, no recognition of the interstrand cross-link by the domain B lacking this region or by the domain A with or without this lysine-rich region attached to its C-terminus is noticed under conditions when these proteins readily bind to 1,2-GG intrastrand adduct. A structural model for the complex formed between the interstrand cross-linked DNA and the domain HMGB1b7 was constructed and refined using molecular mechanics and molecular dynamics techniques. The calculated accessible areas around the deoxyribose protons correlate well with the experimental hydroxyl radical footprint. The model suggests that the only major adaptation necessary for obtaining excellent surface complementarity is extra DNA unwinding (approximately 40 degrees ) at the site of the cross-link. The model structure is consistent with the hypothesis that the enhancement of binding affinity afforded by the basic lysine-rich A/B linker is a consequence of its tight binding to the sugar-phosphate backbone of both DNA strands.

Unrestrained 5 ns molecular dynamics simulation of a cisplatin-DNA 1,2-GG adduct provides a rationale for the NMR features and reveals increased conformational flexibility at the platinum binding site

A 5 ns unrestrained molecular dynamics (MD) simulation of the DNA duplex d(GCCG*G*ATCGC)-d(GCGATCCGGC), bearing a cis-Pt(NH(3))(2)(2+) unit crosslinking the two G* guanine bases, is reported. The MD trajectory was a posteriori correlated with NMR data determined for the same adduct, and it is shown that interproton distances and the characteristic chemical shifts are accounted for by the simulation. The simulation and its confrontation with the NMR data have confirmed the finding derived early from static models that the cytosine complementary to the 5' G*, C17, is mobile with respect to its adjacent bases. However, in contrast to our previous description of this mobility, which included rupture of the Watson-Crick hydrogen bonds and formation of non-Watson-Crick hydrogen bonds, the MD simulation indicated that the G*4-C17 pair moves continuously along a trajectory roughly perpendicular to the local helix axis, with retention of all three Watson-Crick hydrogen bonds. The simulation indicated the reversible formation of a hydrogen bond between the 5' oriented NH(3) ligand of platinum and the C3pG*4 phosphate group, in accord with our former prediction. Furthermore, the simulation has disclosed previously undetected BI <=> BII transitions at the G*5pA6 and A6pT7 steps, connected to formation/rupture of a hydrogen bond between the 3' oriented NH(3) ligand of platinum and the N7 atom of A6. All these conformational equilibria affect the form of the minor groove and increase the conformational flexibility at the platination site, and are thus likely to facilitate recognition by cellular proteins.

Platination of the (T2G4)4 telomeric sequence: a structural and cross-linking study

The telomeric sequence (T(2)G(4))(4) was platinated in aqueous solutions containing 50 mM LiClO(4), NaClO(4), or KClO(4). The identification of the guanines which reacted with [Pt(NH(3))(3)(H(2)O)](2+) revealed that the same type of folding exists in the presence of the three cations and that the latter determine the relative stabilities of the G-quadruplex structures in the order K(+) > Na(+) >> Li(+). The tri-ammine complex yielded ca. 40--90% of adducts, mono- and poly-platinated, bound to 4 guanines out of the 16 guanines in the sequence, in the decreasing amounts G9 > G15 >> G3 > G21. The formation of these adducts was interpreted with a G-quadruplex structure obtained by restrained molecular dynamics (rMD) simulations which confirms the schematic model proposed by Williamson et al. [(1989) Cell 59, 871--880]. The bifunctional complexes cis- and trans-[Pt(NH(3))(2)(H(2)O)(2)](2+) also first reacted with G9 and G15 and gave cross-linked adducts between two guanines, which did not exceed 5% each of the products formed. Both the cis and trans isomers formed a G3-G15 platinum chelate, and the second also formed bis-chelates at both ends of the G-quadruplex structure: G3-G15/G9-G21 and G3-G15/G9-G24. The rMD simulations showed that the cross-linking reactions by the trans complex can occur without disturbing the stacking of the three G-quartets.

Circular dichroism spectra of the individual rotamers of [Pt(N,N-dmen)(GpG)]+ (N,N-dmen=N,N-dimethylethylenediamine) indicate that the base-base oscillator coupling is not the main source of ellipticity in cis-PtG2L2 head-to-head complexes

The asymmetrical platinum complex [PtCl2(N,N-dmen)] (N,N-dmen = N,N-dimethylethylenediamine) reacts with the dinucleotide GpG to form two isomeric chelates of the formula [Pt(N,N-dmen)(GpG)]+ [9]. One of the isomers forms two stable rotamers separable by HPLC, whereas the other apparently prefers one single rotameric form. The favored conformations of these three forms were elucidated by means of molecular mechanics and dynamics techniques. In parallel, we have prepared the adduct, isolated the three rotamers, and recorded their solution circular dichroism (CD) spectra. For the first time we were thus able to correlate the CD features of individual rotamers of a cis-Pt(GpG) chelate with their structures. We show here that the two forms labeled in Inagaki's paper 1'e and 2e have the same right-handed helicoidal arrangement of the guanine bases but display different CD spectra in which the prominent bands have inverted signs. Thus, base-base interactions cannot be the principal cause of the CD of these compounds.