Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Mutation testing”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 145 records · Page 8
All sourcesNASA NTRSDOE OSTI
DOE OSTI2013

Mutational analysis of the rotavirus NSP4 enterotoxic domain that binds to caveolin-1

Background: Rotavirus (RV) nonstructural protein 4 (NSP4) is the first described viral enterotoxin, which induces early secretory diarrhea in neonatal rodents. Our previous data show a direct interaction between RV NSP4 and the structural protein of caveolae, caveolin-1 (cav-1), in yeast and mammalian cells. The binding site of cav-1 mapped to the NSP4 amphipathic helix, and led us to examine which helical face was responsible for the interaction. Methods: A panel of NSP4 mutants were prepared and tested for binding to cav-1 by yeast two hybrid and direct binding assays. The charged residues of the NSP4 amphipathic helix were changed to alanine (NSP4 46-175 -ala6); and three residues in the hydrophobic face were altered to charged amino acids (NSP4 46-175 -HydroMut). In total, twelve mutants of NSP4 were generated to define the cav-1 binding site. Synthetic peptides corresponding to the hydrophobic and charged faces of NSP4 were examined for structural changes by circular dichroism (CD) and diarrhea induction by a neonatal mouse study. Results: Mutations of the hydrophilic face (NSP4 46-175 -Ala6) bound cav-1 akin to wild type NSP4. In contrast, disruption of the hydrophobic face (NSP 446-175 -HydroMut) failed to bind cav-1. These data suggest NSP4 and cav-1 associate via a hydrophobic interaction. Analyses of mutant synthetic peptides in which the hydrophobic residues in the enterotoxic domain of NSP4 were altered suggested a critical hydrophobic residue. Both NSP4 HydroMut112-140 , that contains three charged amino acids (aa113, 124, 131) changed from the original hydrophobic residues and NSP4 AlaAcidic112-140 that contained three alanine residues substituted for negatively charged (aa114, 125, 132) amino acids failed to induce diarrhea. Whereas peptides NSP4wild type 112 -140 and NSP4 AlaBasic112-140 that contained three alanine substituted for positively charged (aa115, 119, 133) amino acids, induced diarrhea. Conclusions: These data show that the cav-1 binding domain is within the hydrophobic face of the NSP4 amphipathic helix. The integrity of the helical structure is important for both cav-1 binding and diarrhea induction implying a connection between NSP4 functional and binding activities.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2026

A Fluorescence‐Based Transient Expression Assay for the Analysis of Upstream Open Reading Frames in Plants

Upstream open reading frames (uORFs) are regulatory elements present in the 5′ leaders of mRNA that can significantly impact downstream gene expression in eukaryotes. In crop engineering, editing of uORFs can provide an avenue to upregulate expression of native genes without the need to add persistent transgenic copies. Even with genome-wide methods to identify translated uORFs such as ribosome profiling, their functional characterization depends on validation through reporter gene assays and mutagenesis studies. Current screening methods for plants use luciferases or protoplasts to measure differential gene expression between wild-type and mutated transcript leaders, which requires tissue processing and/or substrate addition. Here, we present a time- and cost-efficient alternative to investigate transcript leaders by co-expression of two fluorescent proteins in Nicotiana benthamiana leaf tissue and test our assay on genes involved in photoprotection, editing of which could provide a pathway to increase CO 2 assimilation during sun–shade transitions.

Nicotiana benthamiana↗
DOE OSTI2021

Altered Cardiac Energetics and Mitochondrial Dysfunction in Hypertrophic Cardiomyopathy

Background: Hypertrophic cardiomyopathy (HCM) is a complex disease partly explained by the effects of individual gene variants on sarcomeric protein biomechanics. At the cellular level, HCM mutations most commonly enhance force production, leading to higher energy demands. Despite significant advances in elucidating sarcomeric structure–function relationships, there is still much to be learned about the mechanisms that link altered cardiac energetics to HCM phenotypes. Here, we test the hypothesis that changes in cardiac energetics represent a common pathophysiologic pathway in HCM. Methods: We performed a comprehensive multiomics profile of the molecular (transcripts, metabolites, and complex lipids), ultrastructural, and functional components of HCM energetics using myocardial samples from 27 HCM patients and 13 normal controls (donor hearts). Results: Integrated omics analysis revealed alterations in a wide array of biochemical pathways with major dysregulation in fatty acid metabolism, reduction of acylcarnitines, and accumulation of free fatty acids. HCM hearts showed evidence of global energetic decompensation manifested by a decrease in high energy phosphate metabolites (ATP, ADP, and phosphocreatine) and a reduction in mitochondrial genes involved in creatine kinase and ATP synthesis. Accompanying these metabolic derangements, electron microscopy showed an increased fraction of severely damaged mitochondria with reduced cristae density, coinciding with reduced citrate synthase activity and mitochondrial oxidative respiration. These mitochondrial abnormalities were associated with elevated reactive oxygen species and reduced antioxidant defenses. However, despite significant mitochondrial injury, HCM hearts failed to upregulate mitophagic clearance. Conclusions: Overall, our findings suggest that perturbed metabolic signaling and mitochondrial dysfunction are common pathogenic mechanisms in patients with HCM. These results highlight potential new drug targets for attenuation of the clinical disease through improving metabolic function and reducing mitochondrial injury.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2024

Evaluating the mutagenicity of N-nitrosodimethylamine in 2D and 3D HepaRG cell cultures using error-corrected next generation sequencing

Human liver-derived metabolically competent HepaRG cells have been successfully employed in both two-dimensional (2D) and 3D spheroid formats for performing the comet assay and micronucleus (MN) assay. In the present study, we have investigated expanding the genotoxicity endpoints evaluated in HepaRG cells by detecting mutagenesis using two error-corrected next generation sequencing (ecNGS) technologies, Duplex Sequencing (DS) and High-Fidelity (HiFi) Sequencing. Both HepaRG 2D cells and 3D spheroids were exposed for 72 h to N-nitrosodimethylamine (NDMA), followed by an additional incubation for the fixation of induced mutations. NDMA-induced DNA damage, chromosomal damage, and mutagenesis were determined using the comet assay, MN assay, and ecNGS, respectively. The 72-h treatment with NDMA resulted in concentration-dependent increases in cytotoxicity, DNA damage, MN formation, and mutation frequency in both 2D and 3D cultures, with greater responses observed in the 3D spheroids compared to 2D cells. The mutational spectrum analysis showed that NDMA induced predominantly A:T → G:C transitions, along with a lower frequency of G:C → A:T transitions, and exhibited a different trinucleotide signature relative to the negative control. These results demonstrate that the HepaRG 2D cells and 3D spheroid models can be used for mutagenesis assessment using both DS and HiFi Sequencing, with the caveat that severe cytotoxic concentrations should be avoided when conducting DS. With further validation, the HepaRG 2D/3D system may become a powerful human-based metabolically competent platform for genotoxicity testing.

3D spheroids↗
DOE OSTI2014

Mycobacterium tuberculosis Ser/Thr protein kinase B mediates an oxygen-dependent replication switch

In the majority of cases, Mycobacterium tuberculosis (Mtb) infections are clinically latent, characterized by little or no bacterial replication and drug tolerance. Low oxygen tension is a major host factor inducing bacteriostasis, but the molecular mechanisms driving oxygen-dependent replication are poorly understood. Mtb encodes eleven serine/threonine protein kinases, a family of signaling molecules known to regulate similar replicative adaptations in other bacteria. Here, we tested the role of serine/threonine phosphorylation in the Mtb response to altered oxygen status, using an in vitro model of latency (hypoxia) and reactivation (reaeration). Broad kinase inhibition compromised survival of Mtb in hypoxia. Activity-based protein profiling and genetic mutation identified PknB as the kinase critical for surviving hypoxia. Mtb replication was highly sensitive to changes in PknB levels in aerated culture, and even more so in hypoxia. A mutant overexpressing PknB specifically in hypoxia showed a 10-fold loss in viability in low oxygen conditions. In contrast, chemically reducing PknB activity during hypoxia specifically compromised resumption of growth during reaeration. These data support a model in which PknB activity is reduced to achieve bacteriostasis, and elevated when replication resumes. Together, these data show that phosphosignaling controls replicative transitions associated with latency and reactivation, that PknB is a major regulator of these transitions, and that PknB could provide a highly vulnerable therapeutic target at every step of the Mtb life cycle - active disease, latency, and reactivation.

Environmental Molecular Sciences Laboratory↗
DOE OSTI2020

High-Throughput Functional Genomics for Energy Production

Functional genomics remains a foundational field for establishing genotype-phenotype relationships that enable strain engineering. High-throughput (HTP) methods accelerate the Design-Build-Test-Learn cycle that currently drives synthetic biology towards a forward engineering future. Trackable mutagenesis techniques including transposon insertion sequencing and CRISPR-Cas-mediated genome editing allow for rapid fitness profiling of a collection, or library, of mutants to discover beneficial mutations. Due to the relative speed of these experiments compared to adaptive evolution experiments, iterative rounds of mutagenesis can be implemented for next-generation metabolic engineering efforts to design complex production and tolerance phenotypes. Further, the expansion of these mutagenesis techniques to novel bacteria are opening up industrial microbes that show promise for establishing a bio-based economy.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Extending evolutionary forecasts across bacterial species

Improving evolutionary forecasting requires progressing from studying repeated evolution of a single genotype under identical conditions to formulating broad principles. These principles should enable predictions of how similar species will adapt to similar selective pressures. Evolve-and-resequence experiments with multiple species allow testing forecasts on different biological levels and elucidating the causes for failed predictions. Here, we show that forecasts for adaptation to static culture conditions can be extended to multiple species by testing previous predictions for Pseudomonas syringae and Pseudomonas savastanoi. In addition to sequence divergence, these species differ in their repertoire of biofilm regulatory genes and structural components. Consistent with predictions, both species repeatedly produced biofilm mutants with a wrinkly spreader phenotype. Predominantly, mutations occurred in the wsp operon, with less frequent promoter mutations near uncharacterized diguanylate cyclases. However, mutational patterns differed on the gene level, which was explained by a lack of conservation in relative fitness of mutants between more divergent species. The same mutation was the most frequent for both species suggesting that conserved mutation hotspots can increase parallel evolution. This study shows that evolutionary forecasts can be extended across species, but that differences in the genotype–phenotype–fitness map and mutational biases limit predictability on a detailed molecular level.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2019

Variable genome evolution in fungi after transposon-mediated amplification of a housekeeping gene

Transposable elements (TEs) can be key drivers of evolution, but the mechanisms and scope of how they impact gene and genome function are largely unknown. Previous analyses revealed that TE-mediated gene amplifications can have variable effects on fungal genomes, from inactivation of function to production of multiple active copies. For example, a DNA methyltransferase gene in the wheat pathogen Zymoseptoria tritici (synonym Mycosphaerella graminicola) was amplified to tens of copies, all of which were inactivated by Repeat-Induced Point mutation (RIP) including the original, resulting in loss of cytosine methylation. In another wheat pathogen, Pyrenophora tritici-repentis, a histone H3 gene was amplified to tens of copies with little evidence of RIP, leading to many potentially active copies. To further test the effects of transposon-aided gene amplifications on genome evolution and architecture, the repetitive fraction of the significantly expanded genome of the banana pathogen, Pseudocercospora fijiensis, was analyzed in greater detail. These analyses identified a housekeeping gene, histone H3, which was captured and amplified to hundreds of copies by a hAT DNA transposon, all of which were inactivated by RIP, except for the original. In P. fijiensis the original H3 gene probably was not protected from RIP, but most likely was maintained intact due to strong purifying selection. Comparative analyses revealed that a similar event occurred in five additional genomes representing the fungal genera Cercospora, Pseudocercospora and Sphaerulina. These results indicate that the interplay of TEs and RIP can result in different and unpredictable fates of amplified genes, with variable effects on gene and genome evolution.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2014

New Model for Electron Flow for Sulfate Reduction in Desulfovibrio alaskensis G20

To understand the energy conversion activities of the anaerobic sulfate-reducing bacteria, it is necessary to identify the components involved in electron flow. The importance of the abundant type I tetraheme cytochrome c3 (TpIc3) as an electron carrier during sulfate respiration was questioned by the previous isolation of a null mutation in the gene encoding TpIc3, cycA, in Desulfovibrio alaskensis G20. Whereas respiratory growth of the CycA mutant with lactate and sulfate was little affected, growth with pyruvate and sulfate was significantly impaired. We have explored the phenotype of the CycA mutant through physiological tests and transcriptomic and proteomic analyses. Data reported here show that electrons from pyruvate oxidation do not reach adenylyl sulfate reductase, the enzyme catalyzing the first redox reaction during sulfate reduction, in the absence of either CycA or the type I cytochrome c3:menaquinone oxidoreductase transmembrane complex, QrcABCD. In contrast to the wild type, the CycA and QrcA mutants did not grow with H2 or formate and sulfate as the electron acceptor. Transcriptomic and proteomic analyses of the CycA mutant showed that transcripts and enzymes for the pathway from pyruvate to succinate were strongly decreased in the CycA mutant regardless of the growth mode. Neither the CycA nor the QrcA mutant grew on fumarate alone, consistent with the omics results and a redox regulation of gene expression. We conclude that TpIc3 and the Qrc complex are D. alaskensis components essential for the transfer of electrons released in the periplasm to reach the cytoplasmic adenylyl sulfate reductase and present a model that may explain the CycA phenotype through confurcation of electrons.

Environmental Molecular Sciences Laboratory↗
DOE OSTI2014

New Model for Electron Flow for Sulfate Reduction in Desulfovibrio alaskensis G20

To understand the energy conversion activities of the anaerobic sulfate-reducing bacteria, it is necessary to identify the components involved in electron flow. The importance of the abundant type I tetraheme cytochrome c3 (TpIc3) as an electron carrier during sulfate respiration was questioned by the previous isolation of a null mutation in the gene encoding TpIc3, cycA, in Desulfovibrio alaskensis G20. Whereas respiratory growth of the CycA mutant with lactate and sulfate was little affected, growth with pyruvate and sulfate was significantly impaired. We have explored the phenotype of the CycA mutant through physiological tests and transcriptomic and proteomic analyses. Data reported here show that electrons from pyruvate oxidation do not reach adenylyl sulfate reductase, the enzyme catalyzing the first redox reaction during sulfate reduction, in the absence of either CycA or the type I cytochrome c3:menaquinone oxidoreductase transmembrane complex, QrcABCD. In contrast to the wild type, the CycA and QrcA mutants did not grow with H2 or formate and sulfate as the electron acceptor. Transcriptomic and proteomic analyses of the CycA mutant showed that transcripts and enzymes for the pathway from pyruvate to succinate were strongly decreased in the CycA mutant regardless of the growth mode. Neither the CycA nor the QrcA mutant grew on fumarate alone, consistent with the omics results and a redox regulation of gene expression. We conclude that TpIc3 and the Qrc complex are D. alaskensis components essential for the transfer of electrons released in the periplasm to reach the cytoplasmic adenylyl sulfate reductase and present a model that may explain the CycA phenotype through confurcation of electrons.

Rapp-Giles, Barbara J↗
DOE OSTI2014

An Adaptive Unified Differential Evolution Algorithm for Global Optimization

In this paper, we propose a new adaptive unified differential evolution algorithm for single-objective global optimization. Instead of the multiple mutation strate- gies proposed in conventional differential evolution algorithms, this algorithm employs a single equation unifying multiple strategies into one expression. It has the virtue of mathematical simplicity and also provides users the flexibility for broader exploration of the space of mutation operators. By making all control parameters in the proposed algorithm self-adaptively evolve during the process of optimization, it frees the application users from the burden of choosing appro- priate control parameters and also improves the performance of the algorithm. In numerical tests using thirteen basic unimodal and multimodal functions, the proposed adaptive unified algorithm shows promising performance in compari- son to several conventional differential evolution algorithms.

43 PARTICLE ACCELERATORS↗
DOE OSTI2024

Stress-regulated Arabidopsis GAT2 is a low affinity γ-aminobutyric acid transporter

Abstract The four-carbon non-proteinogenic amino acid γ-aminobutyric acid (GABA) accumulates to high levels in plants in response to various abiotic and biotic stress stimuli, and plays a role in C:N balance, signaling, and as a transport regulator. Expression in Xenopus oocytes and voltage-clamping allowed the characterization of Arabidopsis GAT2 (At5g41800) as a low affinity GABA transporter with a K0.5GABA ~8 mM. l-Alanine and butylamine represented additional substrates. GABA-induced currents were strongly dependent on the membrane potential, reaching the highest affinity and highest transport rates at strongly negative membrane potentials. Mutation of Ser17, previously reported to be phosphorylated in planta, did not result in altered affinity. In a short-term stress experiment, AtGAT2 mRNA levels were up-regulated at low water potential and under osmotic stress (polyethylene glycol and mannitol). Furthermore, AtGAT2 promoter activity was detected in vascular tissues, maturating pollen, and the phloem unloading region of young seeds. Even though this suggested a role for AtGAT2 in long-distance transport and loading of sink organs, under the conditions tested neither AtGAT2-overexpressing plants, atgat2 or atgat1 T-DNA insertion lines, nor atgat1 atgat2 doubleknockout mutants differed from wild-type plants in growth on GABA, amino acid levels, or resistance to salt and osmotic stress.

Plant Sciences↗
DOE OSTI2018

Aberrant alternative splicing of RHOA is associated with loss of its expression and activity in diffuse-type gastric carcinoma cells

RhoA is a member of Rho family small GTPases that regulates diverse cellular functions. Recent large-scale sequencing studies have identified recurrent somatic mutations of RHOA in diffuse-type gastric carcinoma (DGC), indicating that RHOA is a driver of DGC. In this study, we investigated the possible abnormalities of RHOA in a panel of gastric carcinoma (GC) cell lines. Pulldown assay and immunoblot analysis showed that the activity and expression of RhoA were detectable in all GC cell lines tested, except for two DGC cell lines, HSC-59 and GSU. RHOA coding region sequencing revealed that aberrant alternative splicing of RHOA occurred in these cell lines. Quantitative real-time PCR analysis showed that the expression of wild-type RHOA was nearly undetectable, whereas splicing variants were almost exclusively expressed in HSC-59 and GSU cell lines. However, the expression levels of RHOA splicing variants were very low and the corresponding proteins were not detected by immunoblotting. Moreover, the splicing isoforms of RhoA protein were neither efficiently expressed nor activated even if ectopically expressed in cells. These results indicate that aberrant alternative splicing of RHOA results in the loss of its activity and expression in DGC cells.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2018

Investigation of TRPV1 loss-of-function phenotypes in TRPV1 Leu206Stop mice generated by N-ethyl-N-nitrosourea mutagenesis

Highlights: • N-ethyl-N-nitrosourea mutagenesis can generate mouse models for pain research. • TRPV1 Leu206Stop mice closely resemble the phenotype of TRPV1 knockout mice. • Mechanical hypersensitivity is not impaired in TRPV1 Leu206Stop mice. N-ethyl-N-nitrosourea (ENU) random mutagenesis was used to generate a mouse model for the analysis of the transient receptor potential vanilloid 1 (TRPV1) cation channel. A transversion from T→A in exon 4 led to a Leu206Stop mutation generating a loss-of-function mutant. The TRPV1 agonist capsaicin was used to analyze functional and nociceptive parameters in vitro and in vivo in TRPV1 Leu206Stop mice and congenic C3HeB/FeJ controls. Capsaicin-induced [Ca{sup 2+}]{sub i} changes in small diameter DRG neurons were significantly diminished in TRPV1 Leu206Stop mice and administration of capsaicin induced neither hypothermia nor nocifensive behaviour in vivo. TRPV1 Leu206Stop mice were tested in the spinal nerve ligation of mononeuropathic pain and developed mechanical hypersensitivity two weeks after nerve injury. In the open field test, a significant increase in spontaneous locomotion was detected in TRPV1 Leu206Stop mice as compared to wildtype controls. TRPV1 knockout mice have been reported to carry a similar phenotype regarding capsaicin-evoked responses in vitro and in vivo. However, in contrast to TRPV1 Leu206Stop mice, TRPV1 knockout mice did not differ in spontaneous locomotion as compared to congenic C57BL/6 mice, suggesting subtle ENU-dependent or independent strain differences between TRPV1 Leu206Stop mice and their wildtype controls. In summary, these data revealed a target-related (i.e. capsaicin-evoked) phenotype of TRPV1 Leu206Stop mice closely resembling that of published TRPV1 knockout mice. However, since ENU-mutant mice are congenic with the mouse strain initially used in random mutagenesis, direct phenotypic comparison with the respective wildtype controls is possible, and the time-consuming backcrossing in lines with targeted mutations is avoided.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

MO-DE-207B-05: Predicting Gene Mutations in Renal Cell Carcinoma Based On CT Imaging Features: Validation Using TCGA-TCIA Datasets

Purpose: The goal of this work is to investigate the use of contrast enhanced computed tomographic (CT) features for the prediction of mutations of BAP1, PBRM1, and VHL genes in renal cell carcinoma (RCC). Methods: For this study, we used two patient databases with renal cell carcinoma (RCC). The first one consisted of 33 patients from our institution (UT Southwestern Medical Center, UTSW). The second one consisted of 24 patients from the Cancer Imaging Archive (TCIA), where each patient is connected by a unique identi?er to the tissue samples from the Cancer Genome Atlas (TCGA). From the contrast enhanced CT image of each patient, tumor contour was first delineated by a physician. Geometry, intensity, and texture features were extracted from the delineated tumor. Based on UTSW dataset, we completed feature selection and trained a support vector machine (SVM) classifier to predict mutations of BAP1, PBRM1 and VHL genes. We then used TCIA-TCGA dataset to validate the predictive model build upon UTSW dataset. Results: The prediction accuracy of gene expression of TCIA-TCGA patients was 0.83 (20 of 24), 0.83 (20 of 24), and 0.75 (18 of 24) for BAP1, PBRM1, and VHL respectively. For BAP1 gene, texture feature was the most prominent feature type. For PBRM1 gene, intensity feature was the most prominent. For VHL gene, geometry, intensity, and texture features were all important. Conclusion: Using our feature selection strategy and models, we achieved predictive accuracy over 0.75 for all three genes under the condition of using patient data from one institution for training and data from other institutions for testing. These results suggest that radiogenomics can be used to aid in prognosis and used as convenient surrogates for expensive and time consuming gene assay procedures.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2025

Combining Reactive Quantum-Mechanical Molecular-Dynamics Simulations with Mutagenesis, Crystallography, and Enzyme Kinetics to Reveal Plausible Steps of Isocyanide Hydratase Catalysis

A complete understanding of enzyme mechanisms requires atomistic details of chemical reactions. Quantum-based molecular dynamics simulations (QMD) are a potential source of this information, but trade-offs between accuracy and computational cost have limited their use. We previously developed extended Lagrangian Born–Oppenheimer molecular dynamics (XL-BOMD) methods that leverage a negligible compromise in accuracy to substantially decrease the cost of QMD simulations. Here, we develop a reactive QMD approach using the latest XL-BOMD formulation, which enables efficient simulations of highly reactive systems, and use it to investigate mechanisms of intermediate formation in isocyanide hydratase (ICH) catalysis. In QMD simulations, molecular analogs of ICH active site residues reacted with para-nitrophenyl isocyanide, forming a thioimidate. Analysis of simulated atomic configurational and charge dynamics revealed a pathway where protonation of the isocyanide carbon occurs prior to thioimidate formation and suggested a possible role of Asp17 as a proton donor in the early phase of ICH catalysis. To test whether the pathway seen using the reactive QMD approach might be relevant to ICH catalysis, we performed X-ray crystallography and pre-steady-state enzyme kinetics studies of wild-type and D17N mutant ICH. Both the structure and kinetics are sensitive to the D17N mutation in a manner that is consistent with the order of the reaction steps seen in the simulations. Mobile protons play essential roles in many enzymes, yet they are difficult to observe experimentally, making the ordering of proton-dependent steps ambiguous in many enzyme mechanisms. The ability to directly simulate model reactions for the design of experiments that provide information about enzyme mechanisms involving mobile protons demonstrates the significance of our reactive QMD approach and motivates further biological applications.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Understanding Proton Movement in [Fe-Fe] Hydrogenases

Nature uses specialized metalloenzymes to carry out small molecule activation reactions, including CO 2 fixation, O 2 activation, and proton reduction, with unparalleled efficiency, rates, and selectivity. The latter reactions are performed by hydrogenases, protein metallo-complexes that interconverts H 2 to protons and electrons (H 2 oxidation) and the reverse reaction (H 2 production) with incredibly low energy input and amazingly fast kinetics. Reproducing both the activity and efficiency of metalloenzymes in sustainable anthropogenic systems remains one of the “holy grails” of inorganic chemistry. However, identifying the precise molecular components responsible for these desirable properties has been challenging in the natural metalloenzymes, hindering efforts to develop analogous processes in synthetic compounds. Considering the inherent complexity of a metalloenzyme and the many interactions, both strong and weak, that contribute to the function of an enzyme, we have elected to model natural metalloenzymes on a biochemical platform. Towards this end, we have a developed structural, functional, and mechanistic mimic of the [Ni-Fe] hydrogenases within a robust protein scaffold, rubredoxin, to understand the influence of the secondary coordination environment on the metal center. This involved preparing a series of three rubredoxin constructs containing a single point mutation at Val positions and making the NMR chemical shift assignments for the paramagnetic (nickel-substituted) and non-paramagnetic (zinc-substituted) form. These physical studies were complemented with in silico molecular dynamic studies on proteins with [Fe-S] clusters to engineer new proton channels to test in vitro . To assist our search for new natural metalloprotein scaffolds within the vast number of sequenced genomes, we created a neural network-based program to identify proteins with specific metal-binding sites. These computational and physical studies with metalloenzymes provide direct insight into the fundamental chemical principles driving the natural systems and offer design principles for developing catalysts that utilize analogous principles.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2019

Expression of Phosphofructokinase Is Not Sufficient to Enable Embden-Meyerhof-Parnas Glycolysis in Zymomonas mobilis ZM4

Zymomonas mobilis is a bacterium that produces ethanol from glucose at up to 97% of theoretical efficiency on a carbon basis. One factor contributing to the high efficiency of ethanol production is that Z. mobilis has a low biomass yield. The low biomass yield may be caused partly by the low ATP yield of the Entner-Doudoroff (ED) glycolytic pathway used by Z. mobilis , which produces only one ATP per glucose consumed. To test the hypothesis that ATP yield limits biomass yield in Z. mobilis , we attempted to introduce the Embden-Meyerhof-Parnas (EMP) glycolytic pathway (with double the ATP yield) by expressing phosphofructokinase (Pfk I) from Escherichia coli . Expression of Pfk I caused growth inhibition and resulted in accumulation of mutations in the pfkA gene. Co-expression of additional EMP enzymes, fructose bisphosphate aldolase (Fba) and triose phosphate isomerase (Tpi), with Pfk I did not enable EMP flux, and resulted in production of glycerol as a side product. Further analysis indicated that heterologous reactions may have operated in the reverse direction because of native metabolite concentrations. This study reveals how the metabolomic context of a chassis organism influences the range of pathways that can be added by heterologous expression.

59 BASIC BIOLOGICAL SCIENCES↗