Docking of 4-oxalocrotonate tautomerase substrates: implications for the catalytic mechanism

New morpholine-containing pyrimidinones act on α-adrenoceptors

Drugs that act on α-adrenoceptors may contain morpholine and pyrimidinone heterocycles. The aim of this study was to synthesize a series of pyrimidinones (S6a-e and S8) and characterize their α-adrenoceptor activity. Cytotoxicity assays (MTT and LDH) were performed in A7r5 and HUVECs. Concentration-effect curves to phenylephrine (Phe) were performed in rat aortic rings in the presence of compounds S6a-e and S8 or vehicle. Nitric oxide (NO) production and NO stable metabolic products, nitrite and nitrate, expressed as total nitrogen oxides (NOx) were assessed in HUVECs by confocal microscopy with the DAF-2DA probe and by the Griess reaction, respectively. Molecular docking simulations were performed using the 6a compound and α2A-adrenoceptor. In the evaluated conditions, the percentage of viable cells and the release of LDH were similar between control cells and cells exposed to the tested pyrimidinones. S6d, S6e, S8, and the positive control prazosin (but not S6a, S6b, and S6c) decreased Phe-induced contractions in endothelium-denuded aortic rings. S6a, S6b, and S6c decreased Phe-induced contractions in endothelium-intact aortic rings. The effect of S6a was abolished by L-NAME. NO production and NOx levels were inhibited in the presence of the α2 receptor antagonist yohimbine and the NOS inhibitor L-NAME. The 6a docking simulation estimated that the mean binding free energy of the compound was lower than the estimated value for yohimbine. These data suggest that S6d, S6e, and S8 may be α1-adrenoceptor antagonists while S6a acts as an agonist of α2-adrenoceptors.

Crucial role of carbon monoxide as a regulator of diarrhea caused by cholera toxin: Evidence of direct interaction with toxin

The present study evaluated the role of heme oxygenase 1 (HO-1)/carbon monoxide (CO) pathway in the cholera toxin-induced diarrhea and its possible action mechanism. The pharmacological modulation with CORM-2 (a CO donor) or Hemin (a HO-1 inducer) decreased the intestinal fluid secretion and Cl- efflux, altered by cholera toxin. In contrast, ZnPP (a HO-1 inhibitor) reversed the antisecretory effect of Hemin and potentiated cholera toxin-induced intestinal secretion. Moreover, CORM-2 also prevented the alteration of intestinal epithelial architecture and local vascular permeability promoted by cholera toxin. The intestinal absorption was not altered by any of the pharmacological modulators. Cholera toxin inoculation also increased HO-1 immunoreactivity and bilirubin levels, a possible protective physiological response. Finally, using fluorometric technique, ELISA assay and molecular docking simulations, we show evidence that CO directly interacts with cholera toxin, forming a complex that affects its binding to GM1 receptor, which help explain the antisecretory effect. Thus, CO is an essential molecule for protection against choleric diarrhea and suggests its use as a possible therapeutic tool.

A second generation of 1,2,4-oxadiazole derivatives with enhanced solubility for inhibition of 3-hydroxykynurenine transaminase (HKT) from Aedes aegypti

The most widely used method for the control of the Aedes aegypti mosquito population is the chemical control method. It represents a time- and cost-effective way to curb several diseases (e.g. dengue, Zika, chikungunya, yellow fever) through vector control. For this reason, the discovery of new compounds with a distinct mode of action from the available ones is essential in order to minimize the rise of insecticide resistance. Detoxification enzymes are an attractive target for the discovery of new insecticides. The kynurenine pathway is an important metabolic pathway, and it leads to the chemically stable xanthurenic acid, biosynthesized from 3-hydroxykynurenine, a precursor of reactive oxygen and nitrogen species, by the enzyme 3-hydroxykynurenine transaminase (HKT). Previously, we have reported the effectiveness of 1,2,4-oxadiazole derivatives acting as larvicides for A. aegypti and AeHKT inhibitors from in vitro and in silico studies. Here, we report the synthesis of new sodium 4-[3-(aryl)-1,2,4-oxadiazol-5-yl] propanoates and the cognate HKT-inhibitory activity. These new derivatives act as competitive inhibitors with IC₅₀ values in the range of 42 to 339 μM. We further performed molecular docking simulations and QSAR analysis for the previously synthesized sodium 4-[3-(aryl)-1,2,4-oxadiazol-5-yl] butanoates reported earlier by our group and the data produced herein. Most of the 1,2,4-oxadiazole derivatives, including the canonical compounds for both series, showed a similar binding mode with HKT. The binding occurs similarly to the co-crystallized inhibitor via anchoring to Arg³⁵⁶ and positioning of the aromatic ring and its substituents outwards at the entry of the active site. QSAR analysis was performed in search of more than 770 molecular descriptors to establish a relationship between the lowest energy conformations and the IC₅₀ values. The five best descriptors were selected to create and validate the model, which exhibited parameters that attested to its robustness and predictability. In summary, we observed that compounds with a para substitution and heavier groups (i.e. CF₃ and NO₂ substituents) had an enhanced HKT-inhibition profile. These compounds comprise a series described as AeHKT inhibitors via enzymatic inhibition experiments, opening the way to further the development of new substances with higher potency against HKT from Aedes aegypti.

Discovery of 1,2,4-oxadiazole derivatives as a novel class of noncompetitive inhibitors of 3-hydroxykynurenine transaminase (HKT) from Aedes aegypti

The mosquito Aedes aegypti is the vector of arboviruses such as Zika, Chikungunya, dengue and yellow fever. These infectious diseases have a major impact on public health. The unavailability of effective vaccines or drugs to prevent or treat most of these diseases makes vector control the main form of prevention. One strategy to promote mosquito population control is the use of synthetic insecticides to inhibit key enzymes in the metabolic pathway of these insects, particularly during larval stages. One of the main targets of the kynurenine detoxification pathway in mosquitoes is the enzyme 3-hydroxykynurenine transaminase (HKT), which catalyzes the conversion of 3-hydroxykynurenine (3-HK) into xanthurenic acid (XA). In this work, we report eleven newly synthesized oxadiazole derivatives and demonstrate that these compounds are potent noncompetitive inhibitors of HKT from Ae. aegypti. The present data provide direct evidence that HKT can be explored as a molecular target for the discovery of novel larvicides against Ae. aegypti. More importantly, it ensures that structural information derived from the HKT 3D-structure can be used to guide the development of more potent inhibitors.

Duffy binding-like 1α adhesin from Plasmodium falciparum recognizes ABH histo-blood group saccharide in a type specific manner

The ability of erythrocytes, infected by Plasmodium falciparum, to adhere to endothelial cells (cytoadherence) and to capture uninfected erythrocyte (rosetting) is the leading cause of death by severe malaria. Evidences link the binding of the adhesin Duffy Binding Like1-α (DBL1α) domain to the ABH histo-blood antigens with formation of rosettes. Inspired by this very close relationship between the disease susceptibility and individual blood type, here we investigate the structural requirements involved in the interaction of DBL1α with A, B and H histo-blood determinants and their subtypes. Our results evidence the high preference of DBL1α to A epitopes, in comparison to B and H epitopes. DBL1α interacts with ABH epitopes in subtype specific manner, presenting a remarkable affinity for type 2 structures, Fucα1-2Galβ1-4GlcNAcβ1, particularly the A2 epitope. The contacts made by DBL1α binding pocket and the ABH histo-blood groups were mapped by theoretical methods and supported by NMR experiments.

Computational protein-ligand docking and virtual drug screening with the AutoDock suite

Computational docking can be used to predict bound conformations and free energies of binding for small-molecule ligands to macromolecular targets. Docking is widely used for the study of biomolecular interactions and mechanisms, and it is applied to structure-based drug design. The methods are fast enough to allow virtual screening of ligand libraries containing tens of thousands of compounds. This protocol covers the docking and virtual screening methods provided by the AutoDock suite of programs, including a basic docking of a drug molecule with an anticancer target, a virtual screen of this target with a small ligand library, docking with selective receptor flexibility, active site prediction and docking with explicit hydration. The entire protocol will require ∼5 h.

A virtual screen discovers novel, fragment-sized inhibitors of Mycobacterium tuberculosis InhA

Isoniazid (INH) is usually administered to treat latent Mycobacterium tuberculosis (Mtb) infections and is used in combination therapy to treat active tuberculosis (TB). Unfortunately, resistance to this drug is hampering its clinical effectiveness. INH is a prodrug that must be activated by Mtb catalase-peroxidase (KatG) before it can inhibit InhA (Mtb enoyl-acyl-carrier-protein reductase). Isoniazid-resistant cases of TB found in clinical settings usually involve mutations in or deletion of katG, which abrogate INH activation. Compounds that inhibit InhA without requiring prior activation by KatG would not be affected by this resistance mechanism and hence would display continued potency against these drug-resistant isolates of Mtb. Virtual screening experiments versus InhA in the GO Fight Against Malaria (GO FAM) project were designed to discover new scaffolds that display base-stacking interactions with the NAD cofactor. GO FAM experiments included targets from other pathogens, including Mtb, when they had structural similarity to a malaria target. Eight of the 16 soluble compounds identified by docking against InhA plus visual inspection were modest inhibitors and did not require prior activation by KatG. The best two inhibitors discovered are both fragment-sized compounds and displayed Ki values of 54 and 59 μM, respectively. Importantly, the novel inhibitors discovered have low structural similarity to known InhA inhibitors and thus help expand the number of chemotypes on which future medicinal chemistry efforts can be focused. These new fragment hits could eventually help advance the fight against INH-resistant Mtb strains, which pose a significant global health threat.

Two isoreticular metal–organic frameworks with CdSO4-like topology: selective gas sorption and drug delivery

In this work, we have also investigated the binding mode of 5-FU to these frameworks, which provides an atomic level view of host–guest interactions.

The enzyme 3-hydroxykynurenine transaminase as potential target for 1,2,4-oxadiazoles with larvicide activity against the dengue vector Aedes aegypti

The mosquito Aedes aegypti is the vector agent responsible for the transmission of yellow fever and dengue fever viruses to over 80 million people in tropical and subtropical regions of the world. Exhaustive efforts have lead to a vaccine candidate with only 30% effectiveness against the dengue virus and failure to protect patients against the serotype 2. Hence, vector control remains the most viable route to dengue fever control programs. We have synthesized a class of 1,2,4-oxadiazole derivatives whose most biologically active compounds exhibit potent activity against Aedes aegypti larvae (ca. of 15 ppm) and low toxicity in mammals. Exposure to these larvicides results in larvae pigmentation in a manner correlated with the LC50 measurements. Structural comparisons of the 1,2,4-oxadiazole nucleus against known inhibitors of insect enzymes allowed the identification of 3-hydroxykynurenine transaminase as a potential target for these synthetic larvicides. Molecular docking calculations indicate that 1,2,4-oxadiazole compounds can bind to 3-hydroxykynurenine transaminase with similar conformation and binding energies as its crystallographic inhibitor 4-(2-aminophenyl)-4-oxobutanoic acid.

Phospholamban modulates the functional coupling between nucleotide domains in Ca-ATPase oligomeric complexes in cardiac sarcoplasmic reticulum

Oligomeric interactions between Ca-ATPase polypeptide chains and their modulation by phospholamban (PLB) were measured in native cardiac sarcoplasmic reticulum (SR) microsomes. Progressive modification of Lys(514) with fluorescein 5-isothiocyanate (FITC), which physically blocks access to the nucleotide binding site by ATP, demonstrates that Ca-ATPase active sites function independently of one another prior to the phosphorylation of PLB. However, upon cAMP-dependent protein kinase (PKA) phosphorylation of PLB, a second-order dependence between residual enzyme activity and the fraction of active sites is observed, consistent with a dimeric functional complex. Complementary distance measurements were made using FITC or 5-iodoacetamidofluorescein (IAF) bound to Cys(674) within the N- or P-domains, respectively, to detect structural coupling within oligomeric complexes. Accompanying the phosphorylation of PLB, neighboring Ca-ATPase polypeptide chains exhibit a 4 +/- 2 A decrease in the proximity between FITC sites within the N-domain and a 9 +/- 3 A increase in the proximity between IAF sites within P-domains. Thus, the phosphorylation of PLB induces spatial rearrangements between the N- and P-domain elements of proximal Ca-ATPase polypeptide chains which restore functional interactions between neighboring polypeptide chains and, in turn, result in increased rates of catalytic turnover. These results are interpreted in terms of a structural model, calculated through optimization of shape complementarity, desolvation, and electrostatic energies, which suggests a dimeric arrangement of Ca-ATPase polypeptide chains through the proximal association of N-domains that accommodates interaction with PLB. We suggest that the phosphorylation of PLB acts to release constraints involving interdomain subunit interactions that enhance catalytically important N-domain motions.

Quantum chemical modeling of enzymatic reactions: the case of 4-oxalocrotonate tautomerase

The reaction mechanism of 4-oxalocrotonate tautomerase (4-OT) is studied using the density functional theory method B3LYP. This enzyme catalyzes the isomerisation of unconjugated alpha-keto acids to their conjugated isomers. Two different quantum chemical models of the active site are devised and the potential energy curves for the reaction are computed. The calculations support the proposed reaction mechanism in which Pro-1 acts as a base to shuttle a proton from the C3 to the C5 position of the substrate. The first step (proton transfer from C3 to proline) is shown to be the rate-limiting step. The energy of the charge-separated intermediate (protonated proline-deprotonated substrate) is calculated to be quite low, in accordance with measured pKa values. The results of the two models are used to evaluate the methodology employed in modeling enzyme active sites using quantum chemical cluster models.

Substrate Orientation in 4-Oxalocrotonate Tautomerase and Its Effect on QM/MM Energy Profiles

The tautomerization of 2-oxo-4E-hexendioate by 4-oxalocrotonate tautomerase has been studied by quantum mechanical/molecular mechanical (QM/MM) methods using three models, A-C, with different substrate orientations. The computed QM/MM energy profiles are rather different. Various energy partitioning analyses indicate the origin of these differences and the role of the active site residues for different substrate orientations. The proposed new model C is preferred over the previously used models A and B because it combines favorable substrate binding geometries with reasonable barriers and is consistent with the experimental evidence from mutation studies concerning the catalytic ability of specific residues in the binding site, especially R11'.

Understanding the Enzymatic Activity of 4-Oxalocrotonate Tautomerase and Its Mutant Analogues: A Computational Study

The effect of replacing arginine residues (Arg) with citrulline residues (Cit) in the binding site of 4-oxalocrotonate tautomerase (4-OT) was investigated with force field molecular dynamics and hybrid quantum mechanics/molecular mechanics studies. It is found that the Arg61Cit mutation has only minor effects on the k(cat) and K(M) values determined experimentally because of the flexibility of this residue. The decrease in k(cat) and increase in K(M) for the Arg11Cit and Arg39Cit mutations are due to the disruption of the binding site, which arises from repulsive interactions with neighboring residues. The results of this investigation shed new light on the effects of mutations in the binding site of 4-OT and consequently on how the enzyme binds the active substrate.

Ab initio QM/MM study shows there is no general acid in the reaction catalyzed by 4-oxalocrotonate tautomerase

The mechanism for the reaction catalyzed by the 4-oxalocrotonate tautomerase (4-OT) enzyme has been studied using a quantum mechanical/molecular mechanical (QM/MM) method developed in our laboratory. Total free energy barriers were obtained for the two steps involved in this reaction. In the first step, Pro-1 acts as a general base to abstract a proton from the third carbon of the substrate, 2-oxo-4-hexenedioate, creating a negative charge on the oxygen at C-2 of this substrate. In the second step, the same hydrogen abstracted by the N-terminal Pro-1 is shuttled back to the fifth carbon of the substrate to form the product, 2-oxo-3-hexenedioate. The calculated total free energy barriers are 14.54 and 16.45 kcal/mol for the first and second steps, respectively. Our calculations clearly show that there is no general acid in the reaction. Arg-39' ', which is hydrogen bonded to the carboxylate group of the substrate, and an ordered water, which moves closer to the site of the charge formed in the transition state and intermediate, play the main role in transition state/intermediate stabilization without acting as general acids in the reaction.

The product of the natural reaction catalyzed by 4-oxalocrotonate tautomerase becomes an affinity label of its mutant

4-Oxalocrotonate tautomerase (4-OT) catalyzes the isomerization of 4-oxalocrotonate, 1, to 2-oxo-3E-hexenedioate, 3, using a general acid/base mechanism that involves a conserved N-terminal proline residue. The P1A and P1G mutants have been shown to catalyze this isomerization but at reduced rates. Analysis of these mutants by mass spectrometry demonstrated that P1A is susceptible to a 1,4-addition of the N-terminal primary amine across the double bond of enone 3 to form a covalent adduct. Although slower than the isomerization reaction, the addition is fast, with 50% of the active sites being alkylated within 12 min. By contrast, the wt4-OT shows no detectable modification over 24 h. These results support the hypothesis that avoidance of nucleophilic reactions, such as the irreversible Michael addition to the product, could be a contributing factor in the evolutionary conservation of N-terminal proline residues in 4OT.

Internal dynamics and ionization states of the macrophage migration inhibitory factor: comparison between wild-type and mutant forms

The macrophage migration inhibitory factor (MIF) is a cytokine that shares a common structural architecture and catalytic strategy with three isomerases: 4-oxalocrotonate tautomerase, 5-carboxymethyl-2-hydroxymuconate isomerase, and D-dopachrome tautomerase. A highly conserved N-terminal proline acts as a base-acid during the proton transfer reaction catalyzed by these enzymes. Such unusual catalytic strategy appears to be possible only due to the N-terminal proline pK(a) shifted to 5.0-6.0 units. Mutations of this residue result in a significant decrease of the catalytic activity of MIF. Two hypotheses have been proposed to explain the catalytic inefficiency of MIF: the lower basicity of primary amines with regard to secondary ones and the increased flexibility resulting from the replacement of a proline by residues like glycine. To investigate that, we have performed molecular dynamics simulations of MIF wild-type and its mutant P1G, as well as calculated the protonation properties of several mutant forms. It was found that the N-terminal glycine does not show larger fluctuations compared to proline, but the former residue is more exposed to the solvent throughout the simulations. The apparent pK(a) of these residues displays very little change (as expected from the structural rigidity of MIF) and is not significantly affected by the surrounding ionizable residues. Instead, the hydrophobic character of the active site seems to be the main factor in determining the pKa of the N-terminal residue and the catalytic efficiency of MIF.

Ionization state and molecular docking studies for the macrophage migration inhibitory factor: the role of lysine 32 in the catalytic mechanism

The macrophage migration inhibitory factor (MIF) is a cytokine that is structurally similar to certain isomerases and for which multiple immune and catalytic roles have been proposed. Different catalytic activities have been reported for MIF, yet the exact mechanism by which MIF acts is not completely known. As a tautomerase, the enzyme uses a general acid-base mechanism of proton transfer in which the amino-terminal proline has been shown to function as the catalytic base. We report the results of molecular docking simulations of macrophage migration inhibitory factor with three substrates, D-dopachrome, L-dopachrome methyl ester and p-(hydroxyphenyl)pyruvate. Electrostatic pK(a) predictions were also performed for the free and complexed forms of the enzyme. The predicted binding mode of p-(hydroxyphenyl)pyruvate is in agreement with the recently published X-ray structure. A model for the binding mode of D-dopachrome and L-dopachrome methyl ester to MIF is proposed which offers insights into the catalytic mechanism of D-dopachrome tautomerase activity of MIF. The proposed catalytic mechanism is further supported by the pK(a) predictions, which suggest that residue Lys32 acts as the general acid for the enzymatic catalysis of D-dopachrome.