Key gene modules and hub genes associated with pyrethroid and organophosphate resistance in Anopheles mosquitoes: a systems biology approach.

Indoor residual spraying (IRS) and insecticide-treated nets (ITNs) are the main methods used to control mosquito populations for malaria prevention. The efficacy of these strategies is threatened by the spread of insecticide resistance (IR), limiting the success of malaria control. Studies of the genetic evolution leading to insecticide resistance could enable the identification of molecular markers that can be used for IR surveillance and an improved understanding of the molecular mechanisms associated with IR. This study used a weighted gene co-expression network analysis (WGCNA) algorithm, a systems biology approach, to identify genes with similar co-expression patterns (modules) and hub genes that are potential molecular markers for insecticide resistance surveillance in Kenya and Benin. A total of 20 and 26 gene co-expression modules were identified via average linkage hierarchical clustering from Anopheles arabiensis and An. gambiae, respectively, and hub genes (highly connected genes) were identified within each module. Three specific genes stood out: serine protease, E3 ubiquitin-protein ligase, and cuticular proteins, which were top hub genes in both species and could serve as potential markers and targets for monitoring IR in these malaria vectors. In addition to the identified markers, we explored molecular mechanisms using enrichment maps that revealed a complex process involving multiple steps, from odorant binding and neuronal signaling to cellular responses, immune modulation, cellular metabolism, and gene regulation. Incorporation of these dynamics into the development of new insecticides and the tracking of insecticide resistance could improve the sustainable and cost-effective deployment of interventions.

An Integrated Computational Analysis of High-Risk SNPs in Angiopoietin-like Proteins (ANGPTL3 and ANGPTL8) Reveals Perturbed Protein Dynamics Associated with Cancer.

Angiopoietin-like proteins (ANGPTL) constitute a family of eight proteins (1-8) which play a pivotal role in the regulation of various pathophysiological processes. The current study sought to identify high-risk, "non-synonymous, single-nucleotide polymorphisms" (nsSNPs) in both ANGPTL3 and ANGPTL8 to evaluate the role that these nsSNPs play in various types of cancer. We retrieved a total of 301 nsSNPs from various databases; 79 of these candidates constitute high-risk nsSNPs. Moreover, we identified eleven high-risk nsSNPs that cause various types of cancer: seven candidates for ANGPTL3 (L57H, F295L, L309F, K329M, R332L, S348C, and G409R) and four candidates for ANGPTL8 (P23L, R85W, R138S, and E148D). Protein-protein interaction analysis revealed a strong association of ANGPTL proteins with several tumor-suppressor proteins such as ITGB3, ITGAV, and RASSF5. 'Gene-expression profiling interactive analysis' (GEPIA) showed that expression of ANGPTL3 is significantly downregulated in five cancers: sarcoma (SARC); cholangio carcinoma (CHOL); kidney chromophobe carcinoma (KICH); kidney renal clear cell carcinoma (KIRC); and kidney renal papillary cell carcinoma (KIRP). GEPIA also showed that expression of ANGPTL8 remains downregulated in three cancers: CHOL; glioblastoma (GBM); and breast invasive carcinoma (BRCA). Survival rate analysis indicated that both upregulation and downregulation of ANGPTL3 and ANGPTL8 leads to low survival rates in various types of cancer. Overall, the current study revealed that both ANGPTL3 and ANGPTL8 constitute potential prognostic biomarkers for cancer; moreover, nsSNPs in these proteins might lead to the progression of cancer. However, further in vivo investigation will be helpful to validate the role of these proteins in the biology of cancer.

Identification of Selective Novel Hits against Plasmodium falciparum Prolyl tRNA Synthetase Active Site and a Predicted Allosteric Site Using in silico Approaches.

Recently, there has been increased interest in aminoacyl tRNA synthetases (aaRSs) as potential malarial drug targets. These enzymes play a key role in protein translation by the addition of amino acids to their cognate tRNA. The aaRSs are present in all Plasmodium life cycle stages, and thus present an attractive malarial drug target. Prolyl tRNA synthetase is a class II aaRS that functions in charging tRNA with proline. Various inhibitors against Plasmodium falciparum ProRS (PfProRS) active site have been designed. However, none have gone through clinical trials as they have been found to be highly toxic to human cells. Recently, a possible allosteric site was reported in PfProRS with two possible allosteric modulators: glyburide and TCMDC-124506. In this study, we sought to identify novel selective inhibitors targeting PfProRS active site and possible novel allosteric modulators of this enzyme. To achieve this, virtual screening of South African natural compounds against PfProRS and the human homologue was carried out using AutoDock Vina. The modulation of protein motions by ligand binding was studied by molecular dynamics (MD) using the GROningen MAchine for Chemical Simulations (GROMACS) tool. To further analyse the protein global motions and energetic changes upon ligand binding, principal component analysis (PCA), and free energy landscape (FEL) calculations were performed. Further, to understand the effect of ligand binding on the protein communication, dynamic residue network (DRN) analysis of the MD trajectories was carried out using the MD-TASK tool. A total of ten potential natural hit compounds were identified with strong binding energy scores. Binding of ligands to the protein caused observable global and residue level changes. Dynamic residue network calculations showed increase in betweenness centrality (BC) metric of residues at the allosteric site implying these residues are important in protein communication. A loop region at the catalytic domain between residues 300 and 350 and the anticodon binding domain showed significant contributions to both PC1 and PC2. Large motions were observed at a loop in the Z-domain between residues 697 and 710 which was also in agreement with RMSF calculations that showed increase in flexibility of residues in this region. Residues in this loop region are implicated in ATP binding and thus a change in dynamics may affect ATP binding affinity. Free energy landscape (FEL) calculations showed that the holo protein (protein-ADN complex) and PfProRS-SANC184 complexes were stable, as shown by the low energy with very few intermediates and hardly distinguishable low energy barriers. In addition, FEL results agreed with backbone RMSD distribution plots where stable complexes showed a normal RMSD distribution while unstable complexes had multimodal RMSD distribution. The betweenness centrality metric showed a loss of functional importance of key ATP binding site residues upon allosteric ligand binding. The deep basins in average L observed at the allosteric region imply that there is high accessibility of residues at this region. To further analyse BC and average L metrics data, we calculated the ΔBC and ΔL values by taking each value in the holo protein BC or L matrix less the corresponding value in the ligand-bound complex BC or L matrix. Interestingly, in allosteric complexes, residues located in a loop region implicated in ATP binding had negative ΔL values while in orthosteric complexes these residues had positive ΔL values. An increase in contact frequency between residues Ser263, Thr267, Tyr285, and Leu707 at the allosteric site and residues Thr397, Pro398, Thr402, and Gln395 at the ATP binding TXE loop was observed. In summary, this study identified five potential orthosteric inhibitors and five allosteric modulators against PfProRS. Allosteric modulators changed ATP binding site dynamics, as shown by RMSF, PCA, and DRN calculations. Changes in dynamics of the ATP binding site and increased contact frequency between residues at the proposed allosteric site and the ATP binding site may explain how allosteric modulators distort the ATP binding site and thus might inhibit PfProRS. The scaffolds of the identified hits in the study can be used as a starting point for antimalarial inhibitor development with low human cytotoxicity.

Aminoacyl tRNA synthetases as malarial drug targets: a comparative bioinformatics study.

BACKGROUND: Treatment of parasitic diseases has been challenging due to evolution of drug resistant parasites, and thus there is need to identify new class of drugs and drug targets. Protein translation is important for survival of malarial parasite, Plasmodium, and the pathway is present in all of its life cycle stages. Aminoacyl tRNA synthetases are primary enzymes in protein translation as they catalyse amino acid addition to the cognate tRNA. This study sought to understand differences between Plasmodium and human aminoacyl tRNA synthetases through bioinformatics analysis. METHODS: Plasmodium berghei, Plasmodium falciparum, Plasmodium fragile, Plasmodium knowlesi, Plasmodium malariae, Plasmodium ovale, Plasmodium vivax, Plasmodium yoelii and human aminoacyl tRNA synthetase sequences were retrieved from UniProt database and grouped into 20 families based on amino acid specificity. These families were further divided into two classes. Both families and classes were analysed. Motif discovery was carried out using the MEME software, sequence identity calculation was done using an in-house Python script, multiple sequence alignments were performed using PROMALS3D and TCOFFEE tools, and phylogenetic tree calculations were performed using MEGA vs 7.0 tool. Possible alternative binding sites were predicted using FTMap webserver and SiteMap tool. RESULTS: Motif discovery revealed Plasmodium-specific motifs while phylogenetic tree calculations showed that Plasmodium proteins have different evolutionary history to the human homologues. Human aaRSs sequences showed low sequence identity (below 40%) compared to Plasmodium sequences. Prediction of alternative binding sites revealed potential druggable sites in PfArgRS, PfMetRS and PfProRS at regions that are weakly conserved when compared to the human homologues. Multiple sequence analysis, motif discovery, pairwise sequence identity calculations and phylogenetic tree analysis showed significant differences between parasite and human aaRSs proteins despite functional and structural conservation. These differences may provide a basis for further exploration of Plasmodium aminoacyl tRNA synthetases as potential drug targets. CONCLUSION: This study showed that, despite, functional and structural conservation, Plasmodium aaRSs have key differences from the human homologues. These differences in Plasmodium aaRSs can be targeted to develop anti-malarial drugs with less toxicity to the host.

Phytochemical analysis, in-vitro and in-silico study of antiproliferative activity of ethyl acetate fraction of Launaea cornuta (Hochst. ex Oliv. & Hiern) C. Jeffrey against human cervical cancer cell line.

Introduction: Cervical cancer is one of the leading causes of death among women globally due to the limitation of current treatment methods and their associated adverse side effects. Launaea cornuta is used as traditional medicine for the treatment of a variety of diseases including cancer. However, there is no scientific validation on the antiproliferative activity of L. cornuta against cervical cancer. Objective: This study aimed to evaluate the selective antiproliferative, cytotoxic and antimigratory effects of L. cornuta and to explore its therapeutical mechanisms in human cervical cancer cell lines (HeLa-229) through a network analysis approach. Materials and methods: The cytotoxic effect of L. cornuta ethyl acetate fraction on the proliferation of cervical cancer cells was evaluated by 3-(4, 5-dimethylthiazol-2-yl)-2, 5-diphenyltetrazolium bromide (MTT) bioassay and the antimigratory effect was assessed by wound healing assays. Compounds were analysed using the qualitative colour method and gas chromatography-mass spectroscopy (GC-MS). Subsequently, bioinformatic analyses, including the protein-protein interaction (PPI) network analysis, Gene Ontology (GO), and Kyoto Encyclopaedia of Genes and Genomes (KEGG) analysis, were performed to screen for potential anticervical cancer therapeutic target genes of L. cornuta. Molecular docking (MD) was performed to predict and understand the molecular interactions of the ligands against cervical cancer. Reverse transcription-quantitative polymerase chain reaction (RT-qPCR) was performed to validate the network analysis results. Results: L. cornuta ethyl acetate fraction exhibited a remarkable cytotoxic effect on HeLa-229 proliferation (IC50 of 20.56 ± 2.83 µg/mL) with a selectivity index (SI) of 2.36 with minimal cytotoxicity on non-cancerous cells (Vero-CCL 81 (IC50 of 48.83 ± 23.02). The preliminary screening revealed the presence of glycosides, phenols, saponins, terpenoids, quinones, and tannins. Thirteen compounds were also identified by GC-MS analysis. 124 potential target genes associated with the effect of L. cornuta ethyl acetate fraction on human cervical cancer were obtained, including AKT1, MDM2, CDK2, MCL1 and MTOR were identified among the top hub genes and PI3K/Akt1, Ras/MAPK, FoxO and EGFR signalling pathways were identified as the significantly enriched pathways. Molecular docking results showed that stigmasteryl methyl ether had a good binding affinity against CDK2, ATK1, BCL2, MDM2, and Casp9, with binding energy ranging from -7.0 to -12.6 kcal/mol. Tremulone showed a good binding affinity against TP53 and P21 with -7.0 and -8.0 kcal/mol, respectively. This suggests a stable molecular interaction of the ethyl acetate fraction of L. cornuta compounds with the selected target genes for cervical cancer. Furthermore, RT-qPCR analysis revealed that CDK2, MDM2 and BCL2 were downregulated, and Casp9 and P21 were upregulated in HeLa-229 cells treated with L. cornuta compared to the negative control (DMSO 0.2%). Conclusion: The findings indicate that L. cornuta ethyl acetate fraction phytochemicals modulates various molecular targets and pathways to exhibit selective antiproliferative and cytotoxic effects against HeLa-229 cells. This study lays a foundation for further research to develop innovative clinical anticervical cancer agents.

Novel dynamic residue network analysis approaches to study allosteric modulation: SARS-CoV-2 Mpro and its evolutionary mutations as a case study.

The rational search for allosteric modulators and the allosteric mechanisms of these modulators in the presence of mutations is a relatively unexplored field. Here, we established novel in silico approaches and applied them to SARS-CoV-2 main protease (Mpro) as a case study. First, we identified six potential allosteric modulators. Then, we focused on understanding the allosteric effects of these modulators on each of its protomers. We introduced a new combinatorial approach and dynamic residue network (DRN) analysis algorithms to examine patterns of change and conservation of critical nodes, according to five independent criteria of network centrality. We observed highly conserved network hubs for each averaged DRN metric on the basis of their existence in both protomers in the absence and presence of all ligands (persistent hubs). We also detected ligand specific signal changes. Using eigencentrality (EC) persistent hubs and ligand introduced hubs we identified a residue communication path connecting the allosteric binding site to the catalytic site. Finally, we examined the effects of the mutations on the behavior of the protein in the presence of selected potential allosteric modulators and investigated the ligand stability. One crucial outcome was to show that EC centrality hubs form an allosteric communication path between the allosteric ligand binding site to the active site going through the interface residues of domains I and II; and this path was either weakened or lost in the presence of some of the mutations. Overall, the results revealed crucial aspects that need to be considered in rational computational drug discovery.