Is the Protein Surrounding the Active Site Critical for Hydrogen Peroxide Reduction by Selenoprotein Glutathione Peroxidase? An ONIOM Study

In silico investigation of falcipain-2 inhibition by hybrid benzimidazole-thiosemicarbazone antiplasmodial agents: A molecular docking, molecular dynamics simulation, and kinetics study

The emergence of artemisinin-resistant variants of Plasmodium falciparum necessitates the urgent search for novel antimalarial drugs. In this regard, an in silico study to screen antimalarial drug candidates from a series of benzimidazole-thiosemicarbazone hybrid molecules with interesting antiplasmodial properties and explore their falcipain-2 (FP2) inhibitory potentials has been undertaken herein. FP2 is a key cysteine protease that degrades hemoglobin in Plasmodium falciparum and is an important biomolecular target in the development of antimalarial drugs. Pharmacokinetic properties, ADMET profiles, MM/GBSA-based binding free energies, reaction mechanisms, and associated barrier heights have been investigated. DFT, molecular dynamics simulation, molecular docking, and ONIOM methods were used. From the results obtained, four 4N-substituted derivatives of the hybrid molecule (E)-2-(1-(5-chloro-1H-benzo[d]imidazol-2-yl)ethylidene)hydrazine-1-carbothioamide (1A) denoted 1B, 1C, 1D, and 1E are drug-like and promising inhibitors of FP2, exhibiting remarkably small inhibitory constants (5.94 × 10-14 - 2.59 × 10-04 n M) and favorable binding free energies (-30.32 to -17.17 kcal/mol). Moreover, the ONIOM results have revealed that 1B and possibly 1C and 1D may act as covalent inhibitors of FP2. The rate-determining step of the thermodynamically favorable covalent binding mechanism occurs across a surmountable barrier height of 24.18 kcal/mol in water and 28.42 kcal/mol in diethyl ether. Our findings are useful for further experimental investigations on the antimalarial activities of the hybrid molecules studied.

Mechanistic Insight into SARS-CoV-2 Mpro Inhibition by Organoselenides: The Ebselen Case Study

The main protease (Mpro) of SARS-CoV-2 is a current target for the inhibition of viral replication. Through a combined Docking and Density Functional Theory (DFT) approach, we investigated in-silico the molecular mechanism by which ebselen (IUPAC: 2-phenyl-1,2-benzoselenazol-3-one), the most famous and pharmacologically active organoselenide, inhibits Mpro. For the first time, we report on a mechanistic investigation in an enzyme for the formation of the covalent -S-Se- bond between ebselen and a key enzymatic cysteine. The results highlight the strengths and weaknesses of ebselen and provide hints for a rational drug design of bioorganic selenium-based inhibitors.

Oxidation of organic diselenides and ditellurides by H2O2 for bioinspired catalyst design

The reactivity of diselenides and ditellurides of general formula (RX)2 (X = Se, Te; R = H, CH3, Ph) toward hydrogen peroxide was studied through a computational approach based on accurate Density Functional Theory (DFT) calculations. The aliphatic and aromatic dichalcogenides have been chosen in light of their activity in glutathione peroxidase (GPx)-like catalytic cycles and their promising features as efficient antioxidant compounds. The reaction products, the energetics and the mechanistic details of these oxidations are discussed. Analogous disulfides are included in our analysis for completeness. We find that the barrier for oxidation of dichalcogenides decreases from disulfides to diselenides to ditellurides. On the other hand, variation of the substituents at the chalcogen nucleus has relatively little effect on the reactivity.

Hydrogen bonding effect between active site and protein environment on catalysis performance in H2-producing [NiFe] hydrogenases

The interaction between the active site and the surrounding protein environment plays a fundamental role in the hydrogen evolution reaction (HER) in [NiFe] hydrogenases. Our density functional theory (DFT) findings demonstrate that the reaction Gibbs free energy required for the rate determining step decreases by 7.1 kcal mol^-1 when the surrounding protein environment is taken into account, which is chiefly due to free energy decreases for the two H⁺/e^- addition steps (the so-called Ni-SIa to I1, and Ni-C to Ni-R), being the largest thermodynamic impediments of the whole reaction. The variety of hydrogen bonds (H-bonds) between the amino acids and the active site is hypothesised to be the main reason for such stability: H-bonds not only work as electrostatic attractive forces that influence the charge redistribution, but more importantly, they act as an electron 'pull' taking electrons from the active site towards the amino acids. Moreover, the electron 'pull' effect through H-bonds via the S^- in cysteine residues shows a larger influence on the energy profile than that via the CN^- ligands on Fe.

Comprehensive modeling of the antioxidant mechanism of ebselen

Three possible catalytic cycles for ebselen have been comprehensively modeled by theoretical calculations using density functional theory (DFT) at a mixed basis set level; the 6-31G(d) basis set for hydrocarbon fragments and the 6-31[Formula: see text]G(d,p) basis set for other atoms. The 2[Formula: see text] cycle is the main pathway in the glutathione peroxidase (GPx) cycle (cycle A), and IM3[Formula: see text]TS3 is the rate controlling process. The 1[Formula: see text]1 cycle is the main pathway for the oxidation cycle (cycle B), and the rate controlling step is the [Formula: see text] step. Ebselen reacts with the selenol 3 to form the diselenide 9, and this is the rate controlling step for cycle C. The extremely high energy barrier for the IM9[Formula: see text]TS9 process indicates that cycle C is unlikely to occur in vivo. Although cycle B is favored based on the energy analysis, with a maximum energy barrier of only 26.68[Formula: see text]kcal/mol at the mixed basis set level, it is generally unlikely to have very high concentrations of peroxides present in vivo. The results indicate that in order to improve the antioxidant activity of ebselen, it would be necessary to suitably modify the molecular structure of ebselen to reduce the energy barrier of the IM3[Formula: see text]TS3 process.

Selenocysteine oxidation in glutathione peroxidase catalysis: an MS-supported quantum mechanics study

Glutathione peroxidases (GPxs) are enzymes working with either selenium or sulfur catalysis. They adopted diverse functions ranging from detoxification of H(2)O(2) to redox signaling and differentiation. The relative stability of the selenoenzymes, however, remained enigmatic in view of the postulated involvement of a highly unstable selenenic acid form during catalysis. Nevertheless, density functional theory calculations obtained with a representative active site model verify the mechanistic concept of GPx catalysis and underscore its efficiency. However, they also allow that the selenenic acid, in the absence of the reducing substrate, reacts with a nitrogen in the active site. MS/MS analysis of oxidized rat GPx4 complies with the predicted structure, an 8-membered ring, in which selenium is bound as selenenylamide to the protein backbone. The intermediate can be re-integrated into the canonical GPx cycle by glutathione, whereas, under denaturing conditions, its selenium moiety undergoes β-cleavage with formation of a dehydro-alanine residue. The selenenylamide bypass prevents destruction of the redox center due to over-oxidation of the selenium or its elimination and likely allows fine-tuning of GPx activity or alternate substrate reactions for regulatory purposes.

Organochalcogen peroxidase mimetics as potential drugs: a long story of a promise still unfulfilled

Organochalcogen compounds have attracted the interest of a multitude of studies to design potential therapeutic agents mimicking the peroxidase activity of selenium-based glutathione peroxidases (GPx's). Starting from the pioneering ebselen, various compounds have been synthesized over the years, which may be traced in three major classes of molecules: cyclic selenenyl amides, diaryl diselenides, and aromatic or aliphatic monoselenides. These compounds share common features and determinants needed to exert an efficient GPx-like activity, such as polarizing groups in close proximity to selenium and steric effects. Nonetheless, the reactivity of selenium, and tellurium as well, poses serious problems for the predictability of the biological effects of these compounds in vivo when used as potential drugs. These molecules, indeed, interfere with thiols of redox-regulated proteins and enzymes, leading to unexpected biological effects. The various chemical aspects of the reaction mechanism of peroxidase mimetics are surveyed here, focusing on experimental evidence and quantum mechanics calculations of organochalcogen representatives of the various classes.

Computational modeling of the direct hydride transfer mechanism for the MAO catalyzed oxidation of phenethylamine and benzylamine: ONIOM (QM/QM) calculations

Monoamine oxidases are two isozymic flavoenzymes which are the important targets for drugs used in the treatment of depression, Parkinson and Alzheimer's diseases. The catalytic reaction taking place between the cofactor FAD and amine substrate is still not completely understood. Herein we employed quantum chemical methods on the recently proposed direct hydride transfer mechanism including full active site residues of MAO isoforms in the calculations. Activation free energy barriers of direct hydride transfer mechanism for MAO-A and MAO-B were calculated by ONIOM (our own n-layered integrated molecular orbital + molecular mechanics) method with QM/QM (quantum mechanics:quantum mechanics) approach employing several density functional theory functionals, B3LYP, WB97XD, CAM-B3LYP and M06-2X, for the high layer. The formation of very recently proposed αC-flavin N5 adduct inside the enzyme has been investigated. ONIOM (M06-2X/6-31+G(d,p):PM6) results revealed that such an adduct may form only in MAO-B suggesting slightly different hydride transfer mechanisms for MAO-A and MAO-B.

Elucidating the catalytic mechanism of β-secretase (BACE1): a quantum mechanics/molecular mechanics (QM/MM) approach

In this quantum mechanics/molecular mechanics (QM/MM) study, the mechanisms of the hydrolytic cleavage of the Met2-Asp3 and Leu2-Asp3 peptide bonds of the amyloid precursor protein (WT-substrate) and its Swedish mutant (SW) respectively catalyzed by β-secretase (BACE1) have been investigated by explicitly including the electrostatic and steric effects of the protein environment in the calculations. BACE1 catalyzes the rate-determining step in the generation of Alzheimer amyloid beta peptides and is widely acknowledged as a promising therapeutic target. The general acid-base mechanism followed by the enzyme proceeds through the following two steps: (1) formation of the gem-diol intermediate and (2) cleavage of the peptide bond. The formation of the gem-diol intermediate occurs with the barriers of 19.6 and 16.1 kcal/mol for the WT- and SW-substrate respectively. The QM/MM energetics predict that with the barriers of 21.9 and 17.2 kcal/mol for the WT- and SW-substrate respectively the cleavage of the peptide bond occurs in the rate-determining step. The computed barriers are in excellent agreement with the measured barrier of ∼18.0 kcal/mol for the SW-substrate and in line with the experimental observation that the cleavage of this substrate is sixty times more efficient than the WT-substrate.

Computational investigation on the structure-activity relationship of the biradical mechanism for monoamine oxidase

Although a considerable amount of mechanistic data has accumulated in literature, the detailed mechanism for amine oxidation by monoamine oxidase is still controversial. The single electron transfer mechanism (SET) has been widely discussed, but not completely understood yet. In the present study, the modified SET mechanism, proposed by Silverman et al., was explored by quantum chemical calculations. The ONIOM method was applied with UDFT/B3LYP/6-31 + G(d,p) for the higher layer and with UHF/6-31G(d) for the lower layer. Isoalloxazin heterocyclic ring and benzylamine were employed in the calculations to represent flavin and the substrate, respectively. The substituents CH(3), OH, OCH(3), H, F, Cl, Br, CF(3) and NO(2) were incorporated at the para position of benzylamine to explore structure-activity relationships. The structures of the reactant complex, transition state and product complex were fully optimized. Activation energies and rate constants of all the reactions were calculated. The results obtained from the linear regression analysis showed that electron-donating groups at the para position of benzylamine increase the reaction rate. A linear but inverse correlation between the log of the calculated rate constants (log k) and the electronic parameter of the substituent was observed (R = 0.93). In accordance with this result, a relatively weak inverse correlation between the calculated log k and the experimental log k was obtained (R = 0.78). The results are contrary to the previous kinetic experiments and the computational study on the effect of p-substituents in the flavin reduction of MAO A by p-substituted benzylamine analogs. Therefore, they present negative evidence for the modeled biradical mechanism.

Matrix metalloproteinase 2 inhibition: combined quantum mechanics and molecular mechanics studies of the inhibition mechanism of (4-phenoxyphenylsulfonyl)methylthiirane and its oxirane analogue

The inhibition mechanism of matrix metalloproteinase 2 (MMP2) by the selective inhibitor (4-phenoxyphenylsulfonyl)methylthiirane (SB-3CT) and its oxirane analogue is investigated computationally. The inhibition mechanism involves C-H deprotonation with concomitant opening of the three-membered heterocycle. SB-3CT was docked into the active site of MMP2, followed by molecular dynamics simulation to prepare the complex for combined quantum mechanics and molecular mechanics (QM/MM) calculations. QM/MM calculations with B3LYP/6-311+G(d,p) for the QM part and the AMBER force field for the MM part were used to examine the reaction of these two inhibitors in the active site of MMP2. The calculations show that the reaction barrier for transformation of SB-3CT is 1.6 kcal/mol lower than its oxirane analogue, and the ring-opening reaction energy of SB-3CT is 8.0 kcal/mol more exothermic than that of its oxirane analogue. Calculations also show that protonation of the ring-opened product by water is thermodynamically much more favorable for the alkoxide obtained from the oxirane than for the thiolate obtained from the thiirane. A six-step partial charge fitting procedure is introduced for the QM/MM calculations to update atomic partial charges of the quantum mechanics region and to ensure consistent electrostatic energies for reactants, transition states, and products.

Quantum mechanics/molecular mechanics investigation of the mechanism of phosphate transfer in human uridine-cytidine kinase 2

The mechanisms of enzyme-catalyzed phosphate transfer and hydrolysis reactions are of great interest due to their importance and abundance in biochemistry. The reaction may proceed in a stepwise fashion, with either a pentavalent phosphorane or a metaphosphate anion intermediate, or by a concerted SN2 mechanism. Despite much theoretical work focused on a few key enzymes, a consensus for the mechanism has not been reached, and examples of all three possibilities have been demonstrated. We have investigated the mechanism of human uridine-cytidine kinase 2 (UCK2, EC 2.7.1.48), which catalyzes the transfer of a phosphate group from ATP to the ribose 5'-hydroxyl of cytidine and uridine. UCK2 is normally expressed in human placenta, but is overexpressed in certain cancer cells, where it is responsible for activating a class of antitumor prodrugs. The UCK2 mechanism was investigated by generating a 2D potential energy surface as a function of the P-O bonds forming and breaking, with energies calculated using a quantum mechanics/molecular mechanics potential (B3LYP/6-31G(d):AMBER). The mechanism of phosphate transfer is shown to be concerted, and is accompanied by concerted proton transfer from the 5'-hydroxyl to a conserved active site aspartic acid that serves as a catalytic base. The calculated barrier for this reaction is 15.1 kcal/mol, in relatively good agreement with the experimental barrier of 17.5 kcal/mol. The interactions of the enzyme active site with the reactant, transition state, and product are examined for their implications on the design of anticancer prodrugs or positron emission tomography (PET) reporter probes for this enzyme.

QM/MM methods for biomolecular systems

Combined quantum-mechanics/molecular-mechanics (QM/MM) approaches have become the method of choice for modeling reactions in biomolecular systems. Quantum-mechanical (QM) methods are required for describing chemical reactions and other electronic processes, such as charge transfer or electronic excitation. However, QM methods are restricted to systems of up to a few hundred atoms. However, the size and conformational complexity of biopolymers calls for methods capable of treating up to several 100,000 atoms and allowing for simulations over time scales of tens of nanoseconds. This is achieved by highly efficient, force-field-based molecular mechanics (MM) methods. Thus to model large biomolecules the logical approach is to combine the two techniques and to use a QM method for the chemically active region (e.g., substrates and co-factors in an enzymatic reaction) and an MM treatment for the surroundings (e.g., protein and solvent). The resulting schemes are commonly referred to as combined or hybrid QM/MM methods. They enable the modeling of reactive biomolecular systems at a reasonable computational effort while providing the necessary accuracy.

General base catalysis for cleavage by the active-site cytosine of the hepatitis delta virus ribozyme: QM/MM calculations establish chemical feasibility

The hepatitis delta virus (HDV) ribozyme is an RNA motif embedded in human pathogenic HDV RNA. Previous experimental studies have established that the active-site nucleotide C75 is essential for self-cleavage of the ribozyme, although its exact catalytic role in the process remains debated. Structural data from X-ray crystallography generally indicate that C75 acts as the general base that initiates catalysis by deprotonating the 2'-OH nucleophile at the cleavage site, while a hydrated magnesium ion likely protonates the 5'-oxygen leaving group. In contrast, some mechanistic studies support the role of C75 acting as general acid and thus being protonated before the reaction. We report combined quantum chemical/molecular mechanical calculations for the C75 general base pathway, utilizing the available structural data for the wild type HDV genomic ribozyme as a starting point. Several starting configurations differing in magnesium ion placement were considered and both one-dimensional and two-dimensional potential energy surface scans were used to explore plausible reaction paths. Our calculations show that C75 is readily capable of acting as the general base, in concert with the hydrated magnesium ion as the general acid. We identify a most likely position for the magnesium ion, which also suggests it acts as a Lewis acid. The calculated energy barrier of the proposed mechanism, approximately 20 kcal/mol, would lower the reaction barrier by approximately 15 kcal/mol compared with the uncatalyzed reaction and is in good agreement with experimental data.

Protein Environment Facilitates O2 Binding in Non-Heme Iron Enzyme. An Insight from ONIOM Calculations on Isopenicillin N Synthase (IPNS)

Binding of dioxygen to a non-heme enzyme has been modeled using the ONIOM combined quantum mechanical/molecular mechanical (QM/MM) method. For the present system, isopenicillin N synthase (IPNS), binding of dioxygen is stabilized by 8-10 kcal/mol for a QM:MM (B3LYP:Amber) protein model compared to a quantum mechanical model of the active site only. In the protein system, the free energy change of O2 binding is close to zero. Two major factors consistently stabilize O2 binding. The first effect, evaluated at the QM level, originates from a change in coordination geometry of the iron center. The active-site model artificially favors the deoxy state (O2 not bound) because it allows too-large rearrangements of the five-coordinate iron site. This error is corrected when the protein is included. The corresponding effect on binding energies is 3-6 kcal/mol, depending on the coordination mode of O2 (side-on or end-on). The second major factor that stabilizes O2 binding is van der Waals interactions between dioxygen and the surrounding enzyme. These interactions, 3-4 kcal/mol at the MM level, are neglected in models that include only the active site. Polarization of the active site by surrounding amino acids does not have a significant effect on the binding energy in the present system.