Search NASA⌕ Search

SEARCH · Search NASA

Results for “data structures”

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 883 records · Page 49

An optimized FM-index library for nucleotide and amino acid search

Abstract Background Pattern matching is a key step in a variety of biological sequence analysis pipelines. The FM-index is a compressed data structure for pattern matching, with search run time that is independent of the length of the database text. Implementation of the FM-index is reasonably complicated, so that increased adoption will be aided by the availability of a fast and flexible FM-index library. Results We present AvxWindowedFMindex (AWFM-index), a lightweight, open-source, thread-parallel FM-index library written in C that is optimized for indexing nucleotide and amino acid sequences. AWFM-index introduces a new approach to storing FM-index data in a strided bit-vector format that enables extremely efficient computation of the FM-index occurrence function via AVX2 bitwise instructions, and combines this with optional on-disk storage of the index’s suffix array and a cache-efficient lookup table for partial k-mer searches. The AWFM-index performs exact match count and locate queries faster than SeqAn3’s FM-index implementation across a range of comparable memory footprints. When optimized for speed, AWFM-index is $$\sim $$ ∼ 2–4x faster than SeqAn3 for nucleotide search, and $$\sim $$ ∼ 2–6x faster for amino acid search; it is also $$\sim $$ ∼ 4x faster with similar memory footprint when storing the suffix array in on-disk SSD storage. Conclusions AWFM-index is easy to incorporate into bioinformatics software, offers run-time performance parameterization, and provides clients with FM-index functionality at both a high-level (count or locate all instances of a query string) and low-level (step-wise control of the FM-index backward-search process). The open-source library is available for download at https://github.com/TravisWheelerLab/AvxWindowFmIndex.

59 BASIC BIOLOGICAL SCIENCES↗

A Taxonomy for Classification and Comparison of Dataflows for GNN Accelerators

Recently, Graph Neural Networks (GNNs) have received a lot of interest because of their success in learning representations from graph structured data. However, GNNs exhibit different compute and memory characteristics compared to traditional Deep Neural Networks (DNNs). Graph convolutions require feature aggregations from neighboring nodes (known as the aggregation phase), which leads to highly irregular data accesses. GNNs also have a very regular compute phase that can be broken down to matrix multiplications (known as the combination phase). All recently proposed GNN accelerators utilize different dataflows and microarchitecture optimizations for these two phases. Different communication strategies between the two phases have been also used. However, as more custom GNN accelerators are proposed, the harder it is to qualitatively classify them and quantitatively contrast them. In this work, we present a taxonomy to describe several diverse dataflows for running GNN inference on accelerators. This provides a structured way to describe and compare the design-space of GNN accelerators.

97 MATHEMATICS AND COMPUTING↗

Adaptively Refined Euler and Navier-Stokes Solutions with a Cartesian-Cell Based Scheme

A Cartesian-cell based scheme with adaptive mesh refinement for solving the Euler and Navier-Stokes equations in two dimensions has been developed and tested. Grids about geometrically complicated bodies were generated automatically, by recursive subdivision of a single Cartesian cell encompassing the entire flow domain. Where the resulting cells intersect bodies, N-sided 'cut' cells were created using polygon-clipping algorithms. The grid was stored in a binary-tree data structure which provided a natural means of obtaining cell-to-cell connectivity and of carrying out solution-adaptive mesh refinement. The Euler and Navier-Stokes equations were solved on the resulting grids using an upwind, finite-volume formulation. The inviscid fluxes were found in an upwinded manner using a linear reconstruction of the cell primitives, providing the input states to an approximate Riemann solver. The viscous fluxes were formed using a Green-Gauss type of reconstruction upon a co-volume surrounding the cell interface. Data at the vertices of this co-volume were found in a linearly K-exact manner, which ensured linear K-exactness of the gradients. Adaptively-refined solutions for the inviscid flow about a four-element airfoil (test case 3) were compared to theory. Laminar, adaptively-refined solutions were compared to accepted computational, experimental and theoretical results.

Coirier, William J.↗

Inputs to GCAM-USA: IM3 Phase 2 Experiments

Overview This dataset contains XML input files for the IM3 Phase 2 version of GCAM-USA. The files are organized into two categories: Scenario-specific inputs represent hydroclimate and socioeconomic effects on water availability, heating and cooling degree-hours, and agricultural productivity. They support eight IM3 canonical scenarios: rcp45cooler_ssp3 rcp45cooler_ssp5 rcp45hotter_ssp3 rcp45hotter_ssp5 rcp85cooler_ssp3 rcp85cooler_ssp5 rcp85hotter_ssp3 rcp85hotter_ssp5 Model-improvement inputs extend GCAM-USA v5.3 with updated representations of coal and nuclear power plant retirements, electricity trade among U.S. interconnections, offshore carbon storage costs, and groundwater depletion constraints. Data structure Scenario-specific inputs rcp45_runoff/ and rcp85_runoff/XML files describing water availability by HUC2 basin under the RCP 4.5 and RCP 8.5 scenarios. rcp45_hdcd/ and rcp85_hdcd/XML files containing monthly-day and monthly-night heating and cooling degree-hours at the U.S. state level for different RCP-SSP combinations. rcp45_agyields/ and rcp85_agyields/XML files describing changes in agricultural productivity at the intersection of GCAM regions and HUC2 water basins for different RCP-SSP combinations. rcp45_emissions_pathway/The emissions-constraint XML file used to represent the RCP 4.5 pathway. Model-improvement inputs core_retire/Updates coal-fired power plant retirement schedules based on New England ISO. GCAMUSA_IM3_elec_trade_interconnect.xmlRestricts electricity trade to occur within the ERCOT, WECC, and IE interconnections. nuclear_USA.xmlUpdates the retirement schedules of the Diablo Canyon and Palisades nuclear power plants. high_cost_offshore_carbon.xmlUpdates the assumed cost of offshore carbon storage. water_supply_constrained_gleeson_5pct.xmlReplaces WaterGAP historical groundwater-depletion estimates with data from the Gleeson dataset and limits groundwater extraction to 5% of the available groundwater in each Superwell grid cell. How to use the data This dataset is designed for use with the IM3 version of GCAM-USA. Download or clone GCAM-USA from the IM3 GCAM GitHub repository at https://github.com/IMMM-SFA/gcam-core and check out the gcam-usa-im3 branch. Place the downloaded folder im3scenarios in the gcam-core/input directory while preserving the provided folder structure.

Energy↗

Cryo-EM of the pathogenic VCP variant R155P reveals long-range conformational changes in the D2 ATPase ring

Single amino acid mutations in valosin containing protein (VCP/p97), a highly conserved member of the ATPases associated with diverse cellular activities (AAA) family of ATPases has been linked to a severe degenerative disease affecting brain, muscle and bone tissue. Previous studies have demonstrated the role of VCP mutations in altering the ATPase activity of the D2 ring; however the structural consequences of these mutations remain unclear. In this study, we report the three-dimensional (3D) map of the pathogenic VCP variant, R155P, as revealed by single-particle Cryo-Electron Microscopy (EM) analysis at 14 Å resolution. We show that the N-terminal R155P mutation induces a large structural reorganisation of the D2 ATPase ring. Results from docking studies using crystal structure data of available wild-type VCP in the EM density maps indicate that the major difference is localized at the interface between two protomers within the D2 ring. Consistent with a conformational change, the VCP R155P variant shifted the isoelectric point of the protein and reduced its interaction with its well-characterized cofactor, nuclear protein localization-4 (Npl4). Together, our results demonstrate that a single amino acid substitution in the N-terminal domain can relay long-range conformational changes to the distal D2 ATPase ring. Our results provide the first structural clues of how VCP mutations may influence the activity and function of the D2 ATPase ring. - Highlights: • p97{sub R155P} and p97{sub A232E} decrease the ability of p97 to bind to its co-factor Npl4. • p97{sub R155P} has a different isoelectric point than that of p97{sub R95G}, p97{sub A232E} and p97{sub WT}. • Mutation R155P changes principally the conformation of the D2 ring. • Mutation R155P modifies the interface between two protomers within the D2 ring.

60 APPLIED LIFE SCIENCES↗

Structural basis of mutant-selectivity and drug-resistance related to CO-1686

Non-small-cell lung cancers (NSCLCs) caused by activating mutations in the kinase domain of epidermal growth factor receptor (EGFR) initially respond to first-generation reversible drugs gefitinib and erlotinib. However, clinical efficacy is limited due to the development of drug-resistance that in more than half of the cases are driven by the secondary T790M mutation. CO-1686 is one of the third generation irreversible inhibitors that inhibits EGFR activating mutants, including those with concurrent T790M, while avoiding the off-target toxicity owing to inhibition of wild-type EGFR in treating EGFR mutation-positive NSCLCs. Despite the remarkable success, the experimentally determined structure of this agent in complex with EGFR T790M remains unknown. In this study, we determined crystal structures of EGFR T790M or L858R mutants covalently bound by CO-1686. Based on these structural data, we can explain why CO-1686 irreversibly inhibits EGFR and selectively prefers T790M, which may help improving this or similar compounds, and explain why EGFR L718Q and L844V mutations incur resistance to this agent.

60 APPLIED LIFE SCIENCES↗

Consistent performance of large language models in rare disease diagnosis across ten languages and 4917 cases

Background Large language models (LLMs) are increasingly used medicine for diverse applications including differential diagnostic support. The training data used to create LLMs such as the Generative Pretrained Transformer (GPT) predominantly consist of English-language texts, but LLMs could be used across the globe to support diagnostics if language barriers could be overcome. Initial pilot studies on the utility of LLMs for differential diagnosis in languages other than English have shown promise, but a large-scale assessment on the relative performance of these models in a variety of European and non-European languages on a comprehensive corpus of challenging rare-disease cases is lacking. Methods We created 4917 clinical vignettes using structured data captured with Human Phenotype Ontology (HPO) terms with the Global Alliance for Genomics and Health (GA4GH) Phenopacket Schema. These clinical vignettes span a total of 360 distinct genetic diseases with 2525 associated phenotypic features. We used translations of the Human Phenotype Ontology together with language-specific templates to generate prompts in English, Chinese, Czech, Dutch, French, German, Italian, Japanese, Spanish, and Turkish. We applied GPT-4o, version gpt-4o-2024-08-06, and the medically fine-tuned Meditron3-70B to the task of delivering a ranked differential diagnosis using a zero-shot prompt. An ontology-based approach with the Mondo disease ontology was used to map synonyms and to map disease subtypes to clinical diagnoses in order to automate evaluation of LLM responses. Findings For English, GPT-4o placed the correct diagnosis at the first rank 19.9% and within the top-3 ranks 27.0% of the time. In comparison, for the nine non-English languages tested here the correct diagnosis was placed at rank 1 between 16.9% and 20.6%, within top-3 between 25.4% and 28.6% of cases. The Meditron3 model placed the correct diagnosis within the first 3 ranks for 20.9% of cases in English and between 19.9% and 24.0% for the other nine languages. Interpretation The differential diagnostic performance of LLMs across a comprehensive corpus of rare-disease cases was largely consistent across the ten languages tested. This suggests that the utility of LLMs in clinical settings may extend to non-English clinical settings.

Artificial intelligence↗

Use of Generalized Fluid System Simulation Program (GFSSP) for Teaching and Performing Senior Design Projects at the Educational Institutions

This paper describes the experience of the authors in using the Generalized Fluid System Simulation Program (GFSSP) in teaching Design of Thermal Systems class at University of Alabama in Huntsville. GFSSP is a finite volume based thermo-fluid system network analysis code, developed at NASA/Marshall Space Flight Center, and is extensively used in NASA, Department of Defense, and aerospace industries for propulsion system design, analysis, and performance evaluation. The educational version of GFSSP is freely available to all US higher education institutions. The main purpose of the paper is to illustrate the utilization of this user-friendly code for the thermal systems design and fluid engineering courses and to encourage the instructors to utilize the code for the class assignments as well as senior design projects. The need for a generalized computer program for thermofluid analysis in a flow network has been felt for a long time in aerospace industries. Designers of thermofluid systems often need to know pressures, temperatures, flow rates, concentrations, and heat transfer rates at different parts of a flow circuit for steady state or transient conditions. Such applications occur in propulsion systems for tank pressurization, internal flow analysis of rocket engine turbopumps, chilldown of cryogenic tanks and transfer lines, and many other applications of gas-liquid systems involving fluid transients and conjugate heat and mass transfer. Computer resource requirements to perform time-dependent, three-dimensional Navier-Stokes computational fluid dynamic (CFD) analysis of such systems are prohibitive and therefore are not practical. Available commercial codes are generally suitable for steady state, single-phase incompressible flow. Because of the proprietary nature of such codes, it is not possible to extend their capability to satisfy the above-mentioned needs. Therefore, the Generalized Fluid System Simulation Program (GFSSP1) has been developed at NASA Marshall Space Flight Center (MSFC) as a general fluid flow system solver capable of handling phase changes, compressibility, mixture thermodynamics and transient operations. It also includes the capability to model external body forces such as gravity and centrifugal effects in a complex flow network. The objectives of GFSSP development are: a) to develop a robust and efficient numerical algorithm to solve a system of equations describing a flow network containing phase changes, mixing, and rotation; and b) to implement the algorithm in a structured, easy-to-use computer program. The analysis of thermofluid dynamics in a complex network requires resolution of the system into fluid nodes and branches, and solid nodes and conductors as shown in Figure 1. Figure 1 shows a schematic and GFSSP flow circuit of a counter-flow heat exchanger. Hot nitrogen gas is flowing through a pipe, colder nitrogen is flowing counter to the hot stream in the annulus pipe and heat transfer occurs through metal tubes. The problem considered is to calculate flowrates and temperature distributions in both streams. GFSSP has a unique data structure, as shown in Figure 2, that allows constructing all possible arrangements of a flow network with no limit on the number of elements. The elements of a flow network are boundary nodes where pressure and temperature are specified, internal nodes where pressure and temperature are calculated, and branches where flowrates are calculated. For conjugate heat transfer problems, there are three additional elements: solid node, ambient node, and conductor. The solid and fluid nodes are connected with solid-fluid conductors. GFSSP solves the conservation equations of mass and energy, and equation of state in internal nodes to calculate pressure, temperature and resident mass. The momentum conservation equation is solved in branches to calculate flowrate. It also solves for energy conservation equations to calculate temperatures of solid nodes. The equations are coupled and nonlinear; therefore, they are solved by an iterative numerical scheme. GFSSP employs a unique numerical scheme known as simultaneous adjustment with successive substitution (SASS), which is a combination of Newton-Raphson and successive substitution methods. The mass and momentum conservation equations and the equation of state are solved by the Newton-Raphson method while the conservation of energy and species are solved by the successive substitution method. GFSSP is linked with two thermodynamic property programs, GASP2 and WASP3 and GASPAK4, that provide thermodynamic and thermophysical properties of selected fluids. Both programs cover a range of pressure and temperature that allows fluid properties to be evaluated for liquid, liquid-vapor (saturation), and vapor region. GASP and WASP provide properties of 12 fluids. GASPAK includes a library of 36 fluids. GFSSP has three major parts. The first part is the graphical user interface (GUI), visual thermofluid analyzer of systems and components (VTASC). VTASC allows users to create a flow circuit by a 'point and click' paradigm. It creates the GFSSP input file after the completion of the model building process. GFSSP's GUI provides the users a platform to build and run their models. It also allows post-processing of results. The network flow circuit is first built using three basic elements: boundary node, internal node, and branch.

Majumdar, A. K.↗

Inversion of the Side-Chain Stereochemistry of Indvidual Thr or Ile Residues in a Protein Molecule: Impact on the Folding, Stability, and Structure of the ShK Toxin

ShK toxin is a cysteine-rich 35-residue protein ion-channel ligand isolated from the sea anemone Stichodactyla helianthus. Here, we studied the effect of inverting the side chain stereochemistry of individual Thr or Ile residues on the properties of the ShK protein. Molecular dynamics simulations were used to calculate the free energy cost of inverting the side-chain stereochemistry of individual Thr or Ile residues. Guided by the computational results, we used chemical protein synthesis to prepare three ShK polypeptide chain analogues, each containing either an allo-Thr or an allo-Ile residue. The three allo-Thr or allo-Ile-containing ShK polypeptides were able to fold into defined protein products, but with different folding propensities. Their relative thermal stabilities were measured and were consistent with the MD simulation data. Structures of the three ShK analogue proteins were determined by quasi-racemic X-ray crystallography and were similar to wild-type ShK. All three ShK analogues retained ion-channel blocking activity.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Phase transformations in oxides above 2000°C: experimental technique development

Here, the oxidation of boride and carbide-based ultra-high-temperature ceramics is the primary limiting factor for their use as aerodynamic surfaces. Understanding the behaviour of the oxides that can result from oxidation of metal borides and carbides at very high temperatures is essential to optimise and tailor the performance of these materials; yet experimental thermodynamic and structural data for refractory oxides above 2000°C are mostly absent. The following techniques that can be applied to fill this gap are discussed: (i) commercial ultra-high-temperature differential thermal analysis for investigation of phase transformations and melting in inert environments to 2500°C, (ii) a combination of laser heating with a splittable nozzle aerodynamic levitator for splat quenching and drop calorimetry from temperatures limited only by sample evaporation, (iii) synchrotron X-ray and neutron diffraction on laser-heated aerodynamically levitated oxide samples for in situ observation of phase transformations in variable atmospheres, refinement of high-temperature structures and thermal expansion. Recent experimental findings include anomalous thermal expansion of the defect fluorite phase of YSZ, thermodynamics of pyrochlore–fluorite transformation from high-temperature structure refinements, and measurement of thermal expansion to the melting temperatures and fusion enthalpies of Zr, Hf, La, Yb and Lu oxides. These methods provide temperatures, enthalpies and volume change for phase transformations above 2000°C, which are required for thermodynamic assessments and calculation of phase diagrams of multicomponent systems.

36 MATERIALS SCIENCE↗

Trace: a high-throughput tomographic reconstruction engine for large-scale datasets

Abstract Background Modern synchrotron light sources and detectors produce data at such scale and complexity that large-scale computation is required to unleash their full power. One of the widely used imaging techniques that generates data at tens of gigabytes per second is computed tomography (CT). Although CT experiments result in rapid data generation, the analysis and reconstruction of the collected data may require hours or even days of computation time with a medium-sized workstation, which hinders the scientific progress that relies on the results of analysis. Methods We present Trace, a data-intensive computing engine that we have developed to enable high-performance implementation of iterative tomographic reconstruction algorithms for parallel computers. Trace provides fine-grained reconstruction of tomography datasets using both (thread-level) shared memory and (process-level) distributed memory parallelization. Trace utilizes a special data structure called replicated reconstruction object to maximize application performance. We also present the optimizations that we apply to the replicated reconstruction objects and evaluate them using tomography datasets collected at the Advanced Photon Source. Results Our experimental evaluations show that our optimizations and parallelization techniques can provide 158× speedup using 32 compute nodes (384 cores) over a single-core configuration and decrease the end-to-end processing time of a large sinogram (with 4501 × 1 × 22,400 dimensions) from 12.5 h to <5 min per iteration. Conclusion The proposed tomographic reconstruction engine can efficiently process large-scale tomographic data using many compute nodes and minimize reconstruction times.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

User alternatives in post-processing for raster-to-vector conversion

A number of Landsat-based coniferous forest stratum maps have been created of the Eldorado National Forest in California. These maps were produced in raster image format which is not directly usable by the U.S. Forest Service's vector-based Wildland Resource Information System (WRIS). As a solution, raster-to-vector conversion software has been developed for processing classified images into polygonal data structures. Before conversion, however, the digital classification images must be simplified to remove high spatial variance ('noise', 'speckle') and meet a USFS ten acre minimum requirement. A post-processing (simplification) strategy different from those commonly used in raster image processing may be desired for preparing maps for conversion to vector format, because simplification routines typically permit diagonal connections in the process of reclassifying pixels and forming new polygons. Diagonal connections are often undesirable when converting to vector format because they permit polygons to effectively cross over each other and occupy a common location. Three alternative methodologies are discussed for simplifying raster data for conversion to vector format.

Logan, T. L.↗

Surface Energy Dynamics and Canopy Structural Properties in Intact and Disturbed Forests in the Southern Amazon

The Brazilian Amazon has been a focus of land development with large swaths of forests converted to agriculture. Forest degradation by selective logging and fires has accompanied the agricultural frontier and has resulted in significant impacts on Amazonian ecosystems. Changes in forest structure resulting from forest disturbances have large impacts on the surface energy balance, including on land surface temperature (LST) and evapotranspiration (ET). This study's objective is to assess the effects of forest disturbances on water fluxes and forest structure in a transitional forest site in the Southern Amazon. We used ET and LST products from MODIS and Landsat 8 and GEDI-derived forest structure data to address our research questions. We found that disturbances induced seasonal water stress, more pronounced in croplands/pastures than in forests (differences up to 20% in the dry season), and more pronounced in second-growth and recently burned areas than in logged and intact forests (differences up to 12% in the dry season). Moreover, ET and LST were negatively related, with more consistent relationships across disturbance classes in the dry season (R 2 : 0.41–0.87) than in the wet season (R 2 : 0.18–0.49). Forest and cropland or pasture classes showed contrasting relationships in the dry season. Finally, we found that forest structure exhibited stronger relationships with ET and LST in the most disturbed forests (R 2 : 0.01–0.43) than in the least disturbed forests (R 2 < 0.05). Our findings help to elucidate degraded forests functioning under a changing climate and to improve estimates of water and energy fluxes in Amazonian degraded forests.

54 ENVIRONMENTAL SCIENCES↗

Influence of pump laser fluence on ultrafast myoglobin structural dynamics

High-intensity femtosecond pulses from an X-ray free-electron laser enable pump–probe experiments for the investigation of electronic and nuclear changes during light-induced reactions. On timescales ranging from femtoseconds to milliseconds and for a variety of biological systems, time-resolved serial femtosecond crystallography (TR-SFX) has provided detailed structural data for light-induced isomerization, breakage or formation of chemical bonds and electron transfer. However, all ultrafast TR-SFX studies to date have employed such high pump laser energies that nominally several photons were absorbed per chromophore. As multiphoton absorption may force the protein response into non-physiological pathways, it is of great concern whether this experimental approach allows valid conclusions to be drawn vis-à-vis biologically relevant single-photon-induced reactions. Here we describe ultrafast pump–probe SFX experiments on the photodissociation of carboxymyoglobin, showing that different pump laser fluences yield markedly different results. In particular, the dynamics of structural changes and observed indicators of the mechanistically important coherent oscillations of the Fe–CO bond distance (predicted by recent quantum wavepacket dynamics) are seen to depend strongly on pump laser energy, in line with quantum chemical analysis. Our results confirm both the feasibility and necessity of performing ultrafast TR-SFX pump–probe experiments in the linear photoexcitation regime. We consider this to be a starting point for reassessing both the design and the interpretation of ultrafast TR-SFX pump–probe experiments such that mechanistically relevant insight emerges.

59 BASIC BIOLOGICAL SCIENCES↗

Structural tunability and origin of two-level systems in amorphous silicon

Amorphous silicon films prepared by electron-beam evaporation have systematically and substantially greater atomic density for higher thickness, higher growth temperature, and slower deposition rate, reaching the density of crystalline Si when films of thickness greater than ~300 nm are grown at 425 °C and at <1 Å/s. Here, a combination of spectroscopic techniques provide insight into atomic disorder, local strains, dangling bonds, and nanovoids. Electron diffraction shows that the short-range order of the amorphous silicon is similar at all growth temperatures, but fluctuation electron microscopy shows that films grown above room temperature show a form of medium-range order not previously observed in amorphous silicon. Atomic disorder and local strain obtained from Raman spectroscopy are reduced with increasing growth temperature and show a nonmonotonic dependence on thickness. Dangling bond density decreases with increasing growth temperature and is only mildly dependent on thickness. Positron annihilation Doppler broadening spectroscopy and electron energy loss spectroscopy show that nanovoids, and not density variations within the network, are responsible for reduced atomic density. Specific heat and mechanical loss measurements, which quantify the density of tunneling two-level systems, in combination with the structural data, suggest that two-level systems in amorphous silicon films are associated with nanovoids and their surroundings, which are in essence loosely bonded regions where atoms are less constrained.

36 MATERIALS SCIENCE↗

Real-time Ridesharing for Transportation Hubs with Demand and Supply Uncertainty

Transportation hubs in major cities generate a significant amount of trips by taxis and for-hail vehicles (FHV), with many of the trips sharing similar destinations. This suggests promising opportunities to leverage the collective travel needs with dedicated ridesharing solutions to reduce the externalities of excessive traffic at transportation hubs. In this study, we develop a novel dynamic ridesharing approach to serve trips from the transportation hub by considering (1) demand (new passengers) and supply (newly available vehicles) in the near future and (2) the uncertainty of future predictions. Our approach consists of two stages. In the first stage, we develop a data structure called hub mobility tree to generate potential combinations of shareable trips as candidate schedules efficiently. Then the generated schedules are used in the second stage to formulate the stochastic hub-based ridesharing problem (SHRP), which is a stochastic integer programming problem with the objective to maximize the total expected ridesharing profit over time. Due to the prohibitive number of shareable trips, we then approximately solve SHRP by the sample average approximate method (SAA), and a dual Lagrangian technique is implemented to further improve the scalability of the solution approach. We demonstrate the performances of the proposed method by simulating the ridesharing service at JFK airport using NYC taxi and FHV data. The results indicate that the proposed method outperforms the myopic ridesharing (maximize profit for a single time step) and the rolling horizon method with point estimation of future demand and supply.

dynamic ridesharing↗

Decoding substrate specificity determining factors in glycosyltransferase-B enzymes – insights from machine learning models

Substrate specificity is an essential characteristic of any enzyme's function and an understanding of the factors that determine this specificity is crucial for enzyme engineering. Unlike the structure of an enzyme which is directly impacted by its sequence, substrate specificity as an enzyme attribute involves a rather indirect relationship with sequence as it also depends on structural aspects that dictate substrate accessibility and active site dynamics. In this study, we explore the performance of classifier-based machine learning models trained on curated sequence and structural data for a class of glycosyltransferases (GTs), namely GT-Bs, to understand their substrate specificity determining factors. GTs enable the transfer of sugar moieties to other biomolecules such as oligosaccharides or proteins and are found in all kingdoms of life. In plants, GTs participate in the biosynthesis of plant cell wall biopolymers (e.g.: hemicelluloses and pectins) and are an integral part of the enzymatic machinery that enables the storage of carbon and energy as plant biomass. To elucidate the substrate specificity of uncharacterized GT-Bs, we constructed multi-label machine learning models (Support Vector Classifier, K-Nearest Neighbors, Gaussian Naïve-Bayes, Random Forest) that incorporate both sequence and structural features. These models achieve good predictive accuracies on test datasets. However, despite our use of structural information, we highlight that there is further scope for improvement in training these models to draw interpretable relationships between sequence, structure and substrate specificity determining motifs in GT-Bs.

97 MATHEMATICS AND COMPUTING↗

Nuclear Resonance Vibrational Spectroscopic Definition of Peroxy Intermediates in Nonheme Iron Sites

Fe III -(hydro)peroxy intermediates have been isolated in two classes of mononuclear nonheme Fe enzymes that are important in bioremediation: the Rieske dioxygenases and the extradiol dioxygenases. The binding mode and protonation state of the peroxide moieties in these intermediates are not well-defined, due to a lack of vibrational structural data. Nuclear resonance vibrational spectroscopy (NRVS) is an important technique for obtaining vibrational information on these and other intermediates, as it is sensitive to all normal modes with Fe displacement. Here in this paper, we present the NRVS spectra of side-on Fe III -peroxy and end-on Fe III -hydroperoxy model complexes and assign these spectra using calibrated DFT calculations. We then use DFT calculations to define and understand the changes in the NRVS spectra that arise from protonation and from opening the Fe–O–O angle. This study identifies four spectroscopic handles that will enable definition of the binding mode and protonation state of Fe III -peroxy intermediates in mononuclear nonheme Fe enzymes. These structural differences are important in determining the frontier molecular orbitals available for reactivity.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗