Molecular recognition of receptor sites using a genetic algorithm with a description of desolvation

In silico evaluation of lapachol derivatives binding to the Nsp9 of SARS-CoV-2

SARS-CoV-2 is the etiological agent of COVID-19, which represents a global health emergency that was rapidly declared a pandemic by the World Health Organization. Currently, there is a dearth of effective targeted therapies against viruses. Natural products isolated from traditional herbal plants have had a huge impact on drug development aimed at various diseases. Lapachol is a 1,4- naphthoquinone compound that has been demonstrated to have therapeutic effects against several diseases. SARS-CoV-2 non-structural proteins (nsps) play an important role in the viral replication cycle. Nsp9 seems to play a key role in transcription of the RNA genome of SARS-CoV-2. Virtual screening by docking and molecular dynamics suggests that lapachol derivatives can interact with Nsp9 from SARS-CoV-2. Complexes of lapachol derivatives V, VI, VIII, IX, and XI with the Nsp9 RNA binding site were subjected to molecular dynamics assays, to assess the stability of the complexes via RMSD. All complexes were stable over the course of 100 ns dynamics assays. Analyses of the hydrogen bonds in the complexes showed that lapachol derivatives VI and IX demonstrated strongest binding, with a stable or increasing number of hydrogen bonds over time. Our results demonstrate that Nsp9 from SARS-CoV-2 could be an important target in prospecting for ligands with antiviral potential. In addition, we showed that lapachol derivatives are potential ligands for SARS-CoV-2 Nsp9.

Communicated by Ramaswamy H. Sarma

Tubulin inhibitory activity of a novel colchicine-binding compounds based on a dinaphthospiropyranran scaffold

Spiropyrans have been investigated for their thermo- and photochromic characteristics, but their biotherapeutic properties have not been addressed. We report anti-proliferative properties of a novel dinaphthospiropyran analogue (1). The compound 1 was synthesized by a simple and expedient method using a one-pot acid-catalyzed aldol condensation of 2-hydroxy-1-naphthaldehyde with 4-piperidone followed by an acetalization reaction. Compound 1 was submitted to anticancer drug screen in the National Cancer Institute's panel of 60 human tumor cell lines. The average concentration of 1 to inhibit 50% cell growth was 5.4 ± 0.23 µM. All cell lines responded at almost the same concentration, suggesting that the action of 1 is not selective for cancer of origin. COMPARE analysis of dose-response data revealed interaction with tubulin as the possible mechanism of action of 1. At molecular level, 1 induced tubulin reorganization in colon cancer HCT-116 cells. Under cell-free conditions, the efficacy of 1 to inhibit tubulin polymerization was comparable to that of paclitaxel and vinblastine. Molecular docking showed that compound 1 binds to the colchicine-binding site of tubulin. We conclude that dinaphthospiropyrans present a novel scaffold for the development of tubulin inhibitors.

Integrated experimental/computational approaches to characterize the systems formed by vanadium with proteins and enzymes

An integrated instrumental/computational approach to characterize metallodrug–protein adducts at the molecular level is reviewed. A series of applications are described, focusing on potential vanadium drugs with a generalization to other metals.

First Synthesis of Racemic Trans Propargylamino-Donepezil, a Pleiotrope Agent Able to Both Inhibit AChE and MAO-B, with Potential Interest against Alzheimer's Disease

Alzheimer’s disease (AD) is a multifactorial neurodegenerative disease towards which pleiotropic approach using Multi-Target Directed Ligands is nowadays recognized as probably convenient. Among the numerous targets which are today validated against AD, acetylcholinesterase (ACh) and Monoamine Oxidase-B (MAO-B) appear as particularly convincing, especially if displayed by a sole agent such as ladostigil, currently in clinical trial in AD. Considering these results, we wanted to take benefit of the structural analogy lying in donepezil (DPZ) and rasagiline, two indane derivatives marketed as AChE and MAO-B inhibitors, respectively, and to propose the synthesis and the preliminary in vitro biological characterization of a structural compromise between these two compounds, we called propargylaminodonepezil (PADPZ). The synthesis of racemic trans PADPZ was achieved and its biological evaluation established its inhibitory activities towards both (h)AChE (IC₅₀ = 0.4 µM) and (h)MAO-B (IC₅₀ = 6.4 µM).

BioMetAll: Identifying Metal-Binding Sites in Proteins from Backbone Preorganization

With a large amount of research dedicated to decoding how metallic species bind to proteins, in silico methods are interesting allies for experimental procedures. To date, computational predictors mostly work by identifying the best possible sequence or structural match of the target protein with metal-binding templates. These approaches are fundamentally focused on the first coordination sphere of the metal. Here, we present the BioMetAll predictor that is based on a different postulate: the formation of a potential metal-binding site is related to the geometric organization of the protein backbone. We first report the set of convenient geometric descriptors of the backbone needed for the algorithm and their parameterization from a statistical analysis. Then, the successful benchmark of BioMetAll on a set of more than 90 metal-binding X-ray structures is presented. Because BioMetAll allows structural predictions regardless of the exact geometry of the side chains, it appears extremely valuable for systems whose structures (either experimental or theoretical) are not optimal for metal-binding sites. We report here its application on three different challenging cases: (i) the modulation of metal-binding sites during conformational transition in human serum albumin, (ii) the identification of possible routes of metal migration in hemocyanins, and (iii) the prediction of mutations to generate convenient metal-binding sites for de novo biocatalysts. This study shows that BioMetAll offers a versatile platform for numerous fields of research at the interface between inorganic chemistry and biology and allows to highlight the role of the preorganization of the protein backbone as a marker for metal binding. BioMetAll is an open-source application available at https://github.com/insilichem/biometall.

Development of 2,5-dihydro-4H-pyrazolo[3,4-d]pyrimidin-4-one inhibitors of aldehyde dehydrogenase 1A (ALDH1A) as potential adjuncts to ovarian cancer chemotherapy

There is strong evidence that inhibition of one or more Aldehyde Dehydrogenase 1A (ALDH1A) isoforms may be beneficial in chemotherapy-resistant ovarian cancer and other tumor types. While many previous efforts have focused on development of ALDH1A1 selective inhibitors, the most deadly ovarian cancer subtype, high-grade serous (HGSOC), exhibits elevated expression of ALDH1A3. Herein, we report continued development of pan-ALDH1A inhibitors to assess whether broad spectrum ALDH1A inhibition is an effective adjunct to chemotherapy in this critical tumor subtype. Optimization of the CM39 scaffold, aided by metabolite ID and several new ALDH1A1 crystal structures, led to improved biochemical potencies, improved cellular ALDH inhibition in HGSOC cell lines, and substantial improvements in microsomal stability culminating in orally bioavailable compounds. We demonstrate that two compounds 68 and 69 are able to synergize with chemotherapy in a resistant cell line and patient-derived HGSOC tumor spheroids, indicating their suitability for future in vivo proof of concept experiments.

Pleiotropic prodrugs: Design of a dual butyrylcholinesterase inhibitor and 5-HT6 receptor antagonist with therapeutic interest in Alzheimer's disease

Beside acetylcholinesterase, butyrylcholinesterase could be considered as a putative target of interest for the symptomatic treatment of Alzheimer's disease (AD). As a result of complexity of AD, no molecule has been approved since 2002. Idalopirdine, a 5-HT₆ receptors antagonist, did not show its effectiveness in clinical trial despite its evaluation as adjunct to cholinesterase inhibitors. Pleiotropic molecules, known as multitarget directed ligands (MTDLs) are currently developed to tackle the multifactorial origin of AD. In this context, we have developed a pleiotropic carbamate 7, that behaves as a covalent inhibitor of BuChE (IC₅₀ = 0.97 μM). The latter will deliver after hydrolysis, compound 6, a potent 5-HT₆ receptors antagonist (K_i = 11.4 nM) related to idalopirdine. In silico and in vitro evaluation proving our concept were performed completed with first in vivo results that demonstrate great promise in restoring working memory.

The determination of Plutella xylostella (L.) GSTs (PxGSTs) involved in the detoxification metabolism of Tolfenpyrad

BACKGROUND

Insect glutathione S-transferases (GSTs) play a crucial role in insecticide detoxification. However, there remains a distinct lack of information regarding the role of GSTs in the detoxification of Tolfenpyrad (TFP) in insects.

RESULTS

Real-time quantitative PCR showed significant upregulation of PxGSTs after exposure to TFP for 6 h. An in vitro inhibition assay showed that TFP could inhibit PxGSTδ, PxGSTε and PxGSTσ, and the most pronounced inhibitory effect was on PxGSTσ. Metabolism assays displayed that PxGSTσ was superior to other test PxGSTs in metabolizing TFP. The molecular docking of TFP and PxGSTσ revealed that the H-bond provided by the sidechains of Tyr107 and Tyr162 were key to the detoxification of TFP by PxGSTσ. Further tests using mutant PxGSTσ proteins at the sites of Tyr107 (PxGSTσY107A) and Tyr162 (PxGSTσY162A) corroborated that the individual replacement of Tyr107 and Tyr162 could greatly weaken the binding and metabolic abilities to TFP.

CONCLUSION

Metabolic interactions between the Plutella xylostella (L.) GSTs (PxGSTs) and TFP were deciphered. This study illustrates the molecular metabolism mechanism of PxGSTσ towards TFP and provides theoretical underpinnings for the design and optimization of novel TFP-like insecticides. © 2020 Society of Chemical Industry.

Integron gene cassettes harboring novel variants of D-alanine-D-alanine ligase confer high-level resistance to D-cycloserine

Antibiotic resistance poses an increasing threat to global health. To tackle this problem, the identification of principal reservoirs of antibiotic resistance genes (ARGs) plus an understanding of drivers for their evolutionary selection are important. During a PCR-based screen of ARGs associated with integrons in saliva-derived metagenomic DNA of healthy human volunteers, two novel variants of genes encoding a d-alanine-d-alanine ligase (ddl6 and ddl7) located within gene cassettes in the first position of a reverse integron were identified. Treponema denticola was identified as the likely host of the ddl cassettes. Both ddl6 and ddl7 conferred high level resistance to d-cycloserine when expressed in Escherichia coli with ddl7 conferring four-fold higher resistance to D-cycloserine compared to ddl6. A SNP was found to be responsible for this difference in resistance phenotype between the two ddl variants. Molecular dynamics simulations were used to explain the mechanism of this phenotypic change at the atomic scale. A hypothesis for the evolutionary selection of ddl containing integron gene cassettes is proposed, based on molecular docking of plant metabolites within the ATP and d-cycloserine binding pockets of Ddl.

Predicting the Substrate Scope of the Flavin-Dependent Halogenase BrvH

The recently described flavin-dependent halogenase BrvH is able to catalyse both the bromination and chlorination of indole, but shows significantly higher bromination activity. BrvH was annotated as a tryptophan halogenase, but does not accept tryptophan as a substrate. Its native substrate remains unknown. A predictive model with the data available for BrvH was analysed. A training set of compounds tested in vitro was docked into the active site of a complete protein model based on the X-ray structure of BrvH. The atoms not resolved experimentally were modelled by using molecular mechanics force fields to obtain this protein model. Furthermore, docking poses for the substrates and known non-substrates have been calculated. Parameters like distance, partial charge and hybridization state were analysed to derive rules for predicting activity. With this model for activity of the BrvH, a virtual screening suggested several structures for potential substrates. Some of the compounds preselected in this way were tested in vitro, and several could be verified as convertible substrates. Based on information on halogenated natural products, a new dataset was created to specifically search for natural products as substrates/products, and virtual screening in this database yielded further hits.

A multi-stage virtual screening of FDA-approved drugs reveals potential inhibitors of SARS-CoV-2 main protease

Coronavirus disease 2019 (COVID-19) caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is an ongoing global health emergency. Repurposing of approved pharmaceutical drugs for COVID-19 treatment represents an attractive approach to quickly identify promising drug candidates. SARS-CoV-2 main protease (M^pro) is responsible for the maturation of viral functional proteins making it a key antiviral target. Based on the recently revealed crystal structures of SARS-CoV-2 M^pro, we herein describe a multi-stage virtual screening protocol including pharmacophore screening, molecular docking and protein-ligand interaction fingerprints (PLIF) post-docking filtration for efficient enrichment of potent SARS-CoV-2 M^pro inhibitors. Potential hits, along with a cocrystallized control were further studied via molecular dynamics. A 150-ns production trajectory was followed by RMSD, free energy calculation, and H-bond analysis for each compound. The applied virtual screening protocol led to identification of five FDA-approved drugs with promising binding modes to key subsites of the substrate-binding pocket of SARS-CoV-2 M^pro. The identified compounds belong to different pharmaceutical classes, including several protease inhibitors, antineoplastic agents and a natural flavonoid. The drug candidates discovered in this study present a potential extension of the recently reported SARS-CoV-2 M^pro inhibitors that have been identified using other virtual screening protocols and may be repurposed for COVID-19 treatment.

Design, synthesis and biological evaluation of novel naphthoquinone-4-aminobenzensulfonamide/carboxamide derivatives as proteasome inhibitors

A series of novel 4-aminobenzensulfonamide/carboxamide derivatives bearing naphthoquinone pharmacophore were designed, sythesized and evaluated for their proteasome inhibitory and antiproliferative activities against human breast cancer cell line (MCF-7). The structures of the synthesized compounds were confirmed by spectral and elemental analyses. The proteasome inhibitory activity studies were carried out using cell-based assay. The antiproteasomal activity results revealed that most of the compounds exhibited inhibitory activity with different percentages against the caspase-like (C-L, β1 subunit), trypsin-like (T-L, β2 subunit) and chymotrypsin-like (ChT-L, β5 subunit) activities of proteasome. Among the tested compounds, compound 14 bearing 5-chloro-2-pyridyl ring on the nitrogen atom of sulfonamide group is the most active compound in the series and displayed higher inhibition with IC₅₀ values of 9.90 ± 0.61, 44.83 ± 4.23 and 22.27 ± 0.15 μM against ChT-L, C-L and T-L activities of proteasome compared to the lead compound PI-083 (IC₅₀ = 12.47 ± 0.21, 53.12 ± 2.56 and 26.37 ± 0.5 μM), respectively. The antiproliferative activity was also determined by MTT (3-(4,5-Dimethyl-2-thiazolyl)-2,5-diphenyl-2H-tetrazolium bromide) assay in vitro. According to the antiproliferative activity results, all of the compounds exhibited cell growth inhibitory activity in a range of IC₅₀ = 1.72 ± 0.14-20.8 ± 0.5 μM and compounds 13 and 28 were found to be the most active compounds with IC₅₀ values of 1.79 ± 0.21 and 1.72 ± 0.14 μM, respectively. Furthermore, molecular modeling studies were carried out for the compounds 13, 14 and 28 to investigate the ligand-enzyme binding interactions.

Structural and Functional Insights into an Archaeal Lipid Synthase

The UbiA superfamily of intramembrane prenyltransferases catalyzes an isoprenyl transfer reaction in the biosynthesis of lipophilic compounds involved in cellular physiological processes. Digeranylgeranylglyceryl phosphate (DGGGP) synthase (DGGGPase) generates unique membrane core lipids for the formation of the ether bond between the glycerol moiety and the alkyl chains in archaea and has been confirmed to be a member of the UbiA superfamily. Here, the crystal structure is reported to exhibit nine transmembrane helices along with a large lateral opening covered by a cytosolic cap domain and a unique substrate-binding central cavity. Notably, the lipid-bound states of this enzyme demonstrate that the putative substrate-binding pocket is occupied by the lipidic molecules used for crystallization, indicating the binding mode of hydrophobic substrates. Collectively, these structural and functional studies provide not only an understanding of lipid biosynthesis by substrate-specific lipid-modifying enzymes but also insights into the mechanisms of lipid membrane remodeling and adaptation.

On single-electron magnesium bonding formation and the effect of methyl substitution

The complexes formed between MgX₂ (X = F, H) molecules and alkyl radicals Y [Y = CH₃, CH₂CH₃, CH(CH₃)₂, and C(CH₃)₃] have been characterized by using quantum chemical methods. The binding distance in all cases is less than the sum of vdW radii of Mg and C, indicating the formation of a non-covalent interaction, namely single-electron magnesium bond. Energy decomposition analysis reveals that electrostatic and polarization contributions are the major components responsible for the stability of the studied complexes. According to interaction energy, atoms in molecules, and independent gradient model analyses, methyl substitution on electron donor Y imposes a positive effect on its complexation with MgX₂. When compared with other nonbonded interactions, the single-electron magnesium bond is found to have strength comparable to those of the single-electron beryllium bond and π-magnesium bond.

The complexes formed between MgX₂ (X = F, H) molecules and alkyl radicals Y [Y = CH₃, CH₂CH₃, CH(CH₃)₂, and C(CH₃)₃] have been characterized by using quantum chemical methods.

Design and synthesis of Coenzyme A analogues as Aurora kinase A inhibitors: An exploration of the roles of the pyrophosphate and pantetheine moieties

Coenzyme A (CoA) is a highly selective inhibitor of the mitotic regulatory enzyme Aurora A kinase, with a novel mode of action. Herein we report the design and synthesis of analogues of CoA as inhibitors of Aurora A kinase. We have designed and synthesised modified CoA structures as potential inhibitors, combining dicarbonyl mimics of the pyrophosphate group with a conserved adenosine headgroup and different length pantetheine-based tail groups. An analogue with a -SH group at the end of the pantotheinate tail showed the best IC50, probably due to the formation of a covalent bond with Aurora A kinase Cys290.

Deciphering the Inhibition Mechanism of under Trial Hsp90 Inhibitors and Their Analogues: A Comparative Molecular Dynamics Simulation

Heat shock protein 90 (Hsp90) performs functions in cellular activities together with other signaling pathways. Hsp90 is evolutionarily conserved and universally articulated as a human cancer-causing agent involved in lung cancer and breast cancer followed by colon and rectum cancers. It has emerged as an effective drug candidate, and inhibition may affect several signaling pathways associated with cancer spread. Therefore, in-silico approaches, molecular docking, molecular dynamics simulation, and binding free energy calculations were applied to create insights into the inhibition mechanism against Hsp90 to identify new cancer therapeutic drugs. Top-docked Hsp90-inhibitor complexes with their analogues were selected as the best complexes based on the GOLD fitness score and orientation. The significant interaction of Hsp90 inhibitors and their analogues were observed to be bound with active site residues as well as residing within the same cavity region. System stability factors RMSD, RMSF, beta-factor, and radius of gyration were analyzed for top-docked complexes and ensure strong binding interaction between inhibitors and the Hsp90 cavity. Cavity bound inhibitors were found to retain consistent hydrogen bonding during the simulation. The radial distribution function (RDF) illustrated that interacting active site residues drive the binding and stability of the inhibitors. Similarly, the axial frequency distribution, which is an indigenously developed analytical tool, produced noteworthy knowledge of the hydrogen-bonding pattern. Results yielded new insights into the design of cancer therapeutic drugs against Hsp90. This finding suggests that under trial Hsp90 inhibitors MPC-3100 could be a potential starting point into the development of potential anticancer agents with the possibility of future directions for the improvement of early existing Hsp90 inhibitors CNF-2024 and SNX-5422 as an anticancer agent.

Unique Domain for a Unique Target: Selective Inhibitors of Host Cell DDX3X to Fight Emerging Viruses

Emerging viruses like dengue, West Nile, chikungunya, and Zika can cause widespread viral epidemics. Developing novel drugs or vaccines against specific targets for each virus is a difficult task. As obligate parasites, all viruses exploit common cellular pathways, providing the possibility to develop broad-spectrum antiviral agents targeting host factors. The human DEAD-box RNA helicase DDX3X is an essential cofactor for viral replication but dispensable for cell viability. Herein, we exploited the presence of a unique structural motif of DDX3X not shared by other cellular enzymes to develop a theoretical model to aid in the design of a novel class of highly selective inhibitors acting against such specific targets, thus limiting off-targeting effects. High-throughput virtual screening led us to identify hit compound 5, endowed with promising antienzymatic activity. To improve its aqueous solubility, 5 and its two enantiomers were synthesized and converted into their corresponding acetate salts (compounds 11, 12, and 13). In vitro mutagenesis and biochemical and cellular assays further confirmed that the developed molecules were selective for DDX3X and were able to suppress replication of West Nile and dengue viruses in infected cells in the micromolar range while showing no toxicity for uninfected cells. These results provide proof of principle for a novel strategy in developing highly selective and broad-spectrum antiviral molecules active against emerging and dangerous viral pathogens. This study paves the way for the development of larger focused libraries targeting such domain to expand SAR studies and fully characterize their mode of interaction.

In silico molecular modelling, structural dynamics simulation and characterization of antifungal nature of β-glucosidase enzyme from Sechium edule

β-glucosidase is an enzyme that has ability to cleave β-glycosidic bonds present in oligosaccharides and glycoconjugates. They are known to be present across all domains of living organism and have important roles in many biological processes including plant defense mechanism. In the present study, a β-glucosidase enzyme identified from seeds of Sechium edule was characterized using various bioinformatics tools. A homology model (SeBG) was generated using a β-glucosidase crystal structure from Oryza sativa (PDB ID: 3PTK) as template. In silico structural binding studies on putative β-glucosidase protein revealed a stable and strong interaction indicative of higher GOLD fitness score with the substrates: p-nitrophenyl-β-d-glucopyranoside (pNPG), laminarin, chitotriose, N-acetylglucosamine and N-acetylmuramic acid suggesting its possible role in broad spectrum antifungal and antimicrobial activity. Assessment of the in vitro enzyme activity with pNPG showed a K_m and V_max values of 2.7 mM and 22 µMmin^-1mL^-1mg^-1, respectively. While, the in vitro enzyme activity with laminarin showed a K_m and V_max values of 0.31 mM and 0.043 µMmin^-1mL^-1mg^-1. The broad spectrum activity of the protein shown in our result indicates SeBG as a promising biocontrol agent against phytopathogens.Communicated by Ramaswamy H. Sarma.

How 'Protein-Docking' Translates into the New Emerging Field of Docking Small Molecules to Nucleic Acids?

In this review, we retraced the '40-year evolution' of molecular docking algorithms. Over the course of the years, their development allowed to progress from the so-called 'rigid-docking' searching methods to the more sophisticated 'semi-flexible' and 'flexible docking' algorithms. Together with the advancement of computing architecture and power, molecular docking's applications also exponentially increased, from a single-ligand binding calculation to large screening and polypharmacology profiles. Recently targeting nucleic acids with small molecules has emerged as a valuable therapeutic strategy especially for cancer treatment, along with bacterial and viral infections. For example, therapeutic intervention at the mRNA level allows to overcome the problematic of undruggable proteins without modifying the genome. Despite the promising therapeutic potential of nucleic acids, molecular docking programs have been optimized mostly for proteins. Here, we have analyzed literature data on nucleic acid to benchmark some of the widely used docking programs. Finally, the comparison between proteins and nucleic acid targets docking highlighted similarity and differences, which are intrinsically related to their chemical and structural nature.

Homology Modeling of the Human P-glycoprotein (ABCB1) and Insights into Ligand Binding through Molecular Docking Studies

The ABCB1 transporter also known as P-glycoprotein (P-gp) is a transmembrane protein belonging to the ATP binding cassette super-family of transporters; it is a xenobiotic efflux pump that limits intracellular drug accumulation by pumping the compounds out of cells. P-gp contributes to a decrease of toxicity and possesses broad substrate specificity. It is involved in the failure of numerous anticancer and antiviral chemotherapies due to the multidrug resistance (MDR) phenomenon, where it removes the chemotherapeutics out of the targeted cells. Understanding the details of the ligand–P-gp interaction is therefore crucial for the development of drugs that might overcome the MRD phenomenon and for obtaining a more effective prediction of the toxicity of certain compounds. In this work, an in silico modeling was performed using homology modeling and molecular docking methods with the aim of better understanding the ligand–P-gp interactions. Based on different mouse P-gp structural templates from the PDB repository, a 3D model of the human P-gp (hP-gp) was constructed by means of protein homology modeling. The homology model was then used to perform molecular docking calculations on a set of thirteen compounds, including some well-known compounds that interact with P-gp as substrates, inhibitors, or both. The sum of ranking differences (SRD) was employed for the comparison of the different scoring functions used in the docking calculations. A consensus-ranking scheme was employed for the selection of the top-ranked pose for each docked ligand. The docking results showed that a high number of π interactions, mainly π–sigma, π–alkyl, and π–π type of interactions, together with the simultaneous presence of hydrogen bond interactions contribute to the stability of the ligand–protein complex in the binding site. It was also observed that some interacting residues in hP-gp are the same when compared to those observed in a co-crystallized ligand (PBDE-100) with mouse P-gp (PDB ID: 4XWK). Our in silico approach is consistent with available experimental results regarding P-gp efflux transport assay; therefore it could be useful in the prediction of the role of new compounds in systemic toxicity.

Extracellular alpha/beta-hydrolase from Paenibacillus species shares structural and functional homology to tobacco salicylic acid binding protein 2

An alpha/ beta hydrolase annotated as a putative salicylate esterase within the genome of a species of Paenibacillus previously identified from differential and selective growth on Kraft lignin was structurally and functionally characterised. Feruloyl esterases are key to the degradation of lignin in several bacterial species and although this activity was investigated, no such activity was observed. The crystal structure of the Paenibacillus esterase, here denoted as PnbE, was determined at 1.32 Å resolution, showing high similarity to Nicotiana tabacum salicylic acid binding protein 2 from the protein database. Structural similarities between these two structures across the core domains and key catalytic residues were observed, with superposition of catalytic residues giving an RMSD of 0.5 Å across equivalent Cα atoms. Conversely, the cap domains of PnbE and Nicotiana tabacum SABP2 showed greater divergence with decreased flexibility in the PnbE cap structure. Activity of PnbE as a putative methyl salicylate esterase was supported with binding studies showing affinity for salicylic acid and functional studies showing methyl salicylate esterase activity. We hypothesise that this activity could enrich Paenibacillus sp. within the rhizosphere by increasing salicylic acid concentrations within the soil.

GM-DockZn: a geometry matching-based docking algorithm for zinc proteins

Motivation

Molecular docking is a widely used technique for large-scale virtual screening of the interactions between small-molecule ligands and their target proteins. However, docking methods often perform poorly for metalloproteins due to additional complexity from the three-way interactions among amino-acid residues, metal ions and ligands. This is a significant problem because zinc proteins alone comprise about 10% of all available protein structures in the protein databank. Here, we developed GM-DockZn that is dedicated for ligand docking to zinc proteins. Unlike the existing docking methods developed specifically for zinc proteins, GM-DockZn samples ligand conformations directly using a geometric grid around the ideal zinc-coordination positions of seven discovered coordination motifs, which were found from the survey of known zinc proteins complexed with a single ligand.

Results

GM-DockZn has the best performance in sampling near-native poses with correct coordination atoms and numbers within the top 50 and top 10 predictions when compared to several state-of-the-art techniques. This is true not only for a non-redundant dataset of zinc proteins but also for a homolog set of different ligand and zinc-coordination systems for the same zinc proteins. Similar superior performance of GM-DockZn for near-native-pose sampling was also observed for docking to apo-structures and cross-docking between different ligand complex structures of the same protein. The highest success rate for sampling nearest near-native poses within top 5 and top 1 was achieved by combining GM-DockZn for conformational sampling with GOLD for ranking. The proposed geometry-based sampling technique will be useful for ligand docking to other metalloproteins.

Availability and implementation

GM-DockZn is freely available at www.qmclab.com/ for academic users.

Supplementary information

Supplementary data are available at Bioinformatics online.

Structure-Based Screening of Non-β-Lactam Inhibitors against Class D β-Lactamases: An Approach of Docking and Molecular Dynamics

The manifestation of class D β-lactamases in the community raises significant concern as they can hydrolyze carbapenem antibiotics. Hence, it is exceptionally alluring to design novel inhibitors. Structure-based virtual screening using docking programs and molecular dynamics simulations was employed to identify two novel non-β-lactam compounds that possess the ability to block different OXA variants. Furthermore, the presence of a nonpolar aliphatic amino acid, valine, near the active site serine, was identified in all OXA variants that can be accounted to block the catalytic activity of OXA enzymes.

Computational insight into the interaction of oxaliplatin with insulin

In an organism, cisplatin and its derivatives are known to interact with proteins besides their principal DNA target. These off-target interactions have major therapeutic consequences including undesired side effects, loss of bioavailability and emergence of resistance. Insulin is one of the prototypical protein targets of platinum drugs as it has been seen to be involved in bioavailability reduction and might also determine resistance in certain cancer lines. However, despite the interest in understanding the nature of the oxaliplatin-insulin adducts, no 3D models have been achieved so far. In this study, we apply our recent computational multiscale protocol optimized for bioinorganic interactions to provide structural insights into these systems. To do so, the initial structures are predicted by blind protein-metalloligand docking calculations optimized to account for a metal-containing species, and then refined using a Molecular Dynamics (MD) and Quantum Mechanics/Molecular Mechanics (QM/MM) integrated protocol. The results are consistent with experimental information obtained from fragment analysis, and also provide novel structural information like conformational changes occurring upon binding and potential effects on the biological functions of the protein. This study opens an avenue towards applying similar strategies to a wide ensemble of metallodrug-protein/peptide systems for which no structural data are available.

Developing an effective polarizable bond method for small molecules with application to optimized molecular docking

Electrostatic interaction plays an essential role in protein-ligand binding. Due to the polarization effect, electrostatic interactions are largely impacted by their local environments. However, traditional force fields use fixed point charge-charge interactions to describe electrostatic interactions but is unable to include the polarization effect. The lack of the polarization effect in the force field representation can result in substantial error in biomolecular studies, such as molecular dynamics and molecular docking. Docking programs usually employ traditional force fields to estimate the binding energy between a ligand and a protein for pose selection or scoring. The intermolecular interaction energy mainly consists of van der Waals and electrostatic interaction in the force field representation. In the current study, we developed an Effective Polarizable Bond (EPB) method for small organic molecules and applied this EPB method to optimize protein-ligand docking in computational tests for a variety of protein-ligand systems. We tested the method on a set of 38 cocrystallized structures taken from the Protein Data Bank (PDB) and found that the maximum error was reduced from 7.98 Å to 2.03 Å when using EPB Dock, providing strong evidence that the use of EPB charges is important. We found that our optimized docking approach with EPB charges could improve the docking performance, sometimes dramatically, and the maximum error was reduced from 12.88 Å to 1.57 Å in Optimized Docking (in the case of 1fqx). The average RMSD decreased from 2.83 Å to 1.85 Å. Further investigations showed that the use of the EBP method could enhance intermolecular hydrogen bonding, which is a major contributing factor to improved docking performance. Developed tools for the calculation of the polarized ligand charge from a protein-ligand complex structure with the EPB method are freely available on GitHub (https://github.com/Xundrug/EPB).

Improving Protein-Peptide Docking Results via Pose-Clustering and Rescoring with a Combined Knowledge-Based and MM-GBSA Scoring Function

Protein-peptide docking, which predicts the complex structure between a protein and a peptide, is a valuable computational tool in peptide therapeutics development and the mechanistic investigation of peptides involved in cellular processes. Although current peptide docking approaches are often able to sample near-native peptide binding modes, correctly identifying those near-native modes from decoys is still challenging because of the extremely high complexity of the peptide binding energy landscape. In this study, we have developed an efficient postdocking rescoring protocol using a combined scoring function of knowledge-based ITScorePP potentials and physics-based MM-GBSA energies. Tested on five benchmark/docking test sets, our postdocking strategy showed an overall significantly better performance in binding mode prediction and score-rmsd correlation than original docking approaches. Specifically, our postdocking protocol outperformed original docking approaches with success rates of 15.8 versus 10.5% for pepATTRACT on the Global_57 benchmark, 5.3 versus 5.3% for CABS-dock on the Global_57 benchmark, 17.0 versus 11.3% for FlexPepDock on the LEADS-PEP data set, 40.3 versus 33.9% for HPEPDOCK on the Local_62 benchmark, and 64.2 versus 52.8% for HPEPDOCK on the LEADS-PEP data set when the top prediction was considered. These results demonstrated the efficacy and robustness of our postdocking protocol.

Combining Molecular Dynamics and Docking Simulations to Develop Targeted Protocols for Performing Optimized Virtual Screening Campaigns on The hTRPM8 Channel

BACKGROUND

There is an increasing interest in TRPM8 ligands of medicinal interest, the rational design of which can be nowadays supported by structure-based in silico studies based on the recently resolved TRPM8 structures. Methods: The study involves the generation of a reliable hTRPM8 homology model, the reliability of which was assessed by a 1.0 μs MD simulation which was also used to generate multiple receptor conformations for the following structure-based virtual screening (VS) campaigns; docking simulations utilized different programs and involved all monomers of the selected frames; the so computed docking scores were combined by consensus approaches based on the EFO algorithm. Results: The obtained models revealed very satisfactory performances; LiGen™ provided the best results among the tested docking programs; the combination of docking results from the four monomers elicited a markedly beneficial effect on the computed consensus models. Conclusions: The generated hTRPM8 model appears to be amenable for successful structure-based VS studies; cross-talk modulating effects between interacting monomers on the binding sites can be accounted for by combining docking simulations as performed on all the monomers; this strategy can have general applicability for docking simulations involving quaternary protein structures with multiple identical binding pockets.

Computational insights into the molecular mechanisms of differentiated allosteric modulation at the mu opioid receptor by structurally similar bitopic modulators

Targeting the mu opioid receptor (MOR) by applying orthosteric ligands is the most frequently employed method to treat opioid use disorder (OUD). Unfortunately, most of MOR orthosteric ligands produce severe side effects, mainly due to their low selectivity over other opioid receptors. In contrast, some G protein-coupled receptor allosteric modulators have been reported to exhibit high subtype selectivity and can effectively modulate the potency and/or efficacy of orthosteric ligands. Recently, NAQ and its analog NCQ were identified as novel MOR bitopic modulators. Interestingly, NAQ and NCQ were similar in structure but exhibited different efficacy profiles to the MOR. NAQ exhibited an antagonism activity to the MOR while NCQ showed a partial agonism activity to the MOR. In the present study, molecular modeling methods were applied to explore the putative molecular mechanisms of their different functional profiles to the MOR. When NAQ binding with the inactive MOR, the 'address' portion of NAQ interacted with the MOR allosteric site but showed no significant allosteric modulation of the efficacy of the 'message' portion of NAQ. However, when NCQ binding with the inactive and active MOR, the 'address' portion of NCQ seemed to be able to positively modulate the efficacy of the 'message' portion of NCQ at varying levels. Evidentially, the substituents at the 1'- and 4'-positions of the isoquinoline ring of NCQ seemed to play a critical role in the modulatory function of the 'address' portion of NCQ. These findings will be invaluable to develop our next generation of MOR bitopic modulators with high affinity and subtype selectivity to potentially treat OUD.

Role of ring positioning and preferential occupation of ligand obtained through molecular dynamics simulation of peptidoglycan associated lipoprotein (Pal)

Peptidoglycan associated lipoprotein (Pal) of Burkholderia cepacia is a highly attractive target as it is essential and conserved across strains of the pathogen and involved significantly in pathogenesis and survival. Herein, we used this target in virtual screening of small drug molecules library with an objective to screen the most appropriate antagonist for (Pal) to block its biological function. Results yielded a pyrimidine with imidazole group of compound "(R)-2-(1-(2, 6-dioxo-1, 2, 3, 6-tetrahydropyrimidine-4-carbonyl) piperidin-3-yl)-5-methlimidazol-3-ium-1-ide" showing the highest affinity for the active pocket of the protein and posed as such that the pyrimidine-4-carbonyl and piperidin-3-yl ring of the molecule is docked deep inside the pocket keeping the (R)-2-(piperidin-3-yl)imidazole-3-ium-1-ide posed outside this cavity. Molecular dynamics (MD) simulations illustrated stable binding of the compound at the docked site throughout the course of 600 ns (ns) with no conformational changes observed in the receptor macromolecule. Higher conformational stability was inferred for the system with mean RMSD of 1.1 Å supported with residual flexibility with mean value of 0.7 Å. Radial distribution function (RDF) showed the profound involvement of Arg56, Asp20 and Asp54 from protein active pocket in hydrogen bonding at the beginning and towards the end of simulation. The binding energy MMPB/GBSA and WaterSwap analyses showed an earnest contract on the affinity of inhibitor in response to hotspot residues of the active pocket. The role of stable positioning of rings inside the cavity is established through simulation, while a relevant comparative analysis indicates the compound as biologically active and potent candidate for therapeutic applications.

Improvement of the novel inhibitor for Mycobacterium enoyl-acyl carrier protein reductase (InhA): a structure-activity relationship study of KES4 assisted by in silico structure-based drug screening

InhA or enoyl-acyl carrier protein reductase of Mycobacterium tuberculosis (mtInhA), which controls mycobacterial cell wall construction, has been targeted in the development of antituberculosis drugs. Previously, our in silico structure-based drug screening study identified a novel class of compounds (designated KES4), which is capable of inhibiting the enzymatic activity of mtInhA, as well as mycobacterial growth. The compounds are composed of four ring structures (A-D), and the MD simulation predicted specific interactions with mtInhA of the D-ring and methylene group between the B-ring and C-ring; however, there is still room for improvement in the A-ring structure. In this study, a structure-activity relationship study of the A-ring was attempted with the assistance of in silico docking simulations. In brief, the virtual chemical library of A-ring-modified KES4 was constructed and subjected to in silico docking simulation against mtInhA using the GOLD program. Among the selected candidates, we achieved synthesis of seven compounds, and the bioactivities (effects on InhA activity and mycobacterial growth and cytotoxicity) of the synthesized molecules were evaluated. Among the compounds tested, two candidates (compounds 3d and 3f) exhibited superior properties as mtInhA-targeted anti-infectives for mycobacteria than the lead compound KES4.

A review of mathematical representations of biomolecular data

Recently, machine learning (ML) has established itself in various worldwide benchmarking competitions in computational biology, including Critical Assessment of Structure Prediction (CASP) and Drug Design Data Resource (D3R) Grand Challenges. However, the intricate structural complexity and high ML dimensionality of biomolecular datasets obstruct the efficient application of ML algorithms in the field. In addition to data and algorithm, an efficient ML machinery for biomolecular predictions must include structural representation as an indispensable component. Mathematical representations that simplify the biomolecular structural complexity and reduce ML dimensionality have emerged as a prime winner in D3R Grand Challenges. This review is devoted to the recent advances in developing low-dimensional and scalable mathematical representations of biomolecules in our laboratory. We discuss three classes of mathematical approaches, including algebraic topology, differential geometry, and graph theory. We elucidate how the physical and biological challenges have guided the evolution and development of these mathematical apparatuses for massive and diverse biomolecular data. We focus the performance analysis on protein-ligand binding predictions in this review although these methods have had tremendous success in many other applications, such as protein classification, virtual screening, and the predictions of solubility, solvation free energies, toxicity, partition coefficients, protein folding stability changes upon mutation, etc.

Rational design of novel benzisoxazole derivatives with acetylcholinesterase inhibitory and serotoninergic 5-HT4 receptors activities for the treatment of Alzheimer's disease

A rigidification strategy was applied to the preclinical candidate donecopride, an acetylcholinesterase inhibitor possessing 5-HT₄R agonist activity. Inspired by promising bioactive benzisoxazole compounds, we have conducted a pharmacomodulation study to generate a novel series of multitarget directed ligands. The chemical synthesis of the ligand was optimized and compounds were evaluated in vitro against each target and in cellulo. Structure-activity relationship was supported by docking analysis in human acetylcholinesterase binding site. Among the synthesized compounds, we have identified a novel hybrid 32a (3-[2-[1-(cyclohexylmethyl)-4-piperidyl]ethyl]-4-methoxy-1,2-benzoxazole) able to display nanomolar acetylcholinesterase inhibitory effects and nanomolar Ki for 5-HT₄R.

Synthesis, biological evaluations and molecular modelling studies of novel indolin-2-ones designing as FGFR inhibitors

A series of novel 3,5-disubstituted indolin-2-ones were designed and synthesized as selective FGFR inhibitors. In the design process of 3,5-disubstituted indolin-2-ones for FGFRs, molecular docking studies were performed to generate and optimize novel compounds which have FGFR inhibitory potency, theoretically. In vitro enzyme inhibitory and selectivity profiles of the synthesized compounds, and their cytotoxicity against NIH-3T3 cells were evaluated. According to enzyme inhibition assay, compound A1 (FGFR1-4; IC₅₀ = 19.82; 5.95; 1419; 37150 nM), compound A5 (FGFR1-4; IC₅₀ = 1890; Nd; 6.50; 18590 nM) and compound A13 (FGFR1-4; IC₅₀ = 6.99; 1022; 17090; 8993 nM) have displayed best inhibitory potency against FGFR2, FGFR3 and FGFR1, respectively. The studied compounds have displayed low affinity to FGFR4 in comparison with other isoforms. Molecular docking study data were used to determine the binding orientations of the synthesized compounds inside FGFRs in accordance with enzyme inhibition assay data. Molecular dynamics simulations and free energy calculations were performed to determine stability, binding modes and dynamics behaviors of compound A1, A5 and A13 inside FGFR-2, FGFR-3 and FGFR-1, respectively. The compounds bearing aromatic groups at the C5 position of indolin-2-one could be lead compounds for the development of more effective and selective FGFR1-3 inhibitors.

Toward a Structural View of hERG Activation by the Small-Molecule Activator ICA-105574

Outward current conducted by human ether-à-go-go-related gene type 1 (hERG1) K+ channels is important for action potential repolarization in the human ventricle. Rapid, voltage-dependent inactivation greatly reduces outward currents conducted by hERG1 channels and involves conformational changes in the ion selectivity filter (SF). Recently, compounds have been found that activate hERG1 channel function by modulating gating mechanisms such as reducing inactivation. Such activating compounds could represent a novel approach to prevent arrhythmias associated with prolonged ventricular repolarization associated with inherited or acquired long QT syndrome. ICA-105574 (ICA), a 3-nitro-n-(4-phenoxyphenyl) benzamide derivative activates hERG1 by strongly attenuating pore-type inactivation. We previously mapped the putative binding site for ICA to a hydrophobic pocket located between two adjacent subunits. Here, we used the recently reported cryoelectron microscopy structures of hERG1 to elucidate the structural mechanisms by which ICA influences the stability of the SF. By combining molecular dynamics simulations, voltage-clamp electrophysiology, and the synthesis of novel ICA derivatives, we provide atomistic insights into SF dynamics and propose a structural link between the SF and S6 segments. Further, our study highlights the importance of the nitro moiety, at the meta position of the benzamide ring, for the activity of ICA and reveals that the (bio)isosteric substitution of this side chain can switch the activity to weak inhibitors. Our findings indicate that ICA increases the stability of the SF to attenuate channel inactivation, and this action requires a fine-tuned compound geometry.

Neutralization of a bothropic PLA2-like protein by caftaric acid, a novel potent inhibitor of ophidian myotoxicity

Envenoming by snakebite is an important global health issue that has received little attention, leading the World Health Organization to naming it as neglected tropical disease. Several snakebites present serious local symptoms manifested on victims that may not be efficiently neutralized by serum therapy. Phospholipase A₂-like (PLA₂-like) toxins are present in Viperidae venoms and are responsible for local myotoxic activity. Herein, we investigated the association between BthTX-I toxin and caftaric acid (CFT), a molecule present in plants. CFT neutralized neuromuscular blocking and muscle-damaging activities promoted by BthTX-I. Calorimetric and light-scattering assays demonstrated that CFT inhibitor interacted with dimeric BthTX-I. Bioinformatics simulations indicated that CFT inhibitor binds to the toxin's hydrophobic channel (HCh). According to the current myotoxic mechanism, three different regions of PLA₂-like toxins have specific tasks: protein allosteric activation (HCh), membrane dockage (MDoS), and membrane rupture (MDiS). We propose CFT inhibitor interferes with the allosteric activation, which is related to the conformation change leading to the exposure/alignment of MDoS/MDiS region. This is the first report of a PLA₂-like toxin fully inhibited by a compound that interacts only with its HCh region. Thus, CFT is a novel candidate to complement serum therapy and improve the treatment of snakebite.

Advances in the Treatment of Explicit Water Molecules in Docking and Binding Free Energy Calculations

Background:

The inclusion of direct effects mediated by water during the ligandreceptor recognition is a hot-topic of modern computational chemistry applied to drug discovery and development. Docking or virtual screening with explicit hydration is still debatable, despite the successful cases that have been presented in the last years. Indeed, how to select the water molecules that will be included in the docking process or how the included waters should be treated remain open questions.

Objective:

In this review, we will discuss some of the most recent methods that can be used in computational drug discovery and drug development when the effect of a single water, or of a small network of interacting waters, needs to be explicitly considered.

Results:

Here, we analyse the software to aid the selection, or to predict the position, of water molecules that are going to be explicitly considered in later docking studies. We also present software and protocols able to efficiently treat flexible water molecules during docking, including examples of applications. Finally, we discuss methods based on molecular dynamics simulations that can be used to integrate docking studies or to reliably and efficiently compute binding energies of ligands in presence of interfacial or bridging water molecules.

Conclusions:

Software applications aiding the design of new drugs that exploit water molecules, either as displaceable residues or as bridges to the receptor, are constantly being developed. Although further validation is needed, workflows that explicitly consider water will probably become a standard for computational drug discovery soon.

Anti-tick effect and cholinesterase inhibition caused by Prosopis juliflora alkaloids: in vitro and in silico studies

Abstract We investigated the in vitro acaricide activity of the methanolic extract (ME) and alkaloid-rich fraction (AF) of Prosopis juliflora on Rhipicephalus microplus and correlated this effect with acetylcholinesterase (AChE) inhibition. The acaricide activity was evaluated using adult and larval immersion tests. Also, we studied the possible interaction mechanism of the major alkaloids present in this fraction via molecular docking at the active site of R. microplus AChE1 (RmAChE1). Higher reproductive inhibitory activity of the AF was recorded, with effective concentration (EC50) four times lower than that of the ME (31.6 versus 121 mg/mL). The AF caused mortality of tick larvae, with lethal concentration 50% (LC50) of 13.8 mg/mL. Both ME and AF were seen to have anticholinesterase activity on AChE of R. microplus larvae, while AF was more active with half-maximal inhibitory concentration (IC50) of 0.041 mg/mL. The LC-MS/MS analyses on the AF led to identification of three alkaloids: prosopine (1), juliprosinine (2) and juliprosopine (3). The molecular docking studies revealed that these alkaloids had interactions at the active site of the RmAChE1, mainly relating to hydrogen bonds and cation-pi interactions. We concluded that the alkaloids of P. juliflora showed acaricide activity on R. microplus and acted through an anticholinesterase mechanism.

Small-Molecule Intervention At The Dimerization Interface Of Survivin By Novel Rigidized Scaffolds

Introduction

Survivin is a nodal protein involved in several cellular pathways. It is a member of the IAP family and an integral component of the chromosomal passenger complex, where it binds to borealin and INCENP through its dimerization interface. By targeting survivin with a small molecule at its dimerization interface, inhibition of the proliferation of cancer cells has been suggested. With Abbott 8, a small-molecule dimerization inhibitor has been recently reported. The structure–activity relationship of this series of inhibitors implied that the middle pyridin-2(1H)-one ring did not tolerate modifications of any kind.

Methods

Based on the synthetic strategy of Abbott 8 using multicomponent reactions, we synthesized a series of small molecules bearing a novel rigidized core scaffold. This rigidization strategy was accomplished by integrating the pyridin-2(1H)-one and its 6-phenyl substituent into a tricyclic structure, linking position 5 of pyridin-2(1H)-one to the phenyl substituent by rings of different sizes. The new scaffolds were designed based on in silico molecular dynamics of survivin.

Results

Binding of these rigidized scaffolds to the recombinant L54M mutant of survivin was evaluated, revealing affinities in the low micromolar range.

Conclusion

This easily accessible, new class of survivin-dimerization modulators is an interesting starting point for further lead optimization.

Screening of different drug design tools to predict the mode of action of steroidal derivatives as anti-cancer agents

There is a pressing need to discover and develop novel drugs against cancer. With the new era of bioinformatics, which integrates different aspects, drug development has been tremendously improved. Recently, extensive research was directed towards the rational modification of steroid molecules against different disease especially cancer. Moreover, heterocyclic steroid derivatives have shown a lot of different biological activities such as antimicrobial, anti-inflammatory, and anti-cancer activities. Molecular docking methods can be used to explore how the steroid derivatives conformations can adopt within the binding sites of specific macromolecular targets involved in cancer progression. We conducted this study to investigate the accuracy of different molecular docking calculations using different steroidal molecular targets, and to define the most accurate one to study the mode of action of steroid derivatives as potential anti-cancer drugs. Our results revealed that the Dock6, PLANTS, AutoDock, GLIDE (SP and XP), and GOLD (ASP, Chemscore, and PLP) software were able to maintain the binding mode of the co-crystallized ligands inside their proteins by achieving RMSD values lower than two. Moreover, molecular docking study revealed that compound 4, and 5 are promising steroidal derivatives as anti-cancer drugs. Further on, the cytotoxic activity of the selected steroidal derivatives were tested against leukemia cell line using MTT assay. The results revealed that compound 4, and 5 were potential cytotoxic agents against THP-1 cells (IC₅₀s were 44.67 µM, and 46.77 µM, respectively), these results are in agreement with the molecular docking study.