Accounting for ligand-bound metal ions in docking small molecules on adenylyl cyclase toxins

Design, Synthesis, Biological Evaluation, and Molecular Docking Study of 4,6-Dimethyl-5-aryl/alkyl-2-[2-hydroxy-3-(4-substituted-1-piperazinyl)propyl]pyrrolo[3,4-c]pyrrole-1,3(2H,5H)-diones as Anti-Inflammatory Agents with Dual Inhibition of COX and LOX

In the present study, we characterize the biological activity of a newly designed and synthesized series of 15 compounds 2-[2-hydroxy-3-(4-substituted-1-piperazinyl)propyl] derivatives of pyrrolo[3,4-c]pyrrole 3a-3o. The compounds were obtained with good yields of pyrrolo[3,4-c]pyrrole scaffold 2a-2c with secondary amines in C₂H₅OH. The chemical structures of the compounds were characterized by ¹H-NMR, ¹³C-NMR, FT-IR, and MS. All the new compounds were investigated for their potencies to inhibit the activity of three enzymes, i.e., COX-1, COX-2, and LOX, by a colorimetric inhibitor screening assay. In order to analyze the structural basis of interactions between the ligands and cyclooxygenase/lipooxygenase, experimental data were supported by the results of molecular docking simulations. The data indicate that all of the tested compounds influence the activity of COX-1, COX-2, and LOX.

Assessing the effects of β-triketone herbicides on HPPD from environmental bacteria using a combination of in silico and microbiological approaches

4-hydroxyphenylpyruvate dioxygenase (HPPD) is the molecular target of β-triketone herbicides in plants. This enzyme, involved in the tyrosine pathway, is also present in a wide range of living organisms, including microorganisms. Previous studies, focusing on a few strains and using high herbicide concentrations, showed that β-triketones are able to inhibit microbial HPPD. Here, we measured the effect of agronomical doses of β-triketone herbicides on soil bacterial strains. The HPPD activity of six bacterial strains was tested with 1× or 10× the recommended field dose of the herbicide sulcotrione. The selected strains were tested with 0.01× to 15× the recommended field dose of sulcotrione, mesotrione, and tembotrione. Molecular docking was also used to measure and model the binding mode of the three herbicides with the different bacterial HPPD. Our results show that responses to herbicides are strain-dependent with Pseudomonas fluorescens F113 HPPD activity not inhibited by any of the herbicide tested, when all three β-triketone herbicides inhibited HPPD in Bacillus cereus ATCC14579 and Shewanella oneidensis MR-1. These responses are also molecule-dependent with tembotrione harboring the strongest inhibitory effect. Molecular docking also reveals different binding potentials. This is the first time that the inhibitory effect of β-triketone herbicides is tested on environmental strains at agronomical doses, showing a potential effect of these molecules on the HPPD enzymatic activity of non-target microorganisms. The whole-cell assay developed in this study, coupled with molecular docking analysis, appears as an interesting way to have a first idea of the effect of herbicides on microbial communities, prior to setting up microcosm or even field experiments. This methodology could then largely be applied to other family of pesticides also targeting an enzyme present in microorganisms.

Multiple Poses and Thermodynamics of Ligands Targeting Protein Surfaces: The Case of Furosemide Binding to mitoNEET in Aqueous Solution

Human NEET proteins, such as NAF-1 and mitoNEET, are homodimeric, redox iron-sulfur proteins characterized by triple cysteine and one histidine-coordinated [2Fe-2S] cluster. They exist in an oxidized and reduced state. Abnormal release of the cluster is implicated in a variety of diseases, including cancer and neurodegeneration. The computer-aided and structure-based design of ligands affecting cluster release is of paramount importance from a pharmaceutical perspective. Unfortunately, experimental structural information so far is limited to only one ligand/protein complex. This is the X-ray structure of furosemide bound to oxidized mitoNEET. Here we employ an enhanced sampling approach, Localized Volume-based Metadynamics, developed by some of us, to identify binding poses of furosemide to human mitoNEET protein in solution. The binding modes show a high variability within the same shallow binding pocket on the protein surface identified in the X-ray structure. Among the different binding conformations, one of them is in agreement with the crystal structure’s one. This conformation might have been overstabilized in the latter because of the presence of crystal packing interactions, absent in solution. The calculated binding affinity is compatible with experimental data. Our protocol can be used in a straightforward manner in drug design campaigns targeting this pharmaceutically important family of proteins.

Role of an orphan substrate-binding protein MhuP in transient heme transfer in Mycobacterium tuberculosis

The redox property of iron makes it an essential cofactor for numerous enzymes involved in various metabolic processes. In vertebrates, iron is attached to either heme molecules or with other circulatory proteins, making its accessibility restricted for bacterial pathogens residing inside the host. Due to this importance, there is always an ongoing battle between the host system and pathogens, known as nutritional immunity. To capture the bound iron from the human hosts, intracellular pathogens like Mycobacterium tuberculosis secrete siderophore molecules which are ultimately uptaken by versatile transport machinery such as ATP-binding cassette (ABC) transporters. Earlier reports have suggested the presence of a heme uptake protein MhuP (ORF id: Rv0265c) in M. tuberculosis, which transiently transfers the bound iron to the protein DppA for further heme transport by utilizing its cognate transport machinery (DppBCD). In the present study, we report the crystal structure of MhuP. The binding experiments of heme with MhuP suggest its specific nature. Molecular docking studies confirm the binding of the protein MhuP with heme as well as to the protein DppA. Thus, the results indicate the binding of heme to MhuP and its probable transient transport via the DppABCD transport system in M. tuberculosis.

Zinc thiotropolone combinations as inhibitors of the SARS-CoV-2 main protease

Numerous organic molecules are known to inhibit the main protease of SARS-CoV-2, (SC2M^pro), a key component in viral replication of the 2019 novel coronavirus. We explore the hypothesis that zinc ions, long used as a medicinal supplement and known to support immune function, bind to the SC2M^pro enzyme in combination with lipophilic tropolone and thiotropolone ligands, L, block substrate docking, and inhibit function. This study combines synthetic inorganic chemistry, in vitro protease activity assays, and computational modeling. While the ligands themselves have half maximal inhibition concentrations, IC₅₀, for SC2M^pro in the 8-34 μM range, the IC₅₀ values are ca. 100 nM for Zn(NO₃)₂ which are further enhanced in Zn-L combinations (59-97 nM). Isolation of the Zn(L)₂ binary complexes and characterization of their ability to undergo ligand displacement is the basis for computational modeling of the chemical features of the enzyme inhibition. Blind docking onto the SC2M^pro enzyme surface using a modified Autodock4 protocol found preferential binding into the active site pocket. Such Zn-L combinations orient so as to permit dative bonding of Zn(L)⁺ to basic active site residues.

In Vitro and In Silico Evaluation of New 1,3,4-Oxadiazole Derivatives of Pyrrolo[3,4-d]pyridazinone as Promising Cyclooxygenase Inhibitors

Since long-term use of classic NSAIDs can cause severe side effects related mainly to the gastroduodenal tract, discovery of novel cyclooxygenase inhibitors with a safe gastric profile still remains a crucial challenge. Based on the most recent literature data and previous own studies, we decided to modify the structure of already reported 1,3,4-oxadiazole based derivatives of pyrrolo[3,4-d]pyridazinone in order to obtain effective COX inhibitors. Herein we present the synthesis, biological evaluation and molecular docking studies of 12 novel compounds with disubstituted arylpiperazine pharmacophore linked in a different way with 1,3,4-oxadiazole ring. None of the obtained molecules show cytotoxicity on NHDF and THP-1 cell lines and, therefore, all were qualified for further investigation. In vitro cyclooxygenase inhibition assay revealed almost equal activity of new derivatives towards both COX-1 and COX-2 isoenzymes. Moreover, all compounds inhibit COX-2 isoform better than Meloxicam which was used as reference. Anti-inflammatory activity was confirmed in biological assays according to which title molecules are able to reduce induced inflammation within cells. Molecular docking studies were performed to describe the binding mode of new structures to cyclooxygenase. Investigated derivatives take place in the active site of COX, very similar to Meloxicam. For some compounds, promising druglikeness was calculated using in silico predictions.

Evaluation of Interactions of Selected Olivacine Derivatives with DNA and Topoisomerase II

Olivacine and ellipticine are model anticancer drugs acting as topoisomerase II inhibitors. Here, we present investigations performed on four olivacine derivatives in light of their antitumor activity. The aim of this study was to identify the best antitumor compound among the four tested olivacine derivatives. The study was performed using CCRF/CEM and MCF-7 cell lines. Comet assay, polarography, inhibition of topoisomerase II activity, histone acetylation, and molecular docking studies were performed. Each tested compound displayed interaction with DNA and topoisomerase II, but did not cause histone acetylation. Compound 2 (9-methoxy-5,6-dimethyl-1-({[1-hydroxy-2-(hydroxymethyl)butan-2-yl]amino}methyl)-6H-pyrido[4,3-b]carbazole) was found to be the best candidate as an anticancer drug because it had the highest affinity for topoisomerase II and caused the least genotoxic damage in cells.

Dinitrosyl iron complexes (DNICs) as inhibitors of the SARS-CoV-2 main protease

By repurposing DNICs designed for other medicinal purposes, the possibility of protease inhibition was investigated in silico using AutoDock 4.2.6 (AD4) and in vitro via a FRET protease assay. AD4 was validated as a predictive computational tool for coordinatively unsaturated DNIC binding using the only known crystal structure of a protein-bound DNIC, PDB- (calculation RMSD = 1.77). From the in silico data the dimeric DNICs TGTA-RRE, [(μ-S-TGTA)Fe(NO)₂]₂ (TGTA = 1-thio-β-d-glucose tetraacetate) and TG-RRE, [(μ-S-TG)Fe(NO)₂]₂ (TG = 1-thio-β-d-glucose) were identified as promising leads for inhibition via coordinative inhibition at Cys-145 of the SARS-CoV-2 Main Protease (SC2M^pro). In vitro studies indicate inhibition of protease activity upon DNIC treatment, with an IC₅₀ of 38 ± 2 μM for TGTA-RRE and 33 ± 2 μM for TG-RRE. This study presents a simple computational method for predicting DNIC-protein interactions; the in vitro study is consistent with in silico leads.

Novel Pyrimidine Derivatives as Potential Anticancer Agents: Synthesis, Biological Evaluation and Molecular Docking Study

In the present paper, new pyrimidine derivatives were designed, synthesized and analyzed in terms of their anticancer properties. The tested compounds were evaluated in vitro for their antitumor activity. The cytotoxic effect on normal human dermal fibroblasts (NHDF) was also determined. According to the results, all the tested compounds exhibited inhibitory activity on the proliferation of all lines of cancer cells (colon adenocarcinoma (LoVo), resistant colon adenocarcinoma (LoVo/DX), breast cancer (MCF-7), lung cancer (A549), cervical cancer (HeLa), human leukemic lymphoblasts (CCRF-CEM) and human monocytic (THP-1)). In particular, their feature stronger influence on the activity of P-glycoprotein of cell cultures resistant to doxorubicin than doxorubicin. Tested compounds have more lipophilic character than doxorubicin, which determines their affinity for the molecular target and passive transport through biological membranes. Moreover, the inhibitory potential against topoisomerase II and DNA intercalating properties of synthesized compounds were analyzed via molecular docking.

Steady-State Kinetics of Enzyme-Catalyzed Hydrolysis of Echothiophate, a P-S Bonded Organophosphorus as Monitored by Spectrofluorimetry

Enzyme-catalyzed hydrolysis of echothiophate, a P–S bonded organophosphorus (OP) model, was spectrofluorimetrically monitored, using Calbiochem Probe IV as the thiol reagent. OP hydrolases were: the G117H mutant of human butyrylcholinesterase capable of hydrolyzing OPs, and a multiple mutant of Brevundimonas diminuta phosphotriesterase, GG1, designed to hydrolyze a large spectrum of OPs at high rate, including V agents. Molecular modeling of interaction between Probe IV and OP hydrolases (G117H butyrylcholinesterase, GG1, wild types of Brevundimonas diminuta and Sulfolobus solfataricus phosphotriesterases, and human paraoxonase-1) was performed. The high sensitivity of the method allowed steady-state kinetic analysis of echothiophate hydrolysis by highly purified G117H butyrylcholinesterase concentration as low as 0.85 nM. Hydrolysis was michaelian with K_m = 0.20 ± 0.03 mM and k_cat = 5.4 ± 1.6 min⁻¹. The GG1 phosphotriesterase hydrolyzed echothiophate with a high efficiency (K_m = 2.6 ± 0.2 mM; k_cat = 53400 min⁻¹). With a k_cat/K_m = (2.6 ± 1.6) × 10⁷ M⁻¹min⁻¹, GG1 fulfills the required condition of potential catalytic bioscavengers. quantum mechanics/molecular mechanics (QM/MM) and molecular docking indicate that Probe IV does not interact significantly with the selected phosphotriesterases. Moreover, results on G117H mutant show that Probe IV does not inhibit butyrylcholinesterase. Therefore, Probe IV can be recommended for monitoring hydrolysis of P–S bonded OPs by thiol-free OP hydrolases.

In silico repositioning of etidronate as a potential inhibitor of the Trypanosoma cruzi enolase

Enolase is a glycolytic enzyme that catalyzes the interconversion between 2-phosphoglycerate and phosphoenolpyruvate. In trypanosomatids, enolase was proposed as a key enzyme after in silico and in vivo analysis and it was validated as a protein essential for the survival of the parasite. Therefore, enolase constitutes an interesting enzyme target for the identification of drugs against Chagas disease. In this work, a combined virtual screening strategy was implemented, employing similarity virtual screening, molecular docking, and molecular dynamics. First, two known enolase inhibitors and the enzyme substrates were used as queries for the similarity screening on the Sweetlead database using five different algorithms. Compounds retrieved in the top 10 of at least three search algorithms were selected for further analysis, resulting in six compounds of medical use (etidronate, pamidronate, fosfomycin, acetohydroxamate, triclofos, and aminohydroxybutyrate). Molecular docking simulations and pose re-scoring predicted that binding with acetohydroxamate and triclofos would be weak, while fosfomycin and aminohydroxybutyrate predicted binding is experimentally implausible. Docking poses obtained for etidronate, pamidronate, and PEP were used for molecular dynamics calculations to describe their mode of binding. From the obtained results, we propose etidronate as a potential TcENO inhibitor and describe molecular motifs to be taken into account in the repurposing or design of drugs targeting this enzyme active site.

Structure and catalytic mechanistic insight into Enterobacter aerogenes acetolactate decarboxylase

α-Acetolactate decarboxylase (ALDC) catalyses α-acetolactate into acetoin (3-hydroxy-2-butanone, AC) and is considered to be the rate-limiting enzyme in the synthesis of 2,3-butanediol. In this work, the enzymatic activity of ALDC from Enterobacter aerogenes ALDC (E.a.-ALDC) was fully characterized with enzyme kinetics, indicating a K_m of 14.83 ± 0.87 mM and a k_cat of 0.81 ± 0.09 s^-1. However, compared with the activities of ALDCs reported from other bacteria, the activity of E.a.-ALDC was determined to present a relatively lower value of 849.08 ± 35.21 U/mg. The enzyme showed maximum activity at pH 5.5. In addition, the activity of E.a.-ALDC was promoted by Mg²⁺. The crystal structure of E.a.-ALDC firstly solved by X-ray crystallography at resolution of 2.4 Å revealed a chelated zinc ion with conserved His199, His201, His212, Glu70 and Glu259. In the active center, the conservative Arg150 was particularly proven to deviate from the zinc ion of the active centre, by adopting a flexible conformational change, resulting in a weak interaction network of the enzyme and the substrate. Further in silico docking of E.a.-ALDC with two enantiomers, (R)-acetolactate and (S)-acetolactate, unaltered the interaction network of E.a.-ALDC from the apo structure, which confirmed the weakened role of Arg150 in the catalytic properties of E.a.-ALDC. Our results reveal a unique structure-function relationship of acetolactate decarboxylase and provide a fundamental basis for the enzymatic synthesis of acetoin.

Designating ligand specificities to metal uptake ABC transporters in Thermus thermophilus HB8

Acquisition of different metal ions by metal uptake ABC transporters of Thermus thermophilus HB8 for accomplishing its various cellular functions.

Screening of novel histone deacetylase 7 inhibitors through molecular docking followed by a combination of molecular dynamics simulations and ligand-based approach

Histone acetylation/deacetylation is a key mechanism for transcription regulation which plays an important role in control of gene expression, tissue growth, and development. In particular, histone deacetylase 7 (HDAC7), a member of class IIa HDACs, is crucial to maintain cell homeostasis, and HDAC7 has emerged as a new target for cancer therapy. In this study, molecular docking was applied to screen candidate inhibitors and 21 compounds were found. Following the 50 ns molecular dynamics simulations and binding free energy calculation, ZINC00156160, ZINC01703144, ZINC04293665, and ZINC13900201 were identified as potential HDAC7 inhibitors, which would provide a sound starting point for further studies involving molecular modeling coupled with biochemical experiments. Meanwhile, similarity computation and substructure search were combined, and then we found that compounds sharing common backbone "CC(=O)N[C@@H](CSc1ccccc1)C(=O)O" could be efficient to inhibit the bioactivity of HDAC7. Then comparative molecular similarity indices analysis (CoMSIA) techniques were implemented to investigate the relationship between properties of the substituent group and bioactivities of small molecules. The CoMSIA model exhibited powerful predictivity, with satisfactory statistical parameters such as q² of 0.659, R² of 0.952, and F of 268.448. Contour maps of the CoMSIA model gave insight into the feature requirements of the common backbone for the HDAC7 inhibitory activity. Finally, details of designing novel HDAC7 inhibitors were confirmed by a combination of receptor-based docking and ligand-based structure-activity relationship. Communicated by Ramaswamy H. Sarma.

Inhibiting effects of common trivalent metal ions on transmembrane-type 2 matrix metalloproteinase

Transmembrane-type 2 matrix metalloproteinase (MT2-MMP) degrades connective extracellular matrix between cells and enables tumor cells to migrate and metastasize, making this substance a potential therapeutic target in various diseases. In this work, the interactions between MT2-MMP and common trivalent metal ions, including aluminum (Al³⁺) and ferrum (Fe³⁺) ions, were investigated. Enzymatic detection revealed that Al³⁺ and Fe³⁺ strongly inhibited the MT2-MMP. Fluorescence spectrography elucidated a static quenching interaction between the negatively charged amino acids on MT2-MMP and the inhibitory trivalent metal ions, indicating that a stable complex was formed between MT2-MMP and metal ions. In addition, fluorescence data and molecular modeling analysis of the binding characteristics revealed that one trivalent metal ion bound with a protein in the stable complex formation process. The potential inhibitory effect of Al³⁺ on MT2-MMP was further examined in an MT2-MMP-overexpressed cell line, HT1080, by using flow cytometry. As a result, Al³⁺ can promote HT1080 cell apoptosis in a micromolar concentration-dependent manner. This work illustrated that common trivalent metal ions can potentially inhibit MT2-MMP-related tumors.

The structural model of Zika virus RNA-dependent RNA polymerase in complex with RNA for rational design of novel nucleotide inhibitors

Zika virus is a global health threat due to significantly elevated risk of fetus malformations in infected pregnant women. Currently, neither an effective therapy nor a prophylactic vaccination is available for clinical use, desperately necessitating novel therapeutics and approaches to obtain them. Here, we present a structural model of the Zika virus RNA-dependent RNA polymerase (ZIKV RdRp) in complex with template and nascent RNAs, Mg²⁺ ions and accessing nucleoside triphosphate. The model allowed for docking studies aimed at effective pre-screening of potential inhibitors of ZIKV RdRp. Applicability of the structural model for docking studies was illustrated with the NITD008 artificial nucleotide that is known to effectively inhibit the function of the ZIKV RdRp. The ZIKV RdRp – RNA structural model is provided for all possible variations of the nascent RNA bases pairs to enhance its general utility in docking and modelling experiments. The developed model makes the rational design of novel nucleosides and nucleotide analogues feasible and thus provides a solid platform for the development of advanced antiviral therapy.

Crystal structure and substrate specificity of ExoY, a unique T3SS mediated secreted nucleotidyl cyclase toxin from Pseudomonas aeruginosa

BACKGROUND

The nucleotidyl cyclase toxin ExoY is an important virulence determinant of Pseudomonas aeruginosa that causes severe acute and chronic infections in immune-compromised individuals. Additionally, this unique T3SS effector shows a striking preference for cUMP, a newly identified non-canonical secondary messenger. Thereby, ExoY is also considered as a potential tool to study unexplored cUMP signaling pathways.

METHODS

The crystal structure of ExoY was determined at 2.2 Å resolutions by in-situ proteolysis assisted crystallization and Rosetta-molecular replacement method. Additionally, isothermal calorimetric (ITC) and molecular dynamic (MD) simulation studies were also carried out to gain molecular insights into its substrate specificity and catalysis.

RESULTS AND CONCLUSION

ExoY is a partially unfolded protein with higher propensity to form soluble higher-order oligomers. However, with meticulous attempts of removing of disordered regions by proteases, the recalcitrant ExoY could be successfully crystallized. The crystal structure of ExoY revealed similar overall structural fold present in other anthrax toxA family of nucleotidyl cyclases, with two-to-three distinctly conserved regions conferring specificity to eukaryotic binding partner. The in-vitro catalytic preference of ExoY is in the following order: cGMP > cUMP > cAMP > cCMP. The substrate specificity of ExoY mainly depends on its ability to bind NTP in proper geometrical orientations. ExoY also seems to prefer one-metal-ion dependent catalysis than two-metal-ion dependent catalysis.

GENERAL SIGNIFICANCE

Our results provide much needed structural insight on ExoY, an important virulence determinant of Pseudomonas aeruginosa and an exciting tool to study non-canonical cNMP signaling pathways.

ACCESSION NUMBERS

The structure factors and coordinate files have been deposited in the Protein Data Bank with accession number 5XNW.

Homology modeling and docking studies of a Δ9-fatty acid desaturase from a Cold-tolerant Pseudomonas sp. AMS8

Membrane-bound fatty acid desaturases perform oxygenated desaturation reactions to insert double bonds within fatty acyl chains in regioselective and stereoselective manners. The Δ9-fatty acid desaturase strictly creates the first double bond between C9 and 10 positions of most saturated substrates. As the three-dimensional structures of the bacterial membrane fatty acid desaturases are not available, relevant information about the enzymes are derived from their amino acid sequences, site-directed mutagenesis and domain swapping in similar membrane-bound desaturases. The cold-tolerant Pseudomonas sp. AMS8 was found to produce high amount of monounsaturated fatty acids at low temperature. Subsequently, an active Δ9-fatty acid desaturase was isolated and functionally expressed in Escherichia coli. In this paper we report homology modeling and docking studies of a Δ9-fatty acid desaturase from a Cold-tolerant Pseudomonas sp. AMS8 for the first time to the best of our knowledge. Three dimensional structure of the enzyme was built using MODELLER version 9.18 using a suitable template. The protein model contained the three conserved-histidine residues typical for all membrane-bound desaturase catalytic activity. The structure was subjected to energy minimization and checked for correctness using Ramachandran plots and ERRAT, which showed a good quality model of 91.6 and 65.0%, respectively. The protein model was used to preform MD simulation and docking of palmitic acid using CHARMM36 force field in GROMACS Version 5 and Autodock tool Version 4.2, respectively. The docking simulation with the lowest binding energy, -6.8 kcal/mol had a number of residues in close contact with the docked palmitic acid namely, Ile26, Tyr95, Val179, Gly180, Pro64, Glu203, His34, His206, His71, Arg182, Thr85, Lys98 and His177. Interestingly, among the binding residues are His34, His71 and His206 from the first, second, and third conserved histidine motif, respectively, which constitute the active site of the enzyme. The results obtained are in compliance with the in vivo activity of the Δ9-fatty acid desaturase on the membrane phospholipids.

Studies of the Interaction of Influenza Virus RNA Polymerase PAN with Endonuclease Inhibitors

Influenza virus is a major causative agent of respiratory viral infections, and RNA polymerase catalyzes its replication and transcription activities in infected cell nuclei. Since it is highly conserved in all virus strains, RNA polymerase becomes a key target of anti-influenza virus agents. Although experimental studies have revealed the good inhibitory activity of endonuclease inhibitors to RNA polymerase, the mechanism is still unclear. In this study, the docking and molecular dynamics simulations have been performed to explore the interaction of three kinds of endonuclease inhibitors with the subunit (PA_N) of RNA polymerase. Our calculations indicate that all these endonuclease inhibitors can bind to the binding pocket of PA_N, in which the electronegative oxygen atoms of the inhibitors form a chelated structure with the two Mn²⁺ cations of the active center. The most important interaction between these inhibitors and PA_N is electrostatic interaction. The electron density of the chelate oxygen atoms determines the magnitude of the electrostatic energy, and the chelated structure and orientation of inhibitors depend largely on the distance between the chelate oxygen atoms.

Captopril/enalapril inhibit promiscuous esterase activity of carbonic anhydrase at micromolar concentrations: An in vitro study

The inhibitory activity of captopril, a thiol-containing competitive inhibitor of the angiotensin-converting enzyme, ACE, against esterase activity of carbonic anhydrase, CA was investigated. This small molecule, as well as enalapril, was selected in order to represents both thiol and carboxylate, as two well-known metal binding functional groups of metalloprotein inhibitors. Since captopril, has also been observed to inhibit other metalloenzymes such as tyrosinase and metallo-beta lactamase through binding to the catalytic metal ions and regarding CA as a zinc-containing metallo-enzyme, in the current study, we set out to determine whether captopril/enalapril inhibit CA esterase activity of the purified human CA II or not? Then, we revealed the inhibitors' potencies (IC₅₀, K_i and K_diss values) and also mode of inhibition. Our results also showed that enalapril is more potent CA inhibitor than captopril. Since enalapril represents no sulfhydryl moiety, thus carboxylate group may have a determinant role in inhibiting of CA esterase activity, the conclusion confirmed by molecular docking studies. Additionally, since CA inhibitory potencies of captopril/enalapril were much lower than those of classic sulfonamide drugs, the findings of the current study may explain why these drugs exhibit no effective CA inhibition at the concentrations reached in vivo and also may shed light on the way of generating new class of inhibitors that will discriminately inhibit various CA isoforms.

Computational modeling and molecular dynamics simulations of mammalian cytoplasmic tyrosyl-tRNA synthetase and its complexes with substrates

Cytoplasmic tyrosyl-tRNA synthetase (TyrRS) is one of the key enzymes of protein biosynthesis. TyrRSs of pathogenic organisms have gained attention as potential targets for drug development. Identifying structural differences between various TyrRSs will facilitate the development of specific inhibitors for the TyrRSs of pathogenic organisms. However, there is a deficiency in structural data for mammalian cytoplasmic TyrRS in complexes with substrates. In this work, we constructed spatial structure of full-length Bos taurus TyrRS (BtTyrRS) and its complexes with substrates using the set of computational modeling techniques. Special attention was paid to BtTyrRS complexes with substrates [L-tyrosine, K⁺ and ATP:Mg²⁺] and intermediate products [tyrosyl-adenylate (Tyr-AMP), K⁺ and PP_i:Mg²⁺] with the different catalytic loop conformations. In order to analyze their dynamical properties, we performed 100 ns of molecular dynamics (MD) simulations. MD simulations revealed new structural data concerning the tyrosine activation reaction in mammalian TyrRS. Formation of strong interaction between Lys154 and γ-phosphate suggests the additional role of CP1 insertion as an important factor for ATP binding. The presence of a potassium-binding pocket within the active site of mammalian TyrRS compensates the absence of the second lysine in the KMSKS motif. Our data provide new details concerning a role of K⁺ ions at different stages of the first step of the tyrosylation reaction, including the coordination of substrates and involvement in the PP_i releasing. The results of this work suggest that differences between ATP-binding sites of mammalian and bacterial TyrRSs are meaningful and could be exploited in the drug design.

Allosteric inhibitors of Coxsackie virus A24 RNA polymerase

Coxsackie virus A24 (CVA24), a causative agent of acute hemorrhagic conjunctivitis, is a prototype of enterovirus (EV) species C. The RNA polymerase (3D(pol)) of CVA24 can uridylylate the viral peptide linked to the genome (VPg) from distantly related EV and is thus, a good model for studying this reaction. Once UMP is bound, VPgpU primes RNA elongation. Structural and mutation data have identified a conserved binding surface for VPg on the RNA polymerase (3D(pol)), located about 20Å from the active site. Here, computational docking of over 60,000 small compounds was used to select those with the lowest (best) specific binding energies (BE) for this allosteric site. Compounds with varying structures and low BE were assayed for their effect on formation of VPgU by CVA24-3D(pol). Two compounds with the lowest specific BE for the site inhibited both uridylylation and formation of VPgpolyU at 10-20μM. These small molecules can be used to probe the role of this allosteric site in polymerase function, and may be the basis for novel antiviral compounds.

Inhibition of sorbitol dehydrogenase by nucleosides and nucleotides

Sorbitol dehydrogenase inhibitors have been found to prevent, or alleviate, various secondary complications of diabetes mellitus. In the present study, the effects of nucleosides and nucleotides on the rate of sorbitol oxidation catalyzed by the sheep liver enzyme were studied by steady-state kinetics at pH 7.4. Various such compounds, including ATP and the 2'-deoxy-analogues of ATP, ADP and AMP, reversibly inhibit enzyme activity by formation of enzyme-coenzyme-inhibitor ternary complexes. In each case, no deviations from linearity were seen in the double-reciprocal plots using sorbitol or NAD(+) as the varied substrate and there was a linear relationship between inhibitor concentration and the observed inhibitory effects. Sorbitol was docked into a model of the sheep SDH-NAD(+) complex based upon the structure of the human SDH-NAD(+) holoenzyme. The resulting structure of the ternary complex of sheep SDH, NAD(+) and sorbitol (PMDB ID code PM 0078068) shows that the reactive C-2 hydroxyl group of sorbitol is oriented toward the 4'-position of the nicotinamide moiety of the coenzyme, and that the adjacent primary hydroxyl group of sorbitol interacts with the catalytic zinc. The results indicate that the ribose moiety of the inhibitor structures is an important determinant for the observed effects. Specifically, the 2'-position of the ribose ring exerts an effect with respect to inhibitor potency.

Pharmacophore selection and redesign of non-nucleotide inhibitors of anthrax edema factor

Antibiotic treatment may fail to protect individuals, if not started early enough, after infection with Bacillus anthracis, due to the continuing activity of toxins that the bacterium produces. Stable and easily stored inhibitors of the edema factor toxin (EF), an adenylyl cyclase, could save lives in the event of an outbreak, due to natural causes or a bioweapon attack. The toxin’s basic activity is to convert ATP to cAMP, and it is thus in principle a simple phosphatase, which means that many mammalian enzymes, including intracellular adenylcyclases, may have a similar activity. While nucleotide based inhibitors, similar to its natural substrate, ATP, were identified early, these compounds had low activity and specificity for EF. We used a combined structural and computational approach to choose small organic molecules in large, web-based compound libraries that would, based on docking scores, bind to residues within the substrate binding pocket of EF. A family of fluorenone-based inhibitors was identified that inhibited the release of cAMP from cells treated with EF. The lead inhibitor was also shown to inhibit the diarrhea caused by enterotoxigenic E. coli (ETEC) in a murine model, perhaps by serving as a quorum sensor. These inhibitors are now being tested for their ability to inhibit Anthrax infection in animal models and may have use against other pathogens that produce toxins similar to EF, such as Bordetella pertussis or Vibrio cholera.

Structure-based redesign of an edema toxin inhibitor

Edema factor (EF) toxin of Bacillus anthracis (NIAID category A), and several other toxins from NIAID category B Biodefense target bacteria are adenylyl cyclases or adenylyl cyclase agonists that catalyze the conversion of ATP to 3',5'-cyclic adenosine monophosphate (cAMP). We previously identified compound 1 (3-[(9-oxo-9H-fluorene-1-carbonyl)-amino]-benzoic acid), that inhibits EF activity in cultured mammalian cells, and reduces diarrhea caused by enterotoxigenic Escherichia coli (ETEC) at an oral dosage of 15μg/mouse. Here, molecular docking was used to predict improvements in potency and solubility of new derivatives of compound 1 in inhibiting edema toxin (ET)-catalyzed stimulation of cyclic AMP production in murine monocyte-macrophage cells (RAW 264.7). Structure-activity relationship (SAR) analysis of the bioassay results for 22 compounds indicated positions important for activity. Several derivatives demonstrated superior pharmacological properties compared to our initial lead compound, and are promising candidates to treat anthrax infections and diarrheal diseases induced by toxin-producing bacteria.

Molecular dynamics simulations of Aβ fibril interactions with β-sheet breaker peptides

Accumulation and aggregation of the 42-residue amyloid-β (Aβ) protein fragment, which originates from the cleavage of amyloid precursor protein by β and γ secretase, correlates with the pathology of Alzheimer's disease (AD). Possible therapies for AD include peptides based on the Aβ sequence, and recently identified small molecular weight compounds designed to mimic these, that interfere with the aggregation of Aβ and prevent its toxic effects on neuronal cells in culture. Here, we use molecular dynamics simulations to compare the mode of interaction of an active (LPFFD) and inactive (LHFFD) β-sheet breaker peptide with an Aβ fibril structure from solid-state NMR studies. We found that LHFFD had a weaker interaction with the fibril than the active peptide, LPFFD, from geometric and energetic considerations, as estimated by the MM/PBSA approach. Cluster analysis and computational alanine scanning identified important ligand-fibril contacts, including a possible difference in the effect of histidine on ligand-fibril π-stacking interactions, and the role of the proline residue in establishing contacts that compete with those essential for maintenance of the inter-monomer β-sheet structure of the fibril. Our results show that molecular dynamics simulations can be a useful way to classify the stability of docking sites. These mechanistic insights into the ability of LPFFD to reverse aggregation of toxic Aβ will guide the redesign of lead compounds, and aid in developing realistic therapies for AD and other diseases of protein aggregation.

ATP conformations and ion binding modes in the active site of anthrax edema factor: a computational analysis

The Edema Factor (EF), one of the virulence factors of anthrax, is an adenylyl cyclase that promotes the overproduction of cyclic-AMP (cAMP) from ATP, and therefore perturbs cell signaling. Crystallographic structures of EF bound to ATP analogs and reaction products, cyclic-AMP, and Pyrophosphate (PPi), revealed different substrate conformations and catalytic-cation binding modes, one or two cations being observed in the active site. To shed light into the biological significance of these crystallographic structures, the energetics, geometry, and dynamics of the active site are analyzed using molecular dynamics simulations. The ATP conformation observed in the one-metal-ion structure allows stronger interactions with the catalytic ion, and ATP is more restrained than in the structure containing two Mg(2+) ions. Therefore, we propose that the conformation observed in the one-ion crystal structure is a more probable starting point for the reaction. The simulations also suggest that a C3'-endo sugar pucker facilitates nucleophilic attack. Additionally, the two-cation binding mode restrains the mobility of the reaction products, and thus their tendency to dissociate.

Testing the efficacy and toxicity of adenylyl cyclase inhibitors against enteric pathogens using in vitro and in vivo models of infection

Enterotoxigenic Escherichia coli (ETEC) produces the ADP-ribosyltransferase toxin known as heat-labile enterotoxin (LT). In addition to the toxic effect of LT resulting in increases of cyclic AMP (cAMP) and disturbance of cellular metabolic processes, this toxin promotes bacterial adherence to intestinal epithelial cells (A. M. Johnson, R. S. Kaushik, D. H. Francis, J. M. Fleckenstein, and P. R. Hardwidge, J. Bacteriol. 191:178-186, 2009). Therefore, we hypothesized that the identification of a compound that inhibits the activity of the toxin would have a suppressive effect on the ETEC colonization capabilities. Using in vivo and in vitro approaches, we present evidence demonstrating that a fluorenone-based compound, DC5, which inhibits the accumulation of cAMP in intoxicated cultured cells, significantly decreases the colonization abilities of adenylyl cyclase toxin-producing bacteria, such as ETEC. These findings established that DC5 is a potent inhibitor both of toxin-induced cAMP accumulation and of ETEC adherence to epithelial cells. Thus, DC5 may be a promising compound for treatment of diarrhea caused by ETEC and other adenylyl cyclase toxin-producing bacteria.

Computational selection of inhibitors of Abeta aggregation and neuronal toxicity

Alzheimer's disease (AD) is characterized by the cerebral accumulation of misfolded and aggregated amyloid-beta protein (Abeta). Disease symptoms can be alleviated, in vitro and in vivo, by 'beta-sheet breaker' pentapeptides that reduce plaque load. However the peptide nature of these compounds, made them biologically unstable and unable to penetrate membranes with high efficiency. The main goal of this study was to use computational methods to identify small molecule mimetics with better drug-like properties. For this purpose, the docked conformations of the active peptides were used to identify compounds with similar activities. A series of related beta-sheet breaker peptides were docked to solid state NMR structures of a fibrillar form of Abeta. The lowest energy conformations of the active peptides were used to design three dimensional (3D)-pharmacophores, suitable for screening the NCI database with Unity. Small molecular weight compounds with physicochemical features and a conformation similar to the active peptides were selected, ranked by docking and biochemical parameters. Of 16 diverse compounds selected for experimental screening, 2 prevented and reversed Abeta aggregation at 2-3microM concentration, as measured by Thioflavin T (ThT) fluorescence and ELISA assays. They also prevented the toxic effects of aggregated Abeta on neuroblastoma cells. Their low molecular weight and aqueous solubility makes them promising lead compounds for treating AD.