Phosphate···π, also called anion···π, contacts occur between nucleobases and anionic phosphate oxygens (OP2) in r(GNRA) and r(UNNN) U-turn motifs (N = A,G,C,U; R = A,G). These contacts were investigated using state-of-the-art quantum-chemical methods (QM) to characterize their physicochemical properties and to serve as a reference to evaluate AMBER force field (AFF) performance. We found that phosphate···π interaction energies calculated with the AFF for dimethyl phosphate···nucleobase model systems are less stabilizing in comparison with double-hybrid DFT and that minimum contact distances are larger for all nucleobases. These distance stretches are also observed in large-scale AFF vs QM/MM computations and classical molecular dynamics (MD) simulations on several r(gcGNRAgc) tetraloop hairpins when compared to experimental data extracted from X-ray/cryo-EM structures (res. ≤ 2.5 Å) using the WebFR3D bioinformatic tool. MD simulations further revealed shifted OP2/nucleobase positions. We propose that discrepancies between the QM and AFF result from a combination of missing polarization in the AFF combined with too large AFF Lennard-Jones (LJ) radii of nucleobase carbon atoms in addition to an exaggerated short-range repulsion of the r-12 LJ repulsive term. We compared these results with earlier data gathered on lone pair···π contacts in CpG Z-steps occurring in r(UNCG) tetraloops. In both instances, charge transfer calculations do not support any significant n → π* donation effects. We also investigated thiophosphate···π contacts that showed reduced stabilizing interaction energies when compared to phosphate···π contacts. Thus, we challenge suggestions that the experimentally observed enhanced thermodynamic stability of phosphorothioated r(GNRA) tetraloops can be explained by larger London dispersion.
- MeSH
- Phosphates MeSH
- RNA * chemistry MeSH
- Molecular Dynamics Simulation * MeSH
- Thermodynamics MeSH
- Computational Biology MeSH
- Publication type
- Journal Article MeSH
Mangiferin is a glycosylated xanthone widely distributed in nature, which exhibits wide pharmacological activities, highlighting its anti-cancer properties. Mangiferin interferes with inflammation, lipid, and calcium signaling, which selectively inhibits multiple NFkB target genes as interleukin-6, tumor necrosis factor, plasminogen, and matrix metalloproteinase, among others. In this work, the interactions of this polyphenol with MMP-9 and NF-κβ are characterized by using computational chemistry methods. The results show MMP-9 inhibition by mangiferina is characterized for the interact with the catalytic Zn atom through a penta-coordinate structure. It is also demonstrated through a strong charge transfer established between mangiferin and Zn in the QM/MM study. Concerning the mangiferin/NF-κβ system, the 92.3% of interactions between p50 sub-unity and DNA are maintained with a binding energy of - 8.04 kcal/mol. These findings indicate that mangiferin blocks the p50-p65/DNA interaction resulting in the loss of the functions of this hetero-dimeric member and suggesting inhibition of the cancer progression. Experimental results concerning the anti-cancer properties of mangiferin show that this natural compound can inhibit selectively MMP-9 and NF-ƙβ. Although the anti-tumor properties of mangiferin are well defined, its molecular mechanisms of actions are not described. In this work, a computational study is carried out to characterize the interactions of mangiferin with these molecular targets. The results obtained corroborate the anti-proliferative and anti-apoptotic activity of mangiferin and provide a depiction of its mechanisms of action.
Inhibition of the human O-linked β-N-acetylglucosaminidase (hOGA, GH84) enzyme is pharmacologically relevant in several diseases such as neurodegenerative and cardiovascular disorders, type 2 diabetes, and cancer. Human lysosomal hexosaminidases (hHexA and hHexB, GH20) are mechanistically related enzymes; therefore, selective inhibition of these enzymes is crucial in terms of potential applications. In order to extend the structure-activity relationships of OGA inhibitors, a series of 2-acetamido-2-deoxy-d-glucono-1,5-lactone sulfonylhydrazones was prepared from d-glucosamine. The synthetic sequence involved condensation of N-acetyl-3,4,6-tri-O-acetyl-d-glucosamine with arenesulfonylhydrazines, followed by MnO2 oxidation to the corresponding glucono-1,5-lactone sulfonylhydrazones. Removal of the O-acetyl protecting groups by NH3/MeOH furnished the test compounds. Evaluation of these compounds by enzyme kinetic methods against hOGA and hHexB revealed potent nanomolar competitive inhibition of both enzymes, with no significant selectivity towards either. The most efficient inhibitor of hOGA was 2-acetamido-2-deoxy-d-glucono-1,5-lactone 1-naphthalenesulfonylhydrazone (5f, Ki = 27 nM). This compound had a Ki of 6.8 nM towards hHexB. To assess the binding mode of these inhibitors to hOGA, computational studies (Prime protein-ligand refinement and QM/MM optimizations) were performed, which suggested the binding preference of the glucono-1,5-lactone sulfonylhydrazones in an s-cis conformation for all test compounds.
- MeSH
- Antigens, Neoplasm chemistry metabolism MeSH
- beta-Hexosaminidase beta Chain chemistry metabolism MeSH
- Histone Acetyltransferases chemistry metabolism MeSH
- Hyaluronoglucosaminidase chemistry metabolism MeSH
- Hydrazones chemical synthesis chemistry pharmacology MeSH
- Enzyme Inhibitors chemical synthesis chemistry pharmacology MeSH
- Lactones chemistry MeSH
- Humans MeSH
- Molecular Conformation MeSH
- Models, Molecular MeSH
- Oxides chemistry MeSH
- Manganese Compounds chemistry MeSH
- Structure-Activity Relationship MeSH
- Check Tag
- Humans MeSH
- Publication type
- Journal Article MeSH
The 3'-5', 3'-5' cyclic dinucleotides (3'3'CDNs) are bacterial second messengers that can also bind to the stimulator of interferon genes (STING) adaptor protein in vertebrates and activate the host innate immunity. Here, we profiled the substrate specificity of four bacterial dinucleotide synthases from Vibrio cholerae (DncV), Bacillus thuringiensis (btDisA), Escherichia coli (dgcZ), and Thermotoga maritima (tDGC) using a library of 33 nucleoside-5'-triphosphate analogues and then employed these enzymes to synthesize 24 3'3'CDNs. The STING affinity of CDNs was evaluated in cell-based and biochemical assays, and their ability to induce cytokines was determined by employing human peripheral blood mononuclear cells. Interestingly, the prepared heterodimeric 3'3'CDNs bound to the STING much better than their homodimeric counterparts and showed similar or better potency than bacterial 3'3'CDNs. We also rationalized the experimental findings by in-depth STING-CDN structure-activity correlations by dissecting computed interaction free energies into a set of well-defined and intuitive terms. To this aim, we employed state-of-the-art methods of computational chemistry, such as quantum mechanics/molecular mechanics (QM/MM) calculations, and complemented the computed results with the {STING:3'3'c-di-ara-AMP} X-ray crystallographic structure. QM/MM identified three outliers (mostly homodimers) for which we have no clear explanation of their impaired binding with respect to their heterodimeric counterparts, whereas the R2 = 0.7 correlation between the computed ΔG'int_rel and experimental ΔTm's for the remaining ligands has been very encouraging.
- MeSH
- Bacillus thuringiensis enzymology ultrastructure MeSH
- Cytokines chemistry genetics MeSH
- Escherichia coli enzymology ultrastructure MeSH
- Crystallography, X-Ray MeSH
- Quantum Theory MeSH
- Leukocytes, Mononuclear chemistry enzymology MeSH
- Humans MeSH
- Membrane Proteins chemistry genetics ultrastructure MeSH
- Nucleotides biosynthesis chemistry genetics MeSH
- Immunity, Innate genetics MeSH
- Substrate Specificity MeSH
- Thermotoga maritima enzymology ultrastructure MeSH
- Vibrio cholerae enzymology ultrastructure MeSH
- Structure-Activity Relationship * MeSH
- Check Tag
- Humans MeSH
- Publication type
- Journal Article MeSH
- Research Support, Non-U.S. Gov't MeSH
Halogenated phenols, such as 2,4-dichlorophenol (2,4-DCP) and 4-bromophenol (4-BP) are pollutants generated by a various industrial sectors like chemical, dye, paper bleaching, pharmaceuticals or in an agriculture as pesticides. The use of Horseradish peroxidase (HRP) in the halogenated phenols treatment has already been mentioned, but it is not well understood how the different phenolic substrates can bind in the peroxidase active site nor how these specific interactions can influence in the bioremediation potential. In this work, different removal efficiencies were obtained for phenolic compounds investigated using HRP as catalyst (93.87 and 59.19% to 4BP and 2,4 DCP, respectively). Thus, to rationalize this result based on the interactions of phenols with active center of HRP, we combine computational and experimental methodologies. The theoretical approaches utilized include density functional theory (DFT) calculations, docking simulation and quantum mechanics/molecular mechanics (QM/MM) technique. Michaelis Menten constant (Km) obtained through experimental methodologies were 2.3 and 0.95 mM to 2,4-DCP and 4-BP, respectively, while the specificity constant (Kcat/Km) found was 1.44 mM-1 s-1 and 0.62 mM-1 s-1 for 4-BP and 2,4-DCP, respectively. The experimental parameters appointed to the highest affinity of HRP to 4-BP. According to the molecular docking calculations, both ligands have shown stabilizing intermolecular interaction energies within the HRP active site, however, the 4-BP showed more stabilizing interaction energy (-53.00 kcal mol-1) than 2,4-dichlorophenol (-49.23 kcal mol-1). Besides that, oxidative mechanism of 4-BP and 2,4-DCP was investigated by the hybrid QM/MM approach. This study showed that the lowest activation energy values for transition states investigated were obtained for 4-BP. Therefore, by theoretical approach, the compound 4-BP showed the more stabilizing interaction and activation energy values related to the interaction within the enzyme and the oxidative reaction mechanism, respectively, which corroborates with experimental parameters obtained. The combination between experimental and theoretical approaches was essential to understand how the degradation potential of the HRP enzyme depends on the interactions between substrate and the active center cavity of the enzyme.
Photosystem II (PSII) is a multi-subunit pigment-protein complex and is one of several protein assemblies that function cooperatively in photosynthesis in plants and cyanobacteria. As more structural data on PSII become available, new questions arise concerning the nature of the charge separation in PSII reaction center (RC). The crystal structure of PSII RC from cyanobacteria Thermosynechococcus vulcanus was selected for the computational study of conformational changes in photosystem II associated to the charge separation process. The parameterization of cofactors and lipids for classical MD simulation with Amber force field was performed. The parametrized complex of PSII was embedded in the lipid membrane for MD simulation with Amber in Gromacs. The conformational behavior of protein and the cofactors directly involved in the charge separation were studied by MD simulations and QM/MM calculations. This study identified the most likely mechanism of the proton-coupled reduction of plastoquinone QB. After the charge separation and the first electron transfer to QB, the system undergoes conformational change allowing the first proton transfer to QB- mediated via Ser264. After the second electron transfer to QBH, the system again adopts conformation allowing the second proton transfer to QBH-. The reduced QBH2 would then leave the binding pocket.
In the present work, we performed a complementary quantum mechanical (QM) study to describe the mechanism by which deprotonated pralidoxime (2-PAM) could reactivate human (Homo sapiens sapiens) acetylcholinesterase (HssAChE) inhibited by the nerve agent VX. Such a reaction is proposed to occur in subsequent addition-elimination steps, starting with a nucleophile bimolecular substitution (SN2) mechanism through the formation of a trigonal bipyramidal transition state (TS). A near attack conformation (NAC), obtained in a former study using molecular mechanics (MM) calculations, was taken as a starting point for this project, where we described the possible formation of the TS. Together, this combined QM/MM study on AChE reactivation shows the feasibility of the reactivation occurring via attack of the deprotonated form of 2-PAM against the Ser203-VX adduct of HssAChE.
- MeSH
- Acetylcholinesterase chemistry drug effects MeSH
- Catalytic Domain MeSH
- Quantum Theory MeSH
- Humans MeSH
- Molecular Conformation MeSH
- Organothiophosphorus Compounds pharmacology MeSH
- Pralidoxime Compounds chemistry pharmacology MeSH
- Protons MeSH
- Serine chemistry MeSH
- Molecular Dynamics Simulation MeSH
- Check Tag
- Humans MeSH
- Publication type
- Journal Article MeSH
- Research Support, Non-U.S. Gov't MeSH
Cyclic dinucleotides are second messengers in the cyclic GMP-AMP synthase (cGAS)-stimulator of interferon genes (STING) pathway, which plays an important role in recognizing tumor cells and viral or bacterial infections. They bind to the STING adaptor protein and trigger expression of cytokines via TANK binding kinase 1 (TBK1)/interferon regulatory factor 3 (IRF3) and inhibitor of nuclear factor-κB (IκB) kinase (IKK)/nuclear factor-κB (NFκB) signaling cascades. In this work, we describe an enzymatic preparation of 2'-5',3'-5'-cyclic dinucleotides (2'3'CDNs) with use of cyclic GMP-AMP synthases (cGAS) from human, mouse, and chicken. We profile substrate specificity of these enzymes by employing a small library of nucleotide-5'-triphosphate (NTP) analogues and use them to prepare 33 2'3'CDNs. We also determine affinity of these CDNs to five different STING haplotypes in cell-based and biochemical assays and describe properties needed for their optimal activity toward all STING haplotypes. Next, we study their effect on cytokine and chemokine induction by human peripheral blood mononuclear cells (PBMCs) and evaluate their cytotoxic effect on monocytes. Additionally, we report X-ray crystal structures of two new CDNs bound to STING protein and discuss structure-activity relationship by using quantum and molecular mechanical (QM/MM) computational modeling.
- MeSH
- Biological Assay MeSH
- Cytokines metabolism MeSH
- HEK293 Cells MeSH
- Protein Conformation MeSH
- Leukocytes, Mononuclear drug effects MeSH
- Humans MeSH
- Membrane Proteins chemistry metabolism MeSH
- Nucleotides, Cyclic chemical synthesis pharmacology MeSH
- Computer Simulation MeSH
- Gene Expression Regulation drug effects MeSH
- Structure-Activity Relationship MeSH
- Check Tag
- Humans MeSH
- Publication type
- Journal Article MeSH
- Research Support, Non-U.S. Gov't MeSH
A series of carbamate-based inhibitors of glutamate carboxypeptidase II (GCPII) were designed and synthesized using ZJ-43, N-[[[(1S)-1-carboxy-3-methylbutyl]amino]carbonyl]-l-glutamic acid, as a molecular template in order to better understand the impact of replacing one of the two nitrogen atoms in the urea-based GCPII inhibitor with an oxygen atom. Compound 7 containing a C-terminal 2-oxypentanedioic acid was more potent than compound 5 containing a C-terminal glutamic acid (2-aminopentanedioic acid) despite GCPII's preference for peptides containing an N-terminal glutamate as substrates. Subsequent crystallographic analysis revealed that ZJ-43 and its two carbamate analogs 5 and 7 with the same (S,S)-stereochemical configuration adopt a nearly identical binding mode while (R,S)-carbamate analog 8 containing a d-leucine forms a less extensive hydrogen bonding network. QM and QM/MM calculations have identified no specific interactions in the GCPII active site that would distinguish ZJ-43 from compounds 5 and 7 and attributed the higher potency of ZJ-43 and compound 7 to the free energy changes associated with the transfer of the ligand from bulk solvent to the protein active site as a result of the lower ligand strain energy and solvation/desolvation energy. Our findings underscore a broader range of factors that need to be taken into account in predicting ligand-protein binding affinity. These insights should be of particular importance in future efforts to design and develop GCPII inhibitors for optimal inhibitory potency.
- MeSH
- Cell Line MeSH
- Drosophila genetics MeSH
- Enzyme Assays MeSH
- Glutamate Carboxypeptidase II antagonists & inhibitors chemistry metabolism MeSH
- Protease Inhibitors chemical synthesis chemistry metabolism MeSH
- Carbamates chemical synthesis chemistry metabolism MeSH
- Catalytic Domain MeSH
- Quantum Theory MeSH
- Humans MeSH
- Urea analogs & derivatives chemical synthesis chemistry metabolism MeSH
- Models, Molecular MeSH
- Stereoisomerism MeSH
- Protein Binding MeSH
- Hydrogen Bonding MeSH
- Animals MeSH
- Check Tag
- Humans MeSH
- Animals MeSH
- Publication type
- Journal Article MeSH
- Research Support, Non-U.S. Gov't MeSH
- Research Support, N.I.H., Extramural MeSH
Organophosphorus compounds (OP) nerve agents are among the most toxic chemical substances known. Their toxicity is due to their ability to bind to acetylcholinesterase. Currently, some enzymes, such as phosphotriesterase, human serum paraoxonase 1 and diisopropyl fluorophosphatase, capable of degrading OP, have been characterized. Regarding the importance of bioremediation methods for detoxication of OP, this work aims to study the interaction modes between the human human deoxyuridine triphosphate nucleotidohydrolase (dUTPase) and Sarin and VX, considering their Rp and Sp enantiomers, to evaluate the asymmetric catalysis of those compounds. In previous work, this enzyme has shown good potential to degrade phosphotriesters, and based on this characteristic, we have applied the human dUTPase to the OP degradation. Molecular docking, chemometrics and mixed quantum and molecular mechanics calculations have been employed, showing a good interaction between dUTPase and OP. Two possible reaction mechanisms were tested, and according to our theoretical results, the catalytic degradation of OP by dUTPase can take place via both mechanisms, beyond being stereoselective, that is, dUTPase cleaves one enantiomer preferentially in relation to other. Chemometric techniques provided excellent assistance for performing this theoretical investigation. The dUTPase study shows importance by the fact of it being a human enzyme. Communicated by Ramaswamy H. Sarma.
- MeSH
- Principal Component Analysis MeSH
- Biodegradation, Environmental MeSH
- Catalytic Domain MeSH
- Quantum Theory * MeSH
- Humans MeSH
- Nerve Agents chemistry metabolism MeSH
- Organophosphorus Compounds chemistry metabolism MeSH
- Organothiophosphorus Compounds chemistry metabolism MeSH
- Pyrophosphatases metabolism MeSH
- Sarin chemistry metabolism MeSH
- Molecular Docking Simulation * MeSH
- Hydrogen Bonding MeSH
- Check Tag
- Humans MeSH
- Publication type
- Journal Article MeSH
