Engineering topics

Smith, Micholas Dean

Publications and source records attributed to Smith, Micholas Dean.

15 records

DOE OSTI2024

Molecular insights into the in situ early-stage assembly of metal–organic frameworks on cellulose nanofibrils

The integration of metal–organic frameworks (MOFs) into renewable matrices, including those derived from cellulose, to create cellulose/MOF hybrids has attracted significant interest due to the synergistic combination of cellulosic biopolymer properties and MOFs' multifunctional features. The interfacial interactions between cellulose and MOF not only ensure stable bonding, but ultimately also determine the bulk physical properties of the macroscopic composites. However, the mechanistic understanding of the in situ assembly between the two materials remains unclear. In this study, we employ a combination study of synthesis, experimental characterization, and molecular dynamics simulations to explore the early-stage assembly of a cellulose/MOF hybrid. Our results revealed that the growth of MOF clusters on the 2,2,6,6-tetramethyl-1-piperidinyloxy oxidized cellulose nanofibrils (TOCNF) follows an inhomogeneous sequential transformation pathway. The carboxylates of TOCNF form coordination-like bonds with the metal ions, while the hydroxyl groups of TOCNF form hydrogen bonds with MOF ligands. These interactions provide the initial nucleation sites that mediate the growth of MOF clusters and also guide the assembly of larger MOF clusters onto the TOCNF substrate. In conclusion, the fundamental insights into the in situ assembly of MOF nanoparticles on cellulosic substrates are essential for the rational design of high-performance materials with tailored morphology and optimized properties.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2024

Macroscale properties and atomic-scale mechanisms of ash removal in low-temperature hydrothermal carbonization

Biogenic ash is a significant impediment to the utilization of agricultural residues in biofuel production. Such challenge can be addressed by various treatments, as demonstrated in this study on the experimental and computational mechanisms involved in the hydrothermal treatment (HT) of wheat straw. A combination of classical (all-atom) molecular dynamics simulations of cellulose carrying silica and calcium species, along with first principles quantum chemical calculations, indicates the dissociation of inorganics from the cellulose with increased HT temperature. Here, this observation is confirmed by experimental evidence of effective ash removal by HT, showing at least 50% removal of sulfur, chlorine, potassium, and calcium, and 12.5% of silica, leading to a reduced total ash content (from 6.7% to 4.2%).Changes in structural features upon HT, such as surface cellular structure and porosity, were revealed, accompanied by an increased specific surface area (from 1.17 to 6.34 m 2 /g). Our simulations suggest that silica binds tightly to the hydrophobic face of cellulose at room temperature, but HT significantly reduces the binding free energy of association with both hydrophobic and hydrophilic surfaces. Most significantly, ash removal leads to an increased calorific value, rising from approximately 16 MJ/kg to about 19 MJ/kg, along with improved thermal behavior. The improved integration combustion index parameter S indicates that the combustion properties improve with ash removal efficiency. The proposed atomic-level mechanism for the observed removal of inorganics during mild HT underscores the potential of such treatment in producing energy-dense wheat straw, a widely available agricultural residue.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2024

Molecular-level design of alternative media for energy-saving pilot-scale fibrillation of nanocellulose

The outstanding mechanical properties, light weight, and biodegradability of cellulose nanofibrils (CNFs) make them promising components of renewable and sustainable next-generation reinforced composite biomaterials and bioplastics. Manufacturing CNFs at a pilot scale requires disc-refining fibrillation of dilute cellulose fibers in aqueous pulp suspensions to shear the fibers apart into their nanodimensional forms, which is, however, an energy-intensive process. Here, in this study, we used atomistic molecular dynamics (MD) simulation to examine media that might facilitate the reduction of interactions between cellulose fibers, thereby reducing energy consumption in fibrillation. The most suitable medium found by the simulations was an aqueous solution with 0.007:0.012 wt.% NaOH:urea, and indeed this was found in pilot-scale experiments to reduce the fibrillation energy by ~21% on average relative to water alone. The NaOH:urea-mediated CNFs have similar crystallinity, morphology, and mechanical strength to those formed in water. The NaOH and urea act synergistically on CNFs to aid fibrillation but at different length scales. NaOH deprotonates hydroxyl groups leading to mesoscale electrostatic repulsion between fibrils, whereas urea forms hydrogen bonds with protonated hydroxyl groups thus disrupting interfibril hydrogen bonds. This suggests a general mechanism in which an aqueous medium that contains a strong base and a small organic molecule acting as a hydrogen-bond acceptor and/or donor may be effectively employed in materials processes where dispersion of deprotonable polymers is required. The study demonstrates how atomic-detail computer simulation can be integrated with pilot-scale experiments in the rational design of materials processes for the circular bioeconomy.

36 MATERIALS SCIENCE↗

DOE OSTI2024

Osteocalcin binds to a GPRC6A Venus fly trap allosteric site to positively modulate GPRC6A signaling

GPRC6A, a member of the Family C G-protein coupled receptors, regulates energy metabolism and sex hormone production and is activated by diverse ligands, including cations, L-amino acids, the osteocalcin (Ocn) peptide and the steroid hormone testosterone. We sought a structural framework for the ability of multiple distinct classes of ligands to active GPRC6A. We created a structural model of GPRC6A using Alphafold2. Using this model we explored a putative orthosteric ligand binding site in the bilobed Venus fly trap (VFT) domain of GPRC6A and two positive allosteric modulator (PAM) sites, one in the VFT and the other in the 7 transmembrane (7TM) domain. We provide evidence that Ocn peptides act as a PAM for GPRC6A by binding to a site in the VFT that is distinct from the orthosteric site for calcium and L-amino acids. In agreement with this prediction, alternatively spliced GPRC6A isoforms 2 and 3, which lack regions of the VFT, and mutations in the computationally predicted Ocn binding site, K352E and H355P, prevent Ocn activation of GPRC6A. These observations explain how dissimilar ligands activate GPRC6A and set the stage to develop novel molecules to activate and inhibit this previously poorly understood receptor.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2024

Amantadine interactions with phase separated lipid membranes

Amantadine, a small amphilphic organic compound that consists of an adamantane backbone and an amino group, was first recognized as an antiviral in 1963 and received approval for prophylaxis against the type A influenza virus in 1976. Since then, it has also been used to treat Parkinson’s disease-related dyskinesia and is being considered as a treatment for corona viruses. Since amantadine usually targets membrane-bound proteins, its interactions with the membrane are also thought to be important. Biological membranes are now widely understood to be laterally heterogeneous and certain proteins are known to preferentially co-localize within specific lipid domains. Does amantadine, therefore, preferentially localize in certain lipid composition domains? To address this question, here, we studied amantadine’s interactions with phase separating membranes composed of cholesterol, DSPC (1,2-distearoyl-sn-glycero-3-phosphocholine), POPC (1-palmitoyl-2-oleoyl-glycero-3-phosphocholine), and DOPC (1,2-dioleoyl-sn-glycero-3-phosphocholine), as well as single-phase DPhPC (1,2-diphytanoyl-sn-glycero-3-phos-phocholine) membranes. From Langmuir trough and differential scanning calorimetry (DSC) measurements, we determined, respectively, that amantadine preferentially binds to disordered lipids, such as POPC, and lowers the phase transition temperature of POPC/DSPC/cholesterol mixtures, implying that amantadine increases membrane disorder. Further, using droplet interface bilayers (DIBs), we observed that amantadine disrupts DPhPC membranes, consistent with its disordering properties. Finally, we carried out molecular dynamics (MD) simulations on POPC/DSPC/cholesterol membranes with varying amounts of amantadine. Consistent with experiment, MD simulations showed that amantadine prefers to associate with disordered POPC-rich domains, domain boundaries, and lipid glycerol backbones. Since different proteins co-localize with different lipid domains, our results have possible implications as to which classes of proteins may be better targets for amantadine.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2024

High-Throughput Screening and Accurate Prediction of Ionic Liquid Viscosities Using Interpretable Machine Learning

Ionic liquids (ILs) are a novel group of green solvents with great promise for various industrial applications, including carbon capture and lignocellulosic biomass deconstruction. However, the use of ILs at the industrial scale remains challenging due to their high viscosities at ambient temperatures. To develop ILs with lower viscosities, a systematic study of their quantitative structure–property relationship (QSPR) is desirable. Here, we developed four machine learning (ML) models to predict viscosity at various temperature and pressure ranges, trained over a wide range of ILs consisting of various cationic and anionic families. ML methods including two-factor polynomial regression (two-factor PR), support vector regression (SVR), feed-forward neural networks (FFNN), and categorical boosting (CATBoost) were developed based on features that have proven useful in previous ML studies: COSMO-RS (conductor-like screening model for real solvents)-derived surface screening charge densities (sigma profiles). FFNN and CATBoost were the most accurate in predicting IL viscosities with lower average absolute relative deviation and higher R2 values on the test set. Tanimoto similarity scores were calculated to characterize the chemical space and structural similarity of the investigated ions. Furthermore, SHapley Additive exPlanation (SHAP) analysis was employed to interpret the ML results. Temperature, the polar area of ILs, and the nonpolar regions of ions are key features that influence the viscosity predictions. Importantly, the IL viscosity prediction here is the most accurate reported to date.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2024

Accurate Machine Learning for Predicting the Viscosities of Deep Eutectic Solvents

Deep eutectic solvents (DESs) are emerging as environmentally friendly designer solvents for mass transport and heat transfer processes in industrial applications; however, the lack of accurate tools to predict and thus control their viscosities under both a range of environmental factors and formulations hinders their general application. While DESs may serve as designer solvents, with nearly unlimited combinations, this unfortunately makes it experimentally infeasible to comprehensively measure the viscosities of all DESs of potential industrial interest. To assist in the design of DESs, we have developed several new machine learning (ML) models that accurately and rapidly predict the viscosities of a diverse group of DESs at different temperatures and molar ratios using, to date, one of the most comprehensive data sets containing the properties of over 670 DESs over a wide range of temperatures (278.15–385.25 K). Three ML models, including support vector regression (SVR), feed forward neural networks (FFNNs), and categorical boosting (CatBoost), were developed to predict DES viscosity as a function of temperature and molar ratio and contrasted with multilinear and two-factor polynomial regression baselines. Further, quantum chemistry-based, COSMO-RS-derived sigma profile (σ-profile) features were used as inputs for the ML models. The CatBoost model is excellent at externally predicting DES viscosity, as indicated by high R 2 (0.99) and low root-mean-square-error (RMSE) and average absolute relative deviations (AARD) (5.22%) values for the testing data sets, and 98% of the data points lie within the 15% of AARD deviations. Furthermore, SHapley additive explanation (SHAP) analysis was employed to interpret the ML results and rationalize the viscosity predictions. The result is an ML approach that accurately predicts viscosity and will aid in accelerating the design of appropriate DESs for industrial applications.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2024

Drugging the entire human proteome: Are we there yet?

Each of the ~20 000 proteins in the human proteome is a potential target for compounds that bind to it and modify its function. The 3D structures of most of these proteins are now available. Here, we discuss the prospects for using these structures to perform proteome-wide virtual HTS (VHTS). Furthermore, we compare physics-based (docking) and AI VHTS approaches, some of which are now being applied with large databases of compounds to thousands of targets. Although preliminary proteome-wide screens are now within our grasp, further methodological developments are expected to improve the accuracy of the results.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2023

Genetic interactions between polycystin-1 and Wwtr1 in osteoblasts define a novel mechanosensing mechanism regulating bone formation in mice

Molecular mechanisms transducing physical forces in the bone microenvironment to regulate bone mass are poorly understood. Here, we used mouse genetics, mechanical loading, and pharmacological approaches to test the possibility that polycystin-1 and Wwtr1 have interdependent mechanosensing functions in osteoblasts. We created and compared the skeletal phenotypes of control Pkd1 flox/+ ;Wwtr1 flox/+ , Pkd1 Oc-cKO , Wwtr1 Oc-cKO , and Pkd1/Wwtr1 Oc-cKO mice to investigate genetic interactions. Consistent with an interaction between polycystins and Wwtr1 in bone in vivo, Pkd1/Wwtr1 Oc-cKO mice exhibited greater reductions of BMD and periosteal MAR than either Wwtr1 Oc-cKO or Pkd1 Oc-cKO mice. Micro-CT 3D image analysis indicated that the reduction in bone mass was due to greater loss in both trabecular bone volume and cortical bone thickness in Pkd1/Wwtr1 Oc-cKO mice compared to either Pkd1 Oc-cKO or Wwtr1 Oc-cKO mice. Pkd1/Wwtr1 Oc-cKO mice also displayed additive reductions in mechanosensing and osteogenic gene expression profiles in bone compared to Pkd1 Oc-cKO or Wwtr1 Oc-cKO mice. Moreover, we found that Pkd1/Wwtr1 Oc-cKO mice exhibited impaired responses to tibia mechanical loading in vivo and attenuation of load-induced mechanosensing gene expression compared to control mice. Finally, control mice treated with a small molecule mechanomimetic, MS2 that activates the polycystin complex resulted in marked increases in femoral BMD and periosteal MAR compared to vehicle control. In contrast, Pkd1/Wwtr1 Oc-cKO mice were resistant to the anabolic effects of MS2. These findings suggest that PC1 and Wwtr1 form an anabolic mechanotransduction signaling complex that mediates mechanical loading responses and serves as a potential novel therapeutic target for treating osteoporosis.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2023

Characterization and molecular simulation of lignin in Cyrene pretreatment of switchgrass

Biomass-derived solvents have been proposed as a novel pathway in biorefining for the realization of biofuels and bioproducts derived from lignocellulosic biomass. Cyrene derived from cellulose has recently been shown to have a high potential as a green organic solvent for pretreating poplar biomass. However, due to its high dynamic viscosity nature, high Cyrene concentration could cause negative effects on the sugar release of the pretreated biomass as well as driving up the operational cost of the lignin recovery. In this study, we combine experimental and computational approaches to examine the impact of Cyrene pretreatment with reduced Cyrene concentration under mild conditions on switchgrass lignin. Our experimental studies indicated correlation between pretreatment condition and recovery and structure modification of lignin. Switchgrass lignin extracted by Cyrene pretreatment possessed high preservation of β-O-4 ether inter-unit linkage, which could provide versatility in the integration of downstream lignin valorization into the modern biorefinery industries. Molecular modeling examining the solvation of switchgrass lignin polymer and the disaggregation of low-molecular weight lignin aggregates under pretreatment conditions indicated that a preferential interaction exists between Cyrene and lignin, which likely drives lignin release, and that the disruption of inter-lignin contacts can be modulated as a non-monotonic function of Cyrene : water ratio. Further, while Cyrene–lignin interactions permit the solubilization of lignin, simulations with proxy reactive-species reveal that changes to the diffusion of these reactive proxies and their localization near linkage sites under Cyrene conditions may inhibit chemical processes. In conclusion, the results indicated that loss of pretreatment efficacy caused by low Cyrene concentration could be compensated by prolonged pretreatment time and high catalyst dosage.

09 BIOMASS FUELS↗

DOE OSTI2023

Predictive understanding of the surface tension and velocity of sound in ionic liquids using machine learning

Knowledge of the physical properties of ionic liquids (ILs), such as the surface tension and speed of sound, is important for both industrial and research applications. Unfortunately, technical challenges and costs limit exhaustive experimental screening efforts of ILs for these critical properties. Previous work has demonstrated that the use of quantum-mechanics-based thermochemical property prediction tools, such as the conductor-like screening model for real solvents, when combined with machine learning (ML) approaches, may provide an alternative pathway to guide the rapid screening and design of ILs for desired physiochemical properties. However, the question of which machine-learning approaches are most appropriate remains. In the present study, we examine how different ML architectures, ranging from tree-based approaches to feed-forward artificial neural networks, perform in generating nonlinear multivariate quantitative structure–property relationship models for the prediction of the temperature- and pressure-dependent surface tension of and speed of sound in ILs over a wide range of surface tensions (16.9–76.2 mN/m) and speeds of sound (1009.7–1992 m/s). The ML models are further interrogated using the powerful interpretation method, shapley additive explanations. We find that several different ML models provide high accuracy, according to traditional statistical metrics. The decision tree-based approaches appear to be the most accurate and precise, with extreme gradient-boosting trees and gradient-boosting trees being the best performers. However, our results also indicate that the promise of using machine-learning to gain deep insights into the underlying physics driving structure–property relationships in ILs may still be somewhat premature.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2023

Quantum Chemistry-Driven Machine Learning Approach for the Prediction of the Surface Tension and Speed of Sound in Ionic Liquids

Ionic liquids (ILs) have unique solvent properties and have thus garnered significant interest. However, exhaustive experimental determination of the physicochemical properties of ILs is unrealistic due to the large structural diversity of anions and cations, their high cost, the requirements of elevated temperature and pressure, and the time required. To circumvent these experimental costs, computational approaches to accurately calculate these properties have emerged. Here in the present study, we present a demonstration of two machine learning (ML) models for the prediction of two critical IL physical properties, the surface tension and the speed of sound, across a wide range of temperatures and pressures. The models make use of molecular descriptors derived from the COSMO-RS, a quantum chemical-based model. The ML models show excellent agreement with experimental observations, with an R2 value of 0.96–0.99 and RMSE of 1.71 mN/m and 16.12 m/s for the surface tension and speed of sound, respectively. This work paves the way for the development of COSMO-RS-informed ML models for the prediction of IL properties which can help to further optimize and accelerate technology development for ILs.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2023

SARS-CoV2 billion-compound docking

Abstract This dataset contains ligand conformations and docking scores for 1.4 billion molecules docked against 6 structural targets from SARS-CoV2, representing 5 unique proteins: MPro, NSP15, PLPro, RDRP, and the Spike protein. Docking was carried out using the AutoDock-GPU platform on the Summit supercomputer and Google Cloud. The docking procedure employed the Solis Wets search method to generate 20 independent ligand binding poses per compound. Each compound geometry was scored using the AutoDock free energy estimate, and rescored using RFScore v3 and DUD-E machine-learned rescoring models. Input protein structures are included, suitable for use by AutoDock-GPU and other docking programs. As the result of an exceptionally large docking campaign, this dataset represents a valuable resource for discovering trends across small molecule and protein binding sites, training AI models, and comparing to inhibitor compounds targeting SARS-CoV-2. The work also gives an example of how to organize and process data from ultra-large docking screens.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2023

Dissecting Phenotype from Genotype with Clinical Isolates of SARS-CoV-2 First Wave Variants

The emergence and availability of closely related clinical isolates of SARS-CoV-2 offers a unique opportunity to identify novel nonsynonymous mutations that may impact phenotype. Global sequencing efforts show that SARS-CoV-2 variants have emerged and then been replaced since the beginning of the pandemic, yet we have limited information regarding the breadth of variant-specific host responses. Using primary cell cultures and the K18-hACE2 mouse, we investigated the replication, innate immune response, and pathology of closely related, clinical variants circulating during the first wave of the pandemic. Mathematical modeling of the lung viral replication of four clinical isolates showed a dichotomy between two B.1. isolates with significantly faster and slower infected cell clearance rates, respectively. While isolates induced several common immune host responses to infection, one B.1 isolate was unique in the promotion of eosinophil-associated proteins IL-5 and CCL11. Moreover, its mortality rate was significantly slower. Lung microscopic histopathology suggested further phenotypic divergence among the five isolates showing three distinct sets of phenotypes: (i) consolidation, alveolar hemorrhage, and inflammation, (ii) interstitial inflammation/septal thickening and peribronchiolar/perivascular lymphoid cells, and (iii) consolidation, alveolar involvement, and endothelial hypertrophy/margination. Together these findings show divergence in the phenotypic outcomes of these clinical isolates and reveal the potential importance of nonsynonymous mutations in nsp2 and ORF8.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2023

Amphiphilic Co-Solvents Modulate the Structure of Membrane Domains

Biofuels are an increasing part of the sustainable energy picture. This makes it a societal and economic imperative to optimize biofuel production. Mitigating the toxic effects of amphiphilic co-solvents is one way to improve the efficiency of biofuel production. Amphiphiles partition into cellular membranes, leading to membrane thinning, destabilization, loss of membrane potential, and, ultimately, cell death. However, this picture of solvent toxicity misses the disruptive impact of co-solvents on lateral membrane organization, which is increasingly recognized as critical for membrane protein sorting and oligomerization. The alteration or disruption of membrane domains has deleterious effects on cellular processes. Here, in this work, we pursue the hypothesis that membrane lateral organization is disrupted by the presence of co-solvents at concentrations lower than those which lead to full membrane destabilization. The disruption occurs due to an increasing interfacial tension between the co-existing phases, resulting in conformational changes to minimize the interfacial length-to-area ratio. This represents an unrecognized mode of solvent-induced stress and a new target for interventions to improve fermentation yields.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗