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 181 records · Page 10

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↗

A novel root gravitropism mutant of Arabidopsis thaliana exhibiting altered auxin physiology

A root gravitropism mutant was isolated from the DuPont Arabidopsis thaliana T-DNA insertional mutagenesis collection. This mutant has reduced root gravitropism, hence the name rgr1. Roots of rgr1 are shorter than those of wild-type, and they have reduced lateral root formation. In addition, roots of rgr1 coil clockwise on inclined agar plates, unlike wild-type roots which grow in a wavy pattern. The rgr1 mutant has increased resistance, as measured by root elongation, to exogenously applied auxins (6-fold to indole-3-acetic acid, 3-fold to 2,4-dichlorophenoxyacetic acid, and 2-fold to napthyleneacetic acid). It is also resistant to polar auxin transport inhibitors (2-fold to triiodobenzoic acid and 3- to 5-fold to napthylphthalamic acid). The rgr1 mutant does not appear to be resistant to other plant hormone classes. When grown in the presence of 10(-7) M 2,4-dichlorophenoxyacetic acid, rgr1 roots have fewer root hairs than wild type. All these rgr1 phenotypes are Mendelian recessives. Complementation tests indicate that rgr1 is not allelic to previously characterized agravitropic or auxin-resistant mutants. The rgr1 locus was mapped using visible markers to 1.4 +/- 0.6 map units from the CH1 locus at 1-65.4. The rgr1 mutation and the T-DNA cosegregate, suggesting that rgr1 was caused by insertional gene inactivation.

NASA Discipline Plant Biology↗

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↗

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↗

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↗

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↗

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↗

DNA Repair Domain Modeling Can Predict Cell Death and Mutation Frequency for Wide Range Spectrum of Radiation

Exploration missions to Mars and other destinations raise many questions about the health of astronauts. The continuous exposure of astronauts to galactic cosmic rays is one of the main concerns for long-term missions. Cosmic ionizing radiations are composed of different ions of various charges and energies notably, highly charged energy (HZE) particles. The HZE particles have been shown to be more carcinogenic than low-LET radiation, suggesting the severity of chromosomal aberrations induced by HZE particles is one possible explanation. However, most mathematical models predicting cell death and mutation frequency are based on directly fitting various HZE dose response and are in essence empirical approaches. In this work, we assume a simple biological mechanism to model DNA repair and use it to simultaneously explain the low- and high-LET response using the exact same fitting parameters. Our work shows that the geometrical position of DNA repair along tracks of heavy ions are sufficient to explain why high-LET particles can induce more death and mutations. Our model is based on assuming DNA double strand breaks (DSBs) are repaired within repair domain, and that any DSBs located within the same repair domain cluster into one repair unit, facilitating chromosomal rearrangements and increasing the probability of cell death. We introduced this model in 2014 using simplified microdosimetry profiles to predict cell death. In this work, we collaborated with NASA Johnson Space Center to generate more accurate microdosimetry profiles derived by Monte Carlo techniques, taking into account track structure of HZE particles and simulating DSBs in realistic cell geometry. We simulated 224 data points (D, A, Z, E) with the BDSTRACKS model, leading to a large coverage of LET from ~10 to 2,400 keV/μm. This model was used to generate theoretical RBE for various particles and energies for both cell death and mutation frequencies. The RBE LET dependence is in agreement with experimental data known in human and murine cells. It suggests that cell shape and its orientation with respect to the HZE particle beam can modify the biological response to radiation. Such discovery will be tested experimentally and, if proven accurate, will be another strong supporting evidence for DNA repair domains and their critical role in interpreting cosmic radiation sensitivity.

Viger, Louise↗

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↗

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↗

NASA's International Space Station: A Testbed for Planetary Protection Protocol Development

Wherever humans go, they inevitably carry along the critters that live in and on them. Conventional wisdom has long held that it is unlikely those critters could survive the space environment, but in 2007 some microscopic aquatic animals called Tardigrades survived exposure to space and in 2008 Cyanobacteria lived for 548 days outside the ISS. Unlike the Mars rovers that were cleaned once and sent on their way, crew members will provide a constantly regenerating contaminant source. Are we prepared to certify that we can meet forward contamination protocols as we search for life at new destinations? What about the organisms we might reasonably expect a crewed spacecraft to leak or vent? Do we even know what they are? How long might our tiny hitch-hikers survive in close proximity to a warm spacecraft that periodically leaks/vents water or oxygen and how might they mutate with long-duration exposure? How will these contaminants migrate from their source in conditions encountered in space or on other planetary surfaces? This project aims to answer some of these questions by bringing together key stakeholder communities to develop a human forward contamination test, analysis, and integration plan. A system engineering approach to identify the experiments, analysis, and modeling needed to develop the contamination control protocols required will be used as a roadmap to integrate the many different parts of this problem - from launch to landing, living, and working on another planetary surface.

Bell, M. S.↗

Nasa's International Space Station: A Testbed for Planetary Protection Protocol Development

Wherever humans go, they inevitably carry along the critters that live in and on them. Conventional wisdom has long held that it is unlikely those critters could survive the space environment, but in 2007 some microscopic aquatic animals called Tardigrades survived exposure to space and in 2008 Cyanobacteria lived for 548 days outside the ISS. Unlike the Mars rovers that were cleaned once and sent on their way, crew members will provide a constantly regenerating contaminant source. Are we prepared to certify that we can meet forward contamination protocols as we search for life at new destinations? What about the organisms we might reasonably expect a crewed spacecraft to leak or vent? Do we even know what they are? How long might our tiny hitch-hikers survive in close proximity to a warm spacecraft that periodically leaks/vents water or oxygen and how might they mutate with long-duration exposure? How will these contaminants migrate from their source in conditions encountered in space or on other planetary surfaces? This project aims to answer some of these questions by bringing together key stakeholder communities to develop a human forward contamination test, analysis, and integration plan. A system engineering approach to identify the experiments, analysis, and modeling needed to develop the contamination control protocols required will be used as a roadmap to integrate the many different parts of this problem - from launch to landing, living, and working on another planetary surface.

Bell, M. S.↗

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↗

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↗

Selective norepinephrine reuptake inhibition as a human model of orthostatic intolerance

BACKGROUND: Observations in patients with functional mutations of the norepinephrine transporter (NET) gene suggest that impaired norepinephrine uptake may contribute to idiopathic orthostatic intolerance. METHODS AND RESULTS: We studied the effect of the selective NET blocker reboxetine and placebo in a randomized, double-blind, crossover fashion on cardiovascular responses to cold pressor testing, handgrip testing, and a graded head-up tilt test (HUT) in 18 healthy subjects. In a subset, we determined isoproterenol and phenylephrine sensitivities. Subjects ingested 8 mg reboxetine or placebo 12 hours and 1 hour before testing. In the supine position, heart rate was 65+/-2 bpm with placebo and 71+/-3 bpm with reboxetine. At 75 degrees HUT, heart rate was 84+/-3 and 119+/-4 bpm with placebo and with reboxetine (P<0.0001). Mean arterial pressure was 85+/-2 with placebo and 91+/-2 mm Hg with reboxetine while supine (P<0.01) and 88+/-2 mm Hg and 90+/-3 mm Hg at 75 degrees HUT. Blood pressure responses to cold pressor and handgrip testing were attenuated with reboxetine. Reboxetine increased the sensitivity to the chronotropic effect of isoproterenol and the pressor effect of phenylephrine. Vasovagal reactions occurred in 9 subjects on placebo and in 1 subject on reboxetine. CONCLUSIONS: Selective NET blockade creates a phenotype that resembles idiopathic orthostatic intolerance. This observation supports the hypothesis that disordered norepinephrine uptake mechanisms can contribute to human cardiovascular disease. Our study also suggests that NET inhibition might be useful in preventing vasovagal reactions.

Non-NASA Center↗

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↗

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↗

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↗