Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “conferences”

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 1,549 records · Page 86
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Preface for the Sherwood Fusion Theory 2023 special collection

Pursuing controlled fusion is a monumental endeavor that promises to unlock a transformational source of clean energy. Here, the International Sherwood Fusion Theory Conference is a crucial platform for researchers to discuss the physics governing plasmas in fusion-relevant configurations. From exploring the fundamental properties of plasma to cutting-edge computer simulation techniques, this conference fosters innovation and collaboration, propelling us closer to realizing the dream of a fusion-powered future.

Beidler, Matthew T. [Oak Ridge National Laboratory↗
DOE OSTI2020

Toward an international practical pressure scale: A proposal for an IPPS ruby gauge (IPPS-Ruby2020)

At the 26 th AIRAPT conference in 2017, a task group was formed to work on an International Practical Pressure Scale (IPPS). This report summarizes the activities of the task group toward an IPPS ruby gauge. We have selected three different approaches to establishing the relation between pressure ( P ) and ruby R1-line shift ( Δλ ) with three groups of optimal reference materials for applying these approaches. Using a polynomial form of the second order, the recommended ruby gauge (referred as Ruby2020) is expressed by: P [ G P a ] = 1.87 ( ± 0.01 ) × 10 3 Δ λ λ 0 1 + 5.63 ( ± 0.03 ) Δ λ λ 0 , where λ 0 is the wavelength of the R1-line near 694.25 nm at ambient condition. In June of 2020, the Executive Committee of AIRAPT endorsed the proposed Ruby2020. Finally, we encourage high-pressure practitioners to utilize Ruby2020 within its applicable pressure range (up to 150 GPa), so that pressure data can be directly compared across laboratories and amended consistently as better scales emerge in the future.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Understanding the Role of Populus ECERIFERUM2-Likes in the Biosynthesis of Very-Long-Chain Fatty Acids for Cuticular Waxes

Cuticular waxes are derived from very-long-chain fatty acid (VLCFA) precursors made by the concerted action of four enzymes that form the fatty acid (FA) elongation complex. Here, the condensing enzyme of the complex confers specificity to substrates of different chain lengths, yet on its own cannot account for the biosynthesis of VLCFAs longer than 28 carbons (C28). Recent evidence from Arabidopsis thaliana points to a synergistic role of clade II BAHD acyltransferases and condensing enzymes in the elongation of VLCFAs beyond C28. In Populus trichocarpa, clade II is composed of seven uncharacterized paralogous genes (PtCER2-like1–7). In the present study, five of these genes were heterologously expressed in yeast and their respective FA profiles were determined. PtCER2-likes differentially altered the accumulation of C28 and C30 FAs when expressed in the presence of the condensing enzyme AtCER6. Among these, PtCER2-like5 produced the highest levels of C28 FAs in yeast and its expression was localized to the epidermis in β-glucuronidase-reporter poplar lines, consistent with a role in cuticular wax biosynthesis. Complementation of the A. thaliana cer2-5 mutant with PtCER2-like5 increased the levels of C28-derived cuticular waxes at the expense of C30-derived components. Together, these results demonstrate that the role of CER2-likes in cuticular wax biosynthesis is conserved in Populus clade II BAHD acyltransferases.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Assessment of Clonal Expansion Using CarcSeq Measurement of Lung Cancer Driver Mutations and Correlation With Mouse Strain- and Sex-Related Incidence of Spontaneous Lung Neoplasia

Quantification of variation in levels of spontaneously occurring cancer driver mutations (CDMs) was developed to assess clonal expansion and predict future risk of neoplasm development. Specifically, an error-corrected next-generation sequencing method, CarcSeq, and a mouse CarcSeq panel (analogous to human and rat panels) were developed and used to quantify low-frequency mutations in a panel of amplicons enriched in hotspot CDMs. Mutations in a subset of panel amplicons, Braf, Egfr, Kras, Stk11, and Tp53, were related to incidence of lung neoplasms at 2 years. This was achieved by correlating median absolute deviation (MAD) from the overall median mutant fraction (MF) measured in the lung DNA of 16-week-old male and female, B6C3F1 and CD-1 mice (10 mice/sex/strain) with percentages of spontaneous alveolar/bronchioloalveolar adenomas and carcinomas reported in bioassay control groups. A total of 1586 mouse lung mutants with MFs > x 10-4 were recovered. The ratio of nonsynonymous to synonymous mutations was used to assess the proportion of recovered mutations conferring a positive selective advantage. The greatest ratio was observed in what is considered the most lung tumor-sensitive model examined, male B6C3F1 mice. Of the recurrent, nonsynonymous mouse mutations recovered, 55.5% have been reported in human tumors, with many located in or around the mouse equivalent of human cancer hotspot codons. MAD for the same subset of amplicons measured in normal human lung DNA samples showed a correlation of moderate strength and borderline significance with age (a cancer risk factor), as well as age-related cumulative lung cancer risk, suggesting MAD may inform species extrapolation.

Toxicology↗
DOE OSTI2025

Braiding for the win: Harnessing braiding statistics in topological states to play quantum games

Nonlocal quantum games provide proof of principle that quantum resources can confer an advantage at certain tasks. They also provide a compelling way to explore the computational utility of phases of matter on quantum hardware. In a recent paper [O. Hart et al., Phys. Rev. Lett. 134, 130602 (2025)], we demonstrated that a toric code resource state conferred advantage at a certain nonlocal game, which remained robust to small deformations of the resource state. In this paper we demonstrate that this robust advantage is a generic property of resource states drawn from topological or fracton ordered phases of quantum matter. To this end, we illustrate how several other states from paradigmatic topological and fracton ordered phases can function as resources for suitably defined nonlocal games, notably the three-dimensional toric-code phase, the X-cube fracton phase, and the double-semion phase. The key in every case is to design a nonlocal game that harnesses the characteristic braiding processes of a quantum phase as a source of contextuality. We unify the strategies that take advantage of mutual statistics by relating the operators to be measured to order and disorder parameters of an underlying generalized symmetry-breaking phase transition. Additionally, by connecting the win probability to twist products, we show that success at the game serves as a many-body entanglement witness. Namely, if the players implement a perfect quantum strategy on large length scales, the quantum state they share cannot be connected to a trivial product state via a constant-depth local unitary circuit. Lastly, we massively generalize the family of games that admit perfect strategies when codewords of homological quantum error-correcting codes are used as resources.

Fractons↗
DOE OSTI2022

Non-Boussinesq subgrid-scale model with dynamic tensorial coefficients

A major drawback of Boussinesq-type subgrid-scale stress models used in large-eddy simulations is the inherent assumption of alignment between large-scale strain rates and filtered subgrid-stresses. A priori analyses using direct numerical simulation (DNS) data have shown that this assumption is invalid locally as subgrid-scale stresses are poorly correlated with the large-scale strain rates [J. Bardina, J. Ferziger, and W. Reynolds, Improved subgrid-scale models for large-eddy simulation, in Proceedings of the 13th Fluid and Plasmadynamics Conference, AIAA (1980); C. Meneveau and K. Katz, Scale-invariance and turbulence models for large-eddy simulation, Ann. Rev. Fluid Mech. 32, 1 (2000)]. In the present work, a new, non-Boussinesq subgrid-scale model is presented where the model coefficients are computed dynamically. Some previous non-Boussinesq models have observed issues in providing adequate dissipation of turbulent kinetic energy [e.g., Bardina et al., Proceedings of the 13th Fluid and Plasmadynamics Conference (1980); R. A. Clark, J. Ferziger, and W.C. Reynolds. Evaluation of subgrid-scale models using an accurately simulated turbulent flow, J. Fluid Mech. 91, 1 (1979); S. Stolz and N. A. Adams, An approximate deconvolution procedure for large-eddy simulation, Phys. Fluids 11, 1699 (1999)]; however, the present model is shown to provide sufficient dissipation using dynamic coefficients. Modeled subgrid-scale Reynolds stresses satisfy the consistency requirements of the governing equations for large-eddy simulation (LES), vanish in laminar flow and at solid boundaries, and have the correct asymptotic behavior in the near-wall region of a turbulent boundary layer. The new model, referred to as the dynamic tensor-coefficient Smagorinsky model (DTCSM), has been tested in simulations of canonical flows: decaying and forced homogeneous isotropic turbulence, and wall-modeled turbulent channel flow at high Reynolds numbers. The results show favorable agreement with DNS data. It has been shown that DTCSM offers similar predictive capabilities as the dynamic Smagorinsky model for canonical flows. In order to assess the performance of DTCSM in more complex flows, wall-modeled simulations of high Reynolds number flow over a Gaussian bump (Boeing speed bump) exhibiting smooth-body flow separation are performed. Predictions of surface pressure and skin friction, compared against DNS and experimental data, show improved accuracy from DTCSM in comparison to existing static coefficient (Vreman) and dynamic Smagorinsky model. The computational cost of performing LES with this model is up to 15% higher than the dynamic Smagorinsky model.

42 ENGINEERING↗
DOE OSTI2021

Planetary normal mode computation: Parallel algorithms, performance and reproducibility

This report is an extension of work entitled “Computing planetary interior normal modes with a highly parallel polynomial filtering eigensolver.” by Shi et al., [1] originally presented at the SC18 conference. A highly parallel polynomial filtered eigensolver was developed and exploited to calculate the planetary normal modes. The proposed method is ideally suited for computing interior eigenpairs for large-scale eigenvalue problems as it greatly enhances memory and computational efficiency. In this article, the second-order finite element method is used to further improve the accuracy as only the first-order finite element method was deployed in the previous work. The parallel algorithm, its parallel performance up to 20k processors, and the great computational accuracy are illustrated. The reproducibility of the previous work was successfully performed on the Student Cluster Competition at the SC19 conference by several participant teams using a completely different Mars-model dataset on different clusters. Both weak and strong scaling performances of the reproducibility by the participant teams were impressive and encouraging. The analysis and reflection of their results are demonstrated and future direction is discussed.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

MemXCT: Design, Optimization, Scaling, and Reproducibility of X-Ray Tomography Imaging

Here, this work extends our previous research entitled "MemXCT: Memory-centric X-ray CT Reconstruction with Massive Parallelization" that was originally published at SC19 conference (Hidayetoglu et al., 2019) with reproducibility of the computational imaging performance. X-ray computed tomography (XCT) is regularly used at synchrotron light sources to study the internal morphology of materials at high resolution. However, experimental constraints, such as radiation sensitivity, can result in noisy or undersampled measurements. Further, depending on the resolution, sample size and data acquisition rates, the resulting noisy dataset can be in the order of terabytes. Advanced iterative reconstruction techniques can produce high-quality images from noisy measurements, but their computational requirements have made their use an exception rather than the rule. We propose a novel memory-centric approach that avoids redundant computations at the expense of additional memory complexity. We develop a memory-centric iterative reconstruction system, MemXCT, that uses an optimized SpMV implementation with two-level pseudo-Hilbert ordering and multi-stage input buffering. We evaluate MemXCT on various supercomputer architectures involving KNL and GPU. MemXCT can reconstruct a large (11Kx11K) mouse brain tomogram in 10 seconds using 4096 KNL nodes (256K cores). The results presented in our original article at the SC19 were based on large-scale supercomputing resources. The MemXCT application was selected for the Student Cluster Competition (SCC) Reproducibility Challenge and evaluated on a variety of cloud computing resources by universities around the world in the SC20 conference. We summarize the results of the top-ranked SCC Reproducibility Challenge teams and identify the most pertinent measures for ensuring the reproducibility of our experiments in this article.

47 OTHER INSTRUMENTATION↗
DOE OSTI2022

MYC transcription factors coordinate tryptophan–dependent defence responses and compromise seed yield in Arabidopsis

Robust plant immunity negatively affects other fitness traits, including growth and seed production. Jasmonate (JA) confers broad-spectrum protection against plant consumers by stimulating the degradation of JASMONATE ZIM-DOMAIN (JAZ) proteins, which in turn relieves repression on transcription factors (TFs) coincident with reduced growth and fecundity. The molecular mechanisms underlying JA-mediated decreases in fitness remain largely unknown. To assess the contribution of MYC TFs to growth and reproductive fitness at high levels of defence, we mutated three MYC genes in a JAZ-deficient mutant (jazD) of Arabidopsis thaliana that exhibits strong defence and low seed yield. Genetic epistasis analysis showed that de-repression of MYC TFs in jazD not only conferred strong resistance to insect herbivory but also reduced shoot and root growth, fruit size and seed yield. We also provided evidence that the JAZ–MYC module coordinates the supply of tryptophan with the production of indole glucosinolates and the proliferation of endoplasmic reticulum bodies that metabolise glucosinolates through the action of b-glucosidases. Our results establish MYCs as major regulators of growth- and reproductive–defence tradeoffs and further indicate that these factors coordinate tryptophan availability with the production of amino acid-derived defence compounds.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

On the Cutting of Metals: A Mechanics Viewpoint

The mechanics of large-strain deformation in cutting of metals is discussed, primarily from viewpoint of recent developments in in situ analysis of plastic flow and microstructure characterization. It is shown that a broad range of deformation parameters can be accessed in chip formation—strains of 1–10, strain rates of 10–10 5 /s, and temperatures up to 0.7T m —and controlled. This range is far wider than achievable by any other single-stage, severe plastic deformation (SPD) process. The resulting extreme deformation conditions produce a rich variety of microstructures in the chip. Four principal types of chip formation—continuous, shear-localized, segmented, and mushroom-type—as elucidated first by Nakayama (1974, “The Formation of ‘Saw-Toothed Chip’ in Metal Cutting,” Proceedings of International Conference on Production Engineering, Tokyo, pp. 572–577) are utilized to emphasize the diverse plastic flow phenomena, especially unsteady deformation modes that prevail in cutting. These chip types are intimately connected with the underlying flow, each arising from a distinct mode and triggered by an instability phenomenon. The role of plastic flow instabilities such as shear banding, buckling, and fracture in mediating unsteady flow modes is expounded, along with consequences of the flow modes and chip types for the cutting. Sinuous flow is shown to be the reason why gummy (highly strain-hardening) metals, although relatively soft, are so difficult to cut. Synthesizing the various observations, a hypothesis is put forth that it is the stability of flow modes that determines the mechanics of cutting. Furthermore, this leads to a flow-stability phase diagram that could provide a framework for predicting chip types and process attributes.

42 ENGINEERING↗
DOE OSTI2024

Bactericidal effectors of the Stenotrophomonas maltophilia type IV secretion system: functional definition of the nuclease TfdA and structural determination of TfcB

ABSTRACT Stenotrophomonas maltophilia expresses a type IV protein secretion system (T4SS) that promotes contact-dependent killing of other bacteria and does so partly by secreting the effector TfcB. Here, we report the structure of TfcB, comprising an N-terminal domain similar to the catalytic domain of glycosyl hydrolase (GH-19) chitinases and a C-terminal domain for recognition and translocation by the T4SS. Utilizing a two-hybrid assay to measure effector interactions with the T4SS coupling protein VirD4, we documented the existence of five more T4SS substrates. One of these was protein 20845, an annotated nuclease. A S. maltophilia mutant lacking the gene for 20845 was impaired for killing Escherichia coli , Klebsiella pneumoniae , and Pseudomonas aeruginosa . Moreover, the cloned 20845 gene conferred robust toxicity, with the recombinant E. coli being rescued when 20845 was co-expressed with its cognate immunity protein. The 20845 effector was an 899 amino-acid protein, comprised of a GHH-nuclease domain in its N-terminus, a large central region of indeterminant function, and a C-terminus for secretion. Engineered variants of the 20845 gene that had mutations in the predicted catalytic site did not impede E. coli , indicating that the antibacterial effect of 20845 involves its nuclease activity. Using flow cytometry with DNA staining, we determined that 20845, but not its mutant variants, confers a loss in DNA content of target bacteria. Database searches revealed that uncharacterized homologs of 20845 occur within a range of bacteria. These data indicate that the S. maltophilia T4SS promotes interbacterial competition through the action of multiple toxic effectors, including a potent, novel DNase. IMPORTANCE Stenotrophomonas maltophilia is a multi-drug-resistant, Gram-negative bacterium that is an emerging pathogen of humans. Patients with cystic fibrosis are particularly susceptible to S. maltophilia infection. In hospital water systems and various types of infections, S. maltophilia co-exists with other bacteria, including other pathogens such as Pseudomonas aeruginosa . We previously demonstrated that S. maltophilia has a functional VirB/D4 type VI protein secretion system (T4SS) that promotes contact-dependent killing of other bacteria. Since most work on antibacterial systems involves the type VI secretion system, this observation remains noteworthy. Moreover, S. maltophilia currently stands alone as a model for a human pathogen expressing an antibacterial T4SS. Using biochemical, genetic, and cell biological approaches, we now report both the discovery of a novel antibacterial nuclease (TfdA) and the first structural determination of a bactericidal T4SS effector (TfcB).

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Symbiont-Mediated Protection of Acromyrmex Leaf-Cutter Ants from the Entomopathogenic Fungus Metarhizium anisopliae

Many fungus-growing ants engage in a defensive symbiosis with antibiotic-producing ectosymbiotic bacteria in the genus Pseudonocardia, which help protect the ants’ fungal mutualist from a specialized mycoparasite, Escovopsis. Here, using germfree ant rearing and experimental pathogen infection treatments, we evaluate if Acromyrmex ants derive higher immunity to the entomopathogenic fungus Metarhizium anisopliae from their Pseudonocardia symbionts. We further examine the ecological dynamics and defensive capacities of Pseudonocardia against M. anisopliae across seven different Acromyrmex species by controlling Pseudonocardia acquisition using ant-nonnative Pseudonocardia switches, in vitro challenges, and in situ mass spectrometry imaging (MSI). We show that Pseudonocardia protects the ants against M. anisopliae across different Acromyrmex species and appears to afford higher protection than metapleural gland (MG) secretions. Although Acromyrmex echinatior ants with nonnative Pseudonocardia symbionts receive protection from M. anisopliae regardless of the strain acquired compared with Pseudonocardia-free conditions, we find significant variation in the degree of protection conferred by different Pseudonocardia strains. Additionally, when ants were reared in Pseudonocardia-free conditions, some species exhibit more susceptibility to M. anisopliae than others, indicating that some ant species depend more on defensive symbionts than others. In vitro challenge experiments indicate that Pseudonocardia reduces Metarhizium conidiospore germination area. Our chemometric analysis using matrix-assisted laser desorption/ionization mass spectrometry imaging (MALDI-MSI) reveals that Pseudonocardia-carrying ants produce more chemical signals than Pseudonocardia-free treatments, indicating that Pseudonocardia produces bioactive metabolites on the Acromyrmex cuticle. Our results indicate that Pseudonocardia can serve as a dual-purpose defensive symbiont, conferring increased immunity for both the obligate fungal mutualist and the ants themselves.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

Structural basis for varying drug resistance of SARS-CoV-2 M pro E166 variants

ABSTRACT Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) main protease (M pro ) has an essential role in the virus lifecycle and, accordingly, it is a target for antiviral drugs. Multiple studies have identified an M pro mutation (E166V) that confers strong resistance to clinically relevant inhibitors, including nirmatrelvir, but the underlying mechanism is not fully understood. Here, we report on crystal structures of SARS-CoV-2 M pro E166V in complex with nirmatrelvir, ensitrelvir, and bofutrelvir. The structures suggest that resistance is caused in part by the loss of a direct hydrogen bond and also, especially for nirmatrelvir, by a steric clash with the substituted valine residue. In comparison, the binding of bofutrelvir shows greater flexibility, which may help alleviate this steric effect and allow bofutrelvir to fit the mutant active site despite the loss of a direct polar contact. Thermal stability analyses also corroborate E166V most severely affecting the binding of nirmatrelvir and, to lesser and different extents, ensitrelvir and bofutrelvir. We further show that E166I causes even more severe nirmatrelvir resistance, whereas E166A and E166L have much milder effects. These studies shed light on the molecular mechanisms of a key M pro drug resistance mutation and may help inform the design of next-generation inhibitors. IMPORTANCE Using a combination of high-resolution X-ray crystallographic and biochemical analyses, we reveal the molecular mechanisms by which a mutation in the severe acute respiratory syndrome coronavirus 2 main protease (M pro ) confers strong resistance against clinically relevant antiviral drugs that inhibit M pro activity. The results presented here may help inform the design of next-generation inhibitors to combat the problem of therapy resistance.

Microbiology↗
DOE OSTI2020

Mycobacterium Phage Butters-Encoded Proteins Contribute to Host Defense against Viral Attack [plus supplemental information]

A diverse set of prophage-mediated mechanisms protecting bacterial hosts from infection has been recently uncovered within cluster N mycobacteriophages isolated on the host, Mycobacterium smegmatis mc 2 155. In that context, we unveil a novel defense mechanism in cluster N prophage Butters. By using bioinformatics analyses, phage plating efficiency experiments, microscopy, and immunoprecipitation assays, we show that Butters genes located in the central region of the genome play a key role in the defense against heterotypic viral attack. Our study suggests that a two-component system, articulated by interactions between protein products of genes 30 and 31, confers defense against heterotypic phage infection by PurpleHaze (cluster A/subcluster A3) or Alma (cluster A/subcluster A9) but is insufficient to confer defense against attack by the heterotypic phage Island3 (cluster I/subcluster I1). Therefore, based on heterotypic phage plating efficiencies on the Butters lysogen, additional prophage genes required for defense are implicated and further show specificity of prophage-encoded defense systems. IMPORTANCE: Many sequenced bacterial genomes, including those of pathogenic bacteria, contain prophages. Some prophages encode defense systems that protect their bacterial host against heterotypic viral attack. Understanding the mechanisms undergirding these defense systems is crucial to appreciate the scope of bacterial immunity against viral infections and will be critical for better implementation of phage therapy that would require evasion of these defenses. Furthermore, such knowledge of prophage-encoded defense mechanisms may be useful for developing novel genetic tools for engineering phage-resistant bacteria of industrial importance.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Fragile Earth: AI for Climate Mitigation, Adaptation, and Environmental Justice

The Fragile EarthWorkshop is a recurring event that gathers the research community to find and explore howdata science can measure and progress climate and social issues, following the framework of the United Nations Sustainable Development Goals (SDGs).Fragile Earth 2022: AI for Climate Mitigation, Adaptation, and Environmental Justice is a workshop taking place as part of the ACM's KDD 2022 Conference on research in knowledge discovery and data mining and their applications. The dates for the Conference are August 14-18, 2022.

Abe, Naoki↗
DOE OSTI2024

Pharmacologic or genetic interference with atrogene signaling protects against glucocorticoid-induced musculoskeletal and cardiac disease

Despite their beneficial actions as immunosuppressants, glucocorticoids (GC) have devastating effects on the musculoskeletal and cardiac systems, as long-term treated patients exhibit high incidence of falls, bone fractures, and cardiovascular events. Herein, we show that GC upregulate simultaneously in bone, skeletal muscle, and the heart the expression of E3 ubiquitin ligases (atrogenes), known to stimulate the proteasomal degradation of proteins. Activation of vitamin D receptor (VDR) signaling with the VDR ligands calcitriol or eldecalcitol prevented GC-induced atrogene upregulation in vivo and ex vivo in bone/muscle organ cultures and preserved tissue structure/mass and function of the 3 tissues in vivo. Direct pharmacologic inhibition of the proteasome with carfilzomib also conferred musculoskeletal protection. Genetic loss of the atrogene MuRF1-mediated protein ubiquitination in ΔRING mice afforded temporary or sustained protection from GC excess in bone or skeletal and heart muscle. We concluded that the atrogene pathway downstream of MuRF1 underlies GC action in bone, muscle, and the heart, and it can be pharmacologically or genetically targeted to confer protection against the damaging actions of GC simultaneously in the 3 tissues.

Research & Experimental Medicine↗
DOE OSTI2022

A numerical approach for detecting switch-like bistability in mass action chemical reaction networks with conservation laws

Theoretical analysis of signaling pathways can provide a substantial amount of insight into their function. One particular area of research considers signaling pathways capable of assuming two or more stable states given the same amount of signaling ligand. This phenomenon of bistability can give rise to switch-like behavior, a mechanism that governs cellular decision making. Investigation of whether or not a signaling pathway can confer bistability and switch-like behavior, without knowledge of specific kinetic rate constant values, is a mathematically challenging problem. Recently a technique based on optimization has been introduced, which is capable of finding example parameter values that confer switch-like behavior for a given pathway. Although this approach has made it possible to analyze moderately sized pathways, it is limited to reaction networks that presume a uniterminal structure. It is this limited structure we address by developing a general technique that applies to any mass action reaction network with conservation laws. In this paper we developed a generalized method for detecting switch-like bistable behavior in any mass action reaction network with conservation laws. The method involves (1) construction of a constrained optimization problem using the determinant of the Jacobian of the underlying rate equations, (2) minimization of the objective function to search for conditions resulting in a zero eigenvalue, (3) computation of a confidence level that describes if the global minimum has been found and (4) evaluation of optimization values, using either numerical continuation or directly simulating the ODE system, to verify that a bistability region exists. The generalized method has been tested on three motifs known to be capable of bistability. We have developed a variation of an optimization-based method for the discovery of bistability, which is not limited to uniterminal chemical reaction networks. Successful completion of the method provides an S-shaped bifurcation diagram, which indicates that the network acts as a bistable switch for the given optimization parameters.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗