Structural dynamics of the box C/D RNA kink-turn and its complex with proteins: the role of the A-minor 0 interaction, long-residency water bridges, and structural ion-binding sites revealed by molecular simulations

Sensitivity of the RNA Structure to Ion Conditions as Probed by Molecular Dynamics Simulations of Common Canonical RNA Duplexes

RNA molecules play a key role in countless biochemical processes. RNA interactions, which are of highly diverse nature, are determined by the fact that RNA is a highly negatively charged polyelectrolyte, which leads to intimate interactions with an ion atmosphere. Although RNA molecules are formally single-stranded, canonical (Watson-Crick) duplexes are key components of folded RNAs. A double-stranded (ds) RNA is also important for the design of RNA-based nanostructures and assemblies. Despite the fact that the description of canonical dsRNA is considered the least problematic part of RNA modeling, the imperfect shape and flexibility of dsRNA can lead to imbalances in the simulations of larger RNAs and RNA-containing assemblies. We present a comprehensive set of molecular dynamics (MD) simulations of four canonical A-RNA duplexes. Our focus was directed toward the characterization of the influence of varying ion concentrations and of the size of the solvation box. We compared several water models and four RNA force fields. The simulations showed that the A-RNA shape was most sensitive to the RNA force field, with some force fields leading to a reduced inclination of the A-RNA duplexes. The ions and water models played a minor role. The effect of the box size was negligible, and even boxes with a small fraction of the bulk solvent outside the RNA hydration sphere were sufficient for the simulation of the dsRNA.

Base-specific RNA force field improving the dynamics conformation of nucleotide

RNA plays a key role in numerous biological processes. Traditional experimental methods have difficulties capturing the structure and dynamic conformation of RNA. Thus, Molecular dynamic simulations (MDs) has become an essential complementary for RNA experiment. However, state-of-the-art RNA force fields have two major limitations of overestimation base stacking propensity and generation of a high ratio of intercalated conformations. Therefore, a two-step strategy was used to optimize the parameters of ff99bsc0χOL3 (named BSFF1) to improve these limitations, which as well adjusted the unbonded parameters of nucleobase heavy atoms and added ζ/α grid-based energy correction map energy term with reweighting. MD simulations of tetranucleotides indicate that BSFF1 can significantly decrease the ratio of intercalated conformations. Tests of single-strand RNA and kink-turn show that BSFF1 force field can reproduce more accurate conformers than ff99bsc0χOL3 force field. BSFF1 can also stabilize the conformers of duplex and riboswitch. The successful ab initio folding of tetraloop further supports the performance of BSFF1. These findings confirm that the newly developed force field BSFF1 can improve the conformer sampling of RNA.

Plasma SNORD83A as a potential biomarker for early diagnosis of non-small-cell lung cancer

Aim: This study aimed to access the efficacy of plasma small nucleolar RNAs in early diagnosis of non-small-cell lung cancer (NSCLC). Methods: SNORD83A was selected based on databases and further verified in 48 paired formalin-fixed, paraffin-embedded tissues, as well as in plasma from 150 NSCLC patients and 150 healthy donors. The diagnostic efficiency of plasma SNORD83A, as well as in combination with carcinoembryonic antigen, was determined by receiver operating characteristic analysis. Results: SNORD83A was significantly increased not only in tissues but also in plasma from NSCLC patients compared with those from healthy donors. Plasma SNORD83A was able to act as a diagnostic biomarker for NSCLC. The diagnostic efficiency of carcinoembryonic antigen was also significantly elevated for early-stage NSCLC when combined with SNORD83A. Conclusion: SNORD83A can serve as a diagnostic biomarker for NSCLC.

Structure and Mechanical Stabilities of the Three-Way Junction Motifs in Prohead RNA

The core structure of phi29 prohead RNA (pRNA) is composed of three major helices organized into three-way junction pRNA (3WJ-pRNA) and has stout structural rigidity along the coaxial helices. Prohead RNAs of the other Bacillus subtilis bacteriophages such as GA1 and SF5 share similar secondary structure and function with phi29; whether these pRNAs have similar mechanical rigidity remains to be elucidated. In this study, we constructed the tertiary structures of GA1 and SF5 3WJ-pRNAs by comparative modeling. Both GA1 and SF5 3WJ-pRNAs adopt a similar structure, in which three helices are organized as the three-way junction and two of the three helices are stacked coaxially. Moreover, detailed structural features of GA1 and SF5 3WJ-pRNAs are also similar to those of phi29 3WJ-pRNA: all of the bases of the coaxial helices are paired, and all of the adenines in the junction region are paired, which eliminates the interference of A-minor tertiary interactions. Hence, the coaxial helices tightly join to each other, and the major groove between them is very narrow. Two Mg²⁺ ions can thus fit into this major groove and form double Mg clamps. A steered molecular dynamics simulation was used to study the mechanical properties of these 3WJ-pRNAs. Both GA1 and SF5 3WJ-pRNAs show strong resistance to applied force in the direction of their coaxial helices. Such mechanical stability can be attributed to the Mg clamps.

Computational Approaches for Cancer-Fighting: From Gene Expression to Functional Foods

It is today widely accepted that a healthy diet is very useful to prevent the risk for cancer or its deleterious effects. Nutrigenomics studies are therefore taking place with the aim to test the effects of nutrients at molecular level and contribute to the search for anti-cancer treatments. These efforts are expanding the precious source of information necessary for the selection of natural compounds useful for the design of novel drugs or functional foods. Here we present a computational study to select new candidate compounds that could play a role in cancer prevention and care. Starting from a dataset of genes that are co-expressed in programmed cell death experiments, we investigated on nutrigenomics treatments inducing apoptosis, and searched for compounds that determine the same expression pattern. Subsequently, we selected cancer types where the genes showed an opposite expression pattern and we confirmed that the apoptotic/nutrigenomics expression trend had a significant positive survival in cancer-affected patients. Furthermore, we considered the functional interactors of the genes as defined by public protein-protein interaction data, and inferred on their involvement in cancers and/or in programmed cell death. We identified 7 genes and, from available nutrigenomics experiments, 6 compounds effective on their expression. These 6 compounds were exploited to identify, by ligand-based virtual screening, additional molecules with similar structure. We checked for ADME criteria and selected 23 natural compounds representing suitable candidates for further testing their efficacy in apoptosis induction. Due to their presence in natural resources, novel drugs and/or the design of functional foods are conceivable from the presented results.

High-resolution structure of eukaryotic Fibrillarin interacting with Nop56 amino-terminal domain

Ribosomal RNA (rRNA) carries extensive 2'-O-methyl marks at functionally important sites. This simple chemical modification is thought to confer stability, promote RNA folding, and contribute to generate a heterogenous ribosome population with a yet-uncharacterized function. 2'-O-methylation occurs both in archaea and eukaryotes and is accomplished by the Box C/D RNP enzyme in an RNA-guided manner. Extensive and partially conflicting structural information exists for the archaeal enzyme, while no structural data is available for the eukaryotic enzyme. The yeast Box C/D RNP consists of a guide RNA, the RNA-primary binding protein Snu13, the two scaffold proteins Nop56 and Nop58, and the enzymatic module Nop1. Here we present the high-resolution structure of the eukaryotic Box C/D methyltransferase Nop1 from Saccharomyces cerevisiae bound to the amino-terminal domain of Nop56. We discuss similarities and differences between the interaction modes of the two proteins in archaea and eukaryotes and demonstrate that eukaryotic Nop56 recruits the methyltransferase to the Box C/D RNP through a protein-protein interface that differs substantially from the archaeal orthologs. This study represents a first achievement in understanding the evolution of the structure and function of these proteins from archaea to eukaryotes.

Small nucleolar RNAs signature (SNORS) identified clinical outcome and prognosis of bladder cancer (BLCA)

Background

Small nucleolar RNAs (snoRNAs) are a new non-coding RNAs (ncRNAs), which have not been widely investigated and are identified to be involved in tumorigenesis. But the function of snoRNAs in BLCA has not been reported yet.

Methods

SnoRNAs signature (SNORS) was constructed through LASSO cox regression analysis. Integrated analysis of candidate snoRNAs was performed to detect the correlation between copy number variation (CNV)/DNA methylation/protein/mRNA/alternative splicing (AS). Then we built a nomogram integrating independent prognostic factors to assist the clinical utility.

Results

We have screened out 15 prognostic differentially expressed snoRNAs (DESs) and constructed SNORS consisting of 5 candidate snoRNAs which could appropriately stratify patients into low or high SNORS groups with distinct prognosis. Then we found 5 candidate snoRNAs might be regulated by their own CNV and DNA methylation. Moreover, 5 candidate snoRNAs were significantly correlated mRNA and alternative splicing (AS), which might regulate diverse biological process in tumorigenesis, such as "extracellular matrix", "epithelial-mesenchymal transition (EMT)", etc. signaling pathways. Furthermore, SNORS was an independent prognostic factor, which was strikingly correlated with clinical outcome. Through inporating with other variables, we have established a predictive nomogram, which was more effectively to predict prognosis than any other variables alone.

Conclusion

Our findings first highlighted an important role of snoRNAs in BLCA and established a potential prognostic model which could serve as a biomarker for BLCA.

Are Small Nucleolar RNAs "CRISPRable"? A Report on Box C/D Small Nucleolar RNA Editing in Human Cells

CRISPR technologies are nowadays widely used for targeted knockout of numerous protein-coding genes and for the study of various processes and metabolic pathways in human cells. Most attention in the genome editing field is now focused on the cleavage of protein-coding genes or genes encoding long non-coding RNAs (lncRNAs), while the studies on targeted knockout of intron-encoded regulatory RNAs are sparse. Small nucleolar RNAs (snoRNAs) present a class of non-coding RNAs encoded within the introns of various host genes and involved in post-transcriptional maturation of ribosomal RNAs (rRNAs) in eukaryotic cells. Box C/D snoRNAs direct 2'-O-methylation of rRNA nucleotides. These short RNAs have specific elements in their structure, namely, boxes C and D, and a target-recognizing region. Here, we present the study devoted to CRISPR/Cas9-mediated editing of box C/D snoRNA genes in Gas5. We obtained monoclonal cell lines carrying mutations in snoRNA genes and analyzed the levels of the mutant box C/D snoRNA as well as the 2'-O-methylation status of the target rRNA nucleotide in the obtained cells. Mutations in SNORD75 in the obtained monoclonal cell line were shown to result in aberrant splicing of Gas5 with exclusion of exons 3 to 5, which was confirmed by RT-PCR and RNA-Seq. The obtained results suggest that SNORD75 contains an element for binding of some factors regulating maturation of Gas5 pre-lncRNA. We suggest that METTL3/METTL14 is among such factors, and m⁶A-methylation pathways are involved in regulation of Gas5 splicing. Our results shell light on the role of SNORDs in regulating splicing of the host gene.

A Pan-cancer Analysis of the Expression and Clinical Relevance of Small Nucleolar RNAs in Human Cancer

Increasing evidence has demonstrated that small nucleolar RNAs (snoRNAs) play important roles in tumorigenesis. We systematically investigated the expression landscape and clinical relevance of snoRNAs in >10,000 samples across 31 cancer types from The Cancer Genome Atlas. We observed overall elevated expression of snoRNAs and their ribonucleoproteins in multiple cancer types. We showed complex regulation of snoRNA expression by their host genes, copy number variation, and DNA methylation. Unsupervised clustering revealed that the snoRNA expression subtype is highly concordant with other molecular/clinical subtypes. We further identified 46 clinically relevant snoRNAs and experimentally demonstrated functional roles of SNORD46 in promoting cell proliferation, migration, and invasion. We developed a user-friendly data portal, SNORic, to benefit the research community. Our study highlights the significant roles of snoRNAs in the development and implementation of biomarkers or therapeutic targets for cancer and provides a valuable resource for cancer research.

RNA Structural Dynamics As Captured by Molecular Simulations: A Comprehensive Overview

With both catalytic and genetic functions, ribonucleic acid (RNA) is perhaps the most pluripotent chemical species in molecular biology, and its functions are intimately linked to its structure and dynamics. Computer simulations, and in particular atomistic molecular dynamics (MD), allow structural dynamics of biomolecular systems to be investigated with unprecedented temporal and spatial resolution. We here provide a comprehensive overview of the fast-developing field of MD simulations of RNA molecules. We begin with an in-depth, evaluatory coverage of the most fundamental methodological challenges that set the basis for the future development of the field, in particular, the current developments and inherent physical limitations of the atomistic force fields and the recent advances in a broad spectrum of enhanced sampling methods. We also survey the closely related field of coarse-grained modeling of RNA systems. After dealing with the methodological aspects, we provide an exhaustive overview of the available RNA simulation literature, ranging from studies of the smallest RNA oligonucleotides to investigations of the entire ribosome. Our review encompasses tetranucleotides, tetraloops, a number of small RNA motifs, A-helix RNA, kissing-loop complexes, the TAR RNA element, the decoding center and other important regions of the ribosome, as well as assorted others systems. Extended sections are devoted to RNA-ion interactions, ribozymes, riboswitches, and protein/RNA complexes. Our overview is written for as broad of an audience as possible, aiming to provide a much-needed interdisciplinary bridge between computation and experiment, together with a perspective on the future of the field.

Structural study of the Fox-1 RRM protein hydration reveals a role for key water molecules in RRM-RNA recognition

The Fox-1 RNA recognition motif (RRM) domain is an important member of the RRM protein family. We report a 1.8 Å X-ray structure of the free Fox-1 containing six distinct monomers. We use this and the nuclear magnetic resonance (NMR) structure of the Fox-1 protein/RNA complex for molecular dynamics (MD) analyses of the structured hydration. The individual monomers of the X-ray structure show diverse hydration patterns, however, MD excellently reproduces the most occupied hydration sites. Simulations of the protein/RNA complex show hydration consistent with the isolated protein complemented by hydration sites specific to the protein/RNA interface. MD predicts intricate hydration sites with water-binding times extending up to hundreds of nanoseconds. We characterize two of them using NMR spectroscopy, RNA binding with switchSENSE and free-energy calculations of mutant proteins. Both hydration sites are experimentally confirmed and their abolishment reduces the binding free-energy. A quantitative agreement between theory and experiment is achieved for the S155A substitution but not for the S122A mutant. The S155 hydration site is evolutionarily conserved within the RRM domains. In conclusion, MD is an effective tool for predicting and interpreting the hydration patterns of protein/RNA complexes. Hydration is not easily detectable in NMR experiments but can affect stability of protein/RNA complexes.

It's a loop world - single strands in RNA as structural and functional elements

Unpaired regions in RNA molecules - loops - are centrally involved in defining the characteristic three-dimensional (3D) architecture of RNAs and are of high interest in RNA engineering and design. Loops adopt diverse, but specific conformations stabilised by complex tertiary structural interactions that provide structural flexibility to RNA structures that would otherwise not be possible if they only consisted of the rigid A-helical shapes usually formed by canonical base pairing. By participating in sequence-non-local contacts, they furthermore contribute to stabilising the overall fold of RNA molecules. Interactions between RNAs and other nucleic acids, proteins, or small molecules are also generally mediated by RNA loop structures. Therefore, the function of an RNA molecule is generally dependent on its loops. Examples include intermolecular interactions between RNAs as part of the microRNA processing pathways, ribozymatic activity, or riboswitch-ligand interactions. Bioinformatics approaches have been successfully applied to the identification of novel RNA structural motifs including loops, local and global RNA 3D structure prediction, and structural and conformational analysis of RNAs and have contributed to a better understanding of the sequence-structure-function relationships in RNA loops.

Regulatory role of small nucleolar RNAs in human diseases

Small nucleolar RNAs (snoRNAs) are appreciable players in gene expression regulation in human cells. The canonical function of box C/D and box H/ACA snoRNAs is posttranscriptional modification of ribosomal RNAs (rRNAs), namely, 2'-O-methylation and pseudouridylation, respectively. A series of independent studies demonstrated that snoRNAs, as well as other noncoding RNAs, serve as the source of various short regulatory RNAs. Some snoRNAs and their fragments can also participate in the regulation of alternative splicing and posttranscriptional modification of mRNA. Alterations in snoRNA expression in human cells can affect numerous vital cellular processes. SnoRNA level in human cells, blood serum, and plasma presents a promising target for diagnostics and treatment of human pathologies. Here we discuss the relation between snoRNAs and oncological, neurodegenerative, and viral diseases and also describe changes in snoRNA level in response to artificial stress and some drugs.

Can We Execute Stable Microsecond-Scale Atomistic Simulations of Protein-RNA Complexes?

We report over 30 μs of unrestrained molecular dynamics simulations of six protein-RNA complexes in explicit solvent. We utilize the AMBER ff99bsc0χ(OL3) RNA force field combined with the ff99SB protein force field and its more recent ff12SB version with reparametrized side-chain dihedrals. The simulations show variable behavior, ranging from systems that are essentially stable to systems with progressive deviations from the experimental structure, which we could not stabilize anywhere close to the starting experimental structure. For some systems, microsecond-scale simulations are necessary to achieve stabilization after initial sizable structural perturbations. The results show that simulations of protein-RNA complexes are challenging and every system should be treated individually. The simulations are affected by numerous factors, including properties of the starting structures (the initially high force field potential energy, resolution limits, conformational averaging, crystal packing, etc.), force field imbalances, and real flexibility of the studied systems. These factors, and thus the simulation behavior, differ from system to system. The structural stability of simulated systems does not correlate with the size of buried interaction surface or experimentally determined binding affinities but reflects the type of protein-RNA recognition. Protein-RNA interfaces involving shape-specific recognition of RNA are more stable than those relying on sequence-specific RNA recognition. The differences between the protein force fields are considerably smaller than the uncertainties caused by sampling and starting structures. The ff12SB improves description of the tyrosine side-chain group, which eliminates some problems associated with tyrosine dynamics.

Molecular dynamic simulations of protein/RNA complexes: CRISPR/Csy4 endoribonuclease

BACKGROUND

Many prokaryotic genomes comprise Clustered Regularly Interspaced Short Palindromic Repeats (CRISPRs) offering defense against foreign nucleic acids. These immune systems are conditioned by the production of small CRISPR-derived RNAs matured from long RNA precursors. This often requires a Csy4 endoribonuclease cleaving the RNA 3'-end.

METHODS

We report extended explicit solvent molecular dynamic (MD) simulations of Csy4/RNA complex in precursor and product states, based on X-ray structures of product and inactivated precursor (55 simulations; ~3.7μs in total).

RESULTS

The simulations identify double-protonated His29 and deprotonated terminal phosphate as the likely dominant protonation states consistent with the product structure. We revealed potential substates consistent with Ser148 and His29 acting as the general base and acid, respectively. The Ser148 could be straightforwardly deprotonated through solvent and could without further structural rearrangements deprotonate the nucleophile, contrasting similar studies investigating the general base role of nucleobases in ribozymes. We could not locate geometries consistent with His29 acting as general base. However, we caution that the X-ray structures do not always capture the catalytically active geometries and then the reactive structures may be unreachable by the simulation technique.

CONCLUSIONS

We identified potential catalytic arrangement of the Csy4/RNA complex but we also report limitations of the simulation technique. Even for the dominant protonation state we could not achieve full agreement between the simulations and the structural data.

GENERAL SIGNIFICANCE

Potential catalytic arrangement of the Csy4/RNA complex is found. Further, we provide unique insights into limitations of simulations of protein/RNA complexes, namely, the influence of the starting experimental structures and force field limitations. This article is part of a Special Issue entitled Recent developments of molecular dynamics.

Interconversion between parallel and antiparallel conformations of a 4H RNA junction in domain 3 of foot-and-mouth disease virus IRES captured by dynamics simulations

RNA junctions are common secondary structural elements present in a wide range of RNA species. They play crucial roles in directing the overall folding of RNA molecules as well as in a variety of biological functions. In particular, there has been great interest in the dynamics of RNA junctions, including conformational pathways of fully base-paired 4-way (4H) RNA junctions. In such constructs, all nucleotides participate in one of the four double-stranded stem regions, with no connecting loops. Dynamical aspects of these 4H RNAs are interesting because frequent interchanges between parallel and antiparallel conformations are thought to occur without binding of other factors. Gel electrophoresis and single-molecule fluorescence resonance energy transfer experiments have suggested two possible pathways: one involves a helical rearrangement via disruption of coaxial stacking, and the other occurs by a rotation between the helical axes of coaxially stacked conformers. Employing molecular dynamics simulations, we explore this conformational variability in a 4H junction derived from domain 3 of the foot-and-mouth disease virus internal ribosome entry site (IRES); this junction contains highly conserved motifs for RNA-RNA and RNA-protein interactions, important for IRES activity. Our simulations capture transitions of the 4H junction between parallel and antiparallel conformations. The interconversion is virtually barrier-free and occurs via a rotation between the axes of coaxially stacked helices with a transient perpendicular intermediate. We characterize this transition, with various interhelical orientations, by pseudodihedral angle and interhelical distance measures. The high flexibility of the junction, as also demonstrated experimentally, is suitable for IRES activity. Because foot-and-mouth disease virus IRES structure depends on long-range interactions involving domain 3, the perpendicular intermediate, which maintains coaxial stacking of helices and thereby consensus primary and secondary structure information, may be beneficial for guiding the overall organization of the RNA system in domain 3.

Kink turn sRNA folding upon L7Ae binding using molecular dynamics simulations

The kink-turn sRNA motif in archaea, whose combination with protein L7Ae initializes the assembly of small ribonucleoprotein particles (sRNPs), plays a key role in ribosome maturation and the translation process. Although many studies have been reported on this motif, the mechanism of sRNA folding coupled with protein binding is still poorly understood. Here, room and high temperature molecular dynamics (MD) simulations were performed on the complex of 25-nt kink-turn sRNA and L7Ae. The average RMSD values between the bound and corresponding apo structures and Kolmogorov-Smirnov P test analysis indicate that sRNA may follow an induced fit mechanism upon binding with L7Ae, both locally and globally. These conclusions are further supported by high-temperature unfolding kinetic analysis. Principal component analysis (PCA) found both closing and opening motions of the kink-turn sRNA. This might play a key role in the sRNP assembly and methylation catalysis. These combined computational methods can be used to study the specific recognition of other sRNAs and proteins.

Bioinformatics and molecular dynamics simulation study of L1 stalk non-canonical rRNA elements: kink-turns, loops, and tetraloops

The L1 stalk is a prominent mobile element of the large ribosomal subunit. We explore the structure and dynamics of its non-canonical rRNA elements, which include two kink-turns, an internal loop, and a tetraloop. We use bioinformatics to identify the L1 stalk RNA conservation patterns and carry out over 11.5 μs of MD simulations for a set of systems ranging from isolated RNA building blocks up to complexes of L1 stalk rRNA with the L1 protein and tRNA fragment. We show that the L1 stalk tetraloop has an unusual GNNA or UNNG conservation pattern deviating from major GNRA and YNMG RNA tetraloop families. We suggest that this deviation is related to a highly conserved tertiary contact within the L1 stalk. The available X-ray structures contain only UCCG tetraloops which in addition differ in orientation (anti vs syn) of the guanine. Our analysis suggests that the anti orientation might be a mis-refinement, although even the anti interaction would be compatible with the sequence pattern and observed tertiary interaction. Alternatively, the anti conformation may be a real substate whose population could be pH-dependent, since the guanine syn orientation requires protonation of cytosine in the tertiary contact. In absence of structural data, we use molecular modeling to explore the GCCA tetraloop that is dominant in bacteria and suggest that the GCCA tetraloop is structurally similar to the YNMG tetraloop. Kink-turn Kt-77 is unusual due to its 11-nucleotide bulge. The simulations indicate that the long bulge is a stalk-specific eight-nucleotide insertion into consensual kink-turn only subtly modifying its structural dynamics. We discuss a possible evolutionary role of helix H78 and a mechanism of L1 stalk interaction with tRNA. We also assess the simulation methodology. The simulations provide a good description of the studied systems with the latest bsc0χOL3 force field showing improved performance. Still, even bsc0χOL3 is unable to fully stabilize an essential sugar-edge H-bond between the bulge and non-canonical stem of the kink-turn. Inclusion of Mg(2+) ions may deteriorate the simulations. On the other hand, monovalent ions can in simulations readily occupy experimental Mg(2+) binding sites.

Hyperphenylalaninemia in the Czech Republic: genotype-phenotype correlations and in silico analysis of novel missense mutations

BACKGROUND

Hyperphenylalaninemia (HPA) is one of the most common inherited metabolic disorders caused by deficiency of the enzyme phenylalanine hydroxylase (PAH). HPA is associated with mutations in the PAH gene, which leads to reduced protein stability and/or impaired catalytic function. Currently, almost 700 different disease-causing mutations have been described. The impact of mutations on enzyme activity varies ranging from classical PKU, mild PKU, to non-PKU HPA phenotype.

METHODS

We provide results of molecular genetic diagnostics of 665 Czech unrelated HPA patients, structural analysis of missense mutations associated with classical PKU and non-PKU HPA phenotype, and prediction of effects of 6 newly discovered HPA missense mutations using bioinformatic approaches and Molecular Dynamics simulations.

RESULTS

Ninety-eight different types of mutations were indentified. Thirteen of these were novel (6 missense, 2 nonsense, 1 splicing, and 4 small gene rearrangements). Structural analysis revealed that classical PKU mutations are more non-conservative compared to non-PKU HPA mutations and that specific sequence and structural characteristics of a mutation might be critical when distinguishing between non-PKU HPA and classical PKU mutations. The greatest impact was predicted for the p.(Phe263Ser) mutation while other novel mutations p.(Asn167Tyr), p.(Thr200Asn), p.(Asp229Gly), p.(Leu358Phe), and p.(Ile406Met) were found to be less deleterious.

A-minor tertiary interactions in RNA kink-turns. Molecular dynamics and quantum chemical analysis

The RNA kink-turn is an important recurrent RNA motif, an internal loop with characteristic consensus sequence forming highly conserved three-dimensional structure. Functional arrangement of RNA kink-turns shows a sharp bend in the phosphodiester backbone. Among other signature interactions, kink-turns form A-minor interaction between their two stems. Most kink-turns possess extended A-minor I (A-I) interaction where adenine of the second A•G base pair of the NC-stem interacts with the first canonical pair of the C-stem (i.e., the receptor pair) via trans-sugar-edge/sugar-edge (tSS) and cis-sugar-edge/sugar-edge (cSS) interactions. The remaining kink-turns have less compact A-minor 0 (A-0) interaction with just one tSS contact. We show that kink-turns with A-I in ribosomal X-ray structures keep G═C receptor base pair during evolution while the inverted pair (C═G) is not realized. In contrast, kink-turns with A-0 in the observed structures alternate G═C and C═G base pairs in sequences. We carried out an extended set (~5 μs) of explicit-solvent molecular dynamics simulations of kink-turns to rationalize this structural/evolutionary pattern. The simulations were done using a net-neutral Na(+) cation atmosphere (with ~0.25 M cation concentration) supplemented by simulations with either excess salt KCl atmosphere or inclusion of Mg(2+). The results do not seem to depend on the treatment of ions. The simulations started with X-ray structures of several kink-turns while we tested the response of the simulated system to base substitutions, modest structural perturbations and constraints. The trends seen in the simulations reveal that the A-I/G═C arrangement is preferred over all three other structures. The A-I/C═G triple appears structurally entirely unstable, consistent with the covariation patterns seen during the evolution. The A-0 arrangements tend to shift toward the A-I pattern in simulations, which suggests that formation of the A-0 interaction is likely supported by the surrounding protein and RNA molecules. A-0 may also be stabilized by additional kink-turn nucleotides not belonging to the kink-turn consensus, as shown for the kink-turn from ribosomal Helix 15. Quantum-chemical calculations on all four A-minor triples suggest that there is a different balance of electrostatic and dispersion stabilization in the A-I/G═C and A-I/C═G triples, which may explain different behavior of these otherwise isosteric triples in the context of kink-turns.

Understanding RNA Flexibility Using Explicit Solvent Simulations: The Ribosomal and Group I Intron Reverse Kink-Turn Motifs

Reverse kink-turn is a recurrent elbow-like RNA building block occurring in the ribosome and in the group I intron. Its sequence signature almost matches that of the conventional kink-turn. However, the reverse and conventional kink-turns have opposite directions of bending. The reverse kink-turn lacks basically any tertiary interaction between its stems. We report unrestrained, explicit solvent molecular dynamics simulations of ribosomal and intron reverse kink-turns (54 simulations with 7.4 μs of data in total) with different variants (ff94, ff99, ff99bsc0, ff99χOL, and ff99bsc0χOL) of the Cornell et al. force field. We test several ion conditions and two water models. The simulations characterize the directional intrinsic flexibility of reverse kink-turns pertinent to their folded functional geometries. The reverse kink-turns are the most flexible RNA motifs studied so far by explicit solvent simulations which are capable at the present simulation time scale to spontaneously and reversibly sample a wide range of geometries from tightly kinked ones through flexible intermediates up to extended, unkinked structures. A possible biochemical role of the flexibility is discussed. Among the tested force fields, the latest χOL variant is essential to obtaining stable trajectories while all force field versions lacking the χ correction are prone to a swift degradation toward senseless ladder-like structures of stems, characterized by high-anti glycosidic torsions. The type of explicit water model affects the simulations considerably more than concentration and the type of ions.