Search NASA⌕ Search

SEARCH · Search NASA

Results for “Structural variant”

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 307 records · Page 17
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Structure and epitope of a neutralizing monoclonal antibody that targets the stem helix of β coronaviruses

Monoclonal antibodies that retain neutralizing activity against multiple coronavirus (CoV) lineages and variants of concern (VoC) must be developed to protect against future pandemics. These broadly neutralizing MAbs (BNMAbs) may be used as therapeutics and/or to assist in the rational design of vaccines that induce BNMAbs. 1249A8 is a BNMAb that targets the stem helix (SH) region of CoV spike (S) protein and neutralizes Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) original strain, delta, and omicron VoC, Severe Acute Respiratory Syndrome CoV (SARS-CoV), and Middle East Respiratory Syndrome CoV (MERS-CoV). To understand its mechanism of action, the crystal structure of 1249A8 bound to a MERS-CoV SH peptide was determined at 2.1 Å resolution. BNMAb 1249A8 mimics the SARS-CoV-2 S loop residues 743–749, which interacts with the N-terminal end of the SH helix in the S post-fusion conformation. The conformation of 1249A8-bound SH is distinct from the SH conformation observed in the post-fusion SARS-CoV-2 S structure, suggesting 1249A8 disrupts the secondary structure and refolding events required for CoV post-fusion S to initiate membrane fusion and ultimately infection. This study provides novel insights into the neutralization mechanisms of SH-targeting CoV BNMAbs that may inform vaccine development and the design of optimal BNMAb therapeutics.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2022

Neutralizing monoclonal antibodies elicited by mosaic RBD nanoparticles bind conserved sarbecovirus epitopes

Increased immune evasion by SARS-CoV-2 variants of concern highlights the need for new therapeutic neutralizing antibodies. Immunization with nanoparticles co-displaying spike receptor-binding domains (RBDs) from eight sarbecoviruses (mosaic-8 RBD-nanoparticles) efficiently elicits cross-reactive polyclonal antibodies against conserved sarbecovirus RBD epitopes. Here, we identified monoclonal antibodies (mAbs) capable of cross-reactive binding and neutralization of animal sarbecoviruses and SARS-CoV-2 variants by screening single mouse B cells secreting IgGs that bind two or more sarbecovirus RBDs. Single-particle cryo-EM structures of antibody-spike complexes, including a Fab-Omicron complex, mapped neutralizing mAbs to conserved class 1/4 RBD epitopes. Structural analyses revealed neutralization mechanisms, potentials for intra-spike trimer cross-linking by IgGs, and induced changes in trimer upon Fab binding. In addition, we identified a mAb-resembling Bebtelovimab, an EUA-approved human class 3 anti-RBD mAb. These results support using mosaic RBD-nanoparticle vaccination to generate and identify therapeutic pan-sarbecovirus and pan-variant mAbs.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Rapid assessment of SARS-CoV-2–evolved variants using virus-like particles

Efforts to determine why new severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) variants demonstrate improved fitness have been limited to analyzing mutations in the spike (S) protein with the use of S-pseudotyped particles. In this study, we show that SARS-CoV-2 virus-like particles (SC2-VLPs) can package and deliver exogenous transcripts, enabling analysis of mutations within all structural proteins and at multiple steps in the viral life cycle. In SC2-VLPs, four nucleocapsid (N) mutations found universally in more-transmissible variants independently increased messenger RNA delivery and expression ~10-fold, and in a reverse genetics model, the serine-202→arginine (S202R) and arginine-203→methionine (R203M) mutations each produced >50 times as much virus. SC2-VLPs provide a platform for rapid testing of viral variants outside of a biosafety level 3 setting and demonstrate N mutations and particle assembly to be mechanisms that could explain the increased spread of variants, including B.1.617.2 (Delta, which contains the R203M mutation).

60 APPLIED LIFE SCIENCES↗
DOE OSTI2023

The complexity of { 11 2 ¯ 1 } twin interface structure and energetics in HCP materials

Here, the quest for a deeper understanding of the mechanical properties of HCP alloys has been an ongoing effort. Due to the lack of a sufficient number of slip systems, twin nucleation and migration originate as a substantial plastic deformation mechanism in HCP metals. The paper uncovers the energetics and atomistic mechanism of the {$11\overline{2}1$} < $11\overline{2}6$ > which is one of the most prominent twin modes in HCP materials that has not been fully understood. In view of corrugated {$11\overline{2}1$} interfaces, the positioning of lattice/motif atoms needs to be clarified. First, we establish the initial twin boundary (TB) structure unambiguously via the lattice offsets resulting in volume invariancy and correct relative positioning of the twin variants. Secondly, the TB migration mechanism is determined from a crystallographic (zero net shuffle, $Σ\vec{s}$ = 0) and energetic standpoints (minimum energy path). The Generalized Planar Fault Energy (GPFE) curves are calculated for Ti, Hf, Mg, and Zr. Such GPFE curves for this twin mode are calculated for the first time, and point to unusually high unstable energies (300-850 mJ/m 2 ), and low migration barriers (<60 mJ/m 2 ) for twinning accompanied with large shuffle magnitudes, |$\vec{s}$|, comparable to the Burger’s vector (|$\vec{s}$| ≅ $b_e$ > 0.44Å)

36 MATERIALS SCIENCE↗
DOE OSTI2023

Computational remodeling of an enzyme conformational landscape for altered substrate selectivity

Structural plasticity of enzymes dictates their function. Yet, our ability to rationally remodel enzyme conformational landscapes to tailor catalytic properties remains limited. Here, we report a computational procedure for tuning conformational landscapes that is based on multistate design of hinge-mediated domain motions. Using this method, we redesign the conformational landscape of a natural aminotransferase to preferentially stabilize a less populated but reactive conformation and thereby increase catalytic efficiency with a non-native substrate, resulting in altered substrate selectivity. Steady-state kinetics of designed variants reveals activity increases with the non-native substrate of approximately 100-fold and selectivity switches of up to 1900-fold. Structural analyses by room-temperature X-ray crystallography and multitemperature nuclear magnetic resonance spectroscopy confirm that conformational equilibria favor the target conformation. Our computational approach opens the door to targeted alterations of conformational states and equilibria, which should facilitate the design of biocatalysts with customized activity and selectivity.

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

Progressive Hedging Decomposition for Solutions of Large-Scale Process Family Design Problems

In previous work, we have introduced a mathematical model for solving a discretized version of the process family design problem. This involves two sets of decision variables. One set selects which unit module designs are included in the process platform out of a candidate set of options; the other set determines which of these unit module designs are assigned to each variant. In this work, we exploit a parallelized Progressive Hedging (PH) algorithm to solve even larger scale design problems. PH is a well-known algorithm traditionally used to solve stochastic programming problems. While our problem is not a two-stage stochastic programming problem, the structure is similar, and it can be directly mapped to the PH approach, which we employ here to solve this deterministic optimization problem. We decompose our problem by process variant. We treat the platform unit module design variables as first-stage and the assignment of unit module designs to variants as second-stage, solving the problem using mpi-sppy. We demonstrate this approach on case studies of CC, water desalination, and refrigeration.

Stinchfield, Georgia↗
DOE OSTI2023

Breaking New Ground: MB ene Route toward Selective Vinyl Double Bond Hydrogenation in Nitroarenes

Doping, or incremental substitution of one element for another, is an effective way to tailor a compound’s structure as well as its physical and chemical properties. Herein, we replaced up to 30% of Ni with Co in members of the family of layered LiNiB compounds, stabilizing the high-temperature polymorph of LiNiB while the room-temperature polymorph does not form. By studying this layered boride with in situ high-temperature powder diffraction, we obtained a distorted variant of LiNi 0.7 Co 0.3 B featuring a perfect interlayer placement of [Ni 0.7 Co 0.3 B] layers on top of each other–a structural motif not seen before in other borides. Because of the Co doping, LiNi 0.7 Co 0.3 B can undergo a nearly complete topochemical Li deintercalation under ambient conditions, resulting in a metastable boride with the formula Li 0.04 Ni 0.7 Co 0.3 B. Heating of Li 0.04 Ni 0.7 Co 0.3 B in anaerobic conditions led to yet another metastable boride, Li 0.01 Ni 0.7 Co 0.3 B, with a CoB-type crystal structure that cannot be obtained by simple annealing of Ni, Co, and B. No significant alterations of magnetic properties were detected upon Co-doping in the temperature-independent paramagnet LiNi 0.7 Co 0.3 B or its Li-deintercalated counterparts. Finally, Li 0.01 Ni 0.7 Co 0.3 B stands out as an exceptional catalyst for the selective hydrogenation of the vinyl C=C bond in 3-nitrostyrene, even in the presence of other competing functional groups. Finally, this research showcases an innovative approach to heterogeneous catalyst design by meticulously synthesizing metastable compounds.

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

Regioselectivity mechanism of the Thunbergia alata Δ 6 -16:0-acyl carrier protein desaturase

Plant plastidial acyl–acyl carrier protein (ACP) desaturases are a soluble class of diiron-containing enzymes that are distinct from the diiron-containing integral membrane desaturases found in plants and other organisms. The archetype of this class is the stearoyl-ACP desaturase which converts stearoyl-ACP into oleoyl (18:1Δ 9 cis )-ACP. Several variants expressing distinct regioselectivity have been described including a Δ 6 -16:0-ACP desaturase from black-eyed Susan vine ( Thunbergia alata ). We solved a crystal structure of the T. alata desaturase at 2.05 Å resolution. Using molecular dynamics (MD) simulations, we identified a low-energy complex between 16:0-ACP and the desaturase that would position C6 and C7 of the acyl chain adjacent to the diiron active site. The model complex was used to identify mutant variants that could convert the T. alata Δ 6 desaturase to Δ 9 regioselectivity. Additional modeling between ACP and the mutant variants confirmed the predicted regioselectivity. Furthermore, to validate the in-silico predictions, we synthesized two variants of the T. alata desaturase and analyzed their reaction products using gas chromatography-coupled mass spectrometry. Assay results confirmed that mutants designed to convert T. alata Δ 6 to Δ 9 selectivity exhibited the predicted changes. In complementary experiments, variants of the castor desaturase designed to convert Δ 9 to Δ 6 selectivity lost some of their Δ9 desaturation ability and gained the ability to desaturate at the Δ 6 position. The computational workflow for revealing the mechanistic understanding of regioselectivity presented herein lays a foundation for designing acyl-ACP desaturases with novel selectivities to increase the diversity of monoenes available for bioproduct applications.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2026

Exploring the binding properties and activities of ancestral expansins

Bacterial expansins are non-lytic proteins capable of loosening cellulose networks, offering promising applications in agriculture, biotechnology, and material science. Their ability to disrupt noncovalent interactions in biopolymer matrices such as cellulose and chitin positions them as valuable tools for upgrading abundant natural materials. However, their industrial use remains limited due to their relatively low wall-loosening activity compared to plant expansins. To address this limitation, we applied Ancestral Sequence Resurrection (ASR) to reconstruct and characterize ancient variants of the Bacillus subtilis expansin BsEXLX1. ASR is a powerful evolutionary tool that enables the inference and synthesis of ancestral proteins, allowing researchers to explore functional traits that may have been lost over time. This approach not only provides insights into protein evolution but also facilitates the design of proteins with enhanced properties, such as improved substrate affinity or structural stability. In this study, we combined biochemical and biophysical assays to evaluate the activity and binding behavior of ancestral expansins. Our results reveal that ancestral variants exhibit increased cellulose affinity, reduced binding to acidic polysaccharides, and greater salt resistance. Furthermore, these traits enhance their wall-loosening activity and demonstrate the utility of ASR in engineering surface-active proteins for industrial applications, particularly in biomass processing and cellulose modification.

09 BIOMASS FUELS↗
DOE OSTI2025

Hybrid Metal Halides: Advancing Optoelectronic Materials (Final Technical Report)

Hybrid metal halides, particularly main group halide perovskites, are a unique class of materials that offer exceptional optoelectronic properties along with a remarkable materials design space. Initial research on this class of materials was driven by the ability of the prototype hybrid 3D perovskite-structured compound methylammonium lead iodide to function as the active layer in thin film solar cells. It is now known that this class of materials comprise a large family with tunable band gaps, relatively high charge carrier mobilities in both single crystal and polycrystalline forms, and relatively low concentrations of (deleterious) electrically active states within the band gap. Beyond applications in solar cells, the potential of these materials has been extended to emitters in light emitting diodes, and active components of radiation detectors. Complementing new functionality, the design space for hybrid metal halides continues to increase with discovery of new structural motifs. In this project, the combination of organic and inorganic functionalities has been employed to open pathways to the design and synthesis new, functional hybrid metal halides. Beyond the simple perovskites, Ruddlesden-Popper and Dion-Jacobsen compounds, and other variants, such as the newly advanced “hollow” perovskites have been studied. The materials have provided routes to understanding the unique electronic properties of hybrid metal halides because of their natural quantum well structures as well as other means of controlling band gaps and band dispersions. Understanding these materials, including the design rules for their formation, their structures and compositions in bulk and in thin films form, and the role of local (non-crystallographic) structure has an important aspect of this endeavor. The goal of advancing new materials and new fundamental understanding within this deceptively simple, yet fascinating class of compounds, so richly endowed with interesting and useful functionality, has been fulfilled.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Effect of a micro-scale dislocation pileup on the atomic-scale multi-variant phase transformation and twinning

In this paper, we perform concurrent atomistic–continuum (CAC) simulations to assess the contribution of the internal stress induced by the microscale dislocation pileup at an atomically structured interface to the atomic-scale phase transformations (PTs), reverse PTs, and twinning. The main novelty of this work is to unify the atomistic description of the interface and the coarse-grained (CG) description of the lagging dislocations away from the interface within one single framework. Our major findings are: (a) the interface dynamically responds to a pileup by forming steps/ledges, the height of which is proportional to the number of dislocations arriving at the interface; (b) the pileup-induced internal stress concentration profile follows neither the classical Eshelby model nor the super-dislocation model alone, but a combination of them; (c) when the pre-sheared sample is compressed, a direct square-to-hexagonal PT occurs ahead of the pileup tip and eventually grows into a wedge shape. The two variants of the hexagonal phases form a twin with respect to each other; (d) upon a further increase of the loading, part of the newly formed hexagonal phase transforms back to the square phase. The square product phase resulting from this reverse PT forms a twin with respect to the initial square phase. All phase boundaries (PBs) and twin boundaries (TBs) are stationary and correspond to zero thermodynamic Eshelby driving forces; and (e) the microscale dislocation pileup-induced internal shear stress and the structural change at the atomic-scale interface reduces the stress required for initiating a PT by a factor of 5.5, comparing with that in the sample containing no dislocations. This work is the first characterization of the behavior of PTs/twinning resulting from the reaction between a microscale dislocation slip and an atomically structured interface. In conclusion, the gained knowledge will advance our understanding of how the multi-phase material behaves in many complex physical processes, such as the synthesis of multi-phase high-entropy alloys or superhard ceramics under high-pressure torsion, deep mantle earthquakes in geophysics, and so on, which all involve dislocation slip, PTs, twinning, and their interactions across from the atomistic to the microscale and beyond.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Optimizing iodine adsorption in functionalized metal-organic frameworks via an unprecedented positional isomerism strategy

Porous metal–organic frameworks (MOFs) have emerged as highly promising adsorbents for capturing radioiodine, a predominant fission product released during nuclear fuel reprocessing. However, systematic investigations into the correlation between MOF structure and iodine uptake capacity remain scarce. Here, we present a novel approach to enhance the iodine adsorption capacity of MOFs by optimizing linker functionalization. Using ligand-functionalized thorium-based MOFs as a structural platform, we demonstrate that ortho-amino-substitution near the node of the dicarboxylate linker significantly increases iodine adsorption capacity compared to meta-amino-substitution, where the amino groups are directed away from the node. Specifically, ortho-substituted Th-UiO-68-3,3”-(NH 2 ) 2 exhibits higher iodine uptake capacities than the meta-substituted Th-UiO-68-2,2”-(NH 2 ) 2 via both vapor diffusion-based (2.042 vs. 1.087 g/g) and solution-based (0.841 vs. 0.784 g/g) processes. Notably, the I 2 vapor adsorption capacity (2.042 g/g) of Th-UiO-68-3,3”-(NH 2 ) 2 represents the second highest among all reported Th-MOFs. Pair distribution function (PDF) studies reveal that the superior iodine uptake performance of ortho-functionalized MOFs can be attributed to the reduced steric hindrance of the amino groups compared with the meta-substituted variants. Finally, this research highlights how positional isomerism and its subtle alterations can significantly influence host–guest interactions, extending beyond simple structural considerations.

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

GPU acceleration of hybrid functional calculations in the SPARC electronic structure code

We present a Graphics Processing Unit (GPU)-accelerated version of the real-space SPARC electronic structure code for performing hybrid functional calculations in generalized Kohn–Sham density functional theory. In particular, we develop a batch variant of the recently formulated Kronecker product-based linear solver for the simultaneous solution of multiple linear systems. We then develop a modular, math kernel based implementation for hybrid functionals on NVIDIA architectures, where computationally intensive operations are offloaded to the GPUs, while the remaining workload is handled by the central processing units (CPUs). Considering bulk and slab examples, we demonstrate that GPUs enable up to 8× speedup in node-hours and 80× in core-hours compared to CPU-only execution, reducing the time to solution on V100 GPUs to around 300 s for a metallic system with over 6000 electrons, and significantly reducing the computational resources required for a given wall time.

Kohn-Sham density functional theory↗
DOE OSTI2025

Dynamic response of 17-4 stainless steel as a function of manufacturing method and heat treatment

We present a series of plate-impact experiments on 17-4 stainless steel to study the effect of manufacturing method and heat treatment on the Hugoniot elastic limit (HEL), Hugoniot, phase transformation stress, and spallation strength. Two traditional manufacturing methods were considered, wrought processing and casting, as well as two additive manufacturing methods, laser powder-bed fusion (LPBF) and wire-fed electron beam (EBAM). For both LPBF and EBAM 17-4 stainless steel variants, two billets were printed, enabling the application of two unique heat treatments. The HEL stress depended heavily on the thermal history, with the HEL increasing after the formation of Cu-rich precipitates via heat treatment. The Hugoniot response both below and above the phase transition was unaffected by the manufacturing method or heat treatment. The phase transition stress depended heavily on the thermal history, with its variation being attributed to the presence of various microstructural features. This is supported by a marked increase in the phase transition stress after precipitation hardening. These results suggest that the notion of the phase transition stress being dictated by bulk composition is an oversimplification and the stress fields generated by the meso-scale structure are a dominant force. The spallation strength was lower in the cast material compared to all other 17-4 stainless steel variants due to the presence of brittle δ-ferrite inclusions. Additionally, a drop in the tensile strain-rate was observed in the spallation response above the phase transition stress, which was hypothesized to stem from the kinetics of the reversion to the low-pressure phase during spall.

Compressive stress↗
DOE OSTI2020

Biophysical and X-ray structural studies of the (GGGTT)3GGG G-quadruplex in complex with N-methyl mesoporphyrin IX

The G-quadruplex (GQ) is a well-studied non-canonical DNA structure formed by G-rich sequences found at telomeres and gene promoters. Biological studies suggest that GQs may play roles in regulating gene expression, DNA replication, and DNA repair. Small molecule ligands were shown to alter GQ structure and stability and thereby serve as novel therapies, particularly against cancer. In this work, we investigate the interaction of a G-rich sequence, 5’-GGGTTGGGTTGGGTTGGG-3’ (T1), with a water-soluble porphyrin, N- methyl mesoporphyrin IX (NMM) via biophysical and X-ray crystallographic studies. UV-vis and fluorescence titrations, as well as a Job plot, revealed a 1:1 binding stoichiometry with an impressively tight binding constant of 30–50 μM -1 and ΔG 298 of -10.3 kcal/mol. Eight extended variants of T1 (named T2 –T9) were fully characterized and T7 was identified as a suitable candidate for crystallographic studies. We solved the crystal structures of the T1- and T7-NMM complexes at 2.39 and 2.34 Å resolution, respectively. Both complexes form a 5’-5’ dimer of parallel GQs capped by NMM at the 3’ G-quartet, supporting the 1:1 binding stoichiometry. Our work provides invaluable details about GQ-ligand binding interactions and informs the design of novel anticancer drugs that selectively recognize specific GQs and modulate their stability for therapeutic purposes.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

Structural basis of differential gene expression at eQTLs loci from high-resolution ensemble models of 3D single-cell chromatin conformations

Abstract Motivation Techniques such as high-throughput chromosome conformation capture (Hi-C) have provided a wealth of information on nucleus organization and genome important for understanding gene expression regulation. Genome-Wide Association Studies have identified numerous loci associated with complex traits. Expression quantitative trait loci (eQTL) studies have further linked the genetic variants to alteration in expression levels of associated target genes across individuals. However, the functional roles of many eQTLs in noncoding regions remain unclear. Current joint analyses of Hi-C and eQTLs data lack advanced computational tools, limiting what can be learned from these data. Results We developed a computational method for simultaneous analysis of Hi-C and eQTL data, capable of identifying a small set of nonrandom interactions from all Hi-C interactions. Using these nonrandom interactions, we reconstructed large ensembles (×105) of high-resolution single-cell 3D chromatin conformations with thorough sampling, accurately replicating Hi-C measurements. Our results revealed many-body interactions in chromatin conformation at the single-cell level within eQTL loci, providing a detailed view of how 3D chromatin structures form the physical foundation for gene regulation, including how genetic variants of eQTLs affect the expression of associated eGenes. Furthermore, our method can deconvolve chromatin heterogeneity and investigate the spatial associations of eQTLs and eGenes at subpopulation level, revealing their regulatory impacts on gene expression. Together, ensemble modeling of thoroughly sampled single-cell chromatin conformations combined with eQTL data, helps decipher how 3D chromatin structures provide the physical basis for gene regulation, expression control, and aid in understanding the overall structure-function relationships of genome organization. Availability and implementation It is available at https://github.com/uic-liang-lab/3DChromFolding-eQTL-Loci.

Du, Lin (ORCID:0009000289869812)↗
DOE OSTI2024

CG-Kit: Code Generation Toolkit for performant and maintainable variants of source code applied to Flash-X hydrodynamics simulations

CG-Kit is a new Code Generation tool-Kit that we have developed as a part of the solution for portability and maintainability for multiphysics computing applications. The development of CG-Kit is rooted in the urgent need created by the shifting landscape of high-performance computing platforms and the algorithmic complexities of a particular large-scale multiphysics application: Flash-X. To efficiently use computing resources on a heterogeneous node, an application must have a map of computation to resources and a mechanism to move the data and computation to the resources according to the map. Most existing performance portability solutions are focussed on abstracting the expression of computations so that a unified source code can be specialized to run on different resources. However, such an approach is insufficient for a code like Flash-X, which has a multitude of code components that can be assembled in various permutations and combinations to form different instances of applications. Similar challenges apply to any code that has composability, where a single specified way of apportioning work among devices may not be optimal. Additionally, use cases arise where the optimal control flow of computation may differ for different devices while the underlying numerics remain identical. This combination leads to unique challenges including handling an existing large code base in Fortran and/or C/C++, subdivision of code into a great variety of units supporting a wide range of physics and numerical methods, different parallelization techniques for distributed and shared memory systems and accelerator devices, and heterogeneity of computing platforms requiring coexisting variants of parallel algorithms. All of these challenges demand that scientific software developers apply existing knowledge about domain applications, algorithms, and computing platforms to determine custom abstractions and granularity for code generation. There is a critical lack of tools to tackle those problems. CG-Kit is designed to fill this gap by providing a user with the ability to express their desired control flow and computation-to-resource map in the form a pseudocode-like recipe. It consists of standalone tools that can be combined into highly specific and, we argue, highly effective portability and maintainability toolchains. Here we present the design of our new tools: parametrized source trees, control flow graphs, and recipes. The tools are implemented in Python. They are agnostic to the programming language of the source code targeted for code generation. In conclusion, we demonstrate the capabilities of the toolkit with two examples, first, multithreaded variants of the basic AXPY operation, and second, variants of parallel algorithms within a hydrodynamics solver, called Spark, from Flash-X that operates on block-structured adaptive meshes.

Algorithmic portability↗
DOE OSTI2022

Applying Crystallography and 19 F NMR to investigate dynamics and partner protein interactions in a long chain Flavodoxin

Flavodoxin (Fld) is a small FMN containing protein that is involved in single electron transfer. The long‐chain flavodoxin in Rhodopseudomonas palustris bacteria replaces ferredoxin as a low‐potential electron carrier when iron is scarce. Thus it is proposed to interact with the bifurcating electron transfer flavoprotein (ETF) that yields low‐potential electrons. A surface loop on Fld interacts with another of Fld's partner proteins, so we hypothesize that it also mediates Fld's interaction with ETF. To monitor interactions with ETF directly and investigate dynamics in this loop, we are using 19 F NMR in solution. 19 F is hyperresponsive to changes in its chemical environment with a chemical shift range of >300 ppm. To provide a static reference point and assess structural heterogeneity, we are also exploiting X‐ray crystallography. In this study we selectively fluorinated the five tyrosine residues in Fld. We obtained resonance assignments from 19 F spectra of Fld variants in which individual tyrosine residues have been replaced. We obtain well resolved signals for each residue, but the resonances' linewidths indicate dynamics that affects some resonances more than others. Y90 residue displays two resonances demonstrating two different conformations that interconvert slowly on an NMR time scale. Meanwhile the crystal structure solved at 2.1Å resolution reveals two molecules per asymmetric unit providing two perspectives on the details of the structure. The crystal symmetry is monoclinic in contrast to most of the other Flds, which are orthorhombic. Interestingly, the long loop bearing Y121 and Y123 is not well resolved in chain B of the crystal structure, and the NMR line of Y123 is exceptionally broad, both indicating that the loop is dynamic and capable of altering its conformation to accomodate binding to a partner protein. Future directions include monitoring the changes in the 19 F NMR of the Fld when titrated with the partner protein, temperature dependence of the NMR spectrum and relaxation studies to evaluate time scales of motions.

Khan, Sharique↗