Performance of HADDOCK and a simple contact-based protein-ligand binding affinity predictor in the D3R Grand Challenge 2

Chitosan-Coated Silver Nanoparticles Inhibit Adherence and Biofilm Formation of Uropathogenic Escherichia coli

Urinary tract infections are commonly caused by uropathogenic Escherichia coli (UPEC), which usually presents multiple virulence and resistance mechanisms, making it difficult to treat. It has been demonstrated that silver and polymeric nanoparticles had potential against these pathogens. In this study, we synthesized thiol chitosan-coated silver nanoparticles (SH-Cs-AgNPs) and evaluated their antibacterial, antibiofilm and antiadherence activity against clinical isolates of UPEC. The SH-Cs-AgNPs showed a spherical shape with a size of 17.80 ± 2.67 nm and zeta potential of 18 ± 2 mV. We observed a potent antibacterial and antibiofilm activity as low as 12.5 μg/mL, as well as a reduction in the adherence of UPEC to mammalian cells at concentrations of 1.06 and 0.53 μg/mL. These findings demonstrate that SH-Cs-AgNPs have potential as a new therapeutic compound against infections caused by UPEC.

'Computational studies on coumestrol-ArlR interaction to target ArlRS signaling cascade involved in MRSA virulence'

Two component signaling system ArlRS (Autolysis-related locus) regulates adhesion, biofilm formation and virulence in methicillin resistant Staphylococcus aureus. It consists of a histidine kinase ArlS and response regulator ArlR. ArlR is composed of a N-terminal receiver domain and DNA-binding effector domain at C-terminal. ArlR receiver domain dimerizes upon signal recognition and activates DNA binding by effector domain and subsequent virulence expression. In silico simulation and structural data suggest that coumestrol, a phytochemical found in Pueraria montana, forges a strong intermolecular interaction with residues involved in dimer formation and destabilizes ArlR dimerization, an essential conformational switch required for downstream effector domain to bind to virulent loci. Structural and energy profiles of simulated ArlR-coumestrol complexes suggest lower affinity between ArlR monomers due to structural rigidity at the dimer interface hindering the conformational rearrangements relevant for dimer formation. These analyses could be an attractive strategy to develop therapeutics and potent leads molecules response regulators of two component systems in which are involved in MRSA virulence as well as other drug-resistant pathogens.Communicated by Ramaswamy H. Sarma.

Butterfly Effect in Cytarabine: Combined NMR-NQR Experiment, Solid-State Computational Modeling, Quantitative Structure-Property Relationships and Molecular Docking Study

Cytarabine (Ara-C) is a synthetic isomer of cytidine that differs from cytidine and deoxycytidine only in the sugar. The use of arabinose instead of deoxyribose hinders the formation of phosphodiester linkages between pentoses, preventing the DNA chain from elongation and interrupting the DNA synthesis. The minor structural alteration (the inversion of hydroxyl at the 2' positions of the sugar) leads to change of the biological activity from anti-depressant and DNA/RNA block builder to powerful anti-cancer. Our study aimed to determine the molecular nature of this phenomenon. Three 1H-14N NMR-NQR experimental techniques, followed by solid-state computational modelling (Quantum Theory of Atoms in Molecules, Reduced Density Gradient and 3D Hirshfeld surfaces), Quantitative Structure-Property Relationships, Spackman's Hirshfeld surfaces and Molecular Docking were used. Multifaceted analysis-combining experiments, computational modeling and molecular docking-provides deep insight into three-dimensional packing at the atomic and molecular levels, but is challenging. A spectrum with nine lines indicating the existence of three chemically inequivalent nitrogen sites in the Ara-C molecule was recorded, and the lines were assigned to them. The influence of the structural alteration on the NQR parameters was modeled in the solid (GGA/RPBE). For the comprehensive description of the nature of these interactions several factors were considered, including relative reactivity and the involvement of heavy atoms in various non-covalent interactions. The binding modes in the solid state and complex with dCK were investigated using the novel approaches: radial plots, heatmaps and root-mean-square deviation of the binding mode. We identified the intramolecular OH···O hydrogen bond as the key factor responsible for forcing the glycone conformation and strengthening NH···O bonds with Gln97, Asp133 and Ara128, and stacking with Phe137. The titular butterfly effect is associated with both the inversion and the presence of this intramolecular hydrogen bond. Our study elucidates the differences in the binding modes of Ara-C and cytidine, which should guide the design of more potent anti-cancer and anti-viral analogues.

Proxy-approach in understanding the bisubstrate activity of strictosidine synthases

Besides tryptamine (1) and secologanin (2), non-cognate substrates also undergo a Pictet-Spengler reaction (PSR) catalyzed by strictosidine synthases (STR) with differing catalytic properties. We characterized the bisubstrate binding aspect of catalysis - order, affinity, and cooperativity - with STR orthologs from Rauvolfia serpentina (RsSTR) and Ophiorrhiza pumila (OpSTR) by an isothermal titration calorimetry (ITC) based 'proxy approach' that employed a non-reactive tryptamine analog (m1) to capture its inert ternary complexes with STRs and (2). ITC studies with OpSTR and (2) revealed 'tryptamine-first' cooperative binding with (1) and a simultaneous cooperative binding with (m1). Binding cooperativity among (m1) and (2) towards OpSTR was higher than RsSTR. Crystallographic study of RsSTR-(m1) complex helped to understand the unreactive binding of (m1) in terms of orientation and interactions in the RsSTR pocket. PSR with (m1) was revealed to be energetically unfeasible by the density functional theory (DFT) scans of the first hydrogen abstraction by RsSTR. The effect of pH on the bisubstrate binding to OpSTR was deciphered by molecular dynamics simulations (MDS), which also provided a molecular basis for the stability of complex of OpSTR with (m1) and (2). Therefore, we investigated STRs from a substrate binding perspective to inform drug-design and rational enzyme engineering efforts.

Binding order of substrate and cofactor in sulfonamide monooxygenase during sulfa drug degradation: in silico studies

For decades, sulfonamide antibiotics have been used across industries such as agriculture and animal husbandry. However, the use and inadvertent misuse of these antibiotics have resulted in the advent of sulfonamide-drug-resistant strains due to antibiotic pollution. Enzymatic bioremediation of antibiotics remains a potential emerging solution to combat antibiotic pollution. Here, we propose an enzymatic model for the degradation of sulfonamides by Microbacterium sp. We have employed a multi-pronged computational strategy involving - protein structure modelling, ligand docking and molecular dynamics simulations to decipher a plausible binding order for the enzymatic degradation of sulfonamides by the bacterial sulfonamide monooxygenase, SulX. Our results enable us to predict that this degradation is achieved through the sequential binding of the antibiotic sulfonamide followed by the reduced flavin cofactor FMNH2, thereby laying the computational foundation for further advancements in enzyme-mediated degradation of the antibiotic. We also provide a list of experiments which may be performed to verify and follow-up on our in-silico studies.Communicated by Ramaswamy H. Sarma.

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.

Exploring the impact of CYP11A1's missense SNPs on the interaction between CYP11A1 and cholesterol: A comprehensive structural analysis and MD simulation study

The process of steroidogenesis plays a vital role in human physiology as it governs the biosynthesis of mineralocorticoids, glucocorticoids, and androgens. These three classes of steroid hormones are primarily produced in the adrenal and gonadal glands through steroidogenesis pathways. Initiated by the side chain cleavage of cholesterol (CLR), this process leads to the conversion of cholesterol into pregnenolone and isocaproic aldehyde. The enzyme CYP11A1, encoded by the CYP11A1 gene, plays a key role in catalyzing the side chain cleavage of CLR. Several single nucleotide polymorphisms (SNPs) have been identified in the CYP11A1 gene, which may predispose carriers to disorders associated with abnormal steroidogenesis. Specifically, missense SNPs in the CYP11A1 gene have the potential to negatively impact the interaction between CYP11A1 and CLR, thus affecting the overall metabolome of steroid hormones. In this computational study, we focused on a specific set of missense SNPs reported in the CYP11A1 gene, aiming to identify variants that directly impact the interaction between CYP11A1 and CLR. The three-dimensional structure of the CYP11A1-CLR complex was obtained from the RCSB Protein Data Bank, while missense SNPs in the CYP11A1 gene were retrieved from Ensembl. To predict the most deleterious variants, we utilized the ConSurf server, SIFT, and PolyPhen. Furthermore, we assessed the impact of induced amino acid (AA) substitutions on the CYP11A1-CLR interaction using the PRODIGY server, PyMol, and Ligplot programs. Additionally, molecular dynamics (MD) simulations were conducted to analyze the effects of deleterious variants on the structural dynamics of the CYP11A1-CLR complex. Among the 8096 retrieved variants, we identified ten missense SNPs (E91K, W147G, R151W, R151Q, S391C, V392M, Q395K, Q416E, R460W, and R460Q) as deleterious for the interaction between CYP11A1 and CLR. MD simulations of the CYP11A1-CLR complexes carrying these deleterious AA substitutions revealed that Q416E, W147G, R460Q, and R460W had the most pronounced impacts on the structural dynamics of the complex. Consequently, these missense SNPs were considered the most deleterious ones. Further functional tests are recommended to assess the impact of these four missense SNPs on the enzymatic activity of CYP11A1. Moreover, Genome-Wide Association Studies (GWAS) should be conducted to determine the significance of their association with abnormal steroidogenesis diseases in various patient groups.

The anti-allergic potential of stingless bee honey from different botanical sources via modulation of mast cell degranulation

BACKGROUND

Allergy is an inflammatory disorder affecting around 20% of the global population. The adverse effects of current conventional treatments give rise to the increased popularity of using natural food products as complementary and alternative medicine against allergic diseases. Stingless bee honey, commonly known as Kelulut honey (KH) in Malaysia, has been used locally as a traditional remedy to relieve cough and asthma. This study evaluated the anti-allergic potential of KH collected from four different botanical sources on phorbol ester 12-myristate-3-acetate and calcium ionophore-activated human mast cells.

METHODS

The present study examined the inhibitory effects of all collected honey on the release of selected inflammatory mediators, such as tumor necrosis factor-α (TNF-α), interleukin (IL)-4, IL-6, IL-8, histamine, and β-hexosaminidase in an activated HMC. Besides that, all honey's total phenolic content (TPC) was also examined, followed by using liquid chromatography with tandem mass spectrometry (LC-MS/MS) to identify the phytochemicals in the honey. Further examination of the identified phytochemicals on their potential interaction with selected signaling molecules in an activated mast cell was conducted using computational methods.

RESULTS

The results indicated that there were significant inhibitory effects on all selected inflammatory mediators' release by KH sourced from bamboo (BH) and rubber tree (RH) at 0.5% and 1%, but not KH sourced from mango (AH) and noni (EH). BH and RH were found to have higher TPC values and were rich in their phytochemical profiles based on the LC-MS/MS results. Computational studies were employed to determine the possible molecular target of KH through molecular docking using HADDOCK and PRODIGY web servers.

CONCLUSIONS

In short, the results indicated that KH possesses anti-allergic effects towards an activated HMC, possibly by targeting downstream MAPKs. However, their anti-allergic effects may vary according to their botanical sources. Nevertheless, the present study has provided insight into the potential application of stingless bee honey as a complementary and alternative medicine to treat various allergic diseases.

Non-target biotransformation enzymes as a target for triazole-zinc mixtures

Triazoles inhibit lanosterol 14α-demethylase and block ergosterol biosynthesis in fungal pathogens. However, they also interact with other cytochrome P450 enzymes and influence non-target metabolic pathways. Disturbingly, triazoles may interact with essential elements. The interaction of penconazole (Pen), cyproconazole (Cyp) and tebuconazole (Teb) with Zn2+ results in the formation of deprotonated ligands in their complexes or in the creation of complexes with Cl- as a counterion or doubly charged complexes. Triazoles, as well as their equimolar cocktails with Zn2+ (10-6 mol/L), decreased the activities of the non-target enzymes CYP19A1 and CYP3A4. Pen most decreased CYP19A1 activity and was best bound to its active centre to block the catalytic cycle in computational analysis. For CYP3A4, Teb was found to be the most effective inhibitor by both, activity assay and interaction with the active centre. Teb/Cyp/Zn2+ and Teb/Pen/Cyp/Zn2+ cocktails also decreased the CYP19A1 activity, which was in correlation with the formation of numerous triazole-Zn2+ complexes.

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.

In Silico Identification of Potential Inhibitors of the SARS-CoV-2 Main Protease among a PubChem Database of Avian Infectious Bronchitis Virus 3CLPro Inhibitors

Remarkable structural homologies between the main proteases of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) and the avian infectious bronchitis virus (IBV) were revealed by comparative amino-acid sequence and 3D structural alignment. Assessing whether reported IBV 3CLPro inhibitors could also interact with SARS-CoV-2 has been undertaken in silico using a PubChem BioAssay database of 388 compounds active on the avian infectious bronchitis virus 3C-like protease. Docking studies of this database on the SARS-CoV-2 protease resulted in the identification of four covalent inhibitors targeting the catalytic cysteine residue and five non-covalent inhibitors for which the binding was further investigated by molecular dynamics (MD) simulations. Predictive ADMET calculations on the nine compounds suggest promising pharmacokinetic properties.

Computational analysis of missense variant CYP4F2*3 (V433M) in association with human CYP4F2 dysfunction: a functional and structural impact

BACKGROUND

Cytochrome P450 4F2 (CYP4F2) enzyme is a member of the CYP4 family responsible for the metabolism of fatty acids, therapeutic drugs, and signaling molecules such as arachidonic acid, tocopherols, and vitamin K. Several reports have demonstrated that the missense variant CYP4F2*3 (V433M) causes decreased activity of CYP4F2 and inter-individual variations in warfarin dose in different ethnic groups. However, the molecular pathogenicity mechanism of missense V433M in CYP4F2 at the atomic level has not yet been completely elucidated.

METHODS AND RESULTS

In the current study, we evaluated the effect of the V433M substitution on CYP4F2 using 14 different bioinformatics tools. Further molecular dynamics (MD) simulations were performed to assess the impact of the V433M mutation on the CYP4F2 protein structure, stability, and dynamics. In addition, molecular docking was used to illustrate the effect of V433M on its interaction with vitamin K1. Based on our results, the CYP4F2*3 variant was a damaging amino acid substitution with a destabilizing nature. The simulation results showed that missense V433M affects the dynamics and stability of CYP4F2 by reducing its compactness and stability, which means that it tends to change the overall structural conformation and flexibility of CYP4F2. The docking results showed that the CYP4F2*3 variant decreased the binding affinity between vitamin K1 and CYP4F2, which reduced the activity of CYP4F2*3 compared to native CYP4F2.

CONCLUSIONS

This study determined the molecular pathogenicity mechanism of the CYP4F2*3 variant on the human CYP4F2 protein and provided new information for understanding the structure-function relationship of CYP4F2 and other CYP4 enzymes. These findings will aid in the development of effective drugs and treatment options.

Discovery of potential FGFR3 inhibitors via QSAR, pharmacophore modeling, virtual screening and molecular docking studies against bladder cancer

BACKGROUND

Fibroblast growth factor receptor 3 is known as a favorable aim in vast range of cancers, particularly in bladder cancer treatment. Pharmacophore and QSAR modeling approaches are broadly utilized for developing novel compounds for the determination of inhibitory activity versus the biological target. In this study, these methods employed to identify FGFR3 potential inhibitors.

METHODS

To find the potential compounds for bladder cancer targeting, ZINC and NCI databases were screened. Pharmacophore and QSAR modeling of FGFR3 inhibitors were utilized for dataset screening. Then, with regard to several factors such as Absorption, Distribution, Metabolism, Excretion and Toxicity (ADMET) properties and Lipinski's Rule of Five, the recognized compounds were filtered. In further step, utilizing the flexible docking technique, the obtained compounds interactions with FGFR3 were analyzed.

RESULTS

The best five compounds, namely ZINC09045651, ZINC08433190, ZINC00702764, ZINC00710252 and ZINC00668789 were selected for Molecular Dynamics (MD) studies. Off-targeting of screened compounds was also investigated through CDD search and molecular docking. MD outcomes confirmed docking investigations and revealed that five selected compounds could make steady interactions with the FGFR3 and might have effective inhibitory potencies on FGFR3.

CONCLUSION

These compounds can be considered as candidates for bladder cancer therapy with improved therapeutic properties and less adverse effects.

SARS-CoV-2 proteases Mpro and PLpro: Design of inhibitors with predicted high potency and low mammalian toxicity using artificial neural networks, ligand-protein docking, molecular dynamics simulations, and ADMET calculations

The main (Mpro) and papain-like (PLpro) proteases are highly conserved viral proteins essential for replication of the COVID-19 virus, SARS-COV-2. Therefore, a logical plan for producing new drugs against this pathogen is to discover inhibitors of these enzymes. Accordingly, the goal of the present work was to devise a computational approach to design, characterize, and select compounds predicted to be potent dual inhibitors - effective against both Mpro and PLpro. The first step employed LigDream, an artificial neural network, to create a virtual ligand library. Ligands with computed ADMET profiles indicating drug-like properties and low mammalian toxicity were selected for further study. Initial docking of these ligands into the active sites of Mpro and PLpro was done with GOLD, and the highest-scoring ligands were redocked with AutoDock Vina to determine binding free energies (ΔG). Compounds 89-00, 89-07, 89-32, and 89-38 exhibited favorable ΔG values for Mpro (-7.6 to -8.7 kcal/mol) and PLpro (-9.1 to -9.7 kcal/mol). Global docking of selected compounds with the Mpro dimer identified prospective allosteric inhibitors 89-00, 89-27, and 89-40 (ΔG -8.2 to -8.9 kcal/mol). Molecular dynamics simulations performed on Mpro and PLpro active site complexes with the four top-scoring ligands from Vina demonstrated that the most stable complexes were formed with compounds 89-32 and 89-38. Overall, the present computational strategy generated new compounds with predicted drug-like characteristics, low mammalian toxicity, and high inhibitory potencies against both target proteases to form stable complexes. Further preclinical studies will be required to validate the in silico findings before the lead compounds could be considered for clinical trials.

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.

Cortisol Interaction with Aquaporin-2 Modulates Its Water Permeability: Perspectives for Non-Genomic Effects of Corticosteroids

Aquaporins (AQPs) are water channels widely distributed in living organisms and involved in many pathophysiologies as well as in cell volume regulations (CVR). In the present study, based on the structural homology existing between mineralocorticoid receptors (MRs), glucocorticoid receptors (GRs), cholesterol consensus motif (CCM) and the extra-cellular vestibules of AQPs, we investigated the binding of corticosteroids on the AQP family through in silico molecular dynamics simulations of AQP2 interactions with cortisol. We propose, for the first time, a putative AQPs corticosteroid binding site (ACBS) and discussed its conservation through structural alignment. Corticosteroids can mediate non-genomic effects; nonetheless, the transduction pathways involved are still misunderstood. Moreover, a growing body of evidence is pointing toward the existence of a novel membrane receptor mediating part of these rapid corticosteroids' effects. Our results suggest that the naturally produced glucocorticoid cortisol inhibits channel water permeability. Based on these results, we propose a detailed description of a putative underlying molecular mechanism. In this process, we also bring new insights on the regulatory function of AQPs extra-cellular loops and on the role of ions in tuning the water permeability. Altogether, this work brings new insights into the non-genomic effects of corticosteroids through the proposition of AQPs as the membrane receptor of this family of regulatory molecules. This original result is the starting point for future investigations to define more in-depth and in vivo the validity of this functional model.

Molecular determinants of acrylamide neurotoxicity through covalent docking

Acrylamide (ACR) is formed during food processing by Maillard reaction between sugars and proteins at high temperatures. It is also used in many industries, from water waste treatment to manufacture of paper, fabrics, dyes and cosmetics. Unfortunately, cumulative exposure to acrylamide, either from diet or at the workplace, may result in neurotoxicity. Such adverse effects arise from covalent adducts formed between acrylamide and cysteine residues of several neuronal proteins via a Michael addition reaction. The molecular determinants of acrylamide reactivity and its impact on protein function are not completely understood. Here we have compiled a list of acrylamide protein targets reported so far in the literature in connection with neurotoxicity and performed a systematic covalent docking study. Our results indicate that acrylamide binding to cysteine is favored in the presence of nearby positively charged amino acids, such as lysines and arginines. For proteins with more than one reactive Cys, docking scores were able to discriminate between the primary ACR modification site and secondary sites modified only at high ACR concentrations. Therefore, docking scores emerge as a potential filter to predict Cys reactivity against acrylamide. Inspection of the ACR-protein complex structures provides insights into the putative functional consequences of ACR modification, especially for non-enzyme proteins. Based on our study, covalent docking is a promising computational tool to predict other potential protein targets mediating acrylamide neurotoxicity.

Remdesivir analogs against SARS-CoV-2 RNA-dependent RNA polymerase

The COVID-19 pandemic has already taken many lives but is still continuing its spread and exerting jeopardizing effects. This study is aimed to find the most potent ligands from 703 analogs of remdesivir against RNA-dependent RNA polymerase (RdRp) protein of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) virus . RdRp is a major part of a multi-subunit transcription complex of the virus, which is essential for viral replication. In clinical trials, it has been found that remdesivir is effective to inhibit viral replication in Ebola and in primary human lung cell cultures; it effectively impedes replication of a broad-spectrum pre-pandemic bat coronaviruses and epidemic human coronaviruses. After virtual screening, 30 most potent ligands and remdesivir were modified with triphosphate. Quantum mechanics-based quantitative structure-activity relationship envisages the binding energy for ligands applying partial least square (PLS) regression. PLS regression remarkably predicts the binding energy of the effective ligands with an accuracy of 80% compared to the value attained from molecular docking. Two ligands (L4:58059550 and L28:126719083), which have more interactions with the target protein than the other ligands including standard remdesivir triphosphate, were selected for further analysis. Molecular dynamics simulation is done to assess the stability and dynamic nature of the drug-protein complex. Binding-free energy results via PRODIGY server and molecular mechanics/Poisson-Boltzmann surface area method depict that the potential and solvation energies play a crucial role. Considering all computational analysis, we recommend the best remdesivir analogs can be utilized for efficacy test through in vitro and in vivo trials against SARS-CoV-2.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.

Structural and energetic analyses of SARS-CoV-2 N-terminal domain characterise sugar binding pockets and suggest putative impacts of variants on COVID-19 transmission

Coronavirus disease 2019 (COVID-19) caused by SARS-CoV-2 is an ongoing pandemic that causes significant health/socioeconomic burden. Variants of concern (VOCs) have emerged affecting transmissibility, disease severity and re-infection risk. Studies suggest that the - N-terminal domain (NTD) of the spike protein may have a role in facilitating virus entry via sialic-acid receptor binding. Furthermore, most VOCs include novel NTD variants. Despite global sequence and structure similarity, most sialic-acid binding pockets in NTD vary across coronaviruses. Our work suggests ongoing evolutionary tuning of the sugar-binding pockets and recent analyses have shown that NTD insertions in VOCs tend to lie close to loops. We extended the structural characterisation of these sugar-binding pockets and explored whether variants could enhance sialic acid-binding. We found that recent NTD insertions in VOCs (i.e., Gamma, Delta and Omicron variants) and emerging variants of interest (VOIs) (i.e., Iota, Lambda and Theta variants) frequently lie close to sugar-binding pockets. For some variants, including the recent Omicron VOC, we find increases in predicted sialic acid-binding energy, compared to the original SARS-CoV-2, which may contribute to increased transmission. These binding observations are supported by molecular dynamics simulations (MD). We examined the similarity of NTD across Betacoronaviruses to determine whether the sugar-binding pockets are sufficiently similar to be exploited in drug design. Whilst most pockets are too structurally variable, we detected a previously unknown highly structurally conserved pocket which can be investigated in pursuit of a generic pan-Betacoronavirus drug. Our structure-based analyses help rationalise the effects of VOCs and provide hypotheses for experiments. Our findings suggest a strong need for experimental monitoring of changes in NTD of VOCs.

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.

FIBCD1 is an endocytic GAG receptor associated with a novel neurodevelopmental disorder

Whole-exome sequencing of two patients with idiopathic complex neurodevelopmental disorder (NDD) identified biallelic variants of unknown significance within FIBCD1, encoding an endocytic acetyl group-binding transmembrane receptor with no known function in the central nervous system. We found that FIBCD1 preferentially binds and endocytoses glycosaminoglycan (GAG) chondroitin sulphate-4S (CS-4S) and regulates GAG content of the brain extracellular matrix (ECM). In silico molecular simulation studies and GAG binding analyses of patient variants determined that such variants are loss-of-function by disrupting FIBCD1-CS-4S association. Gene knockdown in flies resulted in morphological disruption of the neuromuscular junction and motor-related behavioural deficits. In humans and mice, FIBCD1 is expressed in discrete brain regions, including the hippocampus. Fibcd1 KO mice exhibited normal hippocampal neuronal morphology but impaired hippocampal-dependent learning. Further, hippocampal synaptic remodelling in acute slices from Fibcd1 KO mice was deficient but restored upon enzymatically modulating the ECM. Together, we identified FIBCD1 as an endocytic receptor for GAGs in the brain ECM and a novel gene associated with an NDD, revealing a critical role in nervous system structure, function and plasticity.

MDO: A Computational Protocol for Prediction of Flexible Enzyme-Ligand Binding Mode

AIM

Developing a method for use in computer aided drug design Background: Predicting the structure of enzyme-ligand binding mode is essential for understanding the properties, functions, and mechanisms of the bio-complex, but is rather difficult due to the enormous sampling space involved.

OBJECTIVE

Accurate prediction of enzyme-ligand binding mode conformation.

METHOD

A new computational protocol, MDO, is proposed for finding the structure of ligand binding pose. MDO consists of sampling enzyme sidechain conformations via molecular dynamics simulation of enzyme-ligand system and clustering of the enzyme configurations, sampling ligand binding poses via molecular docking and clustering of the ligand conformations, and the optimal ligand binding pose prediction via geometry optimization and ranking by the ONIOM method. MDO is tested on 15 enzyme-ligand complexes with known accurate structures.

RESULTS

The success rate of MDO predictions, with RMSD < 2 Å, is 67%, substantially higher than the 40% success rate of conventional methods. The MDO success rate can be increased to 83% if the ONIOM calculations are applied only for the starting poses with ligands inside the binding cavities.

CONCLUSION

The MDO protocol provides high quality enzyme-ligand binding mode prediction with reasonable computational cost. The MDO protocol is recommended for use in the structure-based drug design.

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.

Hydrazones and Thiosemicarbazones Targeting Protein-Protein-Interactions of SARS-CoV-2 Papain-like Protease

The papain-like protease (PLpro) of SARS-CoV-2 is essential for viral propagation and, additionally, dysregulation of the host innate immune system. Using a library of 40 potential metal-chelating compounds we performed an X-ray crystallographic screening against PLpro. As outcome we identified six compounds binding to the target protein. Here we describe the interaction of one hydrazone (H1) and five thiosemicarbazone (T1-T5) compounds with the two distinct natural substrate binding sites of PLpro for ubiquitin and ISG15. H1 binds to a polar groove at the S1 binding site by forming several hydrogen bonds with PLpro. T1-T5 bind into a deep pocket close to the polyubiquitin and ISG15 binding site S2. Their interactions are mainly mediated by multiple hydrogen bonds and further hydrophobic interactions. In particular compound H1 interferes with natural substrate binding by sterical hindrance and induces conformational changes in protein residues involved in substrate binding, while compounds T1-T5 could have a more indirect effect. Fluorescence based enzyme activity assay and complementary thermal stability analysis reveal only weak inhibition properties in the high micromolar range thereby indicating the need for compound optimization. Nevertheless, the unique binding properties involving strong hydrogen bonding and the various options for structural optimization make the compounds ideal lead structures. In combination with the inexpensive and undemanding synthesis, the reported hydrazone and thiosemicarbazones represent an attractive scaffold for further structure-based development of novel PLpro inhibitors by interrupting protein-protein interactions at the S1 and S2 site.

A Perspective for Ménière’s Disease: In Silico Investigations of Dexamethasone as a Direct Modulator of AQP2

Ménière’s disease is a chronic illness characterized by intermittent episodes of vertigo associated with fluctuating sensorineural hearing loss, tinnitus and aural pressure. This pathology strongly correlates with a dilatation of the fluid compartment of the endolymph, so-called hydrops. Dexamethasone is one of the therapeutic approaches recommended when conventional antivertigo treatments have failed. Several mechanisms of actions have been hypothesized for the mode of action of dexamethasone, such as the anti-inflammatory effect or as a regulator of inner ear water homeostasis. However, none of them have been experimentally confirmed so far. Aquaporins (AQPs) are transmembrane water channels and are hence central in the regulation of transcellular water fluxes. In the present study, we investigated the hypothesis that dexamethasone could impact water fluxes in the inner ear by targeting AQP2. We addressed this question through molecular dynamics simulations approaches and managed to demonstrate a direct interaction between AQP2 and dexamethasone and its significant impact on the channel water permeability. Through compartmentalization of sodium and potassium ions, a significant effect of Na+ upon AQP2 water permeability was highlighted as well. The molecular mechanisms involved in dexamethasone binding and in its regulatory action upon AQP2 function are described.

Biophysical, Biochemical, and Molecular Docking Investigations of Anti-Glycating, Antioxidant, and Protein Structural Stability Potential of Garlic

Garlic has been reported to inhibit protein glycation, a process that underlies several disease processes, including chronic complications of diabetes mellitus. Biophysical, biochemical, and molecular docking investigations were conducted to assess anti-glycating, antioxidant, and protein structural protection activities of garlic. Results from spectral (UV and fluorescence) and circular dichroism (CD) analysis helped ascertain protein conformation and secondary structure protection against glycation to a significant extent. Further, garlic showed heat-induced protein denaturation inhibition activity (52.17%). It also inhibited glycation, advanced glycation end products (AGEs) formation as well as lent human serum albumin (HSA) protein structural stability, as revealed by reduction in browning intensity (65.23%), decrease in protein aggregation index (67.77%), and overall reduction in cross amyloid structure formation (33.26%) compared with positive controls (100%). The significant antioxidant nature of garlic was revealed by FRAP assay (58.23%) and DPPH assay (66.18%). Using molecular docking analysis, some of the important garlic metabolites were investigated for their interactions with the HSA molecule. Molecular docking analysis showed quercetin, a phenolic compound present in garlic, appears to be the most promising inhibitor of glucose interaction with the HSA molecule. Our findings show that garlic can prevent oxidative stress and glycation-induced biomolecular damage and that it can potentially be used in the treatment of several health conditions, including diabetes and other inflammatory diseases.

Benchmarking ensemble docking methods in D3R Grand Challenge 4

The discovery of new drugs is a time consuming and expensive process. Methods such as virtual screening, which can filter out ineffective compounds from drug libraries prior to expensive experimental study, have become popular research topics. As the computational drug discovery community has grown, in order to benchmark the various advances in methodology, organizations such as the Drug Design Data Resource have begun hosting blinded grand challenges seeking to identify the best methods for ligand pose-prediction, ligand affinity ranking, and free energy calculations. Such open challenges offer a unique opportunity for researchers to partner with junior students (e.g., high school and undergraduate) to validate basic yet fundamental hypotheses considered to be uninteresting to domain experts. Here, we, a group of high school-aged students and their mentors, present the results of our participation in Grand Challenge 4 where we predicted ligand affinity rankings for the Cathepsin S protease, an important protein target for autoimmune diseases. To investigate the effect of incorporating receptor dynamics on ligand affinity rankings, we employed the Relaxed Complex Scheme, a molecular docking method paired with molecular dynamics-generated receptor conformations. We found that Cathepsin S is a difficult target for molecular docking and we explore some advanced methods such as distance-restrained docking to try to improve the correlation with experiments. This project has exemplified the capabilities of high school students when supported with a rigorous curriculum, and demonstrates the value of community-driven competitions for beginners in computational drug discovery.

Synthesis, In Vitro, In Vivo and In Silico Antidiabetic Bioassays of 4-Nitro(thio)phenoxyisobutyric Acids Acting as Unexpected PPARγ Modulators: An In Combo Study

Four isobutyric acids (two nitro and two acetamido derivatives) were prepared in two steps and characterized using spectral analysis. The mRNA concentrations of PPARγ and GLUT-4 (two proteins documented as key diabetes targets) were increased by 3T3-L1 adipocytes treated with compounds 1–4, but an absence of in vitro expression of PPARα was observed. Docking and molecular dynamics studies revealed the plausible interaction between the synthesized compounds and PPARγ. In vivo studies established that compounds 1–4 have antihyperglycemic modes of action associated with insulin sensitization. Nitrocompound 2 was the most promising of the series, being orally active, and one of multiple modes of action could be selective PPARγ modulation due to its extra anchoring with Gln-286. In conclusion, we demonstrated that nitrocompound 2 showed strong in vitro and in vivo effects and can be considered as an experimental antidiabetic candidate.

Characterization of recombinant pumpkin 2S albumin and mutation studies to unravel potential DNA/RNA binding site

The native pumpkin 2S albumin, a multifunctional protein, possess a variety of potential biotechnologically exploitable properties. The present study reports the characterization of recombinant pumpkin 2S albumin (rP2SA) and unraveling of its potential DNA/RNA binding site. The purification and characterization of the rP2SA established that it retains the characteristic α-helical structure and exhibited comparable DNase, RNase, antifungal and anti-proliferative activities as native protein. In vitro studies revealed that rP2SA exhibits potent antiviral activity against chikungunya virus (CHIKV) at a non-toxic concentration with an IC₅₀ of 114.5 μg/mL. In silico studies and site-directed mutagenesis were employed to unravel the potential DNA/RNA binding site. A strong positive charge distribution due to presence of many arginine residues in proximity of helix 5 was identified as a potential site. The two of the arginine residues, conserved in some 2S albumins, were selected for the mutation studies. The mutated forms of recombinant protein (R84A and R91A) showed a drastic reduction in DNase and RNase activities suggesting their presence at binding site and involvement in the nuclease activity. A metal binding site was also identified adjacent to DNA/RNA binding site. The present study demonstrated the structural and functional integrity of the rP2SA and reports potential antiviral activity against CHIKV. Further, potential DNA/RNA binding site was unraveled through mutation studies and bioinformatics analysis.

Identification of fluoxetine as a direct NLRP3 inhibitor to treat atrophic macular degeneration

The atrophic form of age-related macular degeneration (dry AMD) affects nearly 200 million people worldwide. There is no Food and Drug Administration (FDA)-approved therapy for this disease, which is the leading cause of irreversible blindness among people over 50 y of age. Vision loss in dry AMD results from degeneration of the retinal pigmented epithelium (RPE). RPE cell death is driven in part by accumulation of Alu RNAs, which are noncoding transcripts of a human retrotransposon. Alu RNA induces RPE degeneration by activating the NLRP3-ASC inflammasome. We report that fluoxetine, an FDA-approved drug for treating clinical depression, binds NLRP3 in silico, in vitro, and in vivo and inhibits activation of the NLRP3-ASC inflammasome and inflammatory cytokine release in RPE cells and macrophages, two critical cell types in dry AMD. We also demonstrate that fluoxetine, unlike several other antidepressant drugs, reduces Alu RNA-induced RPE degeneration in mice. Finally, by analyzing two health insurance databases comprising more than 100 million Americans, we report a reduced hazard of developing dry AMD among patients with depression who were treated with fluoxetine. Collectively, these studies identify fluoxetine as a potential drug-repurposing candidate for dry AMD.

Shape-Restrained Modeling of Protein-Small-Molecule Complexes with High Ambiguity Driven DOCKing

Small-molecule docking remains one of the most valuable computational techniques for the structure prediction of protein-small-molecule complexes. It allows us to study the interactions between compounds and the protein receptors they target at atomic detail in a timely and efficient manner. Here, we present a new protocol in HADDOCK (High Ambiguity Driven DOCKing), our integrative modeling platform, which incorporates homology information for both receptor and compounds. It makes use of HADDOCK's unique ability to integrate information in the simulation to drive it toward conformations, which agree with the provided data. The focal point is the use of shape restraints derived from homologous compounds bound to the target receptors. We have developed two protocols: in the first, the shape is composed of dummy atom beads based on the position of the heavy atoms of the homologous template compound, whereas in the second, the shape is additionally annotated with pharmacophore data for some or all beads. For both protocols, ambiguous distance restraints are subsequently defined between those beads and the heavy atoms of the ligand to be docked. We have benchmarked the performance of these protocols with a fully unbound version of the widely used DUD-E (Database of Useful Decoys-Enhanced) dataset. In this unbound docking scenario, our template/shape-based docking protocol reaches an overall success rate of 81% when a reliable template can be identified (which was the case for 99 out of 102 complexes in the DUD-E dataset), which is close to the best results reported for bound docking on the DUD-E dataset.

Structural Biology in the Clouds: The WeNMR-EOSC Ecosystem

Structural biology aims at characterizing the structural and dynamic properties of biological macromolecules at atomic details. Gaining insight into three dimensional structures of biomolecules and their interactions is critical for understanding the vast majority of cellular processes, with direct applications in health and food sciences. Since 2010, the WeNMR project (www.wenmr.eu) has implemented numerous web-based services to facilitate the use of advanced computational tools by researchers in the field, using the high throughput computing infrastructure provided by EGI. These services have been further developed in subsequent initiatives under H2020 projects and are now operating as Thematic Services in the European Open Science Cloud portal (www.eosc-portal.eu), sending >12 millions of jobs and using around 4,000 CPU-years per year. Here we review 10 years of successful e-infrastructure solutions serving a large worldwide community of over 23,000 users to date, providing them with user-friendly, web-based solutions that run complex workflows in structural biology. The current set of active WeNMR portals are described, together with the complex backend machinery that allows distributed computing resources to be harvested efficiently.

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.

In silico and in vitro assessment of androgen receptor antagonists

There is a growing concern for male reproductive health as studies suggest that there is a sharp increase in prostate cancer and other fertility related problems. Apart from lifestyle, pollutants are also known to negatively affect the reproductive system. In addition to many other compounds that have been shown to alter androgen signaling, several environmental pollutants are known to disrupt androgen signaling via binding to androgen receptor (AR) or indirectly affecting the androgen synthesis. We analyzed here the molecular mechanism of the interaction between the human AR Ligand Binding Domain (hAR-LBD) and two environmental pollutants, linuron (a herbicide) and procymidone (a pesticide), and compared with the steroid agonist dihydrotestosterone (DHT) and well-known hAR antagonists bicalutamide and enzalutamide. Using molecular docking and dynamics simulations, we showed that the co-activator interaction site of the hAR-LBD is disrupted in different ways by different ligands. Binding free energies of the ligands were also ordered in increasing order as follows: linuron, procymidone, DHT, bicalutamide, and enzalutamide. These data were confirmed by in vitro assays. Reporter assay with MDA-kb2 cells showed that linuron, procymidone, bicalutamide and enzalutamide can inhibit androgen mediated activation of luciferase activity. Gene expression analysis further showed that these compounds can inhibit the expression of prostate specific antigen (PSA) and microseminoprotein beta (MSMB) in prostate cell line LNCaP. Comparative analysis showed that procymidone is more potent than linuron in inhibiting AR activity. Furthermore, procymidone at 10 μM dose showed equivalent and higher activity to AR inhibitor enzalutamide and bicalutamide respectively.

Molecular Simplification of Natural Products: Synthesis, Antibacterial Activity, and Molecular Docking Studies of Berberine Open Models

Berberine, the main bioactive component of many medicinal plants belonging to various genera such as Berberis, Coptis, and Hydrastis is a multifunctional compound. Among the numerous interesting biological properties of berberine is broad antimicrobial activity including a range of Gram-positive and Gram-negative bacteria. With the aim of identifying berberine analogues possibly endowed with higher lead-likeness and easier synthetic access, the molecular simplification approach was applied to the secondary metabolite and a series of analogues were prepared and screened for their antimicrobial activity against Gram-positive and Gram-negative bacterial test species. Rewardingly, the berberine simplified analogues displayed 2-20-fold higher potency with respect to berberine. Since our berberine simplified analogues may be easily synthesized and are characterized by lower molecular weight than the parent compound, they are further functionalizable and should be more suitable for oral administration. Molecular docking simulations suggested FtsZ, a well-known protein involved in bacterial cell division, as a possible target.

Unraveling the mechanism of arbidol binding and inhibition of SARS-CoV-2: Insights from atomistic simulations

The COVID-19 pandemic has spread rapidly and poses an unprecedented threat to the global economy and human health. Broad-spectrum antivirals are currently being administered to treat severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). China's prevention and treatment guidelines suggest the use of an anti-influenza drug, arbidol, for the clinical treatment of COVID-19. Reports indicate that arbidol could neutralize SARS-CoV-2. Monotherapy with arbidol is superior to lopinavir-ritonavir or favipiravir for treating COVID-19. In SARS-CoV-2 infection, arbidol acts by interfering with viral binding to host cells. However, the detailed mechanism by which arbidol induces the inhibition of SARS-CoV-2 is not known. Here, we present atomistic insights into the mechanism underlying membrane fusion inhibition of SARS-CoV-2 by arbidol. Molecular dynamics (MD) simulation-based analyses demonstrate that arbidol binds and stabilizes at the receptor-binding domain (RBD)/ACE2 interface with a high affinity. It forms stronger intermolecular interactions with the RBD than ACE2. Analyses of the detailed decomposition of energy components and binding affinities revealed a substantial increase in the affinity between the RBD and ACE2 in the arbidol-bound RBD/ACE2 complex, suggesting that arbidol generates favorable interactions between them. Based on our MD simulation results, we propose that the binding of arbidol induces structural rigidity in the viral glycoprotein, thus restricting the conformational rearrangements associated with membrane fusion and virus entry. Furthermore, key residues of the RBD and ACE2 that interact with arbidol were identified, opening the door for developing therapeutic strategies and higher-efficacy arbidol derivatives or lead drug candidates.

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.

Insight into molecular characteristics of SARS-CoV-2 spike protein following D614G point mutation, a molecular dynamics study

Undoubtedly, the SARS-CoV-2 has become a major concern for all societies due to its catastrophic effects on public health. In addition, mutations and changes in the structure of the virus make it difficult to design effective treatment. Moreover, the amino acid sequence of a protein is a major factor in the formation of the second and tertiary structure in a protein. Amino acid replacement can have noticeable effects on the folding of a protein, especially if an asymmetric change (substitution of polar residue with non-polar, charged with an uncharged, positive charge with a negative charge, or large residue with small residue) occurs. D614G as a spike mutant of SARS-CoV-2 previously identified as an associated risk factor with a high mortality rate of this virus. Using structural bioinformatics, our group determined that D614G mutation could cause extensive changes in SARS-CoV-2 behavior including the secondary structure, receptor binding pattern, 3D conformation, and stability of it.

Communicated by Ramaswamy H. Sarma

How Far Are We from the Rapid Prediction of Drug Resistance Arising Due to Kinase Mutations?

In all living organisms, protein kinases regulate various cell signaling events through phosphorylation. The phosphorylation occurs upon transferring an ATP's terminal phosphate to a target residue. Because of the central role of protein kinases in several proliferative pathways, point mutations occurring within the kinase's ATP-binding site can lead to a constitutively active enzyme, and ultimately, to cancer. A select set of these point mutations can also make the enzyme drug resistant toward the available kinase inhibitors. Because of technical and economical limitations, rapid experimental exploration of the impact of these mutations remains to be a challenge. This underscores the importance of kinase-ligand binding affinity prediction tools that are poised to measure the efficacy of inhibitors in the presence of kinase mutations. To this end, here, we compare the performances of six web-based scoring tools (DSX-ONLINE, KDEEP, HADDOCK2.2, PDBePISA, Pose&Rank, and PRODIGY-LIG) in assessing the impact of kinase mutations on their interactions with their inhibitors. This assessment is carried out on a new structure-based BINDKIN benchmark we compiled. BINDKIN contains wild-type and mutant structure pairs of kinase-inhibitor complexes, together with their corresponding experimental binding affinities (in the form of IC₅₀, K _d, and K _i). The performance of various web servers over BINDKIN shows that they cannot predict the binding affinities (ΔGs) of wild-type and mutant cases directly. Still, they could catch whether a mutation improves or worsens the ligand binding (ΔΔGs) where the highest Pearson's R correlation coefficient is reached by DSX-ONLINE over the K _i dataset. When homology models are used instead of K _i-associated crystal structures, DSX-ONLINE loses its predictive capacity. These results highlight that there is room to improve the available scoring functions to estimate the impact of protein kinase point mutations on inhibitor binding. The BINDKIN benchmark with all related results is freely accessible online (https://github.com/CSB-KaracaLab/BINDKIN).

Enfuvirtide, an HIV-1 fusion inhibitor peptide, can act as a potent SARS-CoV-2 fusion inhibitor: an in silico drug repurposing study

Regarding the urgency of therapeutic measures for coronavirus disease 2019 (COVID-19) pandemic, the use of available drugs with FDA approval is preferred because of the less time and cost required for their development. In silico drug repurposing is an accurate way to speed up the screening of the existing FDA-approved drugs to find a therapeutic option for COVID-19. The similarity in SARS-CoV-2 and HIV-1 fusion mechanism to host cells can be a key point for Inhibit SARS-CoV-2 entry into host cells by HIV fusion inhibitors. Accordingly, in this study, an HIV-1 fusion inhibitor called Enfuvirtide (Enf) was selected. The affinity and essential residues involving in the Enf binding to the S2 protein of SARS-CoV-2, HIV-1 gp41 protein and angiotensin-converting enzyme 2 (ACE-2) as a negative control, was evaluated using molecular docking. Eventually, Enf-S2 and Enf-gp41 protein complexes were simulated by molecular dynamics (MD) in terms of binding affinity and stability. Based on the most important criteria such as docking score, cluster size, energy and dissociation constant, the strongest interaction was observed between Enf with the S2 protein. In addition, MD results confirmed that Enf-S2 protein interaction was remarkably stable and caused the S2 protein residues to undergo the fewest fluctuations. In conclusion, it can be stated that Enf can act as a strong SARS-CoV-2 fusion inhibitor and demonstrates the potential to enter the clinical trial phase of COVID-19.

Communicated by Ramaswamy H. Sarma