Mechanistic insights into the conformational changes and alterations in residual communications due to the mutations in the pncA Gene of Mycobacterium tuberculosis: A computational perspective for effective therapeutic solutions

Due to its emerging resistance to first-line anti-TB medications, tuberculosis (TB) is one of the most contagious illness in the world. According to reports, the effectiveness of treating TB is severely impacted by drug resistance, notably resistance caused by mutations in the pncA gene-encoded pyrazinamidase (PZase) to the front-line drug pyrazinamide (PZA). The present study focused on investigating the resistance mechanism caused by the mutations D12N, T47A, and H137R to better understand the structural and molecular events responsible for the resistance acquired by the pncA gene of Mycobacterium tuberculosis (MTB) at the structural level. Bioinformatics analysis predicted that all three mutations were deleterious and located near the active centre of the pncA, affecting its functional activity. Furthermore, molecular dynamics simulation (MDS) results established that mutations significantly reduced the structural stability and caused the rearrangement of FE2+ in the active centre of pncA. Moreover, essential dynamics analysis, including principal component analysis (PCA) and free energy landscape (FEL), concluded variations in the protein motion and decreased conformational space in the mutants. Additionally, the mutations potentially impacted the network topologies and altered the residual communications in the network. The complex simulation study results established the significant movement of the flap region from the active centre of mutant complexes, further supporting the flap region's significance in developing resistance to the PZA drug. This study advances our knowledge of the primary cause of the mechanism of PZA resistance and the structural dynamics of pncA mutants, which will help us to design new and potent chemical scaffolds to treat drug-resistant TB (DR-TB).

In silico approaches for the identification of potential allergens among hypothetical proteins from Alternaria alternata and its functional annotation

Direct exposure to the fungal species Alternaria alternata is a major risk factor for the development of asthma, allergic rhinitis, and inflammation. As of November 23rd 2020, the NCBI protein database showed 11,227 proteins from A. alternata genome as hypothetical proteins (HPs). Allergens are the main causative of several life-threatening diseases, especially in fungal infections. Therefore, the main aim of the study is to identify the potentially allergenic inducible proteins from the HPs in A. alternata and their associated functional assignment for the complete understanding of the complex biological systems at the molecular level. AlgPred and Structural Database of Allergenic Proteins (SDAP) were used for the prediction of potential allergens from the HPs of A. alternata. While analyzing the proteome data, 29 potential allergens were predicted by AlgPred and further screening in SDAP confirmed the allergic response of 10 proteins. Extensive bioinformatics tools including protein family classification, sequence-function relationship, protein motif discovery, pathway interactions, and intrinsic features from the amino acid sequence were used to successfully predict the probable functions of the 10 HPs. The functions of the HPs are characterized as chitin-binding, ribosomal protein P1, thaumatin, glycosyl hydrolase, and NOB1 proteins. The subcellular localization and signal peptide prediction of these 10 proteins has further provided additional information on localization and function. The allergens prediction and functional annotation of the 10 proteins may facilitate a better understanding of the allergenic mechanism of A. alternata in asthma and other diseases. The functional domain level insights and predicted structural features of the allergenic proteins help to understand the pathogenesis and host immune tolerance. The outcomes of the study would aid in the development of specific drugs to combat A. alternata infections.

Computational insights into potential marine natural products as selective inhibitors of Mycobacterium tuberculosis InhA: A structure-based virtual screening study

Several factors are associated with the emergence of drug resistance mechanisms, such as impermeable cell walls, gene mutations, and drug efflux systems. Consequently, bacteria acquire resistance, leading to a decrease in drug efficacy. A new and innovative strategy is required to combat drug resistance in tuberculosis (TB) effectively. Therefore, targeting the mycolic acid biosynthesis pathway, which is involved in synthesising mycolic acids (MAs), essential structural components responsible for mycobacterial pathogenicity, has garnered interest in TB research and the concept of drug resistance. In this context, InhA, which plays a crucial role in the fatty acid synthase-II (FAS-II) system of the MA biosynthetic pathway, was selected as a druggable target for screening investigation. To identify potential lead molecules against InhA, diverse marine natural products (MNPs) were collected from the comprehensive marine natural products database (CMNPD). Virtual screening studies aided in selecting potential lead molecules that best fit within the substrate-binding pocket (SBP) of InhA, forming crucial hydrogen bond interaction with the catalytic residue Tyr158. Three MNPs, CMNPD30814, CMNPD1702, and CMNPD27355, were chosen as prospective alternative molecules due to their favorable pharmacokinetic properties and lack of toxicity according to ProTox-II predictions. Additionally, improved reactivity of the MNPs was observed in the results of density functional theory (DFT) studies. Furthermore, comparative molecular dynamics simulation (MDS), principal component (PC)-based free energy landscape (FEL) analysis, and molecular mechanics Poisson-Boltzmann surface area (MM-PBSA) were employed to show enhanced structural stability, increased H-bond potential, and high binding affinity toward the target InhA. Moreover, the hot spot residues that contributed to the high binding energy profile and anchored the stability of the complexes were revealed with their individual interaction energy. The computational insights from this study provide potential avenues to combat TB through the multifaceted mode of action of these marine lead molecules, which can be further explored in future experimental investigations.

Exploring Nocardia's ecological spectrum and novel therapeutic frontiers through whole-genome sequencing: unraveling drug resistance and virulence factors

Nocardia farcinica is the leading pathogen responsible for nocardiosis, a life-threatening infection primarily affecting immunocompromised patients. In this study, the genomic sequence of a clinically isolated N. farcinica sample was sequenced. Subsequently, the assembled genome was annotated to identify antimicrobial resistance and virulence genes, as well as plasmid and prophages. The analysis of the entire genome size was 6,021,225 bp, with a GC content of 70.78% and consists of 103 contigs and N50 values of 292,531 bp. The genome analysis revealed the presence of several antimicrobial resistance genes, including RbpA, mtrA, FAR-1, blaFAR-1, blaFAR-1_1, and rox. In addition, virulence genes such as relA, icl, and mbtH were also detected. The present study signifies that N. farcinica genome is pivotal for the understanding of antimicrobial resistance and virulence genes is crucial for comprehending resistance mechanism, and developing effective strategies to combat bacterial infections effectively, especially adhesins and toxins. This study aids in identifying crucial drug targets for combating multidrug-resistant N. farcinica in the future.

A novel search engine for proteins involved in Notch crosstalk signaling pathways

To regulate biological activity in humans, the Notch signaling pathway (NSP) plays an essential role in a wide array of cellular development and differentiation process. In recent years, many studies have reported that aberrant activation of Notch is associated with the tumor process; but no appropriate database exists to fill this significant gap. To address this, we created a pioneering database NCSp, which is open access and comprises intercommunicating pathways and related protein mutations. This allows scientists to understand better the cause of single amino acid mutations in proteins. Therefore, NCSp provides information on the predicted functional effect of human protein mutations, which aids in understanding the importance of mutations linked to the Notch crosstalk signaling pathways in cancerous and non-cancerous systems. This database might be helpful for therapeutic mutation analysis, molecular biology, and structural biology researchers. The NCSp database can be accessed through https://bioserver3.physics.iisc.ac.in/cgi-bin/nccspd/.

Multilayer precision-based screening of potential inhibitors targeting Mycobacterium tuberculosis acetate kinase using in silico approaches

Tuberculosis (TB), caused by Mycobacterium tuberculosis (MTB), remains a major threat to global health, with the emergence of multi-drug and extensively drug-resistant strains posing a serious challenge. Thereby, understanding the molecular basis of MTB virulence and disease pathogenesis is critical for developing effective therapeutic strategies. Targeting proteins involved in central metabolism has been recognized as a promising therapeutic approach to combat MTB. In this regard, the enzyme AckA of the acetate metabolic pathway which produces acetate from acetyl phosphate, is an important drug target for various pathogenic organisms. Therefore, this study aimed to identify potential AckA inhibitors through in silico methods, including molecular modeling, molecular dynamics simulation (MDS), and high-throughput virtual screening (HTVS) followed by ADMETox, MMGBSA, Density Functional Theory (DFT) calculations. HTVS of one million compounds from the ZINC database against AckA resulted in the top five hits (ZINC82048449, ZINC1219737510, ZINC1771921358, ZINC119699567, and ZINC1427100376) with better binding affinity and optimal binding free energy. MDS studies on complexes revealed that key residues, Asn195, Asp266, Phe267, Gly314, and Asn318 played a significant role in stable interactions of the top-ranked compounds to AckA. These outcomes provide insights into the optimal binding of the leads to inhibit the acetate pathway and aid in the rational design of novel therapeutic agents. Thus, the identified leads may act as promising compounds for targeting AckA and may serve as a potential therapeutic modality for treating TB. Our findings offer valuable insights into the inhibition of the acetate pathway, while also serving as a blueprint for rational drug design. The identified leads hold promise as compelling compounds for targeting AckA, thereby offering a potential therapeutic avenue for tackling TB. Thus, our study uncovers a pathway toward promising TB therapeutics by elucidating AckA inhibitors. By leveraging in silico methodologies, potent compounds that hold the potential to thwart AckA's role in MTB's acetate pathway have been unveiled. This breakthrough fosters optimism in the quest for novel and effective TB treatments, addressing a global health challenge with renewed vigor.

Differential gene expression analysis combined with molecular dynamics simulation study to elucidate the novel potential biomarker involved in pulmonary TB

Tuberculosis (TB) is a lethal multisystem disease that attacks the lungs' first line of defense. A substantial threat to public health and a primary cause of death is pulmonary TB. This study aimed to identify and investigate the probable differentially expressed genes (DEGs) primarily involved in Pulmonary TB. Accordingly, three independent gene expression data sets, numbered GSE139825, GSE139871, and GSE54992, were utilized for this purpose. The identified DEGs were used for bioinformatics-based analysis, including physical gene interaction, Gene Ontology (GO), network analysis and pathway studies using the Kyoto Encyclopedia of Genes and Genomes pathway (KEGG). The computational analysis predicted that TNFAIP6 is the significant DEG in the gene expression profiling of TB datasets. According to gene ontology analysis, TNFAIP6 is also essential in injury and inflammation. Further, TNFA1P6 is strongly linked to arsenic poisoning, evident from the results of NetworkAnalyst, a comprehensive and interactive platform for gene expression profiling via network visual analytics. As a result, the TNFAIP6 gene was ultimately chosen as a candidate DEG and subsequently employed for in silico structural characterization studies. The tertiary structure of TNFAIP6 was modelled using the ROBETTA server, followed by validation with SAVES and ProSA webserver. Additionally, structural dynamic studies, including molecular dynamics simulation (MDS) and essential dynamics analysis, including principal component (PC) based free energy landscape (FEL) analysis, was used for checking the stability of TNFAIP6 models. The dynamics result established the structural rigidity of modelled TNFAIP6 through RMSD, RMSF and RoG results. The FEL analysis revealed the restricted conformational flexibility of TNFAIP6 by displaying a single minimum energy basin in the contour plot. The comprehensive computational analysis established that TNFAIP6 could serve as a viable biomarker to assess the severity of pulmonary TB.

Regenerative marine waste towards CaCO3 nanoformulation for Alzheimer's therapy

Alzheimer's disorder (AD) is associated with behavioural and cognitive destruction with due respect to the neurological degeneration. Conventional therapeutic approach for treatment of AD using neuroprotective drugs suffered certain limitations such as poor solubility, insufficient bioavailability, adverse side effects at higher dose and ineffective permeability on blood brain barrier (BBB). Development of nanomaterial based drug delivery system helped to overcome these barriers. Hence the present work focused on encapsulating neuroprotective drug citronellyl acetate within CaCO₃ nanoparticles to develop neuroprotective CaCO₃ nanoformulation (CA@CaCO₃ NFs). CaCO₃ was derived from marine conch shell waste, while the neuroprotective drug citronellyl acetate was scrutinized by in-silico high throughput screening. In-vitro findings revealed that CA@CaCO₃ nanoformulation exhibited enhanced free radical scavenging activity of 92% (IC₅₀ value - 29.27 ± 2.6 μg/ml), AChE inhibition of 95% (IC₅₀ value - 25.6292 ± 1.5 μg/ml) at its maximum dose (100 μg/ml). CA@CaCO₃ NFs attenuated the aggregation of β-amyloid peptide (Aβ) and also disaggregated the preformed mature plaques the major risk factor for AD. Overall, the present study reveals that CaCO₃ nanoformulations exhibits potent neuroprotective potential when compared to the CaCO₃ nanoparticles alone and citronellyl acetate alone due to the sustained drug release and synergistic effect of CaCO₃ nanoparticles and citronellyl acetate depicting the fact that CaCO₃ can act as promising drug delivery system for treatment of neurodegenerative and CNS related disorders.

Computational exploration of molecular flexibility and interaction of meropenem analogs with the active site of oxacillinase-23 in Acinetobacter baumannii

Background: Carbapenem-resistant Acinetobacter baumannii is an opportunistic pathogen responsible for nosocomial infections and is one of the biggest global threats according to the World Health Organization (WHO), particularly causing substantial morbidity and mortality. Objectives: This study aimed at using computational approaches to screen meropenem and its analogs against OXA-23-positive Acinetobacter baumannii, analyzing the correlations between kinetic and phenotypic characteristics. Methods: A total of 5,450 compounds were screened using virtual screening workflow (HTVS, Glide-SP, and Glide-XP) to identify the best compounds based on their binding energy and interactions against OXA-23 and OXA-27 as they had phenotypic data available. Molecular dynamics simulation and density functional theory (DFT) studies were performed from the outcome of molecular docking analysis. Results: During simulations, meropenem and its analogs exhibited high-level stable interactions with Ser79, Ser126, Thr217, Trp219, and Arg259 of OXA-23. Meropenem displayed a CovDock energy of about -3.5 and -1.9 kcal mol^-1 against OXA-23 and OXA-27, respectively. Among the 5,450 compounds, Pubchem_10645796, Pubchem_25224737, and ChEMBL_14 recorded CovDock energy between -6.0 and -9.0 kcal mol^-1. Moreover, the infra-red (IR) spectrophotometric analysis revealed C=O and C-N atoms showing bands at 1,800 and 1,125 cm^-1, respectively. These observed data are in congruence with the experimental observations. Conclusion: The identified compounds showed good agreement with the spectrophotometric analysis using DFT methods. In the earlier studies, meropenem's MIC value was 32 μg mL^-1 in OXA-23-positive isolate A2265 compared to the MIC of 1 μg mL^-1 in Δbla _OXA-23 A2265. Comparing the CovDock energy and hydrogen-bonding interactions, the predicted results are in good agreement with the experimental data reported earlier. Our results highlight the importance of OXA-23 molecular docking studies and their compliance with the phenotypic results. It will help further in developing newer antibiotics for treating severe infections associated with carbapenem-resistant A. baumannii.

Computational investigation of the increased virulence and pathogenesis of SARS-CoV-2 lineage B.1.1.7

New variants of SARS-CoV-2 are being reported worldwide. The World Health Organization has reported Alpha (B.1.1.7), Beta (B.1.351), Gamma (P.1), Delta (B.1.617.2) and Omicron (B.1.1.529) as the variants of concern. There are speculations that the variants might evade the host immune responses induced by currently available vaccines and develop resistance to drugs under consideration. The first step of viral infection in COVID-19 occurs through the interaction of the spike protein's receptor-binding domain (RBD) with the peptidase domain of the human ACE-2 (hACE-2) receptor. This study aims to get a molecular-level understanding of the mechanism behind the increased infection rate in the alpha variant. We have computationally studied the spike protein interaction in both the wild-type and B.1.1.7 variant with the hACE-2 receptor using molecular dynamics and MM-GBSA based binding free energy calculations. The binding free energy difference shows that the mutant variant of the spike protein has increased binding affinity for the hACE-2 receptor (i.e. ΔG(N501Y,A570D) is in the range -7.2 to -7.6 kcal mol^-1) and the results were validated using Density functional theory. We demonstrate that with the use of state-of-the-art computational approaches, we can, in advance, predict the virulent nature of variants of SARS-CoV-2 and alert the world healthcare system.

Exploring the druggability of oxidized low-density lipoprotein (ox-LDL) receptor, LOX-1, a proatherogenic drug target involved in atherosclerosis

Lectin-like oxidized low-density lipoprotein (ox-LDL) receptor 1 (LOX-1) is a vital scavenger receptor involved in ox-LDL binding, internalization, and subsequent proatherogenic signaling leading to cellular dysfunction and atherosclerotic plaque formation. Existing data suggest that modulation of ox-LDL - LOX-1 interaction can prevent or slow down atherosclerosis. Therefore, we utilized computational methods such as multi-solvent simulation and characterized two top-ranked druggable sites. Using systematic molecular docking followed by atomistic molecular dynamics simulation, we have identified and shortlisted small molecules from the NCI library that target two key binding sites. We demonstrate, using surface plasmon resonance (SPR), that four of the shortlisted molecules bind one-on-one to the purified C-terminal domain (CTLD) of LOX-1 receptor with high affinity (K_D), ranging from 4.9 nM to 20.1 μM. Further, we performed WaterMap analysis to understand the role of individual water molecules in small molecule binding and the LOX-1-ligand complex stability. Our data clearly show that LOX-1 is druggable with small molecules. Our study provides strategies to identify novel inhibitors to attenuate ox-LDL - LOX-1 interaction.

Functional significance of mouse seminal vesicle sulfhydryl oxidase on sperm capacitation in vitro

During ejaculation, cauda epididymal spermatozoa are suspended in a protein-rich solution of seminal plasma which is composed of proteins mostly secreted from the seminal vesicle. These seminal proteins interact with the sperm cells and bring about changes in their physiology, so that they can become capacitated in order for the fertilization to take place. Sulfhydryl oxidase (SOX) is a member of the QSOX family and its expression is found to be high in the seminal vesicle secretion of mouse. Previously, it has been reported to cross-link thiol containing amino acids among major seminal vesicle secretion (SVS) proteins. However, its role in male reproduction is unclear. In this study, we determined the role of SOX on epididymal sperm maturation and also disclosed the binding effect of SOX on the sperm fertilizing ability in vitro. In order to achieve the above two objectives, we constructed a Sox clone (1.7 kb) using a pET-30a vector. His-tagged recombinant Sox was over expressed in Shuffle Escherichia coli cells and purified using His-Trap column affinity chromatography along with hydrophobic interaction chromatography. The purified SOX was confirmed by Western blot analysis and by its activity with DTT as a substrate. Results obtained from immunocytochemical staining clearly indicated that SOX possesses a binding site on the sperm acrosome. The influence of SOX on oxidation of sperm sulfhydryl to disulfides during epididymal sperm maturation was evaluated by a thiol labelling agent, mBBr. The SOX protein binds on to the sperm cells and increases their progressive motility. The effect of SOX binding on reducing the [Ca2+]i concentration in sperm head, was determined using a calcium probe, Fluo-3 AM. The inhibitory influence of SOX on sperm acrosome reaction was shown by using calcium ionophore A32187 to induce the acrosome reaction. The acrosome-reacted sperm were examined by staining with FITC-conjugated Arachis hypogaea (peanut) lectin. Furthermore, immunocytochemical analysis revealed that SOX remains bound to the sperm cells in the uterus but disappears in the oviduct during their transit in the female reproductive tract. The results from the above experiment revealed that SOX binding on to the sperm acrosome prevents sperm capacitation by affecting the [Ca2+]i concentration in the sperm head and the ionophore-induced acrosome reaction. Thus, the binding of SOX on to the sperm acrosome may possibly serve as a decapacitation factor in the uterus to prevent premature capacitation and acrosome reaction, thus preserving their fertilizing ability.

Structural Modeling of Drosophila melanogaster Gut Cytochrome P450s and Docking Comparison of Fruit Fly Gut and Human Cytochrome P450s

Drosophila melanogaster is a prominent model organism in developmental biology research and in studies related to pathophysiological conditions like cancer and Alzheimer's disease. The fruit fly gut contains several cytochrome P450s (CYP450s) which have central roles in Drosophila development and in the normal physiology of the gut. Since the crystal structures of these proteins have not been deciphered yet, we modeled the structure of 29 different D. melanogaster gut CYP450s using Prime (Schrödinger). The sequences of chosen D. melanogaster gut CYP450s were compared with that of their human counterparts. The common gut (and liver) microsomal CYP450s in humans were chosen for structural comparison to find the homology and identity % of D. melanogaster CYPs with that of their human counterparts. The modeled structures were validated using PROCHECK and the best fit models were used for docking several known human pharmacological agents/drugs to the modeled D. melanogaster gut CYP450s. Based on the binding affinities (ΔG values) of the selected drug molecules with the modeled fly gut CYPs, the plausible differences in metabolism of the prominent drugs in humans and fly were projected. The gut is involved in absorption of oral drugs/pharmacological agents and hence, upregulation of intestinal CYP450 and their reactions with endobiotics and xenobiotics is envisaged. The insights gleaned from this work can validate D. melanogaster as a model organism for studying intestinal drug metabolism, particularly in the context of a) toxicology of pharmacological agents to the gut cells and b) how gut P450 metabolites/products can influence gut homeostasis. This work can help establish a platform for further in vitro investigations on how intestinal CYP450 metabolism can influence gut health. The data from this work can be used for further in silico studies and this work can serve as a platform for future in vitro investigations on intestinal CYP450-mediated metabolism of endo- and xeno-biotics in D. melanogaster.

Correction: Computational investigation of the increased virulence and pathogenesis of SARS-CoV-2 lineage B.1.1.7

Correction for ‘Computational investigation of the increased virulence and pathogenesis of SARS-CoV-2 lineage B.1.1.7’ by N. Arul Murugan et al., Phys. Chem. Chem. Phys., 2022, 24, 20371–20380, https://doi.org/10.1039/D2CP00469K.

Computational screening of potential inhibitors targeting MurF of Brugia malayi Wolbachia through multi-scale molecular docking, molecular dynamics and MM-GBSA analysis

Lymphatic filariasis is a parasitic disease caused by the worms Wuchereria bancrofti, Brugia malayi and Brugia timori. Three anti-filarial drugs namely Diethylcarbamazine, Ivermectin and Albendazole and their combinations are used as the control strategies for filariasis. The disease has received much attention in drug discovery due to the unavailability of vaccines and the toxic pharmaceutical properties of the existing drugs. In Wolbachia endosymbiont Brugia malayi, the UDP-N-acetylmuramoyl-tripeptide-d-alanyl-d-alanine ligase (MurF) plays a key role in peptidoglycan biosynthesis pathway and therefore can be considered as effective drug target against filariasis disease. Therefore, in the present study, MurF was selected as the therapeutic target to identify specific inhibitors against filariasis. Homology modeling was performed to predict the three-dimensional structure of MurF due to the absence of the experimental structure. Further molecular dynamics simulation and structure-based high throughput virtual screening with three different chemical databases (Zinc, Maybridge and Specs) were carried out to identify potent inhibitors and also to check their conformations inside the binding site of MurF, respectively. Top three compounds with high docking score and high relative binding affinity against MurF were selected. Further, validation studies, including predicted ADME (Absorption, Distribution, Metabolism, Excretion) assessment, binding free energy using MM-GBSA (Molecular Mechanics Generalized Born Surface Area) and DFT (Density Functional Theory) calculations were performed for the top three compounds. From the results, it was observed that all the three compounds were predicted to show high reactivity, acceptable range of pharmacokinetic properties and high binding affinity with the drug target MurF. Overall, the results could provide more understanding on the inhibition of MurF enzyme and the screened compounds could lead to the development of new specific anti-filarial drugs.

WaterMap and Molecular Dynamic Simulation-Guided Discovery of Potential PAK1 Inhibitors Using Repurposing Approaches

p21-Activated kinase 1 (PAK1) is positioned at the nexus of several oncogenic signaling pathways. Currently, there are no approved inhibitors for disabling the transfer of phosphate in the active site directly, as they are limited by lower affinity, and poor kinase selectivity. In this work, a repurposing study utilizing FDA-approved drugs from the DrugBank database was pursued with an initial selection of 27 molecules out of ∼2162 drug molecules, based on their docking energies and molecular interaction patterns. From the molecules that were considered for WaterMap analysis, seven molecules, namely, Mitoxantrone, Labetalol, Acalabrutinib, Sacubitril, Flubendazole, Trazodone, and Niraparib, ascertained the ability to overlap with high-energy hydration sites. Considering many other displaced unfavorable water molecules, only Acalabrutinib, Flubendazole, and Trazodone molecules highlighted their prominence in terms of binding affinity gains through ΔΔG that ranges between 6.44 and 2.59 kcal/mol. Even if Mitoxantrone exhibited the highest docking score and greater interaction strength, it did not comply with the WaterMap and molecular dynamics simulation results. Moreover, detailed MD simulation trajectory analyses suggested that the drug molecules Flubendazole, Niraparib, and Acalabrutinib were highly stable, observed from their RMSD values and consistent interaction pattern with Glu315, Glu345, Leu347, and Asp407 including the hydrophobic interactions maintained in the three replicates. However, the drug molecule Trazodone displayed a loss of crucial interaction with Leu347, which was essential to inhibit the kinase activity of PAK1. The molecular orbital and electrostatic potential analyses elucidated the reactivity and strong complementarity potentials of the drug molecules in the binding pocket of PAK1. Therefore, the CADD-based reposition efforts, reported in this work, helped in the successful identification of new PAK1 inhibitors that requires further investigation by in vitro analysis.

A mouse testis serine protease, TESP1, as the potential SPINK3 receptor protein on mouse sperm acrosome

Serine protease inhibitor Kazal type 3 (SPINK3) from mouse seminal vesicles is a Kazal-type trypsin inhibitor. It has been shown to bind to the sperm acrosome and modify sperm activity by influencing the sub-cellular Ca2+ influx. Previously, SPINK3 was reported to suppress in vitro sperm capacitation. However, under natural coitus, SPINK3 is removed from the mouse acrosome in the female reproductive tract, leading to successful fertilisation. Identification of the SPINK3 binding partner becomes essential to develop a contraceptive that works by prolonging the binding of SPINK3 to the sperm acrosome. We identified the SPINK3 receptor by using recombinant SPINK3 (rSPINK3). Testicular serine protease 1 (TESP1) was identified as the receptor for SPINK3 by 2D gel electrophoresis coupled with western blot analysis. To authenticate TESP1 as the receptor for SPINK3, sperm cells were incubated with TESP1 peptide antibody followed by determining the intracellular [Ca2+]i concentration by flow cytometry using Fluo-3 AM as a calcium probe. Furthermore, the 3D structures of SPINK3 and TESP1 were predicted by homology modelling (Schrodinger suite) using the crystal structure of pancreatic secretory trypsin inhibitor (PDB ID-1TGS) and human prostasin (PDB ID-3DFJ) as templates. The modelled protein structures were validated and subjected to molecular dynamics simulation (MDS) using GROMACS v5.0.5. Protein-protein docking was performed using HDOCK and the complex was validated by MDS. The results predicted that SPINK3 and TESP1 had strong binding affinity, with a dock score of -430.70 and 14 hydrogen bonds as key active site residues. If the binding affinity between SPINK3 and TESP1 could be increased, the SPINK3-TESP1 association will be prolonged, which will be helpful in the development of a male contraceptive.

In silico Screening of Natural Phytocompounds Towards Identification of Potential Lead Compounds to Treat COVID-19

COVID-19 is one of the members of the coronavirus family that can easily assail humans. As of now, 10 million people are infected and above two million people have died from COVID-19 globally. Over the past year, several researchers have made essential advances in discovering potential drugs. Up to now, no efficient drugs are available on the market. The present study aims to identify the potent phytocompounds from different medicinal plants (Zingiber officinale, Cuminum cyminum, Piper nigrum, Curcuma longa, and Allium sativum). In total, 227 phytocompounds were identified and screened against the proteins S-ACE2 and M ^pro through structure-based virtual screening approaches. Based on the binding affinity score, 30 active phytocompounds were selected. Amongst, the binding affinity for beta-sitosterol and beta-elemene against S-ACE2 showed -12.0 and -10.9 kcal/mol, respectively. Meanwhile, the binding affinity for beta-sitosterol and beta-chlorogenin against M ^pro was found to be -9.7 and -8.4 kcal/mol, respectively. Further, the selected compounds proceeded with molecular dynamics simulation, prime MM-GBSA analysis, and ADME/T property checks to understand the stability, interaction, conformational changes, binding free energy, and pharmaceutical relevant parameters. Moreover, the hotspot residues such as Lys31 and Lys353 for S-ACE2 and catalytic dyad His41 and Cys145 for M ^pro were actively involved in the inhibition of viral entry. From the in silico analyses, we anticipate that this work could be valuable to ongoing novel drug discovery with potential treatment for COVID-19.

In silico evaluation of isatin-based derivatives with RNA-dependent RNA polymerase of the novel coronavirus SARS-CoV-2

Isatin (1H-indole-2,3-dione)-containing compounds have been shown to possess several remarkable biological activities. We had previously explored a few isatin-based imidazole derivatives for their predicted dual activity against both inflammation and cancer. We explored 47 different isatin-based derivatives (IBDs) for other potential biological activities using in silico tools and found them to possess anti-viral activity. Using AutoDock tools, the binding site, binding energy, inhibitory constant/K_i and receptor-ligand interactions for each of the compounds were analyzed against SARS-CoV-2 RNA-dependent RNA polymerase (RdRp). The partition coefficient (logP) values were predicted using MedChem Designer tool. Based on the best K_i, binding energy and the ideal range of logP (between 1.0 and 3.0), 10 out of total 47 compounds were deemed to be prospective RdRp inhibitors. Some of these compounds gave better K_i, binding energy and logP values when compared to standard RdRp inhibitors, such as remdesivir (REM) (K_i = 15.61 μM, logP = 2.2; binding energy = -6.95), a clinically approved RdRp inhibitor and nine other RdRp inhibitors. The results showed that the 10 selected IBDs could be further explored. Molecular dynamics simulations (MDSs) showed that the selected RdRp-IBD complexes were highly stable compared to the native RdRp and RdRp-REM complex during 100 ns time periods. DFT studies were performed for the compounds 16a, 24a, 28a, 38a and 40a, to evaluate the charge transfer mechanism for the interactions between the IBDs and the RdRp residues. Among these, ADME profiling revealed that 28a is a possible lead compound which can be explored further for anti-RdRp activity in vitro. Communicated by Ramaswamy H. Sarma.

Structural insights on binding mechanism of CAD complexes (CPSase, ATCase and DHOase)

Pyrimidine biosynthetic pathway enzymes constitute an important target for the development of antitumor drugs. To understand the role of binding mechanisms underlying the inborn errors of pyrimidine biosynthetic pathway, structure and function of enzymes have been analyzed. Pyrimidine biosynthetic pathway is initiated by CAD enzymes that harbor the first three enzymatic activities facilitated by Carbamoyl Phosphate Synthetase (CPSase), Aspartate Transcarbamoylase (ATCase) and Dihydroorotase (DHOase). While being an attractive therapeutic target, the lack of data driven us to study the CPSase (CarA and CarB) and its mode of binding to ATCase and DHOase which are the major limitation for its structural optimization. Understanding the binding mode of CPSase, ATCase and DHOase could help to identify the potential interface hotspot residues that favor the mechanism behind it. The mechanistic insight into the CAD complexes were achieved through Molecular modeling, Protein-Protein docking, Alanine scanning and Molecular dynamics (MD) Studies. The hotspot residues present in the CarB region of carboxy phosphate and carbamoyl phosphate synthetic domains are responsible for the assembly of CAD (CPSase-ATCase-DHOase) complexes. Overall analysis suggests that the identified hotspot residues were confirmed by alanine scanning and important for the regulation of pyrimidine biosynthesis. MD simulations analysis provided the prolonged stability of the interacting complexes. The present study reveals the novel hotspot residues such as Glu134, Glu147, Glu154, Asp266, Lys269, Glu274, Asp333, Trp459, Asp526, Asp528, Glu533, Glu544, Glu546, Glu800, Val855, Asp877, Tyr884 and Gln919 which could be targeted for structure-based inhibitor design to potentiate the CAD mediated regulation of aggressive tumors.Communicated by Ramaswamy H. Sarma.

Identifying dual leucine zipper kinase (DLK) inhibitors using e-pharamacophore screening and molecular docking

Alzheimer's is a neural disorder causing gradual loss in structure and function of nerve cell. To treat such disorders, c-Jun N-terminal Kinase (JNK) Pathway inhibitors were developed by representing chemical compounds that were used to inhibit the JNK signaling pathways. DLK is the stress sensor and implicating as regulatory factor in JNK pathway. Therefore, in the present investigation, pharmacophore screening was tried to identify the chemical compounds that involving inhibition of DLK proteins. To explore the pharmacophore region and mode of binding with DLK protein, N- (I H-pyrazol-3-y l) pyridin-2-aminer inhibitors were docked with DLK. Results reveal the information on the interaction mechanism of protein and ligand with chemical characteristics required to inhibit DLK protein. Such predicted information (AAAARH) was used as query to find out potential novel lead compounds sourced from public database. As an outcome of 65 compounds were listed based on the fitness score (2≥), and were subjected to glide HTVS.SP and XP. Best performing 5 lead compounds were shortlisted for dynamic simulations. This exhibited a constant RMSD over 20 ns of timescale.

Identification of potential inhibitors for Penicillinbinding protein (PBP) from Staphylococcus aureus

Staphylococcus aureus is an infectious agent that causes severe skin and soft tissue infection in hospitalized patients. Therefore, it is of interest to develop potent inhibitors for S. aureus. Penicillin Binding protein (PBP) is a known drug target for inhibition of cell wall biosynthesis in S. aureus. Hence, PBP was screened with compounds from six databases using virtual screening approaches. Results shows that the screened lead compound produced higher docking score (-9.87 kcal/mol) compared to resistant drugs. Antimicrobial activity using screened lead compounds and resistant drugs showed maximum activity in potential screened compounds compared to resistant compounds.

Conformational insights into the inhibitory mechanism of phyto-compounds against Src kinase family members implicated in psoriasis

Psoriasis is a chronic immune mediated disorder of the skin. There is growing evidence that the Src family tyrosine kinases (SFK) are highly upregulated in psoriasis. The SFK are the key components of the signaling pathways triggering cell growth and differentiation in addition to the immune cascades. In the current work, the interactions between SFK and selective phyto-compounds were studied using molecular docking approach. Based on the results of docking and binding energy calculations quercetin was identified as potential lead compound. To get a deeper insight into the binding of quercetin with the SFK, a combined molecular dynamics and binding free energy calculations were performed. The binding of quercetin disrupted the intra-molecular contacts making the SFK compact except Src kinase. The MM/PBSA free energy decomposition analysis highlighted the significance of hydrophobic and polar residues which are involved in the binding of quercetin. An experimental validation was carried out against the activated forms of Fyn, Lyn and Src kinases, the top three proteins which showed high preference for quercetin. The flow cytometry analysis showed that the expression levels of Fyn, Lyn and Src kinases were dramatically increased in HaCaT cells. However, the treatment of quercetin at the concentration of 51.65 µM for 24 h markedly decreased their expression in HaCaT cells. Besides, similar results were also observed when the HaCaT cells were treated with the kinase inhibitor Ponitinib (1.43 µM) for 24 h.Communicated by Ramaswamy H. Sarma.

Identification of Pak1 inhibitors using water thermodynamic analysis

p21-activated kinases (Paks) play an integral component in various cellular diverse processes. The full activation of Pak is dependent upon several serine residues present in the N-terminal region, a threonine present at the activation loop, and finally the phosphorylation of these residues ensure the complete activation of Pak1. The present study deals with the identification of novel potent candidates of Pak1 using computational methods as anti-cancer compounds. A diverse energy based pharmacophore (e-pharmacophore) was developed using four co-crystal inhibitors of Pak1 having pharmacophore features of 5 (DRDRR), 6 (DRHADR), and 7 (RRARDRP and DRRDADH) hypotheses. These models were used for rigorous screening against e-molecule database. The obtained hits were filtered using ADME/T and molecular docking to identify the high affinity binders. These hits were subjected to hierarchical clustering using dendritic fingerprint inorder to identify structurally diverse molecules. The diverse hits were scored against generated water maps to obtain WM/MM ΔG binding energy. Furthermore, molecular dynamics simulation and density functional theory calculations were performed on the final hits to understand the stability of the complexes. Five structurally diverse novel Pak1 inhibitors (4835785, 32198676, 32407813, 76038049, and 32945545) were obtained from virtual screening, water thermodynamics and WM/MM ΔG binding energy. All hits revealed similar mode of binding pattern with the hinge region residues replacing the unstable water molecules in the binding site. The obtained novel hits could be used as a platform to design potent drugs that could be experimentally tested against cancer patients having increased Pak1 expression.

Insights from the Molecular modeling, docking analysis of illicit drugs and Bomb Compounds with Honey Bee Odorant Binding Proteins (OBPs)

Analysis of honeybee PBPs is of interest in the development of Biosensor applications. We described the predicted binding of 19 such compounds with 43-honey bee OBPs using molecular modeling, docking and phylogenetic analysis. Therefore, training the honeybees using preferred compounds formulate the bees to identify the illicit drugs and bomb compounds. Consequently, high docking score produced complex such OBP16-N-Phenyl-2-Napthalamine (-12.25k/mol), 3BJH-Crack Cocaine (-11.75k/mol), OBP10-Methadone (- 11.71k/mol), 1TUJ-Dronobinal Cannabis (-11.66k/mol), OBP13-Plasticizer (-11.27k/mol) and OBP24-Ecstasy (-10.89 k/mol) can be used to identify the compounds using biosensor application. The chemical reaction of the compounds for olfactory sensory was analyzed using DFT (Density Functional Theory) studies. Some of these compounds show high binding OBPs across distant phylogeny.

Green synthesis of silver nanoparticles using Lippia nodiflora aerial extract and evaluation of their antioxidant, antibacterial and cytotoxic effects

Silver nanoparticles biosynthesis has received increasing attention in the field of nanotechnology due totheir antimicrobial and biomedical applications. Green synthesis of metal nanoparticles is anticipated asa cost effective and eco-friendly alternative in the current research scenario. With this aim, the aqueousextracts made from the aerial parts of Lippia nodiflora were used as the reducing agents to synthesizesilver nanoparticles (AgNPs) and their antioxidant, antibacterial and cytotoxic properties have also beenevaluated. The biosynthesized AgNPs were characterized by UV–Visible Spectroscopy, Fourier-TransformInfrared Spectroscopy (FTIR) and X-Ray Diffraction (XRD) analysis. The AgNPs were found to be stable at−25.2 mV through zeta potential study. The morphology and size of synthesized silver nanoparticles wereconfirmed by Scanning Electron Microscope with energy dispersive spectra (SEM-EDX) and TransmissionElectron Microscopy (TEM) analysis with size range from 30 to 60 nm. Biosynthesized AgNPs exhibitedstrong antioxidant activity as well as showed potent antibacterial activity against human pathogenic bacteria.The cytotoxicity study of AgNPs was also revealed against MCF-7 breast cancer cell lines in a dosedependentmanner. The recognized bioactivity confirmed by the synthesized AgNPs directs towards theclinical use as an antioxidant, antibacterial and cytotoxic agent.

Identification of potential inhibitors for AIRS from de novo purine biosynthesis pathway through molecular modeling studies - a computational approach

In cancer, de novo pathway plays an important role in cell proliferation by supplying huge demand of purine nucleotides. Aminoimidazole ribonucleotide synthetase (AIRS) catalyzes the fifth step of de novo purine biosynthesis facilitating in the conversion of formylglycinamidine ribonucleotide to aminoimidazole ribonucleotide. Hence, inhibiting AIRS is crucial due to its involvement in the regulation of uncontrollable cancer cell proliferation. In this study, the three-dimensional structure of AIRS from P. horikoshii OT3 was constructed based on the crystal structure from E. coli and the modeled protein is verified for stability using molecular dynamics for a time frame of 100 ns. Virtual screening and induced fit docking were performed to identify the best antagonists based on their binding mode and affinity. Through mutational studies, the residues necessary for catalytic activity of AIRS were identified and among which the following residues Lys35, Asp103, Glu137, and Thr138 are important in determination of AIRS function. The mutational studies help to understand the structural and energetic characteristics of the specified residues. In addition to Molecular Dynamics, ADME properties, binding free-energy, and density functional theory calculations of the compounds were carried out to find the best lead molecule. Based on these analyses, the compound from the NCI database, NCI_121957 was adjudged as the best molecule and could be suggested as the suitable inhibitor of AIRS. In future studies, experimental validation of these ligands as AIRS inhibitors will be carried out.

P21-activated kinase 1 (Pak1) signaling influences therapeutic outcome in pancreatic cancer

BACKGROUND

Therapeutic resistance to gemcitabine in pancreatic ductal adenocarcinoma (PDAC) is attributed to various cellular mechanisms and signaling molecules that influence as a single factor or in combination.

DESIGN

In this study, utilizing in vitro p21-activated kinase 1 (Pak1) overexpression and knockdown cell line models along with in vivo athymic mouse tumor xenograft models and clinical samples, we demonstrate that Pak1 is a crucial signaling kinase in gemcitabine resistance.

RESULTS

Pak1 kindles resistance via modulation of epithelial-mesenchymal transition and activation of pancreatic stellate cells. Our results from gemcitabine-resistant and -sensitive cell line models showed that elevated Pak1 kinase activity is required to confer gemcitabine resistance. This was substantiated by elevated levels of phosphorylated Pak1 and ribonucleotide reductase M1 levels in the majority of human PDAC tumors when compared with normal. Delineation of the signaling pathway revealed that Pak1 confers resistance to gemcitabine by preventing DNA damage, inhibiting apoptosis and regulating survival signals via NF-κB. Furthermore, we found that Pak1 is an upstream interacting substrate of transforming growth factor β-activated kinase 1-a molecule implicated in gemcitabine resistance. Molecular mechanistic studies revealed that gemcitabine docks with the active site of Pak1; furthermore, gemcitabine treatment induces Pak1 kinase activity both in vivo and in cell-free system. Finally, results from athymic mouse tumor models illustrated that Pak1 inhibition by IPA-3 enhances the cytotoxicity of gemcitabine and brings about pancreatic tumor regression.

CONCLUSION

To our knowledge, this is the first study illustrating the mechanistic role of Pak1 in causing gemcitabine resistance via multiple signaling crosstalks, and hence Pak1-specific inhibitors will prove to be a better adjuvant with existing chemotherapy modality for PDAC.

Structure of Bacillus amyloliquefaciens alpha-amylase at high resolution: implications for thermal stability

The crystal structure of Bacillus amyloliquefaciens alpha-amylase (BAA) at 1.4 A resolution revealed ambiguities in the thermal adaptation of homologous proteins in this family. The final model of BAA is composed of two molecules in a back-to-back orientation, which is likely to be a consequence of crystal packing. Despite a high degree of identity, comparison of the structure of BAA with those of other liquefying-type alpha-amylases indicated moderate discrepancies at the secondary-structural level. Moreover, a domain-displacement survey using anisotropic B-factor and domain-motion analyses implied a significant contribution of domain B to the total flexibility of BAA, while visual inspection of the structure superimposed with that of B. licheniformis alpha-amylase (BLA) indicated higher flexibility of the latter in the central domain A. Therefore, it is suggested that domain B may play an important role in liquefying alpha-amylases, as its rigidity offers a substantial improvement in thermostability in BLA compared with BAA.