Search NASA⌕ Search

SEARCH · Search NASA

Results for “in silico 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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Computationally evaluating high-yield metabolites for sustainable aviation fuel (SAF) using machine learning

The computational tool described in this report helps identify promising biological pathways that produce SAF platform molecules (either a drop-in SAF, or a precursor that can be easily converted to a drop-in SAF). The workflow the computational tool follows first identifies possible biological pathways from a user-defined metabolite. These pathways may, or may not lead to a SAF platform molecule, thus the second step involves insilico testing of the end product of each pathway to assess whether it is, or is not, a SAF platform molecule. The identification of biological pathways performed in the first step is facilitated by linking the metabolite to a biological reaction database. Pathways are found by identifying pathways in the reaction database that include the metabolite. The computational tool includes an alternative way to find pathways. The alternative way develops a Flux Balanced Analysis (FBA), and modifying the FBA to include reactions that transform the metabolite. These modifications serve as a basis for understanding, in a semi-quantitative way, if there is an increase in the flux to desirable products. The second step, in silico testing of the end-products, is accomplished by estimating key physical properties relevant to SAF. When good models are available, we have integrated those models into the computational tool. In a few instances, we have developed our own models. In all instances, we have validated the models against available measured data. Finally, we have evaluated the effectiveness of our computational tool by genetically engineering Rhodosporidium toruloides. Validation occurred without the use of a FBA, and further validation is required.

09 BIOMASS FUELS↗
DOE OSTI2021

A new primer set for Clade I nosZ that recovers genes from a broader range of taxa

Denitrification is an important global N cycle process. The gene encoding NosZ that converts nitrous oxide (N 2 O) to N 2 has been widely used as a biomarker to study denitrifying communities. However, conventional PCR primers target a limited range of the genetically diverse Clade I nosZ, and the amplicons are too long for sequencing on current NGS platforms. To address these issues, here we developed a new PCR primer set that amplifies a 355-bp region of Clade I nosZ and captures broader taxonomic coverage than conventional primers in in silico tests. When compared with the widely used nosZF_nosZR_Rich_2003 set using the same soil samples and the same sequencing depth, the new set retrieved genes from four times more unique species, with consistently higher general diversity-based metrics. The new primer set performed well with different sequencing platforms (Ion Torrent and Illumina), and among a wide variety of soils from polar to tropical, desert to agricultural, and surface to a very low biomass subsoil, with significant differences in denitrifying community diversity and composition. This new primer set for Clade I together with the primers recently reported for Clade II by Chee-Sanford et al. provides a more comprehensive assessment of denitrifier gene hosts, their ecological patterns, and the degree of novelty in retrieved gene sequences.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

Enantioselective ecotoxicity of promethazine in two freshwater organisms: daphnia (Daphnia magna) and zebrafish ( Danio rerio )

Abstract Chiral pharmaceuticals, racemic or enantiomerically pure forms and their metabolites, can reach aquatic ecosystems via wastewater effluents (inefficient treatment operations) or by direct human disposal. They may negatively affect nontarget organisms even at low environmental concentrations. To make an accurate risk evaluation, the (eco)toxicity of both enantiomers needs to be assessed. Promethazine (PMZ) is a chiral antihistamine that has been detected in aquatic ecosystems owing to its high consumption. Promethazine undergoes metabolism in the liver, producing chiral metabolites such as promethazine sulfoxide (PMZSO) and N-desmethylpromethazine (DMPMZ) that reach water bodies. However, knowledge regarding the enantioselective toxicity of PMZ and its metabolites on aquatic organisms is missing. This study aimed to explore the potential enantioselective toxicity of PMZ and its metabolites on two relevant freshwater organisms, daphniid and fish, representing different trophic levels. The half maximal effect concentrations (EC50s) in Daphnia magna of PMZ, DMPMZ, and PMZSO were 2.33, 2.31, > 4 mg L−1, respectively, > 4 and 2.50 mg L−1 for (R) and (S)-PMZ, respectively, and > 4 mg L−1 for the enantiomers of DMPMZ and PMZSO. In studies involving zebrafish, Danio rerio, (R, S)-PMZ showed a median lethal concentration (LC50) of .72 mg L−1, and specific assays revealed that (R)-PMZ exhibited more pronounced adverse effects on larvae at the embryonic, morphological, and biochemical level than the racemate and (S)-PMZ. Toxicity and potential bioaccumulation of these compounds in daphniids and fish were also conducted using in silico tests through proprietary software. The results revealed a concordance between the experimental and predicted EC50 and LC50 values in both species.

Coelho, Maria Miguel↗
DOE OSTI2021

CATMoS: Collaborative Acute Toxicity Modeling Suite

Background: Humans are exposed to tens of thousands of chemical substances that need to be assessed for their potential toxicity. Acute systemic toxicity testing serves as the basis for regulatory hazard classification, labeling, and risk management. However, it is cost- and time-prohibitive to evaluate all new and existing chemicals using traditional rodent acute toxicity tests. In silico models built using existing data facilitate rapid acute toxicity predictions without using animals. Objectives: The U.S. Interagency Coordinating Committee on the Validation of Alternative Methods Acute Toxicity Workgroup organized an international collaboration to develop in silico models for predicting acute oral toxicity based on five different endpoints: LD50 value, U.S. Environmental Protection Agency hazard categories, Globally Harmonized System for Classification and Labelling hazard categories, very toxic chemicals (LD50 =50 mg/kg), and non-toxic chemicals (LD50 >2000 mg/kg). Methods: An acute oral toxicity data inventory for 11,992 chemicals was compiled, split into training and evaluation sets, and made available to 35 participating international research groups that submitted a total of 139 predictive models. Predictions that fell within the applicability domains of the submitted models were evaluated using external validation sets. These were then combined into consensus models to leverage strengths of individual approaches. Results: The resulting consensus predictions, which leverage the collective strengths of each individual model, form the Collaborative Acute Toxicity Modeling Suite (CATMoS). CATMoS demonstrated high performance in terms of accuracy and robustness when compared to in vivo results. Discussion: CATMoS is being evaluated by regulatory agencies for its utility and applicability as a potential replacement for in vivo rat acute oral toxicity studies. CATMoS predictions for over 800,000 chemicals have been made available via the NTP’s Integrated Chemical Environment. The models are also implemented in a free, standalone open-source tool, OPERA, which allows predictions of new and untested chemicals to be made.

63 RADIATION, THERMAL, AND OTHER ENVIRON. POLLUTAN↗
DOE OSTI2026

State of the art, gaps, and prospects in fusion materials theory and modelling

Advancing the theory and simulation of materials for fusion applications remains a key component of global roadmaps aimed at delivering much-needed fusion power. Especially as the drive for commercial application increases, prototypes must be designed against radiation damage before the relevant experimental data can be collected and cost reductions that are possible by testing materials in silico become even more important. Here, we summarise the state of the art as it emerged during the 7 th Fusion Materials Theory & Modelling Workshop that took place in 2024, with the aim to highlight present gaps and future directions for the fusion materials modelling community. Of particular interest were the effects of transmutations, chemical complexity with the development of novel alloys and interatomic potentials, advancements in modelling high-dose microstructures, comparison with experimental data and multiscale models for structural assessment relying on high-performance computing and virtual reality.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Robust deep learning–based protein sequence design using ProteinMPNN

Although deep learning has revolutionized protein structure prediction, almost all experimentally characterized de novo protein designs have been generated using physically based approaches such as Rosetta. Here, we describe a deep learning–based protein sequence design method, ProteinMPNN, that has outstanding performance in both in silico and experimental tests. On native protein backbones, ProteinMPNN has a sequence recovery of 52.4% compared with 32.9% for Rosetta. The amino acid sequence at different positions can be coupled between single or multiple chains, enabling application to a wide range of current protein design challenges. We demonstrate the broad utility and high accuracy of ProteinMPNN using x-ray crystallography, cryo–electron microscopy, and functional studies by rescuing previously failed designs, which were made using Rosetta or AlphaFold, of protein monomers, cyclic homo-oligomers, tetrahedral nanoparticles, and target-binding proteins.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Harnessing the power of ab initio calculations, distributed computing and machine learning to efficiently locate extreme molecules for use in carbon-based solar cells (Final Technical Report)

The use of high-throughput virtual screening (HTVS) tools is a powerful tool to expedite the materials discovery of commercially relevant materials. In previous years, our group has developed a molecular discovery platform to generate libraries in order to obtain suitable candidates for different applications, starting from the Harvard Clean Energy Project. This platform is suitable to test in-silico on traditional supercomputing clusters and shared resources, for example, in the IBM World Community Grid. In this project, we used the molecular discovery platform to create and screen a library of candidates of organic photovoltaic (OPVs) molecules. Based on a set of candidates created with combinations of molecular moieties, we were able to filter, by conformation stability, the energy of electronic orbitals and approximated power conversion efficiencies (PCE). To improve the predictions of orbital energies calculated and the PCEs, we used Gaussian Process regression and two sets of molecules. These sets correspond to electronic structure calculations of a higher level of theory and experimental PCE values, respectively. Finally, we selected a subset of the best candidates (molecules with a PCE higher than 10%) to understand its absorbance properties with TD-DFT. This project has demonstrated the capabilities of our molecular discovery platform for HTVS. Finally, machine learning can help us to introduce more complex effects included in bulk conditions and computational intensive calculations on models.

14 SOLAR ENERGY↗
DOE OSTI2022

A novel approach to build algal consortia for sustainable biomass production

In the last decade, microalgae have reemerged as a feedstock for biofuels and a diverse suite of bioproducts. Yet, considerable challenges must be overcome before algal biofuels and bioproducts become technoeconomically viable. At present, single algal strain selected for particular phenotypic traits, such as maximum specific growth rate or lipid content, are commonly scaled for cultivation in open, outdoor raceway ponds due to the low capital costs of these systems. Although this monoculture approach may maximize the production of end products, monocultures are particularly susceptible to crashes associated with environmental and biological variability. An approach that has been proposed to generate more productive and stable microalgal crops is the use of eco-engineered communities, or consortia. Yet, attempts to construct productive consortia have not been consistently successful. We argue that failures stem from the lack of an eco-engineering approach to design species combinations. Here, we used an in silico method to build consortia before testing their performance against monocultures. Focusing on consortia of Nannochloropsis and Microchloropsis, we measured growth of strains along gradients of light, temperature, and salinity and used a functional dispersion approach to generate over 8000 functionally-diverse consortia combinations. We tested the 50 most functionally diverse consortia in a laboratory experiment and found that consortia overwhelmingly outperformed monocultures. Indeed, overyielding (OY) and a positive net biodiversity effect (NBE) was found respetively in 8%-86% and 88-92% of consortia combinations over the different experimental phases. To our knowledge, this is the first application of an in silico approach to design functionally diverse consortia before laboratory and field testing. Furthermore, our results highlight the importance of employing a functional diversity approach for consortia design.

09 BIOMASS FUELS↗
DOE OSTI2024

A goldilocks computational protocol for inhibitor discovery targeting DNA damage responses including replication-repair functions

While many researchers can design knockdown and knockout methodologies to remove a gene product, this is mainly untrue for new chemical inhibitor designs that empower multifunctional DNA Damage Response (DDR) networks. Here, we present a robust Goldilocks (GL) computational discovery protocol to efficiently innovate inhibitor tools and preclinical drug candidates for cellular and structural biologists without requiring extensive virtual screen (VS) and chemical synthesis expertise. By computationally targeting DDR replication and repair proteins, we exemplify the identification of DDR target sites and compounds to probe cancer biology. Our GL pipeline integrates experimental and predicted structures to efficiently discover leads, allowing early-structure and early-testing (ESET) experiments by many laboratories. By employing an efficient VS protocol to examine protein-protein interfaces (PPIs) and allosteric interactions, we identify ligand binding sites beyond active sites, leveraging in silico advances for molecular docking and modeling to screen PPIs and multiple targets. A diverse 3,174 compound ESET library combines Diamond Light Source DSI-poised, Protein Data Bank fragments, and FDA-approved drugs to span relevant chemotypes and facilitate downstream hit evaluation efficiency for academic laboratories. Two VS per library and multiple ranked ligand binding poses enable target testing for several DDR targets. This GL library and protocol can thus strategically probe multiple DDR network targets and identify readily available compounds for early structural and activity testing to overcome bottlenecks that can limit timely breakthrough drug discoveries. By testing accessible compounds to dissect multi-functional DDRs and suggesting inhibitor mechanisms from initial docking, the GL approach may enable more groups to help accelerate discovery, suggest new sites and compounds for challenging targets including emerging biothreats and advance cancer biology for future precision medicine clinical trials.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

Identification and overexpression of endogenous transcription factors to enhance lipid accumulation in the biotechnologically relevant species Chlamydomonas pacifica

Sustainable low-carbon energy solutions are critical to mitigating global carbon emissions. Algae-based platforms offer potential by converting carbon dioxide into valuable products while aiding carbon sequestration. However, scaling algae cultivation faces challenges like contamination in outdoor systems. Previously, our lab evolved Chlamydomonas pacifica, an extremophile green alga, which tolerates high temperature, pH, salinity, and light, making it ideal for large-scale bioproduct production, including biodiesel. Here, we enhanced lipid accumulation in evolved C. pacifica by identifying and overexpressing key endogenous transcription factors through genome-wide in-silico analysis and in-vivo testing. These factors include Lipid Remodeling Regulator 1 (CpaLRL1), Nitrogen Response Regulator 1 (CpaNRR1), Compromised Hydrolysis of Triacylglycerols 7 (CpaCHT7), and Phosphorus Starvation Response 1 (CpaPSR1). Under nitrogen deprivation, CpaLRL1, CpaNRR1, and CpaCHT7 overexpression enhanced lipid accumulation compared to wild-type. However, CpaPSR1 increased lipid accumulation compared to wild-type in normal media and did not increase further under nitrogen deprivation, highlighting the difference in function based on media conditions. Notably, lipid analysis of CpaPSR1 under normal media conditions revealed a 2.4-fold increase in triglycerides (TAGs) compared to the wild-type, highlighting its potential for biodiesel production. This approach provides a framework for transcription factor-focused metabolic engineering in algae, advancing bioenergy and biomaterial production.

Biofuels↗
DOE OSTI2024

Real-time PCR assays that detect genes for botulinum neurotoxin A–G subtypes

The role of Real-Time PCR assays for surveillance and rapid screening for pathogens is garnering more and more attention because of its versatility and ease of adoption. The goal of this study was to design, test, and evaluate Real-Time TaqMan PCR assays for the detection of botulinum neurotoxin (bont/A-G) genes from currently recognized BoNT subtypes. Assays were computationally designed and then laboratory tested for sensitivity and specificity using DNA preparations containing bont genes from 82 target toxin subtypes, including nine bivalent toxin types; 31 strains representing other clostridial species; and an extensive panel that consisted of DNA from a diverse set of prokaryotic (bacterial) and eukaryotic (fungal, protozoan, plant, and animal) species. In addition to laboratory testing, the assays were computationally evaluated using in silico analysis for their ability to detect bont gene sequences from recently identified toxin subtypes. Seventeen specific assays (two for each of the bont/C, bont/D, bont/E, and bont/G subtypes and three for each of the bont/A, bont/B, and bont/F subtypes) were designed and evaluated for their ability to detect bont genes encoding multiple subtypes from all seven serotypes. These assays could provide an additional tool for the detection of botulinum neurotoxins in clinical, environmental and food samples that can complement other existing methods used in clinical diagnostics, regulatory, public health, and research laboratories.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

Building a predictive model for polycyclic aromatic hydrocarbon dosimetry in organotypically cultured human bronchial epithelial cells using benzo[ a ]pyrene

The airway epithelium is a primary route of exposure for inhaled toxicants, and organotypic culture models represent an important advancement for toxicity testing compared to simple in vitro models that may lack metabolic capability and multicellular structure/communication associated with the bronchial epithelium in vivo. A quantitative understanding of chemical dosimetry is key for interpreting and extrapolating study results; however, dosimetry is understudied in organotypic models limiting ability to predict toxicity. We developed a dosimetry model for primary human bronchial epithelial cells (HBECs) cultured at the air-liquid interface (ALI) using benzo[a]pyrene (BAP), a representative polycyclic aromatic hydrocarbon. Dose and time course evaluation of metabolite formation and enzyme activity and expression were utilized to parameterize a cellular dosimetry model to improve the utility of ALI-HBECs for assessing chemical risk. Dosimetry analysis demonstrated absorption of BAP into cells and an increase in Phase 1 and 2 metabolites over time that correlated with regulation of metabolizing enzymes. BAP was cleared from cells by 48 h after exposure, and the primary metabolites generated in ALI-HBECs were BAP-3-phenol, BAP-4,5-dihydrodiol, BAP-7,8-dihydrodiol, BAP-9,10-dihydrodiol, BAP-7,8,9,10-tetrol, BAP-3-phenol-glucuronide, BAP-4,5-dihydrodiol-glucuronide, and BAP-9,10-dihydrodiol-glucuronide. The resulting dosimetry model described BAP and 7,8-dihydrodiol toxicokinetics in ALI-HBECs and suggested active excretion of 7,8-dihydrodiol. Overall, this study demonstrates metabolic competency of ALI-HBECs for BAP metabolism, demonstrates the usefulness of complex in vitro systems for human-relevant toxicity data, and exhibits how in silico models can be utilized for understanding the dosimetry of test compounds to aid in in vitro to human extrapolation of toxicity data for risk assessments.

Benzo[a]pyrene↗
DOE OSTI2026

Paleotribological models using preserved fossil tissue properties reveal functional significance of Eurasian mammoth dental evolution

Eurasian mammoths (Mammuthus) underwent substantial modifications in molar morphology as later-diverging species evolved progressively thinner enamel and increased enamel crest complexity. These features have been hypothesized to reduce whole-tooth wear and extend dental longevity as increasingly graze-dominated diets evolved within the lineage. This hypothesis has yet to be directly tested. Here, in this study, we developed an in-silico wear model using experimentally derived wear rates from fossil and extant proboscidean dental tissues. The models revealed that shifts in tissue topology do not affect whole-tooth wear rate, as inverse trends in lamellar frequency and enamel thickness preserve a consistent surface enamel area fraction; the determining factor of wear. Rather, topological shifts produce a wear-emergent secondary occlusal surface with greater numbers of triturating crests that create a regular, low-relief, file-like shearing pavement. These changes in occlusal architecture likely directly impacted the mastication capacity of Mammuthus dentitions, facilitating their dietary expansion to incorporate fibrous, lower-nutrient graze.

Dental wear↗
DOE OSTI2020

In Silico Guidance for In Vitro Androgen and Glucocorticoid Receptor ToxCast Assays

Molecular initiating events (MIEs) are key events in adverse outcome pathways (AOPs) that link molecular chemistry to target biology. As they are based in chemistry, these interactions are excellent targets for computational chemistry approaches to in silico modelling. In this work, we aim to link ligand chemical structure to MIEs for androgen receptor (AR) and glucocorticoid receptor (GR) binding using ToxCast data. This has been done using an automated computational algorithm to perform maximal common substructure searches on chemical binders for each target from the ToxCast dataset. The models developed show a high level of accuracy, correctly assigning 87.20% of AR binders and 96.81% of GR binders in a 25% test set using holdout cross-validation. The 2D structural alerts developed can be used as in silico models to predict these MIEs, and as guidance for in vitro ToxCast assays to confirm hits. These models can target such experimental work, reducing the number of assays to be performed to gain required toxicological insight. Development of these models has also allowed some structural alerts to be identified as predictors for agonist or antagonist behavior at the receptor target. This work represents a first step in using computational methods to guide and target experimental approaches.

Allen, Timothy H.↗
DOE OSTI2023

Drug discovery efforts at George Mason University

With over 39,000 students, and research expenditures in excess of $200 million, George Mason University (GMU) is the largest R1 (Carnegie Classification of very high research activity) university in Virginia. Mason scientists have been involved in the discovery and development of novel diagnostics and therapeutics in areas as diverse as infectious diseases and cancer. Below are highlights of the efforts being led by Mason researchers in the drug discovery arena. To enable targeted cellular delivery, and non-biomedical applications, Veneziano and colleagues have developed a synthesis strategy that enables the design of self-assembling DNA nanoparticles (DNA origami) with prescribed shape and size in the 10 to 100 nm range. The nanoparticles can be loaded with molecules of interest such as drugs, proteins and peptides, and are a promising new addition to the drug delivery platforms currently in use. The investigators also recently used the DNA origami nanoparticles to fine tune the spatial presentation of immunogens to study the impact on B cell activation. These studies are an important step towards the rational design of vaccines for a variety of infectious agents. To elucidate the parameters for optimizing the delivery efficiency of lipid nanoparticles (LNPs), Buschmann, Paige and colleagues have devised methods for predicting and experimentally validating the pKa of LNPs based on the structure of the ionizable lipids used to formulate the LNPs. These studies may pave the way for the development of new LNP delivery vehicles that have reduced systemic distribution and improved endosomal release of their cargo post administration. To better understand protein-protein interactions and identify potential drug targets that disrupt such interactions, Luchini and colleagues have developed a methodology that identifies contact points between proteins using small molecule dyes. The dye molecules noncovalently bind to the accessible surfaces of a protein complex with very high affinity, but are excluded from contact regions. When the complex is denatured and digested with trypsin, the exposed regions covered by the dye do not get cleaved by the enzyme, whereas the contact points are digested. The resulting fragments can then be identified using mass spectrometry. The data generated can serve as the basis for designing small molecules and peptides that can disrupt the formation of protein complexes involved in disease processes. For example, using peptides based on the interleukin 1 receptor accessory protein (IL-1RAcP), Luchini, Liotta, Paige and colleagues disrupted the formation of IL-1/IL-R/IL-1RAcP complex and demonstrated that the inhibition of complex formation reduced the inflammatory response to IL-1B. Working on the discovery of novel antimicrobial agents, Bishop, van Hoek and colleagues have discovered a number of antimicrobial peptides from reptiles and other species. DRGN-1, is a synthetic peptide based on a histone H1-derived peptide that they had identified from Komodo Dragon plasma. DRGN-1 was shown to disrupt bacterial biofilms and promote wound healing in an animal model. The peptide, along with others, is being developed and tested in preclinical studies. Other research by van Hoek and colleagues focuses on in silico antimicrobial peptide discovery, screening of small molecules for antibacterial properties, as well as assessment of diffusible signal factors (DFS) as future therapeutics. The above examples provide insight into the cutting-edge studies undertaken by GMU scientists to develop novel methodologies and platform technologies important to drug discovery.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Theoretical modeling of hepatitis C acute infection in liver-humanized mice support pre-clinical assessment of candidate viruses for controlled-human-infection studies

Designing and carrying out a controlled human infection (CHI) model for hepatitis C virus (HCV) is critical for vaccine development. However, key considerations for a CHI model protocol include understanding of the earliest viral-host kinetic events during the acute phase and susceptibility of the viral isolate under consideration for use in the CHI model to antiviral treatment before any infections in human volunteers can take place. Humanized mouse models lack adaptive immune responses but provide a unique opportunity to obtain quantitative understanding of early HCV kinetics and develop mathematical models to further understand viral and innate immune response dynamics during acute HCV infection. We show that the models reproduce the measured HCV kinetics in humanized mice, which are consistent with early acute HCV-host dynamics in immunocompetent chimpanzees. Our findings suggest that humanized mice are well-suited to support development of a CHI model. In-silico and in-vivo modeling estimates provide a starting point to characterize candidate viruses for testing in CHI model studies.

Agent-based modeling↗
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 OSTI2020

Density Functionals for Hydrogen Storage: Defining the H2Bind275 Test Set with Ab Initio Benchmarks and Assessment of 55 Functionals

Efficient and high-capacity storage materials are indispensable for a hydrogen-based economy. In silico tools can accelerate the process of discovery of new adsorbent materials with optimal hydrogen adsorption enthalpies. Density functional theory is well-poised to become a very useful tool for enabling high-throughput screening of potential materials. In this work, we have identified density functional approximations that provide good performance for hydrogen binding applications following a two-pronged approach. First, we have compiled a data set (H2Bind275) that comprehensively represents the hydrogen binding problem capturing the chemical and mechanistic diversity in the binding sites encountered in hydrogen storage materials. We have also computed reference interaction energies for this data set using coupled-cluster theory. Second, we have assessed the performance of 55 density functional approximations for predicting H 2 interaction energies and have identified two hybrid density functionals (ωB97X-V and ωB97M-V), two double hybrid density functionals (DSD-PBEPBE-D3(BJ) and PBE0-DH), and one semilocal density functional (B97M-V) as the best performing ones. We have recommended the addition of empirical dispersion corrections to systematically underbinding density functionals such as revPBE, BLYP, and B3LYP for improvements in performance at negligible additional cost. We have also recommended the usage of the def2-TZVPP basis set as it represents a good compromise between accuracy and cost, limiting the finite basis set errors to less than 1 kJ/mol.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗