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.     

N‐[2‐[(3‐Chlorophenyl)amino]‐phenyl]‐3‐(difluoromethyl)‐1‐methyl‐1H‐pyrazole‐4‐carboxamide: Synthesis,Crystal Structure,Molecular Docking and Biological Activities

In continuation of our previous research on the development of novel pyrazole‐4‐carboxamide with potential antifungal activity, compound SCU2028 , namely N‐[2‐[(3‐chlorophenyl)amino]phenyl]‐3‐(difluoromethyl)‐1‐methyl‐1H‐pyrazole‐4‐carboxamide, was synthesized by new method, structurally characterized by IR, HR‐ESI‐MS, ¹H‐ and ¹³C‐NMR spectra and further identified by single‐crystal X‐ray diffraction. In pot tests, compound SCU2028 showed good in vivo antifungal activity against Rhizoctonia solani (R. solani) and IC₅₀ value of it was 7.48 mg L^?1. In field trials, control efficacy of compound SCU2028 at 200 g.a.i. ha^?1 was 42.30 % on the 7^th day after the first spraying and 68.10 % on the 14^th day after the second spraying, only slightly lower than that of thifluzamide (57.20 % and 71.40 %, respectively). Further in vitro inhibitory activity showed inhibitory ability of compound SCU2028 was 45‐fold higher than that of bixafen and molecular docking of compound SCU2028 to SDH predicted its binding orientation in the active site of the target protein SDH. These results suggested that compound SCU2028 was a potential fungicide for control of rice sheath blight.     

Spectroscopic characterization,DFT studies,molecular docking and cytotoxic evaluation of 4‐nitro‐indole‐3‐carboxaldehyde: A potent lung cancer agent

The 4‐nitro‐1H‐indole‐carboxaldehyde (NICA) molecule was characterized experimentally using FT‐IR, FT‐Raman and UV‐Vis spectra, and it was studied theoretically using DFT calculations. The optimized structure of the NICA molecule was determined by DFT calculations using B3LYP functional with cc‐pVTZ basis set. The electron localization function (ELF) and local orbital localizer (LOL) studies were performed to visualize the electron delocalization in the molecule. The experimental and theoretical wavenumbers of the title molecule were assigned using VEDA 4.0 program. The charge delocalization and stability of the title molecule were investigated using natural bond orbital (NBO) analysis. Frontier molecular orbitals (FMOs) and related molecular properties were calculated. UV‐Vis spectrum was calculated theoretically and validated experimentally. The reactive sites of the molecule were studied from the MEP surface and Fukui function analysis. The molecular docking analysis reveals that the NICA ligand shows better inhibitory activity against RAS, which causes lung cancer. The in vitro cytotoxic activity of the molecule against human lung cancer cell lines (A549) was determined by MTT assay. Thus, the NICA molecule can be used as a potential candidate for the development of the drug against lung cancer.     

1,2,3‐Triazole Tagged 3H‐Pyrano[2,3‐d]pyrimidine‐6‐carboxylate Derivatives: Synthesis,in Vitro Cytotoxicity,Molecular Docking and DNA Interaction Studies

A series of novel ethyl 2,7‐dimethyl‐4‐oxo‐3‐[(1‐phenyl‐1H‐1,2,3‐triazol‐4‐yl)methyl]‐4,5‐dihydro‐3H‐pyrano[2,3‐d]pyrimidine‐6‐carboxylate derivatives 7a – 7m were efficiently synthesized employing click chemistry approach and evaluated for in vitro cytotoxic activity against four tumor cell lines: A549 (human lung adenocarcinoma cell line), HepG2 (human hematoma), MCF‐7 (human breast adenocarcinoma), and SKOV3 (human ovarian carcinoma cell line). Among the compounds tested, the compounds 7a , 7b , 7f , 7l , and 7m have shown potential and selective activity against human lung adenocarcinoma cell line (A549) with IC₅₀ ranging from 0.69 to 6.74 μm . Molecular docking studies revealed that the compounds 7a , 7b , 7f , 7l , and 7m are potent inhibitors of human DNA topoisomerase‐II and also showed compliance with stranded parameters of drug likeness. The calculated binding constants, k_b, from UV/VIS absorptional binding studies of 7a and 7l with CT‐DNA were 10.77 × 10⁴, 6.48 × 10⁴, respectively. Viscosity measurements revealed that the binding could be surface binding mainly due to groove binding. DNA cleavage study showed that 7a and 7l have the potential to cleave pBR322 plasmid DNA without any external agents.     

Challenging a paradigm: Theoretical calculations of the protonation state of the Cys25‐His159 catalytic diad in free papain

A central mechanistic paradigm of cysteine proteases is that the His–Cys catalytic diad forms an ion‐pair NH(+)/S(?) already in the catalytically active free enzyme. Most molecular modeling studies of cysteine proteases refer to this paradigm as their starting point. Nevertheless, several recent kinetics and X‐ray crystallography studies of viral and bacterial cysteine proteases depart from the ion‐pair mechanism, suggesting general base catalysis. We challenge the postulate of the ion‐pair formation in free papain. Applying our QM/SCRF(VS) molecular modeling approach, we analyzed all protonation states of the catalytic diad in free papain and its SMe derivative, comparing the predicted and experimental pK_a data. We conclude that the His–Cys catalytic diad in free papain is fully protonated, NH(+)/SH. The experimental pK_a = 8.62 of His159 imidazole in free papain, obtained by NMR‐controlled titration and originally interpreted as the NH(+)/S(?) ? N/S(?) equilibrium, is now assigned to the NH(+)/SH ? N/SH equilibrium. Proteins 2009. © 2009 Wiley‐Liss, Inc.     

Understanding the disrupting mechanism of the Tau aggregation motif “306VQIVYK311” by phenylthiazolyl‐hydrazides inhibitors

Alzheimer's disease is a progressive neurodegenerative disorder characterized by the abnormal processing of the Tau and the amyloid precursor proteins. The unusual aggregation of Tau is based on the formation of intermolecular β‐sheets through two motifs: ²⁷⁵VQIINK²⁸⁰ and ³⁰⁶VQIVYK³¹¹. Phenylthiazolyl‐hydrazides (PTHs) are capable of inhibiting/disassembling Tau aggregates. However, the disaggregation mechanism of Tau oligomers by PTHs is still unknown. In this work, we studied the disruption of the oligomeric form of the Tau motif ³⁰⁶VQIVYK³¹¹ by PTHs through molecular docking, molecular dynamics, and free energy calculations. We predicted hydrophobic interactions as the major driving forces for the stabilization of Tau oligomer, with V306 and I308 being the major contributors. Nonpolar component of the binding free energy is essential to stabilize Tau‐PTH complexes. PTHs disrupted mainly the van der Waals interactions between the monomers, leading to oligomer destabilization. Destabilization of full Tau filament by PTHs and emodin was not observed in the sampled 20 ns; however, in all cases, the nonpolar component of the binding free energy is essential for the formation of Tau filament‐PTH and Tau filament‐emodin. These results provide useful clues for the design of more effective Tau‐aggregation inhibitors.     

A cytotoxicity,optical spectroscopy and computational binding analysis of 4‐[3‐acetyl‐5‐(acetylamino)‐2‐methyl‐2,3‐dihydro‐1,3,4‐thiadiazole‐2‐yl]phenyl benzoate in calf thymus DNA

In this study the interaction mechanism between newly synthesized 4‐(3‐acetyl‐5‐(acetylamino)‐2‐methyl‐2, 3‐dihydro‐1,3,4‐thiadiazole‐2‐yl) phenyl benzoate (thiadiazole derivative) anticancer active drug with calf thymus DNA was investigated by using various optical spectroscopy techniques along with computational technique. The absorption spectrum shows a clear shift in the lower wavelength region, which may be due to strong hypochromic effect in the ctDNA and the drug. The results of steady state fluorescence spectroscopy show that there is static quenching occurring while increasing the thiadiazole drug concentration in the ethidium bromide‐ctDNA system. Also the binding constant (K), thermo dynamical parameters of enthalpy change (ΔH°), entropy change (ΔS°) Gibbs free energy change (ΔG°) were calculated at different temperature (293 K, 298 K) and the results are in good agreement with theoretically calculated MMGBSA binding analysis. Time resolved emission spectroscopy analysis clearly explains the thiadiazole derivative competitive intercalation in the ethidium bromide‐ctDNA system. Further, molecular docking studies was carried out to understand the hydrogen bonding and hydrophobic interaction between ctDNA and thiadiazole derivative molecule. In addition the docking and molecular dynamics charge distribution analysis was done to understand the internal stability of thiadiazole derivative drug binding sites of ctDNA. The global reactivity of thiadiazole derivative such as electronegativity, electrophilicity and chemical hardness has been calculated.     

New insights into the active site structure and catalytic mechanism of tyrosinase and its related proteins

Tyrosinases are widely distributed in nature. They are copper‐containing oxidases belonging to the type 3 copper protein family, together with catechol oxidases and haemocyanins. Tyrosinases are essential enzymes in melanin biosynthesis and therefore responsible for pigmentation of skin and hair in mammals, where two more enzymes, the tyrosinase‐related proteins (Tyrps), participate in the pathway. The structure and catalytic mechanism of mammalian tyrosinases have been extensively studied but they are not completely understood because of the lack of information on the tertiary structure. The availability of crystallographic data of one plant catechol oxidase and one bacterial tyrosinase has improved the model of the three‐dimensional structure of the active site of the enzyme. Furthermore, sequence comparison of tyrosinase and the Tyrps reveals that the three orthologue proteins share many key structural features, because of their common origin from an ancestral gene, although the specific residues responsible for their different catalytic capabilities have not been identified yet. This review summarizes our current knowledge of tyrosinase and Tyrps structure and function and describes the catalytic mechanism of tyrosinase and Dct/Tyrp2, which are better characterized.     

Structural insight into the catalytic mechanism of gluconate 5‐dehydrogenase from Streptococcus suis: Crystal structures of the substrate‐free and quaternary complex enzymes

Gluconate 5‐dehydrogenase (Ga5DH) is an NADP(H)‐dependent enzyme that catalyzes a reversible oxidoreduction reaction between D ‐gluconate and 5‐keto‐D ‐gluconate, thereby regulating the flux of this important carbon and energy source in bacteria. Despite the considerable amount of physiological and biochemical knowledge of Ga5DH, there is little physical or structural information available for this enzyme. To this end, we herein report the crystal structures of Ga5DH from pathogenic Streptococcus suis serotype 2 in both substrate‐free and liganded (NADP⁺/D ‐gluconate/metal ion) quaternary complex forms at 2.0 Å resolution. Structural analysis reveals that Ga5DH adopts a protein fold similar to that found in members of the short chain dehydrogenase/reductase (SDR) family, while the enzyme itself represents a previously uncharacterized member of this family. In solution, Ga5DH exists as a tetramer that comprised four identical ～29 kDa subunits. The catalytic site of Ga5DH shows considerable architectural similarity to that found in other enzymes of the SDR family, but the S. suis protein contains an additional residue (Arg104) that plays an important role in the binding and orientation of substrate. The quaternary complex structure provides the first clear crystallographic evidence for the role of a catalytically important serine residue and also reveals an amino acid tetrad RSYK that differs from the SYK triad found in the majority of SDR enzymes. Detailed analysis of the crystal structures reveals important contributions of Ca²⁺ ions to active site formation and of specific residues at the C‐termini of subunits to tetramer assembly. Because Ga5DH is a potential target for therapy, our findings provide insight not only of catalytic mechanism, but also suggest a target of structure‐based drug design.     

Design,Synthesis and Bioevaluation of Two Series of 3‐[(1‐Benzyl‐1H‐1,2,3‐triazol‐4‐yl)methyl]quinazolin‐4(3H)‐ones and N‐(1‐Benzylpiperidin‐4‐yl)quinazolin‐4‐amines

Two series of 3‐[(1‐benzyl‐1H‐1,2,3‐triazol‐4‐yl)methyl]quinazolin‐4(3H)‐ones and N‐(1‐benzylpiperidin‐4‐yl)quinazolin‐4‐amines were designed initially as potential acetylcholine esterase inhibitors. Biological evaluation demonstrated that N‐(1‐benzylpiperidin‐4‐yl)quinazolin‐4‐amines significantly inhibited AChE activity. Especially, two compounds of them were found to be the most potent with relative AChE inhibition percentages of 87 % in comparison to donepezil. The docking studies with AChE showed similar interactions between donepezil and four derivatives. N‐(1‐Benzylpiperidin‐4‐yl)quinazolin‐4‐amines also exhibited significant DPPH scavenging effects. The two series of compound also exerted moderate to good cytotoxicity against three human cancer cell lines, including SW620 (human colon cancer), PC‐3 (prostate cancer), and NCI?H23 (lung cancer), with 3‐[(1‐benzyl‐1H‐1,2,3‐triazol‐4‐yl)methyl]quinazolin‐4(3H)‐one being the most cytotoxic agent. 3‐[(1‐Benzyl‐1H‐1,2,3‐triazol‐4‐yl)methyl]quinazolin‐4(3H)‐one significantly induced early apoptosis and arrested the SW620 cells at G2/M phase. From this study, two compounds of N‐(1‐benzylpiperidin‐4‐yl)quinazolin‐4‐amines could serve as new leads for further design and AChE inhibitors, while 3‐[(1‐benzyl‐1H‐1,2,3‐triazol‐4‐yl)methyl]quinazolin‐4(3H)‐one could serve as a new lead for the design and development of more potent anticancer agents.     

Synthesis of 5‐(4′‐carboxyphenyl)‐10,15,20‐tris‐(4 pyridyl)‐porphyrin and its peptidyl phosphonate derivatives

The synthesis and characterization of three new 4‐pyridyl porphyrin‐peptidyl‐phosphonate compounds, containing a diphenyl 3‐pyridylmethyl‐phosphonate moiety, is described in this article. Nitrogen atoms in the pyridine rings of the obtained compounds were alkylated using methyl iodide, to give additional three, water soluble derivatives of these peptidyl‐porphyrin conjugates. All the synthesized compounds could serve as potential photosensitizers for the photodynamic therapy (PDT) method of tumor therapy and displayed activity as inhibitors of aminopeptidase N. Copyright © 2009 European Peptide Society and John Wiley & Sons, Ltd.     

Insights into mechanism of pyrido[2,3‐d]pyrimidines as DYRK1A inhibitors based on molecular dynamic simulations

DYRK1A is characterized by the early development and regulation of neuronal proliferation, and its over expression gives rise to neurological abnormalities. As the promising DYRK1A inhibitors, the binding mechanism between DYRK1A and pyrido[2,3‐d]pyrimidines derivatives at molecular level are still veiled. In this article, it was achieved to get the structural insights into pyrido[2,3‐d]pyrimidines derivatives as DYRK1A inhibitors by means of comprehensive computational approaches involving molecular docking, molecular dynamics simulation, free energy calculation, and energy decomposition analysis. The calculated energy values were highly consistent with the experimental activities. Based on the individual energy terms analysis, the van der Waals interaction was the major leading force in the DYRK1A–ligand interaction. Lys188 was the important residue that formed the hydrogen bond, which improved the inhibitory activity. Furthermore, four novel inhibitors with higher predicted activity were designed based on the obtained findings and confirmed by molecular simulations. Our study is expected to provide significant drug design strategy for the development of more promising DYRK1A inhibitors. Proteins 2016; 84:1108–1123. © 2016 Wiley Periodicals, Inc.     

p‐Coumaric acid decarboxylase from Lactobacillus plantarum: Structural insights into the active site and decarboxylation catalytic mechanism

p‐Coumaric acid decarboxylases (PDCs) catalyze the nonoxidative decarboxylation of hydroxycinnamic acids to generate the corresponding vinyl derivatives. Despite the biotechnological relevance of PDCs in food industry, their catalytic mechanism remains largely unknown. Here, we report insights into the structural basis of catalysis for the homodimeric PDC from Lactobacillus plantarum (LpPDC). The global fold of LpPDC is based on a flattened β‐barrel surrounding an internal cavity. Crystallographic and functional analyses of single‐point mutants of residues located within this cavity have permitted identifying a potential substrate‐binding pocket and also to provide structural evidences for rearrangements of surface loops so that they can modulate the accessibility to the active site. Finally, combination of the structural and functional data with in silico results enables us to propose a two‐step catalytic mechanism for decarboxylation of p‐coumaric acid by PDCs where Glu71 is involved in proton transfer, and Tyr18 and Tyr20 are involved in the proper substrate orientation and in the release of the CO₂ product. Proteins 2010. © 2009 Wiley‐Liss, Inc.     

Computational site‐directed mutagenesis studies of the role of the hydrophobic triad on substrate binding in cholesterol oxidase

Cholesterol oxidase (ChOx) is a flavoenzyme that oxidizes and isomerizes cholesterol (CHL) to form cholest‐4‐en‐3‐one. Molecular docking and molecular dynamics simulations were conducted to predict the binding interactions of CHL in the active site. Several key interactions (E361‐CHL, N485‐FAD, and H447‐CHL) were identified and which are likely to determine the correct positioning of CHL relative to flavin‐adenine dinucleotide (FAD). Binding of CHL also induced changes in key residues of the active site leading to the closure of the oxygen channel. A group of residues, Y107, F444, and Y446, known as the hydrophobic triad, are believed to affect the binding of CHL in the active site. Computational site‐directed mutagenesis of these residues revealed that their mutation affects the conformations of key residues in the active site, leading to non‐optimal binding of CHL and to changes in the structure of the oxygen channel, all of which are likely to reduce the catalytic efficiency of ChOx. Proteins 2017; 85:1645–1655. © 2017 Wiley Periodicals, Inc.     

An improved synthesis of (2S, 4S)‐ and (2S, 4R)‐2‐amino‐4‐methyldecanoic acids: assignment of the stereochemistry of culicinins

An improved synthesis of (2S, 4S)‐ and (2S, 4R)‐2‐amino‐4‐methyldecanoic acids was accomplished using a glutamate derivative as starting material and Evans' asymmetric alkylation as the decisive step. The NMR data of the two diastereomers were measured and compared with those of the natural product. As a result, the stereochemistry of this novel amino acid unit in culicinins was assigned as (2S, 4R). Copyright © 2011 European Peptide Society and John Wiley & Sons, Ltd.