Gamma-radiation induced interstrand cross-links in PNA:DNA heteroduplexes

Parameterization of the miniPEG-Modified γPNA Backbone: Toward Induced γPNA Duplex Dissociation

γ-Modified peptide nucleic acids (γPNAs) serve as potential therapeutic agents against genetic diseases. Miniature poly(ethylene glycol) (miniPEG) has been reported to increase solubility and binding affinity toward genetic targets, yet details of γPNA structure and dynamics are not understood. Within our work, we parameterized missing torsional and electrostatic terms for the miniPEG substituent on the γ-carbon atom of the γPNA backbone in the CHARMM force field. Microsecond timescale molecular dynamics simulations were carried out on six miniPEG-modified γPNA duplexes from NMR structures (PDB ID: 2KVJ). Three NMR models for the γPNA duplex (PDB ID: 2KVJ) were simulated as a reference for structural and dynamic changes captured for the miniPEG-modified γPNA duplex. Principal component analysis performed on the γPNA backbone atoms identified a single isotropic conformational substate (CS) for the NMR simulations, whereas four anisotropic CSs were identified for the ensemble of miniPEG-modified γPNA simulations. The NMR structures were found to have a 23° helical bend toward the major groove, consistent with our simulated CS structure of 19.0°. However, a significant difference between simulated methyl- and miniPEG-modified γPNAs involved the opportunistic invasion of miniPEG through the minor and major groves. Specifically, hydrogen bond fractional analysis showed that the invasion was particularly prone to affect the second G-C base pair, reducing the Watson-Crick base pair hydrogen bond by 60% over the six simulations, whereas the A-T base pairs decreased by only 20%. Ultimately, the invasion led to base stack reshuffling, where the well-ordered base stacking was reduced to segmented nucleobase stacking interactions. Our 6 μs timescale simulations indicate that duplex dissociation suggests the onset toward γPNA single strands, consistent with the experimental observation of decreased aggregation. To complement the insight of miniPEG-modified γPNA structure and dynamics, the new miniPEG force field parameters allow for further exploration of such modified γPNA single strands as potential therapeutic agents against genetic diseases.

Improved Force Fields for Peptide Nucleic Acids with Optimized Backbone Torsion Parameters

Peptide nucleic acids are promising nucleic acid analogs for antisense therapies as they can form stable duplex and triplex structures with DNA and RNA. Computational studies of PNA-containing duplexes and triplexes are an important component for guiding their design, yet existing force fields have not been well validated and parametrized with modern computational capabilities. We present updated CHARMM and Amber force fields for PNA that greatly improve the stability of simulated PNA-containing duplexes and triplexes in comparison with experimental structures and allow such systems to be studied on microsecond time scales. The force field modifications focus on reparametrized PNA backbone torsion angles to match high-level quantum mechanics reference energies for a model compound. The microsecond simulations of PNA-PNA, PNA-DNA, PNA-RNA, and PNA-DNA-PNA complexes also allowed a comprehensive analysis of hydration and ion interactions with such systems.

Synthesis and Improved Cross-Linking Properties of C5-Modified Furan Bearing PNAs

Over the past decades, peptide nucleic acid/DNA (PNA:DNA) duplex stability has been improved via backbone modification, often achieved via introducing an amino acid side chain at the α- or γ-position in the PNA sequence. It was previously shown that interstrand cross-linking can further enhance the binding event. In this work, we combined both strategies to fine-tune PNA crosslinking towards single stranded DNA sequences using a furan oxidation-based crosslinking method; for this purpose, γ-l-lysine and γ-l-arginine furan-PNA monomers were synthesized and incorporated in PNA sequences via solid phase synthesis. It was shown that the l-lysine γ-modification had a beneficial effect on crosslink efficiency due to pre-organization of the PNA helix and a favorable electrostatic interaction between the positively-charged lysine and the negatively-charged DNA backbone. Moreover, the crosslink yield could be optimized by carefully choosing the type of furan PNA monomer. This work is the first to describe a selective and biocompatible furan crosslinking strategy for crosslinking of γ-modified PNA sequences towards single-stranded DNA.

Focus on PNA Flexibility and RNA Binding using Molecular Dynamics and Metadynamics

Peptide Nucleic Acids (PNAs) can efficiently target DNA or RNA acting as chemical tools for gene regulation. Their backbone modification and functionalization is often used to increase the affinity for a particular sequence improving selectivity. The understanding of the trading forces that lead the single strand PNA to bind the DNA or RNA sequence is preparatory for any further rational design, but a clear and unique description of this process is still not complete. In this paper we report further insights into this subject, by a computational investigation aiming at the characterization of the conformations of a single strand PNA and how these can be correlated to its capability in binding DNA/RNA. Employing Metadynamics we were able to better define conformational pre-organizations of the single strand PNA and γ-modified PNA otherwise unrevealed through classical molecular dynamics. Our simulations driven on backbone modified PNAs lead to the conclusion that this γ-functionalization affects the single strand preorganization and targeting properties to the DNA/RNA, in agreement with circular dichroism (CD) spectra obtained for this class of compounds. MD simulations on PNA:RNA dissociation and association mechanisms allowed to reveal the critical role of central bases and preorganization in the binding process.

Synthesis of peptide nucleic acids (PNA) with a crosslinking agent to RNA and effective inhibition of dicer

Peptide nucleic acids (PNAs) are structural mimics of nucleic acids that form stable hybrids with DNA and RNA. Due to these characteristics, PNAs are widely used as biochemical tools, for example, in antisense/antigene therapy. In this study, we have synthesized PNAs incorporating 2-amino-6-vinylpurine (AVP) for the covalent targeting of single-stranded DNA and RNA, and evaluated their reactivities for these targets. PNA containing AVP at the N-terminal position showed a high reactivity to uracil in RNA and thymine in DNA at the complementary site with AVP. In addition, the crosslinking reactions to pre-miR122 with PNA containing AVP increased the inhibition effect for the Dicer processing of pre-miR122 in vitro.

Furan-PNA: a mildly inducible irreversible interstrand crosslinking system targeting single and double stranded DNA

We here report on the design and synthesis of tailor-made furan-modified peptide nucleic acid (PNA) probes for covalent targeting of single stranded DNA through a crosslinking strategy. After introducing furan-containing building blocks into a PNA sequence, hybridization and furan-oxidation based crosslinking to DNA is investigated. The structure of the crosslinked products is characterized and preliminary investigations concerning the application of these systems to double stranded DNA are shown.

Modeling the interactions of the nucleotide excision repair UvrA(2) dimer with DNA

The UvrA protein initiates the DNA damage recognition process by the bacterial nucleotide excision repair (NER) system. Recently, crystallographic structures of holo-UvrA(2) dimers from two different microorganisms have been released (Protein Data Bank entries 2r6f , 2vf7 , and 2vf8 ). However, the details of the DNA binding by UvrA(2) and other peculiarities involved in the damage recognition process remain unknown. We have undertaken a molecular modeling approach to appraise the possible modes of DNA-UvrA(2) interaction using molecular docking and short-scale guided molecular dynamics [continuum field, constrained, and/or unrestricted simulated annealing (SA)], taking into account the three-dimensional location of a series of mutation-identified UvrA residues implicated in DNA binding. The molecular docking was based on the assumptions that the UvrA(2) dimer is preformed prior to DNA binding and that no major protein conformational rearrangements, except moderate domain reorientations, are required for binding of undamaged DNA. As a first approximation, DNA was treated as a rigid ligand. From the electrostatic relief of the ventral surface of UvrA(2), we initially identified three, noncollinear DNA binding paths. Each of the three resulting nucleoprotein complexes (C1, C2, and C3) was analyzed separately, including calculation of binding energies, the number and type of interaction residues (including mutated ones), and the predominant mode of translational and rotational motion of specific protein domains after SA to ensure improved DNA binding. The UvrA(2) dimer can accommodate DNA in all three orientations, albeit with different binding strengths. One of the UvrA(2)-DNA complexes (C1) fulfilled most of the requirements (high interaction energy, proximity of DNA to mutated residues, etc.) expected for a natural, high-affinity DNA binding site. This nucleoprotein presents a structural organization that is designed to clamp and bend double-stranded DNA. We examined the binding site in more detail by docking DNAs of significantly different (AT- vs CG-enriched) sequences and by submitting the complexes to DNA-unrestricted SA. It was found that in a manner independent of the DNA sequence and applied MD protocols, UvrA(2) favors binding of a bent and unwound undamaged DNA, with a kink positioned in the proximity of the Zn3 hairpins, anticollinearly aligned at the bottom of the ventral protein surface. It is further hypothesized that the Zn3 modules play an essential role in the damage recognition process and that the apparent existence of a family of DNA binding sites might be biologically relevant. Our data should prove to be useful in rational (structure-based) mutation studies.

Selection against spermatozoa with fragmented DNA after postovulatory mating depends on the type of damage

BACKGROUND

Before ovulation, sperm-oviduct interaction mechanisms may act as checkpoint for the selection of fertilizing spermatozoa in mammals. Postovulatory mating does not allow the sperm to attach to the oviduct, and spermatozoa may only undergo some selection processes during the transport through the female reproductive tract and/or during the zona pellucida (ZP) binding/penetration.

METHODS

We have induced DNA damage in spermatozoa by two treatments, (a) a scrotal heat treatment (42 degrees C, 30 min) and (b) irradiation with 137Cs gamma-rays (4 Gy, 1.25 Gy/min). The effects of the treatments were analyzed 21-25 days post heat stress or gamma-radiation. Postovulatory females mated either with treated or control males were sacrificed at Day 14 of pregnancy, and numbers of fetuses and resorptions were recorded.

RESULTS

Both treatments decreased significantly implantation rates however, the proportion of fetuses/resorptions was only reduced in those females mated to males exposed to radiation, indicating a selection favoring fertilization of sperm with unfragmented DNA on the heat treatment group. To determine if DNA integrity is one of the keys of spermatozoa selection after postovulatory mating, we analyzed sperm DNA fragmentation by COMET assay in: a) sperm recovered from mouse epididymides; b) sperm recovered from three different regions of female uterine horns after mating; and c) sperm attached to the ZP after in vitro fertilization (IVF). Similar results were found for control and both treatments, COMET values decreased significantly during the transit from the uterine section close to the uterotubal junction to the oviduct, and in the spermatozoa attached to ZP. However, fertilization by IVF and intracytoplasmatic sperm injection (ICSI) showed that during sperm ZP-penetration, a stringent selection against fragmented-DNA sperm is carried out when the damage was induced by heat stress, but not when DNA fragmentation was induced by radiation.

CONCLUSION

Our results indicate that in postovulatory mating there is a preliminary general selection mechanism against spermatozoa with low motility and fragmented-DNA during the transport through the female reproductive tract and in the ZP binding, but the ability of the ZP to prevent fertilization by fragmented-DNA spermatozoa is achieved during sperm-ZP penetration, and depends on the source of damage.