Troy Wymore

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. Proteins 2004.

The Class D β-lactamase family: residues governing the maintenance and diversity of function

Class D β-lactamases, a major source of bacterial resistance to β-lactam antibiotic therapies, represent a distinct subset of the β-lactamase superfamily. They share a serine hydrolase mechanism with Classes A/C vs. Class B. Further understanding of their sequence-structure-function relationships would benefit efforts to design a new generation of antibiotics as well as to predict evolutionary mechanisms in response to such therapies. Here we describe analyses based on our high-resolution multiple sequence alignment and phylogenetic tree of ∼80 Class D β-lactamases that leverage several 3D structures of these enzymes. We observe several sequence clusters on the phylogenetic tree, some that are species specific while others include several species from α-, β- and γ-proteobacteria. Residues characteristic of a specific cluster were identified and shown to be located just outside the active site, possibly modulating the function of the catalytic residues to facilitate reactions with specific types of β-lactams. Most significant was the discovery of a likely disulfide bond in a large group composed of α-, β- and γ-proteobacteria that would contribute to enzyme stability and hence bacterial viability under antibiotic assault. A network of co-evolving residues was identified which suggested the importance of maintaining a surface for binding a highly conserved Phe69.

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.

An algorithm for identification and ranking of family-specific residues, applied to the ALDH3 family

An algorithm for detecting amino acid residues characteristic of individual protein families from within aligned collections of paralogous sequences, and its application to the ALDH3 family versus the rest of the ALDH extended family is described. Residues illuminated by this analysis include a key intramolecular tether, a lysine that makes an intersubunit contact at the dimer interface, three residues in close association with the substrate-binding funnel, and a pair of residues suggested to participate in proton relay during the catalytic cycle.

GTKDynamo: a PyMOL plug-in for QC/MM hybrid potential simulations

Hybrid quantum chemical/molecular mechanical (QCMM) potentials are very powerful tools for molecular simulation. They are especially useful for studying processes in condensed phase systems, such as chemical reactions that involve a relatively localized change in electronic structure and where the surrounding environment contributes to these changes but can be represented with more computationally efficient functional forms. Despite their utility, however, these potentials are not always straightforward to apply since the extent of significant electronic structure changes occurring in the condensed phase process may not be intuitively obvious. To facilitate their use, we have developed an open‐source graphical plug‐in, GTKDynamo that links the PyMOL visualization program and the pDynamo QC/MM simulation library. This article describes the implementation of GTKDynamo and its capabilities and illustrates its application to QC/MM simulations.

Molecular dynamics simulation of class 3 aldehyde dehydrogenase.

Molecular dynamics (MD) simulation of the rat class 3 aldehyde dehydrogenase (ALDH) with nicotinamide dinucleotide (NAD) cofactors and explicit water molecules are reported. Our results demonstrate that MD simulation using the latest methodologies can maintain the crystal structure of the enzyme, as well as closely reproduce the short timescale dynamics of the enzyme. Furthermore, the examination of the distance between the nucleophilic Cys-243 and the NAD cofactor reveal important fluctuations that could be linked to ALDH catalysis. Finally, our quantum mechanical model of benzaldehyde in the active site of ALDH demonstrates that the enzyme requires only minor conformational changes to be poised for nucleophilic attack on the substrate.

A mechanism for evolving novel plant sesquiterpene synthase function

Plant sesquiterpene synthases, a subset of the terpene synthase superfamily, are a mechanistically diverse family of enzymes capable of synthesizing hundreds of complex compounds with high regio‐ and stereospecificity and are of biological importance due to their role in plant defense mechanisms. In the current report we describe a large‐scale, high‐resolution phylogenetic analysis of ∼200 plant sesquiterpene synthases integrated with structural and experimental data that address these issues. We observe that all sequences that cluster together on the phylogenetic tree into well‐defined groups share at least the first reaction in the catalytic mechanism subsequent to the initial ionization step and many share steps beyond this down to proton transfers between the enzyme and substrate. Most significant is the previously unreported high conservation of an Asp‐Tyr‐Asp triad. Due to its high conservation, patterns in the phylogenetic tree as well as experimental and modeling results, we suggest that this Asp‐Tyr‐Asp triad is an important functional element responsible for many proton transfers to and from the substrate and intermediates along the plant sesquiterpene synthase catalytic cycle and whose position can be tuned by residues outside the active site that can lead to the evolution of novel enzyme function.

A Distal Disulfide Bridge in OXA-1 β-Lactamase Stabilizes the Catalytic Center and Alters the Dynamics of the Specificity Determining Ω Loop

Widespread antibiotic resistance, particularly when mediated by broad-spectrum β-lactamases, has major implications for public health. Substitutions in the active site often allow broad-spectrum enzymes to accommodate diverse types of β-lactams. Substitutions observed outside the active site are thought to compensate for the loss of thermal stability. The OXA-1 clade of class D β-lactamases contains a pair of conserved cysteines located outside the active site that forms a disulfide bond in the periplasm. Here, the effect of the distal disulfide bond on the structure and dynamics of OXA-1 was investigated via 4 μs molecular dynamics simulations. The results reveal that the disulfide promotes the preorganized orientation of the catalytic residues and affects the conformation of the functionally important Ω loop. Furthermore, principal component analysis reveals differences in the global dynamics between the oxidized and reduced forms, especially in the motions involving the Ω loop. A dynamical network analysis indicates that, in the oxidized form, in addition to its role in ligand binding, the KTG family motif is a central hub of the global dynamics. As activity of OXA-1 has been measured only in the reduced form, we suggest that accurate assessment of its functional profile would require oxidative conditions mimicking periplasm.