Catalytic Fields as a Tool to Analyze Enzyme Reaction Mechanism Variants and Reaction Steps.

Catalytic fields representing the topology of the optimal molecular environment charge distribution that reduces the activation barrier have been used to examine alternative reaction variants and to determine the role of conserved catalytic residues for two consecutive reactions catalyzed by the same enzyme. Until now, most experimental and conventional top-down theoretical studies employing QM/MM or ONIOM methods have focused on the role of enzyme electric fields acting on broken bonds of reactants. In contrast, our bottom-up approach dealing with a small reactant and transition-state model allows the analysis of the opposite effects: how the catalytic field resulting from the charge redistribution during the enzyme reaction acts on conserved amino acid residues and contributes to the reduction of the activation barrier.

Extreme Catalytic Power of Ketosteroid Isomerase Related to the Reversal of Proton Dislocations in Hydrogen-Bond Network.

Dynamic electrostatic catalytic field (DECF) vectors derived from transition state and reactant wavefunctions for the two-step reaction occurring within ketosteroid isomerase (KSI) have been calculated using MP2/aug-cc-pVTZ and lower theory levels to determine the magnitude of the catalytic effect and the optimal directions of proton transfers in the KSI hydrogen-bond network. The most surprising and meaningful finding is that the KSI catalytic activity is enhanced by proton dislocations proceeding in opposite directions for each of the two consecutive reaction steps in the same hydrogen network. Such a mechanism allows an ultrafast switching of the catalytic proton wire environment, possibly related to the exceptionally high KSI catalytic power.

Bottom-Up Nonempirical Approach To Reducing Search Space in Enzyme Design Guided by Catalytic Fields.

Currently developed protocols of theozyme design still lead to biocatalysts with much lower catalytic activity than enzymes existing in nature, and, so far, the only avenue of improvement was the in vitro laboratory-directed evolution (LDE) experiments. In this paper, we propose a different strategy based on "reversed" methodology of mutation prediction. Instead of common "top-down" approach, requiring numerous assumptions and vast computational effort, we argue for a "bottom-up" approach that is based on the catalytic fields derived directly from transition state and reactant complex wave functions.

Predicting substituent effects on activation energy changes by static catalytic fields.

Catalytic fields illustrate topology of the optimal charge distribution of a molecular environment reducing the activation energy for any process involving barrier crossing, like chemical reaction, bond rotation etc. Until now, this technique has been successfully applied to predict catalytic effects resulting from intermolecular interactions with individual water molecules constituting the first hydration shell, aminoacid mutations in enzymes or Si→Al substitutions in zeolites. In this contribution, hydrogen to fluorine (H→F) substitution effects for two model reactions have been examined indicating qualitative applicability of the catalytic field concept in the case of systems involving intramolecular interactions.

Application of a simple quantum chemical approach to ligand fragment scoring for Trypanosoma brucei pteridine reductase 1 inhibition.

There is a need for improved and generally applicable scoring functions for fragment-based approaches to ligand design. Here, we evaluate the performance of a computationally efficient model for inhibitory activity estimation, which is composed only of multipole electrostatic energy and dispersion energy terms that approximate long-range ab initio quantum mechanical interaction energies. We find that computed energies correlate well with inhibitory activity for a compound series with varying substituents targeting two subpockets of the binding site of Trypanosoma brucei pteridine reductase 1.

Are Various σ-Hole Bonds Steered by the Same Mechanisms?

Representative Lewis acid-Lewis base complexes linked by tetrel, pnicogen, chalcogen, and halogen bonds have been studied within the quantum theory of atoms in molecules (QTAIM) approach and the hybrid variation-perturbation theory (HVPT) to analyze possible relationships between these σ-hole dimers. Results obtained at the MP2/aug-cc-pVTZ level indicate numerous correlations similar to hydrogen-bonded systems.

Rapid Estimation of Catalytic Efficiency by Cumulative Atomic Multipole Moments: Application to Ketosteroid Isomerase Mutants.

We propose a simple atomic multipole electrostatic model to rapidly evaluate the effects of mutation on enzyme activity and test its performance on wild-type and mutant ketosteroid isomerase. The predictions of our atomic multipole model are similar to those obtained with symmetry-adapted perturbation theory at a fraction of the computational cost. We further show that this approach is relatively insensitive to the precise amino acid side chain conformation in mutants and may thus be useful in computational enzyme (re)design.

Universal short-range ab initio atom-atom potentials for interaction energy contributions with an optimal repulsion functional form.

The repulsion term in conventional force fields constitutes a major source of error. Assuming that this could originate from a too simple analytical functional form, we analyzed various analytical functions using ab initio exchange component values as a reference and obtained (α + β R (-1))exp(-γ R) as the optimal form to represent the repulsion term. Universal exchange, delocalization, and electrostatic penetration potentials approximating the corresponding interaction energy components defined within hybrid variation-perturbation theory (HVPT) were derived using as a reference a training set of 660 biomolecular complexes.

Physical Nature of Fatty Acid Amide Hydrolase Interactions with Its Inhibitors: Testing a Simple Nonempirical Scoring Model.

Fatty acid amide hydrolase (FAAH) is an enzyme responsible for the deactivating hydrolysis of fatty acid ethanolamide neuromodulators. FAAH inhibitors have gained considerable interest due to their possible application in the treatment of anxiety, inflammation, and pain. In the context of inhibitor design, the availability of reliable computational tools for predicting binding affinity is still a challenging task, and it is now well understood that empirical scoring functions have several limitations that in principle could be overcome by quantum mechanics.

Alkaline hydrolysis of organophosphorus pesticides: the dependence of the reaction mechanism on the incoming group conformation.

The fundamental mechanism of organophosphate hydrolysis is the subject of a growing interest resulting from the need for safe disposal of phosphoroorganic pesticides. Herein, we present a detailed ab initio study of the gas-phase mechanisms of alkaline hydrolysis of P-O and P-S bonds in a number of organophosphorus pesticides, including paraoxon, methyl parathion, fenitrothion, demeton-S, acephate, phosalone, azinophos-ethyl, and malathion. Our main finding is that the incoming group conformation influences the mechanism of decomposition of organophosphate and organothiophosphate compounds.

Physical nature of intermolecular interactions in [BMIM][PF6] ionic liquid.

The intermolecular interaction energy in a popular ionic liquid, [BMIM][PF6] is analyzed using the Hybrid Variation-Perturbation Theory approach. The analysis is performed on a sample of configurations from molecular dynamics simulation, instead of minimized structures. The interaction energy components are quantified, showing that the electrostatics is the dominating but not the only important term.

Low cost prediction of relative stabilities of hydrogen bonded complexes from atomic multipole moments for overly short intermolecular distances.

The relative stability of biologically relevant, hydrogen bonded complexes with shortened distances can be assessed at low cost by the electrostatic multipole term alone more successfully than by ab initio methods. These results imply that atomic multipole moments may help improve ligand-receptor ranking predictions, particularly in cases where accurate structural data are not available.

Nonempirical energetic analysis of reactivity and covalent inhibition of fatty acid amide hydrolase.

Fatty acid amide hydrolase (FAAH) is a member of the amidase signature family and is responsible for the hydrolytic deactivation of fatty acid amide neuromodulators, such as anandamide. FAAH carries an unusual catalytic triad consisting of Lys-Ser-Ser, which uniquely enables the enzyme to cleave amides and esters at similar rates. The acylation of 9Z-octadecenamide (oleamide, a FAAH reference substrate) has been widely investigated by computational methods, and those have shown that conformational fluctuations of the active site affect the reaction barrier.

The Ethidium-UA/AU Intercalation Site: Effect of Model Fragmentation and Backbone Charge State.

We report a systematic analysis of the intermolecular interactions of cationic ethidium intercalated into a UA/AU step of RNA for a single conformation based on crystallographic coordinates. Interaction energies at the MP2/6-31G** level were partitioned into electrostatic, exchange, delocalization, and correlation components. Various pairwise interaction models built from chemically intuitive fragments reproduce within a few percent values obtained when treating the intercalation site as a whole.

Physical nature of intermolecular interactions within cAMP-dependent protein kinase active site: differential transition state stabilization in phosphoryl transfer reaction.

The origin of enzyme catalytic activity may be effectively explored within the nonempirical theory of intermolecular interactions. The knowledge of electrostatic, exchange, delocalization, and correlation components of the transition state and substrates stabilization energy arising from each enzyme active site residue allows to examine the most essential physical effects involved in enzymatic catalysis. Consequently, one can build approximate models of the catalytic activity in a systematic and legitimate manner.

Gas-phase mechanisms of degradation of hazardous organophosphorus compounds: do they follow a common pattern of alkaline hydrolysis reaction as in phosphotriesterase?

A comprehensive ab initio analysis of the gas-phase mechanisms of alkaline hydrolysis for a number of phosphotriesterase substrates--O,O-diisopropyl phosphorofluoridate (DFP), O-isopropyl methyl phosphonofluoridate, O,O-diethyl p-nitrophenyl phosphate (paraoxon), O,O-diethyl p-nitrophenyl thiophosphate (parathion), N-acetyl phosphoramidothioate (acephate), O,O-diethyl S-2-ethylthioethyl phosphorothioate (demeton-S) and O-ethyl N,N-dimethyl phosphoramidocyanidate--has been presented herein. The results indicate that, although an associative mechanism of alkaline hydrolysis is followed by all these compounds, P-F and P-CN bonds are cleaved according to the multistep addition-elimination scheme, whereas the breakage of P-O and P-S bonds appears to be consistent with the one-step direct-displacement mechanism. Of the two alternative reaction pathways present in all those cases (except of acephate), the most probable one involves the proton from a nucleophilic hydroxide experiencing an additional stabilization by the phosphoryl oxygen atom.

Ab initio multireference study of Hetero-Diels-Alder reaction of buta-1,3-diene with alkyl glyoxylates.

Hetero-Diels-Alder (HDA) reaction of methyl glyoxylate with buta-1,3-diene has been investigated using multireference methods (complete active space SCF and multi-reference perturbation theory) and compared with several single-reference methods (including DFT) often used in calculations of catalysed [4+2] cycloadditions. Concerted and stepwise mechanisms, found in the literature, are compared. It is shown, that the stepwise mechanism may be a result of choosing unbalanced active space.

Intriguing relations of interaction energy components in stacked nucleic acids.

Major components of the interaction energy that define several approximate levels starting from second order Möller-Plesset theory were studied for 58 stacked nucleic acid dimers. They included typical B-DNA and A-DNA structures, and selected published geometries. A survey of the various terms yields an unexpected correlation between the Pauli exchange and dispersion or correlation terms, which holds for each class of similar planar geometries and for various basis sets.

The molecular basis of urokinase inhibition: from the nonempirical analysis of intermolecular interactions to the prediction of binding affinity.

Urokinase-type plasminogen activator (uPA) is a trypsin-like serine protease that plays a crucial role in angiogenesis process. In addition to its physiological role in healthy organisms, angiogenesis is extremely important in cancer growth and metastasis, resulting in numerous attempts to understand its control and to develop new approaches to anticancer therapy. The alpha-aminoalkylphosphonate diphenyl esters are well known as highly efficient serine protease inhibitors.

Non-empirical study of the phosphorylation reaction catalyzed by 4-methyl-5-beta-hydroxyethylthiazole kinase: relevance of the theory of intermolecular interactions.

The subject of this study was an analysis of the role of active site residues in the phosphoryl transfer reaction catalyzed by 4-methyl-5-beta-hydroxyethylthiazole kinase (ThiK). The ThiK-catalyzed reaction is of special interest due to the lack of a highly conserved aspartate residue serving as a catalytic base. ONIOM(B3LYP:PM3) models of stationary points along the reaction pathway consisted of reactants, two magnesium ions and several highly conserved ThiK active site residues.

Oxime-induced reactivation of sarin-inhibited AChE: a theoretical mechanisms study.

Oximes (especially oximate anions) are used as potential reactivators of OP-inhibited AChE due to their unique alpha-effect nucleophilic reactivity. In the present study, by applying the DFT approach at the B3LYP/6-311G(d,p) level and the Møller-Plesset perturbation theory at the MP2/6-311G(d,p) level, the formoximate-induced reactivation patterns of the sarin-AChE adduct and the corresponding reaction mechanism have been investigated. The potential energy surface along the pathway of the reactivation reaction of sarin-inhibited AChE by oxime reveals that the reaction can occur quickly due to the relatively low energy barriers.

The nature of interactions in the ionic crystal of 3-pentenenitrile, 2-nitro-5-oxo, ion(-1), sodium.

The hybrid variation -- perturbation many-body interaction energy decomposition scheme has been applied to analyze the physical nature of interactions in the ionic 3-pentenenitrile, 2-nitro-5-oxo, ion(-1), sodium crystal, which can be regarded as a model for a large group of aromatic quaternary nitrogen salts. In the crystal structure the sodium ions and water molecules of adjacent unit cells form a positively charged "inorganic layer" with the sodium ions clustered together along the ab faces with the organic (negative) part in between. This puzzling crystal packing is due to a strong favorable interaction between the water molecule and the sodium ions and a substantial charge transfer from the carbanions that balances out the destabilizing sodium-sodium ion repulsion.

How short can the H...H intermolecular contact be? New findings that reveal the covalent nature of extremely strong interactions.

Ab initio calculations at the MP2/6-311++G(d,p) and MP2/aug-cc-pVDZ//MP2/aug-cc-pVTZ levels have been performed for the following complexes: H2OH+...