The pseudotorsional space of RNA.

The characterization of the conformational landscape of the RNA backbone is rather complex due to the ability of RNA to assume a large variety of conformations. These backbone conformations can be depicted by pseudotorsional angles linking RNA backbone atoms, from which Ramachandran-like plots can be built. We explore here different definitions of these pseudotorsional angles, finding that the most accurate ones are the traditional η (eta) and θ (theta) angles, which represent the relative position of RNA backbone atoms P and C4'. We explore the distribution of η - θ in known experimental structures, comparing the pseudotorsional space generated with structures determined exclusively by one experimental technique. We found that the complete picture only appears when combining data from different sources. The maps provide a quite comprehensive representation of the RNA accessible space, which can be used in RNA-structural predictions. Finally, our results highlight that protein interactions lead to significant changes in the population of the η - θ space, pointing toward the role of induced-fit mechanisms in protein-RNA recognition.

Histone H3 serine-57 is a CHK1 substrate whose phosphorylation affects DNA repair.

Histone post-translational modifications promote a chromatin environment that controls transcription, DNA replication and repair, but surprisingly few phosphorylations have been documented. We report the discovery of histone H3 serine-57 phosphorylation (H3S57ph) and show that it is implicated in different DNA repair pathways from fungi to vertebrates. We identified CHK1 as a major human H3S57 kinase, and disrupting or constitutively mimicking H3S57ph had opposing effects on rate of recovery from replication stress, 53BP1 chromatin binding, and dependency on RAD52. In fission yeast, mutation of all H3 alleles to S57A abrogated DNA repair by both non-homologous end-joining and homologous recombination, while cells with phospho-mimicking S57D alleles were partly compromised for both repair pathways, presented aberrant Rad52 foci and were strongly sensitised to replication stress. Mechanistically, H3S57ph loosens DNA-histone contacts, increasing nucleosome mobility, and interacts with H3K56. Our results suggest that dynamic phosphorylation of H3S57 is required for DNA repair and recovery from replication stress, opening avenues for investigating the role of this modification in other DNA-related processes.

A Fijivirus Major Viroplasm Protein Shows RNA-Stimulated ATPase Activity by Adopting Pentameric and Hexameric Assemblies of Dimers.

Fijiviruses replicate and package their genomes within viroplasms in a process involving RNA-RNA and RNA-protein interactions. Here, we demonstrate that the 24 C-terminal residues (C-arm) of the P9-1 major viroplasm protein of the mal de Río Cuarto virus (MRCV) are required for its multimerization and the formation of viroplasm-like structures. Using an integrative structural approach, the C-arm was found to be dispensable for P9-1 dimer assembly but essential for the formation of pentamers and hexamers of dimers (decamers and dodecamers), which favored RNA binding. Although both P9-1 and P9-1ΔC-arm catalyzed ATP with similar activities, an RNA-stimulated ATPase activity was only detected in the full-length protein, indicating a C-arm-mediated interaction between the ATP catalytic site and the allosteric RNA binding sites in the (do)decameric assemblies. A stronger preference to bind phosphate moieties in the decamer was predicted, suggesting that the allosteric modulation of ATPase activity by RNA is favored in this structural conformation. Our work reveals the structural versatility of a fijivirus major viroplasm protein and provides clues to its mechanism of action. IMPORTANCE The mal de Río Cuarto virus (MRCV) causes an important maize disease in Argentina. MRCV replicates in several species of Gramineae plants and planthopper vectors. The viral factories, also called viroplasms, have been studied in detail in animal reovirids. This work reveals that a major viroplasm protein of MRCV forms previously unidentified structural arrangements and provides evidence that it may simultaneously adopt two distinct quaternary assemblies. Furthermore, our work uncovers an allosteric communication between the ATP and RNA binding sites that is favored in the multimeric arrangements. Our results contribute to the understanding of plant reovirids viroplasm structure and function and pave the way for the design of antiviral strategies for disease control.

MiOS, an integrated imaging and computational strategy to model gene folding with nucleosome resolution.

The linear sequence of DNA provides invaluable information about genes and their regulatory elements along chromosomes. However, to fully understand gene function and regulation, we need to dissect how genes physically fold in the three-dimensional nuclear space. Here we describe immuno-OligoSTORM, an imaging strategy that reveals the distribution of nucleosomes within specific genes in super-resolution, through the simultaneous visualization of DNA and histones. We combine immuno-OligoSTORM with restraint-based and coarse-grained modeling approaches to integrate super-resolution imaging data with Hi-C contact frequencies and deconvoluted micrococcal nuclease-sequencing information. The resulting method, called Modeling immuno-OligoSTORM, allows quantitative modeling of genes with nucleosome resolution and provides information about chromatin accessibility for regulatory factors, such as RNA polymerase II. With Modeling immuno-OligoSTORM, we explore intercellular variability, transcriptional-dependent gene conformation, and folding of housekeeping and pluripotency-related genes in human pluripotent and differentiated cells, thereby obtaining the highest degree of data integration achieved so far to our knowledge.

Molecular basis of Arginine and Lysine DNA sequence-dependent thermo-stability modulation.

We have used a variety of theoretical and experimental techniques to study the role of four basic amino acids-Arginine, Lysine, Ornithine and L-2,4-Diaminobutyric acid-on the structure, flexibility and sequence-dependent stability of DNA. We found that the presence of organic ions stabilizes the duplexes and significantly reduces the difference in stability between AT- and GC-rich duplexes with respect to the control conditions. This suggests that these amino acids, ingredients of the primordial soup during abiogenesis, could have helped to equalize the stability of AT- and GC-rich DNA oligomers, facilitating a general non-catalysed self-replication of DNA. Experiments and simulations demonstrate that organic ions have an effect that goes beyond the general electrostatic screening, involving specific interactions along the grooves of the double helix. We conclude that organic ions, largely ignored in the DNA world, should be reconsidered as crucial structural elements far from mimics of small inorganic cations.

The Impact of the HydroxyMethylCytosine epigenetic signature on DNA structure and function.

We present a comprehensive, experimental and theoretical study of the impact of 5-hydroxymethylation of DNA cytosine. Using molecular dynamics, biophysical experiments and NMR spectroscopy, we found that Ten-Eleven translocation (TET) dioxygenases generate an epigenetic variant with structural and physical properties similar to those of 5-methylcytosine. Experiments and simulations demonstrate that 5-methylcytosine (mC) and 5-hydroxymethylcytosine (hmC) generally lead to stiffer DNA than normal cytosine, with poorer circularization efficiencies and lower ability to form nucleosomes. In particular, we can rule out the hypothesis that hydroxymethylation reverts to unmodified cytosine physical properties, as hmC is even more rigid than mC. Thus, we do not expect dramatic changes in the chromatin structure induced by differences in physical properties between d(mCpG) and d(hmCpG). Conversely, our simulations suggest that methylated-DNA binding domains (MBDs), associated with repression activities, are sensitive to the substitution d(mCpG) ➔ d(hmCpG), while MBD3 which has a dual activation/repression activity is not sensitive to the d(mCpG) d(hmCpG) change. Overall, while gene activity changes due to cytosine methylation are the result of the combination of stiffness-related chromatin reorganization and MBD binding, those associated to 5-hydroxylation of methylcytosine could be explained by a change in the balance of repression/activation pathways related to differential MBD binding.

Impact of DNA methylation on 3D genome structure.

Determining the effect of DNA methylation on chromatin structure and function in higher organisms is challenging due to the extreme complexity of epigenetic regulation. We studied a simpler model system, budding yeast, that lacks DNA methylation machinery making it a perfect model system to study the intrinsic role of DNA methylation in chromatin structure and function. We expressed the murine DNA methyltransferases in Saccharomyces cerevisiae and analyzed the correlation between DNA methylation, nucleosome positioning, gene expression and 3D genome organization. Despite lacking the machinery for positioning and reading methylation marks, induced DNA methylation follows a conserved pattern with low methylation levels at the 5' end of the gene increasing gradually toward the 3' end, with concentration of methylated DNA in linkers and nucleosome free regions, and with actively expressed genes showing low and high levels of methylation at transcription start and terminating sites respectively, mimicking the patterns seen in mammals. We also see that DNA methylation increases chromatin condensation in peri-centromeric regions, decreases overall DNA flexibility, and favors the heterochromatin state. Taken together, these results demonstrate that methylation intrinsically modulates chromatin structure and function even in the absence of cellular machinery evolved to recognize and process the methylation signal.

Sequence-dependent structural properties of B-DNA: what have we learned in 40 years?

The structure of B-DNA, the physiological form of the DNA molecule, has been a central topic in biology, chemistry and physics. Far from uniform and rigid, the double helix was revealed as a flexible and structurally polymorphic molecule. Conformational changes that lead to local and global changes in the helix geometry are mediated by a complex choreography of base and backbone rearrangements affecting the ability of the B-DNA to recognize ligands and consequently on its functionality. In this sense, the knowledge obtained from the sequence-dependent structural properties of B-DNA has always been thought crucial to rationalize how ligands and, most notably, proteins recognize B-DNA and modulate its activity, i.e. the structural basis of gene regulation. Honouring the anniversary of the first high-resolution X-ray structure of a B-DNA molecule, in this contribution, we present the most important discoveries of the last 40 years on the sequence-dependent structural and dynamical properties of B-DNA, from the early beginnings to the current frontiers in the field.

Antimicrobial peptides in the seedling transcriptome of the tree legume Peltophorum dubium.

Antimicrobial peptides (AMPs) play an essential role in plant defense against invading pathogens. Due to their biological properties, these molecules have been considered useful for drug development, as novel agents in disease therapeutics, applicable to both agriculture and medicine. New technologies of massive sequencing open opportunities to discover novel AMP encoding genes in wild plant species. This work aimed to identify cysteine-rich AMPs from Peltophorum dubium, a legume tree from South America. We performed whole-transcriptome sequencing of P. dubium seedlings followed by de novo transcriptome assembly, uncovering 78 AMP transcripts classified into five families: hevein-like, lipid-transfer proteins (LTPs), alpha hairpinins, defensins, and snakin/GASA (Giberellic Acid Stimulated in Arabidopsis) peptides. No transcripts with similarity to cyclotide or thionin genes were identified. Genomic DNA analysis by PCR confirmed the presence of 18 genes encoding six putative defensins and 12 snakin/GASA peptides and allowed the characterization of their exon-intron structure. The present work demonstrates that AMP prediction from a wild species is possible using RNA sequencing and de novo transcriptome assembly, regarding a starting point for studies focused on AMP gene evolution and expression. Moreover, this study allowed the detection of strong AMP candidates for drug development and novel biotechnological products.

Molecular Determinants for Nitric Oxide Regulation of the Murine Cationic Amino Acid Transporter CAT-2A.

Cationic amino acid transporters (CATs) supply cells with essential and semiessential dibasic amino acids. Among them, l-arginine is the substrate for nitric oxide synthases (NOS) to produce nitric oxide (NO), a key signaling molecule and second messenger. In cardiac preparations, we showed that NO acutely and directly modulates transport activity by noncompetitively inhibiting these CATs. We hypothesize that this NO regulation occurs through modification of cysteine residues in CAT proteins. Homology modeling and a computational chemistry approach identified Cys347 as one of two putative targets for NO binding, of 15 Cys residues present in the low-affinity mouse CAT-2A (mCAT-2A). To test this prediction, mammalian cell lines overexpressing mCAT-2A were used for site-directed mutagenesis and uptake studies. When Cys347 was replaced with alanine (Cys347Ala), mCAT-2A became insensitive to inhibition by NO donors. In addition, the transport capacity of this variant decreased by >50% compared to that of the control, without affecting membrane expression levels or apparent affinities for the transported amino acids. Interestingly, replacing Cys347 with serine (Cys347Ser) restored uptake levels to those of the control while retaining NO insensitivity. Other Cys residues, when replaced with Ala, still produced a NO-sensitive CAT-2A. In cells co-expressing NOS and mCAT-2A, exposure to extracellular l-arginine inhibited the uptake activity of control mCAT-2A, via NO production, but not that of the Cys347Ser variant. Thus, the -SH moiety of Cys347 is largely responsible for mCAT-2A inhibition by NO. Because of the endogenous NO effect, this modulation is likely to be physiologically relevant and a potential intervention point for therapeutics.

A multi-modal coarse grained model of DNA flexibility mappable to the atomistic level.

We present a new coarse grained method for the simulation of duplex DNA. The algorithm uses a generalized multi-harmonic model that can represent any multi-normal distribution of helical parameters, thus avoiding caveats of current mesoscopic models for DNA simulation and representing a breakthrough in the field. The method has been parameterized from accurate parmbsc1 atomistic molecular dynamics simulations of all unique tetranucleotide sequences of DNA embedded in long duplexes and takes advantage of the correlation between helical states and backbone configurations to derive atomistic representations of DNA. The algorithm, which is implemented in a simple web interface and in a standalone package reproduces with high computational efficiency the structural landscape of long segments of DNA untreatable by atomistic molecular dynamics simulations.

The static and dynamic structural heterogeneities of B-DNA: extending Calladine-Dickerson rules.

We present a multi-laboratory effort to describe the structural and dynamical properties of duplex B-DNA under physiological conditions. By processing a large amount of atomistic molecular dynamics simulations, we determine the sequence-dependent structural properties of DNA as expressed in the equilibrium distribution of its stochastic dynamics. Our analysis includes a study of first and second moments of the equilibrium distribution, which can be accurately captured by a harmonic model, but with nonlocal sequence-dependence. We characterize the sequence-dependent choreography of backbone and base movements modulating the non-Gaussian or anharmonic effects manifested in the higher moments of the dynamics of the duplex when sampling the equilibrium distribution. Contrary to prior assumptions, such anharmonic deformations are not rare in DNA and can play a significant role in determining DNA conformation within complexes. Polymorphisms in helical geometries are particularly prevalent for certain tetranucleotide sequence contexts and are always coupled to a complex network of coordinated changes in the backbone. The analysis of our simulations, which contain instances of all tetranucleotide sequences, allow us to extend Calladine-Dickerson rules used for decades to interpret the average geometry of DNA, leading to a set of rules with quantitative predictive power that encompass nonlocal sequence-dependence and anharmonic fluctuations.

Epigenetic loss of RNA-methyltransferase NSUN5 in glioma targets ribosomes to drive a stress adaptive translational program.

Tumors have aberrant proteomes that often do not match their corresponding transcriptome profiles. One possible cause of this discrepancy is the existence of aberrant RNA modification landscapes in the so-called epitranscriptome. Here, we report that human glioma cells undergo DNA methylation-associated epigenetic silencing of NSUN5, a candidate RNA methyltransferase for 5-methylcytosine. In this setting, NSUN5 exhibits tumor-suppressor characteristics in vivo glioma models. We also found that NSUN5 loss generates an unmethylated status at the C3782 position of 28S rRNA that drives an overall depletion of protein synthesis, and leads to the emergence of an adaptive translational program for survival under conditions of cellular stress. Interestingly, NSUN5 epigenetic inactivation also renders these gliomas sensitive to bioactivatable substrates of the stress-related enzyme NQO1. Most importantly, NSUN5 epigenetic inactivation is a hallmark of glioma patients with long-term survival for this otherwise devastating disease.

Gene isolation and structural characterization of a legume tree defensin with a broad spectrum of antimicrobial activity.

MAIN CONCLUSION: The recombinant EcgDf1 defensin has an antimicrobial effect against both plant and human pathogens. In silico analyses predict that EcgDf1 is prone to form dimers capable of interacting with the membranes of microorganisms. Plant defensins comprise a large family of antimicrobial peptides (AMP) with a wide range of biological functions. They are cysteine-rich molecules, highly sequence diverse but with a conserved and stable structure. In this work, a defensin gene (EcgDf1) was isolated from Erythrina crista-galli, a legume tree native from South America. The predicted peptide presents eight cysteines, with a γ-core motif GXCX3-9C and six cysteines distributed like the typical defensin αß motif. The mature EcgDf1 coding sequence was heterologously expressed in Escherichia coli strains and purified by affinity chromatography. Possible dimer and oligomers of EcgDf1 were visible in SDS electrophoresis. Moreover, its 3D structure, determined by homology modeling, docking, and molecular dynamics simulations, was found to be compatible with the formation of homodimers between the ß3 and ß1-loop-α1, leaving the ß2-loop-ß3 free to interact with lipid membranes. The purified recombinant peptide inhibited the growth of several critical plant and human pathogens, like the opportunistic fungi Candida albicans and Aspergillus niger and the plant pathogens Clavibacter michiganensis ssp. michiganensis, Penicillium expansum, Botrytis cinerea, and Alternaria alternata. EcgDf1 is a promising candidate for the development of antimicrobial products for use in agriculture and medicine.

VeriNA3d: an R package for nucleic acids data mining.

SUMMARY: veriNA3d is an R package for the analysis of nucleic acids structural data, with an emphasis in complex RNA structures. In addition to single-structure analyses, veriNA3d also implements functions to handle whole datasets of mmCIF/PDB structures that could be retrieved from public/local repositories. Our package aims to fill a gap in the data mining of nucleic acids structures to produce flexible and high throughput analysis of structural databases. AVAILABILITY AND IMPLEMENTATION: http://mmb.irbbarcelona.org/gitlab/dgallego/veriNA3d. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

How B-DNA Dynamics Decipher Sequence-Selective Protein Recognition.

The rules governing sequence-specific DNA-protein recognition are under a long-standing debate regarding the prevalence of base versus shape readout mechanisms to explain sequence specificity and of the conformational selection versus induced fit binding paradigms to explain binding-related conformational changes in DNA. Using a combination of atomistic simulations on a subset of representative sequences and mesoscopic simulations at the protein-DNA interactome level, we demonstrate the prevalence of the shape readout model in determining sequence-specificity and of the conformational selection paradigm in defining the general mechanism for binding-related conformational changes in DNA. Our results suggest that the DNA uses a double mechanism to adapt its structure to the protein: it moves along the easiest deformation modes to approach the bioactive conformation, while final adjustments require localized rearrangements at the base-pair step and backbone level. Our study highlights the large impact of B-DNA dynamics in modulating DNA-protein binding.

Modulation of the helical properties of DNA: next-to-nearest neighbour effects and beyond.

We used extensive molecular dynamics simulations to study the structural and dynamic properties of the central d(TpA) step in the highly polymorphic d(CpTpApG) tetranucleotide. Contrary to the assumption of the dinucleotide-model and its nearest neighbours (tetranucleotide-model), the properties of the central d(TpA) step change quite significantly dependent on the next-to-nearest (hexanucleotide) sequence context and in a few cases are modulated by even remote neighbours (beyond next-to-nearest from the central TpA). Our results highlight the existence of previously undescribed dynamical mechanisms for the transmission of structural information into the DNA and demonstrate the existence of certain sequences with special physical properties that can impact on the global DNA structure and dynamics.

Glyceraldehyde-3-phosphate dehydrogenase is a chaperone that allocates labile heme in cells.

Cellular heme is thought to be distributed between a pool of sequestered heme that is tightly bound within hemeproteins and a labile heme pool required for signaling and transfer into proteins. A heme chaperone that can hold and allocate labile heme within cells has long been proposed but never been identified. Here, we show that the glycolytic protein glyceraldehyde-3-phosphate dehydrogenase (GAPDH) fulfills this role by acting as an essential repository and allocator of bioavailable heme to downstream protein targets. We identified a conserved histidine in GAPDH that is needed for its robust heme binding both in vitro and in mammalian cells. Substitution of this histidine, and the consequent decreases in GAPDH heme binding, antagonized heme delivery to both cytosolic and nuclear hemeprotein targets, including inducible nitric-oxide synthase (iNOS) in murine macrophages and the nuclear transcription factor Hap1 in yeast, even though this GAPDH variant caused cellular levels of labile heme to rise dramatically. We conclude that by virtue of its heme-binding property, GAPDH binds and chaperones labile heme to create a heme pool that is bioavailable to downstream proteins. Our finding solves a fundamental question in cell biology and provides a new foundation for exploring heme homeostasis in health and disease.

Antimicrobial and structural insights of a new snakin-like peptide isolated from Peltophorum dubium (Fabaceae).

Snakins are antimicrobial peptides (AMPs) found, so far, exclusively in plants, and known to be important in the defense against a wide range of pathogens. Like other plant AMPs, they contain several positively charged amino acids, and an even number of cysteine residues forming disulfide bridges which are considered important for their usual function. Despite its importance, studies on snakin tertiary structure and mode of action are still scarce. In this study, a new snakin-like gene was isolated from the native plant Peltophorum dubium, and its expression was verified in seedlings and adult leaves. The deduced peptide (PdSN1) shows 84% sequence identity with potato snakin-1 mature peptide, with the 12 cysteines characteristic from this peptide family at the GASA domain. The mature PdSN1 coding sequence was successfully expressed in Escherichia coli. The purified recombinant peptide inhibits the growth of important plant and human pathogens, like the economically relevant potato pathogen Streptomyces scabies and the opportunistic fungi Candida albicans and Aspergillus niger. Finally, homology and ab initio modeling techniques coupled to extensive molecular dynamics simulations were used to gain insight on the 3D structure of PdSN1, which exhibited a helix-turn-helix motif conserved in both native and recombinant peptides. We found this motif to be strongly coded in the sequence of PdSN1, as it is stable under different patterns of disulfide bonds connectivity, and even when the 12 cysteines are considered in their reduced form, explaining the previous experimental evidences.