Influence of environment on proton-transfer mechanisms in model triads from theoretical calculations

Substrate dependent reaction channels of the Wolff-Kishner reduction reaction: A theoretical study

Wolff–Kishner reduction reactions were investigated by DFT calculations for the first time. B3LYP/6-311+G(d,p) SCRF=(PCM, solvent = 1,2-ethanediol) optimizations were carried out. To investigate the role of the base catalyst, the base-free reaction was examined by the use of acetone, hydrazine (H₂N–NH₂) and (H₂O)₈. A ready reaction channel of acetone → acetone hydrazine (Me₂C=N–NH₂) was obtained. The channel involves two likely proton-transfer routes. However, it was found that the base-free reaction was unlikely at the N₂ extrusion step from the isopropyl diimine intermediate (Me₂C(H)–N=N–H). Two base-catalyzed reactions were investigated by models of the ketone, H₂N–NH₂ and OH⁻(H₂O)₇. Here, ketones are acetone and acetophenone. While routes of the ketone → hydrazone → diimine are similar, those from the diimines are different. From the isopropyl diimine, the N₂ extrusion and the C–H bond formation takes place concomitantly. The concomitance leads to the propane product concertedly. From the (1-phenyl)ethyl substituted diimine, a carbanion intermediate is formed. The para carbon of the phenyl ring of the anion is subject to the protonation, which leads to a 3-ethylidene-1,4-cyclohexadiene intermediate. Its [1,5]-hydrogen migration gives the ethylbenzene product. For both ketone substrates, the diimines undergoing E2 reactions were found to be key intermediates.

Benchmark Study for the Cysteine-Histidine Proton Transfer Reaction in a Protein Environment: Gas Phase, COSMO, QM/MM Approaches

Proton transfer reactions are of crucial interest for the investigation of proteins. We have investigated the accuracy of commonly used quantum chemical methods for the description of proton transfer reactions in different environments (gas phase, COSMO, QM/MM) using the proton transfer between the catalytic dyad residues cysteine 145 and histidine 41 of SARS coronavirus main protease as a case study. The test includes thermodynamic, kinetic, and structural properties. The study comprises computationally demanding ab initio approaches (HF, CC2, MP2, SCS-CC2, SCS-MP2, CCSD(T)), popular density functional theories (BLYP, B3LYP, M06-2X), and semiempirical methods (MNDO/d, AM1, RM1, PM3, PM6). The approximated coupled cluster approach LCCSD(T) is taken as a reference method. We find that the robustness of the tested methods with respect to the environment correlates well with the level of theory. As an example HF, CC2, MP2, and their SCS variants show similar errors for gas phase, COSMO, or QM/MM computations. In contrast for semiempirical methods, the errors strongly diversify if one goes from gas phase to COSMO or QM/MM. Particular problems are observed for the recent semiempirical methods PM6 and RM1, which show the best performance for gas phase calculations but possess larger errors in conjunction with COSMO. Finally, a combination of SCS-MP2 and B3LYP or M06-2X allows reliable estimates about remaining errors.

First principles study of photostability within hydrogen-bonded amino acids

The photochemistry and photophysics of a two-glycine minimal model is studied at the CASPT2//CASSCF level of theory. Different photoinduced processes are discussed, on the basis of the calculated minimum energy paths and the characterization of the electronic state crossings. Two main processes could provide UV-photostability to the hydrogen-bonded peptide system: (i) forward-backward photoinduced electron/proton transfer involving the H in the hydrogen bond, (ii) singlet-singlet energy transfer between two amino acids, providing ultrafast population of the low-energy n,π* state.

Computational studies on electron and proton transfer in phenol-imidazole-base triads

The electron and proton transfer in phenol-imidazole-base systems (base = NH(2)(-) or OH(-)) were investigated by density-functional theory calculations. In particular, the role of bridge imidazole on the electron and proton transfer was discussed in comparison with the phenol-base systems (base = imidazole, H(2)O, NH(3), OH(-), and NH(2)(-)). In the gas phase phenol-imidazole-base system, the hydrogen bonding between the phenol and the imidazole is classified as short strong hydrogen bonding, whereas that between the imidazole and the base is a conventional hydrogen bonding. The n value in sp(n) hybridization of the oxygen and carbon atoms of the phenolic CO sigma bond was found to be closely related to the CO bond length. From the potential energy surfaces without and with zero point energy correction, it can be concluded that the separated electron and proton transfer mechanism is suitable for the gas-phase phenol-imidazole-base triads, in which the low-barrier hydrogen bond is found and the delocalized phenolic proton can move freely in the single-well potential. For the gas-phase oxidized systems and all of the triads in water solvent, the homogeneous proton-coupled electron transfer mechanism prevails.

Ab initio model study on acetylcholinesterase catalysis: potential energy surfaces of the proton transfer reactions

Ab initio molecular orbital (MO) and hybrid density functional theory (DFT) calculations have been applied to the initial step of the acylation reaction catalyzed by acetylcholinesterase (AChE), which is the nucleophiric addition of Ser200 in catalytic triads to a neurotransmitter acetylcholine (ACh). We focus our attention mainly on the effects of oxyanion hole and Glu327 on the potential energy surfaces (PESs) for the proton transfer reactions in the catalytic triad Ser200-His440-Glu327. The activation barrier for the addition reaction of Ser200 to ACh was calculated to be 23.4 kcal/mol at the B3LYP/6-31G(d)//HF/3-21G(d) level of theory. The barrier height under the existence of oxyanion hole, namely, Ser200-His440-Glu327-ACh-(oxyanion hole) system, decreased significantly to 14.2 kcal/mol, which is in reasonable agreement with recent experimental value (12.0 kcal/mol). Removal of Glu327 from the catalytic triad caused destabilization of both energy of transition state for the reaction and tetrahedral intermediate (product). PESs calculated for the proton transfer reactions showed that the first proton transfer process is the most important in the stabilization of tetrahedral intermediate complex. The mechanism of addition reaction of ACh was discussed on the basis of theoretical results.

Contribution of molecular modeling and site-directed mutagenesis to the identification of two structural residues, Arg-220 and Asp-227, in aminopeptidase A

Aminopeptidase A is a zinc metalloenzyme involved in the formation of brain angiotensin III, which exerts a tonic stimulatory action on the central control of blood pressure. Thus, central inhibitors of aminopeptidase A constitute putative central antihypertensive agents. Mutagenic studies have been performed to investigate organization of the aminopeptidase A active site, with a view to designing such inhibitors. The structure of one monozinc aminopeptidase (leukotriene A(4) hydrolase) was recently resolved and used to construct a three-dimensional model of the aminopeptidase A ectodomain. This new model, highly consistent with the results of mutagenic studies, showed a critical structural interaction between two conserved residues, Arg-220 and Asp-227. Mutagenic replacement of either of these two residues disrupted maturation and subcellular localization and abolished the enzymatic activity of aminopeptidase A, confirming the critical structural role of these residues. In this study, we generated the first three-dimensional model of a strict aminopeptidase, aminopeptidase A. This model constitutes a new tool to probe further the active site of aminopeptidase A and to design new inhibitors of this enzyme.

Computer simulation of ion channel gating: the M(2) channel of influenza A virus in a lipid bilayer

The transmembrane fragment of the influenza virus M(2) protein forms a homotetrameric channel that transports protons. In this paper, we use molecular dynamics simulations to help elucidate the mechanism of channel gating by four histidines that occlude the channel lumen in the closed state. We test two competing hypotheses. In the "shuttle" mechanism, the delta nitrogen atom on the extracellular side of one histidine is protonated by the incoming proton, and, subsequently, the proton on the epsilon nitrogen atom is released on the opposite side. In the "water-wire" mechanism, the gate opens because of electrostatic repulsion between four simultaneously biprotonated histidines. This allows for proton transport along the water wire that penetrates the gate. For each system, composed of the channel embedded in a hydrated phospholipid bilayer, a 1.3-ns trajectory was obtained. It is found that the states involved in the shuttle mechanism, which contain either single-protonated histidines or a mixture of single-protonated histidines plus one biprotonated residue, are stable during the simulations. Furthermore, the orientations and dynamics of water molecules near the gate are conducive to proton transfer. In contrast, the fully biprotonated state is not stable. Additional simulations show that if only two histidines are biprotonated, the channel deforms but the gate remains closed. These results support the shuttle mechanism but not the gate-opening mechanism of proton gating in M(2).