Rational optimization of a transcription factor activation domain inhibitor.

Transcription factors are among the most attractive therapeutic targets but are considered largely 'undruggable' in part due to the intrinsically disordered nature of their activation domains. Here we show that the aromatic character of the activation domain of the androgen receptor, a therapeutic target for castration-resistant prostate cancer, is key for its activity as transcription factor, allowing it to translocate to the nucleus and partition into transcriptional condensates upon activation by androgens. On the basis of our understanding of the interactions stabilizing such condensates and of the structure that the domain adopts upon condensation, we optimized the structure of a small-molecule inhibitor previously identified by phenotypic screening. The optimized compounds had more affinity for their target, inhibited androgen-receptor-dependent transcriptional programs, and had an antitumorigenic effect in models of castration-resistant prostate cancer in cells and in vivo. These results suggest that it is possible to rationally optimize, and potentially even to design, small molecules that target the activation domains of oncogenic transcription factors.

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.

Characterization of p38α autophosphorylation inhibitors that target the non-canonical activation pathway.

p38α is a versatile protein kinase that can control numerous processes and plays important roles in the cellular responses to stress. Dysregulation of p38α signaling has been linked to several diseases including inflammation, immune disorders and cancer, suggesting that targeting p38α could be therapeutically beneficial. Over the last two decades, numerous p38α inhibitors have been developed, which showed promising effects in pre-clinical studies but results from clinical trials have been disappointing, fueling the interest in the generation of alternative mechanisms of p38α modulation. Here, we report the in silico identification of compounds that we refer to as non-canonical p38α inhibitors (NC-p38i). By combining biochemical and structural analyses, we show that NC-p38i efficiently inhibit p38α autophosphorylation but weakly affect the activity of the canonical pathway. Our results demonstrate how the structural plasticity of p38α can be leveraged to develop therapeutic opportunities targeting a subset of the functions regulated by this pathway.

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.

A method for multiplexed full-length single-molecule sequencing of the human mitochondrial genome.

Methods to reconstruct the mitochondrial DNA (mtDNA) sequence using short-read sequencing come with an inherent bias due to amplification and mapping. They can fail to determine the phase of variants, to capture multiple deletions and to cover the mitochondrial genome evenly. Here we describe a method to target, multiplex and sequence at high coverage full-length human mitochondrial genomes as native single-molecules, utilizing the RNA-guided DNA endonuclease Cas9. Combining Cas9 induced breaks, that define the mtDNA beginning and end of the sequencing reads, as barcodes, we achieve high demultiplexing specificity and delineation of the full-length of the mtDNA, regardless of the structural variant pattern. The long-read sequencing data is analysed with a pipeline where our custom-developed software, baldur, efficiently detects single nucleotide heteroplasmy to below 1%, physically determines phase and can accurately disentangle complex deletions. Our workflow is a tool for studying mtDNA variation and will accelerate mitochondrial research.

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.

Mechanism of reaction of RNA-dependent RNA polymerase from SARS-CoV-2.

We combine molecular dynamics, statistical mechanics, and hybrid quantum mechanics/molecular mechanics simulations to describe mechanistically the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) RNA-dependent RNA polymerase (RdRp). Our study analyzes the binding mode of both natural triphosphate substrates as well as remdesivir triphosphate (the active form of drug), which is bound preferentially over ATP by RdRp while being poorly recognized by human RNA polymerase II (RNA Pol II). A comparison of incorporation rates between natural and antiviral nucleotides shows that remdesivir is incorporated more slowly into the nascent RNA compared with ATP, leading to an RNA duplex that is structurally very similar to an unmodified one, arguing against the hypothesis that remdesivir is a competitive inhibitor of ATP. We characterize the entire mechanism of reaction, finding that viral RdRp is highly processive and displays a higher catalytic rate of incorporation than human RNA Pol II. Overall, our study provides the first detailed explanation of the replication mechanism of RdRp.

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.

Bioactivity descriptors for uncharacterized chemical compounds.

Chemical descriptors encode the physicochemical and structural properties of small molecules, and they are at the core of chemoinformatics. The broad release of bioactivity data has prompted enriched representations of compounds, reaching beyond chemical structures and capturing their known biological properties. Unfortunately, bioactivity descriptors are not available for most small molecules, which limits their applicability to a few thousand well characterized compounds. Here we present a collection of deep neural networks able to infer bioactivity signatures for any compound of interest, even when little or no experimental information is available for them. Our signaturizers relate to bioactivities of 25 different types (including target profiles, cellular response and clinical outcomes) and can be used as drop-in replacements for chemical descriptors in day-to-day chemoinformatics tasks. Indeed, we illustrate how inferred bioactivity signatures are useful to navigate the chemical space in a biologically relevant manner, unveiling higher-order organization in natural product collections, and to enrich mostly uncharacterized chemical libraries for activity against the drug-orphan target Snail1. Moreover, we implement a battery of signature-activity relationship (SigAR) models and show a substantial improvement in performance, with respect to chemistry-based classifiers, across a series of biophysics and physiology activity prediction benchmarks.

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.

Bioactivity Profile Similarities to Expand the Repertoire of COVID-19 Drugs.

Until a vaccine becomes available, the current repertoire of drugs is our only therapeutic asset to fight the SARS-CoV-2 outbreak. Indeed, emergency clinical trials have been launched to assess the effectiveness of many marketed drugs, tackling the decrease of viral load through several mechanisms. Here, we present an online resource, based on small-molecule bioactivity signatures and natural language processing, to expand the portfolio of compounds with potential to treat COVID-19. By comparing the set of drugs reported to be potentially active against SARS-CoV-2 to a universe of 1 million bioactive molecules, we identify compounds that display analogous chemical and functional features to the current COVID-19 candidates. Searches can be filtered by level of evidence and mechanism of action, and results can be restricted to drug molecules or include the much broader space of bioactive compounds. Moreover, we allow users to contribute COVID-19 drug candidates, which are automatically incorporated to the pipeline once per day. The computational platform, as well as the source code, is available at https://sbnb.irbbarcelona.org/covid19.

Author Correction: Hsp70 and Hsp40 inhibit an inter-domain interaction necessary for transcriptional activity in the androgen receptor.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Hsp70 and Hsp40 inhibit an inter-domain interaction necessary for transcriptional activity in the androgen receptor.

Molecular chaperones such as Hsp40 and Hsp70 hold the androgen receptor (AR) in an inactive conformation. They are released in the presence of androgens, enabling transactivation and causing the receptor to become aggregation-prone. Here we show that these molecular chaperones recognize a region of the AR N-terminal domain (NTD), including a FQNLF motif, that interacts with the AR ligand-binding domain (LBD) upon activation. This suggests that competition between molecular chaperones and the LBD for the FQNLF motif regulates AR activation. We also show that, while the free NTD oligomerizes, binding to Hsp70 increases its solubility. Stabilizing the NTD-Hsp70 interaction with small molecules reduces AR aggregation and promotes its degradation in cellular and mouse models of the neuromuscular disorder spinal bulbar muscular atrophy. These results help resolve the mechanisms by which molecular chaperones regulate the balance between AR aggregation, activation and quality control.

Oncogenic mutations at the EGFR ectodomain structurally converge to remove a steric hindrance on a kinase-coupled cryptic epitope.

Epidermal growth factor receptor (EGFR) signaling is initiated by a large ligand-favored conformational change of the extracellular domain (ECD) from a closed, self-inhibited tethered monomer, to an open untethered state, which exposes a loop required for strong dimerization and activation. In glioblastomas (GBMs), structurally heterogeneous missense and deletion mutations concentrate at the ECD for unclear reasons. We explore the conformational impact of GBM missense mutations, combining elastic network models (ENMs) with multiple molecular dynamics (MD) trajectories. Our simulations reveal that the main missense class, located at the I-II interface away from the self-inhibitory tether, can unexpectedly favor spontaneous untethering to a compact intermediate state, here validated by small-angle X-ray scattering (SAXS). Significantly, such intermediate is characterized by the rotation of a large ECD fragment (N-TR1), deleted in the most common GBM mutation, EGFRvIII, and that makes accessible a cryptic epitope characteristic of cancer cells. This observation suggested potential structural equivalence of missense and deletion ECD changes in GBMs. Corroborating this hypothesis, our FACS, in vitro, and in vivo data demonstrate that entirely different ECD variants all converge to remove N-TR1 steric hindrance from the 806-epitope, which we show is allosterically coupled to an intermediate kinase and hallmarks increased oncogenicity. Finally, the detected extraintracellular coupling allows for synergistic cotargeting of the intermediate with mAb806 and inhibitors, which is proved herein.

A multifunctional toolkit for target-directed cancer therapy.

Here we present 2shRNA, a shRNA-based nanobinder, which can simultaneously attack two therapeutic targets involved in drug resistance pathways and can additionally bind accessory molecules such as cell targeting peptides or fluorophores. We create 2shRNAs designed to specifically kill HER2+ breast cancer cells in the absence of a transfecting agent.

Regulation of Androgen Receptor Activity by Transient Interactions of Its Transactivation Domain with General Transcription Regulators.

The androgen receptor is a transcription factor that plays a key role in the development of prostate cancer, and its interactions with general transcription regulators are therefore of potential therapeutic interest. The mechanistic basis of these interactions is poorly understood due to the intrinsically disordered nature of the transactivation domain of the androgen receptor and the generally transient nature of the protein-protein interactions that trigger transcription. Here, we identify a motif of the transactivation domain that contributes to transcriptional activity by recruiting the C-terminal domain of subunit 1 of the general transcription regulator TFIIF. These findings provide molecular insights into the regulation of androgen receptor function and suggest strategies for treating castration-resistant prostate cancer.

Efficient siRNA-peptide conjugation for specific targeted delivery into tumor cells.

Despite the broad applicability of the Huisgen cycloaddition reaction, the click functionalization of RNAs with peptides still remains a challenge. Here we describe a straightforward method for the click functionalization of siRNAs with peptides of different sizes and complexities. Among them, a promising peptide carrier for the selective siRNA delivery into HER2+ breast cancer cell lines has been reported.

Rational design of novel N-alkyl-N capped biostable RNA nanostructures for efficient long-term inhibition of gene expression.

Computational techniques have been used to design a novel class of RNA architecture with expected improved resistance to nuclease degradation, while showing interference RNA activity. The in silico designed structure consists of a 24-29 bp duplex RNA region linked on both ends by N-alkyl-N dimeric nucleotides (BCn dimers; n = number of carbon atoms of the alkyl chain). A series of N-alkyl-N capped dumbbell-shaped structures were efficiently synthesized by double ligation of BCn-loop hairpins. The resulting BCn-loop dumbbells displayed experimentally higher biostability than their 3'-N-alkyl-N linear version, and were active against a range of mRNA targets. We studied first the effect of the alkyl chain and stem lengths on RNAi activity in a screen involving two series of dumbbell analogues targeting Renilla and Firefly luciferase genes. The best dumbbell design (containing BC6 loops and 29 bp) was successfully used to silence GRB7 expression in HER2+ breast cancer cells for longer periods of time than natural siRNAs and known biostable dumbbells. This BC6-loop dumbbell-shaped structure displayed greater anti-proliferative activity than natural siRNAs.

Whole-genome bisulfite sequencing of two distinct interconvertible DNA methylomes of mouse embryonic stem cells.

The use of two kinase inhibitors (2i) enables derivation of mouse embryonic stem cells (ESCs) in the pluripotent ground state. Using whole-genome bisulfite sequencing (WGBS), we show that male 2i ESCs are globally hypomethylated compared to conventional ESCs maintained in serum. In serum, female ESCs are hypomethyated similarly to male ESCs in 2i, and DNA methylation is further reduced in 2i. Regions with elevated DNA methylation in 2i strongly correlate with the presence of H3K9me3 on endogenous retroviruses (ERVs) and imprinted loci. The methylome of male ESCs in serum parallels postimplantation blastocyst cells, while 2i stalls ESCs in a hypomethylated, ICM-like state. WGBS analysis during adaptation of 2i ESCs to serum suggests that deposition of DNA methylation is largely random, while loss of DNA methylation during reversion to 2i occurs passively, initiating at TET1 binding sites. Together, our analysis provides insight into DNA methylation dynamics in cultured ESCs paralleling early developmental processes.