Leading coordinate analysis of reaction pathways in proton chain transfer: Application to a two-proton transfer model for the green fluorescent protein

The ‘leading coordinate’ approach to computing an approximate reaction pathway, with subsequent determination of the true minimum energy profile, is applied to a two-proton chain transfer model based on the chromophore and its surrounding moieties within the green fluorescent protein (GFP). Using an ab initio quantum chemical method, a number of different relaxed energy profiles are found for several plausible guesses at leading coordinates. The results obtained for different trial leading coordinates are rationalized through the calculation of a two-dimensional relaxed potential energy surface (PES) for the system. Analysis of the 2-D relaxed PES reveals that two of the trial pathways are entirely spurious, while two others contain useful information and can be used to furnish starting points for successful saddle-point searches. Implications for selection of trial leading coordinates in this class of proton chain transfer reactions are discussed, and a simple diagnostic function is proposed for revealing whether or not a relaxed pathway based on a trial leading coordinate is likely to furnish useful information.

The crystal structures of Klebsiella pneumoniae acetolactate synthase with enzyme-bound cofactor and with an unusual intermediate

Acetohydroxyacid synthase (AHAS) and acetolactate synthase (ALS) are thiamine diphosphate (ThDP)-dependent enzymes that catalyze the decarboxylation of pyruvate to give a cofactor-bound hydroxyethyl group, which is transferred to a second molecule of pyruvate to give 2-acetolactate. AHAS is found in plants, fungi, and bacteria, is involved in the biosynthesis of the branched-chain amino acids, and contains non-catalytic FAD. ALS is found only in some bacteria, is a catabolic enzyme required for the butanediol fermentation, and does not contain FAD. Here we report the 2.3-Angstrom crystal structure of Klebsiella pneumoniae ALS. The overall structure is similar to AHAS except for a groove that accommodates FAD in AHAS, which is filled with amino acid side chains in ALS. The ThDP cofactor has an unusual conformation that is unprecedented among the 26 known three-dimensional structures of nine ThDP-dependent enzymes, including AHAS. This conformation suggests a novel mechanism for ALS. A second structure, at 2.0 Angstrom, is described in which the enzyme is trapped halfway through the catalytic cycle so that it contains the hydroxyethyl intermediate bound to ThDP. The cofactor has a tricyclic structure that has not been observed previously in any ThDP-dependent enzyme, although similar structures are well known for free thiamine. This structure is consistent with our proposed mechanism and probably results from an intramolecular proton transfer within a tricyclic carbanion that is the true reaction intermediate. Modeling of the second molecule of pyruvate into the active site of the enzyme with the bound intermediate is consistent with the stereochemistry and specificity of ALS.

Mechanistic aspects of proton chain transfer: A computational study for the green fluorescent protein chromophore

We explore several models for the ground-state proton chain transfer pathway between the green fluorescent protein chromophore and its surrounding protein matrix, with a view to elucidating mechanistic aspects of this process. We have computed quantum chemically the minimum energy pathways (MEPs) in the ground electronic state for one-, two-, and three-proton models of the chain transfer. There are no stable intermediates for our models, indicating that the proton chain transfer is likely to be a single, concerted kinetic step. However, despite the concerted nature of the overall energy profile, a more detailed analysis of the MEPs reveals clear evidence of sequential movement of protons in the chain. The ground-state proton chain transfer does not appear to be driven by the movement of the phenolic proton off the chromophore onto the neutral water bridge. Rather, this proton is the last of the three protons in the chain to move. We find that the first proton movement is from the bridging Ser205 moiety to the accepting Glu222 group. This is followed by the second proton moving from the bridging water to the Ser205for our model this is where the barrier occurs. The phenolic proton on the chromophore is hence the last in the chain to move, transferring to a bridging “water” that already has substantial negative charge.