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.

The HTLV-1 viral oncoproteins Tax and HBZ reprogram the cellular mRNA splicing landscape.

Viral infections are known to hijack the transcription and translation of the host cell. However, the extent to which viral proteins coordinate these perturbations remains unclear. Here we used a model system, the human T-cell leukemia virus type 1 (HTLV-1), and systematically analyzed the transcriptome and interactome of key effectors oncoviral proteins Tax and HBZ. We showed that Tax and HBZ target distinct but also common transcription factors. Unexpectedly, we also uncovered a large set of interactions with RNA-binding proteins, including the U2 auxiliary factor large subunit (U2AF2), a key cellular regulator of pre-mRNA splicing. We discovered that Tax and HBZ perturb the splicing landscape by altering cassette exons in opposing manners, with Tax inducing exon inclusion while HBZ induces exon exclusion. Among Tax- and HBZ-dependent splicing changes, we identify events that are also altered in Adult T cell leukemia/lymphoma (ATLL) samples from two independent patient cohorts, and in well-known cancer census genes. Our interactome mapping approach, applicable to other viral oncogenes, has identified spliceosome perturbation as a novel mechanism coordinated by Tax and HBZ to reprogram the transcriptome.

Kinetic and thermodynamic characterisation of HIV-protease inhibitors against E35D↑G↑S mutant in the South African HIV-1 subtype C protease.

Herein, we report the effect of nine FDA approved protease inhibitor drugs against a new HIV-1 subtype C mutant protease, E35D↑G↑S. The mutant has five mutations, E35D, two insertions, position 36 (G and S), and D60E. Kinetics, inhibition constants, vitality, Gibbs free binding energies are reported. The variant showed a decreased affinity for substrate and low catalytic efficiency compared to the wild type. There was a significant decrease in the binding of seven FDA approved protease inhibitors against the mutant (p < .0001). Amprenavir and ritonavir showed the least decrease, but still significant reduced activity in comparison to the wildtype (4 and 5 folds, respectively, p = .0021 and .003, respectively). Nelfinavir and atazanavir were the worst inhibitors against the variant as seen from the IC50, with values of 1401 ± 3.0 and 685 ± 3.0 nM, respectively. Thermodynamics data showed less favourable Gibbs free binding energies for the protease inhibitors to the mutant.

An insight to the molecular interactions of the FDA approved HIV PR drugs against L38L↑N↑L PR mutant.

The aspartate protease of the human immune deficiency type-1 virus (HIV-1) has become a crucial antiviral target in which many useful antiretroviral inhibitors have been developed. However, it seems the emergence of new HIV-1 PR mutations enhances drug resistance, hence, the available FDA approved drugs show less activity towards the protease. A mutation and insertion designated L38L↑N↑L PR was recently reported from subtype of C-SA HIV-1. An integrated two-layered ONIOM (QM:MM) method was employed in this study to examine the binding affinities of the nine HIV PR inhibitors against this mutant. The computed binding free energies as well as experimental data revealed a reduced inhibitory activity towards the L38L↑N↑L PR in comparison with subtype C-SA HIV-1 PR. This observation suggests that the insertion and mutations significantly affect the binding affinities or characteristics of the HIV PIs and/or parent PR. The same trend for the computational binding free energies was observed for eight of the nine inhibitors with respect to the experimental binding free energies. The outcome of this study shows that ONIOM method can be used as a reliable computational approach to rationalize lead compounds against specific targets. The nature of the intermolecular interactions in terms of the host-guest hydrogen bond interactions is discussed using the atoms in molecules (AIM) analysis. Natural bond orbital analysis was also used to determine the extent of charge transfer between the QM region of the L38L↑N↑L PR enzyme and FDA approved drugs. AIM analysis showed that the interaction between the QM region of the L38L↑N↑L PR and FDA approved drugs are electrostatic dominant, the bond stability computed from the NBO analysis supports the results from the AIM application. Future studies will focus on the improvement of the computational model by considering explicit water molecules in the active pocket. We believe that this approach has the potential to provide information that will aid in the design of much improved HIV-1 PR antiviral drugs.

I36T↑T mutation in South African subtype C (C-SA) HIV-1 protease significantly alters protease-drug interactions.

The efficacy of HIV-1 protease (PR) inhibition therapies is often compromised by the emergence of mutations in the PR molecule that reduces the binding affinity of inhibitors while maintaining viable catalytic activity and affinity for natural substrates. In the present study, we used a recombinant HIV-1 C-SA PR and a recently reported variant for inhibition (Ki, IC50) and thermodynamic studies against nine clinically used inhibitors. This is the first time that binding free energies for C-SA PR and the mutant are reported. This variant PR harbours a mutation and insertion (I36T↑T) at position 36 of the C-SA HIV-1 PR, and did not show a significant difference in the catalytic effect of the HIV-1 PR. However, the nine clinically approved HIV PR drugs used in this study demonstrated weaker inhibition and lower binding affinities toward the variant when compared to the wild type HIV-1 PR. All the protease inhibitors (PIs), except Amprenavir and Ritonavir exhibited a significant decrease in binding affinity (p<0.0001). Darunavir and Nelfinavir exhibited the weakest binding affinity, 155- and 95-fold decreases respectively, toward the variant. Vitality values for the variant PR, against the seven selected PIs, confirm the impact of the mutation and insertion on the South African HIV-1 subtype C PR. This information has important clinical implications for thousands of patients in Sub-Saharan Africa.

Purification and characterization of naturally occurring HIV-1 (South African subtype C) protease mutants from inclusion bodies.

Human immunodeficiency virus (HIV) infections in sub-Saharan Africa represent about 56% of global infections. Many studies have targeted HIV-1 protease for the development of drugs against AIDS. Recombinant HIV-1 protease is used to screen new drugs from synthetic compounds or natural substances. Along with the wild type (C-SA) we also over-expressed and characterized two mutant forms from patients that had shown resistance to protease inhibitors. Using recombinant DNA technology, we constructed three recombinant plasmids in pGEX-6P-1 and expressed them containing a sequence encoding wild type HIV protease and two mutants (I36T↑T contains 100 amino acids and L38L↑N↑L contains 101 amino acids). These recombinant proteins were isolated from inclusion bodies by using QFF anion exchange and GST trap columns. In SDS-PAGE, we obtained these HIV proteases as single bands of approximately 11.5, 11.6 and 11.7 kDa for the wild type, I36T↑Tand L38L↑N↑L mutants, respectively. The enzyme was recovered efficiently (0.25 mg protein/L of Escherichia coli culture) and had high specific activity of 2.02, 2.20 and 1.33 µmol min(-1) mg(-1) at an optimal pH of 5 and temperature of 37 °C for the wild type, I36T↑T and L38L↑N↑L, respectively. The method employed here provides an easy and rapid purification of the HIV-1(C-SA) protease from the inclusion bodies, with high yield and high specific activities.

Identification of small molecule antivirals against HTLV-1 by targeting the hDLG1-Tax-1 protein-protein interaction.

Human T-cell leukemia virus type-1 (HTLV-1) is the first pathogenic retrovirus discovered in human. Although HTLV-1-induced diseases are well-characterized and linked to the encoded Tax-1 oncoprotein, there is currently no strategy to target Tax-1 functions with small molecules. Here, we analyzed the binding of Tax-1 to the human homolog of the drosophila discs large tumor suppressor (hDLG1/SAP97), a multi-domain scaffolding protein involved in Tax-1-transformation ability. We have solved the structures of the PDZ binding motif (PBM) of Tax-1 in complex with the PDZ1 and PDZ2 domains of hDLG1 and assessed the binding of 10 million molecules by virtual screening. Among the 19 experimentally confirmed compounds, one systematically inhibited the Tax-1-hDLG1 interaction in different biophysical and cellular assays, as well as HTLV-1 cell-to-cell transmission in a T-cell model. Thus, our work demonstrates that interactions involving Tax-1 PDZ-domains are amenable to small-molecule inhibition, which provides a framework for the design of targeted therapies for HTLV-1-induced diseases.

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.

A proteome-scale map of the SARS-CoV-2-human contactome.

Understanding the mechanisms of coronavirus disease 2019 (COVID-19) disease severity to efficiently design therapies for emerging virus variants remains an urgent challenge of the ongoing pandemic. Infection and immune reactions are mediated by direct contacts between viral molecules and the host proteome, and the vast majority of these virus-host contacts (the 'contactome') have not been identified. Here, we present a systematic contactome map of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) with the human host encompassing more than 200 binary virus-host and intraviral protein-protein interactions. We find that host proteins genetically associated with comorbidities of severe illness and long COVID are enriched in SARS-CoV-2 targeted network communities. Evaluating contactome-derived hypotheses, we demonstrate that viral NSP14 activates nuclear factor κB (NF-κB)-dependent transcription, even in the presence of cytokine signaling. Moreover, for several tested host proteins, genetic knock-down substantially reduces viral replication. Additionally, we show for USP25 that this effect is phenocopied by the small-molecule inhibitor AZ1. Our results connect viral proteins to human genetic architecture for COVID-19 severity and offer potential therapeutic targets.

Mechanistic insight on the inhibition of D, D-carboxypeptidase from Mycobacterium tuberculosis by ß-lactam antibiotics: an ONIOM acylation study.

Mycobacterium tuberculosis cell wall is intricate and impermeable to many agents. A D, D-carboxypeptidase (DacB1) is one of the enzymes involved in the biosynthesis of cell wall peptidoglycan and catalyzes the terminal D-alanine cleavage from pentapeptide precursors. Catalytic activity and mechanism by which DacB1 functions is poorly understood. Herein, we investigated the acylation mechanism of DacB1 by ß-lactams using a 6-membered ring transition state model that involves a catalytic water molecule in the reaction pathway. The full transition states (TS) optimization plus frequency were achieved using the ONIOM (B3LYP/6-31 + G(d): AMBER) method. Subsequently, the activation free energies were computed via single-point calculations on fully optimized structures using B3LYP/6-311++(d,p): AMBER and M06-2X/6-311++(d,p): AMBER with an electronic embedding scheme. The 6-membered ring transition state is an effective model to examine the inactivation of DacB1 via acylation by ß-lactams antibiotics (imipenem, meropenem, and faropenem) in the presence of the catalytic water. The ΔG# values obtained suggest that the nucleophilic attack on the carbonyl carbon is the rate-limiting step with 13.62, 19.60 and 30.29 kcal mol-1 for Imi-DacB1, Mero-DacB1 and Faro-DacB1, respectively. The electrostatic potential (ESP) and natural bond orbital (NBO) analysis provided significant electronic details of the electron-rich region and charge delocalization, respectively, based on the concerted 6-membered ring transition state. The stabilization energies of charge transfer within the catalytic reaction pathway concurred with the obtained activation free energies. The outcomes of this study provide important molecular insight into the inactivation of D, D-carboxypeptidase by ß-lactams.Communicated by Ramaswamy H. Sarma.

Optimized Procedure for Recovering HIV-1 Protease (C-SA) from Inclusion Bodies.

HIV-1 is an infectious virus that causes acquired immunodeficiency syndrome (AIDS) and it is one of the major causes of deaths worldwide. The production of HIV-1 protease (PR) on a large scale has been a problem for scientists due to its cytotoxicity, low yield, insolubility, and low activity. HIV-1 C-SA protease has been cloned, expressed, and purified previously, however, with low recovery (0.25 mg/L). Herein we report an optimal expression and solubilisation procedure to recover active HIV-1 C-SA protease enzyme from inclusion bodies. The HIV protease was expressed in seven different vectors (pET11b, pET15b, pET28a pET32a, pET39b, pET41b and pGEX 6P-1). The highest expression was achieved when the vector pET32a (Trx tag) was employed. A total of 19.5 mg of fusion protein was refolded of which 5.5 mg of active protease was obtained after cleavage. The free protease had a high specific activity of 2.81 µmoles/min/mg. Interestingly the Trx-fusion protein also showed activity closer (1.24 µmoles/min/mg) to that of the free protease suggesting that the pET32a vector (Trx tag) expressed in BL21(DE3) pLysS provides a more efficient way to obtain HIV-1 protease.

Identification of potent L,D-transpeptidase 5 inhibitors for Mycobacterium tuberculosis as potential anti-TB leads: virtual screening and molecular dynamics simulations.

Virtual screening is a useful in silico approach to identify potential leads against various targets. It is known that carbapenems (doripenem and faropenem) do not show any reasonable inhibitory activities against L,D-transpeptidase 5 (LdtMt5) and also an adduct of meropenem exhibited slow acylation. Since these drugs are active against L,D-transpeptidase 2 (LdtMt2), understanding the differences between these two enzymes is essential. In this study, a ligand-based virtual screening of 12,766 compounds followed by molecular dynamics (MD) simulations was applied to identify potential leads against LdtMt5. To further validate the obtained virtual screening ranking for LdtMt5, we screened the same libraries of compounds against LdtMt2 which had more experimetal and calculated binding energies reported. The observed consistency between the binding affinities of LdtMt2 validates the obtained virtual screening binding scores for LdtMt5. We subjected 37 compounds with docking scores ranging from - 7.2 to - 9.9 kcal mol-1 obtained from virtual screening for further MD analysis. A set of compounds (n = 12) from four antibiotic classes with ≤ - 30 kcal mol-1 molecular mechanics/generalized born surface area (MM-GBSA) binding free energies (ΔGbind) was characterized. A final set of that, all ß-lactams (n = 4), was considered. The outcome of this study provides insight into the design of potential novel leads for LdtMt5. Graphical abstract.

Exploring the flap dynamics of the South African HIV subtype C protease in presence of FDA-approved inhibitors: MD study.

HIV-1 protease (HIV PR) is considered as one of the most attractive targets for the treatment of HIV and the impact of flap dynamics of HIV PR on the binding affinities of protease inhibitors (PIs) is a crucial ongoing research field. Recently, our research group evaluated the binding affinities of different FDA approved PIs against the South African HIV-1 subtype C (C-SA) protease (PR). The CSA-HIV PR displayed weaker binding affinity for most of the clinical PIs compared to HIV-1 B subtype for West and Central Europe, the Americas. In the current work, the flap dynamics of four different systems of HIV-1 C-SA PR complexed to FDA approved second generation PIs and its impact on binding was explored over the molecular dynamic trajectories. It was observed that the interactions of the selected drugs with the binding site residues of the protease may not be the major contributor for affinity towards PIs. Various post-MD analyses were performed, also entropic contributions, solvation free energies and hydrophobic core formation interactions were studied to assess how the flap dynamics of C-SA PR which is affected by such factors. From these contributions, large van der Waals interactions and low solvation free energies were found to be major factors for the higher activity of ATV against C-SA HIV PR. Furthermore, a comparatively stable hydrophobic core may be responsible for higher stability of the PR flaps of the ATV complex. The outcome of this study provides significant guidance to how the flap dynamics of C-SA PR is affected by various factors as a result of the binding affinity of various protease inhibitors. It will also assist with the design of potent inhibitors against C-SA HIV PR that apart from binding in the active site of PR can interacts with the flaps to prevent opening of the flaps resulting in inactivation of the protease.