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

Structural Insight into the Catalytic Mechanisms of an L-Sorbosone Dehydrogenase

L-Sorbosone dehydrogenase (SNDH) is a key enzyme involved in the biosynthesis of 2-keto-L-gulonic acid , which is a direct precursor for the industrial scale production of vitamin C. Elucidating the structure and the catalytic mechanism is essential for improving SNDH performance. By solving the crystal structures of SNDH from Gluconobacter oxydans WSH-004, a reversible disulfide bond between Cys295 and the catalytic Cys296 residues is discovered. It allowed SNDH to switch between oxidation and reduction states, resulting in opening or closing the substrate pocket. Moreover, the Cys296 is found to affect the NADP+ binding pose with SNDH. Combining the in vitro biochemical and site-directed mutagenesis studies, the redox-based dynamic regulation and the catalytic mechanisms of SNDH are proposed. Moreover, the mutants with enhanced activity are obtained by extending substrate channels. This study not only elucidates the physiological control mechanism of the dehydrogenase, but also provides a theoretical basis for engineering similar enzymes.

Aryl Coenzyme A Ligases, a Subfamily of the Adenylate-Forming Enzyme Superfamily

Aryl coenzyme A (CoA) ligases belong to class I of the adenylate-forming enzyme superfamily (ANL superfamily). They catalyze the formation of thioester bonds between aromatic compounds and CoA and occur in nearly all forms of life. These ligases are involved in various metabolic pathways degrading benzene, toluene, ethylbenzene, and xylene (BTEX) or polycyclic aromatic hydrocarbons (PAHs). They are often necessary to produce the central intermediate benzoyl-CoA that occurs in various anaerobic pathways. The substrate specificity is very diverse between enzymes within the same class, while the dependency on Mg²⁺, ATP, and CoA as well as oxygen insensitivity are characteristics shared by the whole enzyme class. Some organisms employ the same aryl-CoA ligase when growing aerobically and anaerobically, while others induce different enzymes depending on the environmental conditions. Aryl-CoA ligases can be divided into two major groups, benzoate:CoA ligase-like enzymes and phenylacetate:CoA ligase-like enzymes. They are widely distributed between the phylogenetic clades of the ANL superfamily and show closer relationships within the subfamilies than to other aryl-CoA ligases. This, together with residual CoA ligase activity in various other enzymes of the ANL superfamily, leads to the conclusion that CoA ligases might be the ancestral proteins from which all other ANL superfamily enzymes developed.

In silico analysis of class I adenylate-forming enzymes reveals family and group-specific conservations

Luciferases, aryl- and fatty-acyl CoA synthetases, and non-ribosomal peptide synthetase proteins belong to the class I adenylate-forming enzyme superfamily. The reaction catalyzed by the adenylate-forming enzymes is categorized by a two-step process of adenylation and thioesterification. Although all of these proteins perform a similar two-step process, each family may perform the process to yield completely different results. For example, luciferase proteins perform adenylation and oxidation to produce the green fluorescent light found in fireflies, while fatty-acyl CoA synthetases perform adenylation and thioesterification with coenzyme A to assist in metabolic processes involving fatty acids. This study aligned a total of 374 sequences belonging to the adenylate-forming superfamily. Analysis of the sequences revealed five fully conserved residues throughout all sequences, as well as 78 more residues conserved in at least 60% of sequences aligned. Conserved positions are involved in magnesium and AMP binding and maintaining enzyme structure. Also, ten conserved sequence motifs that included most of the conserved residues were identified. A phylogenetic tree was used to assign sequences into nine different groups. Finally, group entropy analysis identified novel conservations unique to each enzyme group. Common group-specific positions identified in multiple groups include positions critical to coordinating AMP and the CoA-bound product, a position that governs active site shape, and positions that help to maintain enzyme structure through hydrogen bonds and hydrophobic interactions. These positions could serve as excellent targets for future research.

In silico analysis of heme oxygenase structural homologues identifies group-specific conservations

Heme oxygenases (HO) catalyze the breakdown of heme, aiding the recycling of its components. Several other enzymes have homologous tertiary structures to HOs, while sharing little sequence homology. These homologues include thiaminases, the hydroxylase component of methane monooxygenases, and the R2 component of Class I ribonucleotide reductases (RNR). This study compared these structural homologues of HO, using a large number of protein sequences for each homologue. Alignment of a total of 472 sequences showed little sequence conservation, with no residues having conservation in more than 80% of aligned sequences and only five residues conserved in at least 60% of the sequences. Fourteen additional positions, most of which were critical for hydrophobic packing, displayed amino acid similarity of 60% or higher. Ten conserved sequence motifs were identified in HOs and RNRs. Phylogenetic analysis verified the existence of the four distinct groups of HO homologues, which were then analyzed by group entropy analysis to identify residues critical to the unique function of each enzyme. Other methods for determining functional residues were also performed. Several common index positions identified represent critical evolutionary changes that resulted in the unique function of each enzyme, suggesting potential targets for site‐directed mutagenesis. These positions included residues that coordinate ligands, form the active sites, and maintain enzyme structure.

Enzymes

Heme oxygenase (EC 1.14.14.18), methane monooxygenase (EC 1.14.13.25), ribonucleotide reductase (EC 1.17.4.1), thiaminase II (EC 3.5.99.2).

Analysis of nucleotide diphosphate sugar dehydrogenases reveals family and group-specific relationships

UDP‐glucose dehydrogenase (UDPGDH), UDP‐N‐acetyl‐mannosamine dehydrogenase (UDPNAMDH) and GDP‐mannose dehydrogenase (GDPMDH) belong to a family of NAD⁺‐linked 4‐electron‐transfering oxidoreductases called nucleotide diphosphate sugar dehydrogenases (NDP‐SDHs). UDPGDH is an enzyme responsible for converting UDP‐d‐glucose to UDP‐d‐glucuronic acid, a product that has different roles depending on the organism in which it is found. UDPNAMDH and GDPMDH convert UDP‐N‐acetyl‐mannosamine to UDP‐N‐acetyl‐mannosaminuronic acid and GDP‐mannose to GDP‐mannuronic acid, respectively, by a similar mechanism to UDPGDH. Their products are used as essential building blocks for the exopolysaccharides found in organisms like Pseudomonas aeruginosa and Staphylococcus aureus. Few studies have investigated the relationships between these enzymes. This study reveals the relationships between the three enzymes by analysing 229 amino acid sequences. Eighteen invariant and several other highly conserved residues were identified, each serving critical roles in maintaining enzyme structure, coenzyme binding or catalytic function. Also, 10 conserved motifs that included most of the conserved residues were identified and their roles proposed. A phylogenetic tree demonstrated relationships between each group and verified group assignment. Finally, group entropy analysis identified novel conservations unique to each NDP‐SDH group, including residue positions critical to NDP‐sugar substrate interaction, enzyme structure and intersubunit contact. These positions may serve as targets for future research.

Enzymes

UDP‐glucose dehydrogenase (UDPGDH, EC 1.1.1.22).

The Class D beta-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.

MPI-PHYLIP: parallelizing computationally intensive phylogenetic analysis routines for the analysis of large protein families

BACKGROUND

Phylogenetic study of protein sequences provides unique and valuable insights into the molecular and genetic basis of important medical and epidemiological problems as well as insights about the origins and development of physiological features in present day organisms. Consensus phylogenies based on the bootstrap and other resampling methods play a crucial part in analyzing the robustness of the trees produced for these analyses.

METHODOLOGY

Our focus was to increase the number of bootstrap replications that can be performed on large protein datasets using the maximum parsimony, distance matrix, and maximum likelihood methods. We have modified the PHYLIP package using MPI to enable large-scale phylogenetic study of protein sequences, using a statistically robust number of bootstrapped datasets, to be performed in a moderate amount of time. This paper discusses the methodology used to parallelize the PHYLIP programs and reports the performance of the parallel PHYLIP programs that are relevant to the study of protein evolution on several protein datasets.

CONCLUSIONS

Calculations that currently take a few days on a state of the art desktop workstation are reduced to calculations that can be performed over lunchtime on a modern parallel computer. Of the three protein methods tested, the maximum likelihood method scales the best, followed by the distance method, and then the maximum parsimony method. However, the maximum likelihood method requires significant memory resources, which limits its application to more moderately sized protein datasets.

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.