In-silico functional and structural annotation of hypothetical protein from Klebsiella pneumonia: A potential drug target

Targeting Matrix Metalloproteinase-3 for Dental Caries Prevention Using Herbal Isolates: MMP3 Inhibition by Cinnamic Acids

Objectives: Dental caries, a prevalent infectious disease affecting teeth, ranks highest among 328 diseases, according to a 2017 Lancet study. In demineralized human dentin, matrix metalloproteinase-3 (MMP3) functions as a proteoglycanase, contributing to the degradation of proteoglycan components. This process exposes collagen fibrils, thereby facilitating the demineralization of the dentin matrix. Inhibiting MMP3 shows potential for preventing dental caries. Methods: The binding affinity of 20 cinnamic acid derivatives, namely cynarin, chlorogenic acid, rosmarinic acid, cinnamyl caffeate, phenethyl caffeate, N-p-coumaroyltyramine, caffeic acid 3-glucoside, caffeic acid phenethyl ester, roscovitine, benzyl caffeate, o-coumaric acid, artepillin C, caffeic acid, methyl caffeate, 2-methylcinnamic acid, ferulic acid, drupanin, p-coumaric acid, cinnamic acid, and sinapinic acid, to the MMP3 catalytic cleft, was assessed utilizing AutoDock 4.0. Molecular dynamics simulation was then employed to analyze the stability of backbone atoms in free MMP3, MMP3-positive control inhibitor, and MMP3 complexed with the top-ranked cinnamic acid over a 100 ns computer simulation. Results: Four cinnamic acids demonstrated ΔG binding scores below -10 kcal/mol, with cynarin emerging as the most potent MMP3 inhibitor, featuring a ΔG binding score and inhibition constant value of -15.57 kcal/mol and 3.83 pM, respectively. The MMP3-cynarin complex exhibited stability after a 50 ns computer simulation, showing a root-mean-square deviation of 8 Å. Conclusions: The inhibition of MMP3 by cynarin, chlorogenic acid, rosmarinic acid, and cinnamyl caffeate holds promise as a potential preventive strategy for dental caries.

Investigating the potential of mono-chalcone compounds in targeting breast cancer receptors through network pharmacology, molecular docking, molecular dynamics simulation, antiproliferative effects, and gene expressions

The study aims to investigate various aspects of synthesized mono-chalcone compounds 5 and 8 concerning breast cancer, including network pharmacology, molecular docking, molecular dynamics (MD) simulations, antiproliferative effects, and gene expressions. Initially, the compounds underwent a network pharmacology analysis targeting breast cancer-related targets, with MalaCards, SwissTargetPrediction, and PharmMapper identifying 70 breast cancer target receptors. Subsequently, protein-protein interaction (PPI) network analysis revealed two distinct target gene clusters. Survival analysis identified seven significant target genes following Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment and Gene Ontology (GO) evaluation. Molecular docking and MD simulations were conducted on these seven target genes (AKT2, BRAF, ESR1, FGFR1, IGF1, IGF1R, and KIT), revealing that compound 8 exhibited the highest binding affinities, as well as better stability and compactness when interacting with the targeted proteins. Next, the compounds underwent cell viability assay and gene expression analysis to validate the in silico findings. Both compounds demonstrated the ability to suppress breast cancer proliferation, with compound 8 showing increased selectivity in targeting breast cancer cells while causing minimal harm to normal breast cells. The suppression of breast cancer cell proliferation was attributed to decreased expression levels of AKT2, BRAF, FGFR1, IGF1, IGF1R, KIT, and ESR1. Hence, the results provide insights into the molecular interaction responsible for the anti-breast cancer capabilities of mono-chalcone compounds.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13205-024-03991-y.

In silico exploration of phenolics as modulators of penicillin binding protein (PBP) 2× of Streptococcus pneumoniae

Infections caused by multidrug-resistant Streptococcus pneumoniae remain the leading cause of pneumonia-related deaths in children < 5 years globally, and mutations in penicillin-binding protein (PBP) 2 × have been identified as the major cause of resistance in the organism to beta-lactams. Thus, the development of new modulators with enhanced binding of PBP2x is highly encouraged. In this study, phenolics, due to their reported antibacterial activities, were screened against the active site of PBP2x using structure-based pharmacophore and molecular docking techniques, and the ability of the top-hit phenolics to inhibit the active and allosteric sites of PBP2x was refined through 120 ns molecular dynamic simulation. Except for gallocatechin gallate and lysidicichin, respectively, at the active and allosteric sites of PBP2x, the top-hit phenolics had higher negative binding free energy (ΔGbind) than amoxicillin [active site (- 19.23 kcal/mol), allosteric site (- 33.75 kcal/mol)]. Although silicristin had the best broad-spectrum effects at the active (- 38.41 kcal/mol) and allosteric (- 50.54 kcal/mol) sites of PBP2x, the high thermodynamic entropy (4.90 Å) of the resulting complex might suggest the need for its possible structural refinement for enhanced potency. Interestingly, silicristin had a predicted synthetic feasibility score of < 5 and quantum calculations using the DFT B3LYP/6-31G+ (dp) revealed that silicristin is less stable and more reactive than amoxicillin. These findings point to the possible benefits of the top-hit phenolics, and most especially silicristin, in the direct and synergistic treatment of infections caused by S. pneumoniae. Accordingly, silicristin is currently the subject of further confirmatory in vitro research.

Discovery of potential quality markers of Fritillariae thunbergii bulbus in pneumonia by combining UPLC-QTOF-MS, network pharmacology, and molecular docking

Fritillariae thunbergii bulbus (FTB) is a popular Chinese herbal medicine with various applications in respiratory diseases. The quality evaluation of FTB has been insufficient to date, as the active ingredients and mechanisms of action of FTB remain unclear. This study proposes a novel strategy for exploring the quality markers (Q-markers) of FTB based on UPLC-QTOF-MS analysis, network pharmacology, molecular docking, and molecular dynamics (MD) simulation. A total of 26 compounds in FTB were identified by UPLC-QTOF-MS. Ten of these compounds were screened as Q-markers based on network pharmacology for their anti-pneumonia effects, including imperialine, peimisine, peiminine, ebeiedinone, zhebeirine, puqiedine, 9-hydroxy-10,12-octadecadienoic acid, (9Z,12Z,15Z)-13-hydroxy-9,12,15-octadecatrienoic acid, 9,12,15-octadecatrienoic acid, and (2E,4Z,7Z,10Z,13Z,16Z,19Z)-2,4,7,10,13,16,19-docosaheptaenoic acid methyl ester (DAME). These Q-markers were predicted to act on multiple targets and pathways associated with pneumonia. Molecular docking results revealed that most of the Q-markers showed high affinity with at least one of the main targets of pneumonia, and the top ten complexes were confirmed with MD simulation. Network pharmacology indicated that FTB may act on the TNF signaling pathway, HIF-1 signaling pathway, JAK-STAT signaling pathway, etc. The results demonstrated that imperialine (P8), peimisine (P9), peiminine (P11), ebeiedine (P15), zhebeirine (P16), and puqiedine (P18) may be potential Q-markers of FTB, and AKT1, IL-6, VEGFA, TP53, EGFR, STAT3, PPARG, MMP9, and CASP3 may be promising therapeutic targets for pneumonia treatment that are worthy of further research.

In silico exploration of hypothetical proteins in Neisseria gonorrhoeae for identification of therapeutic targets

Neisseria gonorrhoeae, a World Health Organization (WHO) declared superbug and the second-most frequent cause of bacterial sexually transmitted infections worldwide is responsible for gonorrhea. Hypothetical proteins are gene products that are predicted to be encoded by a particular gene based on the DNA sequence, but their specific functions and characteristics have not been experimentally determined or verified. In the context of this research, annotating hypothetical proteins is crucial for identifying their potential as therapeutic targets. Without proper annotation, these proteins would remain vague, hindering efforts to understand their roles in disease. The methodology used aims to bridge this gap by employing algorithm-based tools and software to annotate hypothetical proteins and assess their suitability as therapeutic targets based on factors such as essentiality, virulence, subcellular localization, and druggability. Out of 716 N. gonorrhoeae hypothetical proteins reported in UniProt, assessment of crucial pathogenic factors, including essentiality, virulence, subcellular localization, and druggability, effectively filtered and prioritized the hypothetical proteins for further therapeutic exploration and lead to 5 proteins being chosen as targets. The molecular docking studies conducted identified 10 hits targeting the five targets. Conclusively, this study aided in identification of targets and hit compounds for therapeutic targeting of gonorrhea disease.

Graphical abstract

Supplementary Information

The online version contains supplementary material available at 10.1007/s40203-023-00186-w.

Identification and mechanistic exploration of structural and conformational dynamics of NF-kB inhibitors: rationale insights from in silico and in vitro studies

Increased expression of target genes that code for proinflammatory chemical mediators results from a series of intracellular cascades triggered by activation of dysregulated NF-κB signaling pathway. Dysfunctional NF-kB signaling amplifies and perpetuates autoimmune responses in inflammatory diseases, including psoriasis. This study aimed to identify therapeutically relevant NF-kB inhibitors and elucidate the mechanistic aspects behind NF-kB inhibition. After virtual screening and molecular docking, five hit NF-kB inhibitors opted, and their therapeutic efficacy was examined using cell-based assays in TNF-α stimulated human keratinocyte cells. To investigate the conformational changes of target protein and inhibitor-protein interaction mechanisms, molecular dynamics (MD) simulations, binding free energy calculations together with principal component (PC) analysis, dynamics cross-correlation matrix analysis (DCCM), free energy landscape (FEL) analysis and quantum mechanical calculations were carried out. Among identified NF-kB inhibitors, myricetin and hesperidin significantly scavenged intracellular ROS and inhibited NF-kB activation. Analysis of the MD simulation trajectories of ligand-protein complexes revealed that myricetin and hesperidin formed energetically stabilized complexes with the target protein and were able to lock NF-kB in a closed conformation. Myricetin and hesperidin binding to the target protein significantly impacted conformational changes and internal dynamics of amino acid residues in protein domains. Tyr57, Glu60, Lys144 and Asp239 residues majorly contributed to locking the NF-kB in a closed conformation. The combinatorial approach employing in silico tools integrated with cell-based approaches substantiated the binding mechanism and NF-kB active site inhibition by the lead molecule myricetin, which can be explored as a viable antipsoriatic drug candidate associated with dysregulated NF-kB.Communicated by Ramaswamy H. Sarma.

Mechanisms Involved in Therapeutic Effects of Scutellaria baicalensis Georgi in Oral Squamous Cell Carcinoma Based on Systems Biology and Structural Bioinformatics Approaches

Objective

Oral squamous cell carcinoma (OSCC) is the most frequent oral cancer, constituting more than 90% of all oral carcinomas. The 5-year survival rate of OSCC patients is not satisfactory, and therefore, there is an urgent need for new practical therapeutic approaches besides the current therapies to overcome OSCC. Scutellaria baicalensis Georgi (SBG) is a plant of the family Lamiaceae with several pharmaceutical properties such as antioxidant, anti-inflammatory, and anticancer effects. Previous studies have demonstrated the curative effects of SBG in OSCC.

Methods

A systems biology approach was conducted to identify differentially expressed miRNAs (DEMs) in OSCC patients with a dismal prognosis compared to OSCC patients with a favorable prognosis. A protein interaction map (PIM) was built based on DEMs targets, and the hub genes within the PIM were indicated. Subsequently, the prognostic role of the hubs was studied using Kaplan-Meier curves. Next, the binding affinity of SBG's main components, including baicalein, wogonin, oroxylin-A, salvigenin, and norwogonin, to the prognostic markers in OSCC was evaluated using molecular docking analysis.

Results

Survival analysis showed that overexpression of CAV1, SERPINE1, ACTB, SMAD3, HMGA2, MYC, EIF2S1, HSPA4, HSPA5, and IL6 was significantly related to a poor prognosis in OSCC. Besides, molecular docking analysis demonstrated the ΔGbinding and inhibition constant values between SBG's main components and SERPINE1, ACTB, HMGA2, EIF2S1, HSPA4, and HSPA5 were as <-8.00 kcal/mol and nanomolar concentration, respectively. The most salient binding affinity was observed between wogonin and SERPINE1 with a criterion of ΔGbinding < -10.02 kcal/mol.

Conclusion

The present results unraveled potential mechanisms involved in therapeutic effects of SBG in OSCC based on systems biology and structural bioinformatics analyses.

Ginkgolide B attenuates cerebral ischemia-reperfusion injury via inhibition of ferroptosis through disrupting NCOA4-FTH1 interaction

ETHNOPHARMACOLOGICAL RELEVANCE

Cerebral ischemia/reperfusion (I/R) injury is a major cause of neuronal damage and death. Ginkgolide B (GB) has been shown to exhibit neuroprotective effects in various brain injury models.

AIM OF STUDY

The aim of study was to investigate the potential role of GB in protecting against cerebral I/R injury and explore the underlying mechanisms.

MATERIALS AND METHODS

Adult male Sprague-Dawley rats were exposed to transient middle cerebral artery occlusion (tMCAO) followed by reperfusion in order to trigger cerebral I/R injury. The rats were treated with different doses of GB, vehicle control or positive drug. Neurological function, infarct volume, and levels of ferroptosis markers were evaluated. In vitro experiments were performed using OGD/R-induced PC12 cells to further investigate the effects of GB on ferroptosis and its mechanisms. In addition, molecular docking, and microscale thermophoresis (MST) assay were conducted to explore the combination of GB and NCOA4.

RESULTS

Reduced infarct volume and enhanced neurological function were signs of dose-dependent protection from cerebral I/R injury by GB therapy. Additionally, GB treatment had an impact on the levels of oxidative stress and ferroptosis markers, including reactive oxygen species (ROS), malondialdehyde (MDA), superoxide dismutase (SOD), and Fe2+ in the cerebral environment during IR injury. Moreover, relevant ferroptosis key factors such as ACSL4, GPX4, FTH1, and NCOA4 can be regulated by GB. In OGD/R-induced PC12 cells, GB protected against ferroptosis by inhibiting autophagy and disrupting the interaction of NCOA4-FTH1.

CONCLUSION

Our findings suggest that GB may protect against cerebral I/R injury by inhibiting ferroptosis through disrupting NCOA4-FTH1 interaction. GB has potential therapeutic applications for cerebral I/R injury, and further investigation of the underlying mechanisms and clinical trials are warranted.

Ultrasound-assisted extraction of emodin from Rheum officinale Baill and its antibacterial mechanism against Streptococcus suis based on CcpA

Emodin was extracted from Rheum officinale Baill by ultrasound-assisted extraction (UAE), and ethanol was chosen as the suitable solvent through SEM and molecular dynamic simulation. Under the optimum conditions (power 541 W, time 23 min, liquid to material ratio 13:1 mL/g, ethanol concentration 83 %) predicted by RSM, the yield of emodin was 2.18 ± 0.11 mg/g. Moreover, ultrasound power and time displayed the significant effects on the extraction process. Extracting dynamics analysis indicated that the extraction process of emodin by UAE conformed to Fick's second diffusion law. The results of antibacterial experiments suggested that emodin can damage cell membrane and inhibit the expression of cps2A, sao, mrp, epf, neu and the hemolytic activity of S. suis. Biolayer interferometry and FT-IR multi-peak fitting assays demonstrated that emodin induced a secondary conformational shift in CcpA. Molecular docking and molecular dynamics confirmed that emodin bound to CcpA through hydrogen bonding (ALA248, GLU249, GLY129 and ASN196) and π-π T-shaped interaction (TYR225 and TYR130), and the mutation of amino acid residues affected the affinity of CcpA to emodin. Therefore, emodin inhibited the sugar utilization of S. suis through binding to CcpA, and CcpA may be a potential target to inhibit the growth of S. suis.

Molecular Modelling of Resveratrol Derivatives with SIRT1 for the Stimulation of Deacetylase Activity

BACKGROUND

Resveratrol is a polyphenol that is found in plants and has been proposed to have a potential therapeutic effect through the activation of SIRT1, which is a crucial member of the mammalian NAD+ -dependent deacetylases. However, how its activity is enhanced toward specific substrates by resveratrol derivatives has not been studied. This study aimed to evaluate the types of interaction of resveratrol and its derivatives with SIRT1 as the target protein, as well as to find out the best ligand with the strangest interaction with SIRT1.

MATERIALS AND METHODS

In this study, we employed the extensive molecular docking analysis using AutoDock Vina to comparatively evaluate the interactions of resveratrol derivatives (22 molecules from the ZINC database) as ligands with SIRT1 (PDB ID: 5BTR) as a receptor. The ChemDraw and Chem3D tools were used to prepare 3D structures of all ligands and energetically minimize them by the MM2 force field.

RESULTS

The molecular docking and visualizations showed that conformational change in resveratrol derivatives significantly influenced the parameter for docking results. Several types of interactions, including conventional hydrogen bonds, carbon-hydrogen bonds, Pi-donor hydrogen bonds, and Pi-Alkyl, were found via docking analysis of resveratrol derivatives and SIRT1 receptors. The possible activation effect of resveratrol 4'-(6-galloylglucoside) with ZINC ID: ZINC230079516 with higher binding energy score (-46.8608 kJ/mol) to the catalytic domain (CD) of SIRT1 was achieved at the maximum value for SIRT1, as compared to resveratrol and its other derivatives.

CONCLUSION

Finally, resveratrol 4'-(6-galloylglucoside), as a derivative for resveratrol, has stably interacted with the CD of SIRT1 and might be a potential effective activator for SIRT1.

Lipase and Protease Production Ability of Multi-drug Resistant Bacteria Worsens the Outcomes of Wound Infections

BACKGROUND

Surgical site infections are one of the major clinical problems in surgical departments that cost hundreds of millions of dollars to healthcare systems around the world.

AIM

The study aimed to address the pressing issue of surgical site infections, which pose significant clinical and financial burdens on healthcare systems globally. Recognizing the substantial costs incurred due to these infections, the research has focused on understanding the role of lipase and protease production by multi-drug resistant bacteria isolated from surgical wounds in the development of post-surgical wound infections.

METHODS

For these purposes, 153 pus specimens were collected from patients with severe post-surgical wound infections having prolonged hospital stays. The specimens were inoculated on appropriate culture media. Gram staining and biochemical tests were used for the identification of bacterial growth on suitable culture media after 24 hours of incubation. The isolated pathogens were then applied for lipase and protease, key enzymes that could contribute to wound development, on tributyrin and skimmed milk agar, respectively. Following the CSLI guidelines, the Kirby-Bauer disc diffusion method was used to assess antibiotic susceptibility patterns. The results revealed that a significant proportion of the samples (127 out of 153) showed bacterial growth of Gram-negative (n = 66) and Gram-positive (n = 61) bacteria. In total, isolated 37 subjects were declared MDR due to their resistance to three or more than three antimicrobial agents. The most prevalent bacteria were Staphylococcus aureus (29.13%), followed by S. epidermidis (18.89%), Klebsiella pneumoniae (18.89%), Escherichia coli (14.96%), Pseudomonas aeruginosa (10.23%), and Proteus mirabilis (7.87%). Moreover, a considerable number of these bacteria exhibited lipase and protease activity with 70 bacterial strains as lipase positive on tributyrin agar, whereas 74 bacteria showed protease activity on skimmed milk agar with P. aeruginosa as the highest lipase (69.23%) and protease (76.92%) producer, followed by S. aureus (lipase 62.16% and protease 70.27%).

RESULTS

The antimicrobial resistance was evaluated among enzyme producers and non-producers and it was found that the lipase and protease-producing bacteria revealed higher resistance to selected antibiotics than non-producers. Notably, fosfomycin and carbapenem were identified as effective antibiotics against the isolated bacterial strains. However, gram-positive bacteria displayed high resistance to lincomycin and clindamycin, while gram-negative bacteria were more resistant to cefuroxime and gentamicin.

CONCLUSION

In conclusion, the findings suggest that lipases and proteases produced by bacteria could contribute to drug resistance and act as virulence factors in the development of surgical site infections. Understanding the role of these enzymes may inform strategies for preventing and managing post-surgical wound infections more effectively.

Ab initio and comparative 3D modeling of FAM222A-encoded protein and target-driven-based virtual screening for the identification of novel therapeutics against Alzheimer's disease

The complex nature of Alzheimer's disease (AD) makes it difficult to understand the exact molecular processes leading to neuron death. However, two molecular factors - the production of amyloid-beta plaques and tau tangles - are considered to be linked to AD. A genetic marker for brain atrophy, FAM222A, has been identified by the unique cross-phenotype meta-analysis of genetics imaging and the molecular features show an interaction between the protein aggregatin encoded by FAM222A and amyloid beta (Aβ)-peptide (1-42) via its N-terminal Aβ binding domain, thus increasing Aβ aggregation. Function of Aggregatin protein is unclear, and its 3D structure has not been investigated in experimental analysis, so far. Hence, in the present study, first time in literature, 3D models of FAM222A-encoded Aggregatin were systematically constructed by applying diverse homology modeling approaches and they were used as target structures at the virtual screening of FDA-approved drugs and drugs currently under research in clinical trials. Then, the identified hit molecules were chosen for further molecular dynamics (MD) simulations and post-MD analyses. Our integrated ligand-based and protein-driven-based virtual screening results show that Cefpiramide, Diniprofylline, Fostriecin, and Droperidol may target Aggregatin.

Potential inhibitors of FemC to combat Staphylococcus aureus: virtual screening, molecular docking, dynamics simulation, and MM-PBSA analysis

FemC is a methicillin resistance factor involved in the alterations of peptidoglycan and glutamine synthesis in Staphylococcus aureus. To identify the potent antibacterial agents, antibacterial molecules were screened against the predicted and validated FemC model. Based on docking scores, presence of essential interactions with active site residues of FemC, pharmacokinetic, and ADMET properties, six candidates were shortlisted and subjected to molecular dynamics to evaluate the stability of FemC-ligand complexes. Further, per residue decomposition analysis and Molecular Mechanics/Poisson-Boltzmann Surface Area (MMPBSA) analysis confirmed that S15, M16, S17, R31, R43, Q47, K48 and R49 of FemC played a vital role in the formation of lower energy stable FemC-inhibitor(s) complexes. Therefore, in the present study, the reported six molecules (Z317461228, Z92241701, Z30923155, Z30202349, Z2609517102 and Z92470167) may pave the path to design the scaffold of novel potent antimicrobials against S. aureus.Communicated by Ramaswamy H. Sarma.

In silico identification and characterization of potential druggable targets among hypothetical proteins of Leptospira interrogans serovar Copenhageni: a comprehensive bioinformatics approach

Leptospirosis is one of the neglected zoonosis, affecting human and animal populations worldwide. Reliable effective therapeutics and concerns to look for more research into the molecular analysis of its genome is therefore needed. In the genomic pool of the Leptospira interrogans many hypothetical proteins are still uncharacterized. In the current research, we performed extensive in silico analysis to prioritize the potential hypothetical proteins of L. interrogans serovar Copenhageni via stepwise reducing the available hypothetical proteins (Total 3606) of the assembly to only 15, based on non-homologous to homosapien, essential, functional, virulent, cellular localization. Out of them, only two proteins WP_000898918.1 (Hypothetical Protein 1) & WP_001014594.1 (Hypothetical Protein 2) were found druggable and involved in protein-protein interaction network. The 3 D structures of these two target proteins were predicted via ab initio homology modeling followed by structures refinement and validation, as no structures were available till date. The analysis also revealed that the functional domains, families and protein-protein interacting partners identified in both proteins are crucial for the survival of the bacteria. The binding cavities were predicted for both the proteins through blind and specific protein-ligand docking with their respective ligands and inhibitors and were found to be in accordance with the druggable sites predicted by DoGSiteScorer. The docking interactions were found within the active functional domains for both the proteins while for Hypothetical Protein 2, the same residues were involved in interactions with Cytidine-5'-triphosphate in blind and specific docking. Furthermore, the simulations of molecular dynamics and free binding energy revealed the stable substrate binding and efficient binding energies, and were in accordance to our docking results. The work predicted two unique hypothetical proteins of L. interrogans as a potential druggable targets for designing of inhibitors for them.Communicated by Ramaswamy H. Sarma.

Promising antibacterials for LLM of Staphylococcus aureus using virtual screening, molecular docking, dynamics, and MMPBSA

In S. aureus, lipophilic membrane (LLM) protein is a methicillin resistance factor and is an essential role in peptidoglycan metabolism. The virtual screening of antibacterial molecules against the model of LLM was performed to identify the potent antibacterial molecules. Molecular docking results of pharmacokinetic filtered molecules illustrated that five molecules had higher binding affinities than tunicamycin (TUM) and were stabled via non-covalent interactions (hydrogen bond and hydrophobic interactions) at the active site of LLM. Further, molecular dynamics results revealed that binding of identified antibacterial molecules with LLM resulted in stable LLM-inhibitor(s) complexes. Molecular Mechanics/Position-Boltzmann Surface Area (MMPBSA) analysis showed that LLM-inhibitor(s) complexes had high binding affinities in the range of -213.49 ± 2.24 to -227.42 ± 3.05 kJ/mol. The amino acid residues decomposition analysis confirmed that identified antibacterial molecules bound at the active site (Asn148, Leu149, Asp151, Asp208, His269, His271, and His272) of LLM. Noticeably, the current study found five antibacterial molecules (BDE 27575101, BDE 33638168, BDE 33672484, LAS 51502073, and BDE 25098678) were highly potent than TUM and even than earlier reported molecules. Therefore, here reported antibacterial molecules may be used directly or developed to inhibit LLM of S. aureus.Communicated by Ramaswamy H. Sarma.

Exploring and targeting potential druggable antimicrobial resistance targets ArgS, SecY, and MurA in Staphylococcus sciuri with TCM inhibitors through a subtractive genomics strategy

Staphylococcus sciuri (also currently Mammaliicoccus sciuri) are anaerobic facultative and non-motile bacteria that cause significant human pathogenesis such as endocarditis, wound infections, peritonitis, UTI, and septic shock. Methicillin-resistant S. sciuri (MRSS) strains also infects animals that include healthy broilers, cattle, dogs, and pigs. The emergence of MRSS strains thereby poses a serious health threat and thrives the scientific community towards novel treatment options. Herein, we investigated the druggable genome of S. sciuri by employing subtractive genomics that resulted in seven genes/proteins where only three of them were predicted as final targets. Further mining the literature showed that the ArgS (WP_058610923), SecY (WP_058611897), and MurA (WP_058612677) are involved in the multi-drug resistance phenomenon. After constructing and verifying the 3D protein homology models, a screening process was carried out using a library of Traditional Chinese Medicine compounds (consisting of 36,043 compounds). The molecular docking and simulation studies revealed the physicochemical stability parameters of the docked TCM inhibitors in the druggable cavities of each protein target by identifying their druggability potential and maximum hydrogen bonding interactions. The simulated receptor-ligand complexes showed the conformational changes and stability index of the secondary structure elements. The root mean square deviation (RMSD) graph showed fluctuations due to structural changes in the helix-coil-helix and beta-turn-beta changes at specific points where the pattern of the RMSD and root mean square fluctuation (RMSF) (< 1.0 Å) support any major domain shifts within the structural framework of the protein-ligand complex and placement of ligand was well complemented within the binding site. The β-factor values demonstrated instability at few points while the radius of gyration for structural compactness as a time function for the 100-ns simulation of protein-ligand complexes showed favorable average values and denoted the stability of all complexes. It is assumed that such findings might facilitate researchers to robustly discover and develop effective therapeutics against S. sciuri alongside other enteric infections.

Human adenovirus DNA polymerase is evolutionarily and functionally associated with human telomerase reverse transcriptase based on in silico molecular characterization that implicate abacavir and zidovudine

Human adenoviruses (HAdVs) are non-enveloped, small double stranded DNA (dsDNA) viruses that cause asymptomatic infections, clinical syndromes and significant susceptibility to infections in immunocompromised people. The aim of the present study was to identify critical host proteins and HAdV hypothetical proteins that could be developed as potential host-viral targets for antiHAdV therapy. Here, the function of selected hypothetical proteins of HAdV based on phylogenetic relationship with the therapeutic targets of antiretroviral drugs of human immunodeficiency virus (HIV) was predicted computationally, and characterized the molecular dynamics and binding affinity of DNA polymerase of HAdV. Thirty-eight hypothetical proteins (HPs) of human adenovirus (HAdV) were used in this study. The results showed that HAdV DNA polymerase (P03261) is related to Human TERT (O14746) and HLA-B (P01889) genes. The protein-protein interaction of human five molecular targets (PNP, TERT, CCR5, HLA-B, and NR1I2) of ARVDs are well-coordinated/networked with CD4, AHR, FKBP4, NR3C1, HSP90AA1, and STUB1 proteins in the anti-HIV infection mechanism. The results showed that the free energy score of abacavir and zidovudine binding to HAdV DNA polymerase are -5.8 and -5.4 kcal mol^-1 respectively. Also, the control drug, cidofovir and ganciclovir have less binding affinity for DNA polymerase of HAdV when compare to that of abacavir and zidovudine. Similarity was observed in the binding of abacavir and zidovudine to HAdV DNA polymerase (ASP742, ALA743, LEU772, ARG773 and VAL776). In conclusion, combination of abacavir and zidovudine was predicted to be potential therapy for controlling HAdV infection targeting HAdV DNA polymerase.

Molecular modeling study of natural products as potential bioactive compounds against SARS-CoV-2

CONTEXT

The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), the causative agent of the COVID-19 infection and responsible for millions of victims worldwide, remains a significant threat to public health. Even after the development of vaccines, research interest in the emergence of new variants is still prominent. Currently, the focus is on the search for effective and safe drugs, given the limitations and side effects observed for the synthetic drugs administered so far. In this sense, bioactive natural products that are widely used in the pharmaceutical industry due to their effectiveness and low toxicity have emerged as potential options in the search for safe drugs against COVID-19. Following this line, we screened 10 bioactive compounds derived from cholesterol for molecules capable of interacting with the receptor-binding domain (RBD) of the spike protein from SARS-CoV-2 (SC2Spike), responsible for the virus's invasion of human cells. Rounds of docking followed by molecular dynamics simulations and binding energy calculations enabled the selection of three compounds worth being experimentally evaluated against SARS-CoV-2.

METHODS

The 3D structures of the cholesterol derivatives were prepared and optimized using the Spartan 08 software with the semi-empirical method PM3. They were then exported to the Molegro Virtual Docking (MVD®) software, where they were docked onto the RBD of a 3D structure of the SC2Spike protein that was imported from the Protein Data Bank (PDB). The best poses obtained from MVD® were subjected to rounds of molecular dynamics simulations using the GROMACS software, with the OPLS/AA force field. Frames from the MD simulation trajectories were used to calculate the ligand's free binding energies using the molecular mechanics - Poisson-Boltzmann surface area (MM-PBSA) method. All results were analyzed using the xmgrace and Visual Molecular Dynamics (VMD) software.

An integrated computational approach towards novel drugs discovery against polyketide synthase 13 thioesterase domain of Mycobacterium tuberculosis

The acquired drug resistance by Mycobacterium tuberculosis (M. tuberculosis) to antibiotics urges the need for developing novel anti-M. tuberculosis drugs that possess novel mechanism of action. Since traditional drug discovery is a labor-intensive and costly process, computer aided drug design is highly appreciated tool as it speeds up and lower the cost of drug development process. Herein, Asinex antibacterial compounds were virtually screened against thioesterase domain of Polyketide synthase 13, a unique enzyme that forms α-alkyl β-ketoesters as a direct precursor of mycolic acids which are essential components of the lipid-rich cell wall of M. tuberculosis. The study identified three drug-like compounds as the most promising leads; BBB_26582140, BBD_30878599 and BBC_29956160 with binding energy value of - 11.25 kcal/mol, - 9.87 kcal/mol and - 9.33 kcal/mol, respectively. The control molecule binding energy score is -9.25 kcal/mol. Also, the docked complexes were dynamically stable with maximum root mean square deviation (RMSD) value of 3 Å. Similarly, the MM-GB\PBSA method revealed highly stable complexes with mean energy values < - 75 kcal/mol for all three systems. The net binding energy scores are validated by WaterSwap and entropy energy analysis. Furthermore, The in silico druglike and pharmacokinetic investigation revealed that the compounds could be suitable candidates for additional experimentations. In summary, the study findings are significant, and the compounds may be used in experimental validation pipeline to develop potential drugs against drug-resistant tuberculosis.

Molecular Docking of Bacterial Protein Modulators and Pharmacotherapeutics of Carica papaya Leaves as a Promising Therapy for Sepsis: Synchronising In Silico and In Vitro Studies

Sepsis is a serious health concern globally, which necessitates understanding the root cause of infection for the prevention of proliferation inside the host's body. Phytochemicals present in plants exhibit antibacterial and anti-proliferative properties stipulated for sepsis treatment. The aim of the study was to determine the potential role of Carica papaya leaf extract for sepsis treatment in silico and in vitro. We selected two phytochemical compounds, carpaine and quercetin, and docked them with bacterial proteins, heat shock protein (PDB ID: 4PO2), surfactant protein D (PDB ID: 1PW9), and lactobacillus bacterial protein (PDB ID: 4MKS) against imipenem and cyclophosphamide. Quercetin showed the strongest interaction with 1PW9 and 4MKS proteins. The leaves were extracted using ethanol, methanol, and water through Soxhlet extraction. Total flavonoid content, DPPH assay, HPTLC, and FTIR were performed. In vitro cytotoxicity of ethanol extract was screened via MTT assay on the J774 cell line. Ethanol extract (EE) possessed the maximum number of phytocomponents, the highest amount of flavonoid content, and the maximum antioxidant activity compared to other extracts. FTIR analysis confirmed the presence of N-H, O-H, C-H, C=O, C=C, and C-Cl functional groups in ethanol extract. Cell viability was highest (100%) at 25 µg/mL of EE. The present study demonstrated that the papaya leaves possessed antibacterial and cytotoxic activity against sepsis infection.

Molecular modeling, simulation and docking of Rv1250 protein from Mycobacterium tuberculosis

Computational prediction and protein structure modeling have come to the aid of various biological problems in determining the structure of proteins. These technologies have revolutionized the biological world of research, allowing scientists and researchers to gain insights into their biological questions and design experimental research much more efficiently. Pathogenic Mycobacterium spp. is known to stay alive within the macrophages of its host. Mycobacterium tuberculosis is an acid-fast bacterium that is the most common cause of tuberculosis and is considered to be the main cause of resistance of tuberculosis as a leading health issue. The genome of Mycobacterium tuberculosis contains more than 4,000 genes, of which the majority are of unknown function. An attempt has been made to computationally model and dock one of its proteins, Rv1250 (MTV006.22), which is considered as an apparent drug-transporter, integral membrane protein, and member of major facilitator superfamily (MFS). The most widely used techniques, i.e., homology modeling, molecular docking, and molecular dynamics (MD) simulation in the field of structural bioinformatics, have been used in the present work to study the behavior of Rv1250 protein from M. tuberculosis. The structure of unknown TB protein, i.e., Rv1250 was retrived using homology modeling with the help of I-TASSER server. Further, one of the sites responsible for infection was identified and docking was done by using the specific Isoniazid ligand which is an inhibitor of this protein. Finally, the stability of protein model and analysis of stable and static interaction between protein and ligand molecular dynamic simulation was performed at 100 ns The designing of novel Rv1250 enzyme inhibitors is likely achievable with the use of proposed predicted model, which could be helpful in preventing the pathogenesis caused by M. tuberculosis. Finally, the MD simulation was done to evaluate the stability of the ligand for the specific protein.

Drug-Repurposing Approach To Combat Staphylococcus aureus: Biomolecular and Binding Interaction Study

Staphylococcus aureus is considered as one of the most widespread bacterial pathogens and continues to be a prevalent cause of mortality and morbidity across the globe. FmtA is a key factor linked with methicillin resistance in S. aureus. Consequently, new antibacterial compounds are crucial to combat S. aureus resistance. Here, we present the virtual screening of a set of compounds against the available crystal structure of FmtA. The findings indicate that gemifloxacin, paromomycin, streptomycin, and tobramycin were the top-ranked potential drug molecules based on the binding affinity. Furthermore, these drug molecules were analyzed with molecular dynamics simulations, which showed that the identified molecules formed highly stable FmtA-inhibitor(s) complexes. Molecular mechanics Poisson-Boltzmann surface area and quantum mechanics/molecular mechanics calculations suggested that the active site residues (Ser127, Lys130, Tyr211, and Asp213) of FmtA are crucial for the interaction with the inhibitor(s) to form stable protein-inhibitor(s) complexes. Moreover, fluorescence- and isothermal calorimetry-based binding studies showed that all the molecules possess dissociation constant values in the micromolar scale, revealing a strong binding affinity with FmtA^Δ80, leading to stable protein-drug(s) complexes. The findings of this study present potential beginning points for the rational development of advanced, safe, and efficacious antibacterial agents targeting FmtA.