Search NASA⌕ Search

SEARCH · Search NASA

Results for “conformance analysis”

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

A Refinement-by-Superposition -Method for (curl)- and (div)-Conforming Discretizations

Here, we present refinement-by-superposition (RBS) hp-refinement infrastructure for computational electromagnetics (CEMs), which permits exponential rates of convergence. In contrast to dominant approaches to hp-refinement for continuous Galerkin methods, which rely on explicit constraint equations, the multilevel strategy presented drastically reduces the implementation complexity. Through the RBS methodology, enforcement of continuity occurs by construction, enabling arbitrary levels of refinement with ease, and without the practical (but not theoretical) limitations of constrained-node refinement. We outline the construction of the RBS hp-method for refinement with H (curl)- and H (div)-conforming finite cells. Numerical simulations for the 2-D finite element method (FEM) solution of the Maxwell eigenvalue problem demonstrate the effectiveness of RBS hp-refinement. As an additional goal of this work, we aim to promote the use of mixed-order (low- and high-order) elements in practical CEM applications.

42 ENGINEERING↗
DOE OSTI2024

Eliminating chemo-mechanical degradation of lithium solid-state battery cathodes during >4.5 V cycling using amorphous Nb2O5 coatings

Abstract Lithium solid-state batteries offer improved safety and energy density. However, the limited stability of solid electrolytes (SEs), as well as irreversible structural and chemical changes in the cathode active material, can result in inferior electrochemical performance, particularly during high-voltage cycling (>4.3 V vs Li/Li + ). Therefore, new materials and strategies are needed to stabilize the cathode/SE interface and preserve the cathode material structure during high-voltage cycling. Here, we introduce a thin (~5 nm) conformal coating of amorphous Nb 2 O 5 on single-crystal LiNi 0.5 Mn 0.3 Co 0.2 O 2 cathode particles using rotary-bed atomic layer deposition (ALD). Full cells with Li 4 Ti 5 O 12 anodes and Nb 2 O 5 -coated cathodes demonstrate a higher initial Coulombic efficiency of 91.6% ± 0.5% compared to 82.2% ± 0.3% for the uncoated samples, along with improved rate capability (10x higher accessible capacity at 2C rate) and remarkable capacity retention during extended cycling (99.4% after 500 cycles at 4.7 V vs Li/Li + ). These improvements are associated with reduced cell polarization and interfacial impedance for the coated samples. Post-cycling electron microscopy analysis reveals that the Nb 2 O 5 coating remains intact and prevents the formation of spinel and rock-salt phases, which eliminates intra-particle cracking of the single-crystal cathode material. These findings demonstrate a potential pathway towards stable and high-performance solid-state batteries during high-voltage operation.

Science & Technology - Other Topics↗
DOE OSTI2024

A discontinuous piecewise polynomial generalized moving least squares scheme for robust finite element analysis on arbitrary grids

A variational approach is developed with a meshless discretization to enable accurate and robust numerical simulation of partial differential equations for meshes that are of poor quality. Traditional finite element methods use the mesh to both discretize the geometric domain and to define the finite element shape functions. The latter creates a dependence between the quality of the mesh and the properties of the finite element basis that may adversely affect the accuracy of the discretized problem. Here, we propose a new approach for defining finite element shape functions that breaks this dependence and separates mesh quality from the discretization quality, which we call discontinuous piecewise polynomial generalized moving least squares (DPP-GMLS). At the core of the approach is a meshless definition of the shape functions, which limits the purpose of the mesh to representing the geometric domain and integrating the basis functions without having any role in their approximation quality. The resulting non-conforming space can be utilized within a standard discontinuous Galerkin framework, providing a rigorous foundation for solving partial differential equations on low-quality meshes. We present a collection of numerical experiments demonstrating our approach in a wide range of settings: strongly coercive elliptic problems, linear elasticity in the compressible regime, and the stationary Stokes problem. We demonstrate convergence for all problems and stability for element pairs for problems which usually require inf-sup compatibility for conforming methods, also referring to a minor modification possible through the symmetric interior penalty Galerkin framework for stabilizing element pairs that would otherwise be traditionally unstable. Mesh robustness is particularly critical for elasticity, and we provide an example that our approach provides a greater than 5 x improvement in accuracy and allows for taking an 8 x larger stable timestep for a highly deformed mesh, compared to the continuous Galerkin finite element method.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Parametrization of Generalized Parton Distributions from 𝑡-Channel String Exchange in AdS Spaces

We introduce a string-based parametrization for nucleon quark and gluon generalized parton distributions (GPDs) that is valid for all skewness. Our approach leverages conformal moments, representing them as the sum of spin-𝑗 nucleon 𝐴-form factor and skewness-dependent spin-𝑗 nucleon 𝐷-form factor, derived from 𝑡-channel string exchange in AdS spaces consistent with Lorentz invariance and unitarity. This model-independent framework, satisfying the polynomiality condition due to Lorentz invariance, uses Mellin moments from empirical data to estimate these form factors. With just five Regge slope parameters, our method accurately produces various nucleon quark GPD types and symmetric nucleon gluon GPDs through pertinent Mellin-Barnes integrals. Our isovector nucleon quark GPD is in agreement with existing lattice data, promising to improve the empirical extraction and global analysis of nucleon GPDs in exclusive processes, by avoiding the deconvolution problem at any skewness, for the first time.

QCD phenomenology↗
DOE OSTI2026

Crystal structure of Methanococcus jannaschii dihydroorotase with substrate bound

Here, we report the X-ray structural analysis of dihydroorotase from Methanococcus jannaschii co-crystallized with dihydroorotate at pH 6.5. The crystals are isomorphous with the crystals of the apoenzyme, with space group P3 2 21 and unit-cell dimensions a = b = 111.4, c = 101.2 Å. The structure was refined to R = 0.169 and R free = 0.186 at a resolution of 1.87 Å at room temperature. During crystallization the degradative reaction took place, and the electron density in the active site corresponds to the substrate carbamoyl aspartate. Carbamoyl aspartate interacts with the protein in the active site in a manner similar to that observed for Escherichia coli and human dihydroorotases. However, the flexible loop (residues 140–151) adopts both conformations in the crystal, loop-out and loop-in, with the loop-out conformation having higher occupancy. This contrasts with our expectations for the flexible loop to be exclusively in the loop-in conformation, which is the conformation that it adopts to stabilize the binding of the substrate in the known systems. Additional studies with substrate analogs that resemble carbamoyl aspartate and different crystallization conditions would provide further insight into the conformation of the flexible loop in this system.

Methanococcus jannaschii↗
DOE OSTI2025

Accurate and reliable thermochemistry by data analysis of complex thermochemical networks using Active Thermochemical Tables: the case of glycine thermochemistry

Active Thermochemical Tables (ATcT) were successfully used to resolve the existing inconsistencies related to the thermochemistry of glycine, based on statistically analyzing and solving a thermochemical network that includes >3350 chemical species interconnected by nearly 35 000 thermochemically-relevant determinations from experiment and high-level theory. Here, the current ATcT results for the 298.15 K enthalpies of formation are −394.70 ± 0.55 kJ mol −1 for gas phase glycine, −528.37 ± 0.20 kJ mol −1 for solid α-glycine, −528.05 ± 0.22 kJ mol −1 for β-glycine, −528.64 ± 0.23 kJ mol −1 for γ-glycine, −514.22 ± 0.20 kJ mol −1 for aqueous undissociated glycine, and −470.09 ± 0.20 kJ mol −1 for fully dissociated aqueous glycine at infinite dilution. In addition, a new set of thermophysical properties of gas phase glycine was obtained from a fully corrected nonrigid rotor anharmonic oscillator (NRRAO) partition function, which includes all conformers. Corresponding sets of thermophysical properties of α-, β-, and γ-glycine are also presented.

Active Thermochemical Tables↗
DOE OSTI2022

SPARC: Structural properties associated with residue constraints

SPARC facilitates the generation of plausible hypotheses regarding underlying biochemical mechanisms by structurally characterizing protein sequence constraints. Such constraints appear as residues co-conserved in functionally related subgroups, as subtle pairwise correlations (i.e., direct couplings), and as correlations among these sequence features or with structural features. SPARC performs three types of analyses. First, based on pairwise sequence correlations, it estimates the biological relevance of alternative conformations and of homomeric contacts, as illustrated here for death domains. Second, it estimates the statistical significance of the correspondence between directly coupled residue pairs and interactions at heterodimeric interfaces. Third, given molecular dynamics simulated structures, it characterizes interactions among constrained residues or between such residues and ligands that: (a) are stably maintained during the simulation; (b) undergo correlated formation and/or disruption of interactions with other constrained residues; or (c) switch between alternative interactions. We illustrate this for two homohexameric complexes: the bacterial enhancer binding protein (bEBP) NtrC1, which activates transcription by remodeling RNA polymerase (RNAP) containing σ 54 , and for DnaB helicase, which opens DNA at the bacterial replication fork. Based on the NtrC1 analysis, we hypothesize possible mechanisms for inhibiting ATP hydrolysis until ADP is released from an adjacent subunit and for coupling ATP hydrolysis to restructuring of σ 54 binding loops. Based on the DnaB analysis, we hypothesize that DnaB ‘grabs’ ssDNA by flipping every fourth base and inserting it into cavities between subunits and that flipping of a DnaB-specific glutamine residue triggers ATP hydrolysis.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Machine learning-assisted profiling of a kinked ladder polymer structure using scattering

Ladder polymers consisting of fused rings in the backbone have very limited conformational freedom, which results in very different properties from traditional linear polymers. However, accurately determining their size and chain conformations from solution scattering remains a challenge. Their chain conformations of kinked ladder polymers are largely governed by the structures and relative orientations or configurations of the repeat units, unlike conventional polymer chains whose bending angles between repeat units follow a unimodal Gaussian distribution. Meanwhile, traditional scattering models for polymer chains do not account for these unique structural features. This work introduces a novel approach that integrates machine learning with Monte Carlo simulations to construct a model that can describe the geometry of a type of kinked CANAL ladder polymers. We first develop a Monte Carlo simulation model for sampling the configuration space of CANAL ladder polymers, where each repeat unit is modeled as a biaxial segment. Then, we establish a machine learning-assisted scattering analysis framework based on Gaussian Process Regression. Finally, we conduct small-angle neutron scattering experiments on a CANAL ladder polymer solution to apply our approach. Our method uncovers structural features of such ladder polymers that conventional methods fail to capture.

Ding, Lijie [Oak Ridge National Laboratory (ORNL),↗
DOE OSTI2024

Co-hydrotreatment of Bio-oil and Waste Cooking Oil to Produce Transportation Fuels

This paper reports the co-hydrotreatment of the heavy bio-oil fraction with waste cooking oil (WCO) using NiMo/γ-Al 2 O 3 catalyst, followed by the distillation of resulting deoxygenated oil and the characterization of resulting fuel cuts. The heavy BTG bio-oil fraction was obtained by removing the very reactive light-oxygenated compounds via rotary evaporation, subsequently mixed with 1-butanol. The resulting oil was blended with WCO and subjected to a two-step co-hydrotreatment process. The first step, called “stabilization,” is aimed at saturating highly reactive hydrogen-deficient compounds. The second step, called “deoxygenation,” aimed to remove bio-oil oxygen, primarily as H 2 O. This study examined the impact of varying bio-oil concentrations (0, 10, 20, 30, 40 wt.% of WCO) on the upgraded oil's yield, composition, and fuel properties. The resulting hydrotreated oil was distilled into gasoline-range, kerosene-range, and diesel-range hydrocarbons at <150 °C, 150 to 250 °C, and 250 to 350 °C, respectively. The yield of the hydrotreated oil indicates that as the bio-oil concentration increases, the amounts of coke (0.7 to 2.4 %) and water (2 to 10 wt. %) increase while the organic layer yields decrease (80 to 63 %). The coke yield was comparable to the coke yield obtained when co-processing the pyrolytic lignin fraction. This suggests that coke is formed from both the sugar oligomers and the lignin-derived oligomers. The UV-fluorescence analysis on the hydrotreated oil shows that more polycondensed and conjugated ring compounds formed as the bio-oil concentration is increased. These compounds are precursors of coke. FTIR results showed that most raw materials were converted to biofuels after the hydrotreatment. To achieve less than 1 wt. % of coke yield, blends with up to 20 wt. % pyrolysis oil should be used. An increase in bio-oil concentration leads to a slight increase in gasoline yield and a decrease in kerosene and diesel yields. The identified carbon species found in the fuel cuts include n-paraffin, iso-paraffin, cycloparaffin, and aromatics. Further, the jet fuel cut (kerosene) was characterized by density, surface tension, and viscosity. Our product conforms to the standard specifications for sustainable aviation fuels (Jet A-1). Further research is suggested to fine-tune the operating parameters for achieving reduced coke yield and enhanced kerosene yield.

09 BIOMASS FUELS↗
DOE OSTI2024

Ethylene Polymerizations Catalyzed by Fluorinated “Sandwich” Diimine-Nickel and Palladium Complexes

Nickel and palladium complexes bearing “sandwich” diimine ligands with perfluorinated aryl caps have been synthesized, characterized, and explored in ethylene polymerization reactions. The X-ray crystallographic analysis of the precatalysts 16 and 6b show differences from their non-fluorinated analogues 17 and 19, with the perfluorinated aryl caps centered precisely over the nickel and palladium centers which results in higher buried volumes of the metal centers relative to the non-fluorinated analogs. The sandwich diimine palladium complexes 5a and 5b containing perfluorinated aryl caps polymerize ethylene in a controlled fashion with activities that are substantially increased compared with their non-fluorinated analogues. Migratory insertion rates in relevant methyl ethylene complexes agree with activities exhibited in bulk polymerization experiments. DFT studies suggest that facility of ethylene rotation from its preferred orientation perpendicular to the Pd-alkyl bond into a parallel in-plane conformation contributes to the higher polymerization activity for 5b relative to 18a. For these palladium systems polymer molecular weights can be controlled via hydrogen addition (hydrogenolysis) which is unusual for late transition metal-catalyzed olefin polymerizations with no catalyst deactivation occurring. Sandwich diimine nickel complexes 6a and 6b with perfluorinated aryl caps show ethylene polymerization activities that are about half that of classical tetraisopropyl-substituted catalyst 2 but again are more active than the analogous nonfluorinated sandwich complexes. Ethylene polymerizations exhibit living behavior and branched ultra-high molecular weight polyethylenes (UHMWPEs) with very low molecular weight distributions (less than 1.1) are obtained. The activated nickel catalysts are stable in the absence of monomer and show good long-term stability at 25 °C.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Reinterpreting π-stacking

The nature of π–π interactions has long been debated. The term “π-stacking” is considered by some to be a misnomer, in part because overlapping π-electron densities are thought to incur steric repulsion, and the physical origins of the widely-encountered “slip-stacked” motif have variously been attributed to either sterics or electrostatics, in competition with dispersion. Here, we use quantum-mechanical energy decomposition analysis to investigate π–π interactions in supramolecular complexes of polycyclic aromatic hydrocarbons, ranging in size up to realistic models of graphene, and for comparison we perform the same analysis on stacked complexes of polycyclic saturated hydrocarbons, which are cyclohexane-based analogues of graphane. Our results help to explain the short-range structure of liquid hydrocarbons that is inferred from neutron scattering, trends in melting-point data, the interlayer separation of graphene sheets, and finally band gaps and observation of molecular plasmons in graphene nanoribbons. Analysis of intermolecular forces demonstrates that aromatic π–π interactions constitute a unique and fundamentally quantum-mechanical form of non-bonded interaction. Not only do stacked π–π architectures enhance dispersion, but quadrupolar electrostatic interactions that may be repulsive at long range are rendered attractive at the intermolecular distances that characterize π-stacking, as a result of charge penetration effects. The planar geometries of aromatic sp 2 carbon networks lead to attractive interactions that are “served up on a molecular pizza peel”, and adoption of slip-stacked geometries minimizes steric (rather than electrostatic) repulsion. The slip-stacked motif therefore emerges not as a defect induced by electrostatic repulsion but rather as a natural outcome of a conformational landscape that is dominated by van der Waals interactions (dispersion plus Pauli repulsion), and is therefore fundamentally quantum-mechanical in its origins. This reinterpretation of the forces responsible for π-stacking has important implications for the manner in which non-bonded interactions are modeled using classical force fields, and for rationalizing the prevalence of the slip-stacked π–π motif in protein crystal structures.

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

PyPop: a mature open-source software pipeline for population genomics

Python for Population Genomics (PyPop) is a software package that processes genotype and allele data and performs large-scale population genetic analyses on highly polymorphic multi-locus genotype data. In particular, PyPop tests data conformity to Hardy-Weinberg equilibrium expectations, performs Ewens-Watterson tests for selection, estimates haplotype frequencies, measures linkage disequilibrium, and tests significance. Standardized means of performing these tests is key for contemporary studies of evolutionary biology and population genetics, and these tests are central to genetic studies of disease association as well. Here, we present PyPop 1.0.0, a new major release of the package, which implements new features using the more robust infrastructure of GitHub, and is distributed via the industry-standard Python Package Index. New features include implementation of the asymmetric linkage disequilibrium measures and, of particular interest to the immunogenetics research communities, support for modern nomenclature, including colon-delimited allele names, and improvements to meta-analysis features for aggregating outputs for multiple populations.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Role of Bottlebrush Additives on the Structure of Block Copolymers in the Bulk and Thin Films

Blending block copolymers (BCP) with additives is a useful approach for controlling BCP morphology and properties. In athermal systems, blends of BCPs with polymer additives having very high molecular (M n ) mass generally result in macrophase separation. Bottlebrush polymers, which consist of a linear backbone and grafted side chains, present an interesting alternative where the overall M n of the system can be very large but the low M n side chains may drive miscibility with the BCP. Here, in this study, a bottlebrush with a polynorbornene backbone and polystyrene (PS) side chains is blended with PS-b-poly(methyl methacrylate) (PS-b-PMMA) of varying M n , and the resulting morphologies are examined in both the bulk and thin films. Two different M n of PS-b-PMMA were used in the bulk study, and the analysis of small-angle X-ray scattering data shows that the blends were miscible and lamellar at all concentrations. This deviates from reference series of both low and high M n linear polymer additives, which either showed morphological transitions from lamellae to cylinders (low M n ) or were immiscible at all mass fractions studied (high M n ). The relative molecular mass of the side chain (N SC ) and the corresponding component in the BCP (N A ) dictate the distribution of the bottlebrush throughout the BCP, analogous to BCP/linear blends or grafted nanoparticles in a homopolymer matrix. The studies on thin films show a thickness dependence for bottlebrush mass fractions at or above 0.17, a behavior which may be driven by conformational changes of the bottlebrush upon confinement.

36 MATERIALS SCIENCE↗
DOE OSTI2025

The molecular architecture distinctions between compression, opposite and normal wood of Pinus radiata

In gymnosperms compression wood is a specialised type of structural cell wall formed in response to biomechanical stresses. The differences in terms of gross structure, ultrastructure and chemistry are well-known. However, the differences between compression wood, normal wood, and opposite wood regarding the arrangements and interactions of the various polymers and water within their cell walls still needs to be established. The analysis of 13 C-labelled Pinus radiata by solid-state NMR spectroscopy and other complementary techniques revealed several new aspects of compression and opposite wood molecular architecture. Compared to normal wood, compression wood has a lower water content, its overall nanoporosity is reduced, and the water and matrix polymers have a lower molecular mobility. Galactan, which is a specific marker of compression wood, is broadly distributed within the cell wall, disordered, and not aligned with cellulose, and is found to be in close proximity to xylan. Dehydroabietic acid (a resin acid) is immobilised and close to the H-lignin only in compression wood. Although the overall molecular mobility of normal wood and opposite wood are similar, opposite wood has different arabinose conformations, a large increase in the amount of chain ends, contains significantly more galactan and has additional unassigned mobile components highlighting the different molecular arrangement of cell wall polymers in opposite and normal wood.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Computational Tools and Workflows for Quantitative Risk Assessment and Decision Support for Geologic Carbon Storage Sites: Progress and Insights from the U.S. DOE’s National Risk Assessment Partnership

The 2005 Intergovernmental Panel on Climate Change (IPCC) Special Report on CCS raised the profile of CO2 capture and storage (CCS) as an important technology for reducing greenhouse gas (GHG) emissions. CCS is now recognized as a key component of most climate change mitigation scenarios. Since publication of that report the international research, development, and deployment (RD&D) community has advanced key technical aspects, clarified regulatory requirements, explored value chain and infrastructure solutions, and developed incentive paradigms to enable and promote large-scale deployment of CCS. These efforts have included research to better characterize geologic storage resources, to improve injection performance and storage efficiency, to assess and manage subsurface environmental risks, and to advance monitoring technologies to assure system conformance. These efforts have helped to build confidence in the viability of geologic carbon storage (GCS), but stakeholder concerns about long-term risks and liability associated with GCS remain a hurdle to broad acceptance and large-scale deployment of CCS. Since 2010, the U.S. DOE’s National Risk Assessment Partnership (NRAP) – a research collaboration between five contributing national laboratories – has worked to establish and demonstrate methods and tools to quantify and manage the subsurface environmental risks associated with GCS, amidst uncertainty. This work supports the Office of Fossil Energy and Carbon Management Carbon Transport and Storage Program’s goal of advancing safe and secure commercial-scale GCS deployment. To address the technical challenge of simulating the physical response of the GCS site to large-scale CO2 injection, NRAP has adopted an approach that relies on coupling computationally efficient reduced-order and/or data-driven proxy models of important system components (i.e., storage reservoir, sealing caprock, leakage pathways, intermediate formations, overlying groundwater aquifers, and the atmosphere) in integrated assessment framework. That integrated model of the physical system is complemented with fit-for purpose functionality to support site characterization and risk-related decisions. The recently released NRAP Phase II toolset includes the Open-Source Integrated Assessment Model (NRAP-Open-IAM) for evaluation of trends in leakage risk and potential impact, tools to support monitoring design optimization (Designs for Risk Evaluation and Management – DREAM v3.0 and Passive Seismic Monitoring Tool - PSMT), and tools for state of stress evaluation (State-of-Stress Analysis Tool - SOSAT) and forecasting induced seismicity risk. The NRAP team has also released a pair of reports describing conceptual workflows to incorporate physics-based, quantitative risk assessment into many of the design, planning, operation, and closure decisions for GCS projects. An online catalogue highlights published studies where these tools and methods are demonstrated. In this presentation, the utility of these products to assess risks and address key stakeholder questions will be highlighted through examples, and related insights about the safety and security of geologic carbon storage in qualified storage sites will be discussed. The prospect of rapid, large-scale deployment of GCS technology to aggressively reduce anthropogenic CO2 emissions requires careful consideration of interference between multiple commercial-scale storage projects within a basin. Going forward, NRAP is expanding and adapting site-scale risk quantification tools and methods to enable assessment of risks and inform management decisions for basin-scale deployment. Increasingly, this work will leverage next-generation approaches for surrogate modelling, fast prediction, and advanced visualization enabled by machine learning and artificial intelligence to promote virtual learning, scenario evaluation, and augment risk-based decision making.

quantitative risk assessment, geologic carbon stor↗
DOE OSTI2026

From 2D to 4D: a containerized workflow and browser to explore dynamic chromatin architecture

Background Characterizing the physical organization of the genome is essential for understanding long-range gene regulation, chromatin compartmentalization, and epigenetic accessibility. Hi-C experiments generate two-dimensional (2D) genome-wide contact maps of chromatin interactions by capturing the spatial proximity between genomic loci, which reveal interaction frequencies but lack the spatial resolution needed to interpret the three-dimensional (3D) genome structure(s). Emerging evidence suggests that epigenetic regulation is closely linked to 3D genome architecture, and that structural changes over time (4D) drive key biological processes in development, disease, and environmental response. Thus, integrating 3D structure with functional data is critical for a more complete understanding of genome regulation. Previous work, most notably the 4DHiC chromosome modeling framework, has shown that physical multi-dimensional modeling approaches rooted in polymer physics and molecular dynamics can resolve these structures at biologically meaningful resolutions by integrating temporal Hi-C data with physical constraints to uncover dynamic chromosome reorganization. Thus, molecular dynamics simulations, constrained by Hi-C contact matrices, can resolve fine-scale structural changes and reveal functionally significant transitions in chromatin conformation. Results Herein, we present the 4D Genome Browser Workflow (4DGBWorkflow) and the 4D Genome Browser (4DGB). The algorithm is based on the 4DHiC method, and the containerized tool is an end-to-end workflow that can transform, filter, and view 4D epigenomics and chromatin datasets, allowing non-specialists to apply three-dimensional modeling principles to diverse datasets and experimental conditions. The software executes on a laptop running macOS, Linux or Windows. From input Hi-C files (.hic), the 4DGBWorkflow produces 3D reconstructions of chromosomes, integrates the reconstruction with track data (e.g., epigenetic marks, transcriptome profiles), and provides comparative visualization of the results in a single workflow. Conclusions The 4DGBWorkflow and 4D Genome Browser are open-source tools for comparative analysis and visualization of 4D chromosome datasets, including chromatin architecture and epigenomic signals. Automatic integration of Hi-C data with molecular dynamics democratizes the construction of time resolved 3D genome structures, simplifying complex simulations and data integration schemes.

3D Genome Browser↗
DOE OSTI2023

A Bis(imidazole)-based cysteine labeling tool for metalloprotein assembly

Precise metal-protein coordination by design remains a considerable challenge. Polydentate, high-metal-affinity protein modifications, both chemical and recombinant, can enable metal localization. However, these constructs are often bulky, conformationally and stereochemically ill-defined, or coordinately saturated. Here, we expand the biomolecular metal-coordination toolbox with the irreversible attachment to cysteine of bis(1-methylimidazol-2-yl)ethene (“BMIE”), which generates a compact imidazole-based metal-coordinating ligand. Conjugate additions of small-molecule thiols (thiocresol and N-Boc-Cys) with BMIE confirm general thiol reactivity. The BMIE adducts are shown to complex the divalent metal ions Cu ++ and Zn ++ in bidentate (N 2 ) and tridentate (N 2 S*) coordination geometries. Cysteine-targeted BMIE modification (>90% yield at pH 8.0) of a model protein, the S203C variant of carboxypeptidase G2 (CPG2), measured with ESI-MS, confirms its utility as a site-selective bioconjugation method. ICP-MS analysis confirms mono-metallation of the BMIE-modified CPG2 protein with Zn ++ , Cu ++ , and Co ++ . EPR characterization of the BMIE-modified CPG2 protein reveals the structural details of the site selective 1:1 BMIE-Cu ++ coordination and symmetric tetragonal geometry under physiological conditions and in the presence of various competing and exchangeable ligands (H 2 O/HO – , tris, and phenanthroline). An X-ray protein crystal structure of BMIE-modified CPG2-S203C demonstrates that the BMIE modification is minimally disruptive to the overall protein structure, including the carboxypeptidase active sites, although Zn ++ metalation could not be conclusively discerned at the resolution obtained. The carboxypeptidase catalytic activity of BMIE-modified CPG2-S203C was also assayed and found to be minimally affected. Finally, these features, combined with ease of attachment, define the new BMIE-based ligation as a versatile metalloprotein design tool, and enable future catalytic and structural applications.

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

Proton Transfer Dynamics in the Aprotic Proton Accepting Solvent 1-Methylimidazole

The dynamics of proton transfer to the aprotic solvent 1-methylimidazole (MeIm, proton acceptor) from the photoacid 8-hydroxypyrene-1,3,6-trisulfonic acid trisodium salt (HPTS) were investigated using fast fluorescence measurements. The closely related molecule, 8- methoxypyrene-1,3,6-trisulfonic acid trisodium salt (MPTS), which is not a photoacid, was also studied for comparison. Following optical excitation, the wavelength-dependent population dynamics of HPTS in MeIm resulting from the deprotonation process were collected over the entire fluorescence emission window. Analysis of the time-dependent fluorescence spectra revealed four distinct fluorescence bands that appear and decay on different time scales. We label these four states as: protonated (P), associated I (A I ), associated II (A II ), and deprotonated (D). We find that the simple kinetic scheme of P → A I → A II → D is not consistent with the data. Instead, the kinetic scheme that describes the data has P decaying into A I , which mainly goes on to deprotonation (D), but AI can also feed into A II . A II can return to A I or decay to the ground state, but does not deprotonate within experimental error. Finally, quantum chemistry and excitedstate QM/MM Born-Oppenheimer molecular dynamics simulations indicate that A I and A II are two H-bonding conformations of MeIm to the HPTS hydroxyl, axial and equatorial, respectively.

14 SOLAR ENERGY↗