SEARCH · Engineering Papers

Results for “respiratory complex I”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 91 records · Page 5

All sources NASA NTRS DOE OSTI

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 OSTI2022

Resolving Lipid Structures for Threat-Agnostic Signatures

Development of threat-agnostic biosignatures is critical for the detection of and response to biological threats that go beyond the historical list-based approach. Lipids are a class of structurally diverse biomolecules that hold great promise as relevant threat-agnostic biosignatures, yet the fundamental challenge of resolving and decoding structural differences among lipids remains an obstacle in the pursuit of signature development. Structural analysis of complex lipids remains stymied due to spectral complexity and the inability to resolve low abundance lipids. To address these challenges, we integrated ozonolysis with ion mobility mass spectrometry (IMS-MS) analysis to determine the structural information for complex lipid species and (2) reduced lipid sample complexity prior to introduction into the IMS-MS using normal phase high pressure liquid chromatography (HPLC) to enable the ozonolysis of low abundance lipids. We applied these new capabilities to define the lipid structures in antimicrobial resistant (AMR) and antimicrobial susceptible (AMS) biothreat organism Yersinia pestis and cells infected with Severe Acute Respiratory Syndrome coronavirus 2 (SARS-CoV-2) variants.

59 BASIC BIOLOGICAL SCIENCES↗

NASA NTRS2017

NASA Ames Institutional Scientific Collection (ISC)

NASA's current human space flight research is directed towards enabling human space exploration beyond Low Earth Orbit (LEO). The Space Flight Payload Projects; Rodent Research, Cell Science, and Microbial Labs, flown on the International Space Station (ISS), benefit both the global life sciences and commercial space communities. Verified data sets, science results, peer-reviewed publications, and returned biospecimens, collected and analyzed for flight and ground investigations, are all part of the knowledge base within NASA’s Human Exploration and Operations Mission Directorate’s Space Life and Physical Sciences Research and Applications (SLPSRA) Division, specifically the Human Research and Space Biology Programs. These data and biospecimens are made available through the public LSDA website. The Ames Institutional Scientific Collection (ISC), or ARC Biobank, stores flight and ground biospecimens from Space Shuttle and ISS programs. These specimens are curated and managed by the Ames Life Sciences Data Archive (ALSDA), an internal node of NASA's Life Sciences Data Archive (LSDA). The ARC Biolbank stores over 15,000 specimens from experiments dating from 1984 to present. Currently available specimens include tissues from the circulatory, digestive, endocrine, excretory, integumentary, muscular, neurosensory, reproductive, respiratory and skeletal systems. The most recent contributions include RNA, DNA and protein extracts from Rodent Research 1 and tissues from Rodent Research 4. NASA's biospecimen collection represents a unique and limited resource. The use of these biospecimens maximizes utilization and scientific return from these unique spaceflight payload and ground control research subjects. These biospecimens are harvested following complex, costly NASA research activities to meet primary scientific objectives. Once the primary scientific objectives have been met, the remaining specimens are made available to provide secondary opportunities for complementary studies or new investigations to broaden research without large expenditures of time or resources. Innovative ways of sharing this information ultimately advances the frontiers of human space exploration as well as scientific understanding of the effects of gravity on life on earth.

French, Alison J.↗

DOE OSTI2020

Structural basis of RNA cap modification by SARS-CoV-2

The severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2), the causative agent of COVID-19 illness, has caused millions of infections worldwide. In SARS coronaviruses, the non-structural protein 16 (nsp16), in conjunction with nsp10, methylates the 5'-end of virally encoded mRNAs to mimic cellular mRNAs, thus protecting the virus from host innate immune restriction. We report here the high-resolution structure of a ternary complex of SARS-CoV-2 nsp16 and nsp10 in the presence of cognate RNA substrate analogue and methyl donor, S-adenosyl methionine (SAM). The nsp16/nsp10 heterodimer is captured in the act of 2'-O methylation of the ribose sugar of the first nucleotide of SARS-CoV-2 mRNA. We observe large conformational changes associated with substrate binding as the enzyme transitions from a binary to a ternary state. This induced fit model provides mechanistic insights into the 2'-O methylation of the viral mRNA cap. We also discover a distant (25 Å) ligand-binding site unique to SARS-CoV-2, which can alternatively be targeted, in addition to RNA cap and SAM pockets, for antiviral development.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2023

Respiration is essential for aerobic growth of Zymomonas mobilis ZM4

ABSTRACT Zymomonas mobilis is an alpha-proteobacterium that is a promising platform for industrial scale production of biofuels due to its efficient ethanol fermentation and low biomass generation. Z. mobilis is aerotolerant and encodes a complete respiratory electron transport chain, but the benefit of respiration for growth in oxic conditions has never been confirmed, despite decades of research. Growth and ethanol production of wild-type Z. mobilis is poor in oxic conditions indicating that it does not benefit from oxidative phosphorylation. Additionally, in previous studies, aerobic growth improved significantly when respiratory genes were disrupted ( ndh ) or acquired point mutations ( cydA and cydB ), even if respiration was significantly reduced by these changes. Here, we obtained clean deletions of respiratory genes ndh and cydAB , individually and in combination, and showed, for the first time, that deletion of cydAB completely inhibited O 2 respiration and dramatically reduced growth in oxic conditions. Both respiration and aerobic growth were restored by expressing a heterologous, water-forming NADH oxidase, noxE . Oxygen can have many negative effects, including formation of reactive oxygen species (ROS) or directly inactivating oxygen sensitive enzymes. Our results suggest that the effect of molecular oxygen on enzymes had a greater negative impact on Z. mobilis than formation of ROS. This result shows that the main role of the electron transport chain in Z. mobilis is reducing the intracellular concentration of molecular oxygen, helping to explain why it is beneficial for Z. mobilis to use electron transport chain complexes that have little capacity to contribute to oxidative phosphorylation. IMPORTANCE A key to producing next-generation biofuels is to engineer microbes that efficiently convert non-food materials into drop-in fuels, and to engineer microbes effectively, we must understand their metabolism thoroughly. Zymomonas mobilis is a bacterium that is a promising candidate biofuel producer, but its metabolism remains poorly understood, especially its metabolism when exposed to oxygen. Although Z. mobilis respires with oxygen, its aerobic growth is poor, and disruption of genes related to respiration counterintuitively improves aerobic growth. This unusual result has sparked decades of research and debate regarding the function of respiration in Z. mobilis . Here, we used a new set of mutants to determine that respiration is essential for aerobic growth and likely protects the cells from damage caused by oxygen. We conclude that the respiratory pathway of Z. mobilis should not be deleted from chassis strains for industrial production because this would yield a strain that is intolerant of oxygen, which is more difficult to manage in industrial settings.

09 BIOMASS FUELS↗

DOE OSTI2021

2'-O methylation of RNA cap in SARS-CoV-2 captured by serial crystallography

The genome of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) coronavirus has a capping modification at the 5'-untranslated region (UTR) to prevent its degradation by host nucleases. These modifications are performed by the Nsp10/14 and Nsp10/16 heterodimers using S-adenosylmethionine as the methyl donor. Nsp10/16 heterodimer is responsible for the methylation at the ribose 2'-O position of the first nucleotide. To investigate the conformational changes of the complex during 2'-O methyltransferase activity, we used a fixed-target serial synchrotron crystallography method at room temperature. We determined crystal structures of Nsp10/16 with substrates and products that revealed the states before and after methylation, occurring within the crystals during the experiments. Here we report the crystal structure of Nsp10/16 in complex with Cap-1 analog ( m7 GpppA m2'-O ). Inhibition of Nsp16 activity may reduce viral proliferation, making this protein an attractive drug target.

59 BASIC BIOLOGICAL SCIENCES↗

NASA NTRS1983

Human respiration at rest in rapid compression and at high pressures and gas densities

The ventilation (V), end-tidal PCO2 (PACO2), and CO2 elimination rate were determined in men at rest breathing CO2-free gas over the pressure range 1-50 ATA and the gas density range 0.4-25 g/l, during slow and rapid compressions, at stable elevated ambient pressures and during slow decompressions. Progressive increase in pulmonary gas flow resistance due to elevation of ambient pressure and inspired gas density to the He-O2 equivalent of 5000 feet of seawater was found to produce a complex pattern of change in PACO2. It was found that as both ambient pressure and pulmonary gas flow resistance were progressively raised, PACO2 at first increased, went through a maximum, and then declined towards values near the 1 ATA level. It is concluded that this pattern of PACO2 change results from the interaction on ventilation of the increase in pulmonary resistance due to the elevation of gas density with the increase in respiratory drive postulated as due to generalized central nervous system excitation associated with exposure to high hydrostatic pressure. It is suggested that a similar interaction exists between increased gas flow resistance and the increase in respiratory drive related to nitrogen partial pressure and the resulting narcosis.

Gelfand, R.↗

DOE OSTI2025

Structures of respiratory syncytial virus G bound to broadly reactive antibodies provide insights into vaccine design

Abstract Respiratory syncytial virus (RSV) is a leading cause of severe lower respiratory tract disease in infants and older adults. The attachment glycoprotein (RSV G) binds to the chemokine receptor CX3CR1 to promote viral entry and modulate host immunity. Antibodies against RSV G are a known correlate of protection. Previously, several broadly reactive, high-affinity anti-RSV G human monoclonal antibodies were isolated from RSV-exposed individuals and were shown to be protective in vitro and in vivo. Here, we determined the structures of three of these antibodies in complex with RSV G and defined distinct conformational epitopes comprised of highly conserved RSV G residues. Binding competition and structural studies demonstrated that this highly conserved region displays two non-overlapping antigenic sites. Analyses of anti-RSV G antibody sequences reveal that antigenic site flexibility may promote the elicitation of diverse antibody germlines. Together, these findings provide a foundation for next-generation RSV prophylactics, and they expand concepts in vaccine design for the elicitation of germline lineage-diverse, broadly reactive, high-affinity antibodies.

Science & Technology - Other Topics↗

NASA NTRS1979

Early detection of disease: The correlation of the volatile organic profiles from patients with upper respiratory infections with subjects of normal profiles

A method is described whereby a transevaporator is used for sampling 60-100 microns of aqueous sample. Volatiles are stripped from the sample either by a stream of helium and collection on a porous polymer, Tenax, or by 0.8 ml of 2-chloropropane and collected on glass beads. The volatiles are thermally desorbed into a precolumn which is connected to a capillary gas chromatographic column for analysis. The technique is shown to be reproducible and suitable for determining chromatographic profiles for a wide variety of sample types. Using a transevaporator sampling technique, the volatile profiles from 70 microns of serum were obtained by capillary column gas chromatography. The complex chromatograms were interpreted by a combination of manual and computer techniques and a two peak ratio method devised for the classification of normal and virus infected sera. Using the K-Nearest Neighbor approach, 85.7 percent of the unknown samples were classified correctly. Some preliminary results indicate the possible use of the method for the assessment of virus susceptibility.

Zlatkis, A.↗

DOE OSTI2021

Induction of Protection in Mice against a Chlamydia muridarum Respiratory Challenge by a Vaccine Formulated with the Major Outer Membrane Protein in Nanolipoprotein Particles

Chlamydia trachomatis is a sexually transmitted bacterium that infects over 130 million individuals worldwide annually. To implement a vaccine, we developed a cell-free co-translational system to express the Chlamydia muridarum major outer membrane protein (MOMP). This approach uses a nanolipoprotein particles (tNLP) made from ApoA1 protein, amphiphilic telodendrimer and lipids that self-assemble to form 10–25 nm discs. These tNLP provide a protein-encapsulated lipid support to solubilize and fold membrane proteins. The cell-free system co-translated MOMP and ApoA1 in the presence of telodendrimer mixed with lipids. The MOMP-tNLP complex was amenable to CpG and FSL-1 adjuvant addition. To investigate the ability of MOMP-tNLP+CpG+FSL-1 to induce protection against an intranasal (i.n.) C. muridarum challenge, female mice were vaccinated intramuscularly (i.m.) or i.n. and i.m. simultaneously 4 weeks apart. Following vaccination with MOMP-tNLP+CpG+FSL-1, mice mounted significant humoral and cell-mediated immune responses. Following the i.n. challenge, mice vaccinated with MOMP-tNLP+CpG+FSL-1 i.n. + i.m. group were protected as determined by the percentage change in body weight and by the number of C. muridarum inclusion forming units (IFU) recovered from the lungs. To our knowledge, this is the first time a MOMP-based vaccine formulated in tNLP has been shown to protect against C. muridarum.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2021

Structural, Electronic, and Electrostatic Determinants for Inhibitor Binding to Subsites S1 and S2 in SARS-CoV-2 Main Protease

Creating small-molecule antivirals specific for severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) proteins is crucial to battle coronavirus disease 2019 (COVID-19). SARS-CoV-2 main protease (Mpro) is an established drug target for the design of protease inhibitors. We performed a structure–activity relationship (SAR) study of noncovalent compounds that bind in the enzyme’s substrate-binding subsites S1 and S2, revealing structural, electronic, and electrostatic determinants of these sites. The study was guided by the X-ray/neutron structure of Mpro complexed with Mcule-5948770040 (compound 1), in which protonation states were directly visualized. Virtual reality-assisted structure analysis and small-molecule building were employed to generate analogues of 1. In vitro enzyme inhibition assays and room-temperature X-ray structures demonstrated the effect of chemical modifications on Mpro inhibition, showing that (1) maintaining correct geometry of an inhibitor’s P1 group is essential to preserve the hydrogen bond with the protonated His163; (2) a positively charged linker is preferred; and (3) subsite S2 prefers nonbulky modestly electronegative groups.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2022

Molecular insights into the differential dynamics of SARS-CoV-2 variants of concern

Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) has affected the lives and livelihood of millions of individuals around the world. It has mutated several times after its first inception, with an estimated two mutations occurring every month. Although we have been successful in developing vaccines against the virus, the emergence of variants has enabled it to escape therapy. Few of the generated variants are also reported to be more infectious than the wild-type (WT). In this study, we analyze the attributes of all RBD/ACE2 complexes for the reported VOCs, namely, Alpha, Beta, Gamma, and Delta through computer simulations. Results indicate differences in orientation and binding energies of the VOCs from the WT. Overall, it was observed that electrostatic interactions play a major role in the binding of the complexes. Detailed residue level energetics revealed that the most prominent changes in interaction energies were seen particularly at the mutated residues which were present at RBD/ACE2 interface. We found that the Delta variant is one of the most tightly bound variants of SARS-CoV-2 with dynamics similar to WT. The high binding affinity of RBD towards ACE2 is indicative of an increase in viral transmission and infectivity. The details presented in our study provide additional information for the design and development of effective therapeutic strategies for the emerging variants of the virus in the future.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2025

Structural basis of heme scavenging by the ChtA and HtaA hemophores in Corynebacterium diphtheriae

Corynebacterium diphtheriae causes diphtheria, a potentially fatal infectious disease that damages tissues in the upper respiratory tract. In order to proliferate, this pathogen acquires the essential nutrient iron from heme (iron-protoporphyrin IX) primarily found in human hemoglobin (Hb). C. diphtheriae secretes ChtA and HtaA hemophore proteins that bind ferric heme (hemin) via conserved region (CR) domains. Here, we demonstrate that their CR domains scavenge hemin after it is spontaneously released from Hb, and define the structural basis of hemin binding to ChtA and the N-terminal CR domain from HtaA by determining X-ray crystal structures of their protein-hemin complexes. Resonance Raman and electron paramagnetic resonance experiments demonstrate that the CR domains from ChtA and HtaA engage in pentacoordinate hemin binding through a conserved iron-tyrosyl linkage, though variations in their hemin pockets alter the way they stabilize the axial tyrosine and mask hemin’s metal. The importance of these interactions is probed using isothermal titration calorimetry experiments, which represent the first quantitative assessment of CR-hemin affinity and reveal that ChtA binds hemin via an enthalpically driven process. Hemin partitioning experiments using native mass spectrometry demonstrate that the cohort of CR domains within C. diphtheriae ’s hemin-uptake system have dissociation constants for hemin between 0.8 and 22 nM, raising the possibility that affinity differences contribute to the directional flow of hemin into the cell. Collectively, the results of this work provide insight into how C. diphtheriae and other pathogenic and commensal corynebacterium species utilize CR domains to scavenge iron rich hemin from their environment.

Corynebacterium diphtheriae↗

DOE OSTI2020

EBSD seed LDRD project: Does Corona Virus – 2019 (COVID-19) and Seasonal Flu have similar meteorology and air quality controls driving their spread?

Seasonal influenza and Influenza like Illnesses (ILI) pose a serious public health risk and in turn affect the economy. Various factors affect ILI cases and mortality, including the pathogen and its interaction with the host, as well as environmental and socioeconomic factors such as meteorology, household structure, air pollution, urbanization, and population. Despite the growing number of studies on influenza and ILI, challenges remain in forecasting the timing of seasonal onset, outbreak patterns, and key factors affecting transmission. In particular, the impacts of meteorology and air quality on ILI have been challenging to understand, with linear regression studies focused on different geographic regions producing contradictory results. For example, influenza seasonality has been associated with cold-dry conditions in temperate mid-latitudes, but with humid-rainy conditions in tropical climates. These apparently contradictory results imply that the relationships between influenza cases and atmospheric variables may be too complex to be captured by linear regression models. In this seed project, we analyze meteorology and air quality variables from numerical models to determine which atmospheric variables are most helpful in predicting weekly changes in recorded flu cases. In contrast to most previous studies that relied on linear regression analysis to predict the timing of the flu onset or peak, we employ a robust machine learning algorithm to evaluate the contribution of atmospheric variables to weekly changes in recorded ILI cases. These results may also be relevant to the spread of other respiratory illnesses such as Corona Virus Disease – 2019 (COVID-19).

60 APPLIED LIFE SCIENCES↗

NASA NTRS2010

Filovirus Research: The Need for an Integrated Approach in Time and Space

The episodic appearance of Ebola virus (EBOV) and Marburg virus (MARV) across central Africa over the last 15 years not only underscores the importance of filoviruses as uniquely virulent agents to both human and wildlife communities but also implies a very complex transmission scenario that must be understood if we are to prevent or mitigate filovirus outbreaks in the future. Efforts of a global network of scientists and healthcare workers have expanded our knowledge of filoviruses to meet the growing threat of Ebola and Marburg hemorrhagic fevers in Africa. In recent decades, several newly emerging diseases have resulted in major threats to both affected communities and global public health. Viruses from wildlife hosts in particular, have exhibited a capability for cross-species transmission (CST), and have caused high-impact diseases in humans Such as Ebola and Marburg hemorrhagic fevers, Nipah and severe acute respiratory syndrome (SAILS). It has been estimated that about 60.3% (Jones et al. 2008) of human infectious diseases are of animal origin (zoonoses) and even some important viral diseases that are traditionally considered of human origin, for example measles and smallpox, may very well have their prehistoric origins in wildlife (Wolfe et al 2007). It maybe logical and prudent therefore, to anticipate that there are other, new filoviruses out there that will cross into humans at some point in time. If we anticipate that these will happen and wish to be prepared for and mitigate this potential, then an understanding of filoviruses as a biologic system in the environment will be essential to that process. We will need to know how the ecological dynamic of CST interacts with a 'new' viruse's evolutionary factors to overcome environmental, demographic and host-specific barriers to transmission and infectivity to humans.

Pinzon, Jorge E.↗

DOE OSTI2021

Exploring the Role of Glycans in the Interaction of SARS-CoV-2 RBD and Human Receptor ACE2

COVID-19 is a highly infectious respiratory disease caused by the novel coronavirus SARS-CoV-2. It has become a global pandemic and its frequent mutations may pose new challenges for vaccine design. During viral infection, the Spike RBD of SARS-CoV-2 binds the human host cell receptor ACE2, enabling the virus to enter the host cell. Both the Spike and ACE2 are densely glycosylated, and it is unclear how distinctive glycan types may modulate the interaction of RBD and ACE2. Detailed understanding of these determinants is key for the development of novel therapeutic strategies. To this end, we perform extensive all-atom simulations of the (i) RBD-ACE2 complex without glycans, (ii) RBD-ACE2 with oligomannose MAN9 glycans in ACE2, and (iii) RBD-ACE2 with complex FA2 glycans in ACE2. These simulations identify the key residues at the RBD-ACE2 interface that form contacts with higher probabilities, thus providing a quantitative evaluation that complements recent structural studies. Notably, we find that this RBD-ACE2 contact signature is not altered by the presence of different glycoforms, suggesting that RBD-ACE2 interaction is robust. Applying our simulated results, we illustrate how the recently prevalent N501Y mutation may alter specific interactions with host ACE2 that facilitate the virus-host binding. Furthermore, our simulations reveal how the glycan on Asn90 of ACE2 can play a distinct role in the binding and unbinding of RBD. Finally, an energetics analysis shows that MAN9 glycans on ACE2 decrease RBD-ACE2 affinity, while FA2 glycans lead to enhanced binding of the complex. Together, our results provide a more comprehensive picture of the detailed interplay between virus and human receptor, which is much needed for the discovery of effective treatments that aim at modulating the physical-chemical properties of this virus.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2025

Hardware-in-Loop Modules for Testing Automated Ventilator Controllers

Automated ventilator controllers have the potential to simplify oxygen and carbon dioxide management for trauma. In the pre-hospital or military medicine environment, trauma care can be required for prolonged periods by personnel with limited ventilator management training. As such, there is a need for closed-loop control systems that can adapt ventilator management to a complex, ever-changing medical environment. Here, we present a novel hardware-in-loop test platform for the independent troubleshooting and evaluation of oxygen and carbon dioxide automated ventilator management capabilities. The oxygen management system provides an analogue blood oxygen signal that is responsive to the fraction of inspired oxygen and the peak inspiratory pressure ventilator settings. A tested oxygenation controller successfully reached the target oxygen saturation within 5 min. The carbon dioxide removal system integrates with commercial ventilator technology and mimics carbon dioxide generation, lung compliance, and airway resistance while providing an end-tidal carbon dioxide level that is responsive to changes in the tidal volume and respiratory rate settings. A test mechanical ventilator controller was able to regulate EtCO2 regardless of the starting value within 10 min. This highlights the system’s functionality and provides proof-of-concept demonstrations for how the hardware-in-loop test platforms can be used for evaluating closed-loop controller technologies.

Berard, David (ORCID:0000000322863846)↗

DOE OSTI2023

Mathematical solutions in internal dose assessment: A comparison of Python-based differential equation solvers in biokinetic modeling

Abstract In biokinetic modeling systems employed for radiation protection, biological retention and excretion have been modeled as a series of discretized compartments representing the organs and tissues of the human body. Fractional retention and excretion in these organ and tissue systems have been mathematically governed by a series of coupled first-order ordinary differential equations (ODEs). The coupled ODE systems comprising the biokinetic models are usually stiff due to the severe difference between rapid and slow transfers between compartments. In this study, the capabilities of solving a complex coupled system of ODEs for biokinetic modeling were evaluated by comparing different Python programming language solvers and solving methods with the motivation of establishing a framework that enables multi-level analysis. The stability of the solvers was analyzed to select the best performers for solving the biokinetic problems. A Python-based linear algebraic method was also explored to examine how the numerical methods deviated from an analytical or semi-analytical method. Results demonstrated that customized implicit methods resulted in an enhanced stable solution for the inhaled 60 Co (Type M) and 131 I (Type F) exposure scenarios for the inhalation pathway of the International Commission on Radiological Protection (ICRP) Publication 130 Human Respiratory Tract Model (HRTM). The customized implementation of the Python-based implicit solvers resulted in approximately consistent solutions with the Python-based matrix exponential method ( expm ). The differences generally observed between the implicit solvers and expm are attributable to numerical precision and the order of numerical approximation of the numerical solvers. This study provides the first analysis of a list of Python ODE solvers and methods by comparing their usage for solving biokinetic models using the ICRP Publication 130 HRTM and provides a framework for the selection of the most appropriate ODE solvers and methods in Python language to implement for modeling the distribution of internal radioactivity.

61 RADIATION PROTECTION AND DOSIMETRY↗