Search

Abstract Identification of interactions between chemical compounds and proteins is crucial for various applications, including drug discovery, target identification, network pharmacology, and elucidation of protein functions. Deep neural network-based approaches are becoming increasingly popular in efficiently identifying compound-protein interactions with high-throughput capabilities, narrowing down the scope of candidates for traditional labor-intensive, time-consuming and expensive experimental techniques. In this study, we proposed an end-to-end approach termed SPVec-SGCN-CPI, which utilized simplified graph convolutional network (SGCN) model with low-dimensional and continuous features generated from our previously developed model SPVec and graph topology information to predict compound-protein interactions. The SGCN technique, dividing the local neighborhood aggregation and nonlinearity layer-wise propagation steps, effectively aggregates K-order neighbor information while avoiding neighbor explosion and expediting training. The performance of the SPVec-SGCN-CPI method was assessed across three datasets and compared against four machine learning- and deep learning-based methods, as well as six state-of-the-art methods. Experimental results revealed that SPVec-SGCN-CPI outperformed all these competing methods, particularly excelling in unbalanced data scenarios. By propagating node features and topological information to the feature space, SPVec-SGCN-CPI effectively incorporates interactions between compounds and proteins, enabling the fusion of heterogeneity. Furthermore, our method scored all unlabeled data in ChEMBL, confirming the top five ranked compound-protein interactions through molecular docking and existing evidence. These findings suggest that our model can reliably uncover compound-protein interactions within unlabeled compound-protein pairs, carrying substantial implications for drug re-profiling and discovery. In summary, SPVec-SGCN demonstrates its efficacy in accurately predicting compound-protein interactions, showcasing potential to enhance target identification and streamline drug discovery processes. Scientific contributions The methodology presented in this work not only enables the comparatively accurate prediction of compound-protein interactions but also, for the first time, take sample imbalance which is very common in real world and computation efficiency into consideration simultaneously, accelerating the target identification and drug discovery process.

Abstract The COVID-19 pandemic, caused by the SARS-CoV-2 virus, has led to over six million deaths worldwide. In human immune system, the type 1 interferon (IFN) pathway plays a crucial role in fighting viral infections. However, the ORF8 protein of the virus evade the immune system by interacting with IRF3, hindering its nuclear translocation and consequently downregulate the type I IFN signaling pathway. To block the binding of ORF8–IRF3 and inhibit viral pathogenesis a quick discovery of an inhibitor molecule is needed. Therefore, in the present study, the interface between the ORF8 and IRF3 was targeted on a high-affinity carbon nanotube by using computational tools. After analysis of 62 carbon nanotubes by multiple docking with the induced fit model, the top five compounds with high docking scores of − 7.94 kcal/mol, − 7.92 kcal/mol, − 7.28 kcal/mol, − 7.19 kcal/mol and − 7.09 kcal/mol (top hit1-5) were found to have inhibitory activity against the ORF8–IRF3 complex. Molecular dynamics analysis of the complexes revealed the high compactness of residues, stable binding, and strong hydrogen binding network among the ORF8-nanotubes complexes. Moreover, the total binding free energy for top hit1-5 was calculated to be − 43.21 ± 0.90 kcal/mol, − 41.17 ± 0.99 kcal/mol, − 48.85 ± 0.62 kcal/mol, − 43.49 ± 0.77 kcal/mol, and − 31.18 ± 0.78 kcal/mol respectively. These results strongly suggest that the identified top five nanotubes (hit1-5) possess significant potential for advancing and exploring innovative drug therapies. This underscores their suitability for subsequent in vivo and in vitro experiments, marking them as promising candidates worthy of further investigation.

Abstract Microproteins, prevalent across all kingdoms of life, play a crucial role in cell physiology and human health. Although global gene transcription is widely explored and abundantly available, our understanding of microprotein functions using transcriptome data is still limited. To mitigate this problem, we present a database, Mip-mining ( https://weilab.sjtu.edu.cn/mipmining/ ), underpinned by high-quality RNA-sequencing data exclusively aimed at analyzing microprotein functions. The Mip-mining hosts 336 sets of high-quality transcriptome data from 8626 samples and nine representative living organisms, including microorganisms, plants, animals, and humans, in our Mip-mining database. Our database specifically provides a focus on a range of diseases and environmental stress conditions, taking into account chemical, physical, biological, and diseases-related stresses. Comparatively, our platform enables customized analysis by inputting desired data sets with self-determined cutoff values. The practicality of Mip-mining is demonstrated by identifying essential microproteins in different species and revealing the importance of ATP15 in the acetic acid stress tolerance of budding yeast. We believe that Mip-mining will facilitate a greater understanding and application of microproteins in biotechnology. Moreover, it will be beneficial for designing therapeutic strategies under various biological conditions.

The authors would like to make the following correction to this published paper [...]

Spatial transcriptomic (ST) clustering employs spatial and transcription information to group spots spatially coherent and transcriptionally similar together into the same spatial domain. Graph convolution network (GCN) and graph attention network (GAT), fed with spatial coordinates derived adjacency and transcription profile derived feature matrix are often used to solve the problem. Our proposed method STGIC (spatial transcriptomic clustering with graph and image convolution) is designed for techniques with regular lattices on chips. It utilizes an adaptive graph convolution (AGC) to get high quality pseudo-labels and then resorts to dilated convolution framework (DCF) for virtual image converted from gene expression information and spatial coordinates of spots. The dilation rates and kernel sizes are set appropriately and updating of weight values in the kernels is made to be subject to the spatial distance from the position of corresponding elements to kernel centers so that feature extraction of each spot is better guided by spatial distance to neighbor spots. Self-supervision realized by Kullback-Leibler (KL) divergence, spatial continuity loss and cross entropy calculated among spots with high confidence pseudo-labels make up the training objective of DCF. STGIC attains state-of-the-art (SOTA) clustering performance on the benchmark dataset of 10x Visium human dorsolateral prefrontal cortex (DLPFC). Besides, it's capable of depicting fine structures of other tissues from other species as well as guiding the identification of marker genes. Also, STGIC is expandable to Stereo-seq data with high spatial resolution.

Blood-brain barrier (BBB) permeability prediction plays a pivotal role in drug discovery for neurological disorders which is essential for the development of central nervous system (CNS) drugs. Compounds having high permeability must be found to synthesize brain drugs for the treatment of different brain disorders such as Alzheimer's, Parkinson's, and brain tumors. Developing an accurate mathematical computational model to determine the exact brain permeability value for a possible drug is essential for advancing and improving the success rate of the development of drugs for neurological treatment. We developed a combined method capable of forecasting the logBB value of the compound in question for lightBBB and B3DB datasets by using the Convolutional neural network (CNNs), Recurrent Neural Networks (RNNs) and Ordinary Differential Equations (ODEs). The results demonstrate the overall assessment of the prediction ability of BBB permeability. CNN-LSTM model tests the performance for the prediction of logBB values in two different datasets LightBBB and B3DB. In both datasets, CNN-LSTM achieves lower RMSE values as compared to other models, showing that it has better predictive performance. Specifically, in the LightBBB dataset, the CNN-LSTM model achieves an RMSE of 0.59, while in the B3DB datasets, it achieves an even lower RMSE of 0.85. CNN-LSTM model shows highly effective in accurately predicting logBB values in both datasets.

Amid the ongoing monkeypox outbreak, there is an urgent need for the rapid development of effective therapeutic interventions capable of countering the immune evasion mechanisms employed by the monkeypox virus (MPXV). The evasion strategy involves the binding of the F3L protein to dsRNA, resulting in diminished interferon (IFN) production. Consequently, our current research focuses on utilizing virtual drug screening techniques to target the RNA binding domain of the F3L protein. Out of the 954 compounds within the South African natural compound database, only four demonstrated notable docking scores: −6.55, −6.47, −6.37, and −6.35 kcal/mol. The dissociation constant (KD) analysis revealed a stronger binding affinity of the top hits 1-4 (−5.34, −5.32, −5.29, and −5.36 kcal/mol) with the F3L in the MPXV. All-atom simulations of the top-ranked hits 1 to 4 consistently exhibited stable dynamics, suggesting their potential to interact effectively with interface residues. This was further substantiated through analyses of parameters such as radius of gyration (Rg), Root Mean Square Fluctuation, and hydrogen bonding. Cumulative assessments of binding free energy confirmed the top-performing candidates among all the compounds, with values of −35.90, −52.74, −28.17, and −32.11 kcal/mol for top hits 1-4, respectively. These results indicate that compounds top hit 1-4 could hold significant promise for advancing innovative drug therapies, suggesting their suitability for both in vivo and in vitro experiments.

The growing concerns and cases of COVID-19 with the appearance of novel variants i.e., BA.2.75. BA.5 and XBB have prompted demand for more effective treatment options that could overcome the risk of immune evasion. For this purpose, discovering novel small molecules to inhibit druggable proteins such as PLpro required for viral pathogenesis, replication, survival, and spread is the best choice. Compounds from the Dark chemical matter (DCM) database is consistently active in various screening tests and offer intriguing possibilities for finding drugs that are extremely selective or active against uncommon targets. Considering the essential role of PLpro, the current study uses DCMdatabase for the identification of potential hits using in silico virtual molecular screening and simulation approaches to inhibit the current and emerging variants of SARS-CoV-2. Our results revealed the 10 best compounds with docking scores between −7.99 to −7.03 kcal/mol better than the control drug (GRL0617) among which DC 5977–0726, DC 6623–2024, DC C879-0379 and DC D135-0154 were observed as the best hits. Structural-dynamics properties such as dynamic stability, protein packing, and residue flexibility demonstrated the pharmacologically favorable properties of these top hits in contrast to GRL0617. The hydrogen bonding half-life revealed that Asp164, Arg166, Tyr264, and Tyr268 have major contributions to the hydrogen bonding during the simulation. However, some of the important hydrogen bonds were missing in the control drug (GRL0617). Finally, the total binding free energy was reported to be −34.41 kcal/mol for GRL0617 (control), −41.03 kcal/mol for the DC5977-0726-PLpro, for the DC6623-2024-Plpro complex the TBE was −48.87 kcal/mol, for the for DCC879-0379-Plpro complex the TBE was −45.66 kcal/mol while for the DCD135-0154-PLpro complex the TBE was calculated to be −40.09 kcal/mol respectively, which shows the stronger potency of these compounds against PLpro and further in in vivo and in vitro test are required for the possible usage as potential drug against SARS-CoV-2.

Abstract The joint use of multiple drugs may cause unintended drug-drug interactions (DDIs) and result in adverse consequence to the patients. Accurate identification of DDI types can not only provide hints to avoid these accidental events, but also elaborate the underlying mechanisms by how DDIs occur. Several computational methods have been proposed for multi-type DDI prediction, but room remains for improvement in prediction performance. In this study, we propose a supervised contrastive learning based method, MDDI-SCL, implemented by three-level loss functions, to predict multi-type DDIs. MDDI-SCL is mainly composed of three modules: drug feature encoder and mean squared error loss module, drug latent feature fusion and supervised contrastive loss module, multi-type DDI prediction and classification loss module. The drug feature encoder and mean squared error loss module uses self-attention mechanism and autoencoder to learn drug-level latent features. The drug latent feature fusion and supervised contrastive loss module uses multi-scale feature fusion to learn drug pair-level latent features. The prediction and classification loss module predicts DDI types of each drug pair. We evaluate MDDI-SCL on three different tasks of two datasets. Experimental results demonstrate that MDDI-SCL achieves better or comparable performance as the state-of-the-art methods. Furthermore, the effectiveness of supervised contrastive learning is validated by ablation experiment, and the feasibility of MDDI-SCL is supported by case studies. The source codes are available at https://github.com/ShenggengLin/MDDI-SCL .

Abstract Background Previous studies on European (EUR) samples have obtained inconsistent results regarding the genetic correlation between type 2 diabetes mellitus (T2DM) and Schizophrenia (SCZ). A large-scale trans-ethnic genetic analysis may provide additional evidence with enhanced power. Objective We aimed to explore the genetic basis for both T2DM and SCZ based on large-scale genetic analyses of genome-wide association study (GWAS) data from both East Asian (EAS) and EUR subjects. Methods A range of complementary approaches were employed to cross-validate the genetic correlation between T2DM and SCZ at the whole genome, autosomes (linkage disequilibrium score regression, LDSC), loci (Heritability Estimation from Summary Statistics, HESS), and causal variants (MiXeR and Mendelian randomization, MR) levels. Then, genome-wide and transcriptome-wide cross-trait/ethnic meta-analyses were performed separately to explore the effective shared organs, cells and molecular pathways. Results A weak genome-wide negative genetic correlation between SCZ and T2DM was found for the EUR (rg = − 0.098, P = 0.009) and EAS (rg =- 0.053 and P = 0.032) populations, which showed no significant difference between the EUR and EAS populations (P = 0.22). After Bonferroni correction, the rg remained significant only in the EUR population. Similar results were obtained from analyses at the levels of autosomes, loci and causal variants. 25 independent variants were firstly identified as being responsible for both SCZ and T2DM. The variants associated with the two disorders were significantly correlated to the gene expression profiles in the brain (P = 1.1E-9) and pituitary gland (P = 1.9E-6). Then, 61 protein-coding and non-coding genes were identified as effective genes in the pituitary gland (P

Introduction: Hypoxia-inducible factor (HIF) prolyl hydroxylase domain (PHD) enzymes are major therapeutic targets of anemia and ischemic/hypoxia diseases. To overcome safety issues, liver failure, and problems associated with on-/off-targets, natural products due to their novel and unique structures offer promising alternatives as drug targets.Methods: In the current study, the Marine Natural Products, North African, South African, East African, and North-East African chemical space was explored for HIF-PHD inhibitors discovery through molecular search, and the final hits were validated using molecular simulation and free energy calculation approaches.Results: Our results revealed that CMNPD13808 with a docking score of −8.690 kcal/mol, CID15081178 with a docking score of −8.027 kcal/mol, CID71496944 with a docking score of −8.48 kcal/mol and CID11821407 with a docking score of −7.78 kcal/mol possess stronger activity than the control N-[(4-hydroxy-8-iodoisoquinolin-3-yl)carbonyl]glycine, 4HG (−6.87 kcal/mol). Interaction analysis revealed that the target compounds interact with Gln239, Tyr310, Tyr329, Arg383 and Trp389 residues, and chelate the active site iron in a bidentate manner in PHD2. Molecular simulation revealed that these target hits robustly block the PHD2 active site by demonstrating stable dynamics. Furthermore, the half-life of the Arg383 hydrogen bond with the target ligands, which is an important factor for PHD2 inhibition, remained almost constant in all the complexes during the simulation. Finally, the total binding free energy of each complex was calculated as CMNPD13808-PHD2 −72.91 kcal/mol, CID15081178-PHD2 −65.55 kcal/mol, CID71496944-PHD2 −68.47 kcal/mol, and CID11821407-PHD2 −62.06 kcal/mol, respectively.Conclusion: The results show the compounds possess good activity in contrast to the control drug (4HG) and need further in vitro and in vivo validation for possible usage as potential drugs against HIF-PHD2-associated diseases.

Introduction: The current XBB variant of SARS-CoV-2 with the strongest immune escaping properties is currently the most dominant variant circulating around the world. With the emergence of XBB global morbidities and mortalities have raised again. In the current scenario, it was highly required to delineate the binding capabilities of NTD of XBB subvariant towards human neutralizing antibodies and to dig out the binding affinity of RBD of XBB subvariant with ACE2 receptor.Materials and Methods: The current study uses molecular interaction and simulation-based approaches to decipher the binding mechanism of RBD with ACE2 and mAb interaction with NTD of the spike protein.Results: Molecular docking of the Wild type NTD with mAb revealed a docking score of −113.2 ± 0.7 kcal/mol while XBB NTD docking with mAb reported −76.2 ± 2.3 kcal/mol. On the other hand, wild-type RBD and XBB RBD with ACE2 receptor demonstrated docking scores of −115.0 ± 1.5 kcal/mol and −120.8 ± 3.4 kcal/mol respectively. Moreover, the interaction network analysis also revealed significant variations in the number of hydrogen bonds, salt-bridges, and non-bonded contacts. These findings were further validated by computing the dissociation constant (KD). Molecular simulation analysis such as RMSD, RMSF, Rg and hydrogen bonding analysis revealed variation in the dynamics features of the RBD and NTD complexes due to the acquired mutations. Furthermore, the total binding energy for the wild-type RBD in complex with ACE2 reported −50.10 kcal/mol while XBB-RBD coupled with ACE2 reported −52.66 kcal/mol respectively. This shows though the binding of XBB is slightly increased but due to the variation in the bonding network and other factors makes the XBB variant to enter into the host cell efficiently than the wild type. On the other hand, the total binding free energy for the wildtype NTD-mAb was calculated to be −65.94 kcal/mol while for XBB NTD-mAb was reported to be −35.06 kcal/mol respectively. The significant difference in the total binding energy factors explains that the XBB variant possess stronger immune evasion properties than the others variants and wild type.Conclusions: The current study provides structural features for the XBB variant binding and immune evasion which can be used to design novel therapeutics.

Streptococcus gordonii is an oral bacterium colonizing the dental cavity and leading to plaque formation. This pervasive colonizer is also the etiologic agent of bacterial endocarditis and has a major role in infective endocarditis. The bacteria reach the heart through oral bleeding, leading to inflammation of cardiovascular valves. Over the past 50 years, it has shown a significant pathogenic role in immunocompromised and neutropenic patients. Since antibiotic resistance has created prophylaxis failure towards infective endocarditis, a potent therapeutic candidate is needed. Therefore, multi-epitopes vaccine offers advantages over the other approaches. Thus, herein, numerous molecular-omics tools were exploited to mine immunogenic peptides, i.e., T-cell and B-cell epitopes, and construct a vaccine sequence. Our findings revealed a total of 24 epitopes, including CTL, HTL, and B-cell are responsible for imparting immune responses, which were combined with the help of different linkers, and MEVC was constructed. Multifactorial validation of the candidate vaccine was performed to minimize the risk factors. The final sequence was docked with TLR2 to validate its conformation compatibility with receptor and long-term interactions stability. Our analysis revealed that the vaccine construct is immunogenic and non-allergenic. The construct also established various contacts with the immune receptor. Finally, the vaccine sequence was reverse-translated, optimized for codon usage, and analyzed for expression in the Escherichia coli K12 strain. Maximum expression was noted with a CAI score of 0.95. In silico immune simulation revealed that the antigen was neutralized on the 3rd day after injection. In conclusion, the current study warrants validation of the vaccine construct both in in vitro and in vivo models for accurate therapeutic intervention.

The SARS COV-2 and its variants are spreading around the world at an alarming speed, due to its higher transmissibility and the conformational changes caused by mutations. The resulting COVID-19 pandemic has imposed severe health consequences on human health. Several countries of the world including Pakistan have studied its genome extensively and provided productive findings. In the current study, the mCSM, DynaMut2, and I-Mutant servers were used to analyze the effect of identified mutations on the structural stability of spike protein however, the molecular docking and simulations approaches were used to evaluate the dynamics of the bonding network between the wild-type and mutant spike proteins with furin. We addressed the mutational modifications that have occurred in the spike protein of SARS-COV-2 that were found in 215 Pakistani's isolates of COVID-19 patients to study the influence of mutations on the stability of the protein and its interaction with the host cell. We found 7 single amino acid substitute mutations in various domains that reside in spike protein. The H49Y, N74K, G181V, and G446V were found in the S1 domain while the D614A, V622F, and Q677H mutations were found in the central helices of the spike protein. Based on the observation, G181V, G446V, D614A, and V622F mutants were found highly destabilizing and responsible for structural perturbation. Protein-protein docking and molecular simulation analysis with that of furin have predicted that all the mutants enhanced the binding efficiency however, the V622F mutant has greatly altered the binding capacity which is further verified by the KD value (7.1 E−14) and therefore may enhance the spike protein cleavage by Furin and increase the rate of infectivity by SARS-CoV-2. On the other hand, the total binding energy for each complex was calculated which revealed −50.57 kcal/mol for the wild type, for G181V −52.69 kcal/mol, for G446V −56.44 kcal/mol, for D614A −59.78 kcal/mol while for V622F the TBE was calculated to be −85.84 kcal/mol. Overall, the current finding shows that these mutations have increased the binding of Furin for spike protein and shows that D614A and V622F have significant effects on the binding and infectivity.

Among the different drug targets of SARS-CoV-2, a multi-domain protein known as NSP3 is a critical element of the translational and replication machinery. The macrodomain-I, in particular, has been reported to have an essential role in the viral attack on the innate immune response. In this study, we explore natural medicinal compounds and identify potential inhibitors to target the SARS-CoV-2–NSP3 macrodomain-I. Computational modeling and simulation tools were utilized to investigate the structural-dynamic properties using triplicates of 100 ns MD simulations. In addition, the MM/GBSA method was used to calculate the total binding free energy of each inhibitor bound to macrodomain-I. Two significant hits were identified: 3,5,7,4′-tetrahydroxyflavanone 3′-(4-hydroxybenzoic acid) and 2-hydroxy-3-O-beta-glucopyranosyl-benzoic acid. The structural-dynamic investigation of both compounds with macrodomain-I revealed stable dynamics and compact behavior. In addition, the total binding free energy for each complex demonstrated a robust binding affinity, of ΔG −61.98 ± 0.9 kcal/mol for Compound A, while for Compound B, the ΔG was −45.125 ± 2.8 kcal/mol, indicating the inhibitory potential of these compounds. In silico bioactivity and dissociation constant (KD) determination for both complexes further validated the inhibitory potency of each compound. In conclusion, the aforementioned natural products have the potential to inhibit NSP3, to directly rescue the host immune response. The current study provides the basis for novel drug development against SARS-CoV-2 and its variants.

Introduction: Coronavirus disease 2019 (COVID-19), caused by SARS-CoV-2, has had a disastrous effect worldwide during the previous three years due to widespread infections with SARS-CoV-2 and its emerging variations. More than 674 million confirmed cases and over 6.7 million deaths have been attributed to successive waves of SARS-CoV-2 infections as of 29th January 2023. Similar to other RNA viruses, SARS-CoV-2 is more susceptible to genetic evolution and spontaneous mutations over time, resulting in the continual emergence of variants with distinct characteristics. Spontaneous mutations of SARS-CoV-2 variants increase its transmissibility, virulence, and disease severity and diminish the efficacy of therapeutics and vaccines, resulting in vaccine-breakthrough infections and re-infection, leading to high mortality and morbidity rates.Materials and methods: In this study, we evaluated 10,531 whole genome sequences of all reported variants globally through a computational approach to assess the spread and emergence of the mutations in the SARS-CoV-2 genome. The available data sources of NextCladeCLI 2.3.0 (https://clades.nextstrain.org/) and NextStrain (https://nextstrain.org/) were searched for tracking SARS-CoV-2 mutations, analysed using the PROVEAN, Polyphen-2, and Predict SNP mutational analysis tools and validated by Machine Learning models.Result: Compared to the Wuhan-Hu-1 reference strain NC 045512.2, genome-wide annotations showed 16,954 mutations in the SARS-CoV-2 genome. We determined that the Omicron variant had 6,307 mutations (retrieved sequence:1947), including 67.8% unique mutations, more than any other variant evaluated in this study. The spike protein of the Omicron variant harboured 876 mutations, including 443 deleterious mutations. Among these deleterious mutations, 187 were common and 256 were unique non-synonymous mutations. In contrast, after analysing 1,884 sequences of the Delta variant, we discovered 4,468 mutations, of which 66% were unique, and not previously reported in other variants. Mutations affecting spike proteins are mostly found in RBD regions for Omicron, whereas most of the Delta variant mutations drawn to focus on amino acid regions ranging from 911 to 924 in the context of epitope prediction (B cell & T cell) and mutational stability impact analysis protruding that Omicron is more transmissible.Discussion: The pathogenesis of the Omicron variant could be prevented if the deleterious and persistent unique immunosuppressive mutations can be targeted for vaccination or small-molecule inhibitor designing. Thus, our findings will help researchers monitor and track the continuously evolving nature of SARS-CoV-2 strains, the associated genetic variants, and their implications for developing effective control and prophylaxis strategies.

IntroductionDrought is the largest abiotic factor impacting agriculture. Plants are challenged by both natural and artificial stressors because they are immobile. To produce drought-resistant plants, we need to know how plants react to drought. A largescale proteome study of leaf and root tissue found drought-responsive proteins. Tomato as a vegetable is grown worldwide. Agricultural biotechnology focuses on creating drought-resistant cultivars. This is important because tomato drought is so widespread. Breeders have worked to improve tomato quality, production, and stress resistance. Conventional breeding approaches have only increased drought tolerance because of drought’s complexity. Many studies have examined how tomatoes handle drought. With genomics, transcriptomics, proteomics, metabolomics, and modern sequencing technologies, it’s easier to find drought-responsive genes.MethodBiotechnology and in silico studies has helped demonstrate the function of drought-sensitive genes and generate drought-resistant plant types. The latest tomato genome editing technology is another. WRKY genes are plant transcription factors. They help plants grow and respond to both natural and artificial stimuli. To make plants that can handle stress, we need to know how WRKY-proteins, which are transcription factors, work with other proteins and ligands in plant cells by molecular docking and modeling study.ResultAbscisic acid, a plant hormone generated in stressed roots, was used here to make plants drought-resistant. Abscisic acid binds WRKY with binding affinity -7.4kcal/mol and inhibitory concentration (Ki) 0.12 microM.DiscussionThis study aims to modulate the transcription factor so plants can handle drought and stress better. Therefore, polyphenols found to make Solanum lycopersicum more drought-tolerant.

IntroductionThe COVID-19 pandemic has posed a major burden on healthcare and economic systems across the globe for over 3 years. Even though vaccines are available, the pathogenesis is still unclear. Multiple studies have indicated heterogeneity of immune responses to SARS-CoV-2, and potentially distinct patient immune types that might be related to disease features. However, those conclusions are mainly inferred by comparing the differences of pathological features between moderate and severe patients, some immunological features may be subjectively overlooked.MethodsIn this study, the relevance scores(RS), reflecting which features play a more critical role in the decision-making process, between immunological features and the COVID-19 severity are objectively calculated through neural network, where the input features include the immune cell counts and the activation marker concentrations of particular cell, and these quantified characteristic data are robustly generated by processing flow cytometry data sets containing the peripheral blood information of COVID-19 patients through PhenoGraph algorithm.ResultsSpecifically, the RS between immune cell counts and COVID-19 severity with time indicated that the innate immune responses in severe patients are delayed at the early stage, and the continuous decrease of classical monocytes in peripherial blood is significantly associated with the severity of disease. The RS between activation marker concentrations and COVID-19 severity suggested that the down-regulation of IFN-γ in classical monocytes, Treg, CD8 T cells, and the not down-regulation of IL_17a in classical monocytes, Tregs are highly correlated with the occurrence of severe disease. Finally, a concise dynamic model of immune responses in COVID-19 patients was generalized.DiscussionThese results suggest that the delayed innate immune responses in the early stage, and the abnormal expression of IL-17a and IFN-γ in classical monocytes, Tregs, and CD8 T cells are primarily responsible for the severity of COVID-19.

In the past few years, various somatic point mutations of isocitrate dehydrogenase (IDH) encoding genes (IDH1 and IDH2) have been identified in a broad range of cancers, including glioma. Despite the important function of IDH1 in tumorigenesis and its very polymorphic nature, it is not yet clear how different nsSNPs affect the structure and function of IDH1. In the present study, we employed different machine learning algorithms to screen nsSNPs in the IDH1 gene that are highly deleterious. From a total of 207 SNPs, all of the servers classified 80 mutations as deleterious. Among the 80 deleterious mutations, 14 were reported to be highly destabilizing using structure-based prediction methods. Three highly destabilizing mutations G15E, W92G, and I333S were further subjected to molecular docking and simulation validation. The docking results and molecular simulation analysis further displayed variation in dynamics features. The results from molecular docking and binding free energy demonstrated reduced binding of the drug in contrast to the wild type. This, consequently, shows the impact of these deleterious substitutions on the binding of the small molecule. PCA (principal component analysis) and FEL (free energy landscape) analysis revealed that these mutations had caused different arrangements to bind small molecules than the wild type where the total internal motion is decreased, thus consequently producing minimal binding effects. This study is the first extensive in silico analysis of the IDH1 gene that can narrow down the candidate mutations for further validation and targeting for therapeutic purposes.

Abstract Tuberculosis (TB), caused by Mycobacterium tuberculosis, is endemic in Pakistan. Resistance to both firstline rifampicin and isoniazid drugs (multidrug-resistant TB; MDR-TB) is hampering disease control. Rifampicin resistance is attributed to rpoB gene mutations, but rpoA and rpoC loci may also be involved. To characterise underlying rifampicin resistance mutations in the TB endemic province of Khyber Pakhtunkhwa, we sequenced 51 M. tuberculosis isolates collected between 2016 and 2019; predominantly, MDR-TB (n = 44; 86.3%) and lineage 3 (n = 30, 58.8%) strains. We found that known mutations in rpoB (e.g. S405L), katG (e.g. S315T), or inhA promoter loci explain the MDR-TB. There were 24 unique mutations in rpoA, rpoB, and rpoC genes, including four previously unreported. Five instances of within-host resistance diversity were observed, where two were a mixture of MDR-TB strains containing mutations in rpoB, katG, and the inhA promoter region, as well as compensatory mutations in rpoC. Heteroresistance was observed in two isolates with a single lineage. Such complexity may reflect the high transmission nature of the Khyber Pakhtunkhwa setting. Our study reinforces the need to apply sequencing approaches to capture the full-extent of MDR-TB genetic diversity, to understand transmission, and to inform TB control activities in the highly endemic setting of Pakistan.

C99 is the immediate precursor for amyloid beta (Aβ) and therefore is a central intermediate in the pathway that is believed to result in Alzheimer’s disease (AD). It has been suggested that cholesterol is associated with C99, but the dynamic details of how cholesterol affects C99 assembly and the Aβ formation remain unclear. To investigate this question, we employed coarse-grained and all-atom molecular dynamics simulations to study the effect of cholesterol and membrane composition on C99 dimerization. We found that although the existence of cholesterol delays C99 dimerization, there is no direct competition between C99 dimerization and cholesterol association. In contrast, the existence of cholesterol makes the C99 dimer more stable, which presents a cholesterol binding C99 dimer model. Cholesterol and membrane composition change the dimerization rate and conformation distribution of C99, which will subsequently influence the production of Aβ. Our results provide insights into the potential influence of the physiological environment on the C99 dimerization, which will help us understand Aβ formation and AD’s etiology.

Computer-aided drug design is an efficient method to analyze the development of disease-related drugs. However, developed as binding targets, medicines perform well in cell models and animal models but fail in human models. One main reason for this failure is that the human body has natural barriers, such as the blood-brain barrier, to block exogenous macromolecules. Thus, efficient and accurate predictions of drug molecules that can effectively pass the blood-brain barrier is necessary in developing drug treatments for brain tissue diseases. In this study, 7658 molecular structure features were extracted from 2354 drug molecule SMILE strings using computational methods. By integrating three feature selection algorithms of machine learning, 33 chemical structure features with significantly discriminant performance were screened out and used to construct multiple discriminant models. After a comprehensive comparison, the XGBoost model was selected as the final prediction model. After data preprocessing and parameter optimization, the model achieved 95% accuracy on the training set. To verify the model's stability, we introduced an external data set, which reached 96% accuracy of the model. This study applies new resampling methods and machine learning algorithms, and adjusts the application of resampling methods to obtain new chemical features to construct machine learning predictors. The features may contribute to the significant drug development that integrates biological analysis and machine learning algorithms.

The Wnt-activated $\beta $ -catenin and $Ca^{2+}$ ions play a critical role in the regulation of physiology of cardiomyocytes. The dysregulation of both the components causes the replacement of myocardial mass right ventricle with fibrous and adipose tissue which results in the condition of Arrhythmogenic Right ventricular Cardiomyopathy (ARVC). The major hurdle in ARVC treatment is lack of effective therapies targeting the underlying root molecular biomarkers of pathogenesis. Despite major advancements in interpreting the mechanism of ARVC etiology, the dynamics of molecular links in underlying biological machinery are still being delineated. Previously, formal methods based computational modeling techniques including kinetic logic, Petri Nets, hybrid automata and static analysis have greatly contributed in increasing our comprehension to decipher the molecular systems dynamics. It has allowed the identification of biomarkers which can be utilized for target-based therapies owing to meticulous biological abstractions along with implementing reference map that confines together the discrete biological insights. In this study, we have performed the static analysis of the Biological Regulatory Networks (BRNs) of the Wnt/ $\beta $ -catenin and Wnt/Ca 2+ signaling pathways to identify the significant biomarkers for ARVC. The abstracted qualitative models of afore mentioned BRNs are first constructed in GINsim software tool and then these models are converted into Automata Network (ANs) Models using Pint software tool. The fix point analysis is performed which contributed in pinpointing the possible therapeutic strategies for ARVC treatment by identification of drug targets such as Gsk3, Ck1 and Axin in Wnt/ $\beta $ -catenin AN and Bak & Bax, Parp, mCalpain, JNk and CIn in Wnt/Ca 2+ AN. Moreover, the Bcl2 gene is identified as novel therapeutic remedy in both the ANs of ARVC. The Bcl2 gene prevents the cardiac apoptosis via positive regulation of Wnt in Wnt/ $\beta $ -catenin AN and through inhibition of Bak & Bax (apoptotic component) in Wnt/Ca 2+ AN. The current study tends to fulfill the scientific gap between wet lab studies and provides cost effective and time saving computational strategies for an effectual treatment for deadly diseases like ARVC.

The emergence of recent SARS-CoV-2 has become a global health issue. This single-stranded positive-sense RNA virus is continuously spreading with increasing morbidities and mortalities. The proteome of this virus contains four structural and sixteen nonstructural proteins that ensure the replication of the virus in the host cell. However, the role of phosphoprotein (N) in RNA recognition, replicating, transcribing the viral genome, and modulating the host immune response is indispensable. Recently, the NMR structure of the N-terminal domain of the Nucleocapsid Phosphoprotein has been reported, but its precise structural mechanism of how the ssRNA interacts with it is not reported yet. Therefore, here, we have used an integrated computational pipeline to identify the key residues, which play an essential role in RNA recognition. We generated multiple variants by using an alanine scanning strategy and performed an extensive simulation for each system to signify the role of each interfacial residue. Our analyses suggest that residues T57A, H59A, S105A, R107A, F171A, and Y172A significantly affected the dynamics and binding of RNA. Furthermore, per-residue energy decomposition analysis suggests that residues T57, H59, S105 and R107 are the key hotspots for drug discovery. Thus, these residues may be useful as potential pharmacophores in drug designing.

Achromobacter xylosoxidans is a rod-shaped Gram-negative bacterium linked with causing several infections which mostly includes hematological malignancies. It has been recently reported to be associated with the development and progression of lung cancer and is an emerging respiratory disease-causing bacterium. The treatment of individuals infected with A. xylosoxidans bacteremia is difficult due to the fact that this pathogen has both intrinsic and acquired resistance mechanisms, typically resulting in a phenotype of multidrug resistance (MDR). Efforts are needed to design effective therapeutic strategies to curtail the emergence of this bacterium. Computational vaccine designing has proven its effectiveness, specificity, safety, and stability compared to conventional approaches of vaccine development. Therefore, the whole proteome of A. xylosoxidans was screened for the characterization of potential vaccine targets through subtractive proteomics pipeline for therapeutics design. Annotation of the whole proteome confirmed the three immunogenic vaccine targets, such as (E3HHR6), (E3HH04), and (E3HWA2), which were used to map the putative immune epitopes. The shortlisted epitopes, specific against Cytotoxic T Lymphocytes, Helper T-cell Lymphocytes, and linear B-Cell, were used to design the mRNA and multi-epitopes vaccine (MEVC). Initial validations confirmed the antigenic and non-allergenic properties of these constructs, followed by docking with the immune receptor, TLR-5, which resulted in robust interactions. The interaction pattern that followed in the docking complex included formation of 5 hydrogen bonds, 2 salt bridges, and 165 non-bonded contacts. This stronger binding affinity was also assessed through using the mmGBSA approach, showing a total of free binding energy of −34.64 kcal/mol. Further validations based on in silico cloning revealed a CAI score of 0.98 and an optimal percentage of GC contents (54.4%) indicated a putatively higher expression of the vaccine construct in Escherichia coli. Moreover, immune simulation revealed strong antibodies production upon the injection of the designed MEVC that resulted in the highest peaks of IgM+ IgG production (>3,500) between 10 and 15 days. In conclusion the current study provide basis for vaccine designing against the emerging A. xylosoxidans, which demands further experimental studies for in vitro and in vivo validations.

Elicited artificial in vitro cultures are gaining more interest due to their uniform biosynthesis of industrially valuable secondary metabolites. In this study, a unique methodology was applied, in which different ratios of gold (Au) and silver (Ag) nanoparticles (NPs) were supplemented to submerge cultures to investigate sustainable production of biomass and antioxidant secondary metabolites. Cell suspension cultures were exposed to Ag and AuNPs alone or different ratios of AgAuNPs (1:2; 1:3; 2:1; 3:1) in combination with NAA. The combination of AgAuNPs (3:1) with NAA enhanced fresh (9.25 g/100 ml) and dry biomass (0.64 g/100 ml) of suspended cells than control (6.67; 0.233 g/100 ml). AuNPs with NAA-augmented media enhanced biomass accumulation in lag, log and stationary phases in a period of 49 days. Furthermore, AgAu (3:1) and AgAuNPs (2:1; 1:2) with NAA enhanced protein contents, peroxidase and superoxide dismutase enzymes. However, maximum phenolics (TPC; 10.61 mg/g-DW) and flavonoids (7.62 mg/g-DW) were observed in cell cultures exposed to a combination of AgAuNPs (1:3) and NAA than control (6.27, 5.49 mg/g-DW). The combination of AgAuNPs (2:1) with NAA enhanced antioxidant activity (87.85%) in cell cultures. This study will help in illuminating the impact of NPs on cultures development and production of natural antioxidants.

The emergence of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) and its new variants reported in different countries have posed a serious threat to human health and social fabrics worldwide. In addition, these new variants hindered the efforts of vaccines and other therapeutic developments. In this review article, we explained the emergence of new variants of SARS-CoV-2, their transmission risk, mortality rate, and, more importantly, the impact of each new variant on the efficacy of the developed vaccines reported in different literature and findings. The literature reported that with the emergence of new variants, the efficacy of different vaccines is declined, hospitalization and the risk of reinfection is increased. The reports concluded that the emergence of a variant that entirely evades the immune response triggered by the vaccine is improbable. The emergence of new variants and reports of re-infections are creating a more distressing situation and therefore demands further investigation to formulate an effective therapeutic strategy.

Abstract Background Identifying reliable predictive markers is important to make therapeutic decisions, and determine the prognosis for acute myeloid leukemia (AML) patients. However, approximately 50% patients could not be accurately predicted by existing risk factors. It is necessary to identify novel prognostic factors to subdivide the intermediate-risk group or patients without any cytogenetic and molecular abnormalities. Methods Kaplan–Meier and Cox regression were used for survival analyses in three independent AML datasets. Analyses integrating both bioinformatics and ChIP-qPCR experiments were performed to explore the role of CEBPE in regulating the expression of known prognostic factors. Results CEBPE expression was an independent predictor for both overall survival (OS) and event-free survival (EFS) of AML patients. Moreover, low-expression of CEBPE was found to be associated with high relapse rate. We also proved that differential expression of CEBPE stratified the wild-type patients of multiple genes into good and poor outcomes. In addition, the results showed that no obvious improvement was achieved by allogeneic transplantation in CEBPE high-expressed group, while the survival rate (both OS and EFS) was significantly increased in transplanted patients that with low expression of CEBPE. Finally, we found that CEBPE might regulate the expression of known prognostic factors by localizing on their promoters. Conclusion Our findings indicated that CEBPE expression was an independent prognostic factor for AML survival, relapse and allogeneic transplantation, which will provide useful information for outcome prediction and therapeutic decisions.

Enhancers are short DNA segments that play a key role in biological processes, such as accelerating transcription of target genes. Since the enhancer resides anywhere in a genome sequence, it is difficult to precisely identify enhancers. We presented a bi-directional long-short term memory (Bi-LSTM) and attention-based deep learning method (Enhancer-LSTMAtt) for enhancer recognition. Enhancer-LSTMAtt is an end-to-end deep learning model that consists mainly of deep residual neural network, Bi-LSTM, and feed-forward attention. We extensively compared the Enhancer-LSTMAtt with 19 state-of-the-art methods by 5-fold cross validation, 10-fold cross validation and independent test. Enhancer-LSTMAtt achieved competitive performances, especially in the independent test. We realized Enhancer-LSTMAtt into a user-friendly web application. Enhancer-LSTMAtt is applicable not only to recognizing enhancers, but also to distinguishing strong enhancer from weak enhancers. Enhancer-LSTMAtt is believed to become a promising tool for identifying enhancers.

Pyrazinamide (PZA) is the first-line drug commonly used in treating Mycobacterium tuberculosis (Mtb) infections and reduces treatment time by 33%. This prodrug is activated and converted to an active form, Pyrazinoic acid (POA), by Pyrazinamidase (PZase) enzyme. Mtb resistance to PZA is the outcome of mutations frequently reported in pncA, rpsA, and panD genes. Among the mentioned genes, pncA mutations contribute to 72–99% of the total resistance to PZA. Thus, considering the vital importance of this gene in PZA resistance, its frequent mutations (D49N, Y64S, W68G, and F94A) were investigated through in-depth computational techniques to put conclusions that might be useful for new scaffolds design or structure optimization to improve the efficacy of the available drugs. Mutants and wild type PZase were used in extensive and long-run molecular dynamics simulations in triplicate to disclose the resistance mechanism induced by the above-mentioned point mutations. Our analysis suggests that these mutations alter the internal dynamics of PZase and hinder the correct orientation of PZA to the enzyme. Consequently, the PZA has a low binding energy score with the mutants compared with the wild type PZase. These mutations were also reported to affect the binding of Fe2+ ion and its coordinated residues. Conformational dynamics also revealed that β-strand two is flipped, which is significant in Fe2+ binding. MM-GBSA analysis confirmed that these mutations significantly decreased the binding of PZA. In conclusion, these mutations cause conformation alterations and deformities that lead to PZA resistance.

Background: Coronavirus (CoV) is an emerging human pathogen causing severe acute respiratory syndrome (SARS) around the world. Earlier identification of biomarkers for SARS can facilitate detection and reduce the mortality rate of the disease. Thus, by integrated network analysis and structural modeling approach, we aimed to explore the potential drug targets and the candidate drugs for coronavirus medicated SARS. Methods: Differentially expression (DE) analysis of CoV infected host genes (HGs) expression profiles was conducted by using the Limma. Highly integrated DE-CoV-HGs were selected to construct the protein-protein interaction (PPI) network. Results: Using the Walktrap algorithm highly interconnected modules include module 1 (202 nodes); module 2 (126 nodes) and module 3 (121 nodes) modules were retrieved from the PPI network. MYC, HDAC9, NCOA3, CEBPB, VEGFA, BCL3, SMAD3, SMURF1, KLHL12, CBL, ERBB4, and CRKL were identified as potential drug targets (PDTs), which are highly expressed in the human respiratory system after CoV infection. Functional terms growth factor receptor binding, c-type lectin receptor signaling, interleukin-1 mediated signaling, TAP dependent antigen processing and presentation of peptide antigen via MHC class I, stimulatory T cell receptor signaling, and innate immune response signaling pathways, signal transduction and cytokine immune signaling pathways were enriched in the modules. Protein-protein docking results demonstrated the strong binding affinity (-314.57 kcal/mol) of the ERBB4-3cLpro complex which was selected as a drug target. In addition, molecular dynamics simulations indicated the structural stability and flexibility of the ERBB4-3cLpro complex. Further, Wortmannin was proposed as a candidate drug to ERBB4 to control SARS-CoV-2 pathogenesis through inhibit receptor tyrosine kinase-dependent macropinocytosis, MAPK signaling, and NF-kb singling pathways that regulate host cell entry, replication, and modulation of the host immune system. Conclusion: We conclude that CoV drug target “ERBB4” and candidate drug “Wortmannin” provide insights on the possible personalized therapeutics for emerging COVID-19.

Background: Coronavirus (CoV) is an emerging human pathogen causing severe acute respiratory syndrome (SARS) around the world. Earlier identification of biomarkers for SARS can facilitate detection and reduce the mortality rate of the disease. Thus, by integrated network analysis and structural modeling approach, we aimed to explore the potential drug targets and the candidate drugs for coronavirus medicated SARS. Methods: Differentially expression (DE) analysis of CoV infected host genes (HGs) expression profiles was conducted by using the Limma. Highly integrated DE-CoV-HGs were selected to construct the protein-protein interaction (PPI) network. Results: Using the Walktrap algorithm highly interconnected modules include module 1 (202 nodes); module 2 (126 nodes) and module 3 (121 nodes) modules were retrieved from the PPI network. MYC, HDAC9, NCOA3, CEBPB, VEGFA, BCL3, SMAD3, SMURF1, KLHL12, CBL, ERBB4, and CRKL were identified as potential drug targets (PDTs), which are highly expressed in the human respiratory system after CoV infection. Functional terms growth factor receptor binding, c-type lectin receptor signaling, interleukin-1 mediated signaling, TAP dependent antigen processing and presentation of peptide antigen via MHC class I, stimulatory T cell receptor signaling, and innate immune response signaling pathways, signal transduction and cytokine immune signaling pathways were enriched in the modules. Protein-protein docking results demonstrated the strong binding affinity (-314.57 kcal/mol) of the ERBB4-3cLpro complex which was selected as a drug target. In addition, molecular dynamics simulations indicated the structural stability and flexibility of the ERBB4-3cLpro complex. Further, Wortmannin was proposed as a candidate drug to ERBB4 to control SARS-CoV-2 pathogenesis through inhibit receptor tyrosine kinase-dependent macropinocytosis, MAPK signaling, and NF-kb singling pathways that regulate host cell entry, replication, and modulation of the host immune system. Conclusion: We conclude that CoV drug target “ERBB4” and candidate drug “Wortmannin” provide insights on the possible personalized therapeutics for emerging COVID-19.

GPR (G protein receptor) 139 and 142 are novel foundling GPCRs (G protein-coupled receptors) in the class “A” of the GPCRs family and are suitable targets for various biological conditions. To engage these targets, validated pharmacophores and 3D QSAR (Quantitative structure-activity relationship) models are widely used because of their direct fingerprinting capability of the target and an overall accuracy. The current work initially analyzes GPR139 and GPR142 for its genomic alteration via tumor samples. Next to that, the pharmacophore is developed to scan the 3D database for such compounds that can lead to potential agonists. As a result, several compounds have been considered, showing satisfactory performance and a strong association with the target. Additionally, it is gripping to know that the obtained compounds were observed to be responsible for triggering pan-cancer. This suggests the possible role of novel GPR139 and GPR142 as the substances for initiating a physiological response to handle the condition incurred as a result of cancer.

One of the most well-known cancer subtypes worldwide is triple-negative breast cancer (TNBC) which has reduced prediction due to its antagonistic biotic actions and target’s deficiency for the treatment. The current work aims to discover the countenance outlines and possible roles of lncRNAs in the TNBC via computational approaches. Long non-coding RNAs (lncRNAs) exert profound biological functions and are widely applied as prognostic features in cancer. We aim to identify a prognostic lncRNA signature for the TNBC. First, samples were filtered out with inadequate tumor purity and retrieved the lncRNA expression data stored in the TANRIC catalog. TNBC sufferers were divided into two prognostic classes which were dependent on their survival time (shorter or longer than 3 years). Random forest was utilized to select lncRNA features based on the lncRNAs differential expression between shorter and longer groups. The Stochastic gradient boosting method was used to construct the predictive model. As a whole, 353 lncRNAs were differentially transcribed amongst the shorter and longer groups. Using the recursive feature elimination, two lncRNAs were further selected. Trained by stochastic gradient boosting, we reached the highest accuracy of 69.69% and area under the curve of 0.6475. Our findings showed that the two-lncRNA signs can be proved as potential biomarkers for the prognostic grouping of TNBC’s sufferers. Many lncRNAs remained dysregulated in TNBC, while most of them are likely play a role in cancer biology. Some of these lncRNAs were linked to TNBC’s prediction, which makes them likely to be promising biomarkers.

Accurate and fast characterization of the subtype sequences of Avian influenza A virus (AIAV) hemagglutinin (HA) and neuraminidase (NA) depends on expanding diagnostic services and is embedded in molecular epidemiological studies. A new approach for classifying the AIAV sequences of the HA and NA genes into subtypes using DNA sequence data and physicochemical properties is proposed. This method simply requires unaligned, full-length, or partial sequences of HA or NA DNA as input. It allows for quick and highly accurate assignments of HA sequences to subtypes H1–H16 and NA sequences to subtypes N1–N9. For feature extraction, k-gram, discrete wavelet transformation, and multivariate mutual information were used, and different classifiers were trained for prediction. Four different classifiers, Naïve Bayes, Support Vector Machine (SVM), K nearest neighbor (KNN), and Decision Tree, were compared using our feature selection method. This comparison is based on the 30% dataset separated from the original dataset for testing purposes. Among the four classifiers, Decision Tree was the best, and Precision, Recall, F1 score, and Accuracy were 0.9514, 0.9535, 0.9524, and 0.9571, respectively. Decision Tree had considerable improvements over the other three classifiers using our method. Results show that the proposed feature selection method, when trained with a Decision Tree classifier, gives the best results for accurate prediction of the AIAV subtype.

The continued emergence of human coronaviruses (hCoVs) in the last few decades has posed an alarming situation and requires advanced cross-protective strategies against these pandemic viruses. Among these, Middle East Respiratory Syndrome coronavirus (MERS-CoV), Severe Acute Respiratory Syndrome coronavirus (SARS-CoV), and Severe Acute Respiratory Syndrome coronavirus-2 (SARS-CoV-2) have been highly associated with lethality in humans. Despite the challenges posed by these viruses, it is imperative to develop effective antiviral therapeutics and vaccines for these human-infecting viruses. The proteomic similarity between the receptor-binding domains (RBDs) among the three viral species offers a potential target for advanced cross-protective vaccine designs. In this study, putative immunogenic epitopes including Cytotoxic T Lymphocytes (CTLs), Helper T Lymphocytes (HTLs), and Beta-cells (B-cells) were predicted for each RBD-containing region of the three highly pathogenic hCoVs. This was followed by the structural organization of peptide- and mRNA-based prophylactic vaccine designs. The validated 3D structures of these epitope-based vaccine designs were subjected to molecular docking with human TLR4. Furthermore, the CTL and HTL epitopes were processed for binding with respective human Lymphocytes Antigens (HLAs). In silico cloning designs were obtained for the prophylactic vaccine designs and may be useful in further experimental designs. Additionally, the epitope-based vaccine designs were evaluated for immunogenic activity through immune simulation. Further studies may clarify the safety and efficacy of these prophylactic vaccine designs through experimental testing against these human-pathogenic coronaviruses.

Inspired by the recent cocrystallization and theory of energetic materials, we theoretically investigated the intermolecular vibrational energy transfer process and the non-covalent intermolecular interactions between explosive compounds. The intermolecular interactions between 2,4,6-trinitrotoluene (TNT) and 2,4,6,8,10,12-hexanitro-2,4,6,8,10,12-hexaazaisowurtzitane (CL-20) and between 1,3,5,7-tetranitro-1,3,5,7-tetrazocane (HMX) and CL-20 were studied using calculated two-dimensional infrared (2D IR) spectra and the independent gradient model based on the Hirshfeld partition (IGMH) method, respectively. Based on the comparison of the theoretical infrared spectra and optimized geometries with experimental results, the theoretical models can effectively reproduce the experimental geometries. By analyzing cross-peaks in the 2D IR spectra of TNT/CL-20, the intermolecular vibrational energy transfer process between TNT and CL-20 was calculated, and the conclusion was made that the vibrational energy transfer process between CL-20 and TNTII (TNTIII) is relatively slower than between CL-20 and TNTI. As the vibration energy transfer is the bridge of the intermolecular interactions, the weak intermolecular interactions were visualized using the IGMH method, and the results demonstrate that the intermolecular non-covalent interactions of TNT/CL-20 include van der Waals (vdW) interactions and hydrogen bonds, while the intermolecular non-covalent interactions of HMX/CL-20 are mainly comprised of vdW interactions. Further, we determined that the intermolecular interaction can stabilize the trigger bond in TNT/CL-20 and HMX/CL-20 based on Mayer bond order density, and stronger intermolecular interactions generally indicate lower impact sensitivity of energetic materials. We believe that the results obtained in this work are important for a better understanding of the cocrystal mechanism and its application in the field of energetic materials.

Meiotic recombination is the driving force of evolutionary development and an important source of genetic variation. The meiotic recombination does not take place randomly in a chromosome but occurs in some regions of the chromosome. A region in chromosomes with higher rate of meiotic recombination events are considered as hotspots and a region where frequencies of the recombination events are lower are called coldspots. Prediction of meiotic recombination spots provides useful information about the basic functionality of inheritance and genome diversity. This study proposes an intelligent computational predictor called iRSpots-DNN for the identification of recombination spots. The proposed predictor is based on a novel feature extraction method and an optimized deep neural network (DNN). The DNN was employed as a classification engine whereas, the novel features extraction method was developed to extract meaningful features for the identification of hotspots and coldspots across the yeast genome. Unlike previous algorithms, the proposed feature extraction avoids bias among different selected features and preserved the sequence discriminant properties along with the sequence-structure information simultaneously. This study also considered other effective classifiers named support vector machine (SVM), K-nearest neighbor (KNN), and random forest (RF) to predict recombination spots. Experimental results on a benchmark dataset with 10-fold cross-validation showed that iRSpots-DNN achieved the highest accuracy, i.e., 95.81%. Additionally, the performance of the proposed iRSpots-DNN is significantly better than the existing predictors on a benchmark dataset. The relevant benchmark dataset and source code are freely available at: https://github.com/Fatima-Khan12/iRspot_DNN/tree/master/iRspot_DNN.

Type IV secreted effectors (T4SEs) can be translocated into the cytosol of host cells via type IV secretion system (T4SS) and cause diseases. However, experimental approaches to identify T4SEs are time- and resource-consuming, and the existing computational tools based on machine learning techniques have some obvious limitations such as the lack of interpretability in the prediction models. In this study, we proposed a new model, T4SE-XGB, which uses the eXtreme gradient boosting (XGBoost) algorithm for accurate identification of type IV effectors based on optimal features based on protein sequences. After trying 20 different types of features, the best performance was achieved when all features were fed into XGBoost by the 5-fold cross validation in comparison with other machine learning methods. Then, the ReliefF algorithm was adopted to get the optimal feature set on our dataset, which further improved the model performance. T4SE-XGB exhibited highest predictive performance on the independent test set and outperformed other published prediction tools. Furthermore, the SHAP method was used to interpret the contribution of features to model predictions. The identification of key features can contribute to improved understanding of multifactorial contributors to host-pathogen interactions and bacterial pathogenesis. In addition to type IV effector prediction, we believe that the proposed framework can provide instructive guidance for similar studies to construct prediction methods on related biological problems. The data and source code of this study can be freely accessed at https://github.com/CT001002/T4SE-XGB.

This study aims to achieve a clearer and stronger understanding of all the mechanisms involved in the occurrence as well as in the progression of lung cancer along with discovering trustworthy prognostic markers. We combined four gene expression profiles (GSE19188, GSE19804, GSE101929, and GSE18842) from the GEO database and screened the commonly differentially expressed genes (CDEGs). We performed differentially expressed group analysis on CDEGs, alteration and mutational analysis, and expression level verification of core differential genes. Systems biology discoveries in our examination are predictable with past reports. Curiously, our examination revealed that screened biomarker adjustments, for the most part, coexist in lung cancer. After screening 952 CDEGs, we found that the up-regulation of neuromedin U (NMU) and GTSE1 in the case of lung cancer is related to poor prognosis. On the other hand, FOS CDKN1C expression is associated with poor prognosis and is responsible for the down-regulation of CDKN1C and FOS. Changes in these qualities are on free pathways to lung cancer and are not usually of combined quality variety. Even though biomarkers were related to both survival occasions in our examination, it gives us another point of view while playing out the investigation of hereditary changes and clinical highlights employing information mining. Based on our results, we found potential and prospective clinical applications in GTSE1, NMU, FOS, and CDKN1C to act as prognostic markers in case of lung cancer.

LncRNAs, miRNAs, mRNAs, methylation, and proteins exert profound biological functions and are widely applied as prognostic features in liver cancer. This study aims to identify prognostic biomarkers’ signature for liver cancer. Samples with inadequate tumor purity were filtered out and the expression data from different resources were retrieved. The Spares learning approach was applied to select lncRNAs, miRNAs, mRNAs, methylation, and proteins’ features based on their differentially expressed groups. The LASSO boosting technique was employed for the predictive model construction. A total of 200 lncRNAs, 200 miRNAs, 371 mRNAs, 371 methylations, and 184 proteins were observed to be differentially expressed. Five lncRNAs, 11 miRNAs, 30 mRNAs, 4 methylations, and 3 proteins were selected for further evaluation using the feature elimination technique. The highest accuracy of 89.32% is achieved as a result of training and learning by Spares learning methodology. Final outcomes revealed that 5 lncRNA, 11 miRNA, 30 mRNA, 4 methylation, and 3 protein signatures could be potential biomarkers for the prognosis of liver cancer patients.

A central approach for better understanding the forces involved in maintaining protein structures is to investigate the protein folding and thermodynamic properties. The effect of the folding process is often disturbed in mutated states. To explore the dynamic properties behind mutations, molecular dynamic (MD) simulations have been widely performed, especially in unveiling the mechanism of drug failure behind mutation. When comparing wild type (WT) and mutants (MTs), the structural changes along with solvation free energy (SFE), and Gibbs free energy (GFE) are calculated after the MD simulation, to measure the effect of mutations on protein structure. Pyrazinamide (PZA) is one of the first-line drugs, effective against latent Mycobacterium tuberculosis isolates, affecting the global TB control program 2030. Resistance to this drug emerges due to mutations in pncA and rpsA genes, encoding pyrazinamidase (PZase) and ribosomal protein S1 (RpsA) respectively. The question of how the GFE may be a measure of PZase and RpsA stabilities, has been addressed in the current review. The GFE and SFE of MTs have been compared with WT, which were already found to be PZA-resistant. WT structures attained a more stable state in comparison with MTs. The physiological effect of a mutation in PZase and RpsA may be due to the difference in energies. This difference between WT and MTs, depicted through GFE plots, might be useful in predicting the stability and PZA-resistance behind mutation. This study provides useful information for better management of drug resistance, to control the global TB problem.

Nanotechnology utilizes engineered materials and devices which function with biological systems at the molecular level and could transform the management of neurodegenerative diseases (NDs) by provoking, reacting to, and intermingling with target sites to stimulate physiological responses while minimizing side effects. Blood–brain barrier (BBB) protects the brain from harmful agents, and transporting drugs across the BBB is a major challenge for diagnosis, targeting, and treatment of NDs. The BBB provides severe limitations for diagnosis and treatment of Alzheimer’s disease (AD), Parkinson’s disease (PD), and various other neurological diseases. Conventional drug delivery systems generally fail to cross the BBB, thus are inefficient in treatment. Although gradual development through research is ensuring the progress of nanotheranostic approaches from animal to human modeling, aspects of translational applicability and safety are a key concern. This demands a deep understanding of the interaction of body systems with nanomaterials. There are various plant-based nanobioactive compounds which are reported to have applicability in the diagnosis and treatment of these NDs. This review article provides an overview of applications of nanotheranostics in AD and PD. The review also discusses nano-enabled drug delivery systems and their current and potential applications for the treatment of various NDs.