ONIOM (DFT:MM) Study of the Catalytic Mechanism of myo-Inositol Monophosphatase: Essential Role of Water in Enzyme Catalysis in the Two-Metal Mechanism

Computational Insights into Amide Bond Formation Catalyzed by the Condensation Domain of Nonribosomal Peptide Synthetases

Nonribosomal peptide synthetases (NRPSs) are important enzymes that synthesize an array of nongenetically encoded peptides. The latter have diverse physicochemical properties and roles. NRPSs are modular enzymes in which, for example, the condensation (C-) domain catalyzes the formation of amide bonds. The NRPS tyrocidine synthetase from Brevibacillus brevis is responsible for synthesizing the cyclic-peptide antibiotic tyrocidine. The first step is formation of an amide bond between a proline and phenylalanine which is catalyzed by a C-domain. In this study, a multiscale computational approach (molecular dynamics and QM/MM) has been used to investigate substrate binding and catalytic mechanism of the C-domain of tyrocidine synthetase. Overall, the mechanism is found to proceed through three exergonic steps in which an active site Histidine, His222, acts as a base and acid. First, His222 acts as a base to facilitate nucleophilic attack of the prolyl nitrogen at the phenylalanyl's carbonyl carbon. This is also the rate-limiting step with a free energy barrier of 38.8 kJ mol-1. The second step is collapse of the resulting tetrahedral intermediate with cleavage of the S-C bond between the phenylalanyl and its Ppant arm, along with formation of the above amide bond. Meanwhile, the now protonated His222 imidazole has rotated toward the newly formed thiolate of the Ppant arm. In the final step, His222 acts as an acid, protonating the thiolate and regenerating a neutral His222. The overall mechanism is found to be exergonic with the final product complex being 46.3 kJ mol-1 lower in energy than the initial reactant complex.

Role of water in cyclooxygenase catalysis and design of anti-inflammatory agents targeting two sites of the enzyme

While designing the anti-inflammatory agents targeting cyclooxygenase-2 (COX-2), we first identified a water loop around the heme playing critical role in the enzyme catalysis. The results of molecular dynamic studies supported by the strong hydrogen-bonding equilibria of the participating atoms, radical stabilization energies, the pK_a of the H-donor/acceptor sites and the cyclooxygenase activity of pertinent muCOX-2 ravelled the working of the water–peptide channel for coordinating the flow of H·/electron between the heme and Y385. Based on the working of H·/electron transfer channel between the 12.5 Å distant radical generation and the radical disposal sites, a series of molecules was designed and synthesized. Among this category of compounds, an appreciably potent anti-inflammatory agent exhibiting IC₅₀ 0.06 μM against COX-2 and reversing the formalin induced analgesia and carageenan induced inflammation in mice by 90% was identified. Further it was revealed that, justifying its bidentate design, the compound targets water loop (heme bound site) and the arachidonic acid binding pockets of COX-2.

QM/MM Study of the Reaction Catalyzed by Alkyladenine DNA Glycosylase: Examination of the Substrate Specificity of a DNA Repair Enzyme

Human alkyladenine DNA glycosylase (AAG) functions as part of the base excision repair pathway to excise structurally diverse oxidized and alkylated DNA purines. Specifically, AAG uses a water molecule activated by a general base and a nonspecific active site lined with aromatic residues to cleave the N-glycosidic bond. Despite broad substrate specificity, AAG does not target the natural purines (adenine (A) and guanine (G)). Using the ONIOM(QM:MM) methodology, we provide fundamental atomic level details of AAG bound to DNA-containing a neutral substrate (hypoxanthine (Hx)), a nonsubstrate (G), or a cationic substrate (7-methylguanine (7MeG)) and probe changes in the reaction pathway that occur when AAG targets different nucleotides. We reveal that subtle differences in protein-DNA contacts upon binding different substrates within the flexible AAG active site can significantly affect the deglycosylation reaction. Notably, we predict that AAG excises Hx in a concerted mechanism that is facilitated through correct alignment of the (E125) general base due to hydrogen bonding with a neighboring aromatic amino acid (Y127). Hx departure is further stabilized by π-π interactions with aromatic amino acids and hydrogen bonds with active site water. Despite possessing a similar structure to Hx, G is not excised since the additional exocyclic amino group leads to misalignment of the general base due to disruption of the key E125-Y127 hydrogen bond, the catalytically unfavorable placement of water within the active site, and weakened π-contacts between aromatic amino acids and the nucleobase. In contrast, cationic 7MeG does not occupy the same position within the AAG active site as G due to steric clashes with the additional N7 methyl group, which results in the correct alignment of the general base and permits nucleobase excision as observed for neutral Hx. Overall, our structural data rationalizes the observed substrate specificity of AAG and contributes to our fundamental understanding of enzymes with flexible active sites and broad substrate specificities.

Implicit and explicit solvent effects on the selectivity of the cycloaddition reaction of cyclopentadiene and methyl acrylate; a theoretical study

The product selectivity of the Diels–Alder reaction between cyclopentadiene and methyl acrylate has been studied through a solvent effect analysis by implicit and explicit solvent models. Two paths for this reaction have been proposed and thermodynamic and kinetic parameters for each path were calculated by the B3LYP functional with 6-311++G(d,p) basis set at 298.15 K. In the explicit solvent model, hydrogen bonds were investigated by the QM/MM method. According to the results, the proportion of the exo product in the gas phase is more than for the endo one, which is in contrast to the solvent phase. HOMO–LUMO analyses in the gas phase and different solvents have been done and solvent effects on the energies of the HOMO–LUMO orbitals were investigated. According to the results, the band gap of the TSs was increased in the presence of the solvents. Finally, natural population analysis has been performed to calculate the electronic charges of some atoms that are involved in the centre of the reaction. The data obtained showed that different solvents, including polar and nonpolar solvents, have different effects on the electronic charges of the atoms. Topological analysis of the structures during the reaction by the quantum theory of atoms in molecules (QTAIM) method in the gas phase and different solvents confirmed a cyclic structure for the transition state.