AI is a viable alternative to high throughput screening: a 318-target study

High throughput screening (HTS) is routinely used to identify bioactive small molecules. This requires physical compounds, which limits coverage of accessible chemical space. Computational approaches combined with vast on-demand chemical libraries can access far greater chemical space, provided that the predictive accuracy is sufficient to identify useful molecules. Through the largest and most diverse virtual HTS campaign reported to date, comprising 318 individual projects, we demonstrate that our AtomNet® convolutional neural network successfully finds novel hits across every major therapeutic area and protein class. We address historical limitations of computational screening by demonstrating success for target proteins without known binders, high-quality X-ray crystal structures, or manual cherry-picking of compounds. We show that the molecules selected by the AtomNet® model are novel drug-like scaffolds rather than minor modifications to known bioactive compounds. Our empirical results suggest that computational methods can substantially replace HTS as the first step of small-molecule drug discovery.

Case report: Homozygous variants of NEB and KLHL40 in two Arab patients with nemaline myopathy

Objective: Nemaline myopathies are a heterogeneous group of congenital myopathies caused by mutations in different genes associated with the structural and functional proteins of thin muscular filaments. Most patients have congenital onset characterized by hypotonia, respiratory issues, and abnormal deep tendon reflexes, which is a phenotype encountered in a wide spectrum of neuromuscular disorders. Whole-exome sequencing (WES) contributes to a faster diagnosis and facilitates genetic counseling. Methods: Here, we report on two Arab patients from consanguineous families diagnosed with nemaline myopathy of different phenotype spectrum severities. Results: Clinical assessment and particular prenatal history raised suspicion of neuromuscular disease. WES identified homozygous variants in NEB and KLHL40. Muscle biopsy and muscle magnetic resonance imaging studies linked the genetic testing results to the clinical phenotype. The novel variant in the NEB gene resulted in a classical type 2 nemaline myopathy, while the KLHL40 gene variant led to a severe phenotype of nemaline myopathy, type 8. Both patients were identified as having other gene variants with uncertain roles in their complex phenotypes. Conclusions: This study enriches the phenotypic spectrum of nemaline myopathy caused by NEB and KLHL40 variants and highlights the importance of detailed prenatal, neonatal, and infancy assessments of muscular weakness associated with complex systemic features. Variants of uncertain significance in genes associated with nemaline myopathy may be correlated with the phenotype. Early, multidisciplinary intervention can improve the outcome in patients with mild forms of nemaline myopathies. WES is essential for clarifying complex clinical phenotypes encountered in patients from consanguineous families. Targeted carrier screening of extended family members would enable accurate genetic counseling and potential genetic prevention.

How concerning is a SARS-CoV-2 variant of concern? Computational predictions and the variants labeling system

In this study, we evaluated the use of a predictive computational approach for SARS-CoV-2 genetic variations analysis in improving the current variant labeling system. First, we reviewed the basis of the system developed by the World Health Organization (WHO) for the labeling of SARS-CoV-2 genetic variants and the derivative adapted by the United States Centers for Disease Control and Prevention (CDC). Both labeling systems are based on the virus’ major attributes. However, we found that the labeling criteria of the SARS-CoV-2 variants derived from these attributes are not accurately defined and are used differently by the two agencies. Consequently, discrepancies exist between the labels given by WHO and the CDC to the same variants. Our observations suggest that giving the variant of concern (VOC) label to a new variant is premature and might not be appropriate. Therefore, we used a comparative computational approach to predict the effects of the mutations on the virus structure and functions of five VOCs. By linking these data to the criteria used by WHO/CDC for variant labeling, we ascertained that a predictive computational comparative approach of the genetic variations is a good way for rapid and more accurate labeling of SARS-CoV-2 variants. We propose to label all emergent variants, variant under monitoring or variant being monitored (VUM/VBM), and to carry out computational predictive studies with thorough comparison to existing variants, upon which more appropriate and informative labels can be attributed. Furthermore, harmonization of the variant labeling system would be globally beneficial to communicate about and fight the COVID-19 pandemic.

Knocking down Israa, the Zmiz1 intron-nested gene, unveils interrelated T cell activation functions in mouse

We previously reported Israa (immune-system-released activating agent), a novel gene nested in intron 6 of the mouse Zmiz1 gene. Zmiz1 is involved in several functions such as fertility and T cell development and its knockout leads to non-viable embryos. We also reported ISRAA's expression in lymphoid organs, particularly in the thymus CD³⁺ T cells during all developmental stages. In addition, we showed that ISRAA is a binding partner of Fyn and Elf-1 and regulates the expression of T cell activation-related genes in vitro. In this paper, we report the generation and characterization of an Israa ^-/- constitutive knockout mouse. The histological study shows that Israa ^-/- mice exhibit thymus and spleen hyperplasia. Israa ^-/- derived T cells showed increased proliferation compared to the wild-type mice T cells. Moreover, gene expression analysis revealed a set of differentially expressed genes in the knockout and wild-type animals during thymus development (mostly genes of T cell activation pathways). Immunological phenotyping of the thymocytes and splenocytes of Israa ^-/- showed no difference with those of the wild-type. Moreover, we observed that knocking out the Zmiz1 intron embedded Israa gene does not affect mice fertility, thus does not disturb this Zmiz1 function. The characterization of the Israa ^-/- mouse confirms the role ISRAA plays in the expression regulation of genes involved in T cell activation established in vitro. Taken together, our findings point toward a potential functional interrelation between the intron nested Israa gene and the Zmiz1 host gene in regulating T cell activation. This constitutively Israa ^-/- mice can be a good model to study T cell activation and to investigate the relationship between host and intron-nested genes.

A Computational Approach to Evaluate the Combined Effect of SARS-CoV-2 RBD Mutations and ACE2 Receptor Genetic Variants on Infectivity: The COVID-19 Host-Pathogen Nexus

SARS-CoV-2 infectivity is largely determined by the virus Spike protein binding to the ACE2 receptor. Meanwhile, marked infection rate differences were reported between populations and individuals. To understand the disease dynamic, we developed a computational approach to study the implications of both SARS-CoV-2 RBD mutations and ACE2 polymorphism on the stability of the virus-receptor complex. We used the 6LZG PDB RBD/ACE2 3D model, the mCSM platform, the LigPlot+ and PyMol software to analyze the data on SARS-CoV-2 mutations and ACE variants retrieved from GISAID and Ensembl/GnomAD repository. We observed that out of 351 RBD point mutations, 83% destabilizes the complex according to free energy (ΔΔG) differences. We also spotted variations in the patterns of polar and hydrophobic interactions between the mutations occurring in 15 out of 18 contact residues. Similarly, comparison of the effect on the complex stability of different ACE2 variants showed that the pattern of molecular interactions and the complex stability varies also according to ACE2 polymorphism. We infer that it is important to consider both ACE2 variants and circulating SARS-CoV-2 RBD mutations to assess the stability of the virus-receptor association and evaluate infectivity. This approach might offers a good molecular ground to mitigate the virus spreading.

Genetic variants of the EGFR ligand-binding domain and their association with structural alterations in Arab cancer patients

Objective

This study aimed to identify novel genetic variants in the CR2 extracellular domain of the epidermal growth factor receptor (EGFR) in healthy individuals and patients with six different types of adenocarcinoma, in Arabian peninsula populations. It also aimed to investigate the effects of these variants on the EGFR structure and their eventual relevance to tumorigenesis.

Results

We detected seven new EGFR genetic variants in 168 cancer patients and 114 controls. A SNP rs374670788 was more frequent in bladder cancer but not significantly associated to. However, a missense mutation (V550M) was significantly associated to colon, ovary, lung, bladder and thyroid cancer samples (p < 0.05). Three mutations (H590R, E602K and T605T) were found in the heterozygous form only in colon cancer patients. Genomic analysis of the synonymous mutation (G632G) showed that the T/A genotype could be associated to thyroid cancer in Arab patients (p < 0.05). An additional novel SNP rs571064657 was observed in control individuals. Computational analysis of the genetic variants revealed a reduction in the stabilization of the EGFR tethered form for both V550M and the common R521K variant with low energetic state (− ∆∆G). Molecular interactions analysis suggested that these mutations might affect the receptor’s function and promote tumorigenesis.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13104-021-05559-y.

A computer-aided approach to identify novel Leishmania major protein disulfide isomerase inhibitors for treatment of leishmaniasis

Leishmaniasis is an infectious disease caused by parasites of the genus Leishmania and transmitted by the bite of a sand fly. To date, most available drugs for treatment are toxic and beyond the economic means of those affected by the disease. Protein disulfide isomerase (PDI) is a chaperone protein that plays a major role in the folding of newly synthesized proteins, specifically assisting in disulfide bond formation, breakage, or rearrangement in all non-native proteins. In previous work, we demonstrated that Leishmania major PDI (LmPDI) has an essential role in pathogen virulence. Furthermore, inhibition of LmPDI further blocked parasite infection in macrophages. In this study, we utilized a computer-aided approach to design a series of LmPDI inhibitors. Fragment-based virtual screening allowed for the understanding of the inhibitors' modes of action on LmPDI active sites. The generated compounds obtained after multiple rounds of virtual screening were synthesized and significantly inhibited target LmPDI reductase activity and were shown to decrease in vitro parasite growth in human monocyte-derived macrophages. This novel cheminformatics and synthetic approach led to the identification of a new series of compounds that might be optimized into novel drugs, likely more specific and less toxic for the treatment of leishmaniasis.

Vascular endothelial ERp72 is involved in the inflammatory response in a rat model of skeletal muscle injury

The vascular inflammatory response involves the coordinated action of a large network of molecular mediators and culminates in the transmigration of leukocytes into the site of inflammation. Inflammatory mediators include a variety of protein families, including adhesion molecules such as integrins and members of the immunoglobulin superfamily, as well as other cytokines and chemokines. In this study, a rat model of traumatic skeletal muscle injury was used to demonstrate endoplasmic reticulum resident protein 72 (ERp72) overexpression in the early phase of the inflammatory response that follows skeletal muscle injury. Reverse transcription‑quantitative PCR, western blotting, dual‑labeling immunohistochemistry and immunofluorescence experiments confirmed that ERp72 was expressed on the endothelial cells of blood vessels present at the injured area. In addition, a cell‑based neutrophil adhesion assay indicated that a polyclonal antibody specific for ERp72 significantly reduced adhesion of neutrophils to activated human umbilical vein endothelial cells (35% reduction). These data suggested that ERp72 expression on vascular endothelial cells may play a role in skeletal muscle inflammation and could be considered as a target for the modulation of leukocyte‑endothelial cell interactions in an inflammatory setting.

Genetic analysis of TMPRSS6 gene in Saudi female patients with iron deficiency anemia

OBJECTIVE/BACKGROUND

Mutations in transmembrane protease serine 6 (TMPRSS6) gene induce high hepcidin level, which causes iron-refractory iron deficiency anemia (IRIDA) by preventing duodenal iron absorption. This study aims to identify the common genetic variations of the TMPRSS6 gene that affect iron levels among Saudi female patients with iron deficiency anemia (IDA).

METHODS

All study participants were Saudi females (12-49 years old): 32 patients with IDA, 32 patients with IRIDA, and 34 healthy individuals comprising the control group. Hematological investigations, iron profile, serum hepcidin level, and TMPRSS6 gene transcription were determined. The TMPRSS6 gene was amplified, sequenced, and analyzed among all study participants.

RESULTS

The mean hepcidin and TMPRSS6 RNA transcription levels in IDA and IRIDA groups were significantly lower than those in the control group. TMPRSS6 gene sequence analysis detected 41 variants: two in the 5' untranslated region (5'UTR), 17 in introns, and 22 in exons. Thirty-three variants were previously reported in the Single Nucleotide Polymorphism Database, and eight variants were novel; one novel variant was in 5'UTR (g.-2 T > G); five novel variants were detected in exons (p.W73X, p.D479N, p.E523K, p.L674L, and p.I799I). At the time of the sequence analysis of our samples, two variants-p.D479N and p.674L-were novel. However, these variants are present at a very low allele frequency in other populations (L674L, 0.00007761 and D479N, 0.000003980).

CONCLUSION

This is the first study to investigate the genetic variants of TMPRSS6 gene in Saudi female patients with IDA. The generated data will serve as a reference for future studies on IDA in the Arab population.

Four novel mutations in the mitochondrial ND4 gene of complex I in patients with multiple sclerosis

Multiple sclerosis (MS) is an immune-mediated neurological, inflammatory disease of the central nervous system. Recent studies have suggested that genetic variants in mitochondrial DNA (mtDNA)-encoded complexes of respiratory chain, particularly, complex I (NADH dehydrogenase), contribute to the pathogenicity of MS among different ethnicities, and targeting mitochondrial function may represent a novel approach for MS therapy. In this study, we sequenced ND genes (ND1, ND2, ND3, ND4, ND4L, ND5 and ND6) encoding subunits of complex I in 124 subjects, 60 patients with relapsing-remitting MS and 64 healthy individuals, in order to identify potential novel mutations in these patients. We found several variants in ND genes in both the patients and controls, and specific variants only in patients with MS. While the majority of these variants were synonymous, 4 variants in the ND4 gene were identified as missense mutations in patients with MS. Of these, m.11150G>A was observed in one patient, whereas m.11519A>C, m.11523A>C and m.11527C>T were observed in another patient. Functional analysis predicted the mutations, m.11519A>C, m.11523A>C and m.11150G>A, as deleterious with a direct impact on ND4 protein stability and complex I function, whereas m.11527C>T mutation had no effect on ND4 protein stability. However, the 3 mutations, m.11519A>C, m.11523A>C and m.11527C>T, which were observed in the same patient, were predicted to cause a cumulative destabilizing effect on ND4 protein, and could thus disrupt complex I function. On the whole, this study identified 4 novel mutations in the mtDNA-encoded ND4 gene in patients with MS, which could lead to complex I dysfunction, and further confirmed the implication of mtDNA mutations in the pathogenicity of MS. The identified novel mutations in patients with MS may be ethnic-related and may prove to be significant in personalized treatment.

The mouse intron-nested gene, Israa, is expressed in the lymphoid organs and involved in T-cell activation and signaling

We have previously reported Israa, immune-system-released activating agent, as a novel gene nested in intron 8 of the mouse Zmiz1 gene. We have also shown that Israa encodes for a novel FYN-binding protein and might be involved in the regulation of T-cell activation. In this report, we demonstrate that Israa gene product regulates the expression of a pool of genes involved in T-cell activation and signaling. Real time PCR and GFP knock-in expression analysis showed that Israa is transcribed and expressed in the spleen mainly by CD3⁺CD8⁺ cells as well as in the thymus by CD3⁺ (DP and DN), CD4⁺SP and CD8⁺SP cells at different developmental stages. We also showed that Israa is downregulated in T-cells following activation of T-cell receptor. Using yeast two-hybrid analysis, we identified ELF1, a transcription factor involved in T-cell regulation, as an ISRAA-binding partner. Transcriptomic analysis of an EL4 cell line overexpressing ISRAA revealed differential expression of several genes involved in T-cell signaling, activation and development. Among these genes, Prkcb, Mib2, Fos, Ndfip2, Cxxc5, B2m, Gata3 and Cd247 were upregulated whereas Itk, Socs3, Tigit, Ifng, Il2ra and FoxJ1 were downregulated. Our findings support the existence in mouse of a novel FYN-related T-cell regulation pathway involving the product of an intron-nested gene.

Engineering of the upper hinge region of human IgG1 Fc enhances the binding affinity to FcγIIIa (CD16a) receptor isoform

The interaction between antibodies and Immune cells surface FcγRIIIa (CD16a) receptor triggers a variety of immune responses including antibody-dependent cell-mediated cytotoxicity, antibody neutralization, phagocytosis, inflammation and tissue injury. Recent studies showed that IgG1 upper hinge region and FcγRs polymorphism play a major role in the interaction with Fcγ receptors and in the stability of the immune complex hence, in mounting strong inflammatory response. To further investigate this issue, we developed a tool box of IgG1 Fc isoforms to depict the affinity between mutated IgG1 Fc regions and extracellular domain variants (V158F) of CD16a. Our strategy consisted of designing different random upper-hinge mutated variants of IgG1 Fc domain, reproducing the naturally occurring two variants of CD16a and producing all of them as recombinant fusion proteins in Pichia Pastoris. The interactions were assayed using the Surface Plasmon Resonance (Biacore) method along with an in silico analysis to identify the major interaction and key residues that underline the affinity between the Fc region and CD16a variants. Our data showed that the affinity of the Fc region to the CD16a is strongly correlated to polar interactions. This molecular engineering approach yielded an IgG1Fc mutant with enhanced binding affinity to CD16a F158 variant.

A Central Nervous System-Dependent Intron-Embedded Gene Encodes a Novel Murine Fyn Binding Protein

The interplay between the nervous and immune systems is gradually being unraveled. We previously reported in the mouse the novel soluble immune system factor ISRAA, whose activation in the spleen is central nervous system-dependent. We also showed that ISRAA plays a role in modulating anti-infection immunity. Herein, we report the genomic description of the israa locus, along with some insights into the structure-function relationship of the protein. Our findings revealed that israa is nested within intron 6 of the mouse zmiz1 gene. Protein sequence analysis revealed a typical SH2 binding motif (Y₁₀₂TEV), with Fyn being the most likely binding partner. Docking simulation showed a favorable conformation for the ISRAA-Fyn complex, with a specific binding mode for the binding of the YTEV motif to the SH2 domain. Experimental studies showed that in vitro, recombinant ISRAA is phosphorylated by Fyn at tyrosine 102. Cell transfection and pull-down experiments revealed Fyn as a binding partner of ISRAA in the EL4 mouse T-cell line. Indeed, we demonstrated that ISRAA downregulates T-cell activation and the phosphorylation of an activation tyrosine (Y₄₁₆) of Src-family kinases in mouse splenocytes. Our observations highlight ISRAA as a novel Fyn binding protein that is likely to be involved in a signaling pathway driven by the nervous system.

Induction of dissociated cytokine profiles by ISRAA with selective critical involvement of ERK1/2 in its signaling functions

The immune system-released activating agent (ISRAA) is an immune mediator activated as a result of a nerve stimulus initiated by immune challenge. We have previously demonstrated that ISRAA and tumor necrosis factor (TNF) receptor 1 (TNFR1) share an interspecies-conserved motif (72% homology) that induces the apoptosis and proliferation of human peripheral blood mononuclear cells (hPBMCs) in a dose-dependent manner. In the present study, cytokine profiles were examined in response to the stimulation of hPBMCs with ISRAA. Furthermore, the signaling pathways induced by ISRAA were mapped. The results revealed high measurable levels of TNF-α, interleukin (IL)-6, IL-8, IL-10 and interferon (IFN)-γ, but not IL-4, IL-17 (IL-17A) or transforming growth factor (TGF)-β. The analysis of signaling pathways revealed the activation of extracellular-regulated protein kinase (ERK)1/2 as a downstream signal in the mitogen‑activated protein kinase (MAPK) pathway during TNF‑α and IL-6 production and apoptosis, but not during proliferation following stimulation with ISRAA by triggering the Fas-associated protein with death domain (FADD). STAT3 was found to be unphosphorylated in the ISRAA‑stimulated hPBMCs, and STAT3 was ubiquitously expressed in unstimulated cells, suggesting that ISRAA has a protein inhibitor of activated STAT (PIAS)-like activity, by functioning as a negative regulator of the effects of STAT3 on the Janus kinase (JAK)/STAT pathway. The determination of the nature of cytokine responses together with the signaling pathways of cellular activity induced by ISRAA paves the way for the investigation of a potential target of ISRAA and for the development of novel therapeutic approaches for the treatment of immune-regulated disorders.

Leishmania major protein disulfide isomerase as a drug target: enzymatic and functional characterization

Leishmaniasis is a major health problem worldwide and tools available for their control are limited. Effective vaccines are still lacking, drugs are toxic and expensive, and parasites develop resistance to chemotherapy. In this context, new antimicrobials are urgently needed to control the disease in both human and animal. Here, we report the enzymatic and functional characterization of a Leishmania virulence factor, Leishmania major Protein disulfide isomerase (LmPDI) that could constitute a potential drug target. LmPDI possesses domain structure organization similar to other PDI family members (a, a', b, b' and c domains), and it displays the three enzymatic and functional activities specific of PDI family members: isomerase, reductase and chaperone. These results suggest that LmPDI plays a key role in assisting Leishmania protein folding via its capacity to catalyze formation, breakage, and rearrangement of disulfide bonds in nascent polypeptides. Moreover, Bacitracin, a reductase activity inhibitor, and Ribostamycin, a chaperone activity inhibitor, were tested in LmPDI enzymatic assays and versus Leishmania promastigote in vitro cultures and Leishmania amastigote multiplication inside infected THP-1-derived macrophages. Bacitracin inhibited both isomerase and reductase activities, while Ribostamycin had no effect on the chaperone activity. Interestingly, Bacitracin blocked in vitro promastigote growth as well as amastigote multiplication inside macrophages with EC(50) values of 39 μM. These results suggest that LmPDI may constitute an interesting target for the development of new anti-Leishmania drugs.

A high-throughput turbidometric assay for screening inhibitors of Leishmania major protein disulfide isomerase

The use of a high-throughput technique to perform a pilot screen for Leishmania major protein disulfide isomerase (LmPDI) inhibitors identification is reported. In eukaryotic cells, protein disulfide isomerase (PDI) plays a crucial role in protein folding by catalyzing the rearrangement of disulfide bonds in substrate proteins following their synthesis. LmPDI displays similar domain structure organization and functional properties to other PDI family members and is involved in Leishmania virulence. The authors used a method based on the enzyme-catalyzed reduction of insulin in the presence of dithiothreitol. The screen of a small library of 1920 compounds was performed in a 384-well format and led to the identification of 27 compounds with inhibitory activity against LmPDI. The authors further tested the cytotoxicity of these compounds using Jurkat cells as well as their effect on Leishmania donovani amastigotes using high-content analysis. Results show hexachlorophene and a mixture of theaflavin monogallates inhibit Leishmania multiplication in infected macrophages derived from THP-1 cells, although the inhibitory effect on LmPDI enzymatic activity does not necessarily correlate with the antileishmanial activity.