Predicting nearest neighbor free energies of modified RNA with LIE: results for pseudouridine and N1-methylpseudouridine within RNA duplexes.

Pseudouridine (Ψ) and N1-methylpseudouridine (m1Ψ) are among the key modifications in the field of mRNA therapeutics and vaccine research. The accuracy of the design and development of therapeutic RNAs containing such modifications depends on the accuracy of the secondary structure prediction, which in turn depends on the nearest neighbor (NN) thermodynamic parameters for the standard and modified residues. Here, we propose a simple approach based on molecular dynamics simulations and linear interaction energy (LIE) approximation that is able to predict the NN free energy parameters for U-A, Ψ-A and m1Ψ-A pairs in reasonable agreement with the recent experimental reports. We report the NN thermodynamic parameters for different U, Ψ and m1Ψ base pairs, which might be helpful for a deeper understanding of the effect of these modifications in RNA. The predicted NN free energy parameters in this study are able to closely reproduce the folding free energies of duplexes containing internal Ψ for which the thermodynamic data were available. Additionally, we report the predicted folding free energies for the duplexes containing internal m1Ψ.

Analysis of the interaction of a non-canonical twin half-site of Cyclic AMP-Response Element (CRE) with CRE-binding protein.

Differential regulation of a gene having either canonical or non-canonical cyclic AMP response element (CRE) in its promoter is primarily accomplished by its interactions with CREB (cAMP-response element binding protein). The present study aims to delineate the mechanism of the CREB-CRE interactions at the Oncostatin-M (osm) promoter by in vitro and in silico approaches. The non-canonical CREosm consists of two half-CREs separated by a short intervening sequence of 9 base pairs. In this study, in vitro binding assays revealed that out of the two CRE half-sites, the right half-CRE was indispensable for binding of CREB, while the left sequence showed weaker binding ability and specificity. Genome-wide modeling and high throughput free energy calculations for the energy-minimized models containing CREB-CREosm revealed that there was no difference in the binding of CREB to the right half of CREosm site when compared to the entire CREosm. These results were in accordance with the in vitro studies, confirming the indispensable role of the right half-CREosm site in stable complex formation with the CREB protein. Additionally, conversion of the right half-CREosm site to a canonical CRE palindrome showed stronger CREB binding, irrespective of the presence or absence of the left CRE sequence. Thus, the present study establishes an interesting insight into the interaction of CREB with a CRE variant located at the far end of a TATA-less promoter of a cytokine-encoding gene, which in turn could be involved in the regulation of transcription under specific conditions.

LNA-induced dynamic stability in a therapeutic aptamer: insights from molecular dynamics simulations.

Modulation of structural and thermodynamic properties of nucleic acids with synthetic modifications is a promising area of research with possible applications in nanotechnology and nanotherapeutics. Locked nucleic acid (LNA) is one such modification in which the C4' and O2' atoms of the sugar moiety are connected through a methylene bridge. The LNA modified DNA aptamer RNV66, and its unmodified counterpart V7t1, both of which target the vascular endothelial growth factor (VEGF) implicated in oncogenic angiogenesis, have a G-rich tract that can fold into G-quadruplex structures. However, it is not understood why V7t1 has a polymorphic structure while its LNA modified counterpart RNV66 has a unique quadruplex fold with higher nuclease resistance, thermal stability and greater binding affinity for VEGF. In this work, we have performed extensive molecular dynamics simulations of RNV66 and V7t1 to study and compare the structural and dynamic consequences of the insertion of LNAs. It was observed that the increase in dynamic stability was significant in the presence of LNA residues and our protocol for combining different torsional parameters using OL15 for the DNA aptamer and parm99_LNA along with parmbsc0 and ßOL15 for the LNAs nicely reproduced the experimentally observed conformational features of RNV66. Our observations would help in further theoretical studies in understanding the lack of frustration in the folding of the LNA modified aptamer and its higher affinity for VEGF.Communicated by Ramaswamy H. Sarma.

Prevention of inorganic arsenic induced squamous cell carcinoma of the skin in Swiss albino mice by black tea through epigenetic modulation.

Consumption of inorganic Arsenic (iAs) above the safe level may lead to many diseases including cancers of skin. It is known that carcinogenicity of iAs is mediated through generation of excessive reactive oxygen species and polyphenols present in black tea extract (BTE) ameliorate the deleterious effect. Epigenetics also plays vital roles in carcinogenesis. The aim of this paper is to study the influence of iAs on epigenetics and the modulatory effect of BTE. Male Swiss albino mice were divided into three groups, (i) control, (ii) iAs-administered and (iii) iAs + BTE administered. Group (ii) developed invasive squamous cell carcinoma (SCC) of the skin after 330 days, while only hyperplasic and dysplastic changes were observed in group (iii). Expression levels of histone methylation, acetylation marks and several histone methylases, demethylases and acetylases due to iAs were studied; most aberrant expression levels due to iAs were modulated by BTE. JARID1B, a histone demethylase implicated as one of the markers in SCC and a therapeutic target gets upregulated by iAs, but is not influenced by BTE. However, SCC is prevented by BTE. Upregulation of JARID1B by iAs represses H3K4me3; BTE upregulates H3K4me3 without influencing JARID1B expression level. It is known that theaflavin compounds in BTE are transported to the nucleus and interact with histone proteins. in-silico findings in this paper hint that theaflavin compounds present in BTE are very good inhibitors of JARID1B and BTE inhibits its demethylating activity. BTE reverses the epigenetic alterations caused by iAs, thus aids in prevention of SCC.

Data-informed reparameterization of modified RNA and the effect of explicit water models: application to pseudouridine and derivatives.

Pseudouridine is one of the most abundant post-transcriptional modifications in RNA. We have previously shown that the FF99-derived parameters for pseudouridine and some of its naturally occurring derivatives in the AMBER distribution either alone or in combination with the revised γ torsion parameters (parmbsc0) failed to reproduce their conformational characteristics observed experimentally (Deb et al. in J Chem Inf Model 54:1129-1142, 2014; Deb et al. in J Comput Chem 37:1576-1588, 2016; Dutta et al. in J Chem Inf Model 60:4995-5002, 2020). However, the application of the recommended bsc0 correction did lead to an improvement in the description not only of the distribution in the γ torsional space but also of the sugar pucker distributions. In an earlier study, we examined the transferability of the revised glycosidic torsion parameters (χIDRP) for Ψ to its derivatives. We noticed that although these parameters in combination with the AMBER FF99-derived parameters and the revised γ torsional parameters resulted in conformational properties of these residues that were in better agreement with experimental observations, the sugar pucker distributions were still not reproduced accurately. Here we report a new set of partial atomic charges for pseudouridine, 1-methylpseudouridine, 3-methylpseudouridine and 2'-O-methylpseudouridine and a new set of glycosidic torsional parameters (χND) based on chosen glycosidic torsional profiles that most closely corresponded to the NMR data for conformational propensities and studied their effect on the conformational distributions using REMD simulations at the individual nucleoside level. We have also studied the effect of the choice of water model on the conformational characteristics of these modified nucleosides. Our observations suggest that the current revised set of parameters and partial atomic charges describe the sugar pucker distributions for these residues more accurately and that the choice of a suitable water model is important for the accurate description of their conformational properties. We have further validated the revised sets of parameters by studying the effect of substitution of uridine with pseudouridine within single stranded RNA oligonucleotides on their conformational and hydration characteristics.

Inosine and its methyl derivatives: Occurrence, biogenesis, and function in RNA.

Inosine is one of the most common post-transcriptional modifications. Since its discovery, it has been noted for its ability to contribute to non-Watson-Crick interactions within RNA. Rapidly accumulating evidence points to the widespread generation of inosine through hydrolytic deamination of adenosine to inosine by different classes of adenosine deaminases. Three naturally occurring methyl derivatives of inosine, i.e., 1-methylinosine, 2'-O-methylinosine and 1,2'-O-dimethylinosine are currently reported in RNA modification databases. These modifications are expected to lead to changes in the structure, folding, dynamics, stability and functions of RNA. The importance of the modifications is indicated by the strong conservation of the modifying enzymes across organisms. The structure, binding and catalytic mechanism of the adenosine deaminases have been well-studied, but the underlying mechanism of the catalytic reaction is not very clear yet. Here we extensively review the existing data on the occurrence, biogenesis and functions of inosine and its methyl derivatives in RNA. We also included the structural and thermodynamic aspects of these modifications in our review to provide a detailed and integrated discussion on the consequences of A-to-I editing in RNA and the contribution of different structural and thermodynamic studies in understanding its role in RNA. We also highlight the importance of further studies for a better understanding of the mechanisms of the different classes of deamination reactions. Further investigation of the structural and thermodynamic consequences and functions of these modifications in RNA should provide more useful information about their role in different diseases.

Molecular Dynamics Simulation of the Conformational Preferences of Pseudouridine Derivatives: Improving the Distribution in the Glycosidic Torsion Space.

There are only four derivatives of pseudouridine (Ψ) that are known to occur naturally in RNA as post-transcriptional modifications. We have studied the conformational consequences of pseudouridylation and further modifications using replica exchange molecular dynamics simulations at the nucleoside level, and the simulated conformational preferences were compared with the available experimental (NMR) data. We found that the existing AMBER FF99-derived parameters for these nucleosides did not reproduce the observed experimental features and while the recommended bsc0 correction could be combined with these parameters leading to an improvement in the description of sugar pucker distributions, the χOL3 correction could not be applied to these nucleosides as such because of base isomerization. On the other hand, the revised χ torsion parameters (χIDRP) for Ψ developed earlier by us (Deb, I., J. Comput. Chem., 2016, 37, 1576-1588) in combination with the AMBER provided parameters and the revised γ torsion parameters generated conformational distributions, which generally were in better agreement with the experimental data. A significant shift of the distribution of base orientation toward the syn conformation was observed with our revised parameter sets compared to the large excess of anti conformation predicted by the FF99 parameters. Overall, our observations indicated that our revised set of parameters (χIDRP) for Ψ were also able to generate conformational distributions for all of the derivatives of Ψ in better agreement with the experimental data.

Ensemble Allosteric Model for the Modified Wobble Hypothesis.

The residue 2-thiouridine (s2U) provides a remarkable example for the "modified wobble" hypothesis, which postulates that some post-transcriptional modifications at the wobble position of tRNAs restrict recognition of degenerate codons. Through extensive molecular dynamics simulations using our χIDRP force field parameters, we demonstrate how this modification shifts the conformational ensemble from a predominantly disordered, as in the case of an RNA pentamer (GUUUC), to a substantially ordered population in Gs2UUUC. Our simulations clearly showed that the van der Waals interaction of sulfur played a major role in driving the disorder-to-order transition. The conformational redistribution and the slowing down of the transition between the clusters within the population in the presence of s2U suggest ensemble allostery to be a key mechanism that may play a general role in the functioning of the wobble modifications of tRNAs.

Comparative study of the SBP-box gene family in rice siblings.

SBP-box genes are a class of plant-specific transcription factors which have a common DNA-binding domain (SBP-domain) with an unusual zinc-finger architecture. Many of the genes in this class are thought to play a developmental role and a few are involved in the determination of plant architecture. We have made a comparative study of these genes in the genomes of rice (Oryza sativa japonica and Oryza sativa indica) and its nine siblings using a recently proposed hybrid method for orthology and paralogy detection (HyPPO). According to HyPPO, the SBP-box proteins of rice siblings could be divided into twenty primary orthologous groups on the basis of their overall sequence features. This contrasts with a much less number of groups found in earlier work with other plant genomes using phylogenetic analysis of the SBP-domains only. The orthologous groups reported by HyPPO showed close correspondence in exon-intron structure and motif conservation. Comparison between different Oryza species revealed disparity in the maintenance of orthologous genes which may result in their different developmental characteristics. Inclusion of the SBP-box proteins from A. thaliana did not result in any change in the orthologous groups except for the A. thaliana proteins being added to some of the existing groups. The closer correspondence between the proteins in the primary orthologous clusters is expected to help in a more reliable prediction of their functions. It is also expected to provide better insight into the evolutionary history of this class of plant-specific proteins.

Computational and NMR studies of RNA duplexes with an internal pseudouridine-adenosine base pair.

Pseudouridine (Ψ) is the most common chemical modification present in RNA. In general, Ψ increases the thermodynamic stability of RNA. However, the degree of stabilization depends on the sequence and structural context. To explain experimentally observed sequence dependence of the effect of Ψ on the thermodynamic stability of RNA duplexes, we investigated the structure, dynamics and hydration of RNA duplexes with an internal Ψ-A base pair in different nearest-neighbor sequence contexts. The structures of two RNA duplexes containing 5'-GΨC/3'-CAG and 5'-CΨG/3'-GAC motifs were determined using NMR spectroscopy. To gain insight into the effect of Ψ on duplex dynamics and hydration, we performed molecular dynamics (MD) simulations of RNA duplexes with 5'-GΨC/3'-CAG, 5'-CΨG/3'-GAC, 5'-AΨU/3'-UAA and 5'-UΨA/3'-AAU motifs and their unmodified counterparts. Our results showed a subtle impact from Ψ modification on the structure and dynamics of the RNA duplexes studied. The MD simulations confirmed the change in hydration pattern when U is replaced with Ψ. Quantum chemical calculations showed that the replacement of U with Ψ affected the intrinsic stacking energies at the base pair steps depending on the sequence context. The calculated intrinsic stacking energies help to explain the experimentally observed sequence dependent changes in the duplex stability from Ψ modification.

Probing the functional conformations of an atypical proline-rich fusion peptide.

The left-handed polyproline II (PPII) type helical structures are thought to play a very important role in many essential biological processes, particularly in recognition mechanisms. However, reliable characterisation of PPII conformation in solution can be experimentally challenging. Computational simulation of these structures offers an attractive alternative, but the accuracy of the results is dependent on the accuracy of the force field employed. In this report, we present the results of simulation of the structural and dynamical properties of a proline-rich viral fusion peptide for which a solution NMR study reported a substantial stretch of PPII conformation in the central region. The suggested mode of action of the p15 fusion peptide depended on the exposure of the flanking N-terminal hydrophobic residues to solvent thereby facilitating their interaction with the target membrane. Our simulations with a set of four force field and water model combinations consisting of (AMBER ff99SB*-ILDNP + TIP3P), (OPLS-AA + SPC/E), (AMBER ff03ws + TIP4P/2005 water with scaled protein-water interactions) and (CHARMM36m + TIP3P) showed a general agreement with the NMR results for all the four force field and water model combinations. The central region encompassing positions 9-15 showed a large PPII propensity, reduced flexibility and lower conformational entropy. The PPII conformations were stable and satisfied the Burgi-Dunitz criteria without the participation of any significant water bridging interaction. However, comparison of the experimentally observed chemical shifts with the distribution of shifts predicted from the simulated ensembles showed a much better agreement for the CHARMM36m + TIP3P and AMBER ff03ws + TIP4P/2005 combinations. The models based on these two force fields also generated conformations which were in much better agreement with the NMR model than the much more compact structures predicted by the AMBER ff99SB*-ILDNP and OPLS-AA force fields and predicted a substantially larger solvent accessible surface area in accordance with the suggested mechanism of action of the peptide.

Plant Polypeptide Hormone Systemin Prefers Polyproline II Conformation in Solution.

Systemin, an 18 amino-acid-signaling peptide, was the first plant polypeptide hormone to be discovered. Earlier structural studies involving NMR spectroscopy indicated a lack of definite structure in solution while circular dichroism spectroscopy suggested the presence of left-handed polyproline II (PPII) conformation. Here, we report the results of molecular dynamics simulations of the peptide in explicit solvent with two different force fields, namely, ff99SBildn and ff99IDPs, both of which showed a large propensity for PPII-like conformations in spite of showing differing features for other conformational characteristics. More remarkably, the conformations with predicted chemical shifts that agreed better with the NMR observations had a larger than average PPII content, especially for the ff99IDPs force field. An independent docking calculation of the molecule with the putative receptor SR160 also retained this conformational preference for PPII structure. The results suggest PPII to be an important class of conformation for systemin which may have a role in its bioactivity.

Effect of Inactivating Mutations on Peptide Conformational Ensembles: The Plant Polypeptide Hormone Systemin.

As part of their basal immune mechanism against insect/herbivore attacks, plants have evolved systemic response mechanisms. Such a systemic wound response in tomato was found to involve an 18 amino acid polypeptide called systemin, the first polypeptide hormone to be discovered in plants. Systematic alanine scanning and deletion studies showed differential modulation in its activity, particularly a major loss of function due to alanine substitution at positions 13 and 17 and less extentive loss of function due to substitution at position 12. We have studied the conformational ensembles of wild-type systemin along with its 17 variants by carrying out a total of 5.76 µs of replica-exchange molecular dynamics simulation in an implicit solvent environment. In our simulations, wild-type systemin showed a lack of α-helical and ß-sheet structures, in conformity with earlier circular dichroism and NMR data. On the other hand, two regions containing diproline segments showed a tendency to adopt polyproline II structures. Examination of conformational ensembles of the 17 variants revealed a change in the population distributions, suggesting a less flexible structure for alanine substitutions at positions 12 and 13 but not for position 17. Combined with the experimental observations that positions 1-14 of systemin are important for the formation of the peptide-receptor complex, this leads to the hypothesis that loss of conformational flexibility may play a role in the loss of activity of systemin due to the P12A and P13A substitutions, while T17A deactivation probably occurs for a different reason, most likely the loss of the threonine phosphorylation site. We also indicate possible structural reasons why the substitution of the prolines at positions 12 and 13 leads to a loss of conformational freedom in the peptide.

Reparameterizations of the χ Torsion and Lennard-Jones σ Parameters Improve the Conformational Characteristics of Modified Uridines.

UNLABELLED: The currently available force field parameters for modified RNA residues in AMBER show significant deviations in conformational properties from experimental observations. The examination of the transferability of the recently revised torsion parameters revealed that there was an overall improvement in the conformational properties for some of the modifications but the improvements were still insufficient in describing the sugar pucker preferences (J. Chem. Inf. MODEL: 2014, 54, 1129-1142). Here, we report an approach for the development and fine tuning of the AMBER force field parameters for 2-thiouridine, 4-thiouridine, and pseudouridine with diverse conformational preferences. The χ torsion parameters were reparameterized at the individual nucleoside level. The effect of combining the revised γ torsion parameter and modifying the Lennard-Jones σ parameters were also tested by directly comparing the conformational preferences obtained from our extensive molecular dynamics simulations with those from experimental observations. © 2016 Wiley Periodicals, Inc.

Conformational preferences of modified uridines: comparison of AMBER derived force fields.

The widespread occurrence of modified residues in RNA sequences necessitates development of accurate parameters for these modifications for reliable modeling of RNA structure and dynamics. A comprehensive set of parameters for the 107 naturally occurring RNA modifications was proposed by Aduri et al. (J. Chem. Theory Comput. 2007, 3, 1464-1475) for the AMBER FF99 force field. In this work, we tested these parameters on a set of modified uridine residues, namely, dihydrouridine, 2-thiouridine, 4-thiouridine, pseudouridine, and uridine-5-oxyacetic acid, by performing molecular dynamics and replica exchange molecular dynamics simulations of these nucleosides. Although our simulations using the FF99 force field did not, in general, reproduce the experimentally observed conformational characteristics well, combination of the parameter set with recent revisions of the FF99 force field for RNA showed noticeable improvement for some of the nucleosides.

Rapid communication capturing the destabilizing effect of dihydrouridine through molecular simulations.

The structural effects of the commonly occurring modified nucleoside dihydrouridine (D) observed experimentally in model oligonucleotides include a strong destabilization of the C3'-endo sugar conformation of D, the disruption of stacking interactions of neighboring residues with D and a possible destabilization of the C3'-endo sugar pucker of the 5'-neighboring nucleoside. Our simulations with a combination of a set of parameters for modified RNA residues with the recently developed AMBER FF99χ force field having reoptimized glycosidic torsion angle parameters for standard nucleosides was found to reproduce the destabilizing effect of dihydrouridine better than with the AMBER FF99 force field for nucleic acids for which the parameters for the modified residues were originally developed.

Specificity determinants for the abscisic acid response element.

Abscisic acid (ABA) response elements (ABREs) are a group of cis-acting DNA elements that have been identified from promoter analysis of many ABA-regulated genes in plants. We are interested in understanding the mechanism of binding specificity between ABREs and a class of bZIP transcription factors known as ABRE binding factors (ABFs). In this work, we have modeled the homodimeric structure of the bZIP domain of ABRE binding factor 1 from Arabidopsis thaliana (AtABF1) and studied its interaction with ACGT core motif-containing ABRE sequences. We have also examined the variation in the stability of the protein-DNA complex upon mutating ABRE sequences using the protein design algorithm FoldX. The high throughput free energy calculations successfully predicted the ability of ABF1 to bind to alternative core motifs like GCGT or AAGT and also rationalized the role of the flanking sequences in determining the specificity of the protein-DNA interaction.

Genome wide gene expression regulation by HIP1 Protein Interactor, HIPPI: prediction and validation.

BACKGROUND: HIP1 Protein Interactor (HIPPI) is a pro-apoptotic protein that induces Caspase8 mediated apoptosis in cell. We have shown earlier that HIPPI could interact with a specific 9 bp sequence motif, defined as the HIPPI binding site (HBS), present in the upstream promoter of Caspase1 gene and regulate its expression. We also have shown that HIPPI, without any known nuclear localization signal, could be transported to the nucleus by HIP1, a NLS containing nucleo-cytoplasmic shuttling protein. Thus our present work aims at the investigation of the role of HIPPI as a global transcription regulator. RESULTS: We carried out genome wide search for the presence of HBS in the upstream sequences of genes. Our result suggests that HBS was predominantly located within 2 Kb upstream from transcription start site. Transcription factors like CREBP1, TBP, OCT1, EVI1 and P53 half site were significantly enriched in the 100 bp vicinity of HBS indicating that they might co-operate with HIPPI for transcription regulation. To illustrate the role of HIPPI on transcriptome, we performed gene expression profiling by microarray. Exogenous expression of HIPPI in HeLa cells resulted in up-regulation of 580 genes (p < 0.05) while 457 genes were down-regulated. Several transcription factors including CBP, REST, C/EBP beta were altered by HIPPI in this study. HIPPI also interacted with P53 in the protein level. This interaction occurred exclusively in the nuclear compartment and was absent in cells where HIP1 was knocked down. HIPPI-P53 interaction was necessary for HIPPI mediated up-regulation of Caspase1 gene. Finally, we analyzed published microarray data obtained with post mortem brains of Huntington's disease (HD) patients to investigate the possible involvement of HIPPI in HD pathogenesis. We observed that along with the transcription factors like CREB, P300, SREBP1, Sp1 etc. which are already known to be involved in HD, HIPPI binding site was also significantly over-represented in the upstream sequences of genes altered in HD. CONCLUSIONS: Taken together, the results suggest that HIPPI could act as an important transcription regulator in cell regulating a vast array of genes, particularly transcription factors and at least, in part, play a role in transcription deregulation observed in HD.

Arabidopsis thaliana regulatory element analyzer.

UNLABELLED: In the Arabidopsis thaliana regulatory element analyzer (AtREA) server, we have integrated sequence data, genome-wide expression data and functional annotation data in three application modules which will be useful to identify major regulatory targets of a user-provided cis-regulatory element (CRE), study different features of CRE distribution and evaluate the role of a set of CREs in the regulation of gene expression--independently as well as in combination with other user-provided CREs. AVAILABILITY: AtREA is freely available at http://www.bioinformatics.org/grn/atrea.html.