Determination of ensemble-average pairwise root mean-square deviation from experimental B-factors

In-Silico modulation of Interleukin-8 (IL8) for the therapeutic management of endodontic pulpitis

Emerging clinical evidences highlight the association of Interleukin-8 (IL8) with endodontic pulpitis. Relatively higher expression of IL8 has been found in the pulp samples of pulpitis patients with moderate/severe pain. It is speculated that IL8 can be considered as a potential target for therapeutics of endodontic pulpitis. A library consisting of 3072 small molecules from the ZINC database was used to identify potential lead molecules with drug-like properties against the IL8. Based on the in-silico structure-assisted drug designing involving molecular docking, MD simulations, and MMPBSA analyses, we found a small molecule ZINC14613097 inhibits IL8. This study provides a new lead molecule than can be further validated in in-vitro, in-vivo, and ongoing clinical studies for the therapeutic management of endodontic pulpitis.

Targeting Cytotoxin-Associated Antigen A, a Virulent Factor of Helicobacter pylori-Associated Gastric Cancer: Structure-Based In Silico Screening of Natural Compounds

Gastric cancer is the fifth most frequent cancer and the third major cause of mortality worldwide. Helicobacter pylori, a bacterial infection linked with GC, injects the cytotoxin-associated antigen A (CagA; an oncoprotein) into host cells. When the phosphorylated CagA protein enters the cell, it attaches to other cellular components, interfering with normal cellular signaling pathways. CagA plays an important role in the progression of GC by interacting with phosphatidylserine of the host cell membrane. Therefore, disrupting the CagA-phosphatidylserine connection using small molecules appears to be a promising therapeutic approach. In this report, we screened the natural compounds from ZINC database against the CagA protein using the bioinformatics tools. Hits were initially chosen based on their physicochemical, absorption, distribution, metabolism, excretion, and toxicity (ADMET) characteristics, as well as other drug-like characteristics. To locate safe and effective hits, the PAINS filter, binding affinities estimation, and interaction analysis were used. Three compounds with high binding affinity and specificity for the CagA binding pocket were discovered. The final hits, ZINC153731, ZINC69482055, and ZINC164387, were found to bind strongly with CagA protein, with binding energies of -11.53, -10.67, and -9.21 kcal/mol, respectively, which were higher than that of the control compound (-7.25 kcal/mol). Further, based on binding affinity and interaction pattern, two leads (ZINC153731, ZINC69482055) were chosen for molecular dynamics (MD) simulation analysis. MD results showed that they displayed stability in their vicinity at 100 ns. This study suggested that these compounds could be used as possible inhibitors of CagA protein in the fight against GC. However, additional benchwork tests are required to validate them as CagA protein inhibitors.

New Imidazole-Based N-Phenylbenzamide Derivatives as Potential Anticancer Agents: Key Computational Insights

An efficient atom-economical synthetic protocol to access new imidazole-based N-phenylbenzamide derivatives is described. A one-pot three-component reaction was utilized to provide a series of N-phenylbenzamide derivatives in a short reaction time (2–4 h) with an 80–85% yield. The cytotoxic evaluation revealed that derivatives 4e and 4f exhibited good activity, with IC₅₀ values between 7.5 and 11.1 μM against the tested cancer cell lines. Computational studies revealed interesting insights: the docking of the active derivatives (4e and 4f) showed a higher affinity toward the target receptor protein than the control. Molecular dynamic simulations revealed that the active derivatives form stable complexes with the ABL1 kinase protein. Moreover, the ADME and drug-likeness of the derivatives reinforced the potential of the derivatives to be taken up for further development as anticancer agents.

Methyltransferase as Antibiotics Against Foodborne Pathogens: An In Silico Approach for Exploring Enzyme as Enzymobiotics

The development of resistance in microbes against antibiotics and limited choice for the use of chemical preservatives in food lead the urgent need to search for an alternative to antibiotics. The enzymes are catalytic proteins that catalyze digestion of bacterial cell walls and protein requirements for the survival of the cell. To study methyltransferase as antibiotics against foodborne pathogen, the methyltransferase enzyme sequence was modeled and its interactions were analyzed against a membrane protein of the gram-positive and gram-negative bacteria through in silico protein-protein interactions. The methyltransferase interaction with cellular protein was found to be maximum, due to the maximum PatchDock Score (15808), which was followed by colicin (12864) and amoxicillin (4122). The modeled protein has found to be interact more significantly to inhibit the indicator bacteria than the tested antibiotics and antimicrobial colicin protein. Thus, model enzyme methyltransferase could be used as enzymobiotics. Moreover, peptide sequences similar to this enzyme sequence need to be designed and evaluated against the microbial pathogen.

Investigation of anti-diabetic potential and molecular simulation studies of dihydropyrimidinone derivatives

In an attempt to find new targets for α-amylase and α-glucosidase for the treatment of type 2 diabetes mellitus, the present study aims in determining the anti-diabetic potential of synthesized dihydropyrimidinone derivatives. The in vitro α-glucosidase and α-amylase inhibitory activity was performed and the molecular docking analysis of the ligand in the active binding site of target protein was determined. The results revealed significant percent inhibition of α-glucosidase by the compound 6-benzyl-4-(4-hydroxyphenyl)-3,4,6,7-tetrahydro-1H-pyrrolo[3,4-d]pyrimidine-2,5-dione (compound A). The active compound showed 81.99% inhibition when compared to standard ascorbic acid having percent inhibition 81.18%. The IC₅₀ of active compound (A) showed to be 1.02 µg/ml. The molecular docking analysis revealed that the ligand bound to the active binding site of protein with the lowest binding energy of -7.9 kcal/mol that was also significantly similar to standard having -7.8 kcal/mol binding energy. The molecular dynamic simulation studies also revealed stable binding of ligand in the active binding site of protein with low RMSD of 1.7 Å similar to the protein RMSD 1.6Å In conclusion, the study revealed a potential new target against α-glucosidase to treat type 2 diabetes mellitus.

Determination of tyrosinase-cyanidin-3-O-glucoside and (-/+)-catechin binding modes reveal mechanistic differences in tyrosinase inhibition

Tyrosinase, exquisitely catalyzes the phenolic compounds into brown or black pigment, inhibition is used as a treatment for dermatological or neurodegenerative disorders. Natural products, such as cyanidin-3-O-glucoside and (-/+)-catechin, are considered safe and non-toxic food additives in tyrosinase inhibition but their ambiguous inhibitory mechanism against tyrosinase is still elusive. Thus, we presented the mechanistic insights into tyrosinase with cyanidin-3-O-glucoside and (-/+)-catechin using computational simulations and in vitro assessment. Initial molecular docking results predicted ideal docked poses (- 9.346 to - 5.795 kcal/mol) for tyrosinase with selected flavonoids. Furthermore, 100 ns molecular dynamics simulations and post-simulation analysis of docked poses established their stability and oxidation of flavonoids as substrate by tyrosinase. Particularly, metal chelation via catechol group linked with the free 3-OH group on the unconjugated dihydropyran heterocycle chain was elucidated to contribute to tyrosinase inhibition by (-/+)-catechin against cyanidin-3-O-glucoside. Also, predicted binding free energy using molecular mechanics/generalized Born surface area for each docked pose was consistent with in vitro enzyme inhibition for both mushroom and murine tyrosinases. Conclusively, (-/+)-catechin was observed for substantial tyrosinase inhibition and advocated for further investigation for drug development against tyrosinase-associated diseases.

Determination of Novel Anti-Cancer Agents by Targeting OGG1 Enzyme Using Integrated Bioinformatics Methods

The 8-oxoguanine DNA glycosylase (OGG1) enzyme is a key DNA glycosylase mediating the excision of 7,8-dihydro-8-oxoguanine (8-oxoG) from DNA molecule to the start base excision repair pathway. The OGG1 glycosylase function depletion has been seen to obstruct pathological conditions such as inflammation, A3 T-cell lymphoblastic acute leukemia growth, and neurodegenerative diseases, thus warranting OGG1 as an attractive anti-cancer enzyme. Herein, we employed several drug libraries intending to screen non-toxic inhibitory molecules against the active pocket of the enzyme that achieved stable binding mode in dynamics. Two anti-cancer compounds ([O-]C1=C(CC2=CC=CC=C2)SC(=[N+]1CC(=O)NC3=NC=C(CC4=CC=CC=C4)S3)S and CCCN(CCC)[S]-(=O)(=O)C1=CC=C(C=C1)C(=O)NNC2=NC3=CC=C(Br)C=C3C(=N2)C4=CC=CC=C4) from Selleckchem.com were identified to occupy the active pocket of OGG1 and bind with greater affinity than Control TH5487. The binding affinity of Top-1 is -11.6 kcal/mol while that of Top-2 is -10.7 kcal/mol in contrast to TH5487 Control (-9 kcal/mol). During molecular dynamic simulations versus time, the said compounds are tightly held by the enzyme with no minor structural deviations reported except flexible loops in particular those present at the N and C-terminal. Both the compounds produced extensive hydrophobic interactions with the enzyme along with stable hydrogen bonding. The docking and molecular dynamics simulations predictions were further validated by molecular mechanics with generalized Born and surface area solvation (MM/GBSA) and Poisson Boltzmann surface area (MM/PBSA), and WaterSwap binding energies that validated strong binding of the compounds to the enzyme. The MM/GBSA binding free energy for Top-1 complex is -28.10 kcal/mol, Top-2 complex is -50.14 kcal/mol) and Control is -46.91 kcal/mol while MM/PBSA value for Top-1, Top-2 and Control is -23.38 kcal/mol, -35.29 kcal/mol and -38.20 kcal/mol, respectively. Computational pharmacokinetics support good druglike candidacy of the compounds with acceptable profile of pharmacokinetics and very little toxicity. All these findings support the notion that the compounds can be used in experiments to test their anti-cancer activities.

Disruption of Water Networks is the Cause of Human/Mouse Species Selectivity in Urokinase Plasminogen Activator (uPA) Inhibitors Derived from Hexamethylene Amiloride (HMA)

The urokinase plasminogen activator (uPA) plays a critical role in tumor cell invasion and migration and is a promising antimetastasis target. 6-Substituted analogues of 5-N,N-(hexamethylene)amiloride (HMA) are potent and selective uPA inhibitors that lack the diuretic and antikaliuretic properties of the parent drug amiloride. However, the compounds display pronounced selectivity for human over mouse uPA, thus confounding interpretation of data from human xenograft mouse models of cancer. Here, computational and experimental findings reveal that residue 99 is a key contributor to the observed species selectivity, whereby enthalpically unfavorable expulsion of a water molecule by the 5-N,N-hexamethylene ring occurs when residue 99 is Tyr (as in mouse uPA). Analogue 7 lacking the 5-N,N-hexamethylene ring maintained similar water networks when bound to human and mouse uPA and displayed reduced selectivity, thus supporting this conclusion. The study will guide further optimization of dual-potent human/mouse uPA inhibitors from the amiloride class as antimetastasis drugs.

Discovery of Eucalyptin C, derived from the fruits of Eucalyptus globulus Labill., as a novel selective PI3Kγ inhibitor for immunosuppressive treatment

The fruits of Eucalyptus globulus Labill. are known to have a plenty of medicinal properties, such as anti-tumor, anti-inflammatory, and immunosuppressive activity. Our previous study found that the phloroglucinol-sesquiterpene adducts in the fruits of E. globulus were immunosuppressive active constituents, especially Eucalyptin C (EuC). Phosphoinositide 3-kinases-γ (PI3Kγ) plays a pivotal role in T cell mediated excessive immune responses. In this study, EuC was first discovered to be a novel selective PI3Kγ inhibitor with an IC₅₀ value of 0.9 μmol·L^-1 and selectivity over 40-fold towards the other PI3K isoforms. Molecular docking, molecular dynamics simulation, and cellular thermal shift assay showed that EuC bound to PI3Kγ. Furthermore, EuC suppressed the downstream of PI3Kγ to induce the apoptosis and inhibit the activation of primary spleen cells derived from allergic contact dermatitis mice. This work highlights the role of the fruits of E. globulus as a source of bioactive plant with immunosuppressive activity.

Exposing structural variations in SARS-CoV-2 evolution

The mutation of SARS-CoV-2 influences viral function as residue replacements affect both physiochemical properties and folding conformations. Although a large amount of data on SARS-CoV-2 is available, the investigation of how viral functions change in response to mutations is hampered by a lack of effective structural analysis. Here, we exploit the advances of protein structure fingerprint technology to study the folding conformational changes induced by mutations. With integration of both protein sequences and folding conformations, the structures are aligned for SARS-CoV to SARS-CoV-2, including Alpha variant (lineage B.1.1.7) and Delta variant (lineage B.1.617.2). The results showed that the virus evolution with change in mutational positions and physicochemical properties increased the affinity between spike protein and ACE2, which plays a critical role in coronavirus entry into human cells. Additionally, these structural variations impact vaccine effectiveness and drug function over the course of SARS-CoV-2 evolution. The analysis of structural variations revealed how the coronavirus has gradually evolved in both structure and function and how the SARS-CoV-2 variants have contributed to more severe acute disease worldwide.

Identification of Potent Natural Resource Small Molecule Inhibitor to Control Vibrio cholera by Targeting Its Outer Membrane Protein U: An In Silico Approach

Vibrio cholerae causes the diarrheal disease cholera which affects millions of people globally. The outer membrane protein U (OmpU) is the outer membrane protein that is most prevalent in V. cholerae and has already been recognized as a critical component of pathogenicity involved in host cell contact and as being necessary for the survival of pathogenic V. cholerae in the host body. Computational approaches were used in this study to screen a total of 37,709 natural compounds from the traditional Chinese medicine (TCM) database against the active site of OmpU. Following a sequential screening of the TCM database, we report three lead compounds-ZINC06494587, ZINC85510056, and ZINC95910434-that bind strongly to OmpU, with binding affinity values of -8.92, -8.12, and -8.78 kcal/mol, which were higher than the control ligand (-7.0 kcal/mol). To optimize the interaction, several 100 ns molecular dynamics simulations were performed, and the resulting complexes were shown to be stable in their vicinity. Additionally, these compounds were predicted to have good drug-like properties based on physicochemical properties and ADMET assessments. This study suggests that further research be conducted on these compounds to determine their potential use as cholera disease treatment.

A Comprehensive Computational Screening of Phytochemicals Derived from Saudi Medicinal Plants against Human CC Chemokine Receptor 7 to Identify Potential Anti-Cancer Therapeutics

Homeostatic trafficking of immune cells by CC chemokine receptor 7 (CCR7) keeps immune responses and tolerance in a balance. The involvement of this protein in lymph node metastasis in cancer marks CCR7 as a penitential drug target. Using the crystal structure of CCR7, herein, a comprehensive virtual screening study is presented to filter novel strong CCR7 binding phytochemicals from Saudi medicinal plants that have a higher binding affinity for the intracellular allosteric binding pocket. By doing so, three small natural molecules named as Hit-1 (1,8,10-trihydroxy-3-methoxy-6-methylanthracen-9(4H)-one), Hit-2 (4-(3,4-dimethoxybenzyl)-3-(4-hydroxy-3-methoxybenzyl)dihydrofuran-2(3H)-one), and Hit-3 (10-methyl-12,13-dihydro-[1,2]dioxolo[3,4,5-de]furo[3,2-g]isochromeno[4,3-b]chromen-8-ol) are predicted showing strong binding potential for the CC chemokine receptor 7 allosteric pocket. During molecular dynamics simulations, the compounds were observed in the formation of several chemical bonding of short bond distances. Additionally, the molecules remained in strong contact with the active pocket residues and experienced small conformation changes that seemed to be mediated by the CCR7 loops to properly engage the ligands. Two types of binding energy methods (MM/GBPBSA and WaterSwap) were additionally applied to further validate docking and simulation findings. Both analyses complement the good affinity of compounds for CCR7, the electrostatic and van der Waals energies being the most dominant in intermolecular interactions. The active pocket residue's role in compounds binding was further evaluated via alanine scanning, which highlighted their importance in natural compounds binding. Additionally, the compounds fulfilled all drug-like rules: Lipinski, Ghose, Veber, Egan, and Muegge passed many safety parameters, making them excellent anti-cancer candidates for experimental testing.

Thymoquinone and quercetin induce enhanced apoptosis in non-small cell lung cancer in combination through the Bax/Bcl2 cascade

The treatments available for non-small cell lung cancer exert various side effects in patients, and the burden of treatment cost is high. Therefore, exploring the alternative system of medicines, including therapies based on natural compounds, has become inevitable in developing anticancer therapeutics. This study used an integrated approach involving in-silico and in-vitro methods to explore natural compounds targeting Bax and Bcl2 for their apoptotic potential. Molecular docking followed by molecular dynamics (MD) simulation of thymoquinone (Tq) and quercetin (Qu) with Bax and Bcl2 were carried out to explore their interactions and stability under explicit solvent conditions. Tq and Qu showed appreciable binding affinities toward Bax (-6.2 and -7.1 kcal/mol, respectively) and Bcl2 (-5.6 and -6.4 kcal/mol, respectively) with well-organized conformational fitting compatibility. The MD simulation results revealed the development of stable complexes maintained by various noncovalent interactions that were preserved throughout the 100 ns trajectories. Further studies with these compounds were carried out using various in-vitro experimental approaches like MTT assay, apoptotic assay, and Western blot. IC₅₀ values of Tq and Qu alone in A549 cells were found to be 45.78 and 35.69 µM, while in combination, it comes down to 22.49 µM, which is quite impressive. Similarly, in apoptosis assay, a combination of Tq and Qu shows 50.9% early apoptosis compared to Tq (40.6%) and Qu (33.3%) when taken alone. These assays signify their apoptotic induction potential, whereas both compounds significantly reduce the expression of antiapoptotic protein Bcl2 and induce proapoptotic Bax, suggestive of sensitizing NSCLS cells toward apoptosis.

Designing a Multi-Epitope Vaccine against Chlamydia trachomatis by Employing Integrated Core Proteomics, Immuno-Informatics and In Silico Approaches

Chlamydia trachomatis, a Gram-negative bacterium that infects the rectum, urethra, congenital sites, and columnar epithelium of the cervix. It is a major cause of preventable blindness, ectopic pregnancy, and bacterial sexually transmitted infections worldwide. There is currently no licensed multi-epitope vaccination available for this pathogen. This study used core proteomics, immuno-informatics, and subtractive proteomics approaches to identify the best antigenic candidates for the development of a multi-epitope-based vaccine (MEBV). These approaches resulted in six vaccine candidates: Type III secretion system translocon subunit CopD2, SctW family type III secretion system gatekeeper subunit CopN, SycD/LcrH family type III secretion system chaperone Scc2, CT847 family type III secretion system effector, hypothetical protein CTDEC_0668, and CHLPN 76kDa-like protein. A variety of immuno-informatics tools were used to predict B and T cell epitopes from vaccine candidate proteins. An in silico vaccine was developed using carefully selected epitopes (11 CTL, 2 HTL & 10 LBL) and then docked with the MHC molecules (MHC I & MHC II) and human TLR4. The vaccine was coupled with Cholera toxin subunit B (CTB) adjuvant to boost the immune response. Molecular dynamics (MD) simulations, molecular docking, and MMGBSA analysis were carried out to analyze the molecular interactions and binding affinity of MEBV with TLR4 and MHC molecules. To achieve the highest level of vaccine protein expression, the MEBV was cloned and reverse-translated in Escherichia coli. The highest level of expression was achieved, and a CAI score of 0.97 was reported. Further experimental validation of the MEBV is required to prove its efficacy. The vaccine developed will be useful in preventing infections caused by C. trachomatis.

Direct Synthesis of Diamides from Dicarboxylic Acids with Amines Using Nb2O5 as a Lewis Acid Catalyst and Molecular Docking Studies as Anticancer Agents

Several Lewis and Bronsted acid catalysts were tested for the synthesis of some targeted diamides with anticancer activity from dicarboxylic acids and amines under the same reaction condition. Among those catalysts, Nb₂O₅ showed the highest catalytic activity to the corresponding diamides. Nb₂O₅ shows water- and base-tolerant properties for which it gives the highest yield of the synthesized products. Here, we present a novel and sustainable method for the direct synthesis of diamides with anticancer activity using a reusable heterogeneous catalyst Nb₂O₅. A molecular docking study was performed for all of the synthesized compounds with various therapeutical targets of cancer and found that the human epidermal growth factor receptor (HER2) has shown a significant dock score for our synthesized products. After obtaining the best pose from molecular docking, the complex is used for molecular dynamics study by running simulations for 10 ns. The root-mean-square deviations (RMSDs) of α carbon atoms of all systems are analyzed to detect their stability. This method is effective for the direct synthesis of diamides as anticancer agents from dicarboxylic acids and amines using Nb₂O₅ as a base-tolerant heterogeneous catalyst.

Pathogenic genetic variants from highly connected cancer susceptibility genes confer the loss of structural stability

Genetic polymorphisms in DNA damage repair and tumor suppressor genes have been associated with increasing the risk of several types of cancer. Analyses of putative functional single nucleotide polymorphisms (SNP) in such genes can greatly improve human health by guiding choice of therapeutics. In this study, we selected nine genes responsible for various cancer types for gene enrichment analysis and found that BRCA1, ATM, and TP53 were more enriched in connectivity. Therefore, we used different computational algorithms to classify the nonsynonymous SNPs which are deleterious to the structure and/or function of these three proteins. The present study showed that the major pathogenic variants (V1687G and V1736G of BRCA1, I2865T and V2906A of ATM, V216G and L194H of TP53) might have a greater impact on the destabilization of the proteins. To stabilize the high-risk SNPs, we performed mutation site-specific molecular docking analysis and validated using molecular dynamics (MD) simulation and molecular mechanics/Poisson Boltzmann surface area (MM/PBSA) studies. Additionally, SNPs of untranslated regions of these genes affecting miRNA binding were characterized. Hence, this study will assist in developing precision medicines for cancer types related to these polymorphisms.

Perturbation Free-Energy Toolkit: An Automated Alchemical Topology Builder

Free-energy calculations play an important role in the application of computational chemistry to a range of fields, including protein biochemistry, rational drug design, or materials science. Importantly, the free-energy difference is directly related to experimentally measurable quantities such as partition and adsorption coefficients, water activity, and binding affinities. Among several techniques aimed at predicting free-energy differences, perturbation approaches, involving the alchemical transformation of one molecule into another through intermediate states, stand out as rigorous methods based on statistical mechanics. However, despite the importance of free-energy calculations, the applicability of the perturbation approaches is still largely impeded by a number of challenges, including the definition of the perturbation path, i.e., alchemical changes leading to the transformation of one molecule to the other. To address this, an automatic perturbation topology builder based on a graph-matching algorithm is developed, which can identify the maximum common substructure (MCS) of two or multiple molecules and provide the perturbation topologies suitable for free-energy calculations using the GROMOS and the GROMACS simulation packages. Various MCS search options are presented leading to alternative definitions of the perturbation pathway. Moreover, perturbation topologies generated using the default multistate MCS search are used to calculate the changes in free energy between lysine and its two post-translational modifications, 3-methyllysine and acetyllysine. The pairwise free-energy calculations performed on this test system led to a cycle closure of 0.5 ± 0.3 and 0.2 ± 0.2 kJ mol-1, with GROMOS and GROMACS simulation packages, respectively. The same relative free energies between the three states are obtained by employing the enveloping distribution sampling (EDS) approach when compared to the pairwise perturbations. Importantly, this toolkit is made available online as an open-source Python package (https://github.com/drazen-petrov/SMArt).

An In-Silico Approach to Evaluate the Inhibitory Potency of Selected Hydroxamic Acid Derivatives on Zinc-Dependent Histone Deacetylase Enzyme

Histone deacetylase (HDAC) enzymes modify the histone by removing the acetyl group from the lysine residues, known as histone deacetylation. HDACs have been involved in altering gene expressions, resulting in cancer cells in the body. This study focuses on HDAC inhibitors’ impact on histone deacetylase-like protein (HDLP) stability through computational techniques. Molecular dynamics (MD) analyses were used to examine the atomic-level description of drug binding sites and how the HDAC inhibitors change the HDLP enzyme environment. In this study, two hydroxamic acid-derived inhibitors, such as [Formula: see text]-Carboxycinnamic acid bis-hydroxamide (CBHA) and scriptaid (GCK1026), were selected to examine the inhibition ability in terms with suberanilohydroxamic acid (SAHA) as a reference drug. The crystal structure of the HDLP was downloaded from the Protein Data Bank. The structures of inhibitors were optimized using the G09W package. Docking studies were done by AutoDock-Vina, and the resultant complex was used to initiate MD studies. The trajectories obtained from MD simulation were used to perform the structural analysis. Root-mean-square deviation (RMSD), radius of gyration, hydrogen bond, binding free energy and interaction energy studies revealed that the stability of HDLP-SAHA and HDLP-CBHA is higher than the free HDLP enzyme. The HDLP-CBHA complex shows an increased number of hydrogen bonds (5), high MM-PBSA binding free energy ([Formula: see text][Formula: see text]kJ/mol), high interaction energy ([Formula: see text][Formula: see text]kJ/mol), and an increased number of alpha-helical amino acids (130) compared with HDLP-SAHA. It concluded that the CBHA has the relatively same potential as SAHA to inhibit the HDLP. Consequently, the use of CBHA in clinical application is recommended through this in-silico method.

Molecular Mechanism Exploration of Potent Fluorinated PI3K Inhibitors with a Triazine Scaffold: Unveiling the Unusual Synergistic Effect of Pyridine-to-Pyrimidine Ring Interconversion and CF3 Defluorination

The phosphatidylinostitol-3-kinase (PI3K)/AKT/mammalian target of rapamycin signaling pathway is a vital regulator of cell proliferation, growth, and survival, which is frequently overactivated in many human cancers. To this effect, PI3K, which is an important mediator of this pathway, has been pinpointed as a crucial target in cancer therapy and hence the importance of PI3K inhibitors. It was recently reported that defluorination and pyridine-to-pyrimidine ring interconversion increase the potency of specific small-molecule inhibitors of PI3K. Compound 4, an inhibitor with the difluorinated pyrimidine motif, was found to be eight times more potent against PI3K than compound 1, an inhibitor with the trifluorinated pyridine motif. This observation presents the need to rationally resolve the differential inhibitory mechanisms exhibited by both compounds. In this present work, we employed multiple computational approaches to investigate and distinguish the binding modes of 1 and 4 in addition to the effects they mediate on the secondary structure of PI3K. Likewise, we evaluated two other derivatives, compounds 2 with the difluorinated pyridine motif and 3 with the trifluorinated pyrimidine motif, to investigate the cooperativity effect between the defluorination of CF₃ and pyridine-to-pyrimidine ring interconversion. Findings revealed that PI3K, upon interaction with 4, exhibited a series of structural changes that favored the binding of the inhibitor at the active-site region. Furthermore, a positive (synergistic) cooperativity effect was observed between CF₃ defluorination and pyridine-to-pyrimidine ring interconversion. Moreover, there was a good correlation between the binding free energy estimated and the biological activity reported experimentally. Energy decomposition analysis revealed that the major contributing force to binding affinity variations between 1 and 4 is the electrostatic energy. Per-residue energy-based hierarchical clustering analysis further identified four hot-spot residues ASP841, TYR867, ASP964, and LYS833 and four warm-spot residues ASP836, SER806, ASP837, and LYS808, which essentially mediate the optimal and higher-affinity binding of compound 4 to PI3K relative to 1. This study therefore provides rational insights into the mechanisms by which 4 exhibited superior PI3K-inhibitory activities over 1, which is vital for future structure-based drug discovery efforts in PI3K targeting.

Exposing Structural Variations in SARS-CoV-2 Evolution

The mutation of SARS-CoV-2 influences viral function as residue replacements affect both physiochemical properties and folding conformations. Although a large amount of data on SARS-CoV-2 is available, the investigation of how viral functions change in response to mutations is hampered by a lack of effective structural analysis. Here, we exploit advances in protein structure fingerprint technology to study the folding conformational changes induced by mutations. With the integration of both protein sequences and folding conformations and alignments of SARS-CoV to SARS-CoV-2, the UK variant and India variant, we found that structural variations in the spike protein at the binding interface interacting with ACE2 play a critical role in coronavirus entry into human cells. Additionally, the structural variations impact vaccine effectiveness and drug function over the course of SARS-CoV-2 evolution. The analysis of structural variations revealed how the coronavirus has gradually evolved in both structure and function and how the SARS-CoV-2 variants have contributed to more severe acute disease worldwide.

Phytoconstituents and Medicinal Plants for Anticancer Drug Discovery: Computational Identification of Potent Inhibitors of PIM1 Kinase

Natural products, medicinal plants, and phytoconstituents serve as important sources and accelerators for anticancer drug discovery, especially when they are combined with virtual screening and molecular simulations against molecular drug targets. Proto-oncogene serine/threonine-protein kinase Pim1 (PIM1) is involved in cell survival and proliferation, with great relevance for oncogenesis. PIM1 plays a major role in the progression of various common complex human cancers, including prostate cancer, acute myeloid leukemia, and other hematopoietic malignancies. The overexpression of PIM1 leads to cancer progression, and thus it is considered as a potential target for drug design and development purposes. Here, we report original in silico findings by employing structure-based virtual screening to discover potential phytoconstituents from the medicinal plants-based compounds, which could inhibit the PIM1 activity, using the IMPPAT (a curated database of Indian Medicinal Plants, Phytochemistry And Therapeutics) database. The initial hits were selected based on their binding affinity toward PIM1 calculated through the molecular docking approach. Subsequently, interaction analyses and molecular dynamics (MD) simulation for 100 ns was carried out to study the conformational dynamics and complex stability of PIM1 with the identified compounds. Importantly, we found that PIM1 forms stable protein-ligand complexes with the phytoconstituents Dehydrotectol and Nordracorubin in particular. Our findings suggest that identified phytoconstituents Dehydrotectol and Nordracorubin bind to PIM1 in ATP-competitive binding mode. These findings and the compounds reported herein warrant further exploration as promising scaffolds for anticancer drug design, discovery, and development.

Discovery of Novel Apigenin-Piperazine Hybrids as Potent and Selective Poly (ADP-Ribose) Polymerase-1 (PARP-1) Inhibitors for the Treatment of Cancer

Poly (ADP-ribose) polymerase-1 (PARP-1) is a potential target for the discovery of chemosensitizers and anticancer drugs. Amentoflavone (AMF) is reported to be a selective PARP-1 inhibitor. Here, structural modifications and trimming of AMF have led to a series of AMF derivatives (9a-h) and apigenin-piperazine/piperidine hybrids (14a-p, 15a-p, 17a-h, and 19a-f), respectively. Among these compounds, 15l exhibited a potent PARP-1 inhibitory effect (IC50 = 14.7 nM) and possessed high selectivity to PARP-1 over PARP-2 (61.2-fold). Molecular dynamics simulation and the cellular thermal shift assay revealed that 15l directly bound to the PARP-1 structure. In in vitro and in vivo studies, 15l showed a potent chemotherapy sensitizing effect against A549 cells and a selective cytotoxic effect toward SK-OV-3 cells through PARP-1 inhibition. 15l·2HCl also displayed good ADME characteristics, pharmacokinetic parameters, and a desirable safety margin. These findings demonstrated that 15l·2HCl may serve as a lead compound for chemosensitizers and the (BRCA-1)-deficient cancer therapy.

Therapeutic development by repurposing drugs targeting SARS-CoV-2 spike protein interactions by simulation studies

The human-to-human transmitted respiratory illness in COVID-19 affected by the pathogenic Severe Acute Respiratory Syndrome Corona Virus 2 (SARS-CoV-2), which appeared in the last of December 2019 in Wuhan, China, and rapidly spread in many countries. Thereon, based on the urgent need for therapeutic molecules, we conducted in silico based docking and simulation molecular interaction studies on repurposing drugs, targeting SARS-CoV-2 spike protein. Further, the best binding energy of doxorubicin interacting with virus spike protein (PDB: 6VYB) was observed to be -6.38 kcal/mol and it was followed by exemestane and gatifloxacin. The molecular simulation dynamics analysis of doxorubicin, Reference Mean Square Deviation (RMSD), Root Mean Square fluctuation (RMSF), Radius of Gyration (Rg), and formation of hydrogen bonds plot interpretation suggested, a significant deviation and fluctuation of Doxorubicin-Spike RBD complex during the whole simulation period. The Rg analysis has stated that the Doxorubicin-Spike RBD complex was stable during 15,000-35,000 ps MDS. The results have suggested that doxorubicin could inhibit the virus spike protein and prevent the access of the SARS-CoV-2 to the host cell. Thus, in-vitro/in-vivo research on these drugs could be advantageous to evaluate significant molecules that control the COVID-19 disease.

Discovery of Natural Compounds as Potential Inhibitors of Human Carbonic Anhydrase II: An Integrated Virtual Screening, Docking, and Molecular Dynamics Simulation Study

Carbonic anhydrase II (CAII) is one of the zinc metalloenzymes that catalyze the reversible hydration of carbon dioxide, leading to the formation of bicarbonate and proton. CAII plays a significant role in health and disease. For example, CAII helps to maintain eye pressure while regulating the pH of the tumor microenvironment, and by extension, contributing to cancer progression. Owing to its remarkable role in cancer, visual health, and other human diseases, CAII can serve as an attractive therapeutic target. We report an original study based on high-throughput virtual screening of natural compounds from the ZINC database in search of potential inhibitors of CAII. We selected the hits based on the physicochemical, absorption, distribution, metabolism, excretion, and toxicity (ADMET) properties, pan-assay interference compound (PAINS) patterns, and interaction analysis. Importantly, two natural compounds were identified, ZINC08918123 and ZINC00952700, bearing considerable affinity and specific interactions to the residues of the CAII-binding pocket with well-organized conformational fitting compatibility. We investigated the conformational dynamics of CAII in complex with the identified compounds through molecular dynamics simulation, which revealed the formation of a stable complex preserved throughout the 100 ns trajectories. The stability of the protein/ligand complexes is maintained by significant numbers of noncovalent interactions throughout the simulations. In conclusion, natural compounds identified in the present study specifically and computer-assisted drug design broadly offer a reliable resource and strategy to discover potential promising therapeutic inhibitors of CAII to cure various cancers and glaucoma after further experimental validation and clinical studies.

Remdesivir Strongly Binds to RNA-Dependent RNA Polymerase, Membrane Protein, and Main Protease of SARS-CoV-2: Indication From Molecular Modeling and Simulations

Development of new drugs is a time-taking and expensive process. Comprehensive efforts are being made globally toward the search of therapeutics against SARS-CoV-2. Several drugs such as remdesivir, favipiravir, ritonavir, and lopinavir have been included in the treatment regimen and shown effective results in several cases. Among the existing broad-spectrum antiviral drugs, remdesivir is found to be more effective against SARS-CoV-2. Remdesivir has broad-spectrum antiviral action against many single-stranded RNA viruses including pathogenic SARS-CoV and Middle East respiratory syndrome coronavirus (MERS-CoV). In this study, we proposed that remdesivir strongly binds to membrane protein (Mprotein), RNA-dependent RNA polymerase (RDRP), and main protease (Mprotease) of SARS-CoV-2. It might show antiviral activity by inhibiting more than one target. It has been found that remdesivir binds to Mprotease, Mprotein, and RDRP with -7.8, -7.4, and -7.1 kcal/mol, respectively. The structure dynamics study suggested that binding of remdesivir leads to unfolding of RDRP. It has been found that strong binding of remdesivir to Mprotein leads to decrease in structural deviations and gyrations. Additionally, the average solvent-accessible surface area of Mprotein decreases from 127.17 to 112.12 nm², respectively. Furthermore, the eigenvalues and the trace of the covariance matrix were found to be low in case of Mprotease-remdesivir, Mprotein-remdesivir, and RDRP-remdesivir. Binding of remdesivir to Mprotease, Mprotein, and RDRP reduces the average motions in protein due to its strong binding. The MMPBSA calculations also suggested that remdesivir has strong binding affinity with Mprotein, Mprotease, and RDRP. The detailed analysis suggested that remdesivir has more than one target of SARS-CoV-2.

Integrated Core Proteomics, Subtractive Proteomics, and Immunoinformatics Investigation to Unveil a Potential Multi-Epitope Vaccine against Schistosomiasis

Schistosomiasis is a parasitic infection that causes considerable morbidity and mortality in the world. Infections of parasitic blood flukes, known as schistosomes, cause the disease. No vaccine is available yet and thus there is a need to design an effective vaccine against schistosomiasis. Schistosoma japonicum, Schistosoma mansoni, and Schistosoma haematobium are the main pathogenic species that infect humans. In this research, core proteomics was combined with a subtractive proteomics pipeline to identify suitable antigenic proteins for the construction of a multi-epitope vaccine (MEV) against human-infecting Schistosoma species. The pipeline revealed two antigenic proteins-calcium binding and mycosubtilin synthase subunit C-as promising vaccine targets. T and B cell epitopes from the targeted proteins were predicted using multiple bioinformatics and immunoinformatics databases. Seven cytotoxic T cell lymphocytes (CTL), three helper T cell lymphocytes (HTL), and four linear B cell lymphocytes (LBL) epitopes were fused with a suitable adjuvant and linkers to design a 217 amino-acid-long MEV. The vaccine was coupled with a TLR-4 agonist (RS-09; Sequence: APPHALS) adjuvant to enhance the immune responses. The designed MEV was stable, highly antigenic, and non-allergenic to human use. Molecular docking, molecular dynamics (MD) simulations, and molecular mechanics/generalized Born surface area (MMGBSA) analysis were performed to study the binding affinity and molecular interactions of the MEV with human immune receptors (TLR2 and TLR4) and MHC molecules (MHC I and MHC II). The MEV expression capability was tested in an Escherichia coli (strain-K12) plasmid vector pET-28a(+). Findings of these computer assays proved the MEV as highly promising in establishing protective immunity against the pathogens; nevertheless, additional validation by in vivo and in vitro experiments is required to discuss its real immune-protective efficacy.

Insilico drug repurposing using FDA approved drugs against Membrane protein of SARS-CoV-2

The novel coronavirus (SARS-CoV-2) outbreak has started taking away the millions of lives worldwide. Identification of known and approved drugs against novel coronavirus disease (COVID-19) seems to be an urgent need for the repurposing of the existing drugs. So, here we examined a safe strategy of using approved drugs of SuperDRUG2 database against modeled membrane protein (M-protein) of SARS-CoV-2 which is essential for virus assembly by using molecular docking-based virtual screening. A total of 3639 drugs from SuperDRUG2 database and additionally 14 potential drugs reported against COVID-19 proteins were selected. Molecular docking analyses revealed that nine drugs can bind the active site of M-protein with desirable molecular interactions. We therefore applied molecular dynamics simulations and binding free energy calculation using MM-PBSA to analyze the stability of the compounds. The complexes of M-protein with the selected drugs were simulated for 50 ns and ranked according to their binding free energies. The binding mode of the drugs with M-protein was analyzed and it was observed that Colchicine, Remdesivir, Bafilomycin A1 from COVID-19 suggested drugs and Temozolomide from SuperDRUG2 database displayed desirable molecular interactions and higher binding affinity towards M-protein. Interestingly, Colchicine was found as the top most binder among tested drugs against M-protein. We therefore additionally identified four Colchicine derivatives which can bind efficiently with M-protein and have better pharmacokinetic properties. We recommend that these drugs can be tested further through in vitro studies against SARS-CoV-2 M-protein.

Flexible Loop in Carbohydrate-Binding Module 48 Allosterically Modulates Substrate Binding of the 1,4-α-Glucan Branching Enzyme

The 1,4-α-glucan branching enzyme (GBE, EC 2.4.1.18) catalyzes the formation of α-1,6 branching points in starch and plays a key role in synthesis. To obtain mechanistic insights into the catalytic action of the enzyme, we first determined the crystal structure of GBE from Rhodothermus obamensis STB05 (RoGBE) to a resolution of 2.39 Å (PDB ID: 6JOY). The structure consists of three domains: domain A, domain C, and the carbohydrate-binding module 48 (CBM48). An engineered truncated mutant lacking the CBM48 domain (ΔCBM48) showed significantly reduced ligand binding affinity and enzyme activity. Comparison of the structures of RoGBE with other GBEs showed that CBM48 of RoGBE had a longer flexible loop. Truncation of the flexible loops resulted in reduced binding affinity and activity, thereby substantiating the importance of the optimum loop structure for catalysis. In essence, our study shows that CBM48, especially the flexible loop, plays an important role in substrate binding and enzymatic activity of RoGBE. Further, based on the structural analysis, kinetics, and activity assays on wild type and mutants, as well as homology modeling, we proposed a mechanistic model (called the "lid model") to illustrate how the flexible loop triggers substrate binding, ultimately leading to catalysis.

Synthesis, SAR study, and bioactivity evaluation of a series of Quinoline-Indole-Schiff base derivatives: Compound 10E as a new Nur77 exporter and autophagic death inducer

We previously reported 5-((8-methoxy-2-methylquinolin-4-yl)amino)-1H-indole- 2-carbohydrazide derivatives as new Nur77 modulators. In this study, we explored whether the 8-methoxy-2-methylquinoline moiety and bicyclic aromatic rings at the N'-methylene position were critical for their antitumor activity against hepatocellular carcinoma (HCC). For this purpose, a small library of 5-substituted 1H-indole-2-carbohydrazide derivatives was designed and synthesized. We found that the 8-methoxy-2-methylquinoline moiety was a fundamental structure for its biological function, while the introduction of the bicyclic aromatic ring into the N'-methylene greatly improved its anti-tumor effect. We found that the representative compound 10E had a high affinity to Nur77. The K_D values were in the low micromolar (2.25-4.10 μM), which were coincident with its IC₅₀ values against the tumor cell lines (IC₅₀ < 3.78 μM). Compound 10E could induce autophagic cell death of liver cancer cells by targeting Nur77 to mitochondria while knocking down Nur77 greatly impaired anti-tumor effect. These findings provide an insight into the structure-activity relation of Quinoline-Indole-Schiff base derivatives and further demonstrate that antitumor agents targeting Nur77 may be considered as a promising strategy for HCC therapy.

Optimization of N-Phenylpropenoyl-l-amino Acids as Potent and Selective Inducible Nitric Oxide Synthase Inhibitors for Parkinson's Disease

N-Phenylpropenoyl-l-amino acids (NPAs) are inducible nitric oxide synthase (iNOS) inhibitors possessing preventive effects for Parkinson's disease (PD). Here, structural modifications for improving the iNOS inhibitory activity and blood-brain barrier (BBB) permeability of NPAs were conducted, leading to 20 optimized NPA derivatives (1-20). Compound 18, with the most potent activity (IC₅₀ = 74 nM), high BBB permeability (P_e = 19.1 × 10^-6 cm/s), and high selectivity over other NOS isoforms, was selected as the lead compound. Further studies demonstrated that 18 directly binds to iNOS. In the 1-methyl-4-phenyl-1,2,3,6-tetrahydropyridine (MPTP)-induced acute PD model, the oral administration of 18 (1 and 2 mg/kg) exerted preventive effects by alleviating the loss of dopaminergic (DAergic) neurons. Notably, in the MPTP-/probenecid-induced chronic PD model, the same dose of 18 also displayed a therapeutic effect by repairing the damaged DAergic neurons. Finally, good pharmacokinetic properties and low toxicity made 18 a promising candidate for the treatment of PD.

Biomolecular response to hour-long ultralow field microwave radiation: An effective coarse-grained model simulation

Various electronic devices, which we commonly use, radiate microwaves. Such external perturbation influences the functionality of biomolecules. In an ultralow field, the cumulative response of a molecule is expected only over a time scale of hours. To study the structural dynamics of biomolecules over hours, we adopt a simple methodology for constructing the coarse-grained structure of the protein molecule and solve the Langevin equation under different working potentials. In this approach, each amino acid residue of a biomolecule is mapped onto a number of beads, a few for the backbone, and few for the side chain, depending on the complexity of its chemical structure. We choose the force field in such a way that the dynamics of the protein molecule in the presence of ultralow radiation field of microvolt/nm could be followed over the time frame of 2 h. We apply the model to describe a biomolecule, hen egg white lysozyme, and simulate its structural evolution under ultralow strength electromagnetic radiation. The simulation revealed the finer structural details, like the extent of exposure of bioactive residues and the state of the secondary structures of the molecule, further confirmed from spectroscopic measurements [details are available in Phys. Rev. E 97, 052416 (2018)10.1103/PhysRevE.97.052416 and briefly described here]. Though tested for a specific system, the model is quite general. We believe that it harnesses the potential in studying the structural dynamics of any biopolymer under external perturbation over an extended time scale.

Impact of amino acid substitutions on the behavior of a photoactivatable near infrared fluorescent protein PAiRFP1

A photoactivatable near-infrared fluorescent protein (NIR-FP) PAiRFP1 has been developed by 15 amino acid substitutions in its nonfluorescent template Agp2. In our previous communication, we investigated the role of three amino acids in PHY domain distal from BV molecule. The impact of the twelve amino acids in GAF domain, especially five residues near BV-binding pocket is unclear. In this paper, PCR based reverse mutagenesis, spectroscopic methods, molecular modelling and simulations have been employed to explore the roles of these substitutions during the molecular evolution of PAiRFP1. It was found that the residue L163 is important for protein folding in PAiRFP1. The residues F244 and C280 exerted remarkable effects on molar extinction coefficient, NIR fluorescence quantum yield, molecular brightness, fluorescence fold, and dark recovery rate. The residues F244 and V276 modulate the maximum absorption and emission peak position. The reverse mutant L168M exhibited a higher fluorescence fold than PAiRFP1. Additionally, the reverse mutants V203A, V294E, S218G and D127G possessed better spectral properties than PAiRFP1. This study is important for the rational design of a better BphP-based photoactivatable NIR-FPs.

Binding interaction of glyphosate with glyphosate oxidoreductase and C-P lyase: Molecular docking and molecular dynamics simulation studies

Widespread application of glyphosate poses a threat to living organisms. Microbial strains are able to degrade glyphosate via contrasting metabolic pathways with the help of enzymes. Glyphosate oxidoreductase (GOX) and C-P lyase are the key enzymes for the biodegradation of glyphosate and its intermediate metabolite aminomethylphosphonic acid (AMPA) in microbes. The microbial degradation of glyphosate has been reported, but the underlying molecular mechanism is still unclear. Therefore, in this study, the interaction mechanism of GOX and C-P lyase with glyphosate and AMPA were investigated by using molecular docking and molecular dynamics (MD) simulations. The results indicate that glyphosate contacts with the active site of GOX and C-P lyase by hydrogen bonds as well as hydrophobic and van der Waals interactions in aqueous solution to maintain its stability. The presence of glyphosate and AMPA in the active site significantly changes the conformation of GOX and C-P lyase. The results of the MD simulations confirm that GOX and C-P lyase complexes are stable during the catalytic reaction. This study offers a molecular level of understanding of the expression and function of GOX and C-P lyase for the bioremediation of glyphosate.

Design of Inhibitors for Glyceraldehyde-3-phosphate Dehydrogenase (GAPDH) Enzyme of <i>Leishmania mexicana</i>

BACKGROUND

Leishmaniosis is a neglected tropical disease and glyceraldehyde 3- phosphate dehydrogenase (GAPDH) is a key enzyme in the design of new drugs to fight this disease.

OBJECTIVE

The present study aimed to evaluate potential inhibitors of GAPDH enzyme found in Leishmania mexicana (L. mexicana).

METHODS

A search for novel antileishmanial molecules was carried out based on similarities from the pharmacophoric point of view related to the binding site of the crystallographic enzyme using the ZINCPharmer server. The molecules selected in this screening were subjected to molecular docking and molecular dynamics simulations.

RESULTS

Consensual analysis of the docking energy values was performed, resulting in the selection of ten compounds. These ligand-receptor complexes were visually inspected in order to analyze the main interactions and subjected to toxicophoric evaluation, culminating in the selection of three compounds, which were subsequently submitted to molecular dynamics simulations. The docking results showed that the selected compounds interacted with GAPDH from L. mexicana, especially by hydrogen bonds with Cys166, Arg249, His194, Thr167, and Thr226. From the results obtained from molecular dynamics, it was observed that one of the loop regions, corresponding to the residues 195-222, can be related to the fitting of the substrate at the binding site, assisting in the positioning and the molecular recognition via residues responsible for the catalytic activity.

CONCLUSION

The use of molecular modeling techniques enabled the identification of promising compounds as inhibitors of the GAPDH enzyme from L. mexicana, and the results obtained here can serve as a starting point to design new and more effective compounds than those currently available.

Methods for Identifying Microbial Natural Product Compounds that Target Kinetoplastid RNA Structural Motifs by Homology and De Novo Modeled 18S rRNA

The development of novel anti-infectives against Kinetoplastids pathogens targeting proteins is a big problem occasioned by the antigenic variation in these parasites. This is also a global concern due to the zoonosis of these parasites, as they infect both humans and animals. Therefore, we need not only to create novel antibiotics, but also to speed up the development pipeline for these antibiotics. This may be achieved by using novel drug targets for Kinetoplastids drug discovery. In this study, we focused our attention on motifs of rRNA molecules that have been created using homology modeling. The RNA is the most ambiguous biopolymer in the kinetoplatid, which carries many different functions. For instance, tRNAs, rRNAs, and mRNAs are essential for gene expression both in the pro-and eukaryotes. However, all these types of RNAs have sequences with unique 3D structures that are specific for kinetoplastids only and can be used to shut down essential biochemical processes in kinetoplastids only. All these features make RNA very potent targets for antibacterial drug development. Here, we combine in silico methods combined with both computational biology and structure prediction tools to address our hypothesis. In this study, we outline a systematic approach for identifying kinetoplastid rRNA-ligand interactions and, more specifically, techniques that can be used to identify small molecules that target particular RNA. The high-resolution optimized model structures of these kineoplastids were generated using RNA 123, where all the stereochemical conflicts were solved and energies minimized to attain the best biological qualities. The high-resolution optimized model's structures of these kinetoplastids were generated using RNA 123 where all the stereochemical conflicts were solved and energies minimized to attain the best biological qualities. These models were further analyzed to give their docking assessment reliability. Docking strategies, virtual screening, and fishing approaches successfully recognized novel and myriad macromolecular targets for the myxobacterial natural products with high binding affinities to exploit the unmet therapeutic needs. We demonstrate a sensible exploitation of virtual screening strategies to 18S rRNA using natural products interfaced with classical maximization of their efficacy in phamacognosy strategies that are well established. Integration of these virtual screening strategies in natural products chemistry and biochemistry research will spur the development of potential interventions to these tropical neglected diseases.

Designing multi-epitope vaccine against Staphylococcus aureus by employing subtractive proteomics, reverse vaccinology and immuno-informatics approaches

Staphylococcus aureus is a deadly human bacterial pathogen that causes a wide variety of clinical manifestations. Invasive S. aureus infections in hospitals and the community are one of the main causes of mortality and morbidity, as virulent and multi-drug-resistant strains have evolved. There is an unmet and urgent clinical need for immune-based non-antibiotic approaches to treat these infections as the growing antibiotic resistance poses a significant public health danger. Subtractive proteomics assisted reverse vaccinology-based immunoinformatics pipeline was used in this study to target the suitable antigenic proteins for the development of multi-epitope vaccine (MEV). Three essential virulent and antigenic proteins were identified including Glycosyltransferase, Elastin Binding Protein, and Staphylococcal secretory antigen. A variety of immunoinformatics tools have been used to forecast T-cell and B-cell epitopes from target proteins. Seven CTL, five HTL, and eight LBL epitopes, connected through suitable linkers and adjuvant, were employed to design 444 amino acids long MEV construct. The vaccine was paired with the TLR4 agonist 50S ribosomal protein L7/L12 adjuvant to enhance the immune response towards the vaccine. The predicted MEV structure was assessed to be highly antigenic, non-toxic, non-allergenic, flexible, stable, and soluble. Molecular docking simulation of the MEV with the human TLR4 (toll-like receptor 4) and major histocompatibility complex molecules (MHCI and MHCII) was performed to validate the interactions with the receptors. Molecular dynamics (MD) simulation and MMGBSA binding free energy analyses were carried out for the stability evaluation and binding of the MEV docked complexes with TLR4, MHCI and MHCII. To achieve maximal vaccine protein expression with optimal post-translational modifications, MEV was reverse translated, its mRNA structure was analyzed, and finally in silico cloning was performed into E. coli expression host. These rigorous computational analyses supported the effectivity of proposed MEV in protection against infections associated with S. aureus. However, further experimental validations are required to fully evaluate the potential of proposed vaccine candidate.

Circadian clock modulating small molecules repurposing as inhibitors of SARS-CoV-2 Mpro for pharmacological interventions in COVID-19 pandemic

The COVID-19 pandemic caused by SARS-CoV-2 is a global health emergency warranting the development of targeted treatment. The main protease M^pro is considered as a key drug target in coronavirus infections because of its vital role in the proteolytic processing of two essential polyproteins required for the replication and transcription of viral RNA. Targeting and inhibiting the M^pro activity represents a valid approach to prevent the SARS-CoV-2 replication and spread. Based on the structure-assisted drug designing, here we report a circadian clock-modulating small molecule “SRT2183” as a potent inhibitor of M^pro to block the replication of SARS-CoV-2. The findings are expected to pave the way for the development of therapeutics for COVID-19.

Characterization of the binding between anti-tumor drug 5-fluorouracil and human alpha-2-macroglobulin: spectroscopic and molecular docking analyses

5-Fluorouracil (5-FU) is a well-recognized anticancer drug used in the treatment of tumors of head, neck and breast. Drug pharmacokinetics is affected upon binding with protein, thus, making drug-protein interactions imperative to study. Present work investigates the interaction between 5-FU and human major antiproteinase-alpha-2-macroglobulin (α₂M) by multi-spectroscopic, calorimetric and molecular docking techniques. UV/Visible absorption, intrinsic fluorescence and circular dichroism (CD) spectroscopic methods have been employed to unveil the mode and mechanism of 5-FU-α₂M interaction. Synchronous fluorescence showed alteration in the microenvironment of tryptophan and tyrosine residues of protein. Far UV-CD spectra suggest slight alterations in the secondary structure of α₂M by 5-FU. Thermodynamic parameters determined by fluorescence quenching experiments and isothermal titration calorimetry (ITC) suggested the involvement of hydrogen bonds and hydrophobic interactions. Moreover, ITC corroborate the spontaneous and exothermic nature of the interaction process. Molecular docking illustrates that 5-FU binds with moderate affinity and Asp953, Tyr1264, Lys1236, Thr1232, Tyr1323 and Leu951 were the main residues involved. Molecular dynamics simulation studies suggested that 5-FU was stabilizing the α₂M structure and forming a stable complex. It was concluded that 5-FU lower the antiproteolytic activity of α₂M significantly and causes disruption in the native structure and conformation of α₂M.Communicated by Ramaswamy H. Sarma.

Identification of natural compounds as potent inhibitors of SARS-CoV-2 main protease using combined docking and molecular dynamics simulations

Coronavirus disease 2019 (COVID-19) has emerged from China and globally affected the entire population through the human-to-human transmission of a newly emerged virus called severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). The genome of SARS-CoV-2 encodes several proteins that are essential for multiplication and pathogenesis. The main protease (Mpro or 3CLpro) of SARS-CoV-2 plays a central role in its pathogenesis and thus is considered as an attractive drug target for the drug design and development of small-molecule inhibitors. We have employed an extensive structure-based high-throughput virtual screening to discover potential natural compounds from the ZINC database which could inhibit the Mpro of SARS-CoV-2. Initially, the hits were selected on the basis of their physicochemical and drug-like properties. Subsequently, the PAINS filter, estimation of binding affinities using molecular docking, and interaction analyses were performed to find safe and potential inhibitors of SARS-CoV-2 Mpro. We have identified ZINC02123811 (1-(3-(2,5,9-trimethyl-7-oxo-3-phenyl-7H-furo[3,2-g]chromen-6-yl)propanoyl)piperidine-4-carboxamide), a natural compound bearing appreciable affinity, efficiency, and specificity towards the binding pocket of SARS-CoV-2 Mpro. The identified compound showed a set of drug-like properties and preferentially binds to the active site of SARS-CoV-2 Mpro. All-atom molecular dynamics (MD) simulations were performed to evaluate the conformational dynamics, stability and interaction mechanism of Mpro with ZINC02123811. MD simulation results indicated that Mpro with ZINC02123811 forms a stable complex throughout the trajectory of 100 ns. These findings suggest that ZINC02123811 may be further exploited as a promising scaffold for the development of potential inhibitors of SARS-CoV-2 Mpro to address COVID-19.

Small Molecule "Silmitasertib" Repurposed as Inhibitor of Transforming Growth Factor Beta 1 for the Development of Therapeutics for Oral Submucous Fibrosis

Oral submucous fibrosis (OSMF) is considered a premalignant condition characterized by aggressive fibrosis of the submucosal tissues of the oral cavity reflecting its malignant transformation potential. Activation of transforming growth factor beta (TGF-β) signaling has been reported to lead increased collagen production and fibrosis. Recently, significant upregulation of TGF-β1 has been reported in OSMF as compared to normal tissues. Therefore, inhibition of the TGF-β1 may pave for the development of therapeutics of OSMF. Based on the structure-assisted drug designing, we found "silmitasertib" as potent inhibitor of TGF-β1. We suggest that this molecule can be validated and implemented for the treatment of OSMF.

Antifungal activity of amino-alcohols based cationic surfactants and in silico, homology modeling, docking and molecular dynamics studies against lanosterol 14-α-demethylase enzyme

Fungi are being responsible for causing serious infections in humans and animals. The opportunistic microorganisms provoke environmental contaminations in health and storage facilities to represent a serious concern to health security. The present work investigates the antifungal activity of two amino-alcohols based cationic surfactants such as C_nEtOH, C_nPrOH (with n = 14 and 16 are the carbon numbers of alkyl chain and EtOH = Ethanol and PrOH = Propanol) against a collection of different Candida species (Candida tropicalis, Candida albicans, Candida auris, Cyberlindnera jadinii, Candida parapsilosis, Candida glabrata and Candida rugosa) respectively. The amino-alcohols based cationic surfactants exhibited good antifungal activity against all Candida strains tested with minimum inhibitory concentrations (MIC) ranging from 0.002 to 0.30 mM. The MIC evaluation shows an increase as a function of the hydrophobicity of all inhibitors against the majority of the Candida strains tested. The different location of the alcoholic OH function in the polar head shows the influence on the availability of N⁺ responsible for electrostatic interactions with the candidate's cell walls, which remains a very important step in the mode of action of quaternary ammonium cationic surfactants. Hence, a 3D structure of lanosterol 14-α-demethylase enzyme from C. auris was constructed by homology modeling using an online SWISS-MODEL server. The predicted model was analyzed by serval servers. Furthermore, a molecular docking study was carried out to better understand the binding mechanism of lanosterol homologous protein with surfactant ligands. Then, the docked complexes lanosterol-surfactants were refined by the molecular dynamic simulation to analyze their interaction behavior during the simulation.Communicated by Ramaswamy H. Sarma.

Comprehensive in silico analysis and molecular dynamics of the superoxide dismutase 1 (SOD1) variants related to amyotrophic lateral sclerosis

Amyotrophic Lateral Sclerosis (ALS) is the most frequent motor neuron disorder, with a significant social and economic burden. ALS remains incurable, and the only drugs approved for its treatments confers a survival benefit of a few months for the patients. Missense mutations in superoxide dismutase 1 (SOD1), a major cytoplasmic antioxidant enzyme, has been associated with ALS development, accounting for 23% of its familial cases and 7% of all sporadic cases. This work aims to characterize in silico the structural and functional effects of SOD1 protein variants. Missense mutations in SOD1 were compiled from the literature and databases. Twelve algorithms were used to predict the functional and stability effects of these mutations. ConSurf was used to estimate the evolutionary conservation of SOD1 amino-acids. GROMACS was used to perform molecular dynamics (MD) simulations of SOD1 wild-type and variants A4V, D90A, H46R, and I113T, which account for approximately half of all ALS-SOD1 cases in the United States, Europe, Japan, and United Kingdom, respectively. 233 missense mutations in SOD1 protein were compiled from the databases and literature consulted. The predictive analyses pointed to an elevated rate of deleterious and destabilizing predictions for the analyzed variants, indicating their harmful effects. The ConSurf analysis suggested that mutations in SOD1 mainly affect conserved and possibly functionally essential amino acids. The MD analyses pointed to flexibility and essential dynamics alterations at the electrostatic and metal-binding loops of variants A4V, D90A, H46R, and I113T that could lead to aberrant interactions triggering toxic protein aggregation. These alterations may have harmful implications for SOD1 and explain their association with ALS. Understanding the effects of SOD1 mutations on protein structure and function facilitates the design of further experiments and provides relevant information on the molecular mechanism of pathology, which may contribute to improvements in existing treatments for ALS.

Towards a novel peptide vaccine for Middle East respiratory syndrome coronavirus and its possible use against pandemic COVID-19

Middle East respiratory syndrome coronavirus (MERS-CoV) is an emerging health concern due to its high mortality rate of 35%. At present, no vaccine is available to protect against MERS-CoV infections. Therefore, an in silico search for potential antigenic epitopes in the non-redundant proteome of MERS-CoV was performed herein. First, a subtractive proteome-based approach was employed to look for the surface exposed and host non-homologous proteins. Following, immunoinformatics analysis was performed to predict antigenic B and T cell epitopes that were used in the design of a multi-epitopes peptide. Molecular docking study was carried out to predict vaccine construct affinity of binding to Toll-like receptor 3 (TLR3) and understand its binding conformation to extract ideas about its processing by the host immune system. We identified membrane protein, envelope small membrane protein, non-structural protein ORF3, non-structural protein ORF5, and spike glycoprotein as potential candidates for subunit vaccine designing. The designed multi-epitope peptide then linked to β-defensin adjuvant is showing high antigenicity. Further, the sequence of the designed vaccine construct is optimized for maximum expression in the Escherichia coli expression system. A rich pattern of hydrogen and hydrophobic interactions of the construct was observed with the TLR3 allowing stable binding of the construct at the docked site as predicted by the molecular dynamics simulation and MM-PBSA binding energies. We expect that the panel of subunit vaccine candidates and the designed vaccine construct could be highly effective in immunizing populations from infections caused by MERS-CoV and could possible applied on the current pandemic COVID-19.

Molecular level insights into the inhibition of gene expression by thiamine pyrophosphate (TPP) analogs for TPP riboswitch: A well-tempered metadynamics simulations study

Riboswitches are metabolite sensing aptamer domains present in non-coding regions in RNA and act as gene-regulating elements. Thiamine pyrophosphate (TPP) riboswitch is evolved as a new target for developing antibiotics against many pathogenic bacteria. The earlier reports suggest that the modification of the pyrophosphate group in the ligand molecule can enhance gene expression. In this work, we have examined the binding affinity and efficacy of TPP and two recently reported ligands, CH₂-TPP, and CF₂-TPP, using Well-tempered metadynamics (WT-MtD) simulations. The experimental in vitro assays show that both TPP and CH₂-TPP repress the gene expression to the same extent. The calculated binding energies correlate well with the experimental study and show the same trend of binding affinity of ligands for the TPP riboswitch. The root mean square fluctuation profiles suggest that the CH₂-TPP and TPP trigger higher fluctuations in P1 and L3 region, and such fluctuations in the P1 region is involved in the gene regulation process. The metal ion mediated contact of TPP ligand with pyrophosphate binding helix is found to be critical in the gene regulation process. The simulation results corroborate the experimental observations that the role of conformational changes occurring in different riboswitch regions upon ligand binding is essential to repress the gene expression process. This work sheds light on the subtle change in the ligand structure that can induce a more considerable impact on binding affinity and efficacy of ligands with riboswitch.

Biocomputational Prediction Approach Targeting FimH by Natural SGLT2 Inhibitors: A Possible Way to Overcome the Uropathogenic Effect of SGLT2 Inhibitor Drugs

The Food and Drug Administration (FDA) approved a new class of anti-diabetic medication (a sodium-glucose co-transporter 2 (SGLT2) inhibitor) in 2013. However, SGLT2 inhibitor drugs are under evaluation due to their associative side effects, such as urinary tract and genital infection, urinary discomfort, diabetic ketosis, and kidney problems. Even clinicians have difficulty in recommending it to diabetic patients due to the increased probability of urinary tract infection. In our study, we selected natural SGLT2 inhibitors, namely acerogenin B, formononetin, (-)-kurarinone, (+)-pteryxin, and quinidine, to explore their potential against an emerging uropathogenic bacterial therapeutic target, i.e., FimH. FimH plays a critical role in the colonization of uropathogenic bacteria on the urinary tract surface. Thus, FimH antagonists show promising effects against uropathogenic bacterial strains via their targeting of FimH's adherence mechanism with less chance of resistance. The molecular docking results showed that, among natural SGLT2 inhibitors, formononetin, (+)-pteryxin, and quinidine have a strong interaction with FimH proteins, with binding energy (∆G) and inhibition constant (ki) values of -5.65 kcal/mol and 71.95 µM, -5.50 kcal/mol and 92.97 µM, and -5.70 kcal/mol and 66.40 µM, respectively. These interactions were better than those of the positive control heptyl α-d-mannopyranoside and far better than those of the SGLT2 inhibitor drug canagliflozin. Furthermore, a 50 ns molecular dynamics simulation was conducted to optimize the interaction, and the resulting complexes were found to be stable. Physicochemical property assessments predicted little toxicity and good drug-likeness properties for these three compounds. Therefore, formononetin, (+)-pteryxin, and quinidine can be proposed as promising SGLT2 inhibitors drugs, with add-on FimH inhibition potential that might reduce the probability of uropathogenic side effects.