Influence of PNA containing 8-aza-7-deazaadenine on structure stability and binding affinity of PNA·DNA duplex: insights from thermodynamics, counter ion, hydration and molecular dynamics analysis

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.

Conformational Dynamics of Cyanocobalamin and Its Conjugates with Peptide Nucleic Acids

Vitamin B₁₂ also called cobalamin (Cbl) is an important enzymatic cofactor taken up by mammalian and also by many bacterial cells. Peptide nucleic acid (PNA) is a synthetic DNA analogue that has the ability to bind in a complementary manner to natural nucleic acids. Provided that PNA is efficiently delivered to cells, it could act as a steric blocker of functional DNA or RNA and regulate gene expression at the level of transcription or translation. Recently, Cbl has been examined as a transporter of various molecules to cells. Also, PNA, if covalently linked with Cbl, can be delivered to bacterial cells, but it is crucial to verify that Cbl does not change the desired PNA biological properties. We have analyzed the structure and conformational dynamics of conjugates of Cbl with a PNA monomer and oligomer. We synthesized a cyanocobalamin derivative with a PNA monomer C connected via the triazole linker and determined its NMR spectra. Using microsecond-long molecular dynamics simulations, we examined the internal dynamics of cyanocobalamin-C, its conjugate with a 14-mer PNA, and free PNA. The results suggest that all compounds acquire rather compact structures but the PNA oligomer conformations vary. For the Cbl-C conjugate the cross-peaks from the ROESY spectrum corroborated with the clusters from molecular dynamics trajectories. Within PNA the dominant interaction is stacking but the stacking bases are not necessarily neighboring in the PNA sequence. More bases stack in free PNA than in PNA of the conjugate, but stacking is less stable in free PNA. PNA in the conjugate is slightly more exposed to solvent. Overall, cyanocobalamin attached to a PNA oligomer increases the flexibility of PNA in a way that could be beneficial for its hybridization with natural nucleic acid oligomers.

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.

Recognition of 2',5'-linked oligoadenylates by human ribonuclease L: molecular dynamics study

The capability of current MD simulations to be used as a tool in rational design of agonists of medically interesting enzyme RNase L was tested. Dimerization and enzymatic activity of RNase L is stimulated by 2',5'-linked oligoadenylates (pA₂₅A₂₅A; 2-5A). First, it was necessary to ensure that a complex of monomeric human RNase L and 25A was stable in MD simulations. It turned out that Glu131 had to be protonated. The non-protonated Glu131 caused dissociation of 2-5A from RNase L. Because of the atypical 2'-5' internucleotide linkages and a specific spatial arrangement of the 25A trimer, when a single molecule carries all possible conformers of the glycosidic torsion angle, several versions of the AMBER force field were tested. One that best maintained functionally important interactions of 25A and RNase L was selected for subsequent MD simulations. Furthermore, we wonder whether powerful GPUs are able to produce MD trajectories long enough to convincingly demonstrate effects of subtle perturbations of interactions between 25A and RNase L. Detrimental impacts of various point mutations of RNase L (R155A, F126A, W60A, K89A) on 2-5A binding were observed on a time scale of 200 ns. Finally, 2-5A analogues with a bridged 3'--O,4'--C-alkylene linkage (B) introduced into the adenosine units (A) were used to assess ability of MD simulations to distinguish on the time scale of hundreds of nanoseconds between agonists of RNase L (pA₂₅A₂₅B, pB₂₅A₂₅A, pB₂₅A₂₅B) and inactive analogs (pA₂₅B₂₅A, pA₂₅B₂₅B, pB₂₅B₂₅A, pB₂₅B₂₅B). Agonists were potently bound to RNase L during 200 ns MD runs. For inactive 2-5A analogs, by contrast, significant disruptions of their interactions with RNase L already within 100 ns MD runs were found.