Molecular dynamics structures of peptide nucleic acid x DNA hybrid in the wild-type and mutated alleles of Ki-ras proto-oncogene--stereochemical rationale for the low affinity of PNA in the presence of an AC mismatch

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.

Nanoscale structures and mechanics of peptide nucleic acids

Peptide nucleic acids (PNAs) are charge-neutral polyamide oligomers having extremely favorable thermal stability and high affinity to cell membranes when coupled with cationic cell-penetrating peptides (CPPs), as well as the encouraging antisense and antigene activity in cell-free systems. The study of the mechanical properties of short PNA molecules is rare both in experiments and theoretical calculations. Here, we studied the microscopic structures and elastic properties; namely, persistence length, stretch modulus, twist-stretch coupling, and structural crookedness of double-stranded PNA (dsPNA) and their hybrid derivatives using all-atom MD simulation and compared them with those of double-stranded DNA (dsDNA) and double-stranded RNA (dsRNA). The stretch modulus of the dsPNA is found to be ∼160 pN, an order of magnitude lower than that of dsDNA and smaller than dsRNA, respectively. Similarly, the persistence length of dsPNA is found to be ∼35 nm, significantly smaller than those of dsDNA and dsRNA. The PNA-DNA and PNA-RNA hybrid duplexes have elastic properties lying between that of dsPNA and dsDNA/dsRNA. We argue that the neutral backbones of the PNA make it less stiff than dsDNA and dsRNA molecules. Measurement of structural crookedness and principal component analysis additionally support the bending flexibility of dsPNA. Detailed analysis of the helical-rise coupled to helical-twist indicates that the PNA-DNA hybrid over-winds like dsDNA, while PNA-PNA and PNA-RNA unwind like dsRNA upon stretching. Because of the highly flexible nature of PNA, it can bind other biomolecules by adopting a wide range of conformations and is believed to be crucial for future nanobiotechnology research studies.

Conformational distortions induced by periodically recurring A…A in d(CAG).d(CAG) provide stereochemical rationale for the trapping of MSH2.MSH3 in polyQ disorders

CAG repeat instability causes a number of neurodegenerative disorders. The unusual hairpin stem structure formed by the CAG repeats in DNA traps the human mismatch repair MSH2.MSH3 (Mutsβ) complex. To understand the mechanism behind the abnormal binding of Mutsβ with the imperfect hairpin stem structure formed by CAG repeats, molecular dynamics simulations have been carried out for Mutsβ-d(CAG)₂(CAG)(CAG)₂.d(CTG)₂(CAG)(CTG)₂ (1 A…A mismatch) and Mutsβ-d(CAG)₅.d(CAG)₅ (5 mismatches, wherein, A…A occurs periodically) complexes. The interaction of MSH3 residue Tyr₂₄₅ at the minor groove side of A…A, an essential interaction responsible for the recognition by Mutsβ, are retained in both the cases. Nevertheless, the periodic unwinding caused by the nonisostericity of A…A with the flanking canonical base pairs in d(CAG)₅.d(CAG)₅ distorts the regular B-form geometry. Such an unwinding exposes one of the A…A mismatches (that interacts with Tyr₂₄₅) at the major groove side and also facilitates the on and off hydrogen bonding interaction with Lys₅₄₆ sidechain (MSH2-domain-IV). In contrast, kinking of the DNA towards the major groove in Mutsβ-d(CAG)₂(CAG)(CAG)₂.d(CTG)₂(CAG)(CTG)₂ doesn’t facilitate such an exposure of the bases at the major groove. Further, the unwinding of the helix in d(CAG)₅.d(CAG)₅ enhances the tighter binding between MSH2-domain-I and d(CAG)₅.d(CAG)₅ at the major groove side as well as between MSH3-domain-I and MSH3-domain-IV. Markedly, such enhanced interactions are absent in Mutsβ-d(CAG)₂(CAG)(CAG)₂.d(CTG)₂(CAG)(CTG)₂ that has a single A…A mismatch. Thus, the above-mentioned enhancement in intra- and inter- molecular interactions in Mutsβ-d(CAG)₅.d(CAG)₅ provide the stereochemical rationale for the trapping of Mutsβ in CAG repeat expansion disorders.

Spontaneous and frequent conformational dynamics induced by A…A mismatch in d(CAA)·d(TAG) duplex

Base pair mismatches in DNA can erroneously be incorporated during replication, recombination, etc. Here, the influence of A…A mismatch in the context of 5′CAA·5′TAG sequence is explored using molecular dynamics (MD) simulation, umbrella sampling MD, circular dichroism (CD), microscale thermophoresis (MST) and NMR techniques. MD simulations reveal that the A…A mismatch experiences several transient events such as base flipping, base extrusion, etc. facilitating B–Z junction formation. A…A mismatch may assume such conformational transitions to circumvent the effect of nonisostericity with the flanking canonical base pairs so as to get accommodated in the DNA. CD and 1D proton NMR experiments further reveal that the extent of B–Z junction increases when the number of A…A mismatch in d(CAA)·d(T(A/T)G) increases (1–5). CD titration studies of d(CAA)·d(TAG)_n=5 with the hZα_ADAR1 show the passive binding between the two, wherein, the binding of protein commences with B–Z junction recognition. Umbrella sampling simulation indicates that the mismatch samples anti…+ syn/+ syn…anti, anti…anti & + syn…+ syn glycosyl conformations. The concomitant spontaneous transitions are: a variety of hydrogen bonding patterns, stacking and minor or major groove extrahelical movements (with and without the engagement of hydrogen bonds) involving the mismatch adenines. These transitions frequently happen in anti…anti conformational region compared with the other three regions as revealed from the lifetime of these states. Further, 2D-NOESY experiments indicate that the number of cross-peaks diminishes with the increasing number of A…A mismatches implicating its dynamic nature. The spontaneous extrahelical movement seen in A…A mismatch may be a key pre-trapping event in the mismatch repair due to the accessibility of the base(s) to the sophisticated mismatch repair machinery.

Sequence dependent influence of an A…A mismatch in a DNA duplex: An insight into the recognition by hZαADAR1 protein

Base pair mismatches can erroneously be incorporated in the DNA. An adenine pairing with another adenine is one of the eight possible mismatches. The atomistic insights about the structure and dynamics of an A…A mismatch in a DNA (unbound form) is not yet accessible to any experimental technique. Earlier molecular dynamics (MD) simulations have shown that A…A mismatch in the midst of 5'CAG/3'GAC, 5'GAC/3'CAG and 5'CAA/3'GAT (underline represents the mismatch) are highly dynamic in nature. By employing MD simulation, the influence of an A…A mismatch in the midst of 5'GAA/3'CAT, 5'GAG/3'CAC, 5'AAC/3'TAG, 5'AAG/3'TAC, 5'TAA/3'AAT, 5'TAT/3'AAA and 5'AAT/3'TAA sequences have been investigated here. The results indicate that irrespective of the flanking sequences, the mismatch samples a variety of transient conformations, including a B-Z junction. Further, circular dichroism studies have been carried out to explore the ability of these sequences to bind with hZα_ADAR1 which specifically recognizes B-Z junction/Z-DNA. The results indicate that hZα_ADAR1 could not lead to a complete B to Z transition in the above sequences. Notably, a complete transition to Z-form has been reported earlier for 5'GAC/3'CAG upon titrating with hZα_ADAR1. Intriguingly, 5'AAC/3'TAG, 5'AAG/3'TAC and 5'GAA/3'CAT exhibit a B-Z junction formation rather than a complete transition to Z-form, similar to the situation of 5'CAA/3'GAT. These indicate that although A…A mismatch could induce a local B-Z junction transiently, hZα_ADAR1 requires the presence of a G…C/C…G base pair adjacent to the A…A mismatch for the binding. Additionally, the extent of B-Z junction has enhanced upon binding with hZα_ADAR1 in the presence of the A…A mismatch (specifically when CG, CA, AC, GA and AG steps occur), but not in the presence of the canonical base pairs. These confirm the inclination of A…A mismatch towards the B-Z junction.

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.

Predictive Model for the Sequence-Dependent Fluorogenic Response of Forced-Intercalation Peptide Nucleic Acid

The forced-intercalation peptide nucleic acid (FIT-PNA) concept, introduced by Seitz and co-workers, is based on replacing a nucleobase of the PNA sequence with a cyanine dye (such as thiazole orange). The cyanine dye is thus a surrogate base that is forced to intercalate in the duplex (e.g.

, PNA

DNA). This allows single-mismatch sensitivity as the introduction of a mismatch in the vicinity of the dye increases freedom of motion and leads to a significant depletion of its fluorescence because of the free rotation of the monomethine bond separating the two π-systems of the cyanine dye. Herein, we designed and synthesized six FIT-PNA probes, featuring bisquinoline (BisQ), a red-emitting cyanine dye recently developed in our laboratory for FIT-PNAs. By following PNA-DNA duplex fluorescence, we found new sequence-based factors governing the fluorescence response to the mismatched FIT-PNA:DNA duplex. Fluorogenic properties are correlated with the π-stacking energy of three distinctive base pair steps (BPSs) in the PNA:DNA duplex. The first two are the two BPSs opposite BisQ, whereas the third is the BPS of the mismatch position, which presumably becomes unstacked due to the mismatch. We suggest a predictive model for FIT-PNA single-mismatch detection mechanism, a model that can be used in future research to improve FIT-PNA design.

Selective Preference of Parallel DNA Triplexes Is Due to the Disruption of Hoogsteen Hydrogen Bonds Caused by the Severe Nonisostericity between the G*GC and T*AT Triplets

Implications of DNA, RNA and RNA.DNA hybrid triplexes in diverse biological functions, diseases and therapeutic applications call for a thorough understanding of their structure-function relationships. Despite exhaustive studies mechanistic rationale for the discriminatory preference of parallel DNA triplexes with G*GC & T*AT triplets still remains elusive. Here, we show that the highest nonisostericity between the G*GC & T*AT triplets imposes extensive stereochemical rearrangements contributing to context dependent triplex destabilisation through selective disruption of Hoogsteen scheme of hydrogen bonds. MD simulations of nineteen DNA triplexes with an assortment of sequence milieu reveal for the first time fresh insights into the nature and extent of destabilization from a single (non-overlapping), double (overlapping) and multiple pairs of nonisosteric base triplets (NIBTs). It is found that a solitary pair of NIBTs, feasible either at a G*GC/T*AT or T*AT/G*GC triplex junction, does not impinge significantly on triplex stability. But two overlapping pairs of NIBTs resulting from either a T*AT or a G*GC interruption disrupt Hoogsteen pair to a noncanonical mismatch destabilizing the triplex by ~10 to 14 kcal/mol, implying that their frequent incidence in multiples, especially, in short sequences could even hinder triplex formation. The results provide (i) an unambiguous and generalised mechanistic rationale for the discriminatory trait of parallel triplexes, including those studied experimentally (ii) clarity for the prevalence of antiparallel triplexes and (iii) comprehensive perspectives on the sequence dependent influence of nonisosteric base triplets useful in the rational design of TFO's against potential triplex target sites.

Nucleotides containing variously modified sugars: energetics, structure, and mechanical properties

The influence of various sugar residue modifications on intrinsic energetic, conformational, and mechanical properties of 2'-deoxyribonucleotide-5'-monophosphates (dNs) was comprehensively investigated using modern quantum chemical approaches. In total, fourteen sugar modifications, including double bonds and heteroatoms (S and N) inside the sugar ring, as well as fluorination in various positions, were analyzed. Among hundreds of possible conformational states of dNs, only two - AI and BI, corresponding to the most biologically significant forms of a double-helical DNA, were considered for each dN. It was established that the most of the studied modifications tend to strongly stabilize either AI or BI conformation of dNs both in the gas phase and in aqueous solution (modelled by implicit solvent models). Therefore, some of these modifications can be used as a tool for reducing structural polymorphism of nucleic acids in solution as well as for designing oligonucleotides with specific structural features. The evaluation of relaxed force constants (RFC) for glycosidic bonds suggests that the majority of the studied modifications of the sugar residue yield increased strengths of glycosidic bonds in dNs, and can therefore be used for designing modified nucleic acids with an increased resistance to abasic lesions. The most significant reinforcement of the glycosidic bond occurs in dNs containing the CF2 group instead of the O4' oxygen and the fluorine atom at the 2'-α-position. The calculation of the RFC and vibrational root-mean-square (VRMS) deviations for conformational degrees of freedom revealed a strong dependence between mechanical properties of dNs and their energetic characteristics. In particular, electronic energies of AI and BI conformers of dNs calculated in vacuo are closely connected with the values of relaxed force constants (RFC) for the δ angle: the higher RFC(δ) values correspond to more energetically favorable conformers.

Evaluation of the Gibbs Free Energy Changes and Melting Temperatures of DNA/DNA Duplexes Using Hybridization Enthalpy Calculated by Molecular Dynamics Simulation

A molecular dynamics simulation approach was applied for the prediction of the thermal stability of oligonucleotide duplexes. It was shown that the enthalpy of the DNA/DNA complex formation could be calculated using this approach. We have studied the influence of various simulation parameters on the secondary structure and the hybridization enthalpy value of Dickerson-Drew dodecamer. The optimal simulation parameters for the most reliable prediction of the enthalpy values were determined. The thermodynamic parameters (enthalpy and entropy changes) of a duplex formation were obtained experimentally for 305 oligonucleotides of various lengths and GC-content. The resulting database was studied with molecular dynamics (MD) simulation using the optimized simulation parameters. Gibbs free energy changes and the melting temperatures were evaluated using the experimental correlation between enthalpy and entropy changes of the duplex formation and the enthalpy values calculated by the MD simulation. The average errors in the predictions of enthalpy, the Gibbs free energy change, and the melting temperature of oligonucleotide complexes were 11%, 10%, and 4.4 °C, respectively. We have shown that the molecular dynamics simulation gives a possibility to calculate the thermal stability of native DNA/DNA complexes a priori with an unexpectedly high accuracy.

Twisting right to left: A…A mismatch in a CAG trinucleotide repeat overexpansion provokes left-handed Z-DNA conformation

Conformational polymorphism of DNA is a major causative factor behind several incurable trinucleotide repeat expansion disorders that arise from overexpansion of trinucleotide repeats located in coding/non-coding regions of specific genes. Hairpin DNA structures that are formed due to overexpansion of CAG repeat lead to Huntington’s disorder and spinocerebellar ataxias. Nonetheless, DNA hairpin stem structure that generally embraces B-form with canonical base pairs is poorly understood in the context of periodic noncanonical A…A mismatch as found in CAG repeat overexpansion. Molecular dynamics simulations on DNA hairpin stems containing A…A mismatches in a CAG repeat overexpansion show that A…A dictates local Z-form irrespective of starting glycosyl conformation, in sharp contrast to canonical DNA duplex. Transition from B-to-Z is due to the mechanistic effect that originates from its pronounced nonisostericity with flanking canonical base pairs facilitated by base extrusion, backbone and/or base flipping. Based on these structural insights we envisage that such an unusual DNA structure of the CAG hairpin stem may have a role in disease pathogenesis. As this is the first study that delineates the influence of a single A…A mismatch in reversing DNA helicity, it would further have an impact on understanding DNA mismatch repair.

When a set of 3 nucleotides in a DNA sequence repeats beyond a certain number, it leads to incurable neurological or neuromuscular disorders. Such DNA sequences tend to form unusual DNA structures comprising of base pairing schemes different from the canonical A…T/G…C base pairs. Influence of such unusual base pairing on the overall 3-dimensional structure of DNA and its impact on the pathogenesis of disorder is not well understood. CAG repeat overexpansion that leads to Huntington’s disorder and several spinocerebellar ataxias forms noncanonical A…A base pair in between canonical C…G and G…C base pairs. However, no detailed structural information is available on the influence of an A…A mismatch on a DNA structure under any sequence context. Here, we have shown for the first time that A…A base pairing in a CAG repeat provokes the formation of left-handed Z-DNA due to the pronounced structural dissimilarity of A…A base pair with G…C base pair, leading to periodic B-Z junction. Thus, these results suggest that formation of periodic B-Z junction may be one of the molecular bases for CAG repeat instability.

Molecular dynamics simulations of PNA-PNA and PNA-DNA duplexes by the use of new parameters implemented in the GROMACS package: a conformational and dynamics study

Peptide Nucleic Acids (PNAs) still represent a growing research area thanks to their potential applications in many fields of science from chemistry and biology to medicine. In these years, structural investigations by means of either experimental or computational techniques have proved to be very useful for the understanding of the structural organization and the binding properties of PNA. In this context, we here report an all-atoms Molecular Dynamics (MD) study of a PNA-PNA duplex and a PNA-DNA hetero-duplex with the well known GROMACS simulation package, by using new force field parameters properly derived for PNA molecules. The good agreement of our results with the crystallographic and NMR data, available for both the systems under investigation, confirms the validity of our approach. Moreover, our simulations reveal new interesting features related to the conformational-dynamic behavior of the studied systems, thus demonstrating the ability of MD simulations to gain insights into the dynamic properties of biologically relevant systems. This force field parametrization represents a good starting point for the implementation of a computational platform, based on the GROMACS package, useful for the rational design of modified PNA molecules with improved conformational features for selective binding toward DNA or RNA.

Effects of hypoxanthine substitution in peptide nucleic acids targeting KRAS2 oncogenic mRNA molecules: theory and experiment

Genetic disorders can arise from single base substitutions in a single gene. A single base substitution for wild type guanine in the twelfth codon of KRAS2 mRNA occurs frequently to initiate lung, pancreatic, and colon cancer. We have observed single base mismatch specificity in radioimaging of mutant KRAS2 mRNA in tumors in mice by in vivo hybridization with radiolabeled peptide nucleic acid (PNA) dodecamers. We hypothesized that multimutant specificity could be achieved with a PNA dodecamer incorporating hypoxanthine, which can form Watson-Crick base pairs with adenine, cytosine, thymine, and uracil. Using molecular dynamics simulations and free energy calculations, we show that hypoxanthine substitutions in PNAs are tolerated in KRAS2 RNA:PNA duplexes where wild type guanine is replaced by mutant uracil or adenine in RNA. To validate our predictions, we synthesized PNA dodecamers with hypoxanthine, and then measured the thermal stability of RNA:PNA duplexes. Circular dichroism thermal melting results showed that hypoxanthine-containing PNAs are more stable in duplexes where hypoxanthine-adenine and hypoxanthine-uracil base pairs are formed than single mismatch duplexes or duplexes containing hypoxanthine-guanine opposition.

Molecular dynamics of potential rRNA binders: single-stranded nucleic acids and some analogues

By hindering or "silencing" protein translation in vivo, antisense nucleic acid analogues that hybridize to bacterial rRNA could serve as a promising class of antibacterial compounds. Thus, we performed a comparative analysis of the dynamical properties of modified oligonucleotides based upon a sequence (5')r(UGUUACGACU)(3') that is complementary to bacterial ribosomal A-site RNA. In particular, 25 ns explicit solvent molecular dynamics simulations were computed for the following six single-stranded decamers: (1) the above RNA in unmodified form; (2) the 2'-O-methyl-modified RNA; (3) peptide nucleic acid (PNA) analogues of the above sequence, containing either (a) T or (b) U; and (4) two serine-substituted PNAs. Our results show that 2'-O-methylation attenuates RNA backbone dynamics, thereby preventing interconversion between stacked and unstacked conformations. The PNA analogue is rendered less flexible by replacing uracil with thymine; in addition, we found that derivatizing the PNA backbone with serine leads to enhanced base-stacking interactions. Consistent with known solubility properties of these classes of molecules, both RNAs exhibited greater localization of water molecules than did PNA. In terms of counterions, the initially helical conformation of the 2'-O-methyl RNA exhibits the highest Na(+) density among all the simulated decamers, while Na(+) build-up was most negligible for the neutral PNA systems. Further studies of the conformational and physicochemical properties of such modified single-stranded oligomers may facilitate better design of nucleic acid analogues, particularly those capable of serving as specific, high-affinity ribosomal A-site binders.

Molecular dynamics simulations of cyclohexyl modified peptide nucleic acids (PNA)

Peptide Nucleic Acids (PNA) that bind sequence specifically to DNA/RNA are of major interest in the field of molecular biology and could form the basis for gene-targeted drugs. Molecular dynamics simulations are aimed to characterize the structural and dynamical features to understand the effect of backbone modification on the structure and dynamics along with the stability of the resulting 10mer complexes of PNA with DNA/RNA. Twelve Molecular Dynamics (MD) simulations of duplexes and triplexes with and without cyclohexyl modification were carried out for 10ns each. The simulations indicate that the cyclohexyl modification with different stereoisomers has influenced all the PNA-DNA/RNA complexes. Modification has added rigidity to backbone by restricting beta to +60 in case of (1R,2S) cyclohexyl PNA and to -60 in case of (1S,2R) cyclohexyl PNA. The results of MD simulations were able to show the backbone rigidification and preference for RNA complexes over DNA due to presence of cyclohexyl ring in the PNA backbone.

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

Peptide nucleic acids (PNAs) efficiently hybridize with DNA and are promoted as versatile gene-targeting analytical tools and pharmaceuticals. However, PNAs have never been exploited as radiopharmaceuticals, and radiation-induced physicochemical modifications of PNA:DNA heteroduplexes have not been studied. Drug- and radiation-induced creation of covalent cross-links in DNA obstruct crucial cell survival processes such as transcription and replication and are thus considered genotoxic events with a high impact in anticancer therapies. Here we report that gamma-irradiation of complementary PNA:DNA heteroduplexes, wherein the PNA contains l-lysine, free amino, or N-methylmorpholinium N- and C-capping groups, results in the formation of irreversible interstrand cross-links (ICL). The number of detected ICL corresponds to the number of available amino functional groups on the PNA. The effect of DNA sequence on the formation of ICL was studied by modifying the terminal nucleotides of the DNA oligonucleotide to create deletions and overhangs. The involvement of abasic sites (ABS) on the DNA strand in the cross-linking reaction was confirmed by independent experiments with synthetic ABS-containing oligonucleotides. Molecular modeling and molecular dynamics (MD) simulations were applied to elucidate the conformation of the N- and C-capping groups of the PNA oligomer and their interactions with the proximal terminus of the DNA. Good agreement between experimental and modeling results was achieved. Modeling indicated that the presence of positively charged capping groups on the PNA increases the conformational flexibility of the PNA:DNA terminal base pairs and often leads to their melting. This disordered orientation of the duplex ends provides conditions for multiple encounters of the short (amino) and bulky (Lys) side chains with nucleobases and the DNA backbone up to the third base pair along the duplex stem. Dangling duplex ends offer favorable conditions for increased accessibility of the radiation-induced free radicals to terminal nucleotides and their damage. It is suggested that the ICL are produced by initial formation of Schiff base adducts between the PNA amino functions and the opposed DNA oxidation-damaged bases or abasic 2'-deoxyribose-derived aldehydic groups. The subsequent reduction by solvated electrons (e(-)(aq)) or other radiation-produced reducing species results in irreversible covalent interstrand cross-links. The simultaneous involvement of oxidizing, (*)OH, and reducing, e(-)(aq), radicals presents a case in which multiple ionization events along a gamma-particle path lead to DNA injuries that also encompass ICL as part of the multiply damaged sites (MDS). The obtained results may find applications in the development of a new generation of gene-targeted radiosensitizers based on PNA vectors.