Molecular modelling of xylose isomerase catalysis: the role of electrostatics and charge transfer to metals

Metal Dependence of the Xylose Isomerase from Piromyces sp. E2 Explored by Activity Profiling and Protein Crystallography

Xylose isomerase from Piromyces sp. E2 (PirXI) can be used to equip Saccharomyces cerevisiae with the capacity to ferment xylose to ethanol. The biochemical properties and structure of the enzyme have not been described even though its metal content, catalytic parameters, and expression level are critical for rapid xylose utilization. We have isolated the enzyme after high-level expression in Escherichia coli, analyzed the metal dependence of its catalytic properties, and determined 12 crystal structures in the presence of different metals, substrates, and substrate analogues. The activity assays revealed that various bivalent metals can activate PirXI for xylose isomerization. Among these metals, Mn²⁺ is the most favorable for catalytic activity. Furthermore, the enzyme shows the highest affinity for Mn²⁺, which was established by measuring the activation constants (K_act) for different metals. Metal analysis of the purified enzyme showed that in vivo the enzyme binds a mixture of metals that is determined by metal availability as well as affinity, indicating that the native metal composition can influence activity. The crystal structures show the presence of an active site similar to that of other xylose isomerases, with a d-xylose binding site containing two tryptophans and a catalytic histidine, as well as two metal binding sites that are formed by carboxylate groups of conserved aspartates and glutamates. The binding positions and conformations of the metal-coordinating residues varied slightly for different metals, which is hypothesized to contribute to the observed metal dependence of the isomerase activity.

Quantum mechanical modeling: a tool for the understanding of enzyme reactions

Most enzyme reactions involve formation and cleavage of covalent bonds, while electrostatic effects, as well as dynamics of the active site and surrounding protein regions, may also be crucial. Accordingly, special computational methods are needed to provide an adequate description, which combine quantum mechanics for the reactive region with molecular mechanics and molecular dynamics describing the environment and dynamic effects, respectively. In this review we intend to give an overview to non-specialists on various enzyme models as well as established computational methods and describe applications to some specific cases. For the treatment of various enzyme mechanisms, special approaches are often needed to obtain results, which adequately refer to experimental data. As a result of the spectacular progress in the last two decades, most enzyme reactions can be quite precisely treated by various computational methods.

Reintroducing electrostatics into macromolecular crystallographic refinement: application to neutron crystallography and DNA hydration

Most current crystallographic structure refinements augment the diffraction data with a priori information consisting of bond, angle, dihedral, planarity restraints, and atomic repulsion based on the Pauli exclusion principle. Yet, electrostatics and van der Waals attraction are physical forces that provide additional a priori information. Here, we assess the inclusion of electrostatics for the force field used for all-atom (including hydrogen) joint neutron/X-ray refinement. Two DNA and a protein crystal structure were refined against joint neutron/X-ray diffraction data sets using force fields without electrostatics or with electrostatics. Hydrogen-bond orientation/geometry favors the inclusion of electrostatics. Refinement of Z-DNA with electrostatics leads to a hypothesis for the entropic stabilization of Z-DNA that may partly explain the thermodynamics of converting the B form of DNA to its Z form. Thus, inclusion of electrostatics assists joint neutron/X-ray refinements, especially for placing and orienting hydrogen atoms.

Structural and functional analysis of AsbF: origin of the stealth 3,4-dihydroxybenzoic acid subunit for petrobactin biosynthesis

Petrobactin, a virulence-associated siderophore produced by Bacillus anthracis, chelates ferric iron through the rare 3,4-isomer of dihydroxybenzoic acid (3,4-DHBA). Most catechol siderophores, including bacillibactin and enterobactin, use 2,3-DHBA as a biosynthetic subunit. Significantly, siderocalin, a factor involved in human innate immunity, sequesters ferric siderophores bearing the more typical 2,3-DHBA moiety, thereby impeding uptake of iron by the pathogenic bacterial cell. In contrast, the unusual 3,4-DHBA component of petrobactin renders the siderocalin system incapable of obstructing bacterial iron uptake. Although recent genetic and biochemical studies have revealed selected early steps in petrobactin biosynthesis, the origin of 3,4-DHBA as well as the function of the protein encoded by the final gene in the B. anthracis siderophore biosynthetic (asb) operon, asbF (BA1986), has remained unclear. In this study we demonstrate that 3,4-DHBA is produced through conversion of the common bacterial metabolite 3-dehydroshikimate (3-DHS) by AsbF-a 3-DHS dehydratase. Elucidation of the cocrystal structure of AsbF with 3,4-DHBA, in conjunction with a series of biochemical studies, supports a mechanism in which an enolate intermediate is formed through the action of this 3-DHS dehydratase metalloenzyme. Structural and functional parallels are evident between AsbF and other enzymes within the xylose isomerase TIM-barrel family. Overall, these data indicate that microbial species shown to possess homologs of AsbF may, like B. anthracis, also rely on production of the unique 3,4-DHBA metabolite to achieve full viability in the environment or virulence within the host.

Hydride transfer catalyzed by xylose isomerase: mechanism and quantum effects

We have applied molecular dynamics umbrella-sampling simulation and ensemble-averaged variational transition state theory with multidimensional tunneling (EA-VTST/MT) to calculate the reaction rate of xylose-to- xylulose isomerization catalyzed by xylose isomerase in the presence of two Mg2+ ions. The calculations include determination of the free energy of activation profile and ensemble averaging in the transmission coefficient. The potential energy function is approximated by a combined QM/MM/SVB method involving PM3 for the quantum mechanical (QM) subsystem, CHARMM22 and TIP3P for the molecular mechanical (MM) environment, and a simple valence bond (SVB) local function of two bond distances for the hydride transfer reaction. The simulation confirms the essential features of a mechanism postulated on the basis of kinetics and X-ray data by Whitlow et al. (Whitlow, M.; Howard, A. J.; Finzel, B. C.; Poulos, T. L.; Winborne, E.; Gilliland, G. L. Proteins 1991, 9, 153) and Ringe, Petsko, and coworkers (Labie, A.; Allen, K.-N.; Petsko, G. A.; Ringe, D. Biochemistry 1994, 33, 5469). This mechanism involves a rate-determining 1,2-hydride shift with prior and post proton transfers. Inclusion of quantum mechanical vibrational energy is important for computing the free energy of activation, and quantum mechanical tunneling effects are essential for computing kinetic isotope effects (KIEs). It is found that 85% of the reaction proceeds by tunneling and 15% by overbarrier events. The computed KIE for the ratio of hydride to deuteride transfer is in good agreement with the experimental results. The molecular dynamics simulations reveal that proton and hydride transfer reactions are assisted by breathing motions of the mobile Mg2+ ion in the active site, providing evidence for concerted motion of Mg2+ during the hydride transfer step.

Chemical fragmentation in quantum mechanical methods

We give a survey on the application of the chemical fragmentation concept in computer modelling of extended covalent systems. It will be stressed that information on molecular topology, as well as location and composition of the reaction centre allows the construction of a reasonable initial guess for the wave function and thus facilitates the solution of the Schrödinger equation. For systems, where the chemical changes are localised to a few atoms, while others play the role of essentially electrostatic perturbation, a partition into active site and environment is possible providing a background to hybrid quantum mechanical/molecular mechanical (QM/MM) methods. Full molecular orbital treatment of large covalent systems at the minimal basis, semiempirical level becomes possible in the frame of the fragment self-consistent field (FSCF) method which was developed in the past two decades in our laboratory. As an application, we discuss the hydride shift reaction step in xylose isomerase catalysis.

Role of electrostatics at the catalytic metal binding site in xylose isomerase action: Ca(2+)-inhibition and metal competence in the double mutant D254E/D256E

The catalytic metal binding site of xylose isomerase from Arthrobacter B3728 was modified by protein engineering to diminish the inhibitory effect of Ca2+ and to study the competence of metals on catalysis. To exclude Ca2+ from Site 2 a double mutant D254E/D256E was designed with reduced space available for binding. In order to elucidate structural consequences of the mutation the binary complex of the mutant with Mg2+ as well as ternary complexes with bivalent metal ions and the open-chain inhibitor xylitol were crystallized for x-ray studies. We determined the crystal structures of the ternary complexes containing Mg2+, Mn2+, and Ca2+ at 2.2 to 2.5 A resolutions, and refined them to R factors of 16.3, 16.6, and 19.1, respectively. We found that all metals are liganded by both engineered glutamates as well as by atoms O1 and O2 of the inhibitor. The similarity of the coordination of Ca2+ to that of the cofactors as well as results with Be2+ weaken the assumption that geometry differences should account for the catalytic noncompetence of this ion. Kinetic results of the D254E/D256E mutant enzyme showed that the significant decrease in Ca2+ inhibition was accompanied by a similar reduction in the enzymatic activity. Qualitative argumentation, based on the protein electrostatic potential, indicates that the proximity of the negative side chains to the substrate significantly reduces the electrostatic stabilization of the transition state. Furthermore, due to the smaller size of the catalytic metal site, no water molecule, coordinating the metal, could be observed in ternary complexes of the double mutant. Consequently, the proton shuttle step in the overall mechanism should differ from that in the wild type. These effects can account for the observed decrease in catalytic efficiency of the D254E/D256E mutant enzyme.