Structural insights into Entamoeba histolytica arginase and structure‐based identification of novel non‐amino acid based inhibitors as potential antiamoebic molecules

New thiazolidine-2,4-diones as potential anticancer agents and apoptotic inducers targeting VEGFR-2 kinase: Design, synthesis, in silico and in vitro studies

BACKGROUND

VEGFR-2 has emerged as a prominent positive regulator of cancer progression.

AIM

Discovery of new anticancer agents and apoptotic inducers targeting VEGFR-2.

METHODS

Design and synthesis of new thiazolidine-2,4-diones followed by extensive in vitro studies, including VEGFR-2 inhibition assay, MTT assay, apoptosis analysis, and cell migration assay. In silico investigations including docking, MD simulations, ADMET, toxicity, and DFT studies were performed.

RESULTS

Compound 15 showed the strongest VEGFR-2 inhibitory activity with an IC50 value of 0.066 μM. Additionally, most of the synthesized compounds showed anti-proliferative activity against HepG2 and MCF-7 cancer cell lines at the micromolar range with IC50 values ranging from 0.04 to 4.71 μM, relative to sorafenib (IC50 = 2.24 ± 0.06 and 3.17 ± 0.01 μM against HepG2 and MCF-7, respectively). Also, compound 15 showed selectivity indices of 1.36 and 2.08 against HepG2 and MCF-7, respectively. Furthermore, compound 15 showed a significant apoptotic effect and arrested the cell cycle of MCF-7 cells at the S phase. Moreover, compound 15 had a significant inhibitory effect on the ability of MCF-7 cells to heal from. Docking studies revealed that the synthesized thiazolidine-2,4-diones have a binding pattern approaching sorafenib. MD simulations indicated the stability of compound 15 in the active pocket of VEGFR-2 for 200 ns. ADMET and toxicity studies indicated an acceptable pharmacokinetic profile. DFT studies confirmed the ability of compound 15 to interact with VEGFR-2.

CONCLUSION

Compound 15 has promising anticancer activity targeting VEGFR-2 with significant activity as an apoptosis inducer.

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.

Molecular docking, dynamics simulation and pharmacokinetic studies of Cyperus articulatus essential oil metabolites as inhibitors of Staphylococcus aureus

Cyperus articulatus has been extensively studied for its essential oil (EO), active components and antibacterial activities against a wide range of bacteria such as Bacillus megaterium, Streptococcus pyogenes, Staphylococcus epidermidis, Escherichia coli and Staphylococcus aureus. However, knowledge of the biomolecular interaction of the individual EO metabolites responsible for its inhibition activities is lacking. The multi-drug-resistant bacteria S. aureus, which is of prime concern, has been reported to be inhibited by Cyperus articulatus rhizome EO. The present work analyzed the molecular interactions of the major Cyperus articulatus rhizome EO metabolites with the target enzyme TyrRS of S. aureus and studied the conformational dynamics and stability of the protein-ligand complexes. Molecular docking studies of selected EO metabolites such as mustakone, longifolenaldehyde, cyperotundone, α-copaene, β-calacorene, α-calacorene and khusinol were conducted along with standard drug chloramphenicol for comparative analysis of their binding affinity with S. aureus TyrRS. The metabolites khusinol, mustakone, β-calacorene and α-calacorene generated comparable docking scores (-6.4, -6.2, -6.1 and -6.2 kcal/mol, respectively) with that of the drug chloramphenicol (-6.3 kcal/mol). Most EO metabolites did not exhibit H-bonding with the S. aureus TyrRS residues and were stabilized through pi-interactions. The MD simulation study illustrated that compounds like mustakone could effectively bind to the receptors of S. aureus TyrRS with high stability and integrity. Pharmacokinetic, drug-like properties and toxicity analysis of the EO metabolites supported the candidature of mustakone and khusinol as pharmacologically important antibacterial drug ingredients. The study envisaged the structural framework of the EO metabolites for antibacterial drug design.Communicated by Ramaswamy H. Sarma.

Compound effect and mechanism of oxidative damage induced by nanoplastics and benzo [a] pyrene

Nanoplastics and polycyclic aromatic hydrocarbons (PAHs) are ubiquitous in soil environments. In order to objectively evaluate the toxic interaction between polystyrene nanoplastics (PS NPs) and benzo [a] pyrene (BaP), oxidative damage at the level of earthworm cells and biomacromolecules was investigated by experiments combined with molecular dynamics simulation. Studies on cells reveal that PS NPs and BaP had synergistic toxicity when it came to causing oxidative stress. Cellular reactive oxygen species (ROS) levels under combined pollutant exposure were 24% and 19% higher, respectively than when PS NPs and BaP were exposed alone (compared to the blank group). In addition, BaP and PS NPs inhibited the ability of CAT to decompose H2O2 by affecting the structure of the proximal amino acid Tyr 357 in the active center of CAT, which exacerbated oxidative stress to a certain extent. Therefore, the synergistic toxic effect of BaP and PS NPs is due to the mutual complement of the two to the induction of protein structural looseness, and the strengthening of the stability of the conjugate (CAT-BaP-PS) under the weak interaction. This work provides a new perspective and approach on how to talk about the toxicity of combined pollutants.

In vitro and in silico evaluation of new thieno[2,3-d]pyrimidines as anti-cancer agents and apoptosis inducers targeting VEGFR-2

In this study, new thieno[2,3-d]pyrimidine derivatives that could have potential anticancer activity by inhibiting the VEGFR-2 receptor have been designed, synthesized, and investigated. The thieno[2,3-d]pyrimidine derivatives showed strong in vitro abilities to inhibit VEGFR-2 and to prevent cancer cell growth in two different types of cancer cells, MCF-7 and HepG2. Particularly, compound 22 showed the most potent anti-VEGFR-2 activity with an IC50 value of 0.58 µM. Additionally, compound 22 exhibited good anti-proliferative activity against both MCF-7 and HepG2 cancer cell lines, with IC50 values of 11.32 ± 0.32 and 16.66 ± 1.22 µM, respectively. Further investigations revealed that compound 22 induced cell cycle arrest at the G2/M phase and promoted both early and late apoptosis in the MCF-7 cancer cells. Compound 22 also increased the level of BAX (2.8-fold), and reduced the level of Bcl-2 (2.2-fold), hence increasing the rate of apoptosis. Compound 22 also revealed 2.9-fold and 2.8-fold higher levels of caspase-8 and caspase-9, respectively, in the treated MCF-7 cancer cells compared to the control cell lines. The MD simulations showed that the VEGFR-2-22 complex was structurally and energytically stable over 100 ns, while the MM-GBSA study indicated its stable thermodynamic behavior. The bi-dimensional projection analysis confirmed the proper binding of the VEGFR-2-22 complex, while the DFT studies provided optimized geometry, charge distribution, FMO, ESP, the total density of state, and QTAIM maps of compound 22. Finally, computational ADMET studies were performed to assess the drug development potential of the thieno[2,3-d]pyrimidine derivatives. Overall, this study suggests that compound 22 has the potential as an anticancer lead compound by inhibiting VEGFR-2, which may be a guide for future drug design and development.

Illuminating the structure-function landscape of an evolutionary nonconserved motif in the arginases of Helicobacter gastric pathogens

The bimetallic enzyme arginase catalyses the conversion of L-arginine to L-ornithine and urea. In Helicobacter pylori (a known human gastric pathogen), this enzyme is an important virulence factor. In spite of the conservation of the catalytic and the metal-binding residues, the H. pylori homolog possesses a 13-residue motif (-153 ESEEKAWQKLCSL165 -) present in the middle of the protein sequence, whose role was recently elucidated. Despite several reviews available on arginases, no report has thoroughly illustrated the underlying basis for the importance of the above motif of the H. pylori enzyme in structure and function. In this review, we systematically describe a mechanistic basis for its importance in structure and function based on the known data. This motif of the H. pylori enzyme is present exclusively in the arginases of other Helicobacter gastric pathogens, where the critical residues are conserved, implying that the nonconserved stretch has been selected during the evolution of the enzyme in these gastric pathogens in a specific manner to perform its role in the structure and function. The combined information can be useful for understanding the function of arginases in other Helicobacter gastric pathogens. Additionally, this knowledge can be utilised to screen and design new small molecule inhibitors, specific to the arginases of these pathogens.

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.

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.

Computational characteristics of the structure-activity relationship of inhibitors targeting Pks13-TE domain

Multiple studies have established the Pks13-TE domain as a promising target for anti-tuberculosis drug development. However, recent findings have revealed that the lead compound currently in the pipeline for Pks13-TE has significant cardiotoxicity issues. Given the pressing need for new chemical structures for Pks13-TE inhibitors, this study aims to provide a detailed understanding of the Pks13-TE domain binding site through the application of computational chemical biology techniques. Our results highlight the size and shape of the Pks13-TE domain binding pocket, key residues including Asp1644, Asn1640, Phe1670, and Tyr1674 within the pocket, and inhibitor pharmacophore characteristics such as aromatic ring sites, positively charged sites, and hydrogen bond donors. To our knowledge, these simulation results are novel and contribute to the discovery of next-generation Pks13-TE inhibitors without similar prior studies.

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 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.

Sequence Divergence in the Arginase Domain of Ornithine Decarboxylase/Arginase in Fusobacteriacea Leads to Loss of Function in Oral Associated Species

A number of species within the Fusobacteriaceae family of Gram-negative bacteria uniquely encode for an ornithine decarboxylase/arginase (ODA) that ostensibly channels l-ornithine generated by hydrolysis of l-arginine to putrescine formation. However, two aspartate residues required for coordination to a catalytically obligatory manganese cluster of arginases are substituted for a serine and an asparagine. Curiously, these natural substitutions occur only in a clade of Fusobacterium species that inhabit the oral cavity. Herein, we expressed and isolated full-length ODA from the opportunistic oral pathogen Fusobacterium nucleatum along with the individual arginase and ornithine decarboxylase components. The crystal structure of the arginase domain reveals that it adopts the classical α/β arginase-fold, but metal ions are absent in the active site. As expected, the ureohydrolase activity with l-arginine was not detected for wild-type ODA or the isolated arginase domain. However, engineering of the complete metal coordination environment through site-directed mutagenesis restored Mn²⁺ binding capacity and arginase activity, although the catalytic efficiency for l-arginine was low (60-100 M^-1 s^-1). Full-length ODA and the isolated ODC component were able to decarboxylate both l-ornithine and l-arginine to form putrescine and agmatine, respectively, but k_cat/K_M of l-ornithine was ∼20-fold higher compared to l-arginine. We discuss environmental conditions that may have led to the natural selection of an inactive arginase in the oral associated species of Fusobacterium.

In-silico evidences on filarial cystatin as a putative ligand of human TLR4

Cystatin is a small molecular weight immunomodulatory protein of filarial parasite that plays a pivotal role in downregulating the host immune response to prolong the survival of the parasite inside the host body. Hitherto, this protein is familiar as an inhibitor of human proteases. However, growing evidences on the role of cystatin in regulating inflammatory homeostasis prompted us to investigate the molecular reasons behind the explicit anti-inflammatory trait of this protein. We have explored molecular docking and molecular dynamics simulation approaches to explore the interaction of cystatin of Wuchereria bancrofti (causative parasite of human filariasis) with human Toll-like receptors (TLRs). TLRs are the most crucial component of frontline host defence against pathogenic infections including filarial infection. Our in-silico data clearly revealed that cystatin strongly interacts with the extracellular domain of TLR4 (binding energy=-93.5 ± 10 kJ/mol) and this biophysical interaction is mediated by hydrogen bonding and hydrophobic interaction. Molecular dynamics simulation analysis revealed excellent stability of the cystatin-TLR4 complex. Taken together, our data indicated that cystatin appears to be a ligand of TLR4 and we hypothesize that cystatin-TLR4 interaction most likely to play a key role in activating the alternative activation pathways to establish an anti-inflammatory milieu. Thus, the study provokes the development of chemotherapeutics and/or vaccines for targeting the cystatin-TLR4 interaction to disrupt the pathological attributes of human lymphatic filariasis. Our findings are expected to provide a novel dimension to the existing knowledge on filarial immunopathogenesis and it will encourage the scientific communities for experimental validation of the present investigation. Communicated by Ramaswamy H. Sarma.

In-silico screening and identification of potential inhibitors against 2Cys peroxiredoxin of Candidatus Liberibacter asiaticus

Huanglongbing (HLB) is a worldwide citrus plant disease-related to non-culturable and fastidious α-proteobacteria Candidatus Liberibacter asiaticus (CLas). In CLas, Peroxiredoxin (Prx) plays a major role in the reduction of the level of reactive species such as reactive oxygen species (ROS), free radicals and peroxides, etc. Here, we have used structure-based drug designing approach was used to screen and identify the potent molecules against 2Cys Prx. The virtual screening of fragments library was performed against the three-dimensional validated model of Prx. To evaluate the binding affinity, the top four molecules (N-Boc-2-amino isobutyric acid (B2AI), BOC-L-Valine (BLV), 1-(boc-amino) cyclobutane carboxylic acid (1BAC), and N-Benzoyl-DL-alanine (BDLA)) were docked at the active site of Prx. The molecular docking results revealed that all the identified molecules had a higher binding affinity than Tert butyl hydroperoxide (TBHP), a substrate of Prx. Molecular dynamics analysis such as RMSD, Rg, SASA, hydrogen bonds, and PCA results indicated that Prx-inhibitor(s) complexes had lesser fluctuations and were more stable and compact than Prx-TBHP complex. MMPBSA results confirmed that the identified compounds could bind at the active site of Prx to form a lower energy Prx-inhibitor(s) complex than Prx-TBHP complex. The identified potent molecules may pave the path for the development of antimicrobial agents against CLA.Communicated by Ramaswamy H. Sarma.

Designing a natural inhibitor against human kynurenine aminotransferase type II and a comparison with PF-04859989: a computational effort against schizophrenia

The Kynurenine aminotransferase II (KATII) enzyme has an essential role in L-kynurenine transmission to kynurenic acid (KYNA). High concentration of kynurenic acid associates with schizophrenia and some neurocognitive disorders. Decreasing KYNA production via inhibiting KATII would be an effective method for treating and understanding the related central nervous system (CNS) diseases. This study aimed to discover a potent inhibitor against human KATII (hKATII) in comparison with PF-04859989. We utilized the computational methods of molecular dynamics, virtual screening, docking, and binding free-energy calculations. Initially, the 58722 compounds from three drug libraries, including IBS library, DrugBank library, and Analyticon library, were obtained. At the next stage, these sets of compounds were screened by AutoDock Vina software, and a potent inhibitor (ZINC35466084) was selected. Following the screening, molecular dynamics simulations for both ZINC35466084 and PF-04859989 were performed by GROMACS software. MM-PBSA analysis showed that the amount of binding free energy for ZINC35466084 (-61.26 KJ mol^-1) is more potent than PF-04859989 (-43.14 KJ mol^-1). Furthermore, the ADME analysis results revealed that the pharmacokinetic parameters of ZINC35466084 are acceptable for human use. Eventually, our data demonstrated that ZINC35466084 is suitable for hKATII inhibition, and it is an appropriate candidate for further studies in the laboratory. Communicated by Ramaswamy H. Sarma.

Structure-based mimicking of hydroxylated biphenyl congeners (OHPCBs) for human transthyretin, an important enzyme of thyroid hormone system

In humans, transthyretin (hTTR) is a plasma protein act as a transporter of thyroxine (T4) in the blood. Polychlorinated biphenyls (PCBs) are used in coolants, transformers, plasticizers, and pesticide extenders, etc. due to their physical properties, chemical stability, and dielectric properties. Cytochrome P450 can oxidize the PCBs into hydroxylated PCBs (OHPCBs) which can further interact with hTTR results in hepatoxicity, loss of metabolic rate, memory problems, and neurotoxicity. Molecular docking results show that OHPCBs bind at the active site of hTTR with a more binding affinity as compared to T4. Further, molecular dynamics simulation has been done to confirm the stability of hTTR-OHPCBs complexes. Several analysis parameters like RMSD, RMSF, Rg, SASA, hydrogen bonds numbers, PCA, and FEL revealed that binding of OHPCBs with hTTR results in the formation of stable hTTR-OHPCBs complexes. Individual residues decomposition analysis confirms that Lys15, Leu17, Ala108, Ala109, Leu110, Ser117, and Thr119 of hTTR plays a major role in the binding of OHPCBs to form the lower energy hTTR-OHPCBs complexes. Molecular docking and simulations results emphasize that OHPCBs can efficiently bind at the active site of hTTR, which further leads to inhibition of transportation of T4 in human blood.

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

Staphylococcus aureus is resistant to β-lactam antibiotics and causes several skin diseases to life-threatening diseases. FmtA is found to be one of the main factors involved in methicillin resistance in S. aureus. FmtA exhibits an esterase activity that removes the D-Ala from teichoic acid. Teichoic acids played a significant role in cell wall synthesis, cell division, colonization, biofilm formation, virulence, antibiotic resistance, and pathogenesis. The virtual screening of drug molecules against the crystal structure of FmtA was performed and the binding affinities of top three molecules (ofloxacin, roflumilast, and furazolidone) were predicted using molecular docking. The presence of positive potential and electron affinity regions in screened drug molecules by DFT analysis illustrated that these molecules are reactive in nature. The protein-ligand complexes were subjected to molecular dynamics simulation. Molecular dynamics analysis such as RMSD, RMSF, Rg, SASA, PCA, and FEL results suggested that FmtA-drug(s) complexes are stable. MM-GBSA binding affinity and QM/MM results (ΔG, ΔH, and ΔS) revealed that active site residues (Ser127, Lys130, Tyr211, Asp213, and Asn343) of FmtA played an essential for the binding of the drug(s) to form a lower energy stable protein-ligand complexes. FmtA^Δ42 was purified using cation exchange and gel filtration chromatography. Fluorescence spectroscopy and circular dichroism results showed that interactions of drugs with FmtA^Δ42 affect the tertiary structure and increase the thermostability of the protein. The screened molecules need to be tested and could be further modified to develop the antimicrobial compounds against S. aureus.

Computational guided identification of novel potent inhibitors of N-terminal domain of nucleocapsid protein of severe acute respiratory syndrome coronavirus 2

The Coronavirus Disease 2019, caused by the severe acute respiratory syndrome coronavirus 2 is an exceptionally contagious disease that leads to global epidemics with elevated mortality and morbidity. There are currently no efficacious drugs targeting coronavirus disease 2019, therefore, it is an urgent requirement for the development of drugs to control this emerging disease. Owing to the importance of nucleocapsid protein, the present study focuses on targeting the N-terminal domain of nucleocapsid protein from severe acute respiratory syndrome coronavirus 2 to identify the potential compounds by computational approaches such as pharmacophore modeling, virtual screening, docking and molecular dynamics. We found three molecules (ZINC000257324845, ZINC000005169973 and ZINC000009913056), which adopted a similar conformation as guanosine monophosphate (GMP) within the N-terminal domain active site and exhibiting high binding affinity (>−8.0 kcalmol⁻¹). All the identified compounds were stabilized by hydrogen bonding with Arg107, Tyr111 and Arg149 of N-terminal domain. Additionally, the aromatic ring of lead molecules formed π interactions with Tyr109 of N-terminal domain. Molecular dynamics and Molecular mechanic/Poisson–Boltzmann surface area results revealed that N-terminal domain – ligand(s) complexes are less dynamic and more stable than N-terminal domain – GMP complex. As the identified compounds share the same corresponding pharmacophore properties, therefore, the present results may serve as a potential lead for the development of inhibitors against severe acute respiratory syndrome coronavirus 2.

Communicated by Ramaswamy H. Sarma

Structural characterization and in-silico analysis of Momordica charantia 7S globulin for stability and ACE inhibition

Momordica charantia (Mc) seeds are widely used edible crop with high nutritional quality. The food and pharmaceutical industries use it as a natural anti-oxygenic agent. Herein, a ~52 kDa protein, which is a major part of seed proteome has been purified, biochemically characterized and structure has been determined. MALDI-ESI-MS identified peptide fragments and contig-deduced sequence suggested the protein to be homologous to 7S globulins. The crystal structure shows that protein has a bicupin fold similar to 7S globulins and the electron density for a copper and acetate ligand were observed in the C-terminal barrel domain. In silico study reveals that a tripeptide (VFK) from Mc7S possess a higher binding affinity for angiotensin converting enzyme (ACE) than already reported drug Lisinopril (LPR). The protein is a glycoprotein and highly stable under varying thermal and pH conditions due to its secondary structures. The DPPH (2,2-diphenyl-1-picryl-hydrazyl-hydrate) assay showed the protein to have an anti-oxygenic nature and can aid in scavenging free radical from sample. The protein can assist to enhance the nutritional and functional value of food by acting as a food antioxidant. Further, characterization of Mc7S required which might add in importance of Mc7S as antioxidant, anti-diabetic and anti-hypertensive.

Molecular docking and dynamic approach to virtual screen inhibitors against Esbp of Candidatus Liberibacter asiaticus

Citrus greening (huanglongbing) is the most destructive disease of citrus worldwide caused by Candidatus Liberibacter asiaticus (CLA). Currently, no strategies have been developed to manage the Huanglongbing (HLB) disease and to stop the spreading of this disease to new citrus areas. Esbp is an extracellular solute-binding protein, involved in the uptake of iron in CLA. Thus, inhibiting this process may be a promising approach to design a drug against CLA. Thus, the present study focused on the identification of novel effective inhibitors which can inhibit the activity of CLas Esbp. A series of small molecules were screened against the CLas Esbp and the binding affinities were assessed using docking simulation studies. Top scored molecules were screened for different pharmacophore properties and Inhibitory Concentration 50 (IC50) values. Density functional theory was employed to check the chemical properties of the molecules. Further, Molecular Dynamics simulation analysis like RMSD, RMSF, Rg, SASA and MMPBSA results reveal that the identified molecules (ZINC03143779, ZINC05491830, ZINC19210425, ZINC08750867, and ZINC14671545) exhibit a good binding affinity for CLas Esbp and results in the formation of stable CLas Esbp-inhibitor(s) complex. The present study reported that these compounds appeared to be the suitable novel inhibitor of CLas Esbp and pave the way to further development of antimicrobial agents against CLA.