AI-guided pipeline for protein-protein interaction drug discovery identifies a SARS-CoV-2 inhibitor.

Protein-protein interactions (PPIs) offer great opportunities to expand the druggable proteome and therapeutically tackle various diseases, but remain challenging targets for drug discovery. Here, we provide a comprehensive pipeline that combines experimental and computational tools to identify and validate PPI targets and perform early-stage drug discovery. We have developed a machine learning approach that prioritizes interactions by analyzing quantitative data from binary PPI assays or AlphaFold-Multimer predictions. Using the quantitative assay LuTHy together with our machine learning algorithm, we identified high-confidence interactions among SARS-CoV-2 proteins for which we predicted three-dimensional structures using AlphaFold-Multimer. We employed VirtualFlow to target the contact interface of the NSP10-NSP16 SARS-CoV-2 methyltransferase complex by ultra-large virtual drug screening. Thereby, we identified a compound that binds to NSP10 and inhibits its interaction with NSP16, while also disrupting the methyltransferase activity of the complex, and SARS-CoV-2 replication. Overall, this pipeline will help to prioritize PPI targets to accelerate the discovery of early-stage drug candidates targeting protein complexes and pathways.

Despite the odds: formation of the SARS-CoV-2 methylation complex.

Coronaviruses modify their single-stranded RNA genome with a methylated cap during replication to mimic the eukaryotic mRNAs. The capping process is initiated by several nonstructural proteins (nsp) encoded in the viral genome. The methylation is performed by two methyltransferases, nsp14 and nsp16, while nsp10 acts as a co-factor to both. Additionally, nsp14 carries an exonuclease domain which operates in the proofreading system during RNA replication of the viral genome. Both nsp14 and nsp16 were reported to independently bind nsp10, but the available structural information suggests that the concomitant interaction between these three proteins would be impossible due to steric clashes. Here, we show that nsp14, nsp10, and nsp16 can form a heterotrimer complex upon significant allosteric change. This interaction is expected to encourage the formation of mature capped viral mRNA, modulating nsp14's exonuclease activity, and protecting the viral RNA. Our findings show that nsp14 is amenable to allosteric regulation and may serve as a novel target for therapeutic approaches.

AI-guided pipeline for protein-protein interaction drug discovery identifies a SARS-CoV-2 inhibitor.

Protein-protein interactions (PPIs) offer great opportunities to expand the druggable proteome and therapeutically tackle various diseases, but remain challenging targets for drug discovery. Here, we provide a comprehensive pipeline that combines experimental and computational tools to identify and validate PPI targets and perform early-stage drug discovery. We have developed a machine learning approach that prioritizes interactions by analyzing quantitative data from binary PPI assays and AlphaFold-Multimer predictions. Using the quantitative assay LuTHy together with our machine learning algorithm, we identified high-confidence interactions among SARS-CoV-2 proteins for which we predicted three-dimensional structures using AlphaFold Multimer. We employed VirtualFlow to target the contact interface of the NSP10-NSP16 SARS-CoV-2 methyltransferase complex by ultra-large virtual drug screening. Thereby, we identified a compound that binds to NSP10 and inhibits its interaction with NSP16, while also disrupting the methyltransferase activity of the complex, and SARS-CoV-2 replication. Overall, this pipeline will help to prioritize PPI targets to accelerate the discovery of early-stage drug candidates targeting protein complexes and pathways.

SLAM-Drop-seq reveals mRNA kinetic rates throughout the cell cycle.

RNA abundance is tightly regulated in eukaryotic cells by modulating the kinetic rates of RNA production, processing, and degradation. To date, little is known about time­dependent kinetic rates during dynamic processes. Here, we present SLAM­Drop­seq, a method that combines RNA metabolic labeling and alkylation of modified nucleotides in methanol­fixed cells with droplet­based sequencing to detect newly synthesized and preexisting mRNAs in single cells. As a first application, we sequenced 7280 HEK293 cells and calculated gene­specific kinetic rates during the cell cycle using the novel package Eskrate. Of the 377 robust­cycling genes that we identified, only a minor fraction is regulated solely by either dynamic transcription or degradation (6 and 4%, respectively). By contrast, the vast majority (89%) exhibit dynamically regulated transcription and degradation rates during the cell cycle. Our study thus shows that temporally regulated mRNA degradation is fundamental for the correct expression of a majority of cycling genes. SLAM­Drop­seq, combined with Eskrate, is a powerful approach to understanding the underlying mRNA kinetics of single­cell gene expression dynamics in continuous biological processes.

HDLBP binds ER-targeted mRNAs by multivalent interactions to promote protein synthesis of transmembrane and secreted proteins.

The biological role of RNA-binding proteins in the secretory pathway is not well established. Here, we describe that human HDLBP/Vigilin directly interacts with more than 80% of ER-localized mRNAs. PAR-CLIP analysis reveals that these transcripts represent high affinity HDLBP substrates and are specifically bound in their coding sequences (CDS), in contrast to CDS/3'UTR-bound cytosolic mRNAs. HDLBP crosslinks strongly to long CU-rich motifs, which frequently reside in CDS of ER-localized mRNAs and result in high affinity multivalent interactions. In addition to HDLBP-ncRNA interactome, quantification of HDLBP-proximal proteome confirms association with components of the translational apparatus and the signal recognition particle. Absence of HDLBP results in decreased translation efficiency of HDLBP target mRNAs, impaired protein synthesis and secretion in model cell lines, as well as decreased tumor growth in a lung cancer mouse model. These results highlight a general function for HDLBP in the translation of ER-localized mRNAs and its relevance for tumor progression.

Loss of m1acp3Ψ Ribosomal RNA Modification Is a Major Feature of Cancer.

The ribosome is an RNA-protein complex that is essential for translation in all domains of life. The structural and catalytic core of the ribosome is its ribosomal RNA (rRNA). While mutations in ribosomal protein (RP) genes are known drivers of oncogenesis, oncogenic rRNA variants have remained elusive. We identify a cancer-specific single-nucleotide variation in 18S rRNA at nucleotide 1248.U in up to 45.9% of patients with colorectal carcinoma (CRC) and present across >22 cancer types. This is the site of a unique hyper-modified base, 1-methyl-3-α-amino-α-carboxyl-propyl pseudouridine (m1acp3Ψ), a >1-billion-years-conserved RNA modification at the peptidyl decoding site of the ribosome. A subset of CRC tumors we call hypo-m1acp3Ψ shows sub-stoichiometric m1acp3Ψ modification, unlike normal control tissues. An m1acp3Ψ knockout model and hypo-m1acp3Ψ patient tumors share a translational signature characterized by highly abundant ribosomal proteins. Thus, m1acp3Ψ-deficient rRNA forms an uncharacterized class of "onco-ribosome" which may serve as a chemotherapeutic target for treating cancer patients.

The suppressive cap-binding complex factor 4EIP is required for normal differentiation.

Trypanosoma brucei live in mammals as bloodstream forms and in the Tsetse midgut as procyclic forms. Differentiation from one form to the other proceeds via a growth-arrested stumpy form with low messenger RNA (mRNA) content and translation. The parasites have six eIF4Es and five eIF4Gs. EIF4E1 pairs with the mRNA-binding protein 4EIP but not with any EIF4G. EIF4E1 and 4EIP each inhibit expression when tethered to a reporter mRNA, but while tethered EIF4E1 suppresses only when 4EIP is present, suppression by tethered 4EIP does not require the interaction with EIF4E1. In growing bloodstream forms, 4EIP is preferentially associated with unstable mRNAs. Bloodstream- or procyclic-form trypanosomes lacking 4EIP have only a marginal growth disadvantage. Bloodstream forms without 4EIP are, however, defective in translation suppression during stumpy-form differentiation and cannot subsequently convert to growing procyclic forms. Intriguingly, the differentiation defect can be complemented by a truncated 4EIP that does not interact with EIF4E1. In contrast, bloodstream forms lacking EIF4E1 have a growth defect, stumpy formation seems normal, but they appear unable to grow as procyclic forms. We suggest that 4EIP and EIF4E1 fine-tune mRNA levels in growing cells, and that 4EIP contributes to translation suppression during differentiation to the stumpy form.

Translation Regulation and RNA Granule Formation after Heat Shock of Procyclic Form Trypanosoma brucei: Many Heat-Induced mRNAs Are also Increased during Differentiation to Mammalian-Infective Forms.

African trypanosome procyclic forms multiply in the midgut of tsetse flies, and are routinely cultured at 27°C. Heat shocks of 37°C and above result in general inhibition of translation, and severe heat shock (41°C) results in sequestration of mRNA in granules. The mRNAs that are bound by the zinc-finger protein ZC3H11, including those encoding refolding chaperones, escape heat-induced translation inhibition. At 27°C, ZC3H11 mRNA is predominantly present as an untranslated cytosolic messenger ribonucleoprotein particle, but after heat shocks of 37°C-41°C, the ZC3H11 mRNA moves into the polysomal fraction. To investigate the scope and specificities of heat-shock translational regulation and granule formation, we analysed the distributions of mRNAs on polysomes at 27°C and after 1 hour at 39°C, and the mRNA content of 41°C heat shock granules. We found that mRNAs that bind to ZC3H11 remained in polysomes at 39°C and were protected from sequestration in granules at 41°C. As previously seen for starvation stress granules, the mRNAs that encode ribosomal proteins were excluded from heat-shock granules. 70 mRNAs moved towards the polysomal fraction after the 39°C heat shock, and 260 increased in relative abundance. Surprisingly, many of these mRNAs are also increased when trypanosomes migrate to the tsetse salivary glands. It therefore seems possible that in the wild, temperature changes due to diurnal variations and periodic intake of warm blood might influence the efficiency with which procyclic forms develop into mammalian-infective forms.

Integrative analysis of the Trypanosoma brucei gene expression cascade predicts differential regulation of mRNA processing and unusual control of ribosomal protein expression.

BACKGROUND: Trypanosoma brucei is a unicellular parasite which multiplies in mammals (bloodstream form) and Tsetse flies (procyclic form). Trypanosome RNA polymerase II transcription is polycistronic, individual mRNAs being excised by trans splicing and polyadenylation. We previously made detailed measurements of mRNA half-lives in bloodstream and procyclic forms, and developed a mathematical model of gene expression for bloodstream forms. At the whole transcriptome level, many bloodstream-form mRNAs were less abundant than was predicted by the model. RESULTS: We refined the published mathematical model and extended it to the procyclic form. We used the model, together with known mRNA half-lives, to predict the abundances of individual mRNAs, assuming rapid, unregulated mRNA processing; then we compared the results with measured mRNA abundances. Remarkably, the abundances of most mRNAs in procyclic forms are predicted quite well by the model, being largely explained by variations in mRNA decay rates and length. In bloodstream forms substantially more mRNAs are less abundant than predicted. We list mRNAs that are likely to show particularly slow or inefficient processing, either in both forms or with developmental regulation. We also measured ribosome occupancies of all mRNAs in trypanosomes grown in the same conditions as were used to measure mRNA turnover. In procyclic forms there was a weak positive correlation between ribosome density and mRNA half-life, suggesting cross-talk between translation and mRNA decay; ribosome density was related to the proportion of the mRNA on polysomes, indicating control of translation initiation. Ribosomal protein mRNAs in procyclics appeared to be exceptionally rapidly processed but poorly translated. CONCLUSIONS: Levels of mRNAs in procyclic form trypanosomes are determined mainly by length and mRNA decay, with some control of precursor processing. In bloodstream forms variations in nuclear events play a larger role in transcriptome regulation, suggesting aquisition of new control mechanisms during adaptation to mammalian parasitism.

Regulating a Post-Transcriptional Regulator: Protein Phosphorylation, Degradation and Translational Blockage in Control of the Trypanosome Stress-Response RNA-Binding Protein ZC3H11.

The life cycle of the mammalian pathogen Trypanosoma brucei involves commuting between two markedly different environments: the homeothermic mammalian host and the poikilothermic invertebrate vector. The ability to resist temperature and other stresses is essential for trypanosome survival. Trypanosome gene expression is mainly post-transcriptional, but must nevertheless be adjusted in response to environmental cues, including host-specific physical and chemical stresses. We investigate here the control of ZC3H11, a CCCH zinc finger protein which stabilizes stress response mRNAs. ZC3H11 protein levels increase at least 10-fold when trypanosomes are stressed by heat shock, proteasome inhibitors, ethanol, arsenite, and low doses of puromycin, but not by various other stresses. We found that increases in protein stability and translation efficiency both contribute to ZC3H11 accumulation. ZC3H11 is an in vitro substrate for casein kinase 1 isoform 2 (CK1.2), and results from CK1.2 depletion and other experiments suggest that phosphorylation of ZC3H11 can promote its instability in vivo. Results from sucrose density centrifugation indicate that under normal culture conditions translation initiation on the ZC3H11 mRNA is repressed, but after suitable stresses the ZC3H11 mRNA moves to heavy polysomes. The ZC3H11 3'-UTR is sufficient for translation suppression and a region of 71 nucleotides is required for the regulation. Since the control works in both bloodstream forms, where ZC3H11 translation is repressed at 37°C, and in procyclic forms, where ZC3H11 translation is activated at 37°C, we predict that this regulatory RNA sequence is targeted by repressive trans acting factor that is released upon stress.

Trypanosome MKT1 and the RNA-binding protein ZC3H11: interactions and potential roles in post-transcriptional regulatory networks.

The trypanosome zinc finger protein ZC3H11 binds to AU-rich elements in mRNAs. It is essential for survival of the mammalian-infective bloodstream form, where it stabilizes several mRNAs including some encoding chaperones, and is also required for stabilization of chaperone mRNAs during the heat-shock response in the vector-infective procyclic form. When ZC3H11 was artificially 'tethered' to a reporter mRNA in bloodstream forms it increased reporter expression. We here show that ZC3H11 interacts with trypanosome MKT1 and PBP1, and that domains required for both interactions are necessary for function in the bloodstream-form tethering assay. PBP1 interacts with MKT1, LSM12 and poly(A) binding protein, and localizes to granules during parasite starvation. All of these proteins are essential for bloodstream-form trypanosome survival and increase gene expression in the tethering assay. MKT1 is cytosolic and polysome associated. Using a yeast two-hybrid screen and tandem affinity purification we found that trypanosome MKT1 interacts with multiple RNA-binding proteins and other potential RNA regulators, placing it at the centre of a post-transcriptional regulatory network. A consensus interaction sequence, H(E/D/N/Q)PY, was identified. Recruitment of MKT1-containing regulatory complexes to mRNAs via sequence-specific mRNA-binding proteins could thus control several different post-transcriptional regulons.

Post-transcriptional regulation of the trypanosome heat shock response by a zinc finger protein.

In most organisms, the heat-shock response involves increased heat-shock gene transcription. In Kinetoplastid protists, however, virtually all control of gene expression is post-transcriptional. Correspondingly, Trypanosoma brucei heat-shock protein 70 (HSP70) synthesis after heat shock depends on regulation of HSP70 mRNA turnover. We here show that the T. brucei CCCH zinc finger protein ZC3H11 is a post-transcriptional regulator of trypanosome chaperone mRNAs. ZC3H11 is essential in bloodstream-form trypanosomes and for recovery of insect-form trypanosomes from heat shock. ZC3H11 binds to mRNAs encoding heat-shock protein homologues, with clear specificity for the subset of trypanosome chaperones that is required for protein refolding. In procyclic forms, ZC3H11 was required for stabilisation of target chaperone-encoding mRNAs after heat shock, and the HSP70 mRNA was also decreased upon ZC3H11 depletion in bloodstream forms. Many mRNAs bound to ZC3H11 have a consensus AUU repeat motif in the 3'-untranslated region. ZC3H11 bound preferentially to AUU repeats in vitro, and ZC3H11 regulation of HSP70 mRNA in bloodstream forms depended on its AUU repeat region. Tethering of ZC3H11 to a reporter mRNA increased reporter expression, showing that it is capable of actively stabilizing an mRNA. These results show that expression of trypanosome heat-shock genes is controlled by a specific RNA-protein interaction. They also show that heat-shock-induced chaperone expression in procyclic trypanosome enhances parasite survival at elevated temperatures.