Background: In spore-forming bacteria, the molecular mechanisms of accumulation of transfer RNA (tRNA) during sporulation must be a priority as tRNAs play an essential role in protein synthesis during spore germination and outgrowth. However, tRNA processing has not been extensively studied in these conditions, and knowledge of these mechanisms is important to understand long-term stress survival.    Methods:To gain further insight into tRNA processing during spore germination and outgrowth, the expression of the single copy tRNA Cys gene was analyzed in the presence and absence of 1.2 M NaCl in Bacillus subtilis using RNA-Seq data obtained from the Gene Expression Omnibus (GEO) database. The CLC Genomics work bench 12.0.2 (CLC Bio, Aarhus, Denmark, https://www.qiagenbioinformatics.com/) was used to analyze reads from the tRNA Cys gene.  Results:The results show that spores store different populations of tRNA Cys-related molecules.  One such population, representing 60% of total tRNA Cys, was composed of tRNA Cys fragments.  Half of these fragments (3´-tRF) possessed CC, CCA or incorrect additions at the 3´end. tRNA Cys with correct CCA addition at the 3´end represented 23% of total tRNA Cys, while with CC addition represented 9% of the total and with incorrect addition represented 7%. While an accumulation of tRNA Cys precursors was induced by upregulation of the rrnD operon under the control of  σ A -dependent promoters under both conditions investigated, salt stress produced only a modest effect on tRNA Cys expression and the accumulation of tRNA Cys related species. Conclusions:The results demonstrate that tRNA Cys molecules resident in spores undergo dynamic processing to produce functional molecules that may play an essential role during protein synthesis.

The wobble nucleoside 5-methylaminomethyl-2-thio-uridine (mnm5s2U) is present in bacterial tRNAs specific for Lys and Glu and 5-carboxymethylaminomethyl-2-thio-uridine (cmnm5s2U) in tRNA specific for Gln. The sulfur of (c)mnm5s2U may be exchanged by selenium (Se)-a reaction catalyzed by the selenophosphate-dependent tRNA 2-selenouridine synthase encoded by the mnmH (ybbB, selU, sufY) gene. The MnmH protein has a rhodanese domain containing one catalytic Cys (C97) and a P-loop domain containing a Walker A motif, which is a potential nucleotide binding site. We have earlier isolated a mutant of Salmonella enterica, serovar Typhimurium with an alteration in the rhodanese domain of the MnmH protein (G67E) mediating the formation of modified nucleosides having a geranyl (ge)-group (C10H17-fragment) attached to the s2 group of mnm5s2U and of cmnm5s2U in tRNA. To further characterize the structural requirements to increase the geranylation activity, we here report the analysis of 39 independently isolated mutants catalyzing the formation of mnm5ges2U. All these mutants have amino acid substitutions in the rhodanese domain demonstrating that this domain is pivotal to increase the geranylation activity. The wild type form of MnmH+ also possesses geranyltransferase activity in vitro although only a small amount of the geranyl derivatives of (c)mnm5s2U is detected in vivo. The selenation activity in vivo has an absolute requirement for the catalytic Cys97 in the rhodanese domain whereas the geranylation activity does not. Clearly, MnmH has two distinct enzymatic activities for which the rhodanese domain is pivotal. An intact Walker motif in the P-loop domain is required for the geranylation activity implying that it is the binding site for geranylpyrophosphate (GePP), which is the donor molecule in vitro in the geranyltransfer reaction. Purified MnmH from wild type and from the MnmH(G67E) mutant have bound tRNA, which is enriched with geranylated tRNA. This in conjunction with earlier published data, suggests that this bound geranylated tRNA may be an intermediate in the selenation of the tRNA.

Picornaviruses are non-enveloped RNA viruses that enter cells via receptor-mediated endocytosis. Because they lack an envelope, picornaviruses face the challenge of delivering their RNA genomes across the membrane of the endocytic vesicle into the cytoplasm to initiate infection. Currently, the mechanism of genome release and translocation across membranes remains poorly understood. Within the enterovirus genus, poliovirus, rhinovirus 2, and rhinovirus 16 have been proposed to release their genomes across intact endosomal membranes through virally induced pores, whereas one study has proposed that rhinovirus 14 releases its RNA following disruption of endosomal membranes. For the more distantly related aphthovirus genus (e.g. foot-and-mouth disease viruses and equine rhinitis A virus) acidification of endosomes results in the disassembly of the virion into pentamers and in the release of the viral RNA into the lumen of the endosome, but no details have been elucidated as how the RNA crosses the vesicle membrane. However, more recent studies suggest aphthovirus RNA is released from intact particles and the dissociation to pentamers may be a late event. In this study we have investigated the RNase A sensitivity of genome translocation of poliovirus using a receptor-decorated-liposome model and the sensitivity of infection of poliovirus and equine-rhinitis A virus to co-internalized RNase A. We show that poliovirus genome translocation is insensitive to RNase A and results in little or no release into the medium in the liposome model. We also show that infectivity is not reduced by co-internalized RNase A for poliovirus and equine rhinitis A virus. Additionally, we show that all poliovirus genomes that are internalized into cells, not just those resulting in infection, are protected from RNase A. These results support a finely coordinated, directional model of viral RNA delivery that involves viral proteins and cellular membranes.

In plant, non-specific lipid transfer proteins (nsLTPs) are small, basic proteins that have been reported to be involved in numerous biological processes such as transfer of phospholipids, reproductive development, pathogen defence and abiotic stress response. To date, only a tiny fraction of plant nsLTPs have been functionally identified, and even fewer have been identified in maize [Zea mays (Zm)].

Symbiotic mutualisms of bacteria and animals are ubiquitous in nature, running a continuum from facultative to obligate from the perspectives of both partners. The loss of functions required for living independently but not within a host gives rise to reduced genomes in many symbionts. Although the phenomenon of genome reduction can be explained by existing evolutionary models, the initiation of the process is not well understood. Here, we describe the microbiome associated with the eggs of the beetle Lagria villosa, consisting of multiple bacterial symbionts related to Burkholderia gladioli, including a reduced-genome symbiont thought to be the exclusive producer of the defensive compound lagriamide. We show that the putative lagriamide-producing symbiont is the only member of the microbiome undergoing genome reduction and that it has already lost the majority of its primary metabolism and DNA repair pathways. The key step preceding genome reduction in the symbiont was likely the horizontal acquisition of the putative lagriamide lga biosynthetic gene cluster. Unexpectedly, we uncovered evidence of additional horizontal transfers to the symbiont's genome while genome reduction was occurring and despite a current lack of genes needed for homologous recombination. These gene gains may have given the genome-reduced symbiont a selective advantage in the microbiome, especially given the maintenance of the large lga gene cluster despite ongoing genome reduction.IMPORTANCE Associations between microorganisms and an animal, plant, or fungal host can result in increased dependence over time. This process is due partly to the bacterium not needing to produce nutrients that the host provides, leading to loss of genes that it would need to live independently and to a consequent reduction in genome size. It is often thought that genome reduction is aided by genetic isolation-bacteria that live in monocultures in special host organs, or inside host cells, have less access to other bacterial species from which they can obtain genes. Here, we describe exposure of a genome-reduced beetle symbiont to a community of related bacteria with nonreduced genomes. We show that the symbiont has acquired genes from other bacteria despite going through genome reduction, suggesting that isolation has not yet played a major role in this case of genome reduction, with horizontal gene gains still offering a potential route for adaptation.

Fluorescence Resonance Energy Transfer (FRET) is a well-known methodology for detection and quantitation of structural changes of proteins in solution. FRET requires site-specific protein labeling with two fluorophores, one of which functions as an energy donor and the other one as an energy acceptor. However, the site-specific labeling of protein is often complex and difficult, particularly when inserting two fluorophores in specific sites. We have examined several protein labeling approaches with a varying degree of success. Described here is a dual labeling strategy that worked reproducibly in a number of protein targets and we believe will be applicable to a variety of proteins, which have few or no native cysteine (Cys) residues. We have successfully double-labeled DnaA protein of Bacillus anthracis, which lacks intrinsic Cys residues. A cysteine residue was inserted at the N-terminus by in vitro mutagenesis and a Cys-Cys-Phe-Gly-Cys-Cys (CCPGCC) sequence at the C-terminus by PCR. This protein was labeled site-specifically with a fluorescein derivative, FlAsH, at the CCPGCC sequence followed by Alexa568 maleimide at the N-terminus Cys residue. Structural changes of the protein with nucleotide, DNA and an inhibitor protein binding were determined by FRET analysis of the double-labeled protein. This comprehensive novel methodology for site-specific protein labeling with different fluorophores is applicable for understanding different in vitro proteomic structural studies. Here, we describe a verified technique used for FRET spectral analysis and quantitative evaluation of structural changes using fluorophore labeled DnaA protein constructs as an example.

Recent efforts towards the comprehensive identification of RNA-bound proteomes have revealed a large, surprisingly diverse family of candidate RNA-binding proteins (RBPs). Quantitative metrics for characterization and validation of protein-RNA interactions and their dynamic interactions have, however, proven analytically challenging and prone to error. Here we report a method termed LEAP-RBP (Liquid-Emulsion-Assisted-Purification of RNA-Bound Protein) for the selective, quantitative recovery of UV-crosslinked RNA-protein complexes. By virtue of its high specificity and yield, LEAP-RBP distinguishes RNA-bound and RNA-free protein levels and reveals common sources of experimental noise in RNA-centric RBP enrichment methods. We introduce strategies for accurate RBP identification and signal-based metrics for quantifying protein-RNA complex enrichment, relative RNA occupancy, and method specificity. In this work, the utility of our approach is validated by comprehensive identification of RBPs whose association with mRNA is modulated in response to global mRNA translation state changes and through in-depth benchmark comparisons with current methodologies.

Neural stem cells (NSCs) in the adult murine subependymal zone balance their self-renewal capacity and glial identity with the potential to generate neurons during the lifetime. Adult NSCs exhibit lineage priming via pro-neurogenic fate determinants. However, the protein levels of the neural fate determinants are not sufficient to drive direct differentiation of adult NSCs, which raises the question of how cells along the neurogenic lineage avoid different conflicting fate choices, such as self-renewal and differentiation. Here, we identify RNA-binding protein MEX3A as a post-transcriptional regulator of a set of stemness associated transcripts at critical transitions in the subependymal neurogenic lineage. MEX3A regulates a quiescence-related RNA signature in activated NSCs that is needed for their return to quiescence, playing a role in the long-term maintenance of the NSC pool. Furthermore, it is required for the repression of the same program at the onset of neuronal differentiation. Our data indicate that MEX3A is a pivotal regulator of adult murine neurogenesis acting as a translational remodeller.

MicroRNAs (miRNAs) are small non-coding RNA molecules acting as post-transcriptional regulators of gene expression. They are involved in many biological processes, and their dysregulation is implicated in various diseases, including multiple sclerosis (MS). Interferon-beta (IFN-beta) is widely used as a first-line immunomodulatory treatment of MS patients. Here, we present the first longitudinal study on the miRNA expression changes in response to IFN-beta therapy. Peripheral blood mononuclear cells (PBMC) were obtained before treatment initiation as well as after two days, four days, and one month, from patients with clinically isolated syndrome (CIS) and patients with relapsing-remitting MS (RRMS). We measured the expression of 651 mature miRNAs and about 19,000 mRNAs in parallel using real-time PCR arrays and Affymetrix microarrays. We observed that the up-regulation of IFN-beta-responsive genes is accompanied by a down-regulation of several miRNAs, including members of the mir-29 family. These differentially expressed miRNAs were found to be associated with apoptotic processes and IFN feedback loops. A network of miRNA-mRNA target interactions was constructed by integrating the information from different databases. Our results suggest that miRNA-mediated regulation plays an important role in the mechanisms of action of IFN-beta, not only in the treatment of MS but also in normal immune responses. miRNA expression levels in the blood may serve as a biomarker of the biological effects of IFN-beta therapy that may predict individual disease activity and progression.

Several reports describing PIWI-interacting RNAs (piRNAs) in human cancer cells or in the bloodstream are affected by the presence of false positives in piRNA databases. A recent report suggested that piR-36249 regulates testicular cancer progression by engaging with DHX36 to regulate OAS2. However, piR-36249 is a tRNA-Cys 5' half capable of forming intermolecular G-quadruplexes. It is therefore expected that DHX36, a helicase with high affinity for DNA and RNA G-quadruplexes, was pulled down using piR-36249 mimicking probes. The suggestion of using piR-36249 as a therapeutic target for testicular cancer is therefore questionable, due to the consequences that tRNA inhibition could have on healthy cells.

The σ subunit of bacterial RNA polymerase (RNAP) controls recognition of the -10 and -35 promoter elements during transcription initiation. Free σ adopts a "closed," or inactive, conformation incompatible with promoter binding. The conventional two-state model of σ activation proposes that binding to core RNAP induces formation of an "open," active, σ conformation, which is optimal for promoter recognition. Using single-molecule Förster resonance energy transfer, we demonstrate that vegetative-type σ subunits exist in open and closed states even after binding to the RNAP core. As an extreme case, RNAP from Mycobacterium tuberculosis preferentially retains σ in the closed conformation, which is converted to the open conformation only upon binding by the activator protein RbpA and interaction with promoter DNA. These findings reveal that the conformational dynamics of the σ subunit in the RNAP holoenzyme is a target for regulation by transcription factors and plays a critical role in promoter recognition.

Transfer RNA (tRNA) is an adaptor molecule indispensable for assigning amino acids to codons on mRNA during protein synthesis. 2-thiouridine (s2U) derivatives in the anticodons (position 34) of tRNAs for glutamate, glutamine, and lysine are post-transcriptional modifications essential for precise and efficient codon recognition in all organisms. s2U34 is introduced either by (i) bacterial MnmA/eukaryote mitochondrial Mtu1 or (ii) eukaryote cytosolic Ncs6/archaeal NcsA, and the latter enzymes possess iron-sulfur (Fe-S) cluster. Here, we report the identification of novel-type MnmA homologs containing three conserved Cys residues, which could support Fe-S cluster binding and catalysis, in a broad range of bacteria, including thermophiles, Cyanobacteria, Mycobacteria, Actinomyces, Clostridium, and Helicobacter Using EPR spectroscopy, we revealed that Thermus thermophilus MnmA (TtMnmA) contains an oxygen-sensitive [4Fe-4S]-type cluster. Efficient in vitro formation of s2U34 in tRNALys and tRNAGln by holo-TtMnmA occurred only under anaerobic conditions. Mutational analysis of TtMnmA suggested that the Fe-S cluster is coordinated by the three conserved Cys residues (Cys105, Cys108, and Cys200), and is essential for its activity. Evolutionary scenarios for the sulfurtransferases, including the Fe-S cluster containing Ncs6/NcsA s2U thiouridylases and several distantly related sulfurtransferases, are proposed.

Accurate nutrient sensing is important for rapid fungal growth and exploitation of available resources. Sulfur is an important nutrient source found in a number of biological macromolecules, including proteins and lipids. The model filamentous fungus Neurospora crassa is capable of utilizing sulfur found in a variety of sources from amino acids to sulfate. During sulfur starvation, the transcription factor CYS-3 is responsible for upregulation of genes involved in sulfur uptake and assimilation. Using a combination of RNA sequencing and DNA affinity purification sequencing, we performed a global survey of the N. crassa sulfur starvation response and the role of CYS-3 in regulating sulfur-responsive genes. The CYS-3 transcription factor bound the promoters and regulated genes involved in sulfur metabolism. Additionally, CYS-3 directly activated the expression of a number of uncharacterized transporter genes, suggesting that regulation of sulfur import is an important aspect of regulation by CYS-3. CYS-3 also directly regulated the expression of genes involved in mitochondrial electron transfer. During sulfur starvation, genes involved in nitrogen metabolism, such as amino acid and nucleic acid metabolic pathways, along with genes encoding proteases and nucleases that are necessary for scavenging nitrogen, were activated. Sulfur starvation also caused changes in the expression of genes involved in carbohydrate metabolism, such as those encoding glycosyl hydrolases. Thus, our data suggest a connection between sulfur metabolism and other aspects of cellular metabolism. IMPORTANCE Identification of nutrients present in the environment is a challenge common to all organisms. Sulfur is an important nutrient source found in proteins, lipids, and electron carriers that are required for the survival of filamentous fungi such as Neurospora crassa. Here, we transcriptionally profiled the response of N. crassa to characterize the global response to sulfur starvation. We also used DNA affinity purification sequencing to identify the direct downstream targets of the transcription factor responsible for regulating genes involved in sulfur uptake and assimilation. Along with genes involved in sulfur metabolism, this transcription factor regulated a number of uncharacterized transporter genes and genes involved in mitochondrial electron transfer. Our data also suggest a connection between sulfur, nitrogen, and carbon metabolism, indicating that the regulation of a number of metabolic pathways is intertwined.

Interfacing gene delivery vehicles with biomaterials has the potential to play a key role in diversifying gene transfer capabilities, including localized, patterned, and controlled delivery. However, strategies for modifying biomaterials to interact with delivery vectors must be redesigned whenever new delivery vehicles and applications are explored. We have developed a vector-independent biomaterial platform capable of interacting with various adeno-associated viral (AAV) serotypes. A water-soluble, cysteine-tagged, recombinant protein version of the recently discovered multi-AAV serotype receptor (AAVR), referred to as cys-AAVR, was conjugated to maleimide-displaying polycaprolactone (PCL) materials using click chemistry. The resulting cys-AAVR-PCL system bound to a broad range of therapeutically relevant AAV serotypes, thereby providing a platform capable of modulating the delivery of all AAV serotypes. Intramuscular injection of cys-AAVR-PCL microspheres with bound AAV vectors resulted in localized and sustained gene delivery as well as reduced spread to off-target organs compared to a vector solution. This cys-AAVR-PCL system is thus an effective approach for biomaterial-based AAV gene delivery for a broad range of therapeutic applications.

The flavoenzyme glutathione reductase (GR; NADPH+glutathione disulphide+H(+)-->NADP(+)+2 glutathione-SH) of Plasmodium falciparum is a promising drug target against tropical malaria. As P. falciparum genes are assumed to be highly polymorphic we have cloned and expressed the GR cDNA of the chloroquine-sensitive strain 3D7. In comparison to the known GR of the chloroquine-resistant K1 strain there are three base exchanges all of them leading to amino acid substitutions (residues 281, 285 and 335). The catalytic efficiency k(cat)/K(m) of the 3D7 enzyme is 5-fold lower than for the K1 enzyme. In contrast, vis-à-vis the drugs carmustine, methylene blue and fluorophenyliso-alloxazine the two enzyme species exhibited identical inhibition kinetics. Two structural motifs which are specific for P. falciparum GR were studied by mutational deletion analysis of 3D7 GR. Loop 126-138 appears to be important for folding and stability of the enzyme, whereas the subdomain 318-350 was found to be involved in FAD-binding. The subdomain has no major influence on the known functions of the catalytic triad Cys-40, Cys-45 and His-485'. Flavin absorption spectroscopy of inactive point mutants showed that Cys-45 forms a thiolate charge transfer complex and Cys-40 is the interchange thiol, which reduces glutathione disulphide. The mutant His-485-->Gln had a normal K(m) for glutathione disulphide reduction but only 0.8% residual catalytic activity when compared with wild-type GR, which confirms its function as an acid/base catalyst. The parasite-specific domains in combination with the reactive catalytic residues appear to be a suitable target matrix for inhibiting GR in vivo.

Nuclear group I introns are restricted to the ribosomal DNA locus where they interrupt genes for small subunit and large subunit ribosomal RNAs at conserved sites in some eukaryotic microorganisms. Here, the myxomycete protists are a frequent source of nuclear group I introns due to their unique life strategy and a billion years of separate evolution. The ribosomal DNA of the myxomycete Mucilago crustacea was investigated and found to contain seven group I introns, including a direct repeat-containing intron at insertion site S1389 in the small subunit ribosomal RNA gene. We collected, analyzed, and compared 72 S1389 group IC1 introns representing diverse myxomycete taxa. The consensus secondary structure revealed a conserved ribozyme core, but with surprising sequence variations in the guanosine binding site in segment P7. Some S1389 introns harbored large extension sequences in the peripheral region of segment P9 containing direct repeat arrays. These repeats contained up to 52 copies of a putative internal guide sequence motif. Other S1389 introns harbored homing endonuclease genes in segment P1 encoding His-Cys proteins. Homing endonuclease genes were further interrupted by small spliceosomal introns that have to be removed in order to generate the open reading frames. Phylogenetic analyses of S1389 intron and host gene indicated both vertical and horizontal intron transfer during evolution, and revealed sporadic appearances of direct repeats, homing endonuclease genes, and guanosine binding site variants among the myxomycete taxa.

The reactive cysteine residue at position 106 (Cys106) of DJ-1 is preferentially oxidized under oxidative stress, generating oxidized DJ-1 (oxDJ-1). Oxidation of Cys106 to sulfinic acid changes the biologic action of DJ-1 and increases its cytoprotective properties. The similar activation step is known in peroxiredoxins (Prxs), in which oxidation of reactive Cys to sulfinic acid induces polymerization of Prxs and changes its enzyme characteristic from peroxidase to molecular chaperone. In the present study, oxDJ-1 was prepared and its polymerization and related amino acid residues were investigated. We found that oxDJ-1 formed a characteristic polymer with disulfide bonds and with noncovalent and covalent binding other than disulfide. The physiological concentration of glutathione resolved the polymer form of oxDJ-1, and glutathionylation of other two Cys residues, such as Cys 46 and 53, was detected. Mutant analysis indicated the necessity not only of Cys106 but also of Cys46 for the polymer formation. The cellular experiment demonstrated that the electrophilic quinone treatment induced a high-molecular-weight complex containing oxDJ-1. Dynamic polymerization of oxDJ-1 with a ring and a stacked structure was observed by an atomic force microscope. Collectively, these results clearly demonstrated the characteristic polymer formation of oxDJ-1 with a disulfide bond and noncovalent and covalent binding other than disulfide, which might be related to the biologic function of oxDJ-1.

Thioredoxin systems, composed of thioredoxin reductase (TrxR), thioredoxin (Trx) and NADPH, play important roles in maintaining cellular redox homeostasis and redox signaling. Recently the cytosolic Trx1 system has been shown to be a cellular target of arsenic containing compounds. To elucidate the relationship of the structure of arsenic compounds with their ability of inhibiting TrxR1 and Trx1, and cytotoxicity, we have investigated the reaction of Trx1 system with seven arsenic trithiolates: As(Cys)3, As(GS)3, As(Penicillamine)3, As(Mercaptoethanesulfonate)3, As(Mercaptopurine)3, As(2-mercaptopyridine)3 and As(2-mercaptopyridine N-oxide)3. The cytotoxicity of these arsenicals was consistent with their ability to inhibit TrxR1 in vitro and in cells. Unlike other arsenicals, As(Mercaptopurine)3 which did not show inhibitory effects on TrxR1 had very weak cytotoxicity, indicating that TrxR1 is a reliable drug target for arsenicals. Moreover, the two aromatic compounds As(2-mercaptopyridine)3 and As(2-mercaptopyridine N-oxide)3 showed stronger cytotoxicity than the others. As(2-mercaptopyridine)3 which selectively oxidized two structural cysteines (Cys62 and Cys69) in Trx1 showed mild improvement in cytotoxicity. As(2-mercaptopyridine N-oxide)3 oxidized all the Cys residues in Trx1, exhibiting the strongest cytotoxicity. Oxidation of Trx1 by As(2-mercaptopyridine)3 and As(2-mercaptopyridine N-oxide)3 affected electron transfer from NADPH and TrxR1 to peroxiredoxin 1 (Prx1), which could result in the reactive oxygen species elevation and trigger cell death process. These results suggest that oxidation of structural cysteine residues in Trx1 by aromatic group in TrxR1-targeting drugs may sensitize tumor cells to cell death, providing a novel approach to regulate cellular redox signaling and also a basis for rational design of new anticancer agents.

Non-native conformations drive protein misfolding diseases, complicate bioengineering efforts, and fuel molecular evolution. No current experimental technique is well-suited for elucidating them and their phenotypic effects. Especially intractable are the transient conformations populated by intrinsically disordered proteins. We describe an approach to systematically discover, stabilize, and purify native and non-native conformations, generated in vitro or in vivo, and directly link conformations to molecular, organismal, or evolutionary phenotypes. This approach involves high-throughput disulfide scanning (HTDS) of the entire protein. To reveal which disulfides trap which chromatographically resolvable conformers, we devised a deep-sequencing method for double-Cys variant libraries of proteins that precisely and simultaneously locates both Cys residues within each polypeptide. HTDS of the abundant E. coli periplasmic chaperone HdeA revealed distinct classes of disordered hydrophobic conformers with variable cytotoxicity depending on where the backbone was cross-linked. HTDS can bridge conformational and phenotypic landscapes for many proteins that function in disulfide-permissive environments.

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.