Molecular Dynamics Simulation of a PNA·DNA·PNA Triple Helix in Aqueous Solution

Thermodynamic Parameter Estimation for Modified Oligonucleotides Using Molecular Dynamics Simulations

This study investigates the thermodynamic parameters of 1300 RNA/DNA hybrid duplexes, including both natural and chemically modified forms, using molecular dynamics (MD) simulations. Modified duplexes consist of phosphorothioate (PS) backbones and 2'-O-methoxyethyl (MOE) modifications, both commonly used in therapeutic oligonucleotides. Hybridization enthalpy and entropy were calculated from MD trajectories using molecular mechanics Poisson-Boltzmann surface area (MMPBSA) and molecular mechanics generalized Born surface area (MMGBSA) approaches. To address discrepancies with experimental data, we established empirical relationships by comparing calculated values with known experimental results of natural hybrid duplexes, then extended these relationships to the entire data set. The corrected parameters were subsequently used to generate nearest-neighbor (NN) models, allowing for experimentally reliable melting temperature predictions. In this process, MMGBSA demonstrated superior predictive performance with high convergence and consistency for both natural and modified duplexes. Specifically, MMGBSA captured the stabilizing effects of the MOE modification with minimal bias, while MMPBSA exhibited greater variability and limited reliability. These findings highlight the potential of MMGBSA for accurate thermodynamic modeling of both natural and modified nucleic acids, providing a robust framework and experimentally meaningful insights for applications in nucleic acid-based therapeutic design and biotechnology.

Impact of Intracellular Proteins on μ-Opioid Receptor Structure and Ligand Binding

Chronic pain is a prevalent problem affecting approximately one out of every five adults in the U.S. The most effective way to treat chronic pain is with opioids, but they cause dangerous side effects such as tolerance, addiction, and respiratory depression, which makes them quite deadly. Opioids, such as fentanyl, target the μ-opioid receptor (MOR), which can then bind to the intracellular Gi protein or the β-arrestin protein. The Gi pathway is primarily responsible for pain relief and potential side effects, but the β-arrestin pathway is chiefly responsible for the unwanted side effects. Ideally, an effective pain medication without side effects would bind to MOR, which would bias signaling solely through the Gi pathway. We used the Bio3D library to conduct principal component analysis to compare the cryo-electron microscopy MOR structure in complex with the Gi versus an X-ray crystallography MOR structure with a nanobody acting as a Gi mimic. Our results agree with a previous study by Munro, which concluded that nanobody-bound MOR is structurally different than Gi-bound MOR. Furthermore, we investigated the structural diversity of opioids that can bind to MOR. Quantum mechanical calculations show that the low energy solution structures of fentanyl differ from the one bound to MOR in the experimental structure, and pKa calculations reveal that fentanyl is protonated in aqueous solution. Glide docking studies show that higher energy structures of fentanyl in solution form favorable docking complexes with MOR. Our calculations show the relative abundance of each fentanyl conformation in solution as well as the energetic barriers that need to be overcome to bind to MOR. Docking studies confirm that multiple fentanyl conformations can bind to the receptor. Perhaps a variety of conformations of fentanyl can stabilize multiple conformations of the MOR, which can explain why fentanyl can induce different intracellular signaling and multiple physiological effects.

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.

Three's a crowd - stabilisation, structure, and applications of DNA triplexes

DNA is a strikingly flexible molecule and can form a variety of secondary structures, including the triple helix, which is the subject of this review. The DNA triplex may be formed naturally, during homologous recombination, or can be formed by the introduction of a synthetic triplex forming oligonucleotide (TFO) to a DNA duplex. As the TFO will bind to the duplex with sequence specificity, there is significant interest in developing TFOs with potential therapeutic applications, including using TFOs as a delivery mechanism for compounds able to modify or damage DNA. However, to combine triplexes with functionalised compounds, a full understanding of triplex structure and chemical modification strategies, which may increase triplex stability or in vivo degradation, is essential - these areas will be discussed in this review. Ruthenium polypyridyl complexes, which are able to photooxidise DNA and act as luminescent DNA probes, may serve as a suitable photophysical payload for a TFO system and the developments in this area in the context of DNA triplexes will also be reviewed.

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.

Molecularly resolved, label-free nucleic acid sensing at solid-liquid interface using non-ionic DNA analogues

Nucleic acid-based biosensors, where the capture probe is a nucleic acid, e.g., DNA or its synthetic analogue xeno nucleic acid (XNA), offer interesting ways of eliciting clinically relevant information from hybridization/dehybridization signals. In this respect, the application of XNA probes is attractive since the drawbacks of DNA probes might be overcome. Within the XNA probe repertoire, peptide nucleic acid (PNA) and morpholino (MO) are promising since their backbones are non-ionic. Therefore, in the absence of electrostatic charge repulsion between the capture probe and the target nucleic acid, a stable duplex can be formed. In addition, these are nuclease-resistant probes. Herein, we have tested the molecularly resolved nucleic acid sensing capacity of PNA and MO capture probes using a fluorescent label-free single molecule force spectroscopy approach. As far as single nucleobase mismatch discrimination is concerned, both PNA and MO performed better than DNA, while the performance of the MO probe was the best. We propose that the conformationally more rigid backbone of MO, compared to the conformationally flexible PNA, is an advantage for MO, since the probe orientation can be made more upright on the surface and therefore MO can be more effectively accessed by the target sequences. The performance of the XNA probes has been compared to that of the DNA probe, using fixed nucleobase sequences, so that the effect of backbone variation could be investigated. To our knowledge, this is the first report on molecularly resolved nucleic acid sensing by non-ionic capture probes, here, MO and PNA.

Improved nucleic acid sensing in terms of single nucleobase mismatch discrimination, as achieved by the surface-confined non-ionic PNA and MO capture probes, is exemplified by single molecule force spectroscopy.

Nanoscale structures and mechanics of peptide nucleic acids

Peptide nucleic acids are charge-neutral polyamide oligomers with extremely flexible backbones that have a strong affinity for hybridization with complementary DNA or RNA, as well as encouraging antisense and antigene activity in cell-free systems.

Antisense peptide nucleic acid inhibits the growth of KPC-producing Klebsiella pneumoniae strain

Klebsiella pneumoniae causes common and severe hospital- and community-acquired infections with a high incidence of multidrug resistance (MDR) and mortality. In this study, we investigated the ability of the antisense peptide nucleic acids (PNA) conjugated to the (KFF)3K cell-penetrating peptide (CPP) to target the gyrA KPC-producing K. pneumoniae and inhibit bacterial growth in vitro. The inhibitory effect on gyrA gene was evaluated by measuring 16s gene amplification in KPC-producing K. pneumoniae treated with the antisense PNA conjugate. The hemolytic property of the antisense PNA conjugate was accessed toward mice red blood cells. Finally, molecular modeling and dynamics simulations analyses in aqueous solutions were performed to predict the PNA conformation alone in contact with DNA (gyrA gene sequence). PNA was capable of inhibiting bacterial growth at 50 μM, also reducing 16S gene amplification in 96.7%. Besides, PNA presented low hemolytic activity (21.1% hemolysis) at this same concentration. Bioinformatics analysis demonstrated that the structure of the PNA is stable in water without major changes in its secondary structure. The ability of PNA and its conjugated CPP ((KFF)3K) to inhibit bacterial growth demonstrates the potential of this new class of antibacterial agents, encouraging further in vivo studies to confirm its therapeutic efficacy.

Structure-Based Approaches to Classify the Functional Impact of ZBTB18 Missense Variants in Health and Disease

The Cys₂His₂ type zinc finger is a motif found in many eukaryotic transcription factor proteins that facilitates binding to genomic DNA so as to influence cellular gene expression. One such transcription factor is ZBTB18, characterized as a repressor that orchestrates the development of mammalian tissues including skeletal muscle and brain during embryogenesis. In humans, it has been recognized that disease-associated ZBTB18 missense variants mapping to the coding sequence of the zinc finger domain influence sequence-specific DNA binding, disrupt transcriptional regulation, and impair neural circuit formation in the brain. Furthermore, general population ZBTB18 missense variants that influence DNA binding and transcriptional regulation have also been documented within this domain; however, the molecular traits that explain why some variants cause disease while others do not are poorly understood. Here, we have applied five structure-based approaches to evaluate their ability to discriminate between disease-associated and general population ZBTB18 missense variants. We found that thermodynamic integration and Residue Scanning in the Schrodinger Biologics Suite were the best approaches for distinguishing disease-associated variants from general population variants. Our results demonstrate the effectiveness of structure-based approaches for the functional characterization of missense alleles to DNA binding, zinc finger transcription factor protein-coding genes that underlie human health and disease.

The proto-Nucleic Acid Builder: a software tool for constructing nucleic acid analogs

The helical structures of DNA and RNA were originally revealed by experimental data. Likewise, the development of programs for modeling these natural polymers was guided by known structures. These nucleic acid polymers represent only two members of a potentially vast class of polymers with similar structural features, but that differ from DNA and RNA in the backbone or nucleobases. Xeno nucleic acids (XNAs) incorporate alternative backbones that affect the conformational, chemical, and thermodynamic properties of XNAs. Given the vast chemical space of possible XNAs, computational modeling of alternative nucleic acids can accelerate the search for plausible nucleic acid analogs and guide their rational design. Additionally, a tool for the modeling of nucleic acids could help reveal what nucleic acid polymers may have existed before RNA in the early evolution of life. To aid the development of novel XNA polymers and the search for possible pre-RNA candidates, this article presents the proto-Nucleic Acid Builder (https://github.com/GT-NucleicAcids/pnab), an open-source program for modeling nucleic acid analogs with alternative backbones and nucleobases. The torsion-driven conformation search procedure implemented here predicts structures with good accuracy compared to experimental structures, and correctly demonstrates the correlation between the helical structure and the backbone conformation in DNA and RNA.

Molecular insights on the dynamic stability of peptide nucleic acid functionalized carbon and boron nitride nanotubes

The strategic approaches to the design of self-assembled hybrids of biomolecular systems at the nanoscale such as deoxyribonucleic acid (DNA) with single-wall carbon nanotubes (CNTs) and their structural analog, boron nitride nanotubes (BNNTs), rely on understanding how biomolecules recognize and mediate intermolecular interactions with the nanomaterial's surface. In this paper, we consider peptide nucleic acid (PNA), which is a synthetic analog of DNA, and investigate its interaction with a zigzag CNT and BNNT of similar diameter. The results based on the molecular dynamics method find that PNA provides definitive contrasts in the adsorption on the tubular surface in aqueous solution: it prefers to wrap along the circumferential direction on a (11,0) CNT, whereas it binds along the axial direction adopting an extended configuration on a (11,0) BNNT. Moreover, gas-phase Monte Carlo simulations show a dependence of the nanotube diameter on the calculated adsorption energy, with BNNTs exhibiting higher adsorption energy compared to CNTs, and the largest-diameter (25,0) tubular configuration facilitates encapsulation of PNA rather than PNA being adsorbed on its sidewall. The results are expected to be of relevance in the design of novel PNA-based archetypal hybrid materials for nanoscale applications in health-related areas including biosensing.

Vitamin B12-peptide nucleic acids use the BtuB receptor to pass through the Escherichia coli outer membrane

Short modified oligonucleotides that bind in a sequence-specific way to messenger RNA essential for bacterial growth could be useful to fight bacterial infections. One such promising oligonucleotide is peptide nucleic acid (PNA), a synthetic DNA analog with a peptide-like backbone. However, the limitation precluding the use of oligonucleotides, including PNA, is that bacteria do not import them from the environment. We have shown that vitamin B₁₂, which most bacteria need to take up for growth, delivers PNAs to Escherichia coli cells when covalently linked with PNAs. Vitamin B₁₂ enters E. coli via a TonB-dependent transport system and is recognized by the outer-membrane vitamin B₁₂-specific BtuB receptor. We engineered the E. coli ΔbtuB mutant and found that transport of the vitamin B₁₂-PNA conjugate requires BtuB. Thus, the conjugate follows the same route through the outer membrane as taken by free vitamin B₁₂. From enhanced sampling all-atom molecular dynamics simulations, we determined the mechanism of conjugate permeation through BtuB. BtuB is a β-barrel occluded by its luminal domain. The potential of mean force shows that conjugate passage is unidirectional and its movement into the BtuB β-barrel is energetically favorable upon luminal domain unfolding. Inside BtuB, PNA extends making its permeation mechanically feasible. BtuB extracellular loops are actively involved in transport through an induced-fit mechanism. We prove that the vitamin B₁₂ transport system can be hijacked to enable PNA delivery to E. coli cells.

Dynamics Studies of DNA with Non-canonical Structure Using NMR Spectroscopy

The non-canonical structures of nucleic acids are essential for their diverse functions during various biological processes. These non-canonical structures can undergo conformational exchange among multiple structural states. Data on their dynamics can illustrate conformational transitions that play important roles in folding, stability, and biological function. Here, we discuss several examples of the non-canonical structures of DNA focusing on their dynamic characterization by NMR spectroscopy: (1) G-quadruplex structures and their complexes with target proteins; (2) i-motif structures and their complexes with proteins; (3) triplex structures; (4) left-handed Z-DNAs and their complexes with various Z-DNA binding proteins. This review provides insight into how the dynamic features of non-canonical DNA structures contribute to essential biological processes.

Disease-associated missense variants in ZBTB18 disrupt DNA binding and impair the development of neurons within the embryonic cerebral cortex

The activities of DNA-binding transcription factors, such as the multi-zinc-finger protein ZBTB18 (also known as RP58, or ZNF238), are essential to coordinate mammalian neurodevelopment, including the birth and radial migration of newborn neurons within the fetal brain. In humans, the majority of disease-associated missense mutations in ZBTB18 lie within the DNA-binding zinc-finger domain and are associated with brain developmental disorder, yet the molecular mechanisms explaining their role in disease remain unclear. To address this, we developed in silico models of ZBTB18, bound to DNA, and discovered that half of the missense variants map to residues (Asn461, Arg464, Glu486) predicted to be essential to sequence-specific DNA contact, whereas others map to residues (Leu434, Tyr447, Arg495) with limited contributions to DNA binding. We studied pathogenic variants to residues with close (N461S) and limited (R495G) DNA contact and found that each bound DNA promiscuously, displayed altered transcriptional regulatory activity in vitro, and influenced the radial migration of newborn neurons in vivo in different ways. Taken together, our results suggest that altered transcriptional regulation could represent an important pathological mechanism for ZBTB18 missense variants in brain developmental disease.

Synthetic, Structural, and RNA Binding Studies on 2-Aminopyridine-Modified Triplex-Forming Peptide Nucleic Acids

The development of new RNA-binding ligands is attracting increasing interest in fundamental science and the pharmaceutical industry. The goal of this study was to improve the RNA binding properties of triplex-forming peptide nucleic acids (PNAs) by further increasing the pK_a of 2-aminopyridine (M). Protonation of M was the key for enabling triplex formation at physiological pH in earlier studies. Substitution on M by an electron-donating 4-methoxy substituent resulted in slight destabilization of the PNA-dsRNA triplex, contrary to the expected stabilization due to more favorable protonation. To explain this unexpected result, the first NMR structural studies were performed on an M-modified PNA-dsRNA triplex which, combined with computational modeling identified unfavorable steric and electrostatic repulsion between the 4-methoxy group of M and the oxygen of the carbonyl group connecting the adjacent nucleobase to PNA backbone. The structural studies also provided insights into hydrogen-bonding interactions that might be responsible for the high affinity and unusual RNA-binding preference of PNAs.

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.

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.

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.

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

This paper describes the synthesis of a novel 8-aza-7-deazapurin-2,6-diamine (DPP)-containing peptide nucleic acid (PNA) monomer and Boc protecting group-based oligomerization of PNA, replacing adenine (A) with DPP monomers in the PNA strand. The PNA oligomers were synthesized against the biologically relevant SV40 promoter region (2494-AATTTTTTTTATTTA-2508) of pEGFP-N3 plasmid. The DPP-PNA·DNA duplex showed enhanced stability as compared to normal duplex (A-PNA·DNA). The electronic distribution of DPP monomer suggested that DPP had better electron donor properties over 2,6-diamino purine. UV melting and thermodynamic analysis revealed that the PNA oligomer containing a diaminopyrazolo(3,4-d)pyrimidine moiety (DPP) stabilized the PNA·DNA hybrids compared to A-PNA·DNA. DPP-PNA·DNA duplex showed higher water activity (Δnw = 38.5) in comparison to A-PNA·DNA duplex (Δnw = 14.5). The 50 ns molecular dynamics simulations of PNA·DNA duplex containing DPP or unmodified nucleobase-A showed average H-bond distances in the DPP-dT base pair of 2.90 Å (OH-N bond) and 2.91 Å (NH-N bond), which were comparably shorter than in the A-dT base pair, in which the average distances were 3.18 Å (OH-N bond) and 2.97 Å (NH-N bond), and there was one additional H-bond in the DPP-dT base pair of around 2.98 Å (O2H-N2 bond), supporting the higher stability of DPP-PNA·DNA. The analysis of molecular dynamics simulation data showed that the system binding free energy increased at a rate of approximately -4.5 kcal mol(-1) per DPP base of the PNA·DNA duplex. In summary, increased thermal stability, stronger hydrogen bonding and more stable conformation in the DPP-PNA·DNA duplex make it a better candidate as antisense/antigene therapeutic agents.

The DNA and RNA sugar-phosphate backbone emerges as the key player. An overview of quantum-chemical, structural biology and simulation studies

Knowledge of geometrical and physico-chemical properties of the sugar-phosphate backbone substantially contributes to the comprehension of the structural dynamics, function and evolution of nucleic acids. We provide a side by side overview of structural biology/bioinformatics, quantum chemical and molecular mechanical/simulation studies of the nucleic acids backbone. We highlight main features, advantages and limitations of these techniques, with a special emphasis given to their synergy. The present status of the research is then illustrated by selected examples which include classification of DNA and RNA backbone families, benchmark structure-energy quantum chemical calculations, parameterization of the dihedral space of simulation force fields, incorporation of arsenate into DNA, sugar-phosphate backbone self-cleavage in small RNA enzymes, and intricate geometries of the backbone in recurrent RNA building blocks. Although not apparent from the current literature showing limited overlaps between the QM, simulation and bioinformatics studies of the nucleic acids backbone, there in fact should be a major cooperative interaction between these three approaches in studies of the sugar-phosphate backbone.

Long-timescale molecular-dynamics simulations of the major urinary protein provide atomistic interpretations of the unusual thermodynamics of ligand binding

The mouse major urinary protein (MUP) has proved to be an intriguing test bed for detailed studies on protein-ligand recognition. NMR, calorimetric, and modeling investigations have revealed that the thermodynamics of ligand binding involve a complex interplay between competing enthalpic and entropic terms. We performed six independent, 1.2 micros molecular-dynamics simulations on MUP--three replicates on the apo-protein, and three on the complex with the pheromone isobutylmethoxypyrazine. Our findings provide the most comprehensive picture to date of the structure and dynamics of MUP, and how they are modulated by ligand binding. The mechanical pathways by which amino acid side chains can transmit information regarding ligand binding to surface loops and either increase or decrease their flexibility (entropy-entropy compensation) are identified. Dewetting of the highly hydrophobic binding cavity is confirmed, and the results reveal an aspect of ligand binding that was not observed in earlier, shorter simulations: bound ligand retains extensive rotational freedom. Both of these features have significant implications for interpretations of the entropic component of binding. More generally, these simulations test the ability of current molecular simulation methods to produce a reliable and reproducible picture of protein dynamics on the microsecond timescale.

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.

Structural properties of g,t-parallel duplexes

The structure of G,T-parallel-stranded duplexes of DNA carrying similar amounts of adenine and guanine residues is studied by means of molecular dynamics (MD) simulations and UV- and CD spectroscopies. In addition the impact of the substitution of adenine by 8-aminoadenine and guanine by 8-aminoguanine is analyzed. The presence of 8-aminoadenine and 8-aminoguanine stabilizes the parallel duplex structure. Binding of these oligonucleotides to their target polypyrimidine sequences to form the corresponding G,T-parallel triplex was not observed. Instead, when unmodified parallel-stranded duplexes were mixed with their polypyrimidine target, an interstrand Watson-Crick duplex was formed. As predicted by theoretical calculations parallel-stranded duplexes carrying 8-aminopurines did not bind to their target. The preference for the parallel-duplex over the Watson-Crick antiparallel duplex is attributed to the strong stabilization of the parallel duplex produced by the 8-aminopurines. Theoretical studies show that the isomorphism of the triads is crucial for the stability of the parallel triplex.

Molecular renormalization group coarse-graining of polymer chains: application to double-stranded DNA

Coarse-graining of atomistic force fields allows us to investigate complex biological problems, occurring at long timescales and large length scales. In this work, we have developed an accurate coarse-grained model for double-stranded DNA chain, derived systematically from atomistic simulations. Our approach is based on matching correlators obtained from atomistic and coarse-grained simulations, for observables that explicitly enter the coarse-grained Hamiltonian. We show that this requirement leads to equivalency of the corresponding partition functions, resulting in a one-step renormalization. Compared to prior works exploiting similar ideas, the main novelty of this work is the introduction of a highly compact set of Hamiltonian basis functions, based on molecular interaction potentials. We demonstrate that such compactification allows us to reproduce many-body effects, generated by one-step renormalization, at low computational cost. In addition, compact Hamiltonians greatly increase the likelihood of finding unique solutions for the coarse-grained force-field parameter values. By successfully applying our molecular renormalization group coarse-graining technique to double-stranded DNA, we solved, for the first time, a long-standing problem in coarse-graining polymer systems, namely, how to accurately capture the correlations among various polymeric degrees of freedom. Excellent agreement is found among atomistic and coarse-grained distribution functions for various structural observables, including those not included in the Hamiltonian. We also suggest higher-order generalization of this method, which may allow capturing more subtle correlations in biopolymer dynamics.

Computational approaches for the design of peptides with anti-breast cancer properties

Background: Breast cancer is the most common cancer among women. Tamoxifen is the preferred drug for estrogen receptor-positive breast cancer treatment, yet many of these cancers are intrinsically resistant to tamoxifen or acquire resistance during treatment. Therefore, scientists are searching for breast cancer drugs that have different molecular targets. Methodology: Recently, a computational approach was used to successfully design peptides that are new lead compounds against breast cancer. We used replica exchange molecular dynamics to predict the structure and dynamics of active peptides, leading to the discovery of smaller bioactive peptides. Conclusions: These analogs inhibit estrogen-dependent cell growth in a mouse uterine growth assay, a test showing reliable correlation with human breast cancer inhibition. We outline the computational methods that were tried and used along with the experimental information that led to the successful completion of this research.

Destabilization of DNA duplexes by oxidative damage at guanine: implications for lesion recognition and repair

We have used molecular dynamics simulations to study the structure and dynamics of a range of DNA duplexes containing the 2,6-diamino-4-hydroxy-5-formamidopyrimidine (FapydG) lesion that can result from oxidative damage at guanine. Compared to the corresponding undamaged DNA duplexes, FapydG-containing duplexes show little gross structural changes-the damaged base remains stacked in to the DNA double helix and retains hydrogen bonds to its cytosine partner. However, the experimentally observed reduction in DNA stability that accompanies lesion formation can be explained by a careful energetic analysis of the simulation data. Irrespective of the nature of the base pairs on either side of the lesion site, conversion of a guanine to a FapydG base results in increased dynamical flexibility in the base (but not in the DNA as a whole) that significantly weakens its hydrogen-bonding interactions. Surprisingly, the stacking interactions with its neighbours are not greatly altered. The formamido group adopts a non-planar conformation that can interact significantly and in a sequence-dependent manner with its 3'-neighbour. We conclude that the recognition of FapydG lesions by the repair protein formamidopyrimidine-DNA glycosylase probably does not involve the protein capturing an already-extrahelical FapydG base, but rather it relies on detecting alterations to the DNA structure and flexibility created by the lesion site.

Insight into why pyrrolidinyl peptide nucleic acid binding to DNA is more stable than the DNA x DNA duplex

Molecular dynamics (MD) simulations and experimental measurements of the stability of a novel pyrrolidinyl PNA binding to DNA (PNA x DNA) in both parallel and antiparallel configurations were carried out. For comparison, simulations were also performed for the DNA x DNA duplex. The conformations of the three simulated systems were found to retain well-defined base pairing and base stacking as their starting B-like structure. A large gas-phase energy repulsion of the two negatively charged sugar-phosphate backbones of the DNA strands was found to reduce the stability of the DNA x DNA duplex significantly compared with that of the PNA x DNA complexes, especially in the antiparallel binding configuration. In addition, the antiparallel PNA x DNA was observed to be less solvated than that of the other two systems. The simulated binding free energies and the experimental melting temperatures for the three investigated systems are in good agreement, indicating that the antiparallel PNA x DNA is the most stable duplex.

Towards a molecular dynamics consensus view of B-DNA flexibility

We present a systematic study of B-DNA flexibility in aqueous solution using long-scale molecular dynamics simulations with the two more recent versions of nucleic acids force fields (CHARMM27 and parmbsc0) using four long duplexes designed to contain several copies of each individual base pair step. Our study highlights some differences between pambsc0 and CHARMM27 families of simulations, but also extensive agreement in the representation of DNA flexibility. We also performed additional simulations with the older AMBER force fields parm94 and parm99, corrected for non-canonical backbone flips. Taken together, the results allow us to draw for the first time a consensus molecular dynamics picture of B-DNA flexibility.

Refinement of the AMBER force field for nucleic acids: improving the description of alpha/gamma conformers

We present here the parmbsc0 force field, a refinement of the AMBER parm99 force field, where emphasis has been made on the correct representation of the alpha/gamma concerted rotation in nucleic acids (NAs). The modified force field corrects overpopulations of the alpha/gamma = (g+,t) backbone that were seen in long (more than 10 ns) simulations with previous AMBER parameter sets (parm94-99). The force field has been derived by fitting to high-level quantum mechanical data and verified by comparison with very high-level quantum mechanical calculations and by a very extensive comparison between simulations and experimental data. The set of validation simulations includes two of the longest trajectories published to date for the DNA duplex (200 ns each) and the largest variety of NA structures studied to date (15 different NA families and 97 individual structures). The total simulation time used to validate the force field includes near 1 mus of state-of-the-art molecular dynamics simulations in aqueous solution.

Theoretical study of large conformational transitions in DNA: the B<-->A conformational change in water and ethanol/water

We explore here the possibility of determining theoretically the free energy change associated with large conformational transitions in DNA, like the solvent-induced B<-->A conformational change. We find that a combination of targeted molecular dynamics (tMD) and the weighted histogram analysis method (WHAM) can be used to trace this transition in both water and ethanol/water mixture. The pathway of the transition in the A-->B direction mirrors the B-->A pathway, and is dominated by two processes that occur somewhat independently: local changes in sugar puckering and global rearrangements (particularly twist and roll) in the structure. The B-->A transition is found to be a quasi-harmonic process, which follows closely the first spontaneous deformation mode of B-DNA, showing that a physiologically-relevant deformation is in coded in the flexibility pattern of DNA.

Molecular dynamics study of small PNA molecules in lipid-water system

We present the results of molecular dynamics simulations of small peptide nucleic acid (PNA) molecules, synthetic analogs of DNA, at a lipid bilayer in water. At neutral pH, without any salt, and in the NP(n)gammaT ensemble, two similar PNA molecules (6-mers) with the same nucleic base sequence and different terminal groups are investigated at the interface between water and a 1-palmitoyl-2-oleoylphosphatidylcholine lipid bilayer. The results of our simulations suggest that at low ionic strength of the solution, both PNA molecules adsorb at the lipid-water interface. In the case where the PNA molecule has charged terminal groups, the main driving force of adsorption is the electrostatic attraction between the charged groups of PNA and the lipid heads. The main driving force of adsorption of the PNA molecule with neutral terminal groups is the hydrophobic interaction of the nonpolar groups. Our simulations suggest that the system free energy change associated with PNA adsorption at the lipid-water interface is on the order of several tens of kT per PNA molecule in both cases.

Theoretical Examination of Two Opposite Mechanisms Proposed for Hepatitis Delta Virus Ribozyme

Two mechanisms were previously proposed for the hepatitis delta virus (HDV) ribozyme where an active-site cytosine residue (C75) either functioned as a general base to deprotonate the 2'-OH at the rupture site or as a general acid to protonate the O5' leaving group. Here, we reported the first theoretical examination of the two mechanisms using a combination of the quantum mechanics (QM)/molecular mechanics (MM), molecular dynamics (MD), and near-attack-conformation (NAC) techniques. Our theoretical results supported the C75-acid mechanism, which was demonstrated to have an unfavorable starting geometry (in agreement with the crystallographic data) but a significantly lower energy barrier as compared to the C75-base mechanism. Therefore, the chemical details of the transition state in the HDV ribozyme may dramatically differ from those inferred from the structural studies.

Orientation preferences of backbone secondary amide functional groups in peptide nucleic acid complexes: quantum chemical calculations reveal an intrinsic preference of cationic D-amino acid-based chiral PNA analogues for the P-form

Geometric descriptions of nonideal interresidue hydrogen bonding and backbone-base water bridging in the minor groove are established in terms of polyamide backbone carbonyl group orientation from analyses of residue junction conformers in experimentally determined peptide nucleic acid (PNA) complexes. Two types of interresidue hydrogen bonding are identified in PNA conformers in heteroduplexes with nucleic acids that adopt A-like basepair stacking. Quantum chemical calculations on the binding of a water molecule to an O2 base atom in glycine-based PNA thymine dimers indicate that junctions modeled with P-form backbone conformations are lower in energy than a dimer comprising the predominant conformation observed in A-like helices. It is further shown in model systems that PNA analogs based on D-lysine are better able to preorganize in a conformation exclusive to P-form helices than is glycine-based PNA. An intrinsic preference for this conformation is also exhibited by positively charged chiral PNA dimers carrying 3-amino-D-alanine or 4-aza-D-leucine residue units that provide for additional rigidity by side-chain hydrogen bonding to the backbone carbonyl oxygen. Structural modifications stabilizing P-form helices may obviate the need for large heterocycles to target DNA pyrimidine bases via PNA.DNA-PNA triplex formation. Quantum chemical modeling methods are used to propose candidate PNA Hoogsteen strand designs.

Electrostatic, Steric, and Hydration Interactions Favor Na+ Condensation around DNA Compared with K+

Condensation of monovalent counterions around DNA influences polymer properties of the DNA chain. For example, the Na(+) ions show markedly stronger propensity to induce multiple DNA chains to assemble into compact structures compared with the K(+) ions. To investigate the similarities and differences in the sodium and potassium ion condensation around DNA, we carried out a number of extensive all-atom molecular dynamics simulations of a DNA oligomer consisting of 16 base pairs, [d(CGAGGTTTAAACCTCG)](2), in explicit water. We found that the Na(+) ions penetrate the DNA interior and condense around the DNA exterior to a significantly larger degree compared with the K(+) ions. We have provided a microscopic explanation for the larger Na(+) affinity toward DNA that is based on a combination of steric, electrostatic, and hydration effects. Unexpectedly, we found that the Cl(-) co-ions provide more efficient electrostatic screening for the K(+) ions than for the Na(+) ions, contributing to the larger Na(+) condensation around DNA. To examine the importance of the discrete nature of water and ions, we also computed the counterion distributions from the mean-field electrostatic theory, demonstrating significant disagreements with the all-atom simulations. Prior experimental results on the relative extent of the Na(+) and K(+) condensation around DNA were somewhat contradictory. Recent DNA compaction experiments may be interpreted to suggest stronger Na(+) condensation around DNA compared to K(+), which is consistent with our simulations. We also provide a simple interpretation for the experimentally observed increase in DNA electrophoretic mobility in the alkali metal series, Li(+) < Na(+) < K(+) < Rb(+). We compare the DNA segment conformational preferences in various buffers with the proposed NMR models.

G-Quadruplexes Can Maintain Their Structure in the Gas Phase

Several very extended (0.5-1 micros) molecular dynamics (MD) simulations of parallel and antiparallel G-quadruplex DNA strongly suggest that in the presence of suitable cations the quadruplex not only remains stable in the gas phase, but also displays a structure that closely resembles that found in extended (25-ns long) trajectories in aqueous solution. In the absence of the crucial cations, the trajectories become unstable and in general the quadruplex structure is lost. To our knowledge, this is the first physiologically relevant structure of DNA for which very large MD simulations suggest that the structure in water and in the gas phase are indistinguishable.

Thermodynamics and kinetics of PNA-DNA quadruplex-forming chimeras

PNA-DNA chimeras present the interesting properties of PNA, such as the high binding affinity to complementary single-strand (DNA or RNA), and the resistance to nuclease and protease degradation. At the same time, the limitations of an oligomer containing all PNA residues, such as low water solubility, self-aggregation, and low cellular uptake, are effectively overcome. Further, PNA-DNA chimeras possess interesting biological properties as antisense agents. We have explored the ability of PNA-DNA chimeric strands to assemble in quadruplex structures. The rate constant for association of the quadruplexes and their thermodynamic properties have been determined by CD spectroscopy and differential scanning calorimetry (DSC). Thermal denaturation experiments indicated higher thermal and thermodynamic stabilities for chimeric quadruplexes in comparison with the corresponding unmodified DNA quadruplex. Singular value decomposition analysis (SVD) suggests the presence of kinetically stable intermediate species in the quadruplex formation process. The experimental results have been discussed on the basis of molecular dynamic simulations. The ability of PNA-DNA chimeras to form stable quadruplex structures expands their potential utility as therapeutic agents.