Understanding the importance of protein structure to nature's routes for divergent evolution in TIM barrel enzymes

Elucidating the Enzymatic Mechanism of Dihydrocoumarin Degradation: Insight into the Functional Evolution of Methyl-Parathion Hydrolase from QM/MM and MM MD Simulations

Methyl-parathion hydrolase (MPH), which evolved from dihydrocoumarin hydrolase, offers one of the most efficient enzymes for the hydrolysis of methyl-parathion. Interestingly, the substrate preference of MPH shifts from the methyl-parathion to the lactone dihydrocoumarin (DHC) after its mutation of five specific residues (R72L, L273F, L258H, T271I, and S193Δ, m5-MPH). Here, extensive QM/MM calculations and MM MD simulations have been used to delve into the structure-function relationship of MPH enzymes and plausible mechanisms for the chemical and nonchemical steps, including the transportation and binding of the substrate DHC to the active site, the hydrolysis reaction, and the product release. The results reveal that the five mutations remodel the active pocket and reposition DHC within the active site, leading to stronger enzyme-substrate interactions. The MM/GBSA-estimated binding free energies are about -20.7 kcal/mol for m5-MPH and -17.1 kcal/mol for wild-type MPH. Furthermore, this conformational adjustment of the protein may facilitate the chemical step of DHC hydrolysis and the product release, although there is a certain influence on the substrate transport. The hydrolytic reaction begins with the nucleophilic attack of the bridging OH- with the energy barriers of 22.0 and 18.0 kcal/mol for the wild-type and m5-MPH enzymes, respectively, which is rate-determining for the entire process. Unraveling these mechanistic intricacies may help in the understanding of the natural evolution of enzymes for diverse substrates and establish the enzyme structure-function relationship.

Identification of Thermal Conduits That Link the Protein-Water Interface to the Active Site Loop and Catalytic Base in Enolase

We report here on the salient role of protein mobility in accessing conformational landscapes that enable efficient enzyme catalysis. We are focused on yeast enolase, a highly conserved lyase with a TIM barrel domain and catalytic loop, as part of a larger study of the relationship of site selective protein motions to chemical reactivity within superfamilies. Enthalpically hindered variants were developed by replacement of a conserved hydrophobic side chain (Leu 343) with smaller side chains. Leu343 is proximal to the active site base in enolase, and comparative pH rate profiles for the valine and alanine variants indicate a role for side chain hydrophobicity in tuning the pK_a of the catalytic base. However, the magnitude of a substrate deuterium isotope effect is almost identical for wild-type (WT) and Leu343Ala, supporting an unchanged rate-determining proton abstraction step. The introduced hydrophobic side chains at position 343 lead to a discontinuous break in both activity and activation energy as a function of side chain volume. Hydrogen-deuterium exchange mass spectrometry (HDX-MS) experiments were performed as a function of time and temperature for WT and Leu343Ala, and provide a spatially resolved map of changes in protein flexibility following mutation. Impacts on protein flexibility are localized to specific networks that arise at the protein-solvent interface and terminate in a loop that has been shown by X-ray crystallography to close over the active site. These interrelated effects are discussed in the context of long-range, solvent-accessible and thermally activated networks that play key roles in tuning the precise distances and interactions among reactants.

Expression, biochemical characterization, and mutation of a water forming NADH: FMN oxidoreductase from Lactobacillus rhamnosus

Enzyme-catalyzed cofactor regeneration is a significant approach to avoid large quantities consumption of oxidized cofactor, which is vital in a variety of bioconversion reactions. NADH: FMN oxidoreductase is an ideal regenerating enzyme because innocuous molecular oxygen is required as an oxidant. But the by-product H₂O₂ limits its further applications at the industrial scale. Here, novel NADH: FMN oxidoreductase (LrFOR) from Lactobacillus rhamnosus comprised of 1146 bp with a predicted molecular weight of 42 kDa was cloned and overexpressed in Escherichia coli BL21 (DE3). Enzyme assay shows that the purified recombinant LrFOR has both the NADPH and NADH oxidation activity. Biochemical characterizations suggested that LrFOR exhibits the specific activity of 39.8 U·mg^-1 with the optimal pH and temperature of 5.6 and 35 °C and produces H₂O instead of potentially harmful peroxide. To further study its catalytic function, a critical Thr29 residue and its six mutants were investigated. Mutants T29G, T29A, and T29D show notable enhancement in activities compared with the wild type. Molecular docking of NADH into wild type and its mutants reveal that a small size or electronegative of residue in position29 could shorten the distance of NADH and FMN, promoting the electrons transfer and resulting in the increased activity. This work reveals the pivotal role of position 29 in the catalytic function of LrFOR and provides effective catalysts in NAD⁺ regeneration.

Taxonomic Landscape of the Dark Proteomes: Whole-Proteome Scale Interplay Between Structural Darkness, Intrinsic Disorder, and Crystallization Propensity

Growth rate of the protein sequence universe dramatically exceeds the speed of expansion for the protein structure universe, generating an immense dark proteome that includes proteins with unknown structure. A whole-proteome scale analysis of 5.4 million proteins from 987 proteomes in the three domains of life and viruses to systematically dissect an interplay between structural coverage, degree of putative intrinsic disorder, and predicted propensity for structure determination is performed. It has been found that Archaean and Bacterial proteomes have relatively high structural coverage and low amounts of disorder, whereas Eukaryotic and Viral proteomes are characterized by a broad spread of structural coverage and higher disorder levels. The analysis reveals that dark proteomes (i.e., proteomes containing high fractions of proteins with unknown structure) have significantly elevated amounts of intrinsic disorder and are predicted to be difficult to solve structurally. Although the majority of dark proteomes are of viral origin, many dark viral proteomes have at least modest crystallization propensity and only a handful of them are enriched in the intrinsic disorder. The disorder, structural coverage, and propensity are mapped for structural determination onto a novel proteome-level sequence similarity network to analyze the interplay of these characteristics in the taxonomic landscape.

Biomimetic water-collecting materials inspired by nature

Nowadays, water shortage is a severe issue all over the world, especially in some arid and undeveloped areas. Interestingly, a variety of natural creatures can collect water from fog, which can provide a source of inspiration to develop novel and functional water-collecting materials. Recently, as an increasingly hot research topic, bioinspired materials with the water collection ability have captured vast scientific attention in both practical applications and fundamental research studies. In this review, we summarize the mechanisms of water collection in various natural creatures and present the fabrications, functions, applications, and new developments of bioinspired materials in recent years. The theoretical basis related to the phenomenon of water collection containing wetting behaviors and water droplet transportations is described in the beginning, i.e., the Young's equation, Wenzel model, Cassie model, surface energy gradient model and Laplace pressure equation. Then, the water collection mechanisms of three typical and widely researched natural animals and plants are discussed and their corresponding bioinspired materials are simultaneously detailed, which are cactus, spider, and desert beetles, respectively. This is followed by introducing another eight animals and plants (butterfly, shore birds, wheat awns, green bristlegrass, the Cotula fallax plant, Namib grass, green tree frogs and Australian desert lizards) that are rarely reported, exhibiting water collection properties or similar water droplet transportation. Finally, conclusions and outlook concerning the future development of bioinspired fog-collecting materials are presented.

Rapid kinetic characterization of glycosyl hydrolases based on oxime derivatization and nanostructure-initiator mass spectrometry (NIMS)

Glycoside hydrolases (GHs) are critical to cycling of plant biomass in the environment, digestion of complex polysaccharides by the human gut microbiome, and industrial activities such as deployment of cellulosic biofuels. High-throughput sequencing methods show tremendous sequence diversity among GHs, yet relatively few examples from the over 150,000 unique domain arrangements containing GHs have been functionally characterized. Here, we show how cell-free expression, bioconjugate chemistry, and surface-based mass spectrometry can be used to study glycoside hydrolase reactions with plant biomass. Detection of soluble products is achieved by coupling a unique chemical probe to the reducing end of oligosaccharides in a stable oxime linkage, while the use of (13)C-labeled monosaccharide standards (xylose and glucose) allows quantitation of the derivatized glycans. We apply this oxime-based nanostructure-initiator mass spectrometry (NIMS) method to characterize the functional diversity of GHs secreted by Clostridium thermocellum, a model cellulolytic organism. New reaction specificities are identified, and differences in rates and yields of individual enzymes are demonstrated in reactions with biomass substrates. Numerical analyses of time series data suggests that synergistic combinations of mono- and multifunctional GHs can decrease the complexity of enzymes needed for the hydrolysis of plant biomass during the production of biofuels.

CLYBL is a polymorphic human enzyme with malate synthase and β-methylmalate synthase activity

CLYBL is a human mitochondrial enzyme of unknown function that is found in multiple eukaryotic taxa and conserved to bacteria. The protein is expressed in the mitochondria of all mammalian organs, with highest expression in brown fat and kidney. Approximately 5% of all humans harbor a premature stop polymorphism in CLYBL that has been associated with reduced levels of circulating vitamin B12. Using comparative genomics, we now show that CLYBL is strongly co-expressed with and co-evolved specifically with other components of the mitochondrial B12 pathway. We confirm that the premature stop polymorphism in CLYBL leads to a loss of protein expression. To elucidate the molecular function of CLYBL, we used comparative operon analysis, structural modeling and enzyme kinetics. We report that CLYBL encodes a malate/β-methylmalate synthase, converting glyoxylate and acetyl-CoA to malate, or glyoxylate and propionyl-CoA to β-methylmalate. Malate synthases are best known for their established role in the glyoxylate shunt of plants and lower organisms and are traditionally described as not occurring in humans. The broader role of a malate/β-methylmalate synthase in human physiology and its mechanistic link to vitamin B12 metabolism remain unknown.

Protein structure quality assessment based on the distance profiles of consecutive backbone Cα atoms

Predicting the three dimensional native state structure of a protein from its primary sequence is an unsolved grand challenge in molecular biology. Two main computational approaches have evolved to obtain the structure from the protein sequence - ab initio/de novo methods and template-based modeling - both of which typically generate multiple possible native state structures. Model quality assessment programs (MQAP) validate these predicted structures in order to identify the correct native state structure. Here, we propose a MQAP for assessing the quality of protein structures based on the distances of consecutive Cα atoms. We hypothesize that the root-mean-square deviation of the distance of consecutive Cα (RDCC) atoms from the ideal value of 3.8 Å, derived from a statistical analysis of high quality protein structures (top100H database), is minimized in native structures. Based on tests with the top100H set, we propose a RDCC cutoff value of 0.012 Å, above which a structure can be filtered out as a non-native structure. We applied the RDCC discriminator on decoy sets from the Decoys 'R' Us database to show that the native structures in all decoy sets tested have RDCC below the 0.012 Å cutoff. While most decoy sets were either indistinguishable using this discriminator or had very few violations, all the decoy structures in the fisa decoy set were discriminated by applying the RDCC criterion. This highlights the physical non-viability of the fisa decoy set, and possible issues in benchmarking other methods using this set. The source code and manual is made available at https://github.com/sanchak/mqap and permanently available on 10.5281/zenodo.7134.

The dynamic determinants of reaction specificity in the IMPDH/GMPR family of (β/α)(8) barrel enzymes

The inosine monophosphate dehydrogenase (IMPDH)/guanosine monophosphate reductase (GMPR) family of (β/α)(8) enzymes presents an excellent opportunity to investigate how subtle changes in enzyme structure change reaction specificity. IMPDH and GMPR bind the same ligands with similar affinities and share a common set of catalytic residues. Both enzymes catalyze a hydride transfer reaction involving a nicotinamide cofactor hydride, and both reactions proceed via the same covalent intermediate. In the case of IMPDH, this intermediate reacts with water, while in GMPR it reacts with ammonia. In both cases, the two chemical transformations are separated by a conformational change. In IMPDH, the conformational change involves a mobile protein flap while in GMPR, the cofactor moves. Thus reaction specificity is controlled by differences in dynamics, which in turn are controlled by residues outside the active site. These findings have some intriguing implications for the evolution of the IMPDH/GMPR family.

Cofactor mobility determines reaction outcome in the IMPDH and GMPR (β-α)8 barrel enzymes

Inosine monophosphate dehydrogenase (IMPDH) and guanosine monophosphate reductase (GMPR) belong to the same structural family, share a common set of catalytic residues and bind the same ligands. The structural and mechanistic features that determine reaction outcome in the IMPDH and GMPR family have not been identified. Here we show that the GMPR reaction uses the same intermediate E-XMP* as IMPDH, but in this reaction the intermediate reacts with ammonia instead of water. A single crystal structure of human GMPR type 2 with IMP and NADPH fortuitously captures three different states, each of which mimics a distinct step in the catalytic cycle of GMPR. The cofactor is found in two conformations: an 'in' conformation poised for hydride transfer and an 'out' conformation in which the cofactor is 6 Å from IMP. Mutagenesis along with substrate and cofactor analog experiments demonstrate that the out conformation is required for the deamination of GMP. Remarkably, the cofactor is part of the catalytic machinery that activates ammonia.

From amino acid to glucosinolate biosynthesis: protein sequence changes in the evolution of methylthioalkylmalate synthase in Arabidopsis

Methylthioalkylmalate synthase (MAM) catalyzes the committed step in the side chain elongation of Met, yielding important precursors for glucosinolate biosynthesis in Arabidopsis thaliana and other Brassicaceae species. MAM is believed to have evolved from isopropylmalate synthase (IPMS), an enzyme involved in Leu biosynthesis, based on phylogenetic analyses and an overlap of catalytic abilities. Here, we investigated the changes in protein structure that have occurred during the recruitment of IPMS from amino acid to glucosinolate metabolism. The major sequence difference between IPMS and MAM is the absence of 120 amino acids at the C-terminal end of MAM that constitute a regulatory domain for Leu-mediated feedback inhibition. Truncation of this domain in Arabidopsis IPMS2 results in loss of Leu feedback inhibition and quaternary structure, two features common to MAM enzymes, plus an 8.4-fold increase in the k(cat)/K(m) for a MAM substrate. Additional exchange of two amino acids in the active site resulted in a MAM-like enzyme that had little residual IPMS activity. Hence, combination of the loss of the regulatory domain and a few additional amino acid exchanges can explain the evolution of MAM from IPMS during its recruitment from primary to secondary metabolism.

Two independent routes of de novo vitamin B6 biosynthesis: not that different after all

Vitamin B6 is well known in its biochemically active form as pyridoxal 5'-phosphate, an essential cofactor of numerous metabolic enzymes. The vitamin is also implicated in numerous human body functions ranging from modulation of hormone function to its recent discovery as a potent antioxidant. Its de novo biosynthesis occurs only in bacteria, fungi and plants, making it an essential nutrient in the human diet. Despite its paramount importance, its biosynthesis was predominantly investigated in Escherichia coli, where it is synthesized from the condensation of deoxyxylulose 5-phosphate and 4-phosphohydroxy-L-threonine catalysed by the concerted action of PdxA and PdxJ. However, it has now become clear that the majority of organisms capable of producing this vitamin do so via a different route, involving precursors from glycolysis and the pentose phosphate pathway. This alternative pathway is characterized by the presence of two genes, Pdx1 and Pdx2. Their discovery has sparked renewed interest in vitamin B6, and numerous studies have been conducted over the last few years to characterize the new biosynthesis pathway. Indeed, enormous progress has been made in defining the nature of the enzymes involved in both pathways, and important insights have been provided into their mechanisms of action. In the present review, we summarize the recent advances in our knowledge of the biosynthesis of this versatile molecule and compare the two independent routes to the biosynthesis of vitamin B6. Surprisingly, this comparison reveals that the key biosynthetic enzymes of both pathways are, in fact, very similar both structurally and mechanistically.

Template-based recognition of protein fold within the midnight and twilight zones of protein sequence similarity

Most homologous pairs of proteins have no significant sequence similarity to each other and are not identified by direct sequence comparison or profile-based strategies. However, multiple sequence alignments of low similarity homologues typically reveal a limited number of positions that are well conserved despite diversity of function. It may be inferred that conservation at most of these positions is the result of the importance of the contribution of these amino acids to the folding and stability of the protein. As such, these amino acids and their relative positions may define a structural signature. We demonstrate that extraction of this fold template provides the basis for the sequence database to be searched for patterns consistent with the fold, enabling identification of homologs that are not recognized by global sequence analysis. The fold template method was developed to address the need for a tool that could comprehensively search the midnight and twilight zones of protein sequence similarity without reliance on global statistical significance. Manual implementations of the fold template method were performed on three folds--immunoglobulin, c-lectin and TIM barrel. Following proof of concept of the template method, an automated version of the approach was developed. This automated fold template method was used to develop fold templates for 10 of the more populated folds in the SCOP database. The fold template method developed three-dimensional structural motifs or signatures that were able to return a diverse collection of proteins, while maintaining a low false positive rate. Although the results of the manual fold template method were more comprehensive than the automated fold template method, the diversity of the results from the automated fold template method surpassed those of current methods that rely on statistical significance to infer evolutionary relationships among divergent proteins.

Mimicking enzyme evolution by generating new (betaalpha)8-barrels from (betaalpha)4-half-barrels

Gene duplication and fusion events that multiply and link functional protein domains are crucial mechanisms of enzyme evolution. The analysis of amino acid sequences and three-dimensional structures suggested that the (betaalpha)8-barrel, which is the most frequent fold among enzymes, has evolved by the duplication, fusion, and mixing of (betaalpha)4-half-barrel domains. Here, we mimicked this evolutionary strategy by generating in vitro (betaalpha)8-barrels from (betaalpha)4-half-barrels that were deduced from the enzymes imidazole glycerol phosphate synthase (HisF) and N'[(5'-phosphoribosyl)formimino]-5-aminoimidazole-4-carboxamide-ribonucleotide isomerase (HisA). To this end, the gene for the C-terminal (betaalpha)4-half-barrel (HisF-C) of HisF was duplicated and fused in tandem to yield HisF-CC, which is more stable than HisF-C. In the next step, by optimizing side-chain interactions within the center of the beta-barrel of HisF-CC, the monomeric and compact (betaalpha)8-barrel protein HisF-C*C was generated. Moreover, the genes for the N- and C-terminal (betaalpha)4-half-barrels of HisF and HisA were fused crosswise to yield the chimeric proteins HisFA and HisAF. Whereas HisFA contains native secondary structure elements but adopts ill-defined association states, the (betaalpha)8-barrel HisAF is a stable and compact monomer that reversibly unfolds with high cooperativity. The results obtained suggest a previously undescribed dimension for the diversification of enzymatic activities: new (betaalpha)8-barrels with novel functions might have evolved by the exchange of (betaalpha)4-half-barrel domains with distinct functional properties.