Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “compound screening”

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 325 records · Page 18
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Inorganic Halide Perovskitoid TlPbI 3 for Ionizing Radiation Detection

Room temperature semiconductor detector (RTSD) materials for γ-ray and X-ray radiation are in great demand for the nonproliferation of nuclear materials as well as for biomedical imaging applications. Halide perovskites have attracted great attention as emerging and promising RTSD materials. In this contribution, the material synthesis, purification, crystal growth, crystal structure, photoluminescence properties, ionizing radiation detection performance, and electronic structure of the inorganic halide perovskitoid compound TlPbI 3 are reported on. This compound crystallizes in the ABX 3 non-perovskite crystal structure with a high density of d = 6.488 g·cm –3 , has a wide bandgap of 2.25 eV, and melts congruently at a low temperature of 360 °C without phase transitions, which allows for facile growth of high quality crystals with few thermally-activated defects. High-quality TlPbI 3 single crystals of centimeter-size are grown using the vertical Bridgman method using purified raw materials. A high electrical resistivity of ≈10 12 Ω·cm is readily obtainable, and detectors made of TlPbI 3 single crystals are highly photoresponsive to Ag K α X-rays (22.4 keV), and detects 122 keV γ-rays from 57 Co radiation source. The electron mobility-lifetime product µ e τ e was estimated at 1.8 × 10 –5 cm 2 ·V –1 . A high relative static dielectric constant of 35.0 indicates strong capability in screening carrier scattering and charged defects in TlPbI 3 .

36 MATERIALS SCIENCE↗
DOE OSTI2024

Two-dimensional heavy fermions in the van der Waals metal CeSiI

Heavy-fermion metals are prototype systems for observing emergent quantum phases driven by electronic interactions. A long-standing aspiration is the dimensional reduction of these materials to exert control over their quantum phases, which remains a significant challenge because traditional intermetallic heavy-fermion compounds have three-dimensional atomic and electronic structures. Here we report comprehensive thermodynamic and spectroscopic evidence of an antiferromagnetically ordered heavy-fermion ground state in CeSiI, an intermetallic comprising two-dimensional (2D) metallic sheets held together by weak interlayer van der Waals (vdW) interactions. Owing to its vdW nature, CeSiI has a quasi-2D electronic structure, and we can control its physical dimension through exfoliation. The emergence of coherent hybridization of f and conduction electrons at low temperature is supported by the temperature evolution of angle-resolved photoemission and scanning tunnelling spectra near the Fermi level and by heat capacity measurements. Electrical transport measurements on few-layer fakes reveal heavy-fermion behaviour and magnetic order down to the ultra-thin regime. Importantly, our work establishes CeSiI and related materials as a unique platform for studying dimensionally confined heavy fermions in bulk crystals and employing 2D device fabrication techniques and vdW heterostructures to manipulate the interplay between Kondo screening, magnetic order and proximity effects.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Multicode benchmark on simulated Ti K-edge x-ray absorption spectra of Ti-O compounds

X-ray absorption spectroscopy (XAS) is an element-specific materials characterization technique that is sensitive to structural and electronic properties. First-principles simulated XAS has been widely used as a powerful tool to interpret experimental spectra and draw physical insights. Recently, there has also been growing interest in building computational XAS databases to enable data analytics and machine learning applications. However, there are non-trivial differences among commonly used XAS simulation codes, both in underlying theoretical formalism and in technical implementation. Reliable and reproducible computational XAS databases require systematic benchmark studies. In this work, we benchmarked Ti K-edge XAS simulations of ten representative Ti-O binary compounds, which we refer to as the Ti-O-10 dataset, using three state-of-the-art codes: XSPECTRA, OCEAN and exciting. We systematically studied the convergence behavior with respect to the input parameters and developed a workflow to automate and standardize the calculations to ensure converged spectra. Our benchmark comparison considers a 35 eV spectral range starting from the K-edge onset, representative of widely used near-edge spectra. Quantitative comparison over this range is based on Spearman’s rank correlation score (r sp ). Our results show: (1) the two Bethe-Salpeter equation (BSE) codes (OCEAN and exciting) have excellent agreement with an average $r$ $sp$ of 0.998; (2) good agreement is obtained between the core-hole potential code (XSPECTRA) and BSE codes (OCEAN and exciting) with an average $r$ $sp$ of 0.990, and this smaller $r$ $sp$ reflects the noticeable differences in the main edge spectral shape that can be primarily attributed to the difference in the strength of the screened core-hole potential; (3) simulations from both methods overall reproduce well the main experimental spectral features of rutile and anatase, and the different treatments of the screened core-hole potential have visible impact on pre-edge intensities and the peak ratio of the main edge; (4) there exist moderate differences in the relative edge alignment of the three codes with a standard deviation of about 0.2 eV, which arise from multiple contributions including the frozen core approximation, final state effects, and different approximations used for the self-energy correction. In conclusion, our benchmark study provides important standards for first-principles XAS simulations with broad impact in data-driven XAS analysis.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Deep Learning-Enabled MS/MS Spectrum Prediction Facilitates Automated Identification Of Novel Psychoactive Substances

The market for illicit drugs has been reshaped by the emergence of more than 1100 new psychoactive substances (NPS) over the past decade, posing a major challenge to the forensic and toxicological laboratories tasked with detecting and identifying them. Tandem mass spectrometry (MS/MS) is the primary method used to screen for NPS within seized materials or biological samples. The most contemporary workflows necessitate labor-intensive and expensive MS/MS reference standards, which may not be available for recently emerged NPS on the illicit market. Here, we present NPS-MS, a deep learning method capable of accurately predicting the MS/MS spectra of known and hypothesized NPS from their chemical structures alone. NPS-MS is trained by transfer learning from a generic MS/MS prediction model on a large data set of MS/MS spectra. We show that this approach enables a more accurate identification of NPS from experimentally acquired MS/MS spectra than any existing method. We demonstrate the application of NPS-MS to identify a novel derivative of phencyclidine (PCP) within an unknown powder seized in Denmark without the use of any reference standards. We anticipate that NPS-MS will allow forensic laboratories to identify more rapidly both known and newly emerging NPS. NPS-MS is available as a web server at https://nps-ms.ca/, which provides MS/MS spectra prediction capabilities for given NPS compounds. Additionally, it offers MS/MS spectra identification against a vast database comprising approximately 8.7 million predicted NPS compounds from DarkNPS and 24.5 million predicted ESI-QToF-MS/MS spectra for these compounds.

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

Discovery and characterization of bromodomain 2–specific inhibitors of BRDT

Bromodomain testis (BRDT), a member of the bromodomain and extraterminal (BET) subfamily that includes the cancer targets BRD2, BRD3, and BRD4, is a validated contraceptive target. All BET subfamily members have two tandem bromodomains (BD1 and BD2). Knockout mice lacking BRDT-BD1 or both bromodomains are infertile. Treatment of mice with JQ1, a BET BD1/BD2 nonselective inhibitor with the highest affinity for BRD4, disrupts spermatogenesis and reduces sperm number and motility. To assess the contribution of each BRDT bromodomain, we screened our collection of DNA-encoded chemical libraries for BRDT-BD1 and BRDT-BD2 binders. High-enrichment hits were identified and resynthesized off-DNA and examined for their ability to compete with JQ1 in BRDT and BRD4 bromodomain AlphaScreen assays. These studies identified CDD-1102 as a selective BRDT-BD2 inhibitor with low nanomolar potency and >1,000-fold selectivity over BRDT-BD1. Structure–activity relationship studies of CDD-1102 produced a series of additional BRDT-BD2/BRD4-BD2 selective inhibitors, including CDD-1302, a truncated analog of CDD-1102 with similar activity, and CDD-1349, an analog with sixfold selectivity for BRDT-BD2 versus BRD4-BD2. BROMOscan bromodomain profiling confirmed the great affinity and selectivity of CDD-1102 and CDD-1302 on all BET BD2 versus BD1 with the highest affinity for BRDT-BD2. Cocrystals of BRDT-BD2 with CDD-1102 and CDD-1302 were determined at 2.27 and 1.90 Å resolution, respectively, and revealed BRDT-BD2 specific contacts that explain the high affinity and selectivity of these compounds. These BD2-specific compounds and their binding to BRDT-BD2 are unique compared with recent reports and enable further evaluation of their nonhormonal contraceptive potential in vitro and in vivo.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

DNA-encoded chemical libraries yield non-covalent and non-peptidic SARS-CoV-2 main protease inhibitors

The development of SARS-CoV-2 main protease (M pro ) inhibitors for the treatment of COVID-19 has mostly benefitted from X-ray structures and preexisting knowledge of inhibitors; however, an efficient method to generate M pro inhibitors, which circumvents such information would be advantageous. As an alternative approach, we show here that DNA-encoded chemistry technology (DEC-Tec) can be used to discover inhibitors of M pro . An affinity selection of a 4-billion-membered DNA-encoded chemical library (DECL) using M pro as bait produces novel non-covalent and non-peptide-based small molecule inhibitors of M pro with low nanomolar K i values. Furthermore, these compounds demonstrate efficacy against mutant forms of M pro that have shown resistance to the standard-of-care drug nirmatrelvir. Overall, this work demonstrates that DEC-Tec can efficiently generate novel and potent inhibitors without preliminary chemical or structural information.

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

Triacetic Acid Lactone and 4-Hydroxycoumarin as Bioprivileged Molecules for the Development of Performance-Advantaged Organic Corrosion Inhibitors

Biomass conversion, especially the development of bioprivileged molecules utilizing integration of biological and chemical processes, has shown potential to produce novel chemicals with enhanced properties. Here, organic corrosion inhibitors based on triacetic acid lactone and 4-hydroxycoumarin have been synthesized in good yield and tested for corrosion inhibition on mild steel in both sulfuric acid and hydrochloric acid. Sixteen novel corrosion inhibitors derived from triacetic acid lactone and 4-hydroxycoumarin were efficiently synthesized and twelve of them showed high corrosion inhibition efficiency as confirmed by electrochemical impedance spectroscopy (giving values greater than 78%) and polarization analysis. While triacetic acid lactone-based compounds showed good corrosion inhibition performance, the 4-hydroxycoumarin-based compounds showed further improvement in corrosion inhibition performance. Scanning electron microscopy (SEM) analysis of the mild steel coupon in the presence of selected inhibitors showed that corrosion was significantly diminished, with complementary X-ray photoelectron spectroscopy (XPS) analysis suggesting the inhibitor molecules strongly adsorbed on the steel surface. Quantum chemical calculations showed poor correlation between calculated parameters and the performance of these molecules. This large set of corrosion inhibitors containing a dozen promising corrosion inhibitors is useful for future rational design of corrosion inhibitors. In conclusion, the results demonstrate a new direction for the development of bioprivileged molecules.

09 BIOMASS FUELS↗
DOE OSTI2023

Crystal growth of new high light yield halide perovskite scintillator RbSrI 3

X-ray radiography systems are a key national security technology enabling non-invasive cargo screening. However, the scintillators currently used in these systems have limitations in performance that prompt research for new suitable materials. Here, in this work, crystal growth, physical properties, and optical and scintillation properties of a new halide perovskite scintillator RbSrI 3 were investigated for the first time. The melting point was determined to be 453 °C. Divalent europium and ytterbium were investigated as potential luminescent activators. Light yield and energy resolution of RbSrI 3 :Eu 5 mol% was 78,700 ph/MeV and 2.8 % at 662 keV, and that of RbSrI 3 :Yb 0.5 mol% was 24,000 ph/MeV and 4.9 % at 662 keV.

36 MATERIALS SCIENCE↗
DOE OSTI2022

High-throughput platform for yeast morphological profiling predicts the targets of bioactive compounds

Abstract Morphological profiling is an omics-based approach for predicting intracellular targets of chemical compounds in which the dose-dependent morphological changes induced by the compound are systematically compared to the morphological changes in gene-deleted cells. In this study, we developed a reliable high-throughput (HT) platform for yeast morphological profiling using drug-hypersensitive strains to minimize compound use, HT microscopy to speed up data generation and analysis, and a generalized linear model to predict targets with high reliability. We first conducted a proof-of-concept study using six compounds with known targets: bortezomib, hydroxyurea, methyl methanesulfonate, benomyl, tunicamycin, and echinocandin B. Then we applied our platform to predict the mechanism of action of a novel diferulate-derived compound, poacidiene. Morphological profiling of poacidiene implied that it affects the DNA damage response, which genetic analysis confirmed. Furthermore, we found that poacidiene inhibits the growth of phytopathogenic fungi, implying applications as an effective antifungal agent. Thus, our platform is a new whole-cell target prediction tool for drug discovery.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

High-throughput characterization of Lu-doped zirconia

Lu-doped zirconia (LuDZ) was reported in the 1960s to be a poor oxygen ion conductor. However, we identified this system as worthy of reexamination based on the results of a machine learning-based screen for promising oxygen ion conductors, further supported by a physical argument: because Lu is the heaviest and smallest of the lanthanides, it should offer the possibility of high conductivity due to its high polarizability and size compatibility with the zirconia structure. Here the relationship between dopant concentration and conductivity in the LuDZ system was determined using a high throughput experimental methodology in which spatially resolved conductivity was measured across a compositionally graded film, spanning from Zr 0.93 Lu 0.07 O 1.965 to Zr 0.82 Lu 0.18 O 1.91 . The conductivity broadly peaks with composition at 15 cation % Lu, coinciding with the approximate phase boundary of the fully stabilized cubic phase. The composition-dependent conductivity measurements reveal a surprising increase in activation energy for ion transport with increasing Lu content with no change in this trend at the composition of the conductivity maximum. The pre-exponential factor increases sharply with Lu content, but displays a distinct plateau once the conductivity maximum composition is attained. At 500 °C, the conductivity of Zr 0.85 Lu 0.15 O 1.925 is 2 × 10 –3 S cm –1 , approaching that of optimally doped Sc-stabilized zirconia, which has the highest reported conductivity among zirconia compounds. In comparison to other oxide ion conductors, Lu-stabilized zirconia offers the advantage over ceria of negligible electronic conductivity and is easier to process than (La,Sr)(Ga,Mg)O 3-δ .

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

Evaluation of GlassNet for physics-informed machine learning of glass stability and glass-forming ability

Glassy materials form the basis of many modern applications, including nuclear waste immobilization, touch-screen displays, and optical fibers, and also hold great potential for future medical and environmental applications. However, their structural complexity and large composition space make design and optimization challenging for certain applications. Of particular importance for glass processing and design is an estimate of a given composition's glass-forming ability (GFA). However, there remain many open questions regarding the underlying physical mechanisms of glass formation, especially in oxide glasses. It is apparent that a proxy for GFA would be highly useful in glass processing and design, but identifying such a surrogate property has proven itself to be difficult. While glass stability (GS) parameters have historically been used as a GFA surrogate, recent research has demonstrated that most of these parameters are not accurate predictors of the GFA of oxide glasses. Here, in this work, we explore the application of an open-source pre-trained neural network model, GlassNet, that can predict the characteristic temperatures necessary to compute GS with reasonable performance and assess the feasibility of using these physics-informed machine learning (PIML)-predicted GS parameters to estimate GFA. In doing so, we track the uncertainties at each step of the computation—from the original ML prediction errors to the compounding of errors during GS estimation, and finally to the final estimation of GFA. While GlassNet exhibits reasonable accuracy on all individual properties, we observe a large compounding of error in the combination of these individual predictions for the PIML prediction of GS, finding that random forest models offer similar accuracy to GlassNet. We also break down the performance of GlassNet on different glass families and find that the error in GS prediction is correlated with the error in crystallization peak temperature prediction. Lastly, we utilize this finding to assess the relationship between top-performing GS parameters and GFA for two ternary glass systems: sodium borosilicate and sodium iron phosphate glasses. We conclude that to obtain true ML predictive capability of GFA, significantly more data needs to be collected.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Non-Contact Vapor Detection of Illicit Drugs via Atmospheric Flow Tube-Mass Spectrometry

Real-time, non-contact detection of illicit drugs is a desirable goal for the interdiction of these controlled substances, but the relatively low vapor pressures of such species present a challenge for trace vapor detection technologies. The introduction of atmospheric flow tube-mass spectrometry (AFT-MS), which has previously been demonstrated to detect gas-phase analytes at low parts-per-quadrillion levels for explosives and organophosphorus compounds, also enables the potential for non-contact drug detection. With AFT-MS, direct vapor detection of cocaine and methamphetamine from ~5 µg residues at room temperature is demonstrated herein. Furthermore, thermal desorption of low- to sub-picogram levels of cocaine, methamphetamine, fentanyl, and heroin is observed via AFT-MS using a carrier flow rate of several L/min of air. These low levels can permit non-contact sampling through collection of vapor, effectively preconcentrating the analyte before desorption and analysis. Quantitative evaluation of the thermal desorption approach has yielded limits of detection (LODs) on order of 10 fg for cocaine and fentanyl, 100 fg for methamphetamine, and 1.6 pg for heroin. The LOD for heroin was lowered to 200 fg by using tributyl phosphate as a dopant to form a proton-bound heterodimer with heroin. When used with AFT-MS, the intentional formation of specific drug-dopant adducts has the potential to enhance detection limits and selectivity of additional drug species. Species that are prone to form adducts present a challenge to analysis, but that difficulty can be overcome by the intentional addition of a dopant. Molecules unlikely to form adducts will remain essentially unimpacted, but the adduct-forming species will interact with the dopant to compress the analyte signal into a single peak. This approach would be valuable in the application of non-contact screening for illicit substances via vapor collection followed by thermal desorption for analysis.

methamphetamine, cocaine, heroin, fentanyl, atmosp↗
DOE OSTI2025

Technical Background and Validation Report on the Residential Water Inhalation Risk Calculator Presented in the Risk Assessment Information System

Indoor air quality (IAQ) is critical for human health. Poor IAQ is linked to respiratory issues, cardiovascular diseases, and cancer. Indoor pollutants are emitted by typical household items such as cleaning products, personal care items, building materials, and tap water - an understudied volatile organic compound (VOC) source. This document presents the Residential Water Inhalation Risk Calculator (RWIRC), which estimates daily VOC exposure concentrations from various household water uses, such as showering and dishwashing, to assess exposure risks for the most vulnerable occupant. Integrated into the Risk Assessment Information System (RAIS) and sponsored by the US Department of Energy (DOE), the calculator divides a house into three compartments: shower, bathroom, and other spaces, accounting for daily water usage patterns and calculating VOC concentrations. Exposure data generated using the calculator can assist health assessors in estimating excess lifetime cancer risk (ELCR) and hazard index (HI) from VOC inhalation. Unlike traditional exposure models that utilize Andelman’s constant, the RWIRC continuously assesses variability in VOC concentrations and environmental conditions using differential equations to track VOC concentrations and air exchange between compartments. The calculator also provides unique volatilization fractions for each chemical and appliance, enhancing accuracy of the exposure concentration estimation. The RWIRC is accessible online and allows users to customize parameters (i.e., number of bathrooms, water temperature, and exhaust fan conditions) and input VOC characteristics (i.e., tap water and ambient air concentrations). This document provides a step-by-step guide on implementing the calculator. It also provides comparisons with the ATSDR-SHOWER calculator, using eight VOCs with varying physicochemical properties to reveal differences in algorithms and output concentrations. Simulations also assess how bathroom door positions and exhaust fan usage affect VOC exposure. The calculator results can enhance EPA risk screening levels for inhalation exposure to VOCs from tap water, offering a sophisticated tool for assessing inhalation risks and improving public health protection.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

High-dimensional isotomics, part 2: Observations of over 100 constraints on methionine's isotome

The abundances of different isotopic forms of a compound, or isotopologues, will vary based on its physical and chemical history. The number of isotopologues increases combinatorically with the size of a molecule, and even small molecules such as amino acids have thousands of potentially observable isotopic variants. However, due to the analytical challenges of separating and observing isotopologues, only a few dimensions of isotopic diversity are routinely measured. Overcoming these challenges requires both an experimental method to observe many isotopic properties and a theoretical framework for interpreting these experiments. In Part 1, we presented such a theoretical framework; here, we demonstrate an experimental method, which we apply to methionine. Our approach uses a Q Exactive HF Orbitrap to perform several “M + N experiments”, where a sample is ionized, a subset of its isotopologues with cardinal mass N daltons greater than the unsubstituted isotopologue is selected and fragmented, and the proportions of all detectable isotopic forms of those fragment ions are quantified. We perform M + 1, M + 2, M + 3, and M + 4 experiments of a methionine sample and standard where the sample has a 100 ‰ enrichment of 13C at the methyl carbon relative to the natural 13C abundance at that position in the standard, and is otherwise identical to the standard. We observe isotopic forms of 8 fragment ion species for each version of the M + N experiment. With the assistance of a forward model of expected mass spectra, we identify isotopic peaks for each fragment ion based on observed mass and abundance, screen these for data quality, and quantify abundances for 146 unique isotopic peaks at precisions of ≈ 0.3–3 ‰. We present our direct observations and use them to reconstruct the concentrations of 19 individual singly, doubly, and triply-substituted isotopologues; doing so gives fewer constraints and broader error bars than working with the direct observations, but may be more interpretable for some applications. We also examine possibilities for measuring additional peaks, which are primarily limited by the detection limit of the Orbitrap-IRMS method. We then suggest some possible uses of our direct measurements for chemical forensics and hypothesis testing. Furthermore, our results demonstrate the diversity of isotopic constraints currently observable and interpretable for organic molecules.

58 GEOSCIENCES↗
DOE OSTI2023

HD-Bind: Encoding of Molecular Structure with Low Precision, Hyperdimensional Binary Representations

Publicly available collections of drug-like molecules have grown to comprise tens of billions of compounds due to advances in combinatorial chemistry. Traditional methods for identifying "hit" molecules from a large collection of potential drug-like candidates have relied on biophysical theory to compute approximations to the Gibbs free energy of the binding interaction between the drug and its protein target. These approaches have the major drawback that they require exceptional computing capabilities for even relatively small collections of molecules. Hyperdimensional Computing (HDC) is a recently-proposed learning paradigm that represents data with high-dimension binary vectors; this allows the use of low-precision binary vector arithmetic to create models of the data that can be learned without the need for the gradient-based optimization required in many conventional machine learning and deep learning methods. This algorithmic simplicity allows for acceleration in hardware that has been previously demonstrated in a range of application areas. We consider existing HDC approaches for molecular property classification and introduce two novel encodings of a commonly-used molecular representation, the extended connectivity fingerprint (ECFP). We show that HDC-based inference methods are as much as 91 times more efficient than traditional machine learning methods, and achieve an acceleration of nearly nine orders of magnitude compared to molecular docking. Our results show that HDC accelerated methods retain competitive accuracy on a number of well-studied tasks such as molecular property predictions using the MoleculeNet dataset, and bind/no-bind activity classification using the DUD-E and LIT-PCBA datasets. Our work thus motivates further investigation into molecular representation learning to develop ultraefficient pre-screening tools.

Jones, William↗
DOE OSTI2020

Supercomputer-Based Ensemble Docking Drug Discovery Pipeline with Application to Covid-19

In this work, we present a supercomputer-driven pipeline for in silico drug discovery using enhanced sampling molecular dynamics (MD) and ensemble docking. Ensemble docking makes use of MD results by docking compound databases into representative protein binding-site conformations, thus taking into account the dynamic properties of the binding sites. We also describe preliminary results obtained for 24 systems involving eight proteins of the proteome of SARS-CoV-2. The MD involves temperature replica exchange enhanced sampling, making use of massively parallel supercomputing to quickly sample the configurational space of protein drug targets. Using the Summit supercomputer at the Oak Ridge National Laboratory, more than 1 ms of enhanced sampling MD can be generated per day. We have ensemble docked repurposing databases to 10 configurations of each of the 24 SARS-CoV-2 systems using AutoDock Vina. Comparison to experiment demonstrates remarkably high hit rates for the top scoring tranches of compounds identified by our ensemble approach. We also demonstrate that, using Autodock-GPU on Summit, it is possible to perform exhaustive docking of one billion compounds in under 24 h. Finally, we discuss preliminary results and planned improvements to the pipeline, including the use of quantum mechanical (QM), machine learning, and artificial intelligence (AI) methods to cluster MD trajectories and rescore docking poses.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2021

Identification of Known and Novel Monomers for Poly(hydroxyurethanes) from Biobased Materials

Sourcing polymers from biobased materials is desirable for a transition to a more sustainable world. However, finding new monomers sourced from biobased molecules requires extensive experimental effort, due to the fact that most biomolecules lack the functional groups or properties needed to serve as a monomer in a particular polymeric system without chemical modification. In this work, we demonstrate a computational screening method to obtain monomers of a non-isocyanate polyurethane (NIPU), polyhydroxyurethane (PHU), from biomolecules. Our method generates candidate monomers from a list of 15 selected biobased starting molecules by performing chemical transformations that are standard to chemical catalysis and organic synthesis on these starting molecules and then screening the results for molecules with moieties that can be used as monomers for a PHU. Product molecules from two generations of reaction, a total number of 20,057, were analyzed automatically to identify molecules that have two or more epoxy or cyclic carbonate functional groups as candidates for further derivatization. These 21 candidate molecules were then subjected to another two generations of reaction, producing a set of 17,913 additional molecules. Further scrutiny of this set for desired functional groups characterizing PHUs and querying against the PubChem database was carried out, and 15 known molecules that are commercially available were identified, with the desirable features that they are predicted to be derived from biobased starting molecules in four or fewer steps. Moreover, molecule space generated from this process was primarily comprised of molecules that were not known to PubChem, indicating the framework also has the ability to generate novel candidate monomers that originate from biobased starting points.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Machine Learning-Driven Solvent Screening for Biobased 2,3-Butanediol Extraction

Biobased 2,3-butanediol (2,3-BDO) is a valuable biomass-derived chemical due to its versatility in being transformed into a wide variety of products. However, the separation and purification of 2,3-BDO from fermentation broth remain a significant challenge owing to its high boiling point and hydrophilic nature. Herein, we developed a machine learning (ML)-based screening workflow that uses molecular calculations as training data and requires only a small number of experimental measurements for validation to identify alternative solvent candidates for the liquid–liquid extraction (LLE) of 2,3-BDO from aqueous solution. In particular, 130 density functional theory (DFT) calculations with the implicit solvation method not only built a correlation between the computational partition coefficient and the experimental distribution coefficient of 2,3-BDO but also parameterized an Extra-Trees ML model to screen the distribution coefficient for a wider range of 6717 organic solvents. The experimental measurements of only 24 solvents were needed to validate the computational results. A list of 50 prioritized solvents was proposed for 2,3-BDO LLE, and seven additional experimental measurements were conducted to further verify our selected solvents. The impact of the extraction temperature and solvent-to-feed ratio was also investigated for selected solvents in experiments. Furthermore, this work suggested alternative solvents for 2,3-BDO LLE and proposed a versatile workflow that requires fewer experiments and can be applied to a broader range of LLE studies.

Extraction↗