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

Comparative molecular dynamics calculations of duplexation of chemically modified analogs of DNA used for antisense applications

We have subjected several analogs of DNA that have been widely used as antisense oligonucleotide (ASO) inhibitors of gene expression to comparative molecular dynamics (MD) calculations of their ability to form duplexes with DNA and RNA. The analogs included in this study are the phosphorothioate (PS), peptide nucleic acid (PNA), locked nucleic acid (LNA), morpholino nucleic acid (PMO), the 2'-OMe, 2'-F, 2'-methoxyethyl (2'-MOE) and the constrained cET analogs, as well as the natural phosphodiester (PO) as control, for a total of nine structures, in both XNA-DNA and XNA-RNA duplexes. This is intended as an objective criterion for their relative ability to duplex with an RNA complement and their comparative potential for antisense applications. We have found that the constrained furanose ring analogs show increased stability when considering this study's structural and energetic parameters. The 2'-MOE modification, even though energetically stable, has an elevated dynamic range and breathing properties due to the bulkier moiety in the C2' position of the furanose. The smaller modifications in the C2' position, 2'-F, 2'-OMe and PS also form stable and energetically favored duplexes with both DNA and RNA. The morpholino moiety allows for increased tolerance in accommodating either DNA or RNA and the PNA, with the PNA being the most energetically stable, although with a preference for the B-form DNA. In summary, we can rank the overall preference of hybrid strand formations as PNA > cET/LNA > PS/2'-F/2'-OMe > morpholino > 2'-MOE for the efficacy of duplex formation.

Enhanced molecular dynamic simulation studies unravel long-range effects caused by sequence variations and partner binding in RNA aptamers

Intrinsic flexibility and structural modularity are two common features of RNA molecules. Although functionally crucial, RNA plasticity often represents a major complication in high-resolution structural studies. To overcome this problem, RNAs may be rigidified through the complexation with high-affinity partners such as Fab molecules. This approach has been previously used to characterize the DIR2-aptamer. However, possible perturbations induced by the insertion of the Fab binding site on the DIR2-aptamer conformational properties were not investigated. Here, using enhanced molecular dynamics simulations, we compared the dynamics of the DIR2 aptamer holding the Fab binding site with that of the parental sequence. Our results suggest that the L2-loop modification for the Fab recognition leads to a significant increase in local flexibility that also affects the mobility of distant regions. The trajectories provide clear indications of the groups and the interactions mediating the dynamics transfer in DIR2. The effectiveness of our approach in addressing RNA flexibility was further corroborated by showing its ability to reproduce the most important events affecting the NF-κB RNA aptamer upon dissociation from the partner. Therefore, REMD analyses, a rarely adopted technique to unravel the structural/dynamical properties of aptamers, could efficiently complement experimental data guiding the rational design of nucleic acid therapeutics.

Bioactive compounds from Pandanous fascicularis as potential therapeutic candidate to tackle hepatitis a inhibition: Docking and molecular dynamics simulation study

Due to extensive pharmacological research, medicinal plants the underpinning of indigenous herbal serve as a possible source of key compounds for the development of new drugs. Hepatitis A, one of the most widespread infectious diseases associated with global public health issues. The transmission of hepatitis A virus (HAV) occurs, through personal contact, as well as contaminated food/water. The HAV 3C cysteine protease is a non-structural protein, plays pivotal role in proliferation and viral replication. Significant phytochemicals of Pandanous fascicularis include phytosterol, kobusin, epipinoresinol, and ceroptene, which have a wide variety of biological functions. Through ADMET investigation, we have screened fifteen phytochemicals for this study. Additionally, using molecular docking, these phytochemicals were docked with the HAV 3C protease which signifies the phytochemicals phytosterol, kobusin, epipinoresinol, and ceroptene have a significant capability to bind with hepatitis A virus protein.The docking study was further accompanied by analyzes RMSD, RMSF, Rg, SASA, H-bond number, and principal component analysis through 100 ns MD simulations. The molecular dynamics study reveals that, all four phytochemicals possess considerable binding efficacy with hepatitis A virus protein. Based on our computational study and MMGBSA calculations, phytosterol, kobusin and epipinoresinol phytochemicals may be a potential drug candidate for inhibition of hepatitis A. The potential therapeutic characteristics of the phytochemicals against hepatitis A inhibition offer additional support for the in vitro and in vivo studies in future.Communicated by Ramaswamy H. Sarma.

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.

Role of lupeol in chemosensitizing therapy-resistant prostate cancer cells by targeting MYC, β-catenin and c-FLIP: in silico and in vitro studies

Prostate cancer (CaP) is one of the most frequent malignancies amongst men. Enzalutamide is the second-generation potent androgen receptor (AR) antagonist used against metastatic and non-metastatic CaP. Unfortunately, the development of chemoresistance in cancer cells reduces the effectiveness of Enzalutamide. Lupeol is a pentacyclic triterpene found in different fruits, vegetables, and medicinal plants and possesses anti-inflammatory and anti-cancer properties. Here, we report in silico and in vitro studies of Lupeol and Enzalutamide against the β-CATENIN, c-FLIPL, and c-MYC, which play a significant role in chemoresistance. We observed that Lupeol significantly inhibits the cell growth of chemoresistant Du145 cells and cancer stem cells (CSCs) either alone or in combination with Enzalutamide. Lupeol and Enzalutamide were also found to dock with β-CATENIN, c-FLIPL, and c-MYC. The following MD simulation data showed both compounds exerting structural changes in these proteins. Finally, they significantly inhibit the transcriptional activity of all these genes, as observed by luciferase assay. Thus, we infer that Lupeol chemosensitizes the CaP cells for Enzalutamide-resistant CaP cells.

Graphical abstract

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.

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.

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.

The allosteric interplay between S-nitrosylation and glycine binding controls the activity of human serine racemase

Human serine racemase (hSR) catalyzes the biosynthesis of D-serine, an obligatory co-agonist of the NMDA receptors. It was previously found that the reversible S-nitrosylation of Cys113 reduces hSR activity. Here, we show by site-directed mutagenesis, fluorescence spectroscopy, mass spectrometry, and molecular dynamics that S-nitrosylation stabilizes an open, less-active conformation of the enzyme. The reaction of hSR with either NO or nitroso donors is conformation-dependent and occurs only in the conformation stabilized by the allosteric effector ATP, in which the ε-amino group of Lys114 acts as a base toward the thiol group of Cys113. In the closed conformation stabilized by glycine-an active-site ligand of hSR-the side chain of Lys114 moves away from that of Cys113, while the carboxyl side-chain group of Asp318 moves significantly closer, increasing the thiol pK_a and preventing the reaction. We conclude that ATP binding, glycine binding, and S-nitrosylation constitute a three-way regulation mechanism for the tight control of hSR activity. We also show that Cys113 undergoes H₂ O₂ -mediated oxidation, with loss of enzyme activity, a reaction also dependent on hSR conformation.

Structural Insights on Tiny Peptide Nucleic Acid (PNA) Analogues of miRNA-34a: An in silico and Experimental Integrated Approach

In the present work, structural features of the interaction between peptide nucleic acid (PNA)-based analogs of the tumor-suppressor microRNA-34a with both its binding sites on MYCN mRNA were investigated. In particular, the region from base 1 to 8 ("seed" region) of miR-34a was reproduced in the form of an 8-mer PNA fragment (tiny PNA), and binding to target 3'UTR MYCN mRNA, was studied by a seldom reported and detailed NMR characterization, providing evidence for the formation of anti-parallel duplexes with a well-organized structural core. The formation of PNA-3'UTR duplexes was also confirmed by Circular Dichroism, and their melting curves were measured by UV spectroscopy. Nevertheless, this study offered a valuable comparison between molecular dynamics predictions and experimental evidence, which showed great correlation. Preliminary uptake assays were carried out in Neuroblastoma Kelly cells, using short peptide conjugates as carriers and FITC fluorescent tag for subcellular localization. Moderate internalization was observed without the use of transfecting agents. The reported results corroborate the interest toward the design and development of chimeric PNA/RNA sequences as effective RNA-targeting agents.

Development of Cu2+-Based Distance Methods and Force Field Parameters for the Determination of PNA Conformations and Dynamics by EPR and MD Simulations

Peptide nucleic acids (PNAs) are a promising group of synthetic analogues of DNA and RNA that offer several distinct advantages over the naturally occurring nucleic acids for applications in biosensing, drug delivery, and nanoelectronics. Because of its structural differences from DNA/RNA, methods to analyze and assess the structure, conformations, and dynamics are needed. In this work, we develop synergistic techniques for the study of the PNA conformation. We use CuQ₂, a Cu²⁺ complex with 8-hydroxyquinoline (HQ), as an alternative base pair and as a spin label in electron paramagnetic resonance (EPR) distance methods. We use molecular dynamics (MD) simulations with newly developed force field parameters for the spin labels to interpret the distance constraints determined by EPR. We complement these methods by UV-vis and circular dichroism measurements and assess the efficacy of the Cu²⁺ label on a PNA duplex whose backbone is based on aminoethylglycine and a duplex with a hydroxymethyl backbone modification. We show that the Cu²⁺ label functions efficiently within the standard PNA and the hydroxymethyl-modified PNA and that the MD parameters may be used to accurately reproduce our EPR findings. Through the combination of EPR and MD, we gain new insights into the PNA structure and conformations as well as into the mechanism of orientational selectivity in Cu²⁺ EPR at X-band. These results present for the first time a rigid Cu²⁺ spin label used for EPR distance measurements in PNA and the accompanying MD force fields for the spin label. Our studies also reveal that the spin labels have a low impact on the structure of the PNA duplexes. The combined MD and EPR approach represents an important new tool for the characterization of the PNA duplex structure and provides valuable information to aid in the rational application of PNA at large.

Significant improvement in catalytic activity and enantioselectivity of a Phaseolus vulgaris epoxide hydrolase, PvEH3, towards ortho-cresyl glycidyl ether based on the semi-rational design

The investigation of substrate spectrum towards five racemic (rac-) aryl glycidyl ethers (1a-5a) indicated that E. coli/pveh3, an E. coli BL21(DE3) transformant harboring a PvEH3-encoding gene pveh3, showed the highest EH activity and enantiomeric ratio (E) towards rac-3a. For efficiently catalyzing the kinetic resolution of rac-3a, the activity and E value of PvEH3 were further improved by site-directed mutagenesis of selected residues. Based on the semi-rational design of an NC-loop in PvEH3, four single-site variants of pveh3 were amplified by PCR, and intracellularly expressed in E. coli BL21(DE3), respectively. E. coli/pveh3^E134K and /pveh3^T137P had the enhanced EH activities of 15.3 ± 0.4 and 16.1 ± 0.5 U/g wet cell as well as E values of 21.7 ± 1.0 and 21.2 ± 1.1 towards rac-3a. Subsequently, E. coli/pveh3^E134K/T137P harboring a double-site variant gene was also constructed, having the highest EH activity of 22.4 ± 0.6 U/g wet cell and E value of 24.1 ± 1.2. The specific activity of the purified PvEH3^E134K/T137P (14.5 ± 0.5 U/mg protein) towards rac-3a and its catalytic efficiency (k_cat/K_m of 5.67 mM^-1 s^-1) for (S)-3a were 1.7- and 3.54-fold those (8.4 ± 0.3 U/mg and 1.60 mM^-1 s^-1) of PvEH3. The gram-scale kinetic resolution of rac-3a using whole wet cells of E. coli/pveh3^E134K/T137P was performed at 20 °C for 7.0 h, producing (R)-3a with 99.4% ee_s and 38.5 ± 1.2% yield. Additionally, the mechanism of PvEH3^E134K/T137P with remarkably improved E value was analyzed by molecular docking simulation.

Antibacterial Peptide Nucleic Acids-Facts and Perspectives

Antibiotic resistance is an escalating, worldwide problem. Due to excessive use of antibiotics, multidrug-resistant bacteria have become a serious threat and a major global healthcare problem of the 21st century. This fact creates an urgent need for new and effective antimicrobials. The common strategies for antibiotic discovery are based on either modifying existing antibiotics or screening compound libraries, but these strategies have not been successful in recent decades. An alternative approach could be to use gene-specific oligonucleotides, such as peptide nucleic acid (PNA) oligomers, that can specifically target any single pathogen. This approach broadens the range of potential targets to any gene with a known sequence in any bacterium, and could significantly reduce the time required to discover new antimicrobials or their redesign, if resistance arises. We review the potential of PNA as an antibacterial molecule. First, we describe the physicochemical properties of PNA and modifications of the PNA backbone and nucleobases. Second, we review the carriers used to transport PNA to bacterial cells. Furthermore, we discuss the PNA targets in antibacterial studies focusing on antisense PNA targeting bacterial mRNA and rRNA.

Molecular Dynamics Study of the Hybridization between RNA and Modified Oligonucleotides

MicroRNAs (miRNAs) are attractive drug candidates for many diseases as they can modulate the expression of gene networks. Recently, we discovered that DNAs targeting microRNA-22-3p (miR-22-3p) hold the potential for treating obesity and related metabolic disorders (type 2 diabetes mellitus, hyperlipidemia, and nonalcoholic fatty liver disease (NAFLD)) by turning fat-storing white adipocytes into fat-burning adipocytes. In this work, we explored the effects of chemical modifications, including phosphorothioate (PS), locked nucleic acid (LNA), and peptide nucleic acid (PNA), on the structure and energy of DNA analogs by using molecular dynamics (MD) simulations. To achieve a reliable prediction of the hybridization free energy, the AMOEBA polarizable force field and the free energy perturbation technique were employed. The calculated hybridization free energies are generally compatible with previous experiments. For LNA and PNA, the enhanced duplex stability can be explained by the preorganization mechanism, i.e., the single strands adopt stable helical structures similar to those in the duplex. For PS, the S and R isomers (Sp and Rp) have preferences for C2'-endo and C3'-endo sugar puckering conformations, respectively, and therefore Sp is less stable than Rp in DNA/RNA hybrids. In addition, the solvation penalty of Rp accounts for its destabilization effect. PS-LNA is similar to LNA as the sugar puckering is dominated by the locked sugar ring. This work demonstrated that MD simulations with polarizable force fields are useful for the understanding and design of modified nucleic acids.

A combined experimental and computational study on peptide nucleic acid (PNA) analogues of tumor suppressive miRNA-34a

MicroRNAs are a ubiquitous class of non-coding RNAs able to regulate gene expression in diverse biological processes. Widespread miRNAs deregulation was reported in numerous diseases including cancer, with several miRNAs playing oncogenic and/or tumor suppressive role by targeting multiple mRNAs simultaneously. Based on these findings, miRNAs have emerged as promising therapeutic tools for cancer treatment. Herein, for the first time, peptide nucleic acids (PNAs) were studied to develop a new class of molecules able to target 3'UTR on MYCN mRNA without a fully complementary base pairing sequence (as miRNAs). For our proof of concept study we have selected as a model the miRNA-34a, which acts as a tumor suppressor in a number of cancers including neuroblastoma. In particular, miRNA-34a is a direct regulator of MYCN oncogene, whose overexpression is a prominent biomarker for the highly aggressive neuroblastoma phenotype. The design and synthesis of three PNA-based oligomers of different length was described, and their interaction with two binding sites on the target MYCN mRNA was investigated by molecular dynamics simulation, and spectroscopic techniques (CD, UV). Intake assay and confocal microscopy of PNA sequences were also carried out in vitro on neuroblastoma Kelly cells. Despite the presence of multiple mismatches, the PNA/RNA hetero duplexes retain very interesting features in terms of stability, affinity as well as of cellular uptake.

Impact of N-(2-aminoethyl) Glycine Unit on Watson-Crick Base Pairs

We aim to understand the structure and stability of the backbone tailored Watson-Crick base pairs, Guanine-Cytosine (GC), Adenine-Thymine (AT) and Adenine-Uracil (AU) by incorporating N-(2-aminoethyl) glycine units (linked by amide bonds) at the purine and pyrimidine sites of the nucleobases. Density functional theory (DFT) is employed in which B3LYP/6-311++G∗ ∗ level of theory has been used to optimize all the structures. The peptide attached base pairs are compared with the natural deoxyribose nucleic acid (DNA)/ribonucleic acid (RNA) base pairs and the calculations are carried out in both the gas and solution phases. The structural propensities of the optimized base pairs are analyzed using base pair geometries, hydrogen bond distances and stabilization energies and, compared with the standard reference data. The structural parameters were found to correlate well with the available data. The addition of peptide chain at the back bone of the DNA/RNA base pairs results only with a minimal distortion and hence does not alter the structural configuration of the base pairs. Also enhanced stability of the base pairs is spotted while adding peptidic chain at the purine site rather than the pyrimidine site of the nucleobases. The stability of the complexes is further interpreted by considering the hydrogen bonded N–H stretching frequencies of the respective base pairs. The discrimination in the interaction energies observed in both gas and solution phases are resulted due to the existence of distinct lowest unoccupied molecular orbitals (LUMO) in the solution phase. The reactivity of the base pairs is also analyzed through the in-depth examinations on the highest occupied molecular orbital (HOMO)-LUMO orbitals.

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.

The intrinsic flexibility of the aptamer targeting the ribosomal protein S8 is a key factor for the molecular recognition

BACKGROUND

Aptamers are RNA/DNA biomolecules representing an emerging class of protein interactors and regulators. Despite the growing interest in these molecules, current understanding of chemical-physical basis of their target recognition is limited. Recently, the characterization of the aptamer targeting the protein-S8 has suggested that flexibility plays important functional roles. We investigated the structural versatility of the S8-aptamer by molecular dynamics simulations.

METHODS

Five different simulations have been conducted by varying starting structures and temperatures.

RESULTS

The simulation of S8-aptamer complex provides a dynamic view of the contacts occurring at the complex interface. The simulation of the aptamer in ligand-free state indicates that its central region is intrinsically endowed with a remarkable flexibility. Nevertheless, none of the trajectory structures adopts the structure observed in the S8-aptamer complex. The aptamer ligand-bound is very rigid in the simulation carried out at 300 K. A structural transition of this state, providing insights into the aptamer-protein recognition process, is observed in a simulation carried out at 400 K. These data indicate that a key event in the binding is linked to the widening of the central region of the aptamer. Particularly relevant is switch of the A26 base from its ligand-free state to a location that allows the G13-C28 base-pairing.

CONCLUSIONS

Intrinsic flexibility of the aptamer is essential for partner recognition. Present data indicate that S8 recognizes the aptamer through an induced-fit rather than a population-shift mechanism.

GENERAL SIGNIFICANCE

The present study provides deeper understanding of the structural basis of the structural versatility of aptamers.

Peptide Nucleic Acids as miRNA Target Protectors for the Treatment of Cystic Fibrosis

Cystic Fibrosis (CF) is one of the most common life shortening conditions in Caucasians. CF is caused by mutations in the CF Transmembrane Conductance Regulator (CFTR) gene which result in reduced or altered CFTR functionality. Several microRNAs (miRNAs) downregulate the expression of CFTR, thus causing or exacerbating the symptoms of CF. In this context, the design of anti-miRNA agents represents a valid functional tool, but its translation to the clinic might lead to unpredictable side effects because of the interference with the expression of other genes regulated by the same miRNAs. Herein, for the first time, is proposed the use of peptide nucleic acids (PNAs) to protect specific sequences in the 3’UTR (untranslated region) of the CFTR messenger RNA (mRNA) by action of miRNAs. Two PNAs (7 and 13 bases long) carrying the tetrapeptide Gly-SerP-SerP-Gly at their C-end, fully complementary to the 3’UTR sequence recognized by miR-509-3p, have been synthesized and the structural features of target PNA/RNA heteroduplexes have been investigated by spectroscopic and molecular dynamics studies. The co-transfection of the pLuc-CFTR-3´UTR vector with different combinations of PNAs, miR-509-3p, and controls in A549 cells demonstrated the ability of the longer PNA to rescue the luciferase activity by up to 70% of the control, thus supporting the use of suitable PNAs to counteract the reduction in the CFTR expression.

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.

Solid phase synthesis of TyrT, a thymine-tyrosine conjugate with poly(A) RNA-binding ability

TyrT nucleoamino amide interacts with poly(A) RNA.

The present work deals with the synthesis and characterization of a novel nucleoamino acid derivative based on a l-tyrosine moiety to which a thymine nucleobase was anchored by means of an amide bond to the N-alpha group. This derivative, denominated by us TyrT, belongs to the family of the nucleobase–amino acid conjugates that show a wide range of biological activities, frequently associated with their ability to interact with nucleic acids. In this respect, the interaction of TyrT with poly(A), a proposed RNA target for anticancer strategies, was studied by circular dichroism (CD) which suggested its ability to bind this RNA. Moreover, the modification of the morphology of a sample of TyrT in the presence of poly(A) was visualised by scanning electron microscopy (SEM) which was in agreement with the evidence that the thyminyl l-tyrosine interacts with poly(A). Finally, computational analyses have been performed to hypothesize the binding mode from a structural point of view, suggesting that the binding is mainly kept via hydrophobic contacts, reproducing a stacking-like interaction between the thymine ring of TyrT and the two successive adenine rings of a poly(A) model.

Effect of backbone chemistry on hybridization thermodynamics of oligonucleic acids: a coarse-grained molecular dynamics simulation study

In this paper we study how varying oligonucleic acid backbone chemistry affects the hybridization/melting thermodynamics of oligonucleic acids. We first describe the coarse-grained (CG) model with tunable parameters that we developed to enable the study of both naturally occurring oligonucleic acids, such as DNA, and their chemically-modified analogues, such as peptide nucleic acids (PNAs) and locked nucleic acids (LNAs). The DNA melting curves obtained using such a CG model and molecular dynamics simulations in an implicit solvent and with explicit ions match with the melting curves obtained using the empirical nearest-neighbor models. We use these CG simulations to then elucidate the effect of backbone flexibility, charge, and nucleobase spacing along the backbone on the melting curves, potential energy and conformational entropy change upon hybridization and base-pair hydrogen bond residence time. We find that increasing backbone flexibility decreases duplex thermal stability and melting temperature mainly due to increased conformational entropy loss upon hybridization. Removing charges from the backbone enhances duplex thermal stability due to the elimination of electrostatic repulsion and as a result a larger energetic gain upon hybridization. Lastly, increasing nucleobase spacing decreases duplex thermal stability due to decreasing stacking interactions that are important for duplex stability.

Programmable nanoscaffolds that control ligand display to a G-protein-coupled receptor in membranes to allow dissection of multivalent effects

A programmable ligand display system can be used to dissect the multivalent effects of ligand binding to a membrane receptor. An antagonist of the A_2A adenosine receptor, a G-protein-coupled receptor that is a drug target for neurodegenerative conditions, was displayed in 35 different multivalent configurations, and binding to A_2A was determined. A theoretical model based on statistical mechanics was developed to interpret the binding data, suggesting the importance of receptor dimers. Using this model, extended multivalent arrangements of ligands were constructed with progressive improvements in binding to A_2A. The results highlight the ability to use a highly controllable multivalent approach to determine optimal ligand valency and spacing that can be subsequently optimized for binding to a membrane receptor. Models explaining the multivalent binding data are also presented.

Conformational studies of chiral D-Lys-PNA and achiral PNA system in binding with DNA or RNA through a molecular dynamics approach

The growing interest in peptide nucleic acid (PNA) oligomers has led to the development of a very wide variety of PNA derivatives. Among others, the introduction of charged chiral groups on a PNA oligomer has proven effective in improving DNA binding ability, complexation direction and cellular uptake. In particular, the introduction of three adjacent chiral monomers based on D-Lys in the middle of the PNA sequence (D-Lys-PNA) has produced noteworthy results in modulating the directionality of the binding with the DNA complementary strand and in mismatch detection. Here, through a molecular dynamics approach, a comparative study has been carried out to investigate the structural properties that drive the interaction of the chiral D-Lys-PNA and the corresponding achiral PNA system with DNA as well as RNA complementary strands, starting from the crystal structure of D-Lys-PNA in complex with DNA. The results obtained complement experimental data and indicate that the binding with the RNA molecule, compared to DNA, is differently affected by the addition of three D-Lys groups on the PNA backbone, suggesting that this modification could be taken into account for the development of new PNA-based molecules able to discriminate between DNA and RNA.