Structure-Based Identification of Potential Drugs Against FmtA of Staphylococcus aureus: Virtual Screening, Molecular Dynamics, MM-GBSA, and QM/MM

Identification of novel CA IX inhibitor: Pharmacophore modeling, docking, DFT, and dynamic simulation

Human Carbonic anhydrase IX (hCA IX) is found to be an essential biomarker for the treatment of hypoxic tumors in both the early and metastatic stages of cancer. Due to its active function in maintaining pH levels and overexpression in hypoxic conditions, hCA IX inhibitors can be a potential candidate specifically designed to target cancer development at various stages. In search of selective hCA IX inhibitors, we developed a pharmacophore model from the existing natural product inhibitors with IC50 values less than 50 nm. The identified hit molecules were then investigated on protein-ligand interactions using molecular docking experiments followed by molecular dynamics simulations. Among the zinc database 186 hits with an RMSD value less than 1 were obtained, indicating good contact with key residues HIS94, HIS96, HIS119, THR199, and ZN301 required for optimum activity. The top three compounds were subjected to molecular dynamics simulations for 100 ns to know the protein-ligand complex stability. Based on the obtained MD simulation results, binding free energies are calculated. Density Functional Theory (DFT) studies confirmed the energy variation between the Highest Occupied Molecular Orbital (HOMO) and Lowest Unoccupied Molecular Orbital (LUMO). The current study has led to the discovery of lead compounds that show considerable promise as hCA IX inhibitors and suggests that three compounds with special molecular features are more likely to be better-inhibiting hCA IX. Compound S35, characterized by a higher stability margin and a smaller energy gap in quantum studies, is an ideal candidate for selective inhibition of CA IX.

A computational odyssey: uncovering classical β-lactamase inhibitors in dry fruits

In the antibacterial arsenal, β-lactams have held a prominent position, but increasing resistance due to unauthorized use and genetic factors requires new strategies. Combining β-lactamase inhibitors with broad-spectrum β-lactams proves effective in combating this resistance. ESBL producers demand new inhibitors, leading to the exploration of plant-derived secondary metabolites for potent β-lactam antibiotics or alternative inhibitors. Using virtual screening, molecular docking, ADMET analysis, and molecular dynamic simulation, this study actively analyzed the inhibitory activity of figs, cashews, walnuts, and peanuts against SHV-1, NDM-1, KPC-2, and OXA-48 β-lactamases. Using AutoDock Vina, the docking affinities of various compounds for target enzymes were initially screened, revealing 12 bioactive compounds with higher affinities for the target enzymes compared to Avibactam and Tazobactam. Top-scoring metabolites, including Oleanolic acid, Protocatechuic acid, and Tannin, were subjected to MD simulation studies to further analyze the stability of the docked complexes using WebGro. The simulation coordinates, in terms of RMSD, RMSF, SASA, Rg, and hydrogen bonds formed, showed that these phytocompounds are stable enough to retain in the active sites at various orientations. The PCA and FEL analysis also showed the stability of the dynamic motion of Cα residues of phytochemical-bound enzymes. The pharmacokinetic analysis of the top phytochemicals was performed to analyze their bioavailability and toxicity. This study provides new insights into the therapeutic potential of phytochemicals of selected dry fruits and contributes to future experimental studies to identify βL inhibitors from plants.Communicated by Ramaswamy H. Sarma.

Molecular modeling and simulation approaches to characterize potential molecular targets for burdock inulin to instigate protection against autoimmune diseases

In the current study, we utilized molecular modeling and simulation approaches to define putative potential molecular targets for Burdock Inulin, including inflammatory proteins such as iNOS, COX-2, TNF-alpha, IL-6, and IL-1β. Molecular docking results revealed potential interactions and good binding affinity for these targets; however, IL-1β, COX-2, and iNOS were identified as the best targets for Inulin. Molecular simulation-based stability assessment demonstrated that inulin could primarily target iNOS and may also supplementarily target COX-2 and IL-1β during DSS-induced colitis to reduce the role of these inflammatory mechanisms. Furthermore, residual flexibility, hydrogen bonding, and structural packing were reported with uniform trajectories, showing no significant perturbation throughout the simulation. The protein motions within the simulation trajectories were clustered using principal component analysis (PCA). The IL-1β-Inulin complex, approximately 70% of the total motion was attributed to the first three eigenvectors, while the remaining motion was contributed by the remaining eigenvectors. In contrast, for the COX2-Inulin complex, 75% of the total motion was attributed to the eigenvectors. Furthermore, in the iNOS-Inulin complex, the first three eigenvectors contributed to 60% of the total motion. Furthermore, the iNOS-Inulin complex contributed 60% to the total motion through the first three eigenvectors. To explore thermodynamically favorable changes upon mutation, motion mode analysis was carried out. The Free Energy Landscape (FEL) results demonstrated that the IL-1β-Inulin achieved a single conformation with the lowest energy, while COX2-Inulin and iNOS-Inulin exhibited two lowest-energy conformations each. IL-1β-Inulin and COX2-Inulin displayed total binding free energies of - 27.76 kcal/mol and - 37.78 kcal/mol, respectively, while iNOS-Inulin demonstrated the best binding free energy results at - 45.89 kcal/mol. This indicates a stronger pharmacological potential of iNOS than the other two complexes. Thus, further experiments are needed to use inulin to target iNOS and reduce DSS-induced colitis and other autoimmune diseases.

Structural insight into the binding mode of cefotaxime and meropenem to TEM-1, SHV-1, KPC-2, and Amp-C type beta-lactamases

Antimicrobial resistance is an emerging threat to public health around the world. The study employs computational and biophysical methods to investigate the properties of cefotaxime and meropenem's binding to various beta-lactamases like TEM-1, SHV-1, KPC-2, and Amp-C. The enzyme kinetics of purified proteins revealed an increase in Michaelis constant (Km) value in the presence of meropenem and cefotaxime, indicating a decrease in enzyme affinity for nitrocefin. Proteins interact with meropenem/cefotaxime, causing quenching through complex formation. All proteins have one binding site, and binding constant (Kb) values are 104, indicating strong interaction. The study found that meropenem and cefotaxime had high fitness scores for Amp-C, KPC-2,TEM-1 and SHV-1, with binding energy ranging from -7.4 to -7.8, and hydrogen bonds between them. Molecular Dynamic simulation of protein-ligand complexes revealed cefotaxime-binding proteins have slightly lower Root Mean Square Deviation(RMSD) than meropenem-binding proteins, indicating stable association antibiotics with these proteins.

Structure-activity relationship and cytotoxicity of the new thiosemicarbazide derivatives and their Cu(II) complexes against prostate and melanoma cancer cells

In this study, eighteen new ligands (B1-B18) containing a thiosemicarbazide core were synthesized and characterized in terms of physicochemical properties, molecular docking and in vitro biological activity. The structures of eleven ligands were investigated using X-Ray diffraction and Hirschfeld Surface analysis. To study the structure-activity relationship, the organic ligands contained pyridin-2-ylmethyl, pyridin-3-ylmethyl or pyridin-4-ylmethyl moieties and various substituents. Their pharmakokinetic profiles and molecular docking results suggest high potential as new drug candidates. The complexing ability of the selected organic ligands was also evaluated, yielding five new Cu(II) complexes (Cu(B1)Cl2, Cu(B4)Cl2, Cu(B10)Cl2, Cu(B17)Cl2, Cu(B18)Cl2). The obtained results suggest the formation of the polymeric structures. All organic ligands and Cu(II) complexes were tested for anticancer activity against prostate and melanoma cancer cells (PC-3, DU-145, LNCaP, A375, G-361, SK-MEL-28) and normal fibroblasts (BJ), as well as antimicrobial activity against six selected bateria strains. Among B1-B18 compounds, B3, B5, B9, B10, B12 and B14 exhibited cytotoxic activity. The studied Cu(II) complexes were in general more active, with Cu(B1)Cl2 exhibiting antincancer activity agains all three prostate cancer cells and Cu(B10)Cl2 reaching the IC50 value equal to 88 μM against G-361 melanoma cells. Several compounds also exhibited antimicrobial activity against gram-positive and gram-negative bacteria. It was found that the type of specific substituents, especially the presence of -chloro and -dichloro substituents had a greated impact on the cytotoxicity than the position of the nitrogen atom in the pyridylacetyl moiety.

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.

Anticholinesterase activities of novel isoindolin-1,3-dione-based acetohydrazide derivatives: design, synthesis, biological evaluation, molecular dynamic study

In pursuit of developing novel cholinesterase (ChE) inhibitors through molecular hybridization theory, a novel series of isoindolin-1,3-dione-based acetohydrazides (compounds 8a-h) was designed, synthesized, and evaluated as possible acetylcholinesterase (AChE) and butyrylcholinesterase (BChE) inhibitors. In vitro results revealed IC50 values ranging from 0.11 ± 0.05 to 0.86 ± 0.02 µM against AChE and 5.7 ± 0.2 to 30.2 ± 2.8 µM against BChE. A kinetic study was conducted on the most potent compound, 8a, to ascertain its mode of inhibition, revealing its competitive mode against AChE. Furthermore, the binding interaction modes of the most active compound within the AChE active site was elucidated. Molecular dynamics simulations of compound 8a were performed to assess the stability of the 8a-AChE complex. In silico pharmacokinetic predictions for the most potent compounds indicated their potential as promising lead structure for the development of new anti-Alzheimer's disease (anti-AD) agents.

Anti-diabetic drug discovery using the bioactive compounds of Momordica charantia by molecular docking and molecular dynamics analysis

Diabetes mellitus (DM) is a multifactorial life-threatening endocrine disease characterized by abnormalities in glucose metabolism. It is a chronic metabolic disease that involves multiple enzymes such as α-amylase and α-glucosidases. Inhibition of these enzymes has been identified as a promising method for managing diabetes, and researchers are currently focusing on discovering novel α-amylase and α-glucosidase inhibitors for diabetes therapy. Hence, we have selected 12 bioactive compounds from the Momordica charantia (MC) plant and performed a virtual screening and molecular dynamics investigation to identify natural inhibitors of α-amylase and α-glucosidases. Our in silico result revealed that phytocompound Rutin showed the highest binding affinity against α-amylase (1HNY) enzymes at (-11.68 kcal/mol), followed by Karaviloside II (-9.39), Momordicoside F (-9.19), Campesterol (-9.11. While docking against α-glucosidases (4J5T), Rutin again showed the greatest binding affinity (-11.93 kcal/mol), followed by Momordicine (-9.89), and Campesterol (-8.99). Molecular dynamics (MD) simulation research is currently the gold standard for drug design and discovery. Consequently, we conducted simulations of 100 nanoseconds (ns) to assess the stability of protein-ligand complexes based on parameters like RMSD, RMSF, RG, PCA, and FEL. The significance of our findings indicates that rutin from MC might serve as an effective natural therapeutic agent for diabetes management due to its strongest binding affinities with α-amylase and α-glucosidase enzymes. Further research in animals and humans is essential to validate the efficacy of these drug molecules.Communicated by Ramaswamy H. Sarma.

In silico investigation of cholesterol-lowering drugs to find potential inhibitors of dehydrosqualene synthase in Staphylococcus aureus

Staphylococcus aureus is a lethal pathogen that can cause various bacterial infections. This study targets the CrtM enzyme of S. aureus, which is crucial for synthesizing golden carotenoid pigment: staphyloxanthin, which provides anti-oxidant activity to this bacterium for combating antimicrobial resistance inside the host cell. The present investigation quests for human SQS inhibitors against the CrtM enzyme by employing structure-based drug design approaches including induced fit docking (IFD), molecular dynamic (MD) simulations, and binding free energy calculations. Depending upon the docking scores, two compounds, lapaquistat acetate and squalestatin analog 20, were identified as the lead molecules exhibit higher affinity toward the CrtM enzyme. These docked complexes were further subjected to 100 ns MD simulation and several thermodynamics parameters were analyzed. Further, the binding free energies (ΔG) were calculated for each simulated protein-ligand complex to study the stability of molecular contacts using the MM-GBSA approach. Pre-ADMET analysis was conducted for systematic evaluation of physicochemical and medicinal chemistry properties of these compounds. The above study suggested that lapaquistat acetate and squalestatin analog 20 can be selected as potential lead candidates with promising binding affinity for the S. aureus CrtM enzyme. This study might provide insights into the discovery of potential drug candidates for S. aureus with a high therapeutic index.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13205-023-03862-y.

Computational investigation of impact of Pb(II) and Ni(II) ions on hUNG enzyme: insights from molecular dynamics simulations

Human uracil DNA glycosylase (hUNG), a crucial player in the initiation of the base excision repair pathway, is susceptible to alterations in function and conformation induced by the accumulation of toxic metals. Despite the recognized impact of toxic metals on DNA repair enzymes, there exists a notable deficiency in theoretical investigations addressing this phenomenon. This study investigates the impact of toxic heavy metal ions, Pb(II) and Ni(II), on the stability of hUNG through molecular dynamics (MD) simulations. The initial analysis involved the identification of key cavities in the hUNG enzyme. Notably, the active site cavity emerged as a promising site for ligand binding. Subsequently, AutoDockTools software was employed to dock Pb(II) and Ni(II) onto the identified cavities, followed by extensive MD simulations. The MD analysis, encompassing parameters such as root mean square deviation, radius of gyration, solvent accessible surface area, hydrogen bond variations, Ramachandran plot, principal component analysis, and root mean square fluctuations, collectively revealed distinct alterations in the behavior of the enzyme upon complexation with Pb(II) and Ni(II). Interestingly, the enzyme exhibited enhanced structural stability, reduced flexibility, and modified hydrogen bonding patterns in the presence of these toxic metal ions. The observed limitation in structural flexibility implies a more rigid and stable conformation when the enzyme complex with Pb(II) and Ni(II) compared to its free form. This structural alteration may lead to a potential reduction in enzymatic activity, suggesting that toxic metal ions influence the functional dynamics of hUNG. These computational findings offer valuable insights into the molecular interactions between metal ions and enzymes.Communicated by Ramaswamy H. Sarma.

Novel Leflunomide Analog, UTLOH-4e, Ameliorates Gouty Arthritis Induced by Monosodium Urate Via NF-κB/NLRP3 Signaling Pathway

BACKGROUND

Gouty arthritis (GA) is a common form of inflammatory arthritis caused by intra-articular deposition of monosodium urate (MSU) crystals; however, there is a tremendous lack of safe and effective therapy in the clinic.

OBJECTIVE

The goal of this work was to investigate a novel leflunomide analogue, N-(2,4- dihydroxyphenyl)-5-methyl-1,2-oxazole-3-carboxamide (UTLOH-4e), for its potential to prevent/ treat gouty arthritis.

METHODS

In this study, the anti-inflammatory activity of UTLOH-4e was evaluated by MSUinduced GA model in vivo and in vitro, and the molecular docking test was applied to estimate the affinity of UTLOH-4e/UTL-5g/b for MAPKs, NF-κB, and NLRP3.

RESULTS

In vitro, UTLOH-4e (1~100 μM) treatment inhibited the inflammatory reaction with no obvious cytotoxicity in PMA-induced THP-1 macrophages exposed to MSU crystals for 24 h, involving the prominent decreased production and gene expression of IL-1β, TNF-α, and IL-6. Western blot analyses demonstrated that UTLOH-4e (1~100 μM) significantly suppressed the activation of NLRP3 inflammasomes, NF-κB, and MAPK pathways. Furthermore, the data from the experiment on gouty rats induced by intra-articular injection of MSU crystal confirmed that UTLOH-4e markedly ameliorated rat paw swelling, articular synovium inflammation and reduced the concentration of IL-1β and TNF-α in serum through down-regulating NLRP3 protein expression.

CONCLUSION

These results manifested that UTLOH-4e ameliorates GA induced by MSU crystals, which contributes to the modulation of NF-κB/ NLRP3 signaling pathway, suggesting that UTLOH- 4e is a promising and potent drug candidate for the prevention and treatment of gouty arthritis.

Deciphering insights into the binding mechanism and plasticity of Telacebec with M. tuberculosis cytochrome bcc-aa3 supercomplex through an unbiased molecular dynamics simulation, free-energy analysis, and DFT study

The cytochrome bcc-aa3 supercomplex, a key component in the electron transport chain pathway involved in bacterial energy production and homeostasis, is a clinically validated target for tuberculosis (TB), leading to Telacebec (Q203). Telacebec is a potent candidate drug under Phase II clinical development for the treatment of drug-sensitive and drug-resistant TB. Recently, the cryo-electron microscopy structure of this supercomplex from Mycobacterium tuberculosis (Mtb) complexed with Q203 was resolved at 6.9 Å resolution (PDB ID: 7E1W). To understand the binding site (QP site) flexibility and Q203's stability at the QP site of the Mtb cytochrome bcc complex, we conducted molecular dynamics (MD) simulation and free energy analysis on this complex in an explicit hydrated lipid bilayer environment for 500 ns. Through this study, the persistence of a range of direct and indirect interactions was observed over the course of the simulation. The significance of the interactions with His375, Tyr161, Ala178, Ala179, Ile183, His355, Leu356, and Thr313 is underlined. Electrostatic energy was the primary source of the net binding free energy, regardless of the important interacting residues. The overall binding free energy for Q203 was -112.84 ± 7.73 kcal/mol, of which the electrostatic and lipophilic energy contributions were -116.31 ± 1.14 and -21.32 ± 2.35 kcal/mol, respectively. Meanwhile, DFT calculations were utilized to elucidate Q203's molecular properties. Overall, this study deciphers key insights into the cytochrome bcc-aa3 supercomplex with Q203 on the ground of molecular mechanics and quantum mechanics that may facilitate structure-based drug design and optimization for the discovery of the next-generation antitubercular drug(s).Communicated by Ramaswamy H. Sarma.

Exploring the role of cyclodextrins as a cholesterol scavenger: a molecular dynamics investigation of conformational changes and thermodynamics

This study presents a comprehensive analysis of the cholesterol binding mechanism and conformational changes in cyclodextrin (CD) carriers, namely βCD, 2HPβCD, and MβCD. The results revealed that the binding of cholesterol to CDs was spontaneous and thermodynamically favorable, with van der Waals interactions playing a dominant role, while Coulombic interactions have a negligible contribution. The solubility of cholesterol/βCD and cholesterol/MβCD complexes was lower compared to cholesterol/2HPβCD complex due to stronger vdW and Coulombic repulsion between water and CDs. Hydrogen bonding was found to have a minor role in the binding process. The investigation of mechanisms and kinetics of binding demonstrated that cholesterol permeates into the CD cavities completely. Replicas consideration indicated that while the binding to 2HPβCD occurred perpendicularly and solely through positioning cholesterol's oxygen toward the primary hydroxyl rim (PHR), the mechanism of cholesterol binding to βCD and MβCD could take place with the orientation of oxygen towards both rims. Functionalization resulted in decreased cavity polarity, increased constriction tendency, and altered solubility and configuration of the carrier. Upon cholesterol binding, the CDs expanded, increasing the cavity volume in cholesterol-containing systems. The effects of cholesterol on the relative shape anisotropy (κ2) and asphericity parameter (b) in cyclodextrins were investigated. βCD exhibited a spherical structure regardless of cholesterol presence, while 2HPβCD and MβCD displayed more pronounced non-sphericity in the absence of cholesterol. Loading cholesterol transformed 2HPβCD and MβCD into more spherical shapes, with increased probabilities of higher κ2. MβCD showed a higher maximum peak of κ2 compared to 2HPβCD after cholesterol loading, while 2HPβCD maintained a significant maximum peak at 0.2 for b.

Lignin-ultrasound method: Enhancement of antimicrobial capacity of MoS2-containing films

To improve the antimicrobial ability of MoS2-containing films, we used lignin and triple-frequency ultrasound for liquid-phase exfoliation (LPE) to obtain MoS2 nanosheets. Photoresponsive antimicrobial films with MoS2 nanosheets, lignin, polyvinyl alcohol and deep eutectic solvents were subsequently prepared. Lignin functionalized the MoS2 nanosheets by chemically linking with S in MoS2 and significantly improved the exfoliation efficiency. Tri-frequency ultrasound produces beneficial effects on the LPE process by creating a more homogeneous sound field and a stronger degree of cavitation. The concentration of MoS2 nanosheets in the exfoliating solution could reach 1.713 mg/mL under the effect of lignin-ultrasound. The antimicrobial ability of the films was analyzed, and the colony-forming units of E. coli and S. aureus could be reduced from 7 × 106 to 1 × 106 cfu/mL under the irradiation of infrared. The lignin in the film undergoes depolymerization and demethoxylation under the irradiation of infrared to have a more phenolic hydroxyl structure, which confers the growth inhibition ability of the films for bacteria that cannot be in close contact with the film. The method we used has some significance for the preparation of MoS2 nanosheets, and composite films prepared from MoS2, and lignin can be used in food packaging, wound antimicrobials, and other fields.

Temperature dependent molecular dynamics simulation study to understand the stabilizing effect of NADES on the protein β-Lactoglobulin

The thermal stability of a protein is an important concern for its practical application in food processing industries. In this study, we have carried out classical molecular dynamics simulations to systematically investigate the effect of NADES (natural deep eutectic solvent) on the stabilization of the protein β-Lactoglobulin (BLG) at different temperatures. This study sheds light on the very aspects of NADES composed of betaine and sorbitol on the stability of the protein. NADES provides better stability to the protein up to a temperature of 400 K than in water. It is observed that the protein starts to unfold above temperature 400 K in spite of the presence of NADES which is quiet evident from the root mean square deviation (RMSD) and radius of gyration (Rg) plots. The decreasing average solvent accessible surface area (SASA) values and increasing intra-protein hydrogen bonds indicate better stability of the protein in NADES medium than in water at temperatures 300 K and 400 K. At high temperatures viz. 450 K and 500 K the number and distribution of solvent species (betaine and sorbitol) around the protein surface show an increment that are evident from the calculations of solvation shell, radial and spatial distribution functions. Increased number of betaine molecules that interact with the protein through electrostatic interaction may lead to destabilization of the protein at these temperatures. This study suggests that NADES could be used as an ideal medium for thermal stability of the protein BLG up to a temperature of 400 K. Beyond this temperature, NADES used for this study fails to exert stabilization effect on the protein.

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.

Anti-quorum Sensing and Anti-biofilm Effect of Nocardiopsis synnemataformans RMN 4 (MN061002) Compound 2,6-Di-tert-butyl, 1,4-Benzoquinone Against Biofilm-Producing Bacteria

In this study, the anti-biofilm compound of 2,6-Di-tert-butyl, 1,4-benzoquinone was purified from Nocardiopsis synnemataformans (N. synnemataformans) RMN 4 (MN061002). To confirm the compound, various spectroscopy analyses were done including ultraviolet (UV) spectrometer, Fourier transform infrared spectroscopy (FTIR), analytical high-performance liquid chromatography (HPLC), preparative HPLC, gas chromatography-mass spectroscopy (GC-MS), liquid chromatography-mass spectroscopy (LC-MS), and 2D nuclear magnetic resonance (NMR). Furthermore, the purified compound was shown 94% inhibition against biofilm-producing Proteus mirabilis (P. mirabilis) (MN396686) at 70 µg/mL concentrations. Furthermore, the metabolic activity, exopolysaccharide damage, and hydrophobicity degradation results of identified compound exhibited excellent inhibition at 100 µg/mL concentration. Furthermore, the confocal laser scanning electron microscope (CLSM) and scanning electron microscope (SEM) results were shown with intracellular damages and architectural changes in bacteria. Consecutively, the in vivo toxicity effect of the compound against Artemia franciscana (A. franciscana) was shown to have a low mortality rate at 100 µg/mL. Finally, the molecular docking interaction between the quorum sensing (QS) genes and identified compound clearly suggested that the identified compound 2,6-Di-tert-butyl, 1,4-benzoquinone has anti-quorum sensing and anti-biofilm activities against P. mirabilis (MN396686).

The inhibitory mechanism of echinacoside against Staphylococcus aureus Ser/Thr phosphatase Stp1 by virtual screening and molecular modeling

CONTEXT

Stp1 is a new potential target closely related to the pathogenicity of Staphylococcus aureus (S. aureus). In this study, effective Stp1 inhibitors were screened via virtual screening and enzyme activity experiments, and the inhibition mechanism was analyzed using molecular dynamics simulation.

METHODS

AutoDock Vina 4.0 software was used for virtual screening. The molecular structures of Stp1 and ligands were obtained from the RCSB Protein Data Bank and Zinc database, respectively. The molecular dynamics simulation used the Gromacs 4.5.5 software package with the Amberff99sb force field and TIP3P water model. AutoDock Tools was used to add polar hydrogen atoms to Stp1 and distribute part of the charge generated by Kollman's combined atoms. The binding free energies were calculated using the Amber 10 package.

RESULTS

The theoretical calculation results are consistent with the experimental results. We found that echinacoside (ECH) substantially inhibits the hydrolytic activity of Stp1. ECH competes with the substrate by binding to the active center of Stp1, resulting in a decrease in Stp1 activity. In addition, Met39, Gly41, Asp120, Asn162, and Ile163 were identified to play key roles in the binding of Stp1 to ECH. The benzene ring of ECH also plays an important role in complex binding. These findings provide a robust foundation for the development of innovative anti-infection drugs.

Identification of coastal pesticide pollutants as potent inhibitors of Bacillus pasteurii urease mediated calcium carbonate precipitation: a computational approach

Microbially induced calcite precipitation (MICP) through urease enzyme has attained a lot of recognition in various fields of civil engineering and geotechnology for stabilizing the strength of soil and various concrete materials. The activity of urease has been found to be affected by various factors like temperature, substrate concentrations, pH of the medium, presence of inhibitors, etc. Through this study, the outcome of the interaction of pesticides (commonly found in Indian coastal regions) on Bacillus pasteurii urease, a major organism reported for MICP studies has been investigated in silico. The results from the study revealed that the enzyme has higher interactions of -4.1, -3.2, and -3.4 kJ/mol with common pesticides like dichloro diphenyl dichloro ethane(DDD), dichloro diphenyl trichloroe thane (DDT), and methyl parathion of organochlorides and organophosphates class. From the molecular dynamics simulation analysis, complex 1 (DDD -receptor) has been found to have the highest and more compact structure followed by methyl parathion -receptor. Prime MM-GBSA analysis also revealed the highest binding energy of -27.8 kcal/mol with the protein and DDD. Thus, it can be inferred from the current study that pesticides, particularly, DDD, DDT, and methyl parathion present in the coastal areas may have an impact on urease. This interaction can result in the inhibition of the urease activity of B. pasteurii, thus preventing the biomineralization process. This study would be the first report on the computational approach to understanding the interaction of prominent pesticides on the coastal region and B. pasteurii urease.Communicated by Ramaswamy H. Sarma.

In silico prediction of mutation sites for anthranilate synthase from Serratia marcesens to deregulate tryptophan feedback inhibition

Allosteric feedback inhibition of the committed step in amino acid biosynthetic pathways is a major concern for production of amino acids at industrial scale. Anthranilate synthase (AS) catalyzes the first reaction of tryptophan biosynthetic pathway found in microorganisms and is feedback inhibited by its own product i.e. tryptophan. Here, we identified new mutant sites in AS using computational mutagenesis approach. MD simulations (20 ns) followed by MMPBSA and per residue decomposition energy analysis identified seven amino acid residues with best binding affinity for tryptophan. All 19 mutant structures were generated for each identified amino acid residue followed by simulation to evaluate effect of mutation on protein stability. Later, molecular docking studies were employed to generate mutant-tryptophan complex and structures with binding energies (kcal/mol) much higher than wild-type AS were selected. Finally, two mutants i.e., S37W and S37H were identified on the basis of positive binding scores and loss of tryptophan binding inside pocket. Further, MD simulations run for 200 ns were performed over these mutant-tryptophan complexes followed by RMSD, RMSF, radius of gyration , solvent accessible surface area , intra-protein hydrogen bond numbers, principal component analysis, free energy landscape (FEL) and secondary structure analysis to rationale effect of mutations on stability of protein. Cross correlation analysis of mutant site amino acids (S37W) with key residues of catalytic site (G325, T326, H395 and G482) was done to evaluate the effect of mutations on catalytic site conformation. Current computational mutagenesis approach predicted two mutants S37W and S37H with proposed deregulated feedback inhibition by tryptophan and retained catalytic activity.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.

Identification and evaluation of potential inhibitor molecules against TcyA from Candidatus Liberibacter asiaticus

Of the two putative amino acid binding periplasmic receptors of ABC transporter family in Candidatus Liberibacter asiaticus (CLas), cystine binding receptor (CLasTcyA) has been shown to mainly express in phloem of citrus plant and is a target for inhibitor development. The crystal structure of CLasTcyA in complex with substrates has been reported earlier. The present work reports the identification and evaluation of potential candidates for their inhibitory potential against CLasTcyA. Among many compounds, selected through virtual screening, and MD simulation, pimozide, clidinium, sulfasalazine and folic acid showed significantly higher affinities and stability in complex with CLasTcyA. The SPR studies with CLasTcyA revealed significantly higher binding affinities for pimozide and clidinium (Kd, 2.73 nM and 70 nM, respectively) as compared to cystine (Kd, 1.26 µM). The higher binding affinities could be attributed to significantly increased number of interactions in the binding pocket as evident from the crystal structures of CLasTcyA in complex with pimozide and clidinium as compared to cystine. The CLasTcyA possess relatively large binding pocket where bulkier inhibitors fit quite well. In planta studies, carried out to assess the effect of inhibitors on HLB infected Mosambi plants, showed significant reduction in CLas titre in plants treated with inhibitors as compared to control plants. The results showed that pimozide exhibited higher efficiency as compared to clidinium in reducing CLas titre in treated plants. Our results showed that the inhibitor development against critical proteins like CLasTcyA can be an important strategy in management of HLB.

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.

Molecular insights into the inhibitory potential of anthocyanidins on glucokinase regulatory protein

Computational methods were used to investigate six anthocyanidins exhibiting antidiabetic activity by inhibiting glucokinase regulatory protein (GKRP) activity. Density functional theory was used to optimise the geometry of anthocyanidins and calculate their quantum chemical properties. A blind docking method was employed to conduct a molecular docking study, which revealed that delphinidin (Del), cyanidin (Cya), and pelargonidin (Pel) as potential GKRP inhibitors with the lowest binding free energy of -8.7, -8.6, and -8.6 kcal/mol, corresponding to high binding affinity. The molecular dynamics study further verified the blind docking results by showing high GKRP-F1P complex stability and high binding affinity calculated through the MM/GBSA method, upon the binding of pelargonidin. The lower RMSF values of pivotal GK-interacting residues for GKRP-F1P-Pel compared to GKRP-F1P, as a positive control, indicating pelargonidin ability to maintain the inactive conformation of GKRP through the inhibition of GK binding. The key residues that control the binding of the F1P to GKRP and anthocyanidin to GKRP-F1P were also identified in this study. Altogether, pelargonidin is anthocyanidins-derived natural products that have the most potential to act as inhibitors of GKRP and as antidiabetic nutraceuticals.

Tanshinone-I for the treatment of uterine fibroids: Molecular docking, simulation, and density functional theory investigations

Uterine fibroids (UF), most prevalent gynecological disorder, require surgery when symptomatic. It is estimated that between 25 and 35 percent of women wait until the symptoms have worsened like extended heavy menstrual bleeding and severe pelvic pain. These UF may be reduced in size through various methods such as medical or surgical intervention. Progesterone (prog) is a crucial hormone that restores the endometrium and controls uterine function. In the current study, 28 plant-based molecules are identified from previous literature and docked onto the prog receptors with 1E3K and 2OVH. Tanshinone-I has shown the best docking score against both proteins. The synthetic prog inhibitor Norethindrone Acetate is used as a standard to evaluate the docking outcomes. The best compound, tanshinone-I, was analyzed using molecular modeling and DFT. The RMSD for the 1E3K protein-ligand complex ranged from 0.10 to 0.42 Å, with an average of 0.21 Å and a standard deviation (SD) of 0.06, while the RMSD for the 2OVH protein-ligand complex ranged from 0.08 to 0.42 Å, with an average of 0.20 Å and a SD of 0.06 showing stable interaction. In principal component analysis, the observed eigen values of HPR-Tanshinone-I fluctuate between -1.11 to 1.48 and -1.07 to 1.25 for PC1 and PC2, respectively (1E3K), and the prog-tanshinone-I complex shows eigen values of -38.88 to -31.32 and -31.32 to 35.87 for PC1 and PC2, respectively (2OVH), which shows Tanshinone-I forms a stable protein-ligand complex with 1E3K in comparison to 2OVH. The Free Energy Landscape (FEL) analysis shows the Gibbs free energy in the range of 0 to 8 kJ/mol for Tanshinone-I with 1E3K and 0 to 14 kJ/mol for Tanshinone-I with the 2OVH complex. The DFT calculation reveals ΔE value of 2.8070 eV shows tanshinone-I as a stable compound. 1E3K modulates the prog pathway, it may have either an agonistic or antagonistic effect on hPRs. Tanshinone-I can cause ROS, apoptosis, autophagy (p62 accumulation), up-regulation of inositol requiring protein-1, enhancer-binding protein homologous protein, p-c-Jun N-terminal kinase (p-JNK), and suppression of MMPs. Bcl-2 expression can change LC3I to LC3II and cause apoptosis through Beclin-1 expression.

Collaterally Sensitive β-Lactam Drugs as an Effective Therapy against the Pre-Existing Methicillin Resistant Staphylococcus aureus (MRSA) Biofilms

Methicillin-resistant Staphylococcus aureus (MRSA) is among the leading causes of nosocomial infections and forms biofilms, which are difficult to eradicate because of their increasing resistance to antimicrobial agents. This is especially true for pre-existing biofilms. The current study focused on evaluating the efficacy of three β-lactam drugs, meropenem, piperacillin, and tazobactam, alone and in combination against the MRSA biofilms. When used individually, none of the drugs exhibited significant antibacterial activity against MRSA in a planktonic state. At the same time, the combination of meropenem, piperacillin, and tazobactam showed a 41.7 and 41.3% reduction in planktonic bacterial cell growth, respectively. These drugs were further assessed for biofilm inhibition and removal. The combination of meropenem, piperacillin, and tazobactam caused 44.3% biofilm inhibition, while the rest of the combinations did not show any significant effects. Results also revealed that piperacillin and tazobactam exhibited the best synergy against the pre-formed biofilm of MRSA, with 46% removal. However, adding meropenem to the piperacillin and tazobactam combination showed a slightly reduced activity towards the pre-formed biofilm of MRSA and removed 38.7% of it. Although the mechanism of synergism is not fully understood, our findings suggest that these three β-lactam drugs can be used in combination as very effective therapeutic agents for the treatment of pre-existing MRSA biofilms. The in vivo experiments on the antibiofilm activity of these drugs will pave the way for applying such synergistic combinations to clinics.

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.

New insights into the inhibitory effect of phenol carboxylic acid antioxidants on mushroom tyrosinase by molecular dynamic studies and experimental assessment

The inhibitory effects of ferulic and chlorogenic acids on tyrosinase activity were investigated through multi-spectroscopic and molecular docking techniques. Ferulic and chlorogenic acids, flavonoid compounds, demonstrated inhibitory monophenolase activities of tyrosinase. The inhibitor effects against monophenolase activity were in a reversible and competitive manner with ki value equal to 6.8 and 7.5 µM respectively. The affinity between tyrosinase and L-DOPA decreased when fatty acids were added to the solution. The multi-spectroscopic techniques like UV-vis, fluorescence, and isothermal calorimetry are employed to investigate changes. Intrinsic fluorescence quenching and conformational changes of tyrosinase by hydrophobic interaction were confirmed. Tyrosinase had two and three binding sites for ferulic and chlorogenic acids with a binding constant in the order of magnitude of -6.8 and -7.2 kcal/mol. In addition, the secondary structural changes with Circular dichroism (CD) analysis, secondary structure (DSSP), radius of gyration (Rg) and analysis of hydrogen bonds (H-bonds) confirmed. Ferulic acid effect can be observed obviously and also content of α-helix decreased. Thermodynamic parameters indicated that the interaction between enzyme and ferulic and chlorogenic acids followed a spontaneous reaction dynamic manner with ΔG = -14.78 kJ/mol and ΔG = -14.61 kJ/mol (298k). The findings highlighted the potential applications of ferulic acid and chlorogenic acids in food and drug industries as potent inhibitors of tyrosinase.Communicated by Ramaswamy H. Sarma.

Identification of Four New Chemical Series of Small Drug-Like Natural Products as Potential Neuropilin-1 Inhibitors by Structure-Based Virtual Screening: Pharmacophore-Based Molecular Docking and Dynamics Simulation

Neuropilin-1 (NRP-1), a surface transmembrane glycoprotein, is one of the most important co-receptors of VEGF-A165 (vascular endothelial growth factor) responsible for pathological angiogenesis. In general, NRP-1 overexpression in cancer correlates with poor prognosis and more tumor aggressiveness. NRP-1 role in cancer has been mainly explained by mediating VEGF-A165-induced effects on tumor angiogenesis. NRP-1 was recently identified as a co-receptor and an independent gateway for SARS-CoV-2 through binding subunit S2 of Spike protein in the same way as VEGF-A165. Thus, NRP-1 is of particular value as a target for cancer therapy and other angiogenesis-dependent diseases as well as for SARS-CoV-2 antiviral intervention. Herein, The Super Natural II, the largest available database of natural products (∼0.33 M), pre-filtered with drug-likeness criteria (absorption, distribution, metabolism and excretion/toxicity), was screened against NRP-1. NRP-1/VEGF-A165 interaction is one of protein-protein interfaces (PPIs) known to be challenging when approached in-silico. Thus, a PPI-suited multi-step virtual screening protocol, incorporating a derived pharmacophore with molecular docking and followed by MD (molecular dynamics) simulation, was designed. Two stages of pharmacophorically constrained molecular docking (standard and extra precisions), a mixed Torsional/Low-mode conformational search and MM-GBSA ΔG binding affinities calculation, resulted in the selection of 100 hits. These 100 hits were subjected to 20 ns MD simulation, that was extended to 100 ns for top hits (20) and followed by post-dynamics analysis (atomic ligand-protein contacts, RMSD, RMSF, MM-GBSA ΔG, Rg, SASA and H-bonds). Post-MD analysis showed that 19 small drug-like nonpeptide natural molecules, grouped in four chemical scaffolds (purine, thiazole, tetrahydropyrimidine and dihydroxyphenyl), well verified the derived pharmacophore and formed stable and compact complexes with NRP-1. The discovered molecules are promising and can serve as a base for further development of new NRP-1 inhibitors.

Computational investigation of potent inhibitors against YsxC: structure-based pharmacophore modeling, molecular docking, molecular dynamics, and binding free energy

In S. aureus, ribosome biogenesis GTP-binding (YsxC), a GTPase interacts with 50S subunit and 30S subunit of ribosome, and β' subunit of RNA polymerase and played an important role in protein synthesis. For the identification of potent lead molecules, we have conducted pharmacophore modeling by consideration of pharmacophore features of GTP among YsxC-GTP complex. Virtual screening and molecular docking results displayed that five pharmacokinetic and ADMET filtered molecules-ZINC000006424138, ZINC000095502032, ZINC000225415132, ZINC000095475800, and ZINC000012990761-had higher binding affinities than GTP with YsxC. All the identified molecules shared similar pharmacophore features of GTP and were stabilized via hydrogen bonds and hydrophobic interactions with YsxC. Molecular dynamics analysis revealed that YsxC-inhibitor(s) complexes were lesser dynamics and higher stable than YsxC-GTP complex. Molecular Mechanics/Poisson-Boltzmann Surface Area (MMPBSA) results confirmed that identified molecules bound at the active site (Arg33, Ser34, Asn35, Val36, Lys38, Ser39, Thr40, Thr54, Ser55, Pro58, Lys60, Thr61, Thr144, Lys145, Ser178, and Ile179) of YsxC and formed the lower energy (-190.32 ± 3.46 to -217.03 ± 2.55 kJ/mol) complexes than YsxC-GTP (-157.16 ± 2.89 kJ/mol) complex. The identified molecules in this study can be further tested and utilized to design novel antimicrobial agents for S. aureus.Communicated by Ramaswamy H. Sarma.

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.

Molecular modeling and dynamics simulation of alcohol dehydrogenase enzyme from high efficacy cellulosic ethanol-producing yeast mutant strain Pichia kudriavzevii BGY1-γm

One of the major constraints limiting the use of abundantly available lignocellulosic biomass as potential feedstock for alcohol industry is the lack of C₆/C₅ co-sugar fermenting yeast. The present study explores a mutant yeast Pichia kudriavzevii BGY1-γm as a potential strain for bioconversion of glucose/xylose sugars of green biomass into ethanol under batch fermentation. The mutant strain having higher alcohol dehydrogenase activity (11.31%) showed significantly higher ethanol concentration during co-fermentation of glucose/xylose sugars (14.2%) as compared to the native strain. Based on 99% sequence similarity of ADH encoding gene from the mutant with the gene sequences from other yeast strains, the ADH enzyme was identified as ADH-1 type. The study reveals first three-dimensional model of ADH-1 utilizing glucose/xylose sugars from P. kudriavzevii BGY1-γm (PkADH mutant). The refined and validated model of PkADH mutant was used for molecular docking against the substrate (acetaldehyde) and product (ethanol). Molecular docking results showed that substrate and product exhibited a binding affinity of -4.55 and -4.5 kcal/mol with PkADH mutant. Acetaldehyde and ethanol interacted at the active site of PkADH mutant via hydrogen bonds (Ser42, His69 and Asp163) and hydrophobic interactions (Cys40, Ser42, His69, Cys95, Trp123 and Asp163) to form the stable protein-ligand complex. Molecular dynamics analysis revealed that PkADH-mutant acetaldehyde and PkADH-mutant ethanol complexes were more stable than PkADH mutant. MMPBSA binding energy confirmed that binding of substrate and product results in the formation of a lower energy stable protein-ligand complex.Communicated by Ramaswamy H. Sarma.

Identification of potential inhibitors for LLM of Staphylococcus aureus: structure-based pharmacophore modeling, molecular dynamics, and binding free energy studies

Staphylococcus aureus causes various life-threatening diseases in humans and developed resistance to several antibiotics. Lipophilic membrane (LLM) protein regulates bacterial lysis rate and methicillin resistance level in S. aureus. To identify potential lead molecules, we performed a structure-based pharmacophore modeling by consideration of pharmacophore properties from LLM-tunicamycin complex. Further, virtual screening of ZINC database against the LLM was conducted and compounds were assessed for Lipinski and ADMET properties. Based on pharmacokinetic, and molecular docking, five potential inhibitors (ZINC000072380005, ZINC000257219974, ZINC000176045471, ZINC000035296288, and ZINC000008789934) were identified. Molecular dynamics simulation (MDS) of these five molecules was performed to evaluate the dynamics and stability of protein after binding of the ligands. Several MDS analysis like RMSD, RMSF, Rg, SASA, and PCA confirm that identified compounds exhibit higher binding affinity as compared to tunicamycin for LLM. The binding free energy analysis reveals that five compounds exhibit higher binding energy in the range of -218.76 to -159.52 kJ/mol, which is higher as compared to tunicamycin (-116.13 kJ/mol). Individual residue decomposition analysis concludes that Asn148, Asp151, Asp208, His271, and His272 of LLM play a significant role in the formation of lower energy LLM-inhibitor(s) complexes. These predicted molecules displayed pharmacological and structural properties and may be further used to develop novel antimicrobial compounds against S. aureus.Communicated by Ramaswamy H. Sarma.

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.

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

Klebsiella pneumonia is known to cause several nosocomial infections in immunocompromised patients. It has developed resistance against a broad range of presently available antibiotics, resulting in high mortality rates in patients and declared an urgent threat. Therefore, exploration of possible novel drug targets against this opportunistic bacteria needs to be undertaken. In the present study, we performed an extensive in-silico analysis for functional and structural annotation and characterized HP CP995_08280 from K. pneumonia as a drug target and aimed to identify potent drug candidates. The functional and structural studies using several bioinformatics tools and databases predicted that HP CP995_08280 is a cytosolic protein that belongs to the β-lactamase family and shares structural similarity with FmtA protein from Staphylococcus aureus (PDB ID: 5ZH8). The structure of HP CP995_08280 was successfully modeled followed by structure-based virtual screening, docking, molecular dynamics, and Molecular mechanic/Poisson-Boltzmann surface area (MMPBSA) were performed to identify the potential compounds. We have found five potent antibacterial molecules, namely BDD 24083171, BDD 24085737, BDE 25098678, BDE 33638819, and BDE 33672484, which exhibited high binding affinity (>-7.5 kcal/mol) and were stabilized by hydrogen bonding and hydrophobic interactions with active site residues (Ser42, Lys45, Tyr126, and Asp128) of protein. Molecular dynamics and MMPBSA revealed that HP CP995_08280 - ligand(s) complexes were less dynamic and more stable than native HP CP995_08280. Hence, the present study may serve as a potential lead for developing inhibitors against drug-resistant Klebsiella pneumonia.

A comprehensive review on genomics, systems biology and structural biology approaches for combating antimicrobial resistance in ESKAPE pathogens: computational tools and recent advancements

In recent decades, antimicrobial resistance has been augmented as a global concern to public health owing to the global spread of multidrug-resistant strains from different ESKAPE pathogens. This alarming trend and the lack of new antibiotics with novel modes of action in the pipeline necessitate the development of non-antibiotic ways to treat illnesses caused by these isolates. In molecular biology, computational approaches have become crucial tools, particularly in one of the most challenging areas of multidrug resistance. The rapid advancements in bioinformatics have led to a plethora of computational approaches involving genomics, systems biology, and structural biology currently gaining momentum among molecular biologists since they can be useful and provide valuable information on the complex mechanisms of AMR research in ESKAPE pathogens. These computational approaches would be helpful in elucidating the AMR mechanisms, identifying important hub genes/proteins, and their promising targets together with their interactions with important drug targets, which is a crucial step in drug discovery. Therefore, the present review aims to provide holistic information on currently employed bioinformatic tools and their application in the discovery of multifunctional novel therapeutic drugs to combat the current problem of AMR in ESKAPE pathogens. The review also summarizes the recent advancement in the AMR research in ESKAPE pathogens utilizing the in silico approaches.

Identification of bio-active food compounds as potential SARS-CoV-2 PLpro inhibitors-modulators via negative image-based screening and computational simulations

Despite significant studies on the COVID-19 pandemic, scientists around the world are still battling to find a definitive therapy against the ongoing severe global health crisis. In this study, advanced computational approaches have been employed to identify bioactive food constituents as potential SARS-CoV-2 PLpro inhibitors-modulators. As a validated antiviral drug target, PLpro has gained tremendous attention for therapeutics developments. Therefore, targeting the multifunctional SARS-CoV-2 PLpro protein, ∼1039 bioactive dietary compounds have been screened extensively through novel techniques like negative image-based (NIB) screening and molecular docking approaches. In particular, the three different models of NIB screening have been generated and used to re-score the dietary compounds based on the negative image which is created by reversing the shape and electrostatics features of PLpro protein's ligand-binding cavity. Further, 100 ns molecular dynamics simulation has been performed and MM-GBSA based binding free energies have been estimated for the final proposed four dietary compounds (PC000550, PC000361, PC000558, and PC000573) as potential inhibitors/modulators of SARS-CoV-2 PLpro protein. Employed computational study outcome also has been compared with respect to the earlier experimentally investigated compound GRL0617 against SARS-CoV-2 PLpro protein, which suggests much greater interaction potential in terms of binding affinity and other energetic contributions for the proposed dietary compounds. Hence, the present study suggests that proposed dietary compounds can be suitable chemical entities for modulating the activity of PLpro protein or can be further utilized for optimizing or screening of novel chemical surrogates, however also needs experimental evaluation for entry in clinical studies for better assessment.

Promising antivirals for PLpro of SARS-CoV-2 using virtual screening, molecular docking, dynamics, and MMPBSA

The recent pandemic of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection (COVID-19) is a viral respiratory disease that has been spread all over the globe. Therefore, it is an urgent requirement to identify and develop drugs for this contagious infection. The papain-like protease (PLpro) of SARS-CoV-2 performs critical functions in virus replication and immune evasion, making it an enticing therapeutic target. SARS-CoV-2 and SARS-CoV PLpro proteases have significant similarities, and an inhibitor discovered for SARS-CoV PLpro is an exciting first step toward therapeutic development. Here, a set of antiviral molecules were screened at the catalytic and S-binding allosteric sites of papain-like protease (PLpro). Molecular docking results suggested that five molecules (44560613, 136277567, S5652, SC75741, and S3833) had good binding affinities at both sites of PLpro. Molecular dynamics analysis like root mean square deviation (RMSD), root mean square fluctuation (RMSF), radius of gyration (Rg), solvent accessible surface area (SASA), and hydrogen bond results showed that identified molecules with PLpro tend to form stable PLpro-inhibitor(s) complexes. Molecular Mechanics/Position-Boltzmann Surface Area (MMPBSA) analysis confirmed that antiviral molecules bound PLpro complex had lower energy (-184.72 ± 7.81 to -215.67 ± 6.73 kJ/mol) complexes. Noticeably, computational approaches revealed promising antivirals candidates for PLpro, which may be further tested by biochemical and cell-based assays to assess their potential for SARS-CoV-2.Communicated by Ramaswamy H. Sarma.

Quantum Mechanics/Molecular Mechanics Studies on the Catalytic Mechanism of a Novel Esterase (FmtA) of Staphylococcus aureus

FmtA is a novel esterase that shares the penicillin-binding protein (PBP) core structural folding but found to hydrolyze the removal of d-Ala from teichoic acids. Molecular docking, dynamics, and MM-GBSA of FmtA and its variants S127A, K130A, Y211A, D213A, and K130AY211A, in the presence or absence of wall teichoic acid (WTA), suggest that active site residues S127, K130, Y211, D213, N343, and G344 play a role in substrate binding. Quantum mechanics (QM)/molecular mechanics (MM) calculations reveal that during WTA catalysis, K130 deprotonates S127, and the nucleophilic S127 attacks the carbonyl carbon of d-Ala bound to WTA. The tetrahedral intermediate (TI) complex is stabilized by hydrogen bonding to the oxyanion holes. The TI complex displays a high energy gap and collapses to an energetically favorable acyl-enzyme complex.

Identification of an AXL kinase inhibitor in triple-negative breast cancer by structure-based virtual screening and bioactivity test

Breast cancer is a malignant tumor that occurs in the glandular epithelium of the breast, and more than 15% of the patients are triple-negative breast cancer (TNBC). Therefore, finding new targets and targeted therapeutic drugs for TNBC is urgent. Overexpression of the AXL is associated with motility and invasiveness of the TNBC cells, which is a potential target for breast cancer therapy. A compound Y041-5921 (IC₅₀  = 6.069 μm for AXL kinase and IC₅₀  = 4.1 μm for MDA-MB-231 cell line) was identified through structure-based virtual screening and bioassay test for the first time. The compound Y041-5921 could significantly inhibit the proliferation and invasion of the TNBC cells and the toxicity of Y041-5921 to normal immortalized breast epithelial cells was far lower than that of commonly used clinical chemotherapy drugs. Besides, it also had well inhibitory effect on the proliferation of many other malignant tumor cell lines (the IC₅₀  value are 10.0 m, 7.1 m, 10.3 m, 11.4 m and 5.8 m for U251 cell, COLO cell, PC-9 cell, CAKI-1 cell and MG63 cell, respectively). The interaction mechanism between Y041-5921 and AXL was studied by molecular dynamics (MD) simulations and binding free energy calculation, and the key residues whose energy contribution mainly comes from non-polar solvation interaction (such as Ala565, Lys567, Met598, Leu620, Pro621, Met623, Lys624, Arg676, Asn677 and Met679) were identified. The small molecule inhibitors Y041-5921 targeting AXL reported in this work will lay a foundation and provide a theoretical basis for the development of the TNBC.

Binding of hydroxychloroquine and chloroquine dimers to palmitoyl-protein thioesterase 1 (PPT1) and its glycosylated forms: a computational approach

The lysosomal enzyme palmitoyl-protein thioesterase 1 (PPT1) removes thioester-linked fatty acid groups from membrane-bound proteins to facilitate their proteolysis. A lack of PPT1 (due to gene mutations) causes the progressive death of cortical neurons and is responsible for infantile neural ceroid lipofuscinosis (INCL), a severe neurodegenerative disorder in children. Conversely, PPT1 is often over-expressed in cancer, and considered as a valid target to control tumor growth. Potent and selective inhibitors of PPT1 have been designed, in particular 4-amino-7-chloro-quinoline derivatives such as hydroxychloroquine (HCQ) and the dimeric analogues Lys05 and DC661. We have modeled the interaction of these three compounds with the enzyme, taking advantage of the PPT1 crystallographic structure. The molecules can fit into the palmitate site of the protein, with the dimeric compounds forming more stable complexes than the monomer. But the molecular modeling suggests that the most favorable binding sites are located outside the active site. Two sites centered on residues Met112 and Gln144 were identified, offering suitable cavities for drug binding. According to the calculated empirical energies of interaction (ΔE), the dimer DC661 forms the most stable complex at site Met112 of palmitate-bound PPT1. N-glycosylated forms of PPT1 were elaborated. Paucimannosidic glycans (M2FA and M3F) and a bulkier tetra-antennary complex glycan were introduced at asparagine residues N197, N212 and N232. These N-glycans do not impede drug binding, thus suggesting that all glycoforms of PPT1 can be targeted with these compounds.Communicated by Ramaswamy H. Sarma.