A quantum chemical study of the mechanism of action of Vitamin K carboxylase (VKC) III. Intermediates and transition states.

A reaction path including transition states is generated for the Dowd mechanism [P. Dowd, R. Hershlne, S.W. Ham, S. Naganathan. Vitamin K and energy transduction: a base strength amplification mechanism. Science 269 (2005) 1684-1691] of action for Vitamin K carboxylase (VKC) using quantum chemical methods (B3LYP/6-311G**). VKC, an essential enzyme in mammalian systems, catalyzes the conversion of hydroquinone form of Vitamin K to the epoxide form in the presence of oxygen. An intermediate species of the oxidation of Vitamin K, an alkoxide, acts apparently to abstract the gamma hydrogen from specifically located glutamate residues. We are able to follow the Dowd proposed path to generate this alkoxide species. The geometries of the proposed model intermediates and transition states in the mechanism are energy optimized. We find that the most energetic step in the mechanism is the uni-deprotonation of the hydroquinone - once this occurs, there is only a small barrier of 3.5kcal/mol for the interaction of oxygen with the carbon to be attacked - and then the reaction proceeds downhill in free energy to form the critical alkoxide species. The results are consistent with the idea that the enzyme probably acts to facilitate the formation of the epoxide by reducing the energy required to deprotonate the hydroquinone form.

A quantum chemical study of the mechanism of action of Vitamin K epoxide reductase (VKOR) II. Transition states.

A reaction path including transition states is generated for the Silverman mechanism [R.B. Silverman, Chemical model studies for the mechanism of Vitamin K epoxide reductase, J. Am. Chem. Soc. 103 (1981) 5939-5941] of action for Vitamin K epoxide reductase (VKOR) using quantum mechanical methods (B3LYP/6-311G**). VKOR, an essential enzyme in mammalian systems, acts to convert Vitamin K epoxide, formed by Vitamin K carboxylase, to its (initial) quinone form for cellular reuse. This study elaborates on a prior work that focused on the thermodynamics of VKOR [D.W. Deerfield II, C.H. Davis, T. Wymore, D.W. Stafford, L.G. Pedersen, Int. J. Quant. Chem. 106 (2006) 2944-2952]. The geometries of proposed model intermediates and transition states in the mechanism are energy optimized. We find that once a key disulfide bond is broken, the reaction proceeds largely downhill. An important step in the conversion of the epoxide back to the quinone form involves initial protonation of the epoxide oxygen. We find that the source of this proton is likely a free mercapto group rather than a water molecule. The results are consistent with the current view that the widely used drug Warfarin likely acts by blocking binding of Vitamin K at the VKOR active site and thereby effectively blocking the initiating step. These results will be useful for designing more complete QM/MM studies of the enzymatic pathway once three-dimensional structural data is determined and available for VKOR.

Molecular recognition of aldehydes by aldehyde dehydrogenase and mechanism of nucleophile activation.

Experimental structural data on the state of substrates bound to class 3 Aldehyde Dehydrogenases (ALDH3A1) is currently unknown. We have utilized molecular mechanics (MM) simulations, in conjunction with new force field parameters for aldehydes, to study the atomic details of benzaldehyde binding to ALDH3A1. Our results indicate that while the nucleophilic Cys243 must be in the neutral state to form what are commonly called near-attack conformers (NACs), these structures do not correlate with increased complexation energy calculated with the MM-Generalized Born Molecular Volume (GBMV) method. The negatively charged Cys243 (thiolate form) of ALDH3A1 also binds benzaldehyde in a stable conformation but in this complex the sulfur of Cys243 is oriented away from benzaldehyde yet yields the most favorable MM-GBMV complexation energy. The identity of the general base, Glu209 or Glu333, in ALDHs remains uncertain. The MM simulations reveal structural and possible functional roles for both Glu209 and Glu333. Structures from the MM simulations that would support either glutamate residue as the general base were further examined with Hybrid Quantum Mechanical (QM)/MM simulations. These simulations show that, with the PM3/OPLS potential, Glu209 must go through a step-wise mechanism to activate Cys243 through an intervening water molecule while Glu333 can go through a more favorable concerted mechanism for the same activation process.

Strategies for multiple sequence alignment.

We present an overview of multiple sequence alignments to outline the practical consequences for the choices among different techniques and parameters. We begin with a discussion of the scoring methods for quantifying the quality of a multiple sequence alignment, followed by a discussion of the algorithms implemented within a variety of multiple sequence alignment programs. We also discuss additional alignment details such as gap penalty and distance metrics. The paper concludes with a discussion on how to improve alignment quality and the limitations of the techniques described in this paper

Initial catalytic events in class 3 aldehyde dehydrogenase: MM and QM/MM simulations.

A novel enzyme mechanism has been predicted by computer simulations for formation of the thiohemiacetal intermediate in the rat ALDH3A1 enzyme. We used molecular mechanics simulations to study the atomic details of substrate binding and quantum mechanical/molecular mechanical methods to study the Cys-243 thiolate attack on benzaldehyde (BA) substrate. BA was found to produce more reactive conformers when aligned for formation of the tetrahedral thiohemiacetal in the R-configuration. In addition, the sulfhydryl proton was seen to be important for initial binding of the substrate. Finally, the free energy differences between forming a thiohemiacetal oxyanion intermediate versus forming a neutral thiohemiacetal intermediate where a proton is donated to the intermediate from the surroundings strongly favor the latter. Our results suggest that the proton donor is the amide proton from the Cys-243 backbone supported by interactions with Lys-235.

Mechanistic implications of the cysteine-nicotinamide adduct in aldehyde dehydrogenase based on quantum mechanical/molecular mechanical simulations.

Recent computer simulations of the cysteine nucleophilic attack on propanal in human mitochondrial aldehyde dehydrogenase (ALDH2) yielded an unexpected result: the chemically reasonable formation of a dead-end cysteine-cofactor adduct when NAD+ was in the "hydride transfer" position. More recently, this adduct found independent crystallographic support in work on formyltetrahydrofolate dehydrogenase, work which further found evidence of the same adduct on re-examination of deposited electron densities of ALDH2. Although the experimental data showed that this adduct was reversible, several mechanistic questions arise from the fact that it forms at all. Here, we present results from further quantum mechanical/molecular mechanical (QM/MM) simulations toward understanding the mechanistic implications of adduct formation. These simulations revealed formation of the oxyanion thiohemiacetal intermediate only when the nicotinamide ring of NAD+ is oriented away from the active site, contrary to prior arguments. In contrast, and in seeming paradox, when NAD is oriented to receive the hydride, disassociation of the oxyanion intermediate to form the dead-end adduct is more thermodynamically favored than maintaining the oxyanion intermediate necessary for catalysis to proceed. However, this disassociation to the adduct could be avoided through proton transfer from the enzyme to the intermediate. Our results continue to indicate that the unlikely source of this proton is the Cys302 main chain amide.

Mathematically complete nucleotide and protein sequence searching using Ssearch.

In this unit a protocol is described for predicting the structure of simple transmembrane a-helical bundles. The protocol is based on a global molecular dynamics search (GMDS) of the configuration space of the helical bundle, yielding several candidate structures. The correct structure among these candidates is selected using information from silent amino acid substitutions, employing the premise that only the correct structure must (by definition) accept all of the silent amino acid substitutions. Thus, the correct structure is found by repeating the GMDS for several close homologs and selecting the structure that persists in all of the trials.

PM3-compatible zinc parameters optimized for metalloenzyme active sites.

Recent studies have shown that semiempirical methods (e.g., PM3 and AM1) for zinc-containing compounds are unreliable for modeling structures containing zinc ions with ligand environments similar to those observed in zinc metalloenzymes. To correct these deficiencies a reparameterization of zinc at the PM3 level was undertaken. In this effort we included frequency corrected B3LYP/6-311G* zinc metalloenzyme ligand environments along with previously utilized experimental data. Average errors for the heats of formation have been reduced from 46.9 kcal/mol (PM3) to 14.2 kcal/mol for this new parameter set, termed ZnB for "Zinc, Biological." In addition, the new parameter sets predict geometries for the Bacillus fragilis active site model and other zinc metalloenzyme mimics that are qualitatively in agreement with high-level ab initio results, something existing parameter sets failed to do.