Search NASA⌕ Search

SEARCH · Search NASA

Results for “in silico”

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 235 records · Page 13
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

High-throughput and data-driven search for stable optoelectronic AMSe 3 materials

The rapid advancement in emerging optoelectronic technologies demands highly efficient, affordable, and ecofriendly materials. In this context, ternary chalcogenides, especially ternary selenides, show early promise as a material class due to their stability and remarkable electronic, optical, and transport properties. In this work, we integrate first-principles-based high-throughput computations with machine learning (ML) techniques to predict the thermodynamic stability and optoelectronic properties of 920 valency-satisfied selenide compounds. Through investigating polymorphism, our study reveals the edge-sharing orthorhombic Pnma phase (NH 4 CdCl 3 -type) as the most stable structure for most ternary selenides. High-fidelity supervised ML models are trained and tested to accelerate stability and band gap predictions. These data-driven models pin down the most influential features that dominantly control key material characteristics. The multistep high-throughput computations identify the ternary selenides with optimal direct band gaps, light carrier masses, and strong optical absorption edges. The extensive materials screening considering phase stability, toxicity, and defect tolerance, finally identifies the seven most suitable candidates for photovoltaic applications. Two of these final compounds, SrZrSe 3 and SrHfSe 3 , have already been synthesized in a single-phase form, with the latter showing an optically suitable band gap, aligning well with our findings. The non-adiabatic molecular dynamics reveal sufficiently long photoexcited charge carrier lifetimes (on the order of nanoseconds) in some of these selected selenide materials, indicating their exciting characteristics. Overall, our study suggests a robust in silico framework that can be extended to screen large datasets of various material classes for identifying promising photoactive candidates.

36 MATERIALS SCIENCE↗
DOE OSTI2025

The influence of protein electrostatics on potential inversion in flavoproteins

Biology uses relatively few electron-transfer cofactors, tuning their potentials, electronic couplings, and reorganization energies to carry out the required chemistry. It is remarkable that the potential ordering of two-electron transfer active flavins can be normal (first oxidation at low potential and second oxidation at high potential) or inverted, and the gap between the potentials can be as large as one volt. Analysis based on structural bioinformatics and electrostatics indicates that the ordering of the flavin redox potential is influenced by protein electrostatics. In all 36 flavoproteins examined, the introduction of a negative charge near the flavin in silico increases the extent of potential inversion (by lowering the electrochemical potential of the second electron-transfer step); the introduction of a positive charge near the flavin favors normally ordered potentials. We also find that the addition of positive charges increases the electrochemical potential for the naturally occurring one-electron transition in flavodoxins (between deprotonated hydroquinone and neutral semiquinone) and also increases the second one-electron transition in bifurcating flavins (between anionic semiquinone and fully oxidized flavin). Finally, we find that proximity of a proton acceptor, notably conserved arginine, supports proton-coupled electron transfer because it may act as a proton acceptor, promoting potential inversion. This key arginine residue may enable two-electron transfer chemistry by promoting the proton-coupled electron transfer process over the pure electron transfer process, suggesting how a protein's flavin environment may influence one- or two-electron chemistry in flavoproteins.

Singh, Niven [Duke Univ., Durham, NC (United State↗
DOE OSTI2022

Insights into pathological mechanisms and interventions revealed by analyzing a mathematical model for cone metabolism

This work analyzes a mathematical model for the metabolic dynamics of a cone photoreceptor, which is the first model to account for energy generation from fatty acids oxidation of shed photoreceptor outer segments (POS). Multiple parameter bifurcation analysis shows that joint variations in external glucose, the efficiency of glucose transporter 1 (GLUT1), lipid utilization for POS renewal, and oxidation of fatty acids affect the cone’s metabolic vitality and its capability to adapt under glucose-deficient conditions. The analysis further reveals that when glucose is scarce, cone viability cannot be sustained by only fueling energy production in the mitochondria, but it also requires supporting anabolic processes to create lipids necessary for cell maintenance and repair. In silico experiments are used to investigate how the duration of glucose deprivation impacts the cell without and with a potential GLUT1 or oxidation of fatty acids intervention as well as a dual intervention. The results show that for prolonged duration of glucose deprivation, the cone metabolic system does not recover with higher oxidation of fatty acids and requires greater effectiveness of GLUT1 to recover. Finally, time-varying global sensitivity analysis (GSA) is applied to assess the sensitivity of the model outputs of interest to changes and uncertainty in the parameters at specific times. The results reveal a critical temporal window where there would be more flexibility for interventions to rescue a cone cell from the detrimental consequences of glucose shortage.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Inverse methods for design of soft materials

Functional soft materials, comprising colloidal and molecular building blocks that self-organize into complex structures as a result of their tunable interactions, enable a wide array of technological applications. Inverse methods provide a systematic means for navigating their inherently high-dimensional design spaces to create materials with targeted properties. Furthermore, while multiple physically motivated inverse strategies have been successfully implemented in silico, their translation to guiding experimental materials discovery has thus far been limited to a handful of proof-of-concept studies. In this perspective, we discuss recent advances in inverse methods for design of soft materials that address two challenges: (1) methodological limitations that prevent such approaches from satisfying design constraints and (2) computational challenges that limit the size and complexity of systems that can be addressed. Strategies that leverage machine learning have proven particularly effective, including methods to discover order parameters that characterize complex structural motifs and schemes to efficiently compute macroscopic properties from the underlying structure. We also highlight promising opportunities to improve the experimental realizability of materials designed computationally, including discovery of materials with functionality at multiple thermodynamic states, design of externally directed assembly protocols that are simple to implement in experiments, and strategies to improve the accuracy and computational efficiency of experimentally relevant models.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Deterministic modeling of hybrid nonlinear effects in epsilon-near-zero thin films

In nonlinear optics, significant effort is concentrated on improving the strength and efficiency of interactions; however, experimentally investigating nonlinear materials is a complex, time-consuming, and costly investment. Moreover, it is often challenging to isolate, study, and optimize material parameters in an experiment due to complexities in the growth process. Recently, epsilon-near-zero materials have received a great deal of attention as promising nonlinear optical materials, but like many up-and-coming materials, the ability to explore and optimize their properties has been challenging. Here, we establish a framework to rapidly evaluate the performance of nonlinear epsilon-near-zero materials for both inter- and intraband effects in silico, requiring only an energy-momentum (E-k) diagram, linear optical properties, and experimental conditions. Measured nonlinear reflection and transmission in gallium-doped zinc oxide films are compared to the numerical framework for both intra- and interband excitation to verify accuracy across wavelength and irradiance while two figures of merit (FoMs) are introduced to quickly evaluate the performance of films without a full numerical framework. This capability is used to predict the performance of highly doped gallium nitride, cadmium oxide, zinc oxide, and indium tin oxide films, and efficient intra- and interband operation conditions are identified. Through this numerical framework and the FoMs, the exploration of unstudied epsilon-near-zero materials is enabled without the need for a nonlinear experiment, thereby accelerating the search for more efficient nonlinear materials and excitation conditions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

High-throughput screening of tribological properties of monolayer films using molecular dynamics and machine learning

Monolayer films have shown promise as a lubricating layer to reduce friction and wear of mechanical devices with separations on the nanoscale. These films have a vast design space with many tunable properties that can affect their tribological effectiveness. For example, terminal group chemistry, film composition, and backbone chemistry can all lead to films with significantly different tribological properties. This design space, however, is very difficult to explore without a combinatorial approach and an automatable, reproducible, and extensible workflow to screen for promising candidate films. Here, using the Molecular Simulation Design Framework (MoSDeF), a combinatorial screening study was performed to explore 9747 unique monolayer films (116 964 total simulations) and a machine learning (ML) model using a random forest regressor, an ensemble learning technique, to explore the role of terminal group chemistry and its effect on tribological effectiveness. The most promising films were found to contain small terminal groups such as cyano and ethylene. The ML model was subsequently applied to screen terminal group candidates identified from the ChEMBL small molecule library. Approximately 193 131 unique film candidates were screened with approximately a five order of magnitude speed-up in analysis compared to simulation alone. The ML model was thus able to be used as a predictive tool to greatly speed up the initial screening of promising candidate films for future simulation studies, suggesting that computational screening in combination with ML can greatly increase the throughput in combinatorial approaches to generate in silico data and then train ML models in a controlled, self-consistent fashion.

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

Effects of elevated-temperature deposition on the atomic structure of amorphous Ta 2 O 5 films

Brownian thermal noise as a result of mechanical loss in optical coatings will become the dominant source of noise at the most sensitive frequencies of ground-based gravitational-wave detectors. Experiments found, however, that a candidate material, amorphous Ta 2 O 5 , is unable to form an ultrastable glass and, consequently, to yield a film with significantly reduced mechanical loss through elevated-temperature deposition alone. X-ray scattering PDF measurements are carried out on films deposited and subsequently annealed at various temperatures. Inverse atomic modeling is used to analyze the short and medium range features in the atomic structure of these films. Furthermore, in silico deposition simulations of Ta 2 O 5 are carried out at various substrate temperatures and an atomic level analysis of the growth at high temperatures is presented. It is observed that upon elevated-temperature deposition, short range features remain identical, whereas medium range order increases. After annealing, however, both the short and medium range orders of films deposited at different substrate temperatures are nearly identical. A discussion on the surface diffusion and glass transition temperatures indicates that future pursuits of ultrastable low-mechanical-loss films through elevated temperature deposition should focus on materials with a high surface mobility, and/or lower glass transition temperatures in the range of achievable deposition temperatures.

36 MATERIALS SCIENCE↗
DOE OSTI2024

A transferable classical force field to describe glyme based lithium solvate ionic liquids

A non-polarizable force field for lithium (Li + ) and bis(trifluoromethanesulfonyl)imide (TFSI – ) ions solvated in diglyme at around 0.2 mol fraction salt concentration was developed based on ab initio molecular dynamics (AIMD) simulations and a modified polymer consistent force field model. A force–torque matching based scheme, in conjunction with a genetic algorithm, was used to determine the Lennard-Jones (LJ) parameters of the ion–ion and ion–solvent interactions. This force field includes a partial charge scaling factor and a scaling factor for the 1–4 interactions. The resulting force field successfully reproduces the radial distribution function of the AIMD simulations and shows better agreement compared to the unmodified force field. The new force field was then used to simulate salt solutions with glymes of increasing chain lengths and different salt concentrations. The comparison of the MD simulations, using the new force field, with experimental data at different salt concentrations and AIMD simulations on equimolar concentrations of the triglyme system demonstrates the transferability of the force field parameters to longer glymes and higher salt concentrations. Furthermore, the force field appears to reproduce the features of the experimental x-ray structure factors, suggesting accuracy beyond the first solvation shell, for equimolar salt solutions using both triglyme and tetraglyme as the solvent. Altogether, the new force field was found to accurately reproduce the molecular descriptions of LiTFSI-glyme systems not only at various salt concentrations but also with glymes of different chain lengths. Thus, the new force field provides a useful and accurate tool to perform in silico studies of this family of systems at the atomistic level.

25 ENERGY STORAGE↗
DOE OSTI2026

Many roads to the seam: How conformational flexibility drives nonadiabatic relaxation in a prototypical tetrapyrrolic chromophore

Large and structurally flexible chromophores pose challenges for in silico modeling of photodeactivation due to the many vibrational modes that can funnel the system toward energy degeneracy. In this work, we examine how the multiple degrees of freedom in biliverdin, a prototypical tetrapyrrolic chromophore, cooperate to drive access to the S 1 /S 0 intersection seam in vacuo. We begin by mapping the ground-state potential energy surface to identify representative biliverdin conformers relevant to photoexcitation. We then use a CASSCF-based framework to map the excited-state landscape and characterize the intersection seam, identifying distinct conical-intersection types. Finally, we employ ab initio multiple spawning to resolve the dynamical pathways by which the system accesses these regions. DFT potential-energy and free-energy mappings indicate that, although several conformers are relevant, the “locked-helix” ZsZsZs conformer predominates in the ground state. The intersection seam comprises numerous geometrically distinct regions characterized by varying degrees and combinations of dihedral torsion, pyramidalization, and bond-length alternation. Yet only select regions lie within energetic reach, and moderate barriers separate them from the S 1 minimum. Nonadiabatic dynamics combined with multivariate analyses show that, despite extensive mode coupling during deactivation that guides the system toward multiple regions of the seam, a single dihedral torsion, together with bond-length alternation, predominantly drives energy degeneracy. This work offers new insight into biliverdin’s intrinsic photochemical response and underscores a general feature of flexible chromophores: many modes may participate during photorelaxation, but only a limited subset ultimately dictates seam accessibility.

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

Molecular characterization of a fungal gasdermin-like protein

Programmed cell death (PCD) in filamentous fungi prevents cytoplasmic mixing following fusion between conspecific genetically distinct individuals (allorecognition) and serves as a defense mechanism against mycoparasitism, genome exploitation, and deleterious cytoplasmic elements (i.e., senescence plasmids). Recently, we identified regulator of cell death-1 ( rcd-1 ), a gene controlling PCD in germinated asexual spores in the filamentous fungus Neurospora crassa . rcd-1 alleles are highly polymorphic and fall into two haplogroups in N. crassa populations. Coexpression of alleles from the two haplogroups, rcd-1–1 and rcd-1–2 , is necessary and sufficient to trigger a cell death reaction. Here, we investigated the molecular bases of rcd-1 -dependent cell death. Based on in silico analyses, we found that RCD-1 is a remote homolog of the N-terminal pore-forming domain of gasdermin, the executioner protein of a highly inflammatory cell death reaction termed pyroptosis, which plays a key role in mammalian innate immunity. We show that RCD-1 localizes to the cell periphery and that cellular localization of RCD-1 was correlated with conserved positively charged residues on predicted amphipathic α-helices, as shown for murine gasdermin-D. Similar to gasdermin, RCD-1 binds acidic phospholipids in vitro, notably, cardiolipin and phosphatidylserine, and interacts with liposomes containing such lipids. The RCD-1 incompatibility system was reconstituted in human 293T cells, where coexpression of incompatible rcd-1–1/rcd-1–2 alleles triggered pyroptotic-like cell death. Oligomers of RCD-1 were associated with the cell death reaction, further supporting the evolutionary relationship between gasdermin and rcd-1 . This report documents an ancient transkingdom relationship of cell death execution modules involved in organismal defense.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Entropic formation of a thermodynamically stable colloidal quasicrystal with negligible phason strain

Quasicrystals have been discovered in a variety of materials ranging from metals to polymers. Yet, why and how they form is incompletely understood. In situ transmission electron microscopy of alloy quasicrystal formation in metals suggests an error-and-repair mechanism, whereby quasiperiodic crystals grow imperfectly with phason strain present, and only perfect themselves later into a high-quality quasicrystal with negligible phason strain. The growth mechanism has not been investigated for other types of quasicrystals, such as dendrimeric, polymeric, or colloidal quasicrystals. Soft-matter quasicrystals typically result from entropic, rather than energetic, interactions, and are not usually grown (either in laboratories or in silico) into large-volume quasicrystals. Consequently, it is unknown whether soft-matter quasicrystals form with the high degree of structural quality found in metal alloy quasicrystals. Here, we investigate the entropically driven growth of colloidal dodecagonal quasicrystals (DQCs) via computer simulation of systems of hard tetrahedra, which are simple models for anisotropic colloidal particles that form a quasicrystal. Using a pattern recognition algorithm applied to particle trajectories during DQC growth, we analyze phason strain to follow the evolution of quasiperiodic order. As in alloys, we observe high structural quality; DQCs with low phason strain crystallize directly from the melt and only require minimal further reduction of phason strain. Finally, we also observe transformation from a denser approximant to the DQC via continuous phason strain relaxation. Our results demonstrate that soft-matter quasicrystals dominated by entropy can be thermodynamically stable and grown with high structural quality––just like their alloy quasicrystal counterparts.

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

Closing the loop between microstructure and charge transport in conjugated polymers by combining microscopy and simulation

A grand challenge in materials science is to identify the impact of molecular composition and structure across a range of length scales on macroscopic properties. We demonstrate a unified experimental–theoretical framework that coordinates experimental measurements of mesoscale structure with molecular-level physical modeling to bridge multiple scales of physical behavior. Here we apply this framework to understand charge transport in a semiconducting polymer. Spatially-resolved nanodiffraction in a transmission electron microscope is combined with a self-consistent framework of the polymer chain statistics to yield a detailed picture of the polymer microstructure ranging from the molecular to device relevant scale. Using these data as inputs for charge transport calculations, the combined multiscale approach highlights the underrepresented role of defects in existing transport models. Short-range transport is shown to be more chaotic than is often pictured, with the drift velocity accounting for a small portion of overall charge motion. Local transport is sensitive to the alignment and geometry of polymer chains. At longer length scales, large domains and gradual grain boundaries funnel charges preferentially to certain regions, creating inhomogeneous charge distributions. While alignment generally improves mobility, these funneling effects negatively impact mobility. The microstructure is modified in silico to explore possible design rules, showing chain stiffness and alignment to be beneficial while local homogeneity has no positive effect. This combined approach creates a flexible and extensible pipeline for analyzing multiscale functional properties and a general strategy for extending the accesible length scales of experimental and theoretical probes by harnessing their combined strengths.

42 ENGINEERING↗
DOE OSTI2022

De novo designed protein inhibitors of amyloid aggregation and seeding

Neurodegenerative diseases are characterized by the pathologic accumulation of aggregated proteins. Known as amyloid, these fibrillar aggregates include proteins such as tau and amyloid-β (Aβ) in Alzheimer’s disease (AD) and alpha-synuclein (αSyn) in Parkinson’s disease (PD). The development and spread of amyloid fibrils within the brain correlates with disease onset and progression, and inhibiting amyloid formation is a possible route toward therapeutic development. Recent advances have enabled the determination of amyloid fibril structures to atomic-level resolution, improving the possibility of structure-based inhibitor design. In this work, we use these amyloid structures to design inhibitors that bind to the ends of fibrils, “capping” them so as to prevent further growth. Using de novo protein design, we develop a library of miniprotein inhibitors of 35 to 48 residues that target the amyloid structures of tau, Aβ, and αSyn. Biophysical characterization of top in silico designed inhibitors shows they form stable folds, have no sequence similarity to naturally occurring proteins, and specifically prevent the aggregation of their targeted amyloid-prone proteins in vitro. The inhibitors also prevent the seeded aggregation and toxicity of fibrils in cells. In vivo evaluation reveals their ability to reduce aggregation and rescue motor deficits in Caenorhabditis elegans models of PD and AD.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Iterative computational design and crystallographic screening identifies potent inhibitors targeting the Nsp3 macrodomain of SARS-CoV-2

The nonstructural protein 3 (NSP3) of the severe acute respiratory syndrome-coronavirus-2 (SARS-CoV-2) contains a conserved macrodomain enzyme (Mac1) that is critical for pathogenesis and lethality. While small-molecule inhibitors of Mac1 have great therapeutic potential, at the outset of the COVID-19 pandemic, there were no well-validated inhibitors for this protein nor, indeed, the macrodomain enzyme family, making this target a pharmacological orphan. Here, we report the structure-based discovery and development of several different chemical scaffolds exhibiting low- to sub-micromolar affinity for Mac1 through iterations of computer-aided design, structural characterization by ultra-high-resolution protein crystallography, and binding evaluation. Potent scaffolds were designed with in silico fragment linkage and by ultra-large library docking of over 450 million molecules. Both techniques leverage the computational exploration of tangible chemical space and are applicable to other pharmacological orphans. Overall, 160 ligands in 119 different scaffolds were discovered, and 153 Mac1-ligand complex crystal structures were determined, typically to 1 Å resolution or better. Our analyses discovered selective and cell-permeable molecules, unexpected ligand-mediated conformational changes within the active site, and key inhibitor motifs that will template future drug development against Mac1.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2023

Identification of a broadly conserved family of enzymes that hydrolyze (p)ppApp

Bacteria produce a variety of nucleotide second messengers to adapt to their surroundings. Although chemically similar, the nucleotides guanosine penta- and tetraphosphate [(p)ppGpp] and adenosine penta- and tetraphosphate [(p)ppApp] have distinct functions in bacteria. (p)ppGpp mediates survival under nutrient-limiting conditions and its intracellular levels are regulated by synthetases and hydrolases belonging to the RelA-SpoT homolog (RSH) family of enzymes. By contrast, (p)ppApp is not known to be involved in nutrient stress responses and is synthesized by RSH-resembling toxins that inhibit the growth of bacterial cells. However, it remains unclear whether there exists a family of hydrolases that specifically act on (p)ppApp to reverse its toxic effects. Here, we present the structure and biochemical characterization of adenosine 3'-pyrophosphohydrolase 1 (Aph1), the founding member of a monofunctional (p)ppApp hydrolase family of enzymes. Our work reveals that Aph1 adopts a histidine-aspartate (HD)-domain fold characteristic of phosphohydrolase metalloenzymes and its activity mitigates the growth inhibitory effects of (p)ppApp-synthesizing toxins. Using an informatic approach, we identify over 2,000 putative (p)ppApp hydrolases that are widely distributed across bacterial phyla and found in diverse genomic contexts, and we demonstrate that 12 representative members hydrolyze ppApp. In addition, our in silico analyses reveal a unique molecular signature that is specific to (p)ppApp hydrolases, and we show that mutation of two residues within this signature broadens the specificity of Aph1 to promiscuously hydrolyze (p)ppGpp in vitro. Overall, our findings indicate that like (p)ppGpp hydrolases, (p)ppApp hydrolases are widespread in bacteria and may play important and underappreciated role(s) in bacterial physiology.

(p)ppApp↗
DOE OSTI2024

Integrating chromosome conformation and DNA repair in a computational framework to assess cell radiosensitivity

Objective. The arrangement of chromosomes in the cell nucleus has implications for cell radiosensitivity. The development of new tools to utilize Hi-C chromosome conformation data in nanoscale radiation track structure simulations allows for in silico investigation of this phenomenon. We have developed a framework employing Hi-C-based cell nucleus models in Monte Carlo radiation simulations, in conjunction with mechanistic models of DNA repair, to predict not only the initial radiation-induced DNA damage, but also the repair outcomes resulting from this damage, allowing us to investigate the role chromosome conformation plays in the biological outcome of radiation exposure. Approach. In this study, we used this framework to generate cell nucleus models based on Hi-C data from fibroblast and lymphoblastoid cells and explore the effects of cell type-specific chromosome structure on radiation response. The models were used to simulate external beam irradiation including DNA damage and subsequent DNA repair. The kinetics of the simulated DNA repair were compared with previous results. Main results. We found that the fibroblast models resulted in a higher rate of inter-chromosome misrepair than the lymphoblastoid model, despite having similar amounts of initial DNA damage and total misrepairs for each irradiation scenario. Significance. This framework represents a step forward in radiobiological modeling and simulation allowing for more realistic investigation of radiosensitivity in different types of cells.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Skeletides: A Modular, Simplified Physical Model of Protein Secondary Structure

Three-dimensional (3D) models are essential for visualization and conceptual understanding of complex architectures such as protein structure. Although there is a plethora of software platforms that allow digital depictions of protein structure at an atomic level in silico, physical models are needed to convey an intuitive understanding of biomolecular architecture. However, it is a challenge to represent all the relevant features of proteins in a single physical model due to their sheer structural complexity. Here, we describe a modular protein model that focuses only on representation of the secondary structure—the underlying structural skeleton. The simplified model consists of amino acid units, which can be linked together to reproduce the two most fundamental structural features of protein secondary structure: the relative positions of the amino acid alpha carbon atoms, and the intra-main chain hydrogen bonding pattern. We use 3D printing and magnets to create a set of three modular amino acid building blocks, which when linked together into a chain, can faithfully represent alpha helices, beta sheets, and turns. These simple to make models can be used to quickly assemble the alpha carbon trace of an entire protein domain, which conveys a tactile experience of the complexity of protein skeletal architecture. These models highlight the modularity of protein structure: where a single structural unit, the amino acid, can be linked together to form a larger, regular secondary structure. These models also have a propensity to spontaneously organize into alpha helices and beta sheets, as demonstrated by their ability to autonomously assemble when placed in a circulating water tank. These models provide a missing educational tool to expand knowledge of protein structure, foster deeper insight into protein folding, and inspire greater interest in biomacromolecular architecture.

36 MATERIALS SCIENCE↗
DOE OSTI2022

MINE 2.0: enhanced biochemical coverage for peak identification in untargeted metabolomics

Abstract Summary Although advances in untargeted metabolomics have made it possible to gather data on thousands of cellular metabolites in parallel, identification of novel metabolites from these datasets remains challenging. To address this need, Metabolic in silico Network Expansions (MINEs) were developed. A MINE is an expansion of known biochemistry which can be used as a list of potential structures for unannotated metabolomics peaks. Here, we present MINE 2.0, which utilizes a new set of biochemical transformation rules that covers 93% of MetaCyc reactions (compared to 25% in MINE 1.0). This results in a 17-fold increase in database size and a 40% increase in MINE database compounds matching unannotated peaks from an untargeted metabolomics dataset. MINE 2.0 is thus a significant improvement to this community resource. Availability and implementation The MINE 2.0 website can be accessed at https://minedatabase.ci.northwestern.edu. The MINE 2.0 web API documentation can be accessed at https://mine-api.readthedocs.io/en/latest/. The data and code underlying this article are available in the MINE-2.0-Paper repository at https://github.com/tyo-nu/MINE-2.0-Paper. MINE 2.0 source code can be accessed at https://github.com/tyo-nu/MINE-Database (MINE construction), https://github.com/tyo-nu/MINE-Server (backend web API) and https://github.com/tyo-nu/MINE-app (web app). Supplementary information Supplementary data are available at Bioinformatics online.

59 BASIC BIOLOGICAL SCIENCES↗