Deficient Endoplasmic Reticulum-Mitochondrial Phosphatidylserine Transfer Causes Liver Disease.

Non-alcoholic fatty liver is the most common liver disease worldwide. Here, we show that the mitochondrial protein mitofusin 2 (Mfn2) protects against liver disease. Reduced Mfn2 expression was detected in liver biopsies from patients with non-alcoholic steatohepatitis (NASH). Moreover, reduced Mfn2 levels were detected in mouse models of steatosis or NASH, and its re-expression in a NASH mouse model ameliorated the disease. Liver-specific ablation of Mfn2 in mice provoked inflammation, triglyceride accumulation, fibrosis, and liver cancer. We demonstrate that Mfn2 binds phosphatidylserine (PS) and can specifically extract PS into membrane domains, favoring PS transfer to mitochondria and mitochondrial phosphatidylethanolamine (PE) synthesis. Consequently, hepatic Mfn2 deficiency reduces PS transfer and phospholipid synthesis, leading to endoplasmic reticulum (ER) stress and the development of a NASH-like phenotype and liver cancer. Ablation of Mfn2 in liver reveals that disruption of ER-mitochondrial PS transfer is a new mechanism involved in the development of liver disease.

4D Genome Rewiring during Oncogene-Induced and Replicative Senescence.

To understand the role of the extensive senescence-associated 3D genome reorganization, we generated genome-wide chromatin interaction maps, epigenome, replication-timing, whole-genome bisulfite sequencing, and gene expression profiles from cells entering replicative senescence (RS) or upon oncogene-induced senescence (OIS). We identify senescence-associated heterochromatin domains (SAHDs). Differential intra- versus inter-SAHD interactions lead to the formation of senescence-associated heterochromatin foci (SAHFs) in OIS but not in RS. This OIS-specific configuration brings active genes located in genomic regions adjacent to SAHDs in close spatial proximity and favors their expression. We also identify DNMT1 as a factor that induces SAHFs by promoting HMGA2 expression. Upon DNMT1 depletion, OIS cells transition to a 3D genome conformation akin to that of cells in replicative senescence. These data show how multi-omics and imaging can identify critical features of RS and OIS and discover determinants of acute senescence and SAHF formation.

A new paradigm for molecular dynamics databases: the COVID-19 database, the legacy of a titanic community effort.

Molecular dynamics (MD) simulations are keeping computers busy around the world, generating a huge amount of data that is typically not open to the scientific community. Pioneering efforts to ensure the safety and reusability of MD data have been based on the use of simple databases providing a limited set of standard analyses on single-short trajectories. Despite their value, these databases do not offer a true solution for the current community of MD users, who want a flexible analysis pipeline and the possibility to address huge non-Markovian ensembles of large systems. Here we present a new paradigm for MD databases, resilient to large systems and long trajectories, and designed to be compatible with modern MD simulations. The data are offered to the community through a web-based graphical user interface (GUI), implemented with state-of-the-art technology, which incorporates system-specific analysis designed by the trajectory providers. A REST API and associated Jupyter Notebooks are integrated into the platform, allowing fully customized meta-analysis by final users. The new technology is illustrated using a collection of trajectories obtained by the community in the context of the effort to fight the COVID-19 pandemic. The server is accessible at https://bioexcel-cv19.bsc.es/#/. It is free and open to all users and there are no login requirements. It is also integrated into the simulations section of the BioExcel-MolSSI COVID-19 Molecular Structure and Therapeutics Hub: https://covid.molssi.org/simulations/ and is part of the MDDB effort (https://mddbr.eu).

CGeNArate: a sequence-dependent coarse-grained model of DNA for accurate atomistic MD simulations of kb-long duplexes.

We present CGeNArate, a new model for molecular dynamics simulations of very long segments of B-DNA in the context of biotechnological or chromatin studies. The developed method uses a coarse-grained Hamiltonian with trajectories that are back-mapped to the atomistic resolution level with extreme accuracy by means of Machine Learning Approaches. The method is sequence-dependent and reproduces very well not only local, but also global physical properties of DNA. The efficiency of the method allows us to recover with a reduced computational effort high-quality atomic-resolution ensembles of segments containing many kilobases of DNA, entering into the gene range or even the entire DNA of certain cellular organelles.

An integrated machine-learning model to predict nucleosome architecture.

We demonstrate that nucleosomes placed in the gene body can be accurately located from signal decay theory assuming two emitters located at the beginning and at the end of genes. These generated wave signals can be in phase (leading to well defined nucleosome arrays) or in antiphase (leading to fuzzy nucleosome architectures). We found that the first (+1) and the last (-last) nucleosomes are contiguous to regions signaled by transcription factor binding sites and unusual DNA physical properties that hinder nucleosome wrapping. Based on these analyses, we developed a method that combines Machine Learning and signal transmission theory able to predict the basal locations of the nucleosomes with an accuracy similar to that of experimental MNase-seq based methods.

Formation of left-handed helices by C2'-fluorinated nucleic acids under physiological salt conditions.

Recent findings in cell biology have rekindled interest in Z-DNA, the left-handed helical form of DNA. We report here that two minimally modified nucleosides, 2'F-araC and 2'F-riboG, induce the formation of the Z-form under low ionic strength. We show that oligomers entirely made of these two nucleosides exclusively produce left-handed duplexes that bind to the Zα domain of ADAR1. The effect of the two nucleotides is so dramatic that Z-form duplexes are the only species observed in 10 mM sodium phosphate buffer and neutral pH, and no B-form is observed at any temperature. Hence, in contrast to other studies reporting formation of Z/B-form equilibria by a preference for purine glycosidic angles in syn, our NMR and computational work revealed that sequential 2'F…H2N and intramolecular 3'H…N3' interactions stabilize the left-handed helix. The equilibrium between B- and Z- forms is slow in the 19F NMR time scale (≥ms), and each conformation exhibited unprecedented chemical shift differences in the 19F signals. This observation led to a reliable estimation of the relative population of B and Z species and enabled us to monitor B-Z transitions under different conditions. The unique features of 2'F-modified DNA should thus be a valuable addition to existing techniques for specific detection of new Z-binding proteins and ligands.

Site-specific incorporation of a fluorescent nucleobase analog enhances i-motif stability and allows monitoring of i-motif folding inside cells.

The i-motif is an intriguing non-canonical DNA structure, whose role in the cell is still controversial. Development of methods to study i-motif formation under physiological conditions in living cells is necessary to study its potential biological functions. The cytosine analog 1,3-diaza-2-oxophenoxazine (tCO) is a fluorescent nucleobase able to form either hemiprotonated base pairs with cytosine residues, or neutral base pairs with guanines. We show here that when tCO is incorporated in the proximity of a G:C:G:C minor groove tetrad, it induces a strong thermal and pH stabilization, resulting in i-motifs with Tm of 39ºC at neutral pH. The structural determination by NMR methods reveals that the enhanced stability is due to a large stacking interaction between the guanines of the tetrad with the tCO nucleobase, which forms a tCO:C+ in the folded structure at unusually-high pHs, leading to an increased quenching in its fluorescence at neutral conditions. This quenching is much lower when tCO is base-paired to guanines and totally disappears when the oligonucleotide is unfolded. By taking profit of this property, we have been able to monitor i-motif folding in cells.

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.

Using interactive Jupyter Notebooks and BioConda for FAIR and reproducible biomolecular simulation workflows.

Interactive Jupyter Notebooks in combination with Conda environments can be used to generate FAIR (Findable, Accessible, Interoperable and Reusable/Reproducible) biomolecular simulation workflows. The interactive programming code accompanied by documentation and the possibility to inspect intermediate results with versatile graphical charts and data visualization is very helpful, especially in iterative processes, where parameters might be adjusted to a particular system of interest. This work presents a collection of FAIR notebooks covering various areas of the biomolecular simulation field, such as molecular dynamics (MD), protein-ligand docking, molecular checking/modeling, molecular interactions, and free energy perturbations. Workflows can be launched with myBinder or easily installed in a local system. The collection of notebooks aims to provide a compilation of demonstration workflows, and it is continuously updated and expanded with examples using new methodologies and tools.

Repair of UV-Induced DNA Damage Independent of Nucleotide Excision Repair Is Masked by MUTYH.

DNA lesions caused by UV damage are thought to be repaired solely by the nucleotide excision repair (NER) pathway in human cells. Patients carrying mutations within genes functioning in this pathway display a range of pathologies, including an increased susceptibility to cancer, premature aging, and neurological defects. There are currently no curative therapies available. Here we performed a high-throughput chemical screen for agents that could alleviate the cellular sensitivity of NER-deficient cells to UV-induced DNA damage. This led to the identification of the clinically approved anti-diabetic drug acetohexamide, which promoted clearance of UV-induced DNA damage without the accumulation of chromosomal aberrations, hence promoting cellular survival. Acetohexamide exerted this protective function by antagonizing expression of the DNA glycosylase, MUTYH. Together, our data reveal the existence of an NER-independent mechanism to remove UV-induced DNA damage and prevent cell death.

Correlated motions in DNA: beyond base-pair step models of DNA flexibility.

Traditional mesoscopic models of DNA flexibility use a reductionist-local approach, which assumes that the flexibility of DNA can be expressed as local harmonic movements (at the base-pair step level) in the helical space, ignoring multimodality and correlations in DNA movements, which have in reality a large impact in modulating DNA movements. We present a new multimodal-harmonic correlated model, which takes both contributions into account, providing, with a small computational cost, results of an unprecedented local and global quality. The accuracy of this method and its computational efficiency make it an alternative to explore the dynamics of long segments of DNA, approaching the chromatin range.

Structural analysis of the Candida albicans mitochondrial DNA maintenance factor Gcf1p reveals a dynamic DNA-bridging mechanism.

The compaction of mitochondrial DNA (mtDNA) is regulated by architectural HMG-box proteins whose limited cross-species similarity suggests diverse underlying mechanisms. Viability of Candida albicans, a human antibiotic-resistant mucosal pathogen, is compromised by altering mtDNA regulators. Among them, there is the mtDNA maintenance factor Gcf1p, which differs in sequence and structure from its human and Saccharomyces cerevisiae counterparts, TFAM and Abf2p. Our crystallographic, biophysical, biochemical and computational analysis showed that Gcf1p forms dynamic protein/DNA multimers by a combined action of an N-terminal unstructured tail and a long helix. Furthermore, an HMG-box domain canonically binds the minor groove and dramatically bends the DNA while, unprecedentedly, a second HMG-box binds the major groove without imposing distortions. This architectural protein thus uses its multiple domains to bridge co-aligned DNA segments without altering the DNA topology, revealing a new mechanism of mtDNA condensation.

Controlled sulfur-based engineering confers mouldability to phosphorothioate antisense oligonucleotides.

Phosphorothioates (PS) have proven their effectiveness in the area of therapeutic oligonucleotides with applications spanning from cancer treatment to neurodegenerative disorders. Initially, PS substitution was introduced for the antisense oligonucleotides (PS ASOs) because it confers an increased nuclease resistance meanwhile ameliorates cellular uptake and in-vivo bioavailability. Thus, PS oligonucleotides have been elevated to a fundamental asset in the realm of gene silencing therapeutic methodologies. But, despite their wide use, little is known on the possibly different structural changes PS-substitutions may provoke in DNA·RNA hybrids. Additionally, scarce information and significant controversy exists on the role of phosphorothioate chirality in modulating PS properties. Here, through comprehensive computational investigations and experimental measurements, we shed light on the impact of PS chirality in DNA-based antisense oligonucleotides; how the different phosphorothioate diastereomers impact DNA topology, stability and flexibility to ultimately disclose pro-Sp S and pro-Rp S roles at the catalytic core of DNA Exonuclease and Human Ribonuclease H; two major obstacles in ASOs-based therapies. Altogether, our results provide full-atom and mechanistic insights on the structural aberrations PS-substitutions provoke and explain the origin of nuclease resistance PS-linkages confer to DNA·RNA hybrids; crucial information to improve current ASOs-based therapies.

Correction: The structural role of SARS-CoV-2 genetic background in the emergence and success of spike mutations: The case of the spike A222V mutation.

[This corrects the article DOI: 10.1371/journal.ppat.1010631.].

The structural role of SARS-CoV-2 genetic background in the emergence and success of spike mutations: The case of the spike A222V mutation.

The S:A222V point mutation, within the G clade, was characteristic of the 20E (EU1) SARS-CoV-2 variant identified in Spain in early summer 2020. This mutation has since reappeared in the Delta subvariant AY.4.2, raising questions about its specific effect on viral infection. We report combined serological, functional, structural and computational studies characterizing the impact of this mutation. Our results reveal that S:A222V promotes an increased RBD opening and slightly increases ACE2 binding as compared to the parent S:D614G clade. Finally, S:A222V does not reduce sera neutralization capacity, suggesting it does not affect vaccine effectiveness.

DNAffinity: a machine-learning approach to predict DNA binding affinities of transcription factors.

We present a physics-based machine learning approach to predict in vitro transcription factor binding affinities from structural and mechanical DNA properties directly derived from atomistic molecular dynamics simulations. The method is able to predict affinities obtained with techniques as different as uPBM, gcPBM and HT-SELEX with an excellent performance, much better than existing algorithms. Due to its nature, the method can be extended to epigenetic variants, mismatches, mutations, or any non-coding nucleobases. When complemented with chromatin structure information, our in vitro trained method provides also good estimates of in vivo binding sites in yeast.

BioExcel Building Blocks Workflows (BioBB-Wfs), an integrated web-based platform for biomolecular simulations.

We present BioExcel Building Blocks Workflows, a web-based graphical user interface (GUI) offering access to a collection of transversal pre-configured biomolecular simulation workflows assembled with the BioExcel Building Blocks library. Available workflows include Molecular Dynamics setup, protein-ligand docking, trajectory analyses and small molecule parameterization. Workflows can be launched in the platform or downloaded to be run in the users' own premises. Remote launching of long executions to user's available High-Performance computers is possible, only requiring configuration of the appropriate access credentials. The web-based graphical user interface offers a high level of interactivity, with integration with the NGL viewer to visualize and check 3D structures, MDsrv to visualize trajectories, and Plotly to explore 2D plots. The server requires no login but is recommended to store the users' projects and manage sensitive information such as remote credentials. Private projects can be made public and shared with colleagues with a simple URL. The tool will help biomolecular simulation users with the most common and repetitive processes by means of a very intuitive and interactive graphical user interface. The server is accessible at https://mmb.irbbarcelona.org/biobb-wfs.

The Rhodium Analogue of Coenzyme B12 as an Anti-Photoregulatory Ligand Inhibiting Bacterial CarH Photoreceptors.

Coenzyme B12 (AdoCbl; 5'-deoxy-5'-adenosylcobalamin), the quintessential biological organometallic radical catalyst, has a formerly unanticipated, yet extensive, role in photoregulation in bacteria. The light-responsive cobalt-corrin AdoCbl performs this nonenzymatic role by facilitating the assembly of CarH photoreceptors into DNA-binding tetramers in the dark, suppressing gene expression. Conversely, exposure to light triggers the decomposition of this AdoCbl-bound complex by a still elusive photochemical mechanism, activating gene expression. Here, we have examined AdoRhbl, the non-natural rhodium analogue of AdoCbl, as a photostable isostructural surrogate for AdoCbl. We show that AdoRhbl closely emulates AdoCbl in its uptake by bacterial cells and structural functionality as a regulatory ligand for CarH tetramerization, DNA binding, and repressor activity. Remarkably, we find AdoRhbl is photostable even when bound "base-off/His-on" to CarH in vitro and in vivo. Thus, AdoRhbl, an antivitamin B12, also represents an unprecedented anti-photoregulatory ligand, opening a pathway to precisely target biomimetic inhibition of AdoCbl-based photoregulation, with new possibilities for selective antibacterial applications. Computational biomolecular analysis of AdoRhbl binding to CarH yields detailed structural insights into this complex, which suggest that the adenosyl group of photoexcited AdoCbl bound to CarH may specifically undergo a concerted non-radical syn-1,2-elimination mechanism, an aspect not previously considered for this photoreceptor.

pH-Dependent Capping Interactions Induce Large-Scale Structural Transitions in i-Motifs.

We study here a DNA oligonucleotide having the ability to form two different i-motif structures whose relative stability depends on pH and temperature. The major species at neutral pH is stabilized by two C:C+ base pairs capped by two minor groove G:C:G:C tetrads. The high pH and thermal stability of this structure are mainly due to the favorable effect of the minor groove tetrads on their adjacent positively charged C:C+ base pairs. At pH 5, we observe a more elongated i-motif structure consisting of four C:C+ base pairs capped by two G:T:G:T tetrads. Molecular dynamics calculations show that the conformational transition between the two structures is driven by the protonation state of key cytosines. In spite of large conformational differences, the transition between the acidic and neutral structures can occur without unfolding of the i-motif. These results represent the first case of a conformational switch between two different i-motif structures and illustrate the dramatic pH-dependent plasticity of this fascinating DNA motif.

BioExcel Building Blocks REST API (BioBB REST API), programmatic access to interoperable biomolecular simulation tools.

MOTIVATION: The BioExcel Building Blocks (BioBB) library offers a broad collection of wrappers on top of common biomolecular simulation and bioinformatics tools. The possibility to access the library remotely and programmatically increases its usability, allowing individual and sporadic executions and enabling remote workflows. RESULTS: BioBB REST API extends and complements the BioBB library offering programmatic access to the collection of biomolecular simulation tools included in the BioExcel Building Blocks library. Molecular Dynamics setup, docking, structure modeling, free energy simulations and flexibility analyses are examples of functionalities included in the endpoints collection. All functionalities are accessible through standard REST API calls, voiding the need for tool installation. AVAILABILITY AND IMPLEMENTATION: All the information related to the BioBB REST API endpoints is accessible from https://mmb.irbbarcelona.org/biobb-api/. Links to extended documentation, including OpenAPI endpoints specification and examples, Read-The-Docs documentation and a complete workflow tutorial can be found in the Supplementary Table S1. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.