Exploring the selectivity of auto-inducer complex with LuxR using molecular docking, mutational studies and molecular dynamics simulations

Anti-quorum sensing effects of SidA protein on Escherichia coli receptors: in silico analysis

Quorum sensing enables cell-cell communication in bacteria and regulates biofilm formation. Biofilm production promotes pathogenicity of Escherichia coli and causes infections. However, antibiotic resistance limits conventional treatment efficacy against biofilm infections. Quorum quenching offers an alternative by disrupting quorum sensing signals. Allicin, extracted from garlic, possesses antimicrobial and anti-quorum sensing properties. This study employed molecular docking and dynamics simulations to investigate allicin's interaction with the E. coli quorum sensing system, specifically targeting the cytoplasmic SidA receptor protein. SidA binds to N-acyl-homoserine lactone ligands and regulates quorum sensing in E. coli. The crystal structure of SidA was obtained from the PDB. Molecular docking revealed that allicin competitively binds to the ligand-binding pocket of SidA. Simulations analyzed the effects of allicin binding on SidA stability and affinity for N-acyl-homoserine lactones over 200 ns. Parameters like RMSD, RMSF, and hydrogen bonding indicated SidA was more stable when bound to allicin compared to unbound. Binding free energies suggested allicin reduced SidA's affinity for native ligands. Therefore, allicin binding to SidA alters its conformation and inhibits interaction with N-acyl-homoserine lactones, disrupting quorum sensing signaling and biofilm production in E. coli.Communicated by Ramaswamy H. Sarma.

Computational identification and experimental validation of anti-filarial lead molecules targeting metal binding/substrate channel residues of Cu/Zn SOD1 from Wuchereria bancrofti

Lymphatic filariasis (LF) is a neglected mosquito-borne parasitic disease, widely caused by Wuchereria bancrofti (Wb) in tropical and sub-tropical countries. During a blood meal, the filarial nematodes are transmitted to humans by the infected mosquito. To counter attack the invaded nematodes, the human immune system produces reactive oxygen species. However, the anti-oxidant enzymes of nematodes counteract the host oxidative cytotoxicity. Cu/Zn Superoxide dismutase (SOD1), a member of antioxidant enzymes and are widely used by the nematodes to sustain the host oxidative stress across its lifecycle, hence targeting SOD1 to develop suitable drug molecules would help to overcome the problems related to efficacy and activity of drugs upon different stages of nematodes. In order to find the potent inhibitors, a three-dimensional structure of Cu/Zn WbSOD1 was modelled and the structural stability was analysed through simulation studies. The structure-guided virtual screening approach has been used to identify lead molecules from the ChemBridge based on the docking score, ADMET properties and protein-ligand complex stability analysis. The identified compounds were observed to interact with the copper, metal binding residues (His48, His63, His80 and His120) and catalytically important residue Arg146, which play a crucial role in the disproportionation of incoming superoxide radicals of Cu/Zn WbSOD1. Further, in vitro validation of the selected leads in the filarial worm Setaria digitata exhibited higher inhibition and better IC50 compared to the standard drug ivermectin. Thus, the identified leads could potentially inhibit enzyme activity, which could subsequently act as drug candidates to control LF.Communicated by Ramaswamy H. Sarma.

Potential inhibitors for peroxiredoxin 6 of W. bancrofti: A combined study of modelling, structure-based drug design and MD simulation

Lymphatic filariasis (LF), a mosquito-borne parasitic disease caused by nematode Wuchereria bancrofti in tropical and sub-tropical countries. These nematodes are transferred into the human host when the infected mosquito carrying L3 larvae is released into the bloodstream during the blood ingestion process. The host immune system produces ROS (Reactive Oxygen Species) as a primary defence mechanism to remove the invading filarial worms. However, well-defined antioxidant enzymes of the nematodes scavenge the host-produced ROS to escape from oxidative stress. The enzyme peroxiredoxin 6 (Prx6) belongs to the peroxiredoxin family, catalyses hydrogen peroxide (H₂O₂) into water (H₂O). In order to find the inhibitors that inhibit the activity of peroxiredoxin 6 of W. bancrofti. We performed the homology modelling to predict the WbPrx6 three-dimensional structure using the Schrödinger-Prime and the dynamic stability of the modelled WbPrx6 was analyzed by carrying out the molecular dynamic (MD) simulation for the time scale of 200ns. Further, the structure-based virtual screening shortlisted the hit molecules from the ChemBridge database based on the glide score. The potential lead molecules (ID: 10239274, 11112883, 79879205, 58160895, and 42133744) that have better binding and satisfied the ADMET properties were selected for further complex simulation and DFT calculations. The identified compounds interact with the N-terminal region of the thioredoxin domain, which plays a key role in reducing phospholipase A2 activity. Interestingly, upon binding the lead molecule, the fluctuation of the loop region that connects α-IV with the β-VI plays a vital role in affecting the geometry of the active site, which in turn affects the activity WbPrx6. The outcomes of the present computational studies could help in future drug development and designing of the effective candidate to control Lymphatic filariasis.

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

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

Investigating the interaction between three perfluorinated carboxylic acids and the G protein-coupled estrogen receptor: spectroscopic analyses and computational simulations

In this paper, perfluorinated compounds (PFCs), such as perfluorobutyric acid (PFBA), perfluorooctanoic acid (PFOA) and perfluorododecanoic acid (PFDoA), were selected as typical representatives of perfluorinated carboxylic acids (PFCAs) to study the effects of PFCAs on the G protein-coupled estrogen receptor (GPER). The interaction mechanism of the three types of PFCAs with the GPER was investigated using steady-state fluorescence spectroscopy, ultraviolet-visible spectroscopy, three-dimensional fluorescence spectroscopy, and Fourier transform infrared spectroscopy combined with molecular docking and molecular dynamics simulations. Among these techniques, steady-state fluorescence and ultraviolet-visible spectroscopic analyses showed that PFBA, PFOA and PFDoA quenched the endogenous GPER fluorescence by combined dynamic and static quenching and non-radiative energy transfer. The binding constants (K_a) of PFCAs on the GPER were all larger than 10⁵ L mol^-1, indicating that their affinity for the GPER was strong. Fourier transform infrared spectroscopy and three-dimensional fluorescence showed that the secondary structure of the GPER changed after binding to PFCAs. Thermodynamic analysis showed ΔG < 0, which indicated that the interaction between the GPER and PFCAs was spontaneous. For the binding of PFBA and PFOA to the GPER, ΔH > 0 and ΔS > 0, indicating that the interaction was mainly driven by hydrophobic forces; for the binding of PFDoA to the GPER, ΔH < 0 and ΔS < 0, suggesting that van der Waals force and hydrogen bonding were the main interaction forces. Molecular dynamics simulations suggested that the stability of the GPER-PFCA complexes was higher than that of the free GPER, and also that the structure and hydrophobicity of the GPER changed after binding to PFCAs. Molecular docking analysis showed that all three PFCAs could form hydrogen bonds with the GPER, which improved the stability of the complex.

Evolutionary significance and functional characterization of streptomycin adenylyltransferase from Serratia marcescens

Complete functional annotations of proteins are essential to understand the role and mechanisms in pathogenesis. Aminoglycoside nucleotidyltransferases are the subclasses of aminoglycosides modifying enzymes conferring resistance to organisms. Insight into the structural and functional understanding of nucleotidyltransferase family protein provides vital information to combat pathogenesis. Phylogenetic analysis is employed to identify the evolutionary significance and common motif's present among the homologs of nucleotidyltransferase family protein. Structure, sequence based approaches and molecular docking were implemented to predict the exact function of the protein. Wide distribution of the nucleotidyltransferase family protein in gram-positive and gram-negative organisms are evidenced from phylogenetic analysis. Five common motifs were present in all the homolog's of nucleotidyltransferase family protein. Sequence-structure based functional annotations predicts that the targeted protein function as ATP-Mg dependent streptomycin adenylyltransferase. Structural comparisons and docking studies correlate well with the identified function. The complete function of nucleotidyltransferase family protein was identified as Streptomycin adenylyltransferase and it could be targeted as a potential therapeutic target to overcome antibiotic resistance.Communicated by Ramaswamy H. SarmaAbbreviationsAACaminoglycoside acetyltransferasesAMEaminoglycoside modifying enzymeANTaminoglycoside nucleotidyltransferasesAPHaminoglycoside phosphotransferasesATPadenosine triphosphateCASTpcomputer atlas and surface topography of proteinsDUFdomains of unknown functionGlidegrid-based ligand docking with energeticHMMhidden Markov modelMASTmotif alignment and search toolMEGAmolecular evolutionary genetics analysisMEMEmultiple Em for motif elicitationMSAmultiple sequence alignmentNMPnucleoside monophosphateNTPnucleoside triphosphateNTnucleotidyltransferaseOPLSoptimized potential for liquid simulationXPextra precision.

Phthalic acid esters disturbed the genetic information processing and improved the carbon metabolism in black soils

Phthalic acid esters (PAEs), such as dimethyl phthalate (DMP) and dibutyl phthalate (DBP), are widely distributed as environmental pollutants. In this study, the effects of these chemicals were investigated in black soils using a metagenomics approach. The results clearly showed that DMP or DBP increased the abundance of genes involved in transcription, replication and repair in black soils. In addition, the abundances of genes associated with metabolic functions was improved following treatment with DMP or DBP, including those involved in lipid transport and metabolism, carbohydrate transport and metabolism, and energy production and conversion. There could be many reasons for these observed changes. First, the DMP or DBP treatments increased the abundances of genes associated with the LuxR family, the UvrABC repair system, DNA replication pathways, the RNA polymerase complex and base excision repair. Second, the abundances of genes associated with isocitrate lyase regulator (IclR) family transcriptional regulators, lipid metabolism and carbohydrate active enzymes (CAZys) were altered by the DMP or DBP treatments. Finally, the DMP or DBP treatments also increased the emission load of CO₂ and altered the fluorescence intensity of humic acid. Therefore, the results of this study suggested that DMP and DBP contamination altered the abundances of genes associated with genetic information processing and improved the carbon metabolism in black soils.

The role of LasR active site amino acids in the interaction with the Acyl Homoserine Lactones (AHLs) analogues: A computational study

The present work combines molecular docking calculations, 3D-QSAR, molecular dynamics simulations and free binding energy calculations (MM/PBSA and MM/GBSA) in a set of 28 structural analogues of acyl homoserine lactones with Quorum Sensing antagonist activity. The aim of this work is to understand how ligand binds and is affected by the molecular microenvironment in the active site of the LasR receptor for pseudomonas aeruginosa. We also study the stability of the interaction to find key structural characteristics that explain the antagonist activities of this set of ligands. This information is relevant for the rational modification or design of molecules and their identification as powerful LasR modulators. The analysis of molecular docking simulations shows that the 28 analogues have a similar binding mode compared to the native ligand. The carbonyl groups belonging to the lactone ring and the amide group of the acyl chain are oriented towards the amino acids forming hydrogen bond like interactions. The difference in antagonist activity is due to location and orientation of the LasR side chains within the hydrophobic pocket in its binding site. Additionally, we carried out molecular dynamics simulations to understand the conformational changes in the ligand-receptor interaction and the stability of each complex. Results show a direct relationship among the interaction energies of the ligands and the activities as an antagonist of the LasR receptor.