Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “biophysical models”

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 577 records · Page 32

Computational design of non-porous pH-responsive antibody nanoparticles

Abstract Programming protein nanomaterials to respond to changes in environmental conditions is a current challenge for protein design and is important for targeted delivery of biologics. Here we describe the design of octahedral non-porous nanoparticles with a targeting antibody on the two-fold symmetry axis, a designed trimer programmed to disassemble below a tunable pH transition point on the three-fold axis, and a designed tetramer on the four-fold symmetry axis. Designed non-covalent interfaces guide cooperative nanoparticle assembly from independently purified components, and a cryo-EM density map closely matches the computational design model. The designed nanoparticles can package protein and nucleic acid payloads, are endocytosed following antibody-mediated targeting of cell surface receptors, and undergo tunable pH-dependent disassembly at pH values ranging between 5.9 and 6.7. The ability to incorporate almost any antibody into a non-porous pH-dependent nanoparticle opens up new routes to antibody-directed targeted delivery.

Biochemistry & Molecular Biology↗

Towards the Structure Determination of a Modulated Protein Crystal: The Semicrystalline State of Profilin:Actin

One of the remaining challenges to structural biology is the solution of modulated structures. While small molecule crystallographers have championed this type of structure, to date, no modulated macromolecular structures have been determined. Modulation of the molecular structures within the crystal can produce satellite reflections or a superlattice of reflections in reciprocal space. We have developed the data collection methods and strategies that are needed to collect and analyze these data. If the macromolecule's crystal lattice is composed of physiologically relevant packing contacts, structural changes induced under physiological conditions can cause distortion relevant to the function and biophysical processes of the molecule making up the crystal. By careful measurement of the distortion, and the corresponding three-dimensional structure of the distorted molecule, we will visualize the motion and mechanism of the biological macromolecule(s). We have measured the modulated diffraction pattern produced by the semicrystalline state of profilin:actin crystals using highly parallel and highly monochromatic synchrotron radiation coupled with fine phi slicing (0.001-0.010 degrees) for structure determination. These crystals present these crystals present a unique opportunity to address an important question in structural biology. The modulation is believed to be due to the formation of actin helical filaments from the actin beta ribbon upon the pH-induced dissociation of profilin. To date, the filamentous state of actin has resisted crystallization and no detailed structures are available. The semicrystalline state profilin:actin crystals provides a unique opportunity to understand the many conformational states of actin. This knowledge is essential for understanding the dynamics underlying shape changes and motility of eukaryotic cells. Many essential processes, such as cytokinesis, phagocytosis, and cellular migration depend upon the capacity of the actin microfilament system to be restructured in a controlled manner via polymerization, depolymerization, severing, cross-linking, and anchorage. The structure the semicrystalline state of profilin:actin will challenge and validate current models of muscle contraction and cell motility. The methodology and theory under development will be easily extendable to other systems.

Borgstahl, G.↗

Interplay of vibration- and environment-assisted energy transfer

Abstract We study the interplay between two environmental influences on excited state energy transfer in photosynthetic light harvesting complexes, namely, vibrationally assisted energy transfer (VAET) and environment-assisted quantum transport (ENAQT), considering a dimeric chromophore donor–acceptor model as a prototype for larger systems. We demonstrate how the basic features of the excitonic energy transfer are influenced by these two environments, both separately and together, with the environment being fully quantum in the case of VAET and treated in the Haken–Strobl–Reineker classical limit in the case of ENAQT. Our results reveal that in the weak noise regime, the presence of a classical noise source is detrimental to the energy transfer that is resonantly assisted by the exciton-vibration interactions intrinsic to VAET. In the strong noise regime we reproduce all the features of ENAQT, including the turnover into a Zeno regime where energy transfer is suppressed, and VAET is insignificant.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Mechanism of the cadherin–catenin F-actin catch bond interaction

Mechanotransduction at cell–cell adhesions is crucial for the structural integrity, organization, and morphogenesis of epithelia. At cell–cell junctions, ternary E-cadherin/β-catenin/αE-catenin complexes sense and transmit mechanical load by binding to F-actin. The interaction with F-actin, described as a two-state catch bond, is weak in solution but is strengthened by applied force due to force-dependent transitions between weak and strong actin-binding states. Here, we provide direct evidence from optical trapping experiments that the catch bond property principally resides in the αE-catenin actin-binding domain (ABD). Consistent with our previously proposed model, the deletion of the first helix of the five-helix ABD bundle enables stable interactions with F-actin under minimal load that are well described by a single-state slip bond, even when αE-catenin is complexed with β-catenin and E-cadherin. Our data argue for a conserved catch bond mechanism for adhesion proteins with structurally similar ABDs. We also demonstrate that a stably bound ABD strengthens load-dependent binding interactions between a neighboring complex and F-actin, but the presence of the other αE-catenin domains weakens this effect. These results provide mechanistic insight to the cooperative binding of the cadherin–catenin complex to F-actin, which regulate dynamic cytoskeletal linkages in epithelial tissues.

59 BASIC BIOLOGICAL SCIENCES↗

Predicting permeation of compounds across the outer membrane of P. aeruginosa using molecular descriptors

The ability Gram-negative pathogens have at adapting and protecting themselves against antibiotics has increasingly become a public health threat. Data-driven models identifying molecular properties that correlate with outer membrane (OM) permeation and growth inhibition while avoiding efflux could guide the discovery of novel classes of antibiotics. Here we evaluate 174 molecular descriptors in 1260 antimicrobial compounds and study their correlations with antibacterial activity in Gram-negative Pseudomonas aeruginosa. The descriptors are derived from traditional approaches quantifying the compounds’ intrinsic physicochemical properties, together with, bacterium-specific from ensemble docking of compounds targeting specific MexB binding pockets, and all-atom molecular dynamics simulations in different subregions of the OM model. Using these descriptors and the measured inhibitory concentrations, we design a statistical protocol to identify predictors of OM permeation/inhibition. We find consistent rules across most of our data highlighting the role of the interaction between the compounds and the OM. An implementation of the rules uncovered in our study is shown, and it demonstrates the accuracy of our approach in a set of previously unseen compounds. Our analysis sheds new light on the key properties drug candidates need to effectively permeate/inhibit P. aeruginosa, and opens the gate to similar data-driven studies in other Gram-negative pathogens.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Carbon nanotube porin diffusion in mixed composition supported lipid bilayers

Carbon nanotube porins (CNTPs), short pieces of carbon nanotubes capable of self-inserting into a lipid bilayer, represent a simplified model of biological membrane channels. We have used high-speed atomic force microscopy (HS-AFM) and all-atom molecular dynamics (MD) simulations to study the behavior of CNTPs in a mixed lipid membrane consisting of DOPC lipid with a variable percentage of DMPC lipid added to it. HS-AFM data reveal that the CNTPs undergo diffusive motion in the bilayer plane. Motion trajectories extracted from the HS-AFM movies indicate that CNTPs exhibit diffusion coefficient values broadly similar to values reported for membrane proteins in supported lipid bilayers. The data also indicate that increasing the percentage of DMPC leads to a marked slowing of CNTP diffusion. MD simulations reveal a CNTP-lipid assembly that diffuses in the membrane and show trends that are consistent with the experimental observations.

77 NANOSCIENCE AND NANOTECHNOLOGY↗

Distinct Secondary Structures of the Leucine-Rich Repeat Proteoglycans Decorin and Biglycan: Glycosylation-Dependent Conformational Stability

Biglycan and decorin, closely related small leucine-rich repeat proteoglycans, have been overexpressed in eukaryotic cers and two major glycoforms isolated under native conditions: a proteoglycan substituted with glycosaminoglycan chains; and a core protein form secreted devoid of glycosaminoglycans. A comparative biophysical study of these glycoforms has revealed that the overall secondary structures of biglycan and decorin are different. Far-UV Circular Dichroism (CD) spectroscopy of decorin and biglycan proteoglycans indicates that, although they are predominantly Beta-sheet, biglycan has a significantly higher content of alpha-helical structure. Decorin proteoglycan and core protein are very similar, whereas the biglycan core protein exhibits closer similarity to the decorin glycoforms than to. the biglycan proteoglycan form. However, enzymatic removal of the chondroitin sulfate chains from biglycan proteoglycan does not induce a shift to the core protein structure, suggesting that the fmal form is influenced by polysaccharide addition only during biosynthesis. Fluorescence emission spectroscopy demonstrated that the single tryptophan residue, which is at a conserved position at the C-terminal domain of both biglycan and decorin, is found in similar microenvironments. This indicates that at least in this specific domain, the different glycoforms do exhibit apparent conservation of structure. Exposure of decorin and biglycan to 10 M urea resulted in an increase in fluorescent intensity, which indicates that the emission from tryptophan in the native state is quenched. Comparison of urea-induced protein unfolding curves provided further evidence that decorin and biglycan assume different structures in solution. Decorin proteoglycan and core protein unfold in a manner similar to a classic two-state model, in which there is a steep transition to an unfolded state between 1-2 M urea. The biglycan core protein also shows a similar steep transition. However, biglycan proteoglycan shows a broad unfolding transition between 1-6 M urea, probably indicating the presence of stable unfolding intermediates.

Krishnan, Priya↗

The homogenous alternative to biomineralization: Zn- and Mn-rich materials enable sharp organismal “tools” that reduce force requirements

We measured hardness, modulus of elasticity, and, for the first time, loss tangent, energy of fracture, abrasion resistance, and impact resistance of zinc- and manganese-enriched materials from fangs, stings and other “tools” of an ant, spider, scorpion and nereid worm. The mechanical properties of the Zn- and Mn-materials tended to cluster together between plain and biomineralized “tool” materials, with the hardness reaching, and most abrasion resistance values exceeding, those of calcified salmon teeth and crab claws. Atom probe tomography indicated that Zn was distributed homogeneously on a nanometer scale and likely bound as individual atoms to more than ¼ of the protein residues in ant mandibular teeth. This homogeneity appears to enable sharper, more precisely sculpted “tools” than materials with biomineral inclusions do, and also eliminates interfaces with the inclusions that could be susceptible to fracture. Based on contact mechanics and simplified models, we hypothesize that, relative to plain materials, the higher elastic modulus, hardness and abrasion resistance minimize temporary or permanent tool blunting, resulting in a roughly 2/3 reduction in the force, energy, and muscle mass required to initiate puncture of stiff materials, and even greater force reductions when the cumulative effects of abrasion are considered. We suggest that the sharpness-related force reductions lead to significant energy savings, and can also enable organisms, especially smaller ones, to puncture, cut, and grasp objects that would not be accessible with plain or biomineralized “tools”.

59 BASIC BIOLOGICAL SCIENCES↗

Accelerated cryo-EM-guided determination of three-dimensional RNA-only structures

The discovery and design of biologically important RNA molecules is outpacing three-dimensional structural characterization. In this work, we demonstrate that cryo-electron microscopy can routinely resolve maps of RNA-only systems and that these maps enable subnanometer-resolution coordinate estimation when complemented with multidimensional chemical mapping and Rosetta DRRAFTER computational modeling. Furthermore ,this hybrid ‘Ribosolve’ pipeline detects and falsifies homologies and conformational rearrangements in 11 previously unknown 119- to 338-nucleotide protein-free RNA structures: full-length Tetrahymena ribozyme, hc16 ligase with and without substrate, full-length Vibrio cholerae and Fusobacterium nucleatum glycine riboswitch aptamers with and without glycine, Mycobacterium SAM-IV riboswitch with and without S-adenosylmethionine, and the computer-designed ATP-TTR-3 aptamer with and without AMP. Simulation benchmarks, blind challenges, compensatory mutagenesis, cross-RNA homologies and internal controls demonstrate that Ribosolve can accurately resolve the global architectures of RNA molecules but does not resolve atomic details. These tests offer guidelines for making inferences in future RNA structural studies with similarly accelerated throughput.

59 BASIC BIOLOGICAL SCIENCES↗

Mitocans induce lipid flip-flop and permeabilize the membrane to signal apoptosis

Pancratistatin (PST) and narciclasine (NRC) are natural therapeutic agents that exhibit specificity toward the mitochondria of cancerous cells and initiate apoptosis. Unlike traditional cancer therapeutic agents, PST and NRC are effective, targeted, and have limited adverse effects on neighboring healthy, noncancerous cells. Currently, the mechanistic pathway of action for PST and NRC remains elusive, which in part inhibits PST and NRC from becoming efficacious therapeutic alternatives. Herein, we use neutron and x-ray scattering in combination with calcein leakage assays to characterize the effects of PST, NRC, and tamoxifen (TAM) on a biomimetic model membrane. We report an increase in lipid flip-flop half-times (t 1/2 ) (≈12.0%, ≈35.1%, and a decrease of ≈45.7%) with 2 mol percent PST, NRC, and TAM respectively. An increase in bilayer thickness (≈6.3%, ≈7.8%, and ≈7.8%) with 2 mol percent PST, NRC, and TAM, respectively, was also observed. Lastly, increases in membrane leakage (≈31.7%, ≈37.0%, and ≈34.4%) with 2 mol percent PST, NRC, and TAM, respectively, were seen. Considering the maintenance of an asymmetric lipid composition across the outer mitochondrial membrane (OMM) is crucial to eukaryotic cellular homeostasis and survival, our results suggest PST and NRC may play a role in disrupting the native distribution of lipids within the OMM. A possible mechanism of action for PST- and NRC-induced mitochondrial apoptosis is proposed via the redistribution of the native OMM lipid organization and through OMM permeabilization.

59 BASIC BIOLOGICAL SCIENCES↗

Storage-Induced Collapse of Lignin Macromolecular Structure and Its Impacts on the Biorefinery

Lignin plays a vital role in the economics of biorefineries, serving as a source of process energy and a feedstock for sustainable fuels and chemical production. While understanding lignin’s chemical composition is crucial, emerging evidence suggests that a more comprehensive understanding of its macromolecular structure is critical to explaining its complex behavior in the biorefinery. This study investigated the partial collapse of the lignin network in corn stover feedstock after harvest and storage as a result of the microbial digestion of hemicellulose. Fluorescence microscopy was used to detect the collapse of lignin in terms of lignin’s inter-molecular interaction and the re-orientation of lignin’s chromophores, by the changes in lignin’s fluorescence lifetime, anisotropy, and the number of effective emitters. With minimal sample perturbation, our in-situ microscopic results revealed lignin's coil-globule transition phenomena, which was only previously predicted by molecular dynamics modeling extracted lignin in solvent. This collapse of lignin macromolecular structure was confirmed by results from NMR, IR, Raman, and powder X-ray diffraction. We also investigated the impact of this storage-induced collapse on the downstream biorefinery processes. Our study revealed that the two major approaches for lignin valorization in the lignin-first biorefinery model, namely monomer extraction and milled wood lignin extraction, were negatively impacted by the lignin collapse. As changes during storage are a source of feedstock variability, our study highlights the importance of understanding the effect of feedstock handling on biorefinery operations and economics.

09 BIOMASS FUELS↗

The J-elongated conformation of β 2 -glycoprotein I predominates in solution: Implications for our understanding of antiphospholipid syndrome

β 2 -Glycoprotein I (β 2 GPI) is an abundant plasma protein displaying phospholipid-binding properties. Because it binds phospholipids, it is a target of antiphospholipid antibodies (aPLs) in antiphospholipid syndrome (APS), a life-threatening autoimmune thrombotic disease. Indeed, aPLs prefer membrane-bound β 2 GPI to that in solution. β 2 GPI exists in two almost equally populated redox states: oxidized, in which all the disulfide bonds are formed, and reduced, in which one or more disulfide bonds are broken. Furthermore, β 2 GPI can adopt multiple conformations (i.e. J-elongated, S-twisted, and O-circular). While strong evidence indicates that the J-form is the structure bound to aPLs, which conformation exists and predominates in solution remains controversial, and so is the conformational pathway leading to the bound state. Here, we report that human recombinant β 2 GPI purified under native conditions is oxidized. Moreover, under physiological pH and salt concentrations, this oxidized form adopts a J-elongated, flexible conformation, not circular or twisted, in which the N-terminal domain I (DI) and the C-terminal domain V (DV) are exposed to the solvent. Consistent with this model, binding kinetics and mutagenesis experiments revealed that in solution the J-form interacts with negatively charged liposomes and with MBB2, a monoclonal anti-DI antibody that recapitulates most of the features of pathogenic aPLs. We conclude that the preferential binding of aPLs to phospholipid-bound β 2 GPI arises from the ability of its preexisting J-form to accumulate on the membranes, thereby offering an ideal environment for aPL binding. We propose that targeting the J-form of β 2 GPI provides a strategy to block pathogenic aPLs in APS.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Structure and dynamics of SARS-CoV-2 proofreading exoribonuclease ExoN

High-fidelity replication of the large RNA genome of coronaviruses (CoVs) is mediated by a 3'-to-5' exoribonuclease (ExoN) in nonstructural protein 14 (nsp14), which excises nucleotides including antiviral drugs misincorporated by the low-fidelity viral RNA-dependent RNA polymerase (RdRp) and has also been implicated in viral RNA recombination and resistance to innate immunity. Here, we determined a 1.6-Å resolution crystal structure of severe acute respiratory syndrome CoV 2 (SARS-CoV-2) ExoN in complex with its essential cofactor, nsp10. The structure shows a highly basic and concave surface flanking the active site, comprising several Lys residues of nsp14 and the N-terminal amino group of nsp10. Modeling suggests that this basic patch binds to the template strand of double-stranded RNA substrates to position the 3' end of the nascent strand in the ExoN active site, which is corroborated by mutational and computational analyses. We also show that the ExoN activity can rescue a stalled RNA primer poisoned with sofosbuvir and allow RdRp to continue its extension in the presence of the chain-terminating drug, biochemically recapitulating proofreading in SARS-CoV-2 replication. Molecular dynamics simulations further show remarkable flexibility of multidomain nsp14 and suggest that nsp10 stabilizes ExoN for substrate RNA binding to support its exonuclease activity. Our high-resolution structure of the SARS-CoV-2 ExoN–nsp10 complex serves as a platform for future development of anticoronaviral drugs or strategies to attenuate the viral virulence.

60 APPLIED LIFE SCIENCES↗

Zinc binding alters the conformational dynamics and drives the transport cycle of the cation diffusion facilitator YiiP

YiiP is a secondary transporter that couples Zn 2+ transport to the proton motive force. Structural studies of YiiP from prokaryotes and Znt8 from humans have revealed three different Zn 2+ sites and a conserved homodimeric architecture. These structures define the inward-facing and outward-facing states that characterize the archetypal alternating access mechanism of transport. To study the effects of Zn 2+ binding on the conformational transition, we use cryo-EM together with molecular dynamics simulation to compare structures of YiiP from Shewanella oneidensis in the presence and absence of Zn 2+ . To enable single-particle cryo-EM, we used a phage-display library to develop a Fab antibody fragment with high affinity for YiiP, thus producing a YiiP/Fab complex. To perform MD simulations, we developed a nonbonded dummy model for Zn 2+ and validated its performance with known Zn 2+ -binding proteins. Using these tools, we find that, in the presence of Zn 2+ , YiiP adopts an inward-facing conformation consistent with that previously seen in tubular crystals. After removal of Zn 2+ with high-affinity chelators, YiiP exhibits enhanced flexibility and adopts a novel conformation that appears to be intermediate between inward-facing and outward-facing states. This conformation involves closure of a hydrophobic gate that has been postulated to control access to the primary transport site. Comparison of several independent cryo-EM maps suggests that the transition from the inward-facing state is controlled by occupancy of a secondary Zn 2+ site at the cytoplasmic membrane interface. This work enhances our understanding of individual Zn 2+ binding sites and their role in the conformational dynamics that govern the transport cycle.

59 BASIC BIOLOGICAL SCIENCES↗

Polyubiquitin ligand-induced phase transitions are optimized by spacing between ubiquitin units

Biomolecular condensates form via multivalent interactions among key macromolecules and are regulated through ligand binding and/or posttranslational modifications. One such modification is ubiquitination, the covalent addition of ubiquitin (Ub) or polyubiquitin chains to target macromolecules. Specific interactions between polyubiquitin chains and partner proteins, including hHR23B, NEMO, and UBQLN2, regulate condensate assembly or disassembly. Here, we used a library of designed polyubiquitin hubs and UBQLN2 as model systems for determining the driving forces of ligand-mediated phase transitions. Perturbations to either the UBQLN2-binding surface of Ub or the spacing between Ub units reduce the ability of hubs to modulate UBQLN2 phase behavior. By developing an analytical model based on polyphasic linkage principles that accurately described the effects of different hubs on UBQLN2 phase separation, we determined that introduction of Ub to UBQLN2 condensates incurs a significant inclusion energetic penalty. This penalty antagonizes the ability of polyUb hubs to scaffold multiple UBQLN2 molecules and cooperatively amplify phase separation. The extent to which polyubiquitin hubs promote UBQLN2 phase separation is encoded in the spacings between Ub units. This spacing is modulated by chains of different linkages and designed chains of different architectures, thus illustrating how the ubiquitin code regulates functionality via the emergent properties of the condensate. The spacing in naturally occurring linear polyubiquitin chains is already optimized to promote phase separation with UBQLN2. We expect our findings to extend to other condensates, emphasizing the importance of ligand properties, including concentration, valency, affinity, and spacing between binding sites in studies and designs of condensates.

97 MATHEMATICS AND COMPUTING↗

Functional Relevance of CASP16 Nucleic Acid Predictions as Evaluated by Structure Providers

ABSTRACT Accurate biomolecular structure prediction enables the prediction of mutational effects, the speculation of function based on predicted structural homology, the analysis of ligand binding modes, experimental model building, and many other applications. Such algorithms to predict essential functional and structural features remain out of reach for biomolecular complexes containing nucleic acids. Here, we report a quantitative and qualitative evaluation of nucleic acid structures for the CASP16 blind prediction challenge by 12 of the experimental groups who provided nucleic acid targets. Blind predictions accurately model secondary structure and some aspects of tertiary structure, including reasonable global folds for some complex RNAs; however, predictions often lack accuracy in the regions of highest functional importance. All models have inaccuracies in non‐canonical regions where, for example, the nucleic‐acid backbone bends, deviating from an A‐form helix geometry, or a base forms a non‐standard hydrogen bond (not a Watson‐Crick base pair). These bends and non‐canonical interactions are integral to forming functionally important regions such as RNA enzymatic active sites. Additionally, the modeling of conserved and functional interfaces between nucleic acids and ligands, proteins, or other nucleic acids remains poor. For some targets, the experimental structures may not represent the only structure the biomolecular complex occupies in solution or in its functional life cycle, posing a future challenge for the community.

Biochemistry & Molecular Biology↗

Molecular level insight into non-bilayer structure formation in thylakoid membranes: a molecular dynamics study

In oxygenic photosynthetic organisms, the light reactions are performed by protein complexes embedded in the lipid bilayer of thylakoid membranes (TMs). The organization of the bulk lipid molecules into bilayer structures provide optimal conditions for the build-up of the proton motive force (pmf) and its utilization for ATP synthesis. However, the lipid composition of TMs is dominated by the non-bilayer lipid species monogalactosyl diacylglycerol (MGDG), and functional plant TMs, besides the bilayer, contain large amounts of non-bilayer lipid phases. Bulk lipids have been shown to be associated with lumenal, stromal-side and marginal-region proteins and proposed to play roles in the self-assembly and photoprotection of the photosynthetic machinery. Furthermore, it has recently been pointed out that the generation and utilization of pmf for ATP synthesis according to the ‘protet’ or protonic charge transfer model Kell (Biochim Biophys Acta Bioenerg 1865(4):149504, 2024), requires high MGDG content Garab (Physiol Plant 177(2):e70230, 2025). In this study, to gain better insight into the structural and functional roles of MGDG, we employed all atom and coarse-grained molecular dynamics simulations to explore how temperature, hydration levels and varying MGDG concentrations affect the structural and dynamic properties of bilayer membranes constituted of plant thylakoid lipids. Our findings reveal that MGDG promotes increased membrane fluidity and dynamic fluctuations in membrane thickness. MGDG-rich stacked bilayers spontaneously formed inverted hexagonal phases; these transitions were enhanced at low hydration levels and at elevated but physiologically relevant temperatures. It can thus be inferred that MGDG plays important roles in heat and drought stress mechanisms.

59 BASIC BIOLOGICAL SCIENCES↗

Computational multiphysics modeling of radioactive aerosol deposition in diverse human respiratory tract geometries

The evaluation of aerosol exposure relies on generic mathematical models that assume uniform particle deposition profiles over the human respiratory tract and do not account for subject-specific characteristics. Here we introduce a hybrid-automated computational workflow that generates personalized particle deposition profiles in 3D reconstructed human airways from computed tomography scans using Computational Fluid and Particle Dynamics simulations. This is the first large-scale study to consider realistic airways variability, where 380 lower and 40 upper human respiratory tract 3D geometries are reconstructed and parameterized. The data is clustered into nine groups using random forest regression. Computational fluid and particle dynamics simulations are conducted on these representative geometries using a realistic heavy-breathing respiratory cycle and radioactive iodine-131 as a source term. Monte Carlo radiation transport simulations are performed to obtain detailed energy deposition maps. Our findings emphasize the importance of personalized studies, as minor respiratory tract variations notably influence deposition patterns rather than global parameters of the lower airways, observing more than 30% variance in the mass deposition fraction.

62 RADIOLOGY AND NUCLEAR MEDICINE↗