Conformational features of benzo‐homologated yDNA duplexes by molecular dynamics simulation

yDNA is a base‐modified nucleic acid duplex containing size‐expanded nucleobases. Base‐modified nucleic acids could expand the genetic alphabet and thereby enhance the functional potential of DNA. Unrestrained 100 ns MD simulations were performed in explicit solvent on the yDNA NMR sequence [5′(yA T yA yA T yA T T yA T)₂] and two modeled yDNA duplexes, [5′(yC yC G yC yC G G yC G G)₂] and [(yT5′ G yT A yC yG C yA yG T3′)?(yA5′ C T C yG C G yT A yC A3′)]. The force field parameters for the yDNA bases were derived in consistent with the well‐established AMBER force field. Our results show that DNA backbone can withstand the stretched size of the bases retaining the Watson‐Crick base pairing in the duplexes. The duplexes retained their double helical structure throughout the simulations accommodating the strain due to expanded bases in the backbone torsion angles, sugar pucker and helical parameters. The effect of the benzo‐expansion is clearly reflected in the extended C1′‐C1′ distances and enlarged groove widths. The size expanded base modification leads to reduction in base pair twist resulting in larger overlapping area between the stacked bases, enhancing inter and intra strand stacking interactions in yDNA in comparison with BDNA. This geometry could favour enhanced interactions with the groove binders and DNA binding proteins., 2016. © 2015 Wiley Periodicals, Inc. Biopolymers 105: 55–64, 2016     

Conformational geometry and vibrational frequencies of nucleic acid chains.

Normal coordinate analysis of diethyl phosphate has been made, which predicts all observed Raman frequencies in the range 170–1300 cm^?1. The force constants from this calculation have been transferred to a vibrational calculation for a simplified model of the backbone of nucleic acids, which also involves the ? O? PO₂^?? O phosphate group and the ? C₅′? C₄′? C₃′? linkage of the ribose. The coordinates of these atoms are those recently given by Arnott and Hukins, which place the ribose ring of B-DNA in a C₃′-exo conformation. This simple polymer model appears to be able to describe adequately the frequency-dependent changes observed in the Raman spectra arising from the backbone vibrations of nucleic acid in going from the B- to A-form. The symmetric ? O? P? O? diester stretch increases in frequency from about 787 cm^?1 in the B-form to 807 cm^?1 in the A-form. The increased frequency characteristic of the A-form is due to the combining of the diester stretch with vibrations involving the C₅′, C₄′, and C₃′ nuclei. The frequency of the symmetric ? O? P? O? diester stretch is shown to be very dependent on the conformation of the ribose ring, indicating that in polynucleotides the ribose ring takes on one of two rigid conformations: C₃′-endo for A-form or C₃′-exo for B-form and “disordered” polynucleotides. The calculation lends confirmation to the atomic coordinates of Arnott and Hukins since the use of other geometries with the same force constants failed to give results in agreement with experimental evidence. The calculations also demonstrate the lowering effect of hydration on the anionic PO stretching frequencies. Experimental results show that the 814-cm^?1 band observed in the spectra of 5′GMP gel arises from a different vibrational mode than that of the 814-cm^?1 band of A-DNA.     

Chemical shift assignments for human apurinic/apyrimidinic endonuclease 1

Apurinic/apyrimidinic endonuclease 1 (APE1 or Ref-1) is the major enzyme in mammals for processing abasic sites in DNA. These cytotoxic and mutagenic lesions arise via spontaneous rupture of the base-sugar bond or the removal of damaged bases by a DNA glycosylase. APE1 cleaves the DNA backbone 5′ to an abasic site, giving a 3′-OH primer for repair synthesis, and mediates other key repair activities. The DNA repair functions are essential for embryogenesis and cell viability. APE1-deficient cells are hypersensitive to DNA-damaging agents, and APE1 is considered an attractive target for inhibitors that could potentially enhance the efficacy of some anti-cancer agents. To enable an important new method for studying the structure, dynamics, catalytic mechanism, and inhibition of APE1, we assigned the chemical shifts (backbone and ¹³C^β) of APE1 residues 39-318. We also report a protocol for refolding APE1, which was essential for achieving complete exchange of backbone amide sites for the perdeuterated protein.     

Revealing structural peculiarities of homopurine GA repetition stuck by i-motif clip

Non-canonical forms of nucleic acids represent challenging objects for both structure-determination and investigation of their potential role in living systems. In this work, we uncover a structure adopted by GA repetition locked in a parallel homoduplex by an i-motif. A series of DNA oligonucleotides comprising GAGA segment and C₃ clip is analyzed by NMR and CD spectroscopies to understand the sequence–structure–stability relationships. We demonstrate how the relative position of the homopurine GAGA segment and the C₃ clip as well as single-base mutations (guanine deamination and cytosine methylation) affect base pairing arrangement of purines, i-motif topology and overall stability. We focus on oligonucleotides C₃GAGA and methylated GAGAC₃ exhibiting the highest stability and structural uniformity which allowed determination of high-resolution structures further analyzed by unbiased molecular dynamics simulation. We describe sequence-specific supramolecular interactions on the junction between homoduplex and i-motif blocks that contribute to the overall stability of the structures. The results show that the distinct structural motifs can not only coexist in the tight neighborhood within the same molecule but even mutually support their formation. Our findings are expected to have general validity and could serve as guides in future structure and stability investigations of nucleic acids.     

Conformational characteristics of polypeptides containing isolated L-proline residues with cis peptide bonds

The conformational characteristics of the peptide sequence X-l-Pro, where X  Gly or l-Ala and the peptide bond joining X and l-Pro is cis, are evaluated. Semi-empirical potential functions are used to estimate the contributions to the conformational energy made by the non-bonded van der Waals' and electrostatic interactions and the intrinsic torsional potentials about the NC^a and C^aC′ bonds. Rotations φ₁ and ψ₁ about the NC^a and C^aC′ bonds in residue X and rotation ψ₂ about the C^aC′ bond in l-Pro are permitted, while the angle of rotation φ₂ about the NC^a bond in l-Pro is fixed at 120 ° by the pyrrolidine ring. The presence of the cis peptide bond connecting X and l-Pro renders the backbone rotations φ₁, ψ₁ in X dependent upon the rotation ψ₂ about the C^aC′ bond in l-Pro. (Interdependence of rotations in neighboring residues joined by a cis peptide bond was previously observed in l-alanine oligomers.) The number of energetically allowed conformations for the Gly and l-Ala residues preceding a cis peptide bond l-Pro residue are found to be substantially reduced from those permitted when the peptide bond is trans or when l-Pro is replaced by an amino acid residue. On the other hand, ψ₂ = 100 to 160 ° (cis′) and 300 to 0 ° (trans′) are found to be the lowest energy conformations of the l-Pro residue irrespective of the cis or trans conformation of the X-l-Pro peptide bond.     

An ancient anion‐binding structural module in RNA and DNA helicases

RNA and DNA helicases manipulate or translocate along single strands of nucleic acids by grasping them using a conserved structural motif. We have examined the available crystal structures of helicases of the two principal superfamilies, SF1 and SF2, and observed that the most conserved interactions with the nucleic acid occur between the phosphosugar backbone of a trinucleotide and the three strand‐helix loops within a (β‐strand/α‐helix)₃ structural module. At the first and third loops is a conserved hydrogen‐bonded feature called a thr‐motif, often seen at α‐helical N‐termini, with the threonine as the N‐cap residue. These loops can be aligned with few insertions or deletions, and their main chain atoms are structurally congruent amongst the family members and between the two modules found as tandem pairs in all SF1 and SF2 proteins. The other highly conserved interactions with nucleic acid involve mainchain NH groups, often at the helical N‐termini, interacting with phosphate groups. We comment on how the sequence motifs that are commonly used to identify helicases map to locations on the module and discuss the implications of the conserved orientation of nucleic acid on the surface of the module for directional stepping along DNA or RNA. Proteins 2010. © 2010 Wiley‐Liss, Inc.     

Inhibition of Topoisomerases by Fatty Acids

The inhibitory effects of various fatty acids on topoisomerases were examined, and their structure-activity relationships and mechanism of action were studied. Saturated fatty acids (C_6:0 to C_22:0) did not inhibit topoisomerase I, but cis-unsaturated fatty acids (C_16:1 to C_22:1) with one double bond showed strong inhibition of the enzyme. The inhibitory potency depended on the carbon chain length and the position of the double bond in the fatty acid molecule. The trans-isomer, methyl ester and hydroxyl derivative of oleic acid had no or little inhibitory effect on topoisomerases I and II. Among the compounds studied petroselinic acid and vaccenic acid (C_18:1) with a cis-double bond were the potent inhibitors. Petroselinic acid was a topoisomerase inhibitor of the cleavable complex-nonforming type and acted directly on the enzyme molecule in a noncompetitive manner without DNA intercalation.     

Germanium‐ and Silicon‐Substituted Donor–Acceptor Type Copolymers: Effect of the Bridging Heteroatom on Molecular Packing and Photovoltaic Device Performance

The effects of heteroatom substitution from a silicon atom to a germanium atom in donor‐acceptor type low band gap copolymers, poly[(4,4′‐bis(2‐ethylhexyl)dithieno[3,2‐b:2′,3′‐d]silole)‐2,6‐diyl‐alt‐(2,1,3‐benzothiadiazole)‐4,7‐diyl] (PSiBTBT) and poly[(4,4′‐bis(2‐ethylhexyl)dithieno[3,2‐b:2′,3′‐d]germole)‐2,6‐diyl‐alt‐(2,1,3‐benzothiadiazole)‐4,7‐diyl] (PGeBTBT), are studied. The optoelectronic and charge transport properties of these polymers are investigated with a particular focus on their use for organic photovoltaic (OPV) devices in blends with phenyl‐C₇₀‐butyric acid methyl ester (PC₇₀BM). It is found that the longer C‐Ge bond length, in comparison to C‐Si, modifies the molecular conformation and leads to a more planar chain conformation in PGeBTBT than PSiBTBT. This increase in molecular planarity leads to enhanced crystallinity and an increased preference for a face‐on backbone orientation, thus leading to higher charge carrier mobility in the diode configuration. These results provide important insight into the impact of the heavy atom substitution on the molecular packing and device performance of polymers based on the poly[2,6‐(4,4‐bis‐(2‐ethylhexyl)‐4H‐cyclopenta[2,1‐b;3,4‐b]‐dithiophene)‐alt‐4,7‐(2,1,3‐benzothiadiazole) (PCPDTBT) backbone.     

High resolution nuclear magnetic resonance spectroscopy of glycosphingolipids. I: 360 MHz 1H and 90.5 MHz 13C NMR analysis of galactosylceramides

The high resolution ¹H and ¹³C nuclear magnetic resonance (NMR) spectra of galactosylceramides containing n-fatty acids and α-hydroxy fatty acids were recorded in dimethylsulfoxide solution with and without addition of D₂O. From the coupling constants of the sugar ring protons, a ⁴C₁ conformation can be deduced. In contrast to the conformation in aqueous solution, the C⁶ hydroxymethylene group is freely rotating around the C⁶C⁵ bond. In the ceramide residue all signals produced by protons linked to carbons bearing electronegative substituents could be attributed. The large difference in coupling constants of the methylene protons of C^1′ to the C^2′ methine proton of the sphingosine indicates a restricted rotation around the C^1′C^2′ bond. The assignments of the hydroxy and amino protons follow from the decoupling of the corresponding methine protons.     

Prediction of Xaa-Pro peptide bond conformation from sequence and chemical shifts

We present a program, named Promega, to predict the Xaa-Pro peptide bond conformation on the basis of backbone chemical shifts and the amino acid sequence. Using a chemical shift database of proteins of known structure together with the PDB-extracted amino acid preference of cis Xaa-Pro peptide bonds, a cis/trans probability score is calculated from the backbone and ¹³C^β chemical shifts of the proline and its neighboring residues. For an arbitrary number of input chemical shifts, which may include Pro-¹³C^γ, Promega calculates the statistical probability that a Xaa-Pro peptide bond is cis. Besides its potential as a validation tool, Promega is particularly useful for studies of larger proteins where Pro-¹³C^γ assignments can be challenging, and for on-going efforts to determine protein structures exclusively on the basis of backbone and ¹³C^β chemical shifts.     

Shape readout of AT‐rich DNA by carbohydrates

Gene expression can be altered by small molecules that target DNA; sequence as well as shape selectivities are both extremely important for DNA recognition by intercalating and groove‐binding ligands. We have characterized a carbohydrate scaffold (1) exhibiting DNA “shape readout” properties. Thermodynamic studies with 1 and model duplex DNAs demonstrate the molecule's high affinity and selectivity towards B* form (continuous AT‐rich) DNA. Isothermal titration calorimetry (ITC), circular dichroism (CD) titration, ultraviolet (UV) thermal denaturation, and Differential Scanning Calorimetry were used to characterize the binding of 1 with a B* form AT‐rich DNA duplex d[5′‐G₂A₆T₆C₂‐3′]. The binding constant was determined using ITC at various temperatures, salt concentrations, and pH. ITC titrations were fit using a two‐binding site model. The first binding event was shown to have a 1:1 binding stoichiometry and was predominantly entropy‐driven with a binding constant of approximately 10⁸ M^?1. ITC‐derived binding enthalpies were used to obtain the binding‐induced change in heat capacity (ΔC_p) of ?225 ± 19 cal/mol·K. The ionic strength dependence of the binding constant indicated a significant electrolytic contribution in ligand:DNA binding, with approximately four to five ion pairs involved in binding. Ligand 1 displayed a significantly higher affinity towards AT‐tract DNA over sequences containing GC inserts, and binding experiments revealed the order of binding affinity for 1 with DNA duplexes: contiguous B* form AT‐rich DNA (d[5′‐G₂A₆T₆C₂‐3′]) >B form alternate AT‐rich DNA (d[5′‐G₂(AT)₆C_2‐3′]) > A form GC‐rich DNA (d[5′‐A₂G₆C₆T₂‐3′]), demonstrating the preference of ligand 1 for B* form DNA. © 2013 Wiley Periodicals, Inc. Biopolymers 101: 720–732, 2014.     

Multinuclear magnetic resonance studies of monomers and dimers containing 2'-fluoro-2'-deoxyadenosine

A series of 2′-fluorinated adenosine compounds, dAfl, dAflp, pdAfl, dAfl-A, A-dAfl, and dAfl-dAfl, have been investigated by nmr spectroscopies. The ¹H-, ¹⁹F-, and ³¹P-nmr data provide structural information from different parts of these moleucles. The pK_a of the phosphate group of these two 2′-fluoro-2′-deoxyadenosine monophosphates was found to be the same as that of hte parent adenosine monophosphate. As for the pentose conformation, the ³E population is greatly increased as a result of the fluorine substitution at the C_2′ position. However, the populations of conformers of gg (C_4′-C_5′) and g′g′ (C_5′-O_5′) and the average angle ?′(C_3′-O_3′) of the 2′-fluoro compounds remain unchanged as compared to the natural riboadenosine monomer and dimer (A-A). Thefefore, the backbone conformation of the 2′-fluoro-2′-deoxy-adenosine, its monophosphates and dimers, resembles that of RNA. The extent of base-base overlapping in these 2′-fluoro-2′-deoxy-adenosine-containing dimers is also found to be similar to or even greater than A-A. Thus, the conformations of these compounds can be considered as those in the RNA family. These fluorocompounds also serve as models for a careful study on the ¹⁹F-nmr in nucleic acid. The ¹⁹F chemical-shift values are sensitive to the environment of the fluorine atom such as ionic structure of the neighboring group(s) (phosphate of base), solvation, and ring-ruccent anisotropic effect from the base(s). Qualitatively, the change of the ¹⁹F chemical-shift values (up to 2 ppm) is much larger than that of ¹H-nmr (up to 0.5 ppm) in the dimers. Using dAfl·poly(U), poly(dAfl)·poly(dAfl), and poly(dAfl)·poly(U) helix–coil transition as model systems, the linewidth of ¹⁹F in dAfl- residues reflects effectively the mobility of the unit in the nucleic acid complex as calibrated by uv data and by ¹H-nmr. Therefore, application of ¹⁹F-nmr spectroscopy on fluorine-substituted nucleic acid can also be used to detect nucleic acid-nucleic acid interaction in complicated systems.     

Probing mechanism of metal catalyzed hydrolysis of Thymidylyl (3′-O, 5′-S) thymidine phosphodiester derivatives

Hydrolysis of nucleic acids is of fundamental importance in biological sciences. Kinetic and theoretical studies on different substrates wherein the phosphodiester bond combined with alkyl or aryl groups and sugar moiety have been the focus of attention in recent literature. The present work focuses on understanding the mechanism and energetics of alkali metal (Li, Na, and K) catalyzed hydrolysis of phosphodiester bond in modeled substrates including Thymidylyl (3′-O, 5′-S) thymidine phosphodiester (Tp-ST) (1), 3′-Thymidylyl (1-trifluoroethyl) phosphodiester (Tp-OCH₂CF₃) (2), 3′-Thymidylyl (o-cholorophenyl) phosphodiester (Tp-OPh(o-Cl)) (3) and 3′-Thymidylyl(p-nitrophenyl) phosphodiester (Tp-OPh(p-NO₂)) (4) employing density functional theory. Theoretical calculations reveal that the reaction follows a single-step (A_ND_N) mechanism where nucleophile attack and leaving group departure take place simultaneously. Activation barrier for potassium catalyzed Tp-ST hydrolysis (12.0 kcal mol^?1) has been nearly twice as large compared to that for hydrolysis incorporating lithium or sodium. Effect of solvent (water) on activation energies has further been analyzed by adding a water molecule to each metal ion of the substrate. It has been shown that activation barrier of phosphodiester hydrolysis correlates well with basicity of leaving group.

2′,4′-BNA/LNA aptamers: CE-SELEX using a DNA-based library of full-length 2′-O,4′-C-methylene-bridged/linked bicyclic ribonucleotides

DNA-based aptamers that contain 2′-O,4′-C-methylene-bridged/linked bicyclic ribonucleotides (B/L nucleotides) over the entire length were successfully obtained using a capillary electrophoresis systematic evolution of ligands by exponential enrichment (CE-SELEX) method. A modified DNA library was prepared with an enzyme mix of KOD Dash and KOD mutant DNA polymerases. Forty 2′-O,4′-C-methylene bridged/locked nucleic acid (2′,4′-BNA/LNA) aptamers were isolated from an enriched pool and classified into six groups according to their sequence. 2′,4′-BNA/LNA aptamers of groups V and VI bound human thrombin with K_d values in the range of several 10 nanomolar levels.     

Conformational studies on polynucleotide chains. IV. Backbone structure of ribonucleic acids.

In preceding papers the energies associated with the internal rotations in the sugar–phosphate–sugar complex were described with an analytical potential consisting of a Lennard-Jones 6–12 term and an intrinsic torsional term and representing the best fit to a large number of energies computed with a quantum mechanical ab initio technique. The complex considered there (of 37 atoms and with the chemical formula C₁₀H₁₈O₈P) is repesentative of deoxyribonucleic acids. In this paper we apply our potential to evaluating the intramolecular energies of the 39-atom complex C₁₀H₁₈O₁₀P, representative of the ribonucleic acids. The potential energies for the internal rotations (considered independent from one another) and the energy maps for rotations about consecutive bonds of the backbone chain are critically compared, both with those obtained for the deoxy system and with those obtained from different theoretical approaches as available from literature. It is shown that, at least for certain combinations of the internal rotation angles, the choice of the starting geometry for the sugarphosphate–sugar molecule (bond lengths and valence angles) strongly affects the value of the computed energy. If a proper geometry is used, very low energies are predicted by our potential in correspondence of the sets of torsional angles found in various RNAs by x-ray crystallography.     

Molecular simulation of cyclohexanyl nucleic acid (CNA) duplexes with CNA,DNA and RNA and CNA triloop and tetraloop hairpin structures

As part of a selection strategy for artificial nucleic acids (XNA) (to be considered as potential new information systems in vivo), we have carried out a modelling study on cyclohexanyl nucleic acids (CNA) duplexes and hairpins. CNA may form a duplex as well as hairpin structures, having the carbocyclic nucleosides in the ⁴C₁ conformation (with equatorial basis). The geometry of ds CNA is close to that of a HNA:RNA duplex. We demonstrated that CNA triphosphates function as a substrate for polymerases. Modelling experiments indicate that the monomers are probably presented to the polymerase in the ¹C₄ conformation.     

What stabilizes the 314‐helix in β3‐peptides? A conformational analysis using molecular simulation

β‐Peptides are analogs of natural α‐peptides and form a variety of remarkably stable structures. Having an additional carbon atom in the backbone of each residue, their folded conformation is not only influenced by the side‐chain sequence but also and foremost by their substitution pattern. The precise mechanism by which the side chains interact with the backbone is, however, hitherto not completely known. To unravel the various effects by which the side chains influence the backbone conformation, we quantify to which extent the dihedral angles of a β³‐substited peptide with an additional methyl group on the central C_α‐atom can be regarded as independent degrees of freedom and analyze the distributions of these dihedral angles. We also selectively capture the steric effect of substituents on the C_α‐ and C_β‐atoms of the central residue by alchemically changing them into dummy atoms, which have no nonbonded interactions. We find that the folded state of the β³‐peptide is primarily stabilized by a steric exclusion of large parts of the unfolded state (entropic effect) and only subsequently by mutual dependence of the ψ‐dihedral angles (enthalpic effect). The folded state of β‐peptides is stabilized by a different mechanism than that of α‐peptides. Proteins 2010. © 2010 Wiley‐Liss, Inc.     

Backbone Conformation in Nucleic Acids: An Analysis of Local Helicity through Heminucleotide Scheme and a Proposal for a Unified Conformational Plot

Abstract

A relationship has been established to express the local helicity of a polynucleotide backbone directly in terms of the virtual bonds spanning the conformationally equivalent heminucleotide repeats, with a view to provide a better understanding of the cumulative effects of all the chemical bond rotational variations on local helicity. Using this, an analysis made with a few oligodeoxynucleotide crystal structures clearly brings forth that it is the concerted movements manifested in the near neighbour correlations between the pair of chemical bonds C4′—C5′ and P—05′ and C4′-C3′ and P-03′ of the 5′ and 3′ heminucleotides respectively that are primarily responsible for the observed non-uniform helical twists both in A and B type helical backbones. That these need not be restricted to oligodeoxynucleotides but may be a feature of oligoribonucleotides backbone also is shown from an analysis of helical segments of yeast tRNA^Phe. A proposal of a unified or a grand two dimensional conformational plot which would help visualise succinctly the overall effect of the variations in all the repeating six chemical bonds of a polynucleotide backbone is made. Apart from considerable simplification, the plot affords identification on it regions characteristic of helical, and loop and bend conformations of nucleic acid backbone chain.