Search

Kink-turns are highly bent internal loop motifs commonly found in the ribosome and other RNA complexes. They frequently act as binding sites for proteins and mediate tertiary interactions in larger RNA structures. Kink-turns have been a topic of intense research, but their elastic properties in the folded state are still poorly understood. Here we use extensive all-atom molecular dynamics simulations to parameterize a model of kink-turn in which the two flanking helical stems are represented by effective rigid bodies. Time series of the full set of six interhelical coordinates enable us to extract minimum energy shapes and harmonic stiffness constants for kink-turns from different RNA functional classes. The analysis suggests that kink-turns exhibit isotropic bending stiffness but are highly anisotropic with respect to lateral displacement of the stems. The most flexible lateral displacement mode is perpendicular to the plane of the static bend. These results may help understand the structural adaptation and mechanical signal transmission by kink-turns in complex natural and artificial RNA structures.

The double-crossover (DX) motif is a key building block of DNA nanostructures. It connects two double helices by two closely spaced Holliday junctions. Despite its prominent importance, the structure and elasticity of the DX motif is not fully understood. Here we employ extensive all-atom molecular dynamics (MD) simulations of an antiparallel DNA DX motif with two full turns between crossovers to infer its global structure and deformability. We quantitatively reproduce the experimentally known two-fold increase of bending stiffness upon incorporation of a DNA duplex into the DX motif, and find out that its stretching and twisting stiffness are only slightly influenced. To further describe the motif, we define four effective rigid bodies, each involving several base pairs flanking the Holliday junctions from the outside, and consider internal coordinates capturing relative displacement and orientation of the bodies. Time series of the coordinates from MD then yield the global structure of the DX motif and its stiffness in the harmonic approximation.

Mismatch repair is a highly conserved cellular pathway responsible for repairing mismatched dsDNA. Errors are detected by the MutS enzyme, which most likely senses altered mechanical property of damaged dsDNA rather than a specific molecular pattern. While the curved shape of dsDNA in crystallographic MutS/DNA structures suggests the role of DNA bending, the theoretical support is not fully convincing. Here, we present a computational study focused on a base-pair opening into the minor groove, a specific base-pair motion observed upon interaction with MutS. Propensities for the opening were evaluated in terms of two base-pair parameters: Opening and Shear. We tested all possible base pairs in anti/anti, anti/syn and syn/anti orientations and found clear discrimination between mismatches and canonical base-pairs only for the opening into the minor groove. Besides, the discrimination gap was also confirmed in hotspot and coldspot sequences, indicating that the opening could play a more significant role in the mismatch recognition than previously recognized. Our findings can be helpful for a better understanding of sequence-dependent mutability. Further, detailed structural characterization of mismatches can serve for designing anti-cancer drugs targeting mismatched base pairs.

Germline mutations in ETV6 gene cause inherited thrombocytopenia with leukemia predisposition. Here, we report on functional validation of ETV6 W380R mutation segregating with thrombocytopenia in a family where two family members also suffered from acute lymphoblastic leukemia (ALL) or essential thrombocythemia (ET). In-silico analysis predicted impaired DNA binding due to W380R mutation. Functional analysis showed that this mutation prevents the ETV6 protein from localizing into the cell nucleus and impairs the transcriptional repression activity of ETV6. Based on the germline ETV6 mutation, ET probably started with somatic JAK2 V617F mutation, whereas ALL could be caused by diverse mechanisms: high-hyperdiploidity; somatic deletion of exon 1 IKZF1 gene; or somatic mutations of other genes found by exome sequencing of the ALL sample taken at the diagnosis.

Bernard-Soulier syndrome (BSS) is a rare inherited disorder characterized by unusually large platelets, low platelet count, and prolonged bleeding time. BSS is usually inherited in an autosomal recessive (AR) mode of inheritance due to a deficiency of the GPIb-IX-V complex also known as the von Willebrand factor (VWF) receptor. We investigated a family with macrothrombocytopenia, a mild bleeding tendency, slightly lowered platelet aggregation tests, and suspected autosomal dominant (AD) inheritance. We have detected a heterozygous GP1BA likely pathogenic variant, causing monoallelic BSS. A germline GP1BA gene variant (NM_000173:c.98G > A:p.C33Y), segregating with the macrothrombocytopenia, was detected by whole-exome sequencing. In silico analysis of the protein structure of the novel GPIbα variant revealed a potential structural defect, which could impact proper protein folding and subsequent binding to VWF. Flow cytometry, immunoblot, and electron microscopy demonstrated further differences between p.C33Y GP1BA carriers and healthy controls. Here, we provide a detailed insight into its clinical presentation and phenotype. Moreover, the here described case first presents an mBSS patient with two previous ischemic strokes.

The primary goal was to determine the yield of next-generation sequencing (NGS) epilepsy gene panels used for epilepsy etiology diagnosing using a multidisciplinary approach and to demonstrate the importance of genotype-phenotype correlations. The secondary goal was to evaluate the application of precision medicine in selected patients.This single-center retrospective study included a total of 175 patients (95 males and 80 females) aged 0-19 years. They were examined between 2015 and 2020 using an NGS epilepsy gene panel (270 genes). A bioinformatic analysis was performed including copy number variation identification. Thorough genotype-phenotype correlation was performed.Out of 175 patients, described pathogenic variants or novel variants with clear pathogenic impact were identified in 30 patients (17.14%). Genotype-phenotype correlations and parental DNA analysis were performed, and genetic diagnosis was confirmed on the basis of the results in another 16 out of 175 patients (9.14%). The diagnostic yield of our study increased from 30 to 46 patients (by 53.33%) by the precise genotype-phenotype correlation.We emphasize a complex genotype-phenotype correlation and a multidisciplinary approach in evaluating the results of the NGS epilepsy gene panel, which enables the most accurate genetic diagnosis and correct interpretation of results.

Here we report a C-terminal RUNX1 mutation in a family with platelet disorder and predisposition to myeloid malignancies. We identified the mutation c.866delG:p.Gly289Aspfs*22 (NM_001754) (RUNX1 b-isoform NM_001001890; c.785delG:p.Gly262Aspfs*22) using exome sequencing of samples obtained from eight members of a single family. The mutation found in our pedigree is within exon eight and the transactivation domain of RUNX1. One of the affected individuals developed myelodysplastic syndrome (MDS), which progressed to acute myelogenous leukemia (AML). A search for the second hit which led to the development of MDS and later AML in this individual revealed the PHF6 gene variant (exon9:c.872G > A:p.G291E; NM_001015877), BCORL1 (exon3:c.1111A > C:p.T371P; NM_001184772) and BCOR gene variant (exon4:c.2076dupT:p.P693fs; NM_001123383), which appear to be very likely second hits participating in the progression to myeloid malignancy.

The genus Ceratostomella has a long history of taxonomic confusion. While species with evanescent asci have been transferred to the Microascales and Ophiostomatales, the taxonomic status of species with persistent asci has not been completely resolved. In previous studies using DNA sequence data, cultures and morphology, several Ceratostomella spp. were allocated in 13 genera in the Eurotiomycetes and Sordariomycetes. In our study, the systematics of the remaining Ceratostomella spp. with persistent asci is revisited with new collection data, cultures and phylogeny based on novel DNA sequences from six nuclear loci. Bayesian inference and Maximum Likelihood analyses support the monophyly of several wood-inhabiting species formerly classified in Ceratostomella and other unknown morphologically similar taxa and their division into four genera, i.e. Lentomitella, Spadicoides, Torrentispora and the newly described Calyptosphaeria. This robust clade represents the order Xenospadicoidales in the Sordariomycetidae. Comparative analysis of the ITS2 secondary structure revealed a genetic variation among Lentomitella isolates; 11 species were recognised, of which five are newly introduced and two are new combinations. Other taxonomic novelties include four new species and eight new combinations in Calyptosphaeria, Spadicoides, and Torrentispora. Molecular data suggest that Spadicoides is polyphyletic. The core of the genus is positioned in the Xenospadicoidales; Spadicoides s. str. is experimentally linked with sexual morphs for the first time. Based on DNA sequence data, the monotypic genera Xenospadicoides and Pseudodiplococcium are reduced to synonymy under Spadicoides, while Fusoidispora and Pseudoannulatascus are synonymised with Torrentispora. Members of the Xenospadicoidales inhabit decaying wood in terrestrial and freshwater environments and share a few morphological characters such as the absence of stromatic tissue, ascomata with a cylindrical or rostrate neck, similar anatomies of the ascomatal walls, thin-walled unitunicate asci with a non-amyloid apical annulus, disintegrating paraphyses, usually ellipsoidal to fusiform ascospores and holoblastic-denticulate or tretic conidiogenesis. Revised Ceratostomella spp. with persistent asci are listed and the taxonomic status of each species is re-evaluated based on revision of the holotype and other representative material, published details and available phylogenetic data. Key words: Ceratostomella, Conidiogenesis, Holoblastic-denticulate, Molecular systematics, New taxa, Phaeoisaria-like, Selenosporella-like, Tretic, Taxonomy, Xenospadicoidales, Taxonomic novelties: New genus: Calyptosphaeria Réblová & A.N. Mill, New species: Calyptosphaeria collapsa Réblová & A.N. Mill., C. tenebrosa Réblová & A.N. Mill., Lentomitella magna Réblová, L. obscura Réblová, L. striatella Réblová, L. sulcata Réblová, L. tenuirostris Réblová, Torrentispora calembola Réblová & A.N. Mill., T. novae-zelandiae Réblová & A.N. Mill, New combinations: Calyptosphaeria subdenudata (Peck) Réblová & A.N. Mill., C. tropica (Huhndorf et al.) Réblová & A.N. Mill., Lentomitella conoidea (Feltg.) Réblová, L. investita (Schw.) Réblová, Spadicoides fuscolutea (Rehm) Réblová, S. hyalostoma (Munk) Réblová, Spadicoides iberica (Hern.-Restr. et al.) Réblová & A.N. Mill., Torrentispora aquatica (Vijaykr. et al.) Réblová & A.N. Mill., T. biatriispora (K.D. Hyde) Réblová & A.N. Mill., T. dubia (Sacc.) Réblová & A.N. Mill

In recent decades, it has become clear that most of human proteins are glycosylated and that protein glycosylation plays an important role in health and diseases. At present, simple, fast and inexpensive methods are sought for clinical applications and particularly for improved diagnostics of various diseases, including cancer. We propose a label- and reagent-free electrochemical method based on chronopotentiometric stripping (CPS) analysis and a hanging mercury drop electrode for the detection of interaction of sialylated protein biomarker a prostate specific antigen (PSA) with two important lectins: Sambucus nigra agglutinin (SNA) and Maackia amurensis agglutinin (MAA). Incubation of PSA-modified electrode with specific SNA lectin resulted in an increase of CPS peak H of the complex as compared to this peak of individual PSA. By adjusting polarization current and temperature, PSA-MAA interaction can be either eliminated or distinguished from the more abundant PSA-SNA complex. CPS data were in a good agreement with the data obtained by complementary methods, namely surface plasmon resonance and fluorescent lectin microarray. It can be anticipated that CPS will find application in glycomics and proteomics.

C-loop is an internal loop motif found in the ribosome and used in artificial nanostructures. While its geometry has been partially characterized, its mechanical properties remain elusive. Here we propose a method to evaluate global shape and stiffness of an internal loop. The loop is flanked by short A-RNA helices modeled as rigid bodies. Their relative rotation and displacement are fully described by six interhelical coordinates. The deformation energy of the loop is assumed to be a general quadratic function of the interhelical coordinates. The model parameters for isolated C-loops are inferred from unrestrained all-atom molecular dynamics simulations. C-loops exhibit high twist as reported earlier, but also a bend and a lateral displacement of the flanking helices. Their bending stiffness and lateral displacement stiffness are nearly isotropic and similar to the control A-RNA duplexes. Nevertheless, we found systematic variations with the C-loop position in the ribosome and the organism of origin. The results characterize global properties of C-loops in the full six-dimensional interhelical space and enable one to choose an optimally stiff C-loop for use in a nanostructure. Our approach can be readily applied to other internal loops and extended to more complex structural motifs.

Splicing-affecting mutations can disrupt gene function by altering the transcript assembly. To ascertain splicing dysregulation principles, we modified a minigene assay for the parallel high-throughput evaluation of different mutations by next-generation sequencing. In our model system, all exonic and six intronic positions of the SMN1 gene's exon 7 were mutated to all possible nucleotide variants, which amounted to 180 unique single-nucleotide mutants and 470 double mutants. The mutations resulted in a wide range of splicing aberrations. Exonic splicing-affecting mutations resulted either in substantial exon skipping, supposedly driven by predicted exonic splicing silencer or cryptic donor splice site (5'ss) and de novo 5'ss strengthening and use. On the other hand, a single disruption of exonic splicing enhancer was not sufficient to cause major exon skipping, suggesting these elements can be substituted during exon recognition. While disrupting the acceptor splice site led only to exon skipping, some 5'ss mutations potentiated the use of three different cryptic 5'ss. Generally, single mutations supporting cryptic 5'ss use displayed better pre-mRNA/U1 snRNA duplex stability and increased splicing regulatory element strength across the original 5'ss. Analyzing double mutants supported the predominating splicing regulatory elements' effect, but U1 snRNA binding could contribute to the global balance of splicing isoforms. Based on these findings, we suggest that creating a new splicing enhancer across the mutated 5'ss can be one of the main factors driving cryptic 5'ss use.

The molecular genetics of well-characterized inherited diseases, such as phenylketonuria (PKU) and hyperphenylalaninemia (HPA) predominantly caused by mutations in the phenylalanine hydroxylase (PAH) gene, is often complicated by the identification of many novel variants, often with no obvious impact on the associated disorder. To date, more than 1100 PAH variants have been identified of which a substantial portion have unknown clinical significance. In this work, we study the functionality of seven yet uncharacterized PAH missense variants p.Asn167Tyr, p.Thr200Asn, p.Asp229Gly, p.Gly239Ala, p.Phe263Ser, p.Ala342Pro, and p.Ile406Met first identified in the Czech PKU/HPA patients. From all tested variants, three of them, namely p.Asn167Tyr, p.Thr200Asn, and p.Ile406Met, exerted residual enzymatic activity in vitro similar to wild type (WT) PAH, however, when expressed in HepG2 cells, their protein level reached a maximum of 72.1% &plusmn, 4.9%, 11.2% &plusmn, 4.2%, and 36.6% &plusmn, 7.3% compared to WT PAH, respectively. Remaining variants were null with no enzyme activity and decreased protein levels in HepG2 cells. The chaperone-like effect of applied BH4 precursor increased protein level significantly for p.Asn167Tyr, p.Asp229Gly, p.Ala342Pro, and p.Ile406Met. Taken together, our results of functional characterization in combination with in silico prediction suggest that while p.Asn167Tyr, p.Thr200Asn, and p.Ile406Met PAH variants have a mild impact on the protein, p.Asp229Gly, p.Gly239Ala, p.Phe263Ser, and p.Ala342Pro severely affect protein structure and function.

Different genes related to alveolar stability have been associated with familial interstitial pneumonia (FIP). Here, we report a novel, rare SFTPA1 variant in a family with idiopathic interstitial pneumonia (IIP). We performed whole-exome sequencing on germline DNA samples from four members of one family; three of them showed signs of pulmonary fibrosis (idiopathic interstitial pneumonia) with autosomal-dominant inheritance. A heterozygous single nucleotide variant c.532 G > A in the SFTPA1 gene has been identified. This variant encodes the substitution p.(Val178Met), localized within the carbohydrate recognition domain of surfactant protein A and segregates with the genes causing idiopathic interstitial pneumonia. This rare variant has not been previously reported. We also analyzed the detected sequence variant in the protein structure in silico. The replacement of valine by the larger methionine inside the protein may cause a disruption in the protein structure. The c.532 G > A variant was further validated using Sanger sequencing of the amplicons, confirming the diagnosis in all symptomatic family members. Moreover, this variant was also found by Sanger sequencing in one other symptomatic family member and one young asymptomatic family member. The autosomal-dominant inheritance, the family history of IIP, and the evidence of a mutation occurring in part of the SFTPA1 gene all suggest a novel variant that causes FIP., Lung disease: A different type of interstitial pneumonia Researchers have identified a novel mutation that causes inheritable and ultimately deadly scarring of the tiny air sacs in the lungs, the alveoli. In familial interstitial pneumonia (FIP), alveoli are scarred and stiffened by inflammation, not by bacterial infection as the word ‘pneumonia’ usually indicates. Michael Doubek at University Hospital and Central European Institute of Technology, Brno, Czech Republic and co-workers investigated a suspected case of FIP following the early death of one family member. Sequencing the genomes of other family members revealed that they shared a mutation in a protein that keeps alveoli moist, aiding oxygen absorption. Computer analysis showed that the mutation probably changed the protein’s shape, preventing it from functioning. Identifying mutations that cause FIP will help provide proactive treatment for family members who are at risk but not yet showing symptoms.

The low density lipoprotein receptor (LDLR) is a transmembrane protein that plays a key role in cholesterol metabolism. It contains 860 amino acids including a 21 amino acid long signal sequence, which directs the protein into the endoplasmic reticulum. Mutations in the LDLR gene lead to cholesterol accumulation in the plasma and results in familial hypercholesterolemia (FH). Knowledge of the impact of a mutation on the LDLR protein structure and function is very important for the diagnosis and management of FH. Unfortunately, for a large proportion of mutations this information is still missing. In this study, we focused on the LDLR signal sequence and carried out functional and in silico analyses of two sequence changes, p.(Gly20Arg) and p.(Leu15Pro), localized in this part of the LDLR. Our results revealed that the p.(Gly20Arg) change, previously described as disease causing, has no detrimental effect on protein expression or LDL particle binding. In silico analysis supports this observation, showing that both the wt and p.(Gly20Arg) signal sequences adopt an expected α-helix structure. In contrast, the mutation p.(Leu15Pro) is not associated with functional protein expression and exhibits a structure with disrupted a α-helical arrangement in the signal sequence, which most likely affects protein folding in the endoplasmic reticulum.

Mutations in the GP1BA gene have been associated with platelet-type von Willebrand disease and Bernard-Soulier syndrome. Here, we report a novel GP1BA mutation in a family with autosomal dominant macrothrombocytopenia and mild bleeding. We performed analyses of seven family members. Using whole-exome sequencing of germline DNA samples, we identified a heterozygous single-nucleotide change in GP1BA (exone2:c.176T>G), encoding a p.Leu59Arg substitution in the N-terminal domain, segregating with macrothrombocytopenia. This variant has not been previously reported. We also analysed the structure of the detected sequence variant in silico. In particular, we used the crystal structure of the human platelet receptor GP Ibα N-terminal domain. Replacement of aliphatic amino-acid Leu 59 with charged, polar and larger arginine probably disrupts the protein structure. An autosomal dominant mode of inheritance, a family history of mild bleeding episodes, aggregation pattern in affected individuals together with evidence of mutation occurring in part of the GP1BA gene encoding the leucine-rich repeat region suggest a novel variant causing monoallelic Bernard-Soulier syndrome.

Activation-induced cytidine deaminase (AID) is a mutator enzyme essential for somatic hypermutation (SHM) and class switch recombination (CSR) during effective adaptive immune responses. Its aberrant expression and activity have been detected in lymphomas, leukemias, and solid tumors. In chronic lymphocytic leukemia (CLL) increased expression of alternatively spliced AID variants has been documented. We used real-time RT-PCR to quantify the expression of AID and its alternatively spliced transcripts (AIDΔE4a, AIDΔE4, AIDivs3, and AIDΔE3E4) in 149 CLL patients and correlated this expression to prognostic markers including recurrent chromosomal aberrations, the presence of complex karyotype, mutation status of the immunoglobulin heavy chain variable gene, and recurrent mutations. We report a previously unappreciated association between higher AID transcript levels and trisomy of chromosome 12. Functional analysis of AID splice variants revealed loss of their activity with respect to SHM, CSR, and induction of double-strand DNA breaks. In silico modeling provided insight into the molecular interactions and structural dynamics of wild-type AID and a shortened AID variant closely resembling AIDΔE4, confirming its loss-of-function phenotype.

Adenosine to inosine (A&ndash, I) editing is the most common modification of double-stranded RNA (dsRNA). This change is mediated by adenosine deaminases acting on RNA (ADARs) enzymes with a preference of U&gt, A&gt, C&gt, G for 5&prime, neighbor and G&gt, C=A&gt, U or G&gt, U=A for 3&prime, neighbor. A&ndash, I editing occurs most frequently in the non-coding regions containing repetitive elements such as ALUs. It leads to disruption of RNA duplex structure, which prevents induction of innate immune response. We employed standard and biased molecular dynamics (MD) simulations to analyze the behavior of RNA duplexes with single and tandem inosine&ndash, uracil (I&ndash, U) base pairs in different sequence context. Our analysis showed that the I&ndash, U pairs induce changes in base pair and base pair step parameters and have different dynamics when compared with standard canonical base pairs. In particular, the first I&ndash, U pair from tandem I&ndash, U/I&ndash, U systems exhibited increased dynamics depending on its neighboring 5&prime, base. We discovered that UII sequence, which is frequently edited, has lower flexibility compared with other sequences (AII, GII, CII), hence it only modestly disrupts dsRNA. This might indicate that the UAA motifs in ALUs do not have to be sufficiently effective in preventing immune signaling.

Thirteen morphologically similar strains of barbatosphaeria- and tectonidula-like fungi were studied based on the comparison of cultural and morphological features of sexual and asexual morphs and phylogenetic analyses of five nuclear loci, i.e. internal transcribed spacer rDNA operon (ITS), large and small subunit nuclear ribosomal DNA, β-tubulin, and second largest subunit of RNA polymerase II. Phylogenetic results were supported by in-depth comparative analyses of common core secondary structure of ITS1 and ITS2 in all strains and the identification of non-conserved, co-evolving nucleotides that maintain base pairing in the RNA transcript. Barbatosphaeria is defined as a well-supported monophyletic clade comprising several lineages and is placed in the Sordariomycetes incertae sedis. The genus is expanded to encompass nine species with both septate and non-septate ascospores in clavate, stipitate asci with a non-amyloid apical annulus and non-stromatic ascomata with a long decumbent neck and carbonised wall often covered by pubescence. The asexual morphs are dematiaceous hyphomycetes with holoblastic conidiogenesis belonging to Ramichloridium and Sporothrix types. The morphologically similar Tectonidula, represented by the type species T. hippocrepida, grouped with members of Barbatosphaeria and is transferred to that genus. Four new species are introduced and three new combinations in Barbatosphaeria are proposed. A dichotomous key to species accepted in the genus is provided.

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.

Phenylketonuria (PKU) and hyperphenylalaninemia (HPA) are a group of genetic disorders predominantly caused by mutations in the phenylalanine hydroxylase (PAH) gene. To date, more than 950 variants have been identified, however the pathogenic mechanism of many variants remains unknown. In this study, in silico prediction and in vitro prokaryotic and eukaryotic expression systems were used to functionally characterize five PAH missense variants (p.F233I, p.R270I, p.F331S, p.S350Y, and p.L358F) previously identified in Slovak and Czech patients. p.F233I, p.R270I, and p.S350Y were classified as deleterious mutations since they showed no specific activity in functional assay and no response to chaperone co-expression. Protein levels of these PAH variants were very low when expressed in HepG2 cells, and only p.S350Y responded to BH4 precursor overload by significant increase in PAH monomer, probably due to reduced rate of protein degradation as the result of proper protein folding. Variants p.F331S and p.L358F exerted residual enzymatic activity in vitro. While the first can be classified as probably pathogenic due to its very low protein levels in HepG2 cells, the latter is considered to be mild mutation with protein levels of approximately 17.85% compared to wt PAH. Our findings contribute to better understanding of structure and function of PAH mutated enzymes and optimal treatment of PKU patients carrying these mutations using BH4 supplementation.

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.

In a case study of fungi of the class Sordariomycetes, we evaluated the effect of multiple sequence alignment (MSA) on the reliability of the phylogenetic trees, topology and confidence of major phylogenetic clades. We compared two main approaches for constructing MSA based on (1) the knowledge of the secondary (2D) structure of ribosomal RNA (rRNA) genes, and (2) automatic construction of MSA by four alignment programs characterized by different algorithms and evaluation methods, CLUSTAL, MAFFT, MUSCLE, and SAM. In the primary fungal sequences of the two functional rRNA genes, the nuclear small and large ribosomal subunits (18 S and 28 S), we identified four and six, respectively, highly variable regions, which correspond mainly to hairpin loops in the 2D structure. These loops are often positioned in expansion segments, which are missing or are not completely developed in the Archaeal and Eubacterial kingdoms. Proper sorting of these sites was a key for constructing an accurate MSA. We utilized DNA sequences from 28 S as an example for one-gene analysis. Five different MSAs were created and analyzed with maximum parsimony and maximum likelihood methods. The phylogenies inferred from the alignments improved with 2D structure with identified homologous segments, and those constructed using the MAFFT alignment program, with all highly variable regions included, provided the most reliable phylograms with higher bootstrap support for the majority of clades. We illustrate and provide examples demonstrating that re-evaluating ambiguous positions in the consensus sequences using 2D structure and covariance is a promising means in order to improve the quality and reliability of sequence alignments.

The L1 stalk is a key mobile element of the large ribosomal subunit which interacts with tRNA during translocation. Here, we investigate the structure and mechanical properties of the rRNA H76/H75/H79 three-way junction at the base of the L1 stalk from four different prokaryotic organisms. We propose a coarse-grained elastic model and parameterize it using large-scale atomistic molecular dynamics simulations. Global properties of the junction are well described by a model in which the H76 helix is represented by a straight, isotropically flexible elastic rod, while the junction core is represented by an isotropically flexible spherical hinge. Both the core and the helix contribute substantially to the overall H76 bending fluctuations. The presence of wobble pairs in H76 does not induce any increased flexibility or anisotropy to the helix. The half-closed conformation of the L1 stalk seems to be accessible by thermal fluctuations of the junction itself, without any long-range allosteric effects. Bending fluctuations of H76 with a bulge introduced in it suggest a rationale for the precise position of the bulge in eukaryotes. Our elastic model can be generalized to other RNA junctions found in biological systems or in nanotechnology.

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.

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.

Kink-turns (K-turns) are recurrent elbow-like RNA motifs that participate in protein-assisted RNA folding and contribute to RNA dynamics. We carried out a set of molecular dynamics (MD) simulations using parm99 and parmbsc0 force fields to investigate structural dynamics of the box C/D RNA and its complexes with two proteins: native archaeal L7ae protein and human 15.5 kDa protein, originally bound to very similar structure of U4 snRNA. The box C/D RNA forms K-turn with A-minor 0 tertiary interaction between its canonical (C) and noncanonical (NC) stems. The local K-turn architecture is thus different from the previously studied ribosomal K-turns 38 and 42 having A-minor I interaction. The simulations reveal visible structural dynamics of this tertiary interaction involving altogether six substates which substantially contribute to the elbow-like flexibility of the K-turn. The interaction can even temporarily shift to the A-minor I type pattern; however, this is associated with distortion of the G/A base pair in the NC-stem of the K-turn. The simulations show reduction of the K-turn flexibility upon protein binding. The protein interacts with the apex of the K-turn and with the NC-stem. The protein-RNA interface includes long-residency hydration sites. We have also found long-residency hydration sites and major ion-binding sites associated with the K-turn itself. The overall topology of the K-turn remains stable in all simulations. However, in simulations of free K-turn, we observed instability of the key C16(O2')-A7(N1) H-bond, which is a signature interaction of K-turns and which was visibly more stable in simulations of K-turns possessing A-minor I interaction. It may reflect either some imbalance of the force field or it may be a correct indication of early stages of unfolding since this K-turn requires protein binding for its stabilization. Interestingly, the 16(O2')-7(N1) H- bond is usually not fully lost since it is replaced by a water bridge with a tightly bound water, which is adenine-specific similarly as the original interaction. The 16(O2')-7(N1) H-bond is stabilized by protein binding. The stabilizing effect is more visible with the human 15.5 kDa protein, which is attributed to valine to arginine substitution in the binding site. The behavior of the A-minor interaction is force-field-dependent because the parmbsc0 force field attenuates the A-minor fluctuations compared to parm99 simulations. Behavior of other regions of the box C/D RNA is not sensitive to the force field choice. Simulation with net-neutralizing Na(+) and 0.2 M excess salt conditions appear in all aspects equivalent. The simulations show loss of a hairpin tetraloop, which is not part of the K-turn. This was attributed to force field limitations.

We present extensive explicit solvent molecular dynamics analysis of three RNA three-way junctions (3WJs) from the large ribosomal subunit: the 3WJ formed by Helices 90-92 (H90-H92) of 23S rRNA; the 3WJ formed by H42-H44 organizing the GTPase associated center (GAC) of 23S rRNA; and the 3WJ of 5S rRNA. H92 near the peptidyl transferase center binds the 3'-CCA end of amino-acylated tRNA. The GAC binds protein factors and stimulates GTP hydrolysis driving protein synthesis. The 5S rRNA binds the central protuberance and A-site finger (ASF) involved in bridges with the 30S subunit. The simulations reveal that all three 3WJs possess significant anisotropic hinge-like flexibility between their stacked stems and dynamics within the compact regions of their adjacent stems. The A-site 3WJ dynamics may facilitate accommodation of tRNA, while the 5S 3WJ flexibility appears to be essential for coordinated movements of ASF and 5S rRNA. The GAC 3WJ may support large-scale dynamics of the L7/L12-stalk region. The simulations reveal that H42-H44 rRNA segments are not fully relaxed and in the X-ray structures they are bent towards the large subunit. The bending may be related to L10 binding and is distributed between the 3WJ and the H42-H97 contact.

The hairpin ribozyme is a prominent member of the group of small catalytic RNAs (RNA enzymes or ribozymes) because it does not require metal ions to achieve catalysis. Biochemical and structural data have implicated guanine 8 (G8) and adenine 38 (A38) as catalytic participants in cleavage and ligation catalyzed by the hairpin ribozyme, yet their exact role in catalysis remains disputed. To gain insight into dynamics in the active site of a minimal self-cleaving hairpin ribozyme, we have performed extensive classical, explicit-solvent molecular dynamics (MD) simulations on time scales of 50-150 ns. Starting from the available X-ray crystal structures, we investigated the structural impact of the protonation states of G8 and A38, and the inactivating A-1(2'-methoxy) substitution employed in crystallography. Our simulations reveal that a canonical G8 agrees well with the crystal structures while a deprotonated G8 profoundly distorts the active site. Thus MD simulations do not support a straightforward participation of the deprotonated G8 in catalysis. By comparison, the G8 enol tautomer is structurally well tolerated, causing only local rearrangements in the active site. Furthermore, a protonated A38H(+) is more consistent with the crystallography data than a canonical A38. The simulations thus support the notion that A38H(+) is the dominant form in the crystals, grown at pH 6. In most simulations, the canonical A38 departs from the scissile phosphate and substantially perturbs the structures of the active site and S-turn. Yet, we occasionally also observe formation of a stable A-1(2'-OH)...A38(N1) hydrogen bond, which documents the ability of the ribozyme to form this hydrogen bond, consistent with a potential role of A38 as general base catalyst. The presence of this hydrogen bond is, however, incompatible with the expected in-line attack angle necessary for self-cleavage, requiring a rapid transition of the deprotonated 2'-oxyanion to a position more favorable for in-line attack after proton transfer from A-1(2'-OH) to A38(N1). The simulations revealed a potential force field artifact, occasional but irreversible formation of "ladder-like", underwound A-RNA structure in one of the external helices. Although it does not affect the catalytic center of the hairpin ribozyme, further studies are under way to better assess possible influence of such force field behavior on long RNA simulations.

Helix 38 (H38) of the large ribosomal subunit, with a length of 110 Å, reaches the small subunit through intersubunit bridge B1a. Previous cryo-EM studies revealed that the tip of H38 moves by more than 10 Å from the non-ratcheted to the ratcheted state of the ribosome while mutational studies implicated a key role of flexible H38 in attenuation of translocation and in dynamical signaling between ribosomal functional centers. We investigate a region including the elbow-shaped kink-turn (Kt-38) in the Haloarcula marismortui archaeal ribosome, and equivalently positioned elbows in three eubacterial species, located at the H38 base. We performed explicit solvent molecular dynamics simulations on the H38 elbows in all four species. They are formed by at first sight unrelated sequences resulting in diverse base interactions but built with the same overall topology, as shown by X-ray crystallography. The elbows display similar fluctuations and intrinsic flexibilities in simulations indicating that the eubacterial H38 elbows are structural and dynamical analogs of archaeal Kt-38. We suggest that this structural element plays a pivotal role in the large motions of H38 and may act as fulcrum for the abovementioned tip motion. The directional flexibility inferred from simulations correlates well with the cryo-EM results.

Dimerization of HIV-1 genomic RNA is initiated by kissing loop interactions at the Dimerization Initiation Site (DIS). Dynamics of purines that flank the 5' ends of the loop-loop helix in HIV-1 DIS kissing complex were explored using explicit solvent molecular dynamics (MD) simulations with the CHARMM force field. Multiple MD simulations (200 ns in total) of X-ray structures for HIV-1 DIS Subtypes A, B, and F revealed conformational variability of flanking purines. In particular, the flanking purines, which in the starting X-ray structures are bulged-out and stack in pairs, formed a consecutive stack of four bulged-out adenines at the beginning of several simulations. This conformation is seen in the crystal structure of DIS Subtype F with no interference from crystal packing, and was frequently reported in our preceding MD studies performed with the AMBER force field. However, as CHARMM simulations progressed, the four continuously stacked adenines showed conformational transitions from the bulged-out into the bulged-in geometries. Although such an arrangement has not been seen in any X-ray structure, it has been suggested by a recent NMR investigation. In CHARMM simulations, in the longer time scale, the flanking purines display the tendency to move to bulged-in conformations. This is in contrast with the AMBER simulations, which indicate a modest prevalence for bulged-out flanking base positions in line with the X-ray data. The simulations also suggest that the intermolecular stacking between purines from the opposite hairpins can additionally stabilize the kissing complex.

The Calosphaeriales is revisited with new collection data, living cultures, morphological studies of ascoma centrum, secondary structures of the internal transcribed spacer (ITS) rDNA and phylogeny based on novel DNA sequences of five nuclear ribosomal and protein-coding loci. Morphological features, molecular evidence and information from predicted RNA secondary structures of ITS converged upon robust phylogenies of the Calosphaeriales and Togniniales. The current concept of the Calosphaeriales includes the Calosphaeriaceae and Pleurostomataceae encompassing five monophyletic genera, Calosphaeria, Flabellascus gen. nov., Jattaea, Pleurostoma and Togniniella, strongly supported by Bayesian and Maximum Likelihood methods. The structural elements of ITS1 form characteristic patterns that are phylogenetically conserved, corroborate observations based on morphology and have a high predictive value at the generic level. Three major clades containing 44 species of Phaeoacremonium were recovered in the closely related Togniniales based on ITS, actin and β-tubulin sequences. They are newly characterized by sexual and RNA structural characters and ecology. This approach is a first step towards understanding of the molecular systematics of Phaeoacremonium and possibly its new classification. In the Calosphaeriales, Jattaea aphanospora sp. nov. and J. ribicola sp. nov. are introduced, Calosphaeria taediosa is combined in Jattaea and epitypified. The sexual morph of Phaeoacremonium cinereum was encountered for the first time on decaying wood and obtained in vitro. In order to achieve a single nomenclature, the genera of asexual morphs linked with the Calosphaeriales are transferred to synonymy of their sexual morphs following the principle of priority, i.e. Calosphaeriophora to Calosphaeria, Phaeocrella to Togniniella and Pleurostomophora to Pleurostoma. Three new combinations are proposed, i.e. Pleurostoma ochraceum comb. nov., P. repens comb. nov. and P. richardsiae comb. nov. The morphology-based key is provided to facilitate identification of genera accepted in the Calosphaeriales.

Autosomal recessive congenital ichthyosis (ARCI) represents a heterogeneous group of disorders of epidermal cornification. Nine genes have been identified to be causative of ARCI, including TGM1 1,2 , ABCA12 3 , NIPAL4 4 , CYP4F22 5 , ALOX12B, ALOXE3 6 , PNPLA1 7 , LIPN 8 , and CERS3 9 . ARCI is rare, with a reported prevalence 1 in 200,000 in European and northern American populations 10 . We started DNA analysis of ARCI in 2012 when PCR-direct sequencing of the TGM1, ALOX12B, ALOXE3, NIPAL4, and CYP4F22 genes was introduced and patients were analysed step by step for mutations in these genes.

Hyperphenylalaninemia (HPA) is one of the most common metabolic disorders. HPA, which is transmitted by an autosomal recessive mode of inheritance, is caused by mutations of the phenylalanine hydroxylase gene. Most mutations are missense and lead to reduced protein stability and/or impaired catalytic function. The impact of such mutations varies, ranging from classical phenylketonuria (PKU), mild PKU, to non-PKU HPA phenotypes. Despite the fact that HPA is a monogenic disease, clinical data show that one PKU genotype can be associated with more in vivo phenotypes, which indicates the role of other (still unknown) factors. To better understand the phenotype–genotype relationships, we analyzed computationally the impact of missense mutations in homozygotes stored in the BIOPKU database. A total of 34 selected homozygous genotypes was divided into two main groups according to their phenotypes: (A) genotypes leading to non-PKU HPA or combined phenotype non-PKU HPA/mild PKU and (B) genotypes leading to classical PKU, mild PKU or combined phenotype mild PKU/classical PKU. Combining in silico analysis and molecular dynamics simulations (in total 3000 ns) we described the structural impact of the mutations, which allowed us to separate 32 out of 34 mutations between groups A and B. Testing the simulation conditions revealed that the outcome of mutant simulations can be modulated by the ionic strength.We also employed programs SNPs3D, Polyphen-2, and SIFT but based on the predictions performed we were not able to discriminate mutations with mild and severe PKU phenotypes.

Ribozymes are catalytically competent examples of highly structured noncoding RNAs, which are ubiquitous in the processing and regulation of genetic information. Combining explicit-solvent molecular dynamics simulation and single molecule fluorescence spectroscopy approaches, we find that a ribozyme from a subviral plant pathogen exhibits a coupled hydrogen bonding network that communicates dynamic structural rearrangements throughout the catalytic core in response to site-specific chemical modification. Trapped long-residency water molecules are critical for this network and only occasionally exchange with bulk solvent as they pass through a breathing interdomain base stack. These highly structured water molecules line up in a string that may potentially also be involved in specific base catalysis. Our observations suggest important, still underappreciated roles for specifically bound water molecules in the structural dynamics and function of noncoding RNAs.

Molecular dynamics (MD) simulations were employed to investigate the structure, dynamics, and local base-pair step deformability of the free 16S ribosomal helix 44 from Thermus thermophilus and of a canonical A-RNA double helix. While helix 44 is bent in the crystal structure of the small ribosomal subunit, the simulated helix 44 is intrinsically straight. It shows, however, substantial instantaneous bends that are isotropic. The spontaneous motions seen in simulations achieve large degrees of bending seen in the X-ray structure and would be entirely sufficient to allow the dynamics of the upper part of helix 44 evidenced by cryo-electron microscopic studies. Analysis of local base-pair step deformability reveals a patch of flexible steps in the upper part of helix 44 and in the area proximal to the bulge bases, suggesting that the upper part of helix 44 has enhanced flexibility. The simulations identify two conformational substates of the second bulge area (bottom part of the helix) with distinct base pairing. In agreement with nuclear magnetic resonance (NMR) and X-ray studies, a flipped out conformational substate of conserved 1492A is seen in the first bulge area. Molecular dynamics (MD) simulations reveal a number of reversible alpha-gamma backbone flips that correspond to transitions between two known A-RNA backbone families. The flipped substates do not cumulate along the trajectory and lead to a modest transient reduction of helical twist with no significant influence on the overall geometry of the duplexes. Despite their considerable flexibility, the simulated structures are very stable with no indication of substantial force field inaccuracies.

Ribosomal RNA K-turn motifs are asymmetric internal loops characterized by a sharp bend in the phosphodiester backbone resulting in "V" shaped structures, recurrently observed in ribosomes and showing a high degree of sequence conservation. We have carried out extended explicit solvent molecular dynamics simulations of selected K-turns, in order to investigate their intrinsic structural and dynamical properties. The simulations reveal an unprecedented dynamical flexibility of the K-turns around their X-ray geometries. The K-turns sample, on the nanosecond timescale, different conformational substates. The overall behavior of the simulations suggests that the sampled geometries are essentially isoenergetic and separated by minimal energy barriers. The nanosecond dynamics of isolated K-turns can be qualitatively considered as motion of two rigid helix stems controlled by a very flexible internal loop which then leads to substantial hinge-like motions between the two stems. This internal dynamics of K-turns is strikingly different for example from the bacterial 5S rRNA Loop E motif or BWYV frameshifting pseudoknot which appear to be rigid in the same type of simulations. Bistability and flexibility of K-turns was also suggested by several recent biochemical studies. Although the results of MD simulations should be considered as a qualitative picture of the K-turn dynamics due to force field and sampling limitations, the main advantage of the MD technique is its ability to investigate the region close to K-turn ribosomal-like geometries. This part of the conformational space is not well characterized by the solution experiments due to large-scale conformational changes seen in the experiments. We suggest that K-turns are well suited to act as flexible structural elements of ribosomal RNA. They can for example be involved in mediation of large-scale motions or they can allow a smooth assembling of the other parts of the ribosome.