Dependence of enzyme reaction mechanism on protonation state of titratable residues and QM level description: lactate dehydrogenase

Dynamic and Electrostatic Effects on the Reaction Catalyzed by HIV-1 Protease

HIV-1 Protease (HIV-1 PR) is one of the three enzymes essential for the replication process of HIV-1 virus, which explains why it has been the main target for design of drugs against acquired immunodeficiency syndrome (AIDS). This work is focused on exploring the proteolysis reaction catalyzed by HIV-1 PR, with special attention to the dynamic and electrostatic effects governing its catalytic power. Free energy surfaces for all possible mechanisms have been computed in terms of potentials of mean force (PMFs) within hybrid QM/MM potentials, with the QM subset of atoms described at semiempirical (AM1) and DFT (M06-2X) level. The results suggest that the most favorable reaction mechanism involves formation of a gem-diol intermediate, whose decomposition into the product complex would correspond to the rate-limiting step. The agreement between the activation free energy of this step with experimental data, as well as kinetic isotope effects (KIEs), supports this prediction. The role of the protein dynamic was studied by protein isotope labeling in the framework of the Variational Transition State Theory. The predicted enzyme KIEs, also very close to the values measured experimentally, reveal a measurable but small dynamic effect. Our calculations show how the contribution of dynamic effects to the effective activation free energy appears to be below 1 kcal·mol^-1. On the contrary, the electric field created by the protein in the active site of the enzyme emerges as being critical for the electronic reorganization required during the reaction. These electrostatic properties of the active site could be used as a mold for future drug design.

The protonation site of para-dimethylaminobenzoic acid using atmospheric pressure ionization methods

The protonation site of para-dimethylaminobenzoic acid (p-DMABA) was investigated using atmospheric pressure ionization methods (ESI and APCI) coupled with collision-induced dissociation (CID), nuclear magnetic resonance (NMR), and computational chemistry. Theoretical calculations and NMR experiments indicate that the dimethyl amino group is the preferred site of protonation both in the gas phase and aqueous solution. Protonation of p-DMABA occurs at the nitrogen atom by ESI independent of the solvents and other operation conditions under typical thermodynamic control. However, APCI produces a mixture of the nitrogen- and carbonyl oxygen-protonated p-DMABA when aprotic organic solvents (acetonitrile, acetone, and tetrahydrofuran) are used, exhibiting evident kinetic characteristics of protonation. But using protic organic solvents (methanol, ethanol, and isopropanol) in APCI still leads to the formation of thermodynamically stable N-protonated p-DMABA. These structural assignments were based on the different CID behavior of the N- and O-protonated p-DMABA. The losses of methyl radical and water are the diagnostic fragmentations of the N- and O-protonated p-DMABA, respectively. In addition, the N-protonated p-DMABA is more stable than the O-protonated p-DMABA in CID revealed by energy resolved experiments and theoretical calculations.

Theoretical studies of the hydrolysis of antibiotics catalyzed by a metallo-β-lactamase

In this paper, hybrid QM/MM molecular dynamics (MD) simulations have been performed to explore the mechanisms of hydrolysis of two antibiotics, Imipenen (IMI), an antibiotic belonging to the subgroup of carbapenems, and the Cefotaxime (CEF), a third-generation cephalosporin antibiotic, in the active site of a mono-nuclear β-lactamase, CphA from Aeromonas hydrophila. Significant different transition state structures are obtained for the hydrolysis of both antibiotics: while the TS of the CEF is an ionic species with negative charge on nitrogen, the IMI TS presents a tetrahedral-like character with negative charge on oxygen atom of the carbonyl group of the lactam ring. Thus, dramatic conformational changes can take place in the cavity of CphA to accommodate different substrates, which would be the origin of its substrate promiscuity. Since CphA shows only activity against carbapenem antibiotic, this study sheds some light into the origin of the selectivity of the different MbL and, as a consequence, into the discovery of specific and potent MβL inhibitors against a broad spectrum of bacterial pathogens. We have finally probed that a re-parametrization of semiempirical methods should be done to properly describe the behavior the metal cation in active site, Zn(2+), when used in QM/MM calculations.

Catalytic mechanism of porphobilinogen synthase: the chemical step revisited by QM/MM calculations

Porphobilinogen synthase (PBGS) catalyzes the asymmetric condensation and cyclization of two 5-aminolevulinic acid (5-ALA) substrate molecules to give porphobilinogen (PBG). The chemical step of PBGS is herein revisited using QM/MM (ONIOM) calculations. Two different protonation states and several different mechanisms are considered. Previous mechanisms based on DFT-only calculations are shown unlikely to occur. According to these new calculations, the deprotonation step rather than ring closure is rate-limiting. Both the C-C bond formation first mechanism and the C-N bond formation first mechanism are possible, depending on how the A-site ALA binds to the enzyme. We furthermore propose that future work should focus on the substrate binding step rather than the enzymatic mechanism.

A hybrid elastic band string algorithm for studies of enzymatic reactions

A common challenge in theoretical biophysics is the identification of a minimum energy path (MEP) for the rearrangement of a group of atoms from one stable configuration to another. The structure with maximum energy along the MEP approximates the transition state for the process and the energy profile itself permits estimation of the transition rates. In this work we describe a computationally efficient algorithm for the identification of minimum energy paths in complicated biosystems. The algorithm is a hybrid of the nudged elastic band (NEB) and string methods. It has been implemented in the pDynamo simulation program and tested by examining elementary steps in the reaction mechanisms of three enzymes: citrate synthase, RasGAP, and lactate dehydrogenase. Good agreement is found for the energies and geometries of the species along the reaction profiles calculated using the new algorithm and previous versions of the NEB and string techniques, and also those obtained by the common method of adiabatic exploration of the potential energy surface as a function of predefined reaction coordinates. Precisely refined structures of the saddle points along the paths may be subsequently obtained with the climbing image variant of the NEB algorithm. Directions in which the utility of the methods that we have implemented can be further improved are discussed.

Dynamic mechanism of proton transfer in mannitol 2-dehydrogenase from Pseudomonas fluorescens: mobile GLU292 controls proton relay through a water channel that connects the active site with bulk solvent

The active site of mannitol 2-dehydrogenase from Pseudomonas fluorescens (PfM2DH) is connected with bulk solvent through a narrow protein channel that shows structural resemblance to proton channels utilized by redox-driven proton pumps. A key element of the PfM2DH channel is the "mobile" Glu(292), which was seen crystallographically to adopt distinct positions up and down the channel. It was suggested that the "down → up" conformational change of Glu(292) could play a proton relay function in enzymatic catalysis, through direct proton shuttling by the Glu or because the channel is opened for water molecules forming a chain along which the protons flow. We report evidence from site-directed mutagenesis (Glu(292) → Ala) substantiated by data from molecular dynamics simulations that support a role for Glu(292) as a gate in a water chain (von Grotthuss-type) mechanism of proton translocation. Occupancy of the up and down position of Glu(292) is influenced by the bonding and charge state of the catalytic acid base Lys(295), suggesting that channel opening/closing motions of the Glu are synchronized to the reaction progress. Removal of gatekeeper control in the E292A mutant resulted in a selective, up to 120-fold slowing down of microscopic steps immediately preceding catalytic oxidation of mannitol, consistent with the notion that formation of the productive enzyme-NAD(+)-mannitol complex is promoted by a corresponding position change of Glu(292), which at physiological pH is associated with obligatory deprotonation of Lys(295) to solvent. These results underscore the important role of conformational dynamics in the proton transfer steps of alcohol dehydrogenase catalysis.

Importance of the lactate dehydrogenase quaternary structure in theoretical calculations

Using the example of lactate dehydrogenase, we show that enzyme quaternary structure has an important influence on the structure of the active site and that models that comprise all amino acids in the vicinity of an active site, but are missing this structural information, can lead to incorrect results. We also show that binding isotope effects are very sensitive to the geometric parameters, and thus one should be very cautious when interpreting results obtained with models that are too coarse. In terms of the type of hydrogen bonds, our results indicate that binding isotope effects are pronounced only when a hydrogen bond exhibits some covalent character.

Theoretical site-directed mutagenesis: Asp168Ala mutant of lactate dehydrogenase

Molecular simulations based on the use of hybrid quantum mechanics/molecular mechanics methods are able to provide detailed information about the complex enzymatic reactions and the consequences of specific mutations on the activity of the enzyme. In this work, the reduction of pyruvate to lactate catalysed by wild-type and Asp168Ala mutant lactate dehydrogenase (LDH) has been studied by means of simulations using a very flexible molecular model consisting of the full tetramer of the enzyme, together with the cofactor NADH, the substrate and solvent water molecules. Our results indicate that the Asp168Ala mutation provokes a shift in the pKa value of Glu199 that becomes unprotonated at neutral pH in the mutant enzyme. This change compensates the loss of the negative charge of Asp168, rendering a still active enzyme. Thus, our methodology gives a calculated barrier height for the Asp168Ala mutant 3 kcal mol-1 higher than that for wild-type LDH, which is in very good agreement with the experiment. The computed potential energy surfaces reveal the reaction pathways and transition structures for the wild-type and mutant enzymes. Hydride transfer is less advanced and the proton transfer is more advanced in the Asp168Ala mutant than in the wild type. This approach provides a very powerful tool for the analysis of the roles of key active-site residues.

A theoretical study of the catalytic mechanism of formate dehydrogenase

A theoretical study of the hydride transfer between formate anion and nicotinamide adenine dinucleotide (NAD(+)) catalyzed by the enzyme formate dehydrogenase (FDH) has been carried out by a combination of two hybrid quantum mechanics/molecular mechanics techniques: statistical simulation methods and internal energy minimizations. Free energy profiles, obtained for the reaction in the enzyme active site and in solution, allow obtaining a comparative analysis of the behavior of both condensed media. Moreover, calculations of the reaction in aqueous media can be used to probe the dramatic differences between reactants state in the enzyme active site and in solution. The results suggest that the enzyme compresses the substrate and the cofactor into a conformation close to the transition structure by means of favorable interactions with the amino acid residues of the active site, thus facilitating the relative orientation of donor and acceptor atoms to favor the hydride transfer. Moreover, a permanent field created by the protein reduces the work required to reach the transition state (TS) with a concomitant polarization of the cofactor that would favor the hydride transfer. In contrast, in water the TS is destabilized with respect to the reactant species because the polarity of the solute diminishes as the reaction proceeds, and consequently the reaction field, which is created as a response to the change in the solute polarity, is also decreased. Therefore protein structure is responsible of both effects; substrate preorganization and TS stabilization thus diminishing the activation barrier. Because of the electrostatic features of the catalyzed reaction, both media preferentially stabilize the ground-state, thus explaining the small rate constant enhancement of this enzyme, but FDH does so to a much lower extent than aqueous solution. Finally, a good agreement between experimental and theoretical kinetic isotope effects is found, thus giving some credit to our results.

A theoretical analysis of rate constants and kinetic isotope effects corresponding to different reactant valleys in lactate dehydrogenase

In some enzymatic systems large conformational changes are coupled to the chemical step, in such a way that dispersion of rate constants can be observed in single-molecule experiments, each corresponding to the reaction from a different reactant valley. Under this perspective here we present a computational study of pyruvate to lactate transformation catalyzed by lactate dehydrogenase. The reaction consists of a hydride transfer and a proton transfer that seem to take place concertedly. The degree of asynchronicity and the energy barrier depend on the particular starting reactant valley. In order to estimate rate constants we used a free energy perturbation technique adapted to follow the intrinsic reaction coordinate for several possible reaction paths. Tunneling effects are also obtained with a slightly modified version of the ensemble-averaged variational transition state theory with multidimensional tunneling contributions. According to our results the closure of the active site by means of a flexible loop can lead to the formation of different reactant complexes displaying different features in the disposition of some key residues (such as Arg109), interactions with the substrate and number of water molecules in the active site. The chemical step of the reaction takes place with a different reaction rate from each of these complexes. Finally, primary kinetic isotope effects for replacement of the transferring hydrogen of the cofactor with a deuteride are in good agreement with experimental observations, thus validating our methodology.