Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “in silico design”

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 91 records · Page 5

Demonstration of Decentralized Physics-Driven Learning

In typical artificial neural networks, neurons adjust according to global calculations of a central processor, but in the brain, neurons and synapses self-adjust based on local information. Contrastive learning algorithms have recently been proposed to train physical systems, such as fluidic, mechanical, or electrical networks, to perform machine-learning tasks from local evolution rules. However, to date, such systems have only been implemented in silico due to the engineering challenge of creating elements that autonomously evolve based on their own response to two sets of global boundary conditions. Here, we introduce and implement a physics-driven contrastive learning scheme for a network of variable resistors, using circuitry to locally compare the response of two identical networks subjected to the two different sets of boundary conditions. Using this method, our system effectively trains itself, optimizing its resistance values without the use of a central processor or external information storage. Once the system is trained for a specified allostery, regression, or classification task, the task is subsequently performed rapidly and automatically by the physical imperative to minimize power dissipation in response to the given voltage inputs. Here, we demonstrate that, unlike typical computers, such learning systems are robust to extreme damage (and thus manufacturing defects) due to their decentralized learning. Our twin-network approach is therefore readily scalable to extremely large or nonlinear networks, where its distributed nature will be an enormous advantage; a laboratory network of only 500 edges will already outpace its in silico counterpart.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Identifying Bayesian optimal experiments for uncertain biochemical pathway models

Abstract Pharmacodynamic (PD) models are mathematical models of cellular reaction networks that include drug mechanisms of action. These models are useful for studying predictive therapeutic outcomes of novel drug therapies in silico. However, PD models are known to possess significant uncertainty with respect to constituent parameter data, leading to uncertainty in the model predictions. Furthermore, experimental data to calibrate these models is often limited or unavailable for novel pathways. In this study, we present a Bayesian optimal experimental design approach for improving PD model prediction accuracy. We then apply our method using simulated experimental data to account for uncertainty in hypothetical laboratory measurements. This leads to a probabilistic prediction of drug performance and a quantitative measure of which prospective laboratory experiment will optimally reduce prediction uncertainty in the PD model. The methods proposed here provide a way forward for uncertainty quantification and guided experimental design for models of novel biological pathways.

97 MATHEMATICS AND COMPUTING↗

Active Learning Guided Computational Discovery of Plant-Based Redoxmers for Organic Nonaqueous Redox Flow Batteries

Organic nonaqueous redox flow batteries (ONRFBs) are promising energy storage devices due to their scalability and reliance on sourceable materials. However, finding suitable redox-active organic molecules (redoxmers) for these batteries remains a challenge. Using plant-based compounds as precursors for these redoxmers can decrease their costs and environmental toxicity. In this computational study, flavonoid molecules have been examined as potential redoxmers for ONRFBs. Flavone and isoflavone derivatives were selected as catholyte (positive charge carrier) and anolyte (negative charge carrier) molecules, respectively. To drive their redox potentials to the opposite extremes, in silico derivatization was performed using a novel algorithm to generate a library of > 40000 candidate molecules that penalizes overly complex structures. A multiobjective Bayesian optimization based active learning algorithm was then used to identify best redoxmer candidates in these search spaces. Furthermore, our study provides methodologies for molecular design and optimization of natural scaffolds and highlights the need of incorporating expert chemistry awareness of the natural products and the basic rules of synthetic chemistry in machine learning.

25 ENERGY STORAGE↗

Ligand- and Structure-Based Analysis of Deep Learning-Generated Potential α2a Adrenoceptor Agonists

The α2a adrenoceptor is a particularly medically relevant subtype of the G protein-couple receptor family. Unfortunately, high-throughput techniques aimed at producing novel drug leads targeted to this receptor have been largely unsuccessful due to the complex pharmacology of adrenergic receptors. As such, cutting-edge in silico ligand- and structure-based assessment and de novo deep learning methods are well-positioned to provide new insights into protein-ligand interactions and potential active compounds. In this work, we (i) collect a dataset of α2a adrenoceptor agonists and provide it as a resource for the drug design community; (ii) use the dataset as a basis to generate candidate active structures via deep learning; and (iii) apply computational ligand- and structure-based analysis techniques to gain new insight into α2a adrenoceptor agonists and assess the quality of the computer-generated compounds. Additionally, we further describe how such assessment techniques can be applied to putative chemical probes with a case study involved proposed medetomidine-based probes.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Genomic surveillance indicates clonal replacement of hypervirulent Klebsiella pneumoniae ST881 and ST29 lineage strains in vivo

The emergence of hypervirulent Klebsiella pneumoniae (hvKp) poses a significant public health threat, particularly regarding its carriage in the healthy population. However, the genomic epidemiological characteristics and population dynamics of hvKp within a single patient across distinct infection episodes remain largely unknown. This study aimed to investigate the clonal replacement of hvKp K2-ST881 and K54-ST29 lineage strains in a single patient experiencing multiple-site infections during two independent episodes. Two strains, designated EDhvKp-1 and EDhvKp-2, were obtained from blood and cerebrospinal fluid during the first admission, and the strain isolated from blood on the second admission was named EDhvKp-3. Whole-genome sequencing, utilizing both short-read Illumina and long-read Oxford Nanopore platforms, was conducted. In silico multilocus sequence typing (MLST), identification of antimicrobial resistance and virulence genes, and the phylogenetic relationship between our strains and other K. pneumoniae ST881 and ST29 genomes retrieved from the public database were performed. Virulence potentials were assessed through a mouse lethality assay. Our study indicated that the strains were highly susceptible to multiple antimicrobial agents. Plasmid sequence analysis confirmed that both virulence plasmids, pEDhvKp-1 (166,008 bp) and pEDhvKp-3 (210,948 bp), belonged to IncFIB type. Multiple virulence genes, including rmpA , rmpA2 , rmpC , rmpD , iroBCDN , iucABCD , and iutA , were identified. EDhvKp-1 and EDhvKp-2 showed the closest relationship to strain 502 (differing by 51 SNPs), while EDhvKp-3 exhibited 69 SNPs differences compared to strain TAKPN-1, which all recovered from Chinese patients in 2020. In the mouse infection experiment, both ST881 EDhvKp-1 and ST29 EDhvKp-3 displayed similar virulence traits, causing 90 and 100% of the mice to die within 72 h after intraperitoneal infection, respectively. Our study expands the spectrum of hvKp lineages and highlights genomic alterations associated with clonal switching between two distinct lineages of hvKP that successively replaced each other in vivo . The development of novel strategies for the surveillance, diagnosis, and treatment of high-risk hvKp is urgently needed.

Liu, Ning↗

Predicting partner fitness based on spatial structuring in a light-driven microbial community

Microbial communities have vital roles in systems essential to human health and agriculture, such as gut and soil microbiomes, and there is growing interest in engineering designer consortia for applications in biotechnology ( e . g ., personalized probiotics, bioproduction of high-value products, biosensing). The capacity to monitor and model metabolite exchange in dynamic microbial consortia can provide foundational information important to understand the community level behaviors that emerge, a requirement for building novel consortia. Where experimental approaches for monitoring metabolic exchange are technologically challenging, computational tools can enable greater access to the fate of both chemicals and microbes within a consortium. In this study, we developed an in-silico model of a synthetic microbial consortia of sucrose-secreting Synechococcus elongatus PCC 7942 and Escherichia coli W. Our model was built on the NUFEB framework for Individual-based Modeling (IbM) and optimized for biological accuracy using experimental data. We showed that the relative level of sucrose secretion regulates not only the steady-state support for heterotrophic biomass, but also the temporal dynamics of consortia growth. In order to determine the importance of spatial organization within the consortium, we fit a regression model to spatial data and used it to accurately predict colony fitness. We found that some of the critical parameters for fitness prediction were inter-colony distance, initial biomass, induction level, and distance from the center of the simulation volume. We anticipate that the synergy between experimental and computational approaches will improve our ability to design consortia with novel function.

59 BASIC BIOLOGICAL SCIENCES↗

Role of the Mineral in the Self-Healing of Cracks in Human Enamel

Human enamel is an incredibly resilient biological material, withstanding repeated daily stresses for decades. The mechanisms behind this resilience remain an open question, with recent studies demonstrating a crack-deflection mechanism contributing to enamel toughness and other studies detailing the roles of the organic matrix and remineralization. Here, we focus on the mineral and hypothesize that self-healing of cracks in enamel nanocrystals may be an additional mechanism acting to prevent catastrophic failure. To test this hypothesis, we used a molecular dynamics (MD) approach to compare the fracture behavior of hydroxyapatite (HAP) and calcite, the main minerals in human enamel and sea urchin teeth, respectively. In this work, we find that cracks heal under pressures typical of mastication by fusion of crystals in HAP but not in calcite, which is consistent with the resilience of HAP enamel that calcite teeth lack. Scanning transmission electron microscopy (STEM) images of structurally intact ("sound") human enamel show dashed-line nanocracks that resemble and therefore might be the cracks healed by fusion of crystals produced in silico . The fast, self-healing mechanism shown here is common in soft materials and ceramics but has not been observed in single crystalline materials at room temperature. The crack self-healing in sound enamel nanocrystals, therefore, is unique in the human body and unique in materials science, with potential applications in designing bioinspired materials.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Knowledge graph-aided Bayesian active learning for top- K genetic interaction discovery

In silico methods for predicting the effects of multi-gene perturbations hold great promise for advancing functional genomics, computational drug discovery, and disease modeling. However, the development of these predictive algorithms for mammalian systems has been hampered by limited datasets and high experimental costs. In this study, we present a Bayesian active learning framework designed to discover pairwise host gene knockdowns that effectively inhibit viral proliferation in an in vitro HIV-1 infection model. Our method leverages a biological knowledge graph as side information and employs a computationally efficient batch diversification approach. We evaluated this framework using a dataset of viral load measurements obtained from multi-day dual-gene depletion experiments, encompassing all possible pairwise knockdowns of over 350 host genes associated with HIV infection. We demonstrate that our framework rapidly identifies the most effective gene knockdown pairs for reducing viral load. Furthermore, we show that incorporating side information enhances performance during the early stages of active learning (low data regime), while our batch diversification strategy significantly boosts performance in later stages (high data regime). This framework is general and can be adapted to explore gene interactions in other contexts, such as synthetic lethality prediction and mapping epistatic effects across quantitative trait loci.

Computational biology and bioinformatics↗

Acetyl-CoA Synthesis Through a Bicyclic Carbon-Fixing Pathway in Gas-Fermenting Bacteria

Gas-fermenting acetogens can upgrade one-carbon (C1) compounds (such as CO2 and CO) to the two-carbon (C2) metabolite acetyl coenzyme A (CoA) and convert sugar feedstocks to acetyl-CoA with minimal CO2 emissions. Fulfilling the biosynthetic potential of these microbes requires overcoming challenges in pathway engineering. Here we design a synthetic acetyl-CoA bi-cycle-in addition to the natural carbon-fixing pathways-for C2 metabolite synthesis. This pathway produces an acetyl-CoA by fixation of two CO2 equivalents via three functional modules acting in sequence: carbon fixation, gluconeogenesis and non-oxidative glycolysis. The pathway was examined by in silico thermodynamic and kinetic analyses. The prototypic pathway was implemented in a syngas-fermenting organism, Clostridium ljungdahlii DSM 13528, by expressing a heterologous phosphoketolase that can work with other native enzymes in the host acetogen. The carbon conversion pathway is possible under various growth conditions and is independent of the Wood-Ljungdahl pathway for the valorization of H2 and CO2. This study reports the improvement of carbon conversion using a reductive acetyl-CoA bi-cycle and the potential impact of redox homoeostasis in the acetogenic host for industrial applications of gas fermentation.

Acetyl-CoA↗

GPU-Accelerated Drug Discovery with Docking on the Summit Supercomputer: Porting, Optimization, and Application to COVID-19 Research

Protein-ligand docking is an in silico tool used to screen potential drug compounds for their ability to bind to a given protein receptor within a drug-discovery campaign. Experimental drug screening is expensive and time consuming, and it is desirable to carry out large scale docking calculations in a high-throughput manner to narrow the experimental search space. Few of the existing computational docking tools were designed with high performance computing in mind. Therefore, optimizations to maximize use of high-performance computational resources available at leadership-class computing facilities enables these facilities to be leveraged for drug discovery. Here we present the porting, optimization, and validation of the AutoDock-GPU program for the Summit supercomputer, and its application to initial compound screening efforts to target proteins of the SARS-CoV-2 virus responsible for the current COVID-19 pandemic.

LeGrand, Scott↗

Understanding interactions between biomolecules and two-dimensional nanomaterials using in silico microscopes

We report two-dimensional (2D) nanomaterials such as graphene are increasingly used in research and industry for various biomedical applications. Extensive experimental and theoretical studies have revealed that 2D nanomaterials are promising drug delivery vehicles, yet certain materials exhibit toxicity under biological conditions. So far, it is known that 2D nanomaterials possess strong adsorption propensities for biomolecules. To mitigate potential toxicity and retain favorable physical and chemical properties of 2D nanomaterials, it is necessary to explore the underlying mechanisms of interactions between biomolecules and nanomaterials for the subsequent design of biocompatible 2D nanomaterials for nanomedicine. The purpose of this review is to integrate experimental findings with theoretical observations and facilitate the study of 2D nanomaterial interaction with biomolecules at the molecular level. We discuss the current understanding and progress of 2D nanomaterial interaction with proteins, lipid membranes, and DNA based on molecular dynamics (MD) simulation. In this review, we focus on the 2D graphene nanosheet and briefly discuss other 2D nanomaterials. With the ever-growing computing power, we can image nanoscale processes using MD simulation that are otherwise not observable in experiment. We expect that molecular characterization of the complex behavior between 2D nanomaterials and biomolecules will help fulfill the goal of designing effective 2D nanomaterials as drug delivery platforms.

2D nanomaterials↗

Molecular mechanism of secreted amyloid-β precursor protein in binding and modulating GABABR1a

A recent phenomenal study discovered that the extension domain of secreted amyloid-β precursor protein (sAPP) can bind to the intrinsically disordered sushi 1 domain of the γ-aminobutyric acid type B receptor subunit 1a (GABA B R1a) and modulate its synaptic transmission. The work provided an important structural foundation for the modulation of GABA B R1a; however, the detailed molecular interaction mechanism, crucial for future drug design, remains elusive. Here, we further investigated the dynamical interactions between sAPP peptides and the natively unstructured sushi 1 domain using all-atom molecular dynamics simulations, for both the 17-residue sAPP peptide (APP 17-mer) and its minimally active 9 residue segment (APP 9-mer). We then explored mutations of the APP 9-mer with rigorous free energy perturbation (FEP) calculations. Our in silico mutagenesis studies revealed key residues (D4, W6, and W7) responsible for the binding with the sushi 1 domain. More importantly, one double mutation based on different vertebrate APP sequences from evolution exhibited a stronger binding (ΔΔG = -1.91 ± 0.66 kcal mol -1 ), indicating a potentially enhanced GABA B R1a modulator. These large-scale simulations may provide new insights into the binding mechanism between sAPP and the sushi 1 domain, which could open new avenues in the development of future GABA B R1a-specific therapeutics.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

In Silico Discovery of Multistep Chemistry Initiated by a Conical Intersection: The Challenging Case of Donor–Acceptor Stenhouse Adducts

Detailed mechanistic understanding of multistep chemical reactions triggered by internal conversion via a conical intersection is a challenging task that emphasizes limitations in theoretical and experimental techniques. We present a discovery-based, hypothesis-free computational approach based on first-principles molecular dynamics to discover and refine the switching mechanism of donor–acceptor Stenhouse adducts (DASAs). We simulate the photochemical experiment in silico, following the “hot” ground state dynamics for 10 ps after photoexcitation. Using state-of-the-art graphical processing units-enabled electronic structure calculations we performed in total ~2 ns of nonadiabatic ab initio molecular dynamics discovering (a) critical intermediates that are involved in the open-to-closed transformation, (b) several competing pathways which lower the overall switching yield, and (c) key elements for future design strategies. Finally, our dynamics describe the natural evolution of both the nuclear and electronic degrees of freedom that govern the interconversion between DASA ground-state intermediates, exposing significant elements for future design strategies of molecular switches.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Building and experimenting with an agent-based model to study the population-level impact of CommunityRx, a clinic-based community resource referral intervention

CommunityRx (CRx), an information technology intervention, provides patients with a personalized list of healthful community resources (HealtheRx). In repeated clinical studies, nearly half of those who received clinical “doses” of the HealtheRx shared their information with others (“social doses”). Clinical trial design cannot fully capture the impact of information diffusion, which can act as a force multiplier for the intervention. Furthermore, experimentation is needed to understand how intervention delivery can optimize social spread under varying circumstances. To study information diffusion from CRx under varying conditions, we built an agent-based model (ABM). This study describes the model building process and illustrates how an ABM provides insight about information diffusion through in silico experimentation. To build the ABM, we constructed a synthetic population (“agents”) using publicly-available data sources. Using clinical trial data, we developed empirically-informed processes simulating agent activities, resource knowledge evolution and information sharing. Using RepastHPC and chiSIM software, we replicated the intervention in silico, simulated information diffusion processes, and generated emergent information diffusion networks. The CRx ABM was calibrated using empirical data to replicate the CRx intervention in silico. We used the ABM to quantify information spread via social versus clinical dosing then conducted information diffusion experiments, comparing the social dosing effect of the intervention when delivered by physicians, nurses or clinical clerks. The synthetic population (N = 802,191) exhibited diverse behavioral characteristics, including activity and knowledge evolution patterns. In silico delivery of the intervention was replicated with high fidelity. Large-scale information diffusion networks emerged among agents exchanging resource information. Varying the propensity for information exchange resulted in networks with different topological characteristics. Community resource information spread via social dosing was nearly 4 fold that from clinical dosing alone and did not vary by delivery mode. This study, using CRx as an example, demonstrates the process of building and experimenting with an ABM to study information diffusion from, and the population-level impact of, a clinical information-based intervention. While the focus of the CRx ABM is to recreate the CRx intervention in silico, the general process of model building, and computational experimentation presented is generalizable to other large-scale ABMs of information diffusion.

59 BASIC BIOLOGICAL SCIENCES↗

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↗

Deciphering the structure of Arabidopsis thaliana 5-enol-pyruvyl-shikimate-3-phosphate synthase: An essential step toward the discovery of novel inhibitors to supersede glyphosate

Glyphosate interferes with plant aromatic metabolism through the inhibition of 5-enol-pyruvyl-shikimate-3-phosphate (EPSP) synthase [EPSPS, EC 2.5.1.19]. For this reason, EPSPS has been extensively studied in a vast array of organisms. This notwithstanding, up to date, the crystal structure of the protein has been solved exclusively in a few prokaryotes, while that of the plant enzyme has been only deduced in silico by similarity. This study aimed at determining the structure of EPSPS from the plant model species Arabidopsis thaliana, which has been cloned, heterologously expressed and affinity-purified. The kinetic properties of the enzyme have been determined, as well as its susceptibility to the inhibition brought about by glyphosate. The crystal structure of the protein has been resolved at high resolution (1.4 Å), showing open conformation of the enzyme, which is the state ready for substrate/inhibitor binding. This provides a framework for the structure-based design of novel EPSPS inhibitors. Surface regions near the active-site cleft entrance or at the interdomain hinge appear promising for inhibitor selectivity, while bound chloride near the active site is a potential placeholder for anionic moieties of future herbicides.

59 BASIC BIOLOGICAL SCIENCES↗

Modeling the impact of structure and coverage on the reactivity of realistic heterogeneous catalysts

Adsorbates often cover the surfaces of catalysts densely as they carry out reactions, dynamically altering their structure and reactivity. Understanding adsorbate-induced phenomena and harnessing them in our broader quest for improved catalysts is a substantial challenge that is only beginning to be addressed. Here, in this work, we chart a path toward a deeper understanding of such phenomena by focusing on emerging in silico modeling methodologies, which will increasingly incorporate machine learning techniques. We first examine how adsorption on catalyst surfaces can lead to local and even global structural changes spanning entire nanoparticles, and how this affects their reactivity. We then evaluate current efforts and the remaining challenges in developing robust and predictive simulations for modeling such behavior. Last, we provide our perspectives in four critical areas—integration of artificial intelligence, building robust catalysis informatics infrastructure, synergism with experimental characterization, and adaptive modeling frameworks—that we believe can help surmount the remaining challenges in rationally designing catalysts in light of these complex phenomena.

catalytic mechanisms↗