Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “biochemicals”

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 739 records · Page 41
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

How “hot” are hotspots: Statistically localizing the high-activity areas on soil and rhizosphere images

The topic of microbial hotspots in soil requires not only visualizing their spatial distribution and biochemical analyses, but also statistical approaches to identify these hotspots and separate them from the surrounding activities (background). We hypothesized that each hotspot type (e.g. enzyme activities in the rhizosphere, root exudation, localization of herbicide accumulation) is a result of local process driven by biotic and/or abiotic factors, and the process rates in the hotspots are much faster than those in the soil background. We further hypothesized that the background and hotspot activities in soil belong to different statistical distributions. Consequently, hotspot determination should be based on statistical separation of activities significantly higher than the background. We analyzed for the statistical distributions of grey values on three groups of published images: 1) 14 C images of carbon input by roots into the rhizosphere, 2) 14 C glyphosate accumulation in the plant, and 3) zymogram of leucine aminopeptidase activity in rooted soil. The two Gaussian distributions were fit (the first representing the background, the second the hotspots) to the distribution of grey values in the images, the parameters (means and standard deviations, SD) of the fitted distributions were calculated, and the background was removed. Thus, we identified hotspots as areas outside of the Mean+2SD image intensity (corresponding to the upper ~ 2.5% of activity, being over 97.5% of background values) and finally, visualized images of solely hotspot locations. Finally, these results were compared with previously used decisions on hotspot intensity thresholding (i.e. Top-25% and 17 standard thresholding approaches in ImageJ) and discussed the advantages of the Mean+2SD as well as Mean+3SD approaches. These advantages include: i) simple unification of the thresholding approach for several imaging methods with various principles of activity distribution, ii) identification of hotspots with various activity levels, iii) analysis of “time-specific” hotspots in temporal sequences of images. Compared with 17 standard thresholding methods, we concluded that objectively elucidating and separating the hotspots should be based on statistical distribution analysis, e.g. using the Mean+2SD or Mean+3SD approaches. Furthermore, this simple Mean+2SD approach delivered suitable results for three groups of images and so, helps to understand the processes responsible for the highest activities and elucidate hotspots.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

Characterizing leaf-scale fluorescence with spectral invariants

Sun-induced chlorophyll fluorescence (SIF) is increasingly recognized as a non-destructive probe for tracking terrestrial photosynthesis. Emerging developments in spectral invariants theory provide an innovative and efficient approach for representing SIF radiative transfer processes at the canopy scale. However, modeling leaf-scale fluorescence based on the spectral invariants properties (SIP) remains underexplored. In this study, the spectral invariants theory is employed for the first time to model the leaf-scale total, backward and forward fluorescence (leaf-SIP SIF). The leaf-SIP SIF model separates the leaf-scale radiative transfer process into two distinct components: the wavelength-dependent one associated with leaf biochemical properties, and the wavelength-independent component linked to leaf structural characteristics. The leaf structure-related effects are characterized by two spectrally invariant parameters: the photon recollision probability (p) and the scattering asymmetry parameter (q), which are parameterized using the directly measurable leaf dry matter. Evaluation against field measurements shows that the proposed leaf-SIP SIF model has a good performance, with coefficient of determination (R 2 ) of 0.89, 0.89, 0.90 and root mean squared errors (RMSE) of 1.28, 0.69, 0.74 Wm -2 µm -1 sr -1 , respectively for the total, backward, and forward fluorescence (660–800 nm). The leaf-SIP SIF model with a more concise formulation demonstrates comparable performance with the widely used Fluspect model. Further, the leaf-SIP SIF model provides a simple and efficient approach for simulating leaf-scale fluorescence, with the potential to be integrated into a unified SIP-based model framework for simulating the radiative transfer processes across the soil-leaf-canopy-atmosphere continuum.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

Combining geometric-optical and spectral invariants theories for modeling canopy fluorescence anisotropy

The spectral invariants theory ( p -theory) has received much attention in the field of quantitative remote sensing over the past few decades and has been adopted for modeling of canopy solar-induced chlorophyll fluorescence (SIF). However, the spectral invariant properties (SIP) in simple analytical formulae have not been applied for modeling canopy fluorescence anisotropy primarily because they are parameterized in terms of leaf total scattering, which precludes the differentiation between forward and backward leaf SIF emissions. In this study, we have developed the canopy-SIP SIF model by combining geometric-optical (GO) theory to account for asymmetric leaf SIF forward and backward emissions at the first-order scattering and by modeling multiple scattering based on the p-theory, thus avoiding the dependence on radiative transfer models. The applicability of the model simulations especially over 3D heterogeneous canopies was improved by incorporating canopy structure through multi-angular clumping index, and by modeling single scattering from the four components of the scene in view according to the GO approach. The results show good consistency with both the state-of-the-art SIF models and multi-angular field SIF observations over grass and chickpea canopies. Further, the coefficient of determination (R²) between the simulated SIF and field measurements was 0.75 (red) and 0.74 (far-red) for chickpea, and 0.65 (both red and far-red) for grass. The average relative error was approximately 3% for 1D homogeneous scenes when comparing the canopy-SIP SIF model simulations to the SCOPE model simulations, and around 4% for the 3D heterogeneous scene when comparing to the LESS model simulations. The results indicate that the proposed approach for separating asymmetric leaf SIF emissions is a robust way to keep a balance between satisfactory simulation accuracy and efficiency. Model simulations suggest that neglecting the leaf SIF asymmetry can lead to an underestimation of canopy red SIF by 6.3% to 42.6% for various leaf biochemical and canopy structural parameters. This study presents a simple but efficient analytical approach for canopy fluorescence modeling, with potential for large-scale canopy fluorescence simulations.

3D heterogeneous structure↗
DOE OSTI2025

RSPECT: A PROSPECT-based model incorporating the real structure of rice leaves

Radiative transfer models (RTMs) describe how light is absorbed, scattered, and transmitted within leaves by simulating mechanistic light propagation processes. The PROSPECT model is based on measurable parameters (the leaf biochemical content) and a non-measurable parameter (the leaf anatomical structure represented by the leaf structure parameter (N)). The effect of N on the optical properties of leaves has been investigated through a number of local and global sensitivity analyses. Other studies have directly evaluated the effect of the leaf anatomical structure on spectral reflectance, particularly in the near infrared region. However, the relationship between N and the anatomical structure is unclear. Here, in this study, we leveraged eLeaf, a ray tracing-based 3D rice leaf simulator, to establish relationships between leaf anatomical features and spectral properties, enabling us to replace N in the PROSPECT-4 model with measurable leaf anatomical parameters and develop the RSPECT model. The leaf thickness at minor vein, leaf thickness at bulliform cells, mesophyll thickness at minor vein, and distance between two minor veins could be used to predict N effectively. The RSPECT model achieved spectral simulation accuracy comparable to PROSPECT-4 and was more suitable for parameter inversion of the physical and chemical properties of rice leaves, with relative root mean square errors of 7.4% for chlorophyll content, 5.6% for equivalent water thickness, and 7.5% for dry matter content. In conclusion, RSPECT improves radiative transfer modeling by integrating measurable anatomical features and provides a framework for extending this approach to other vegetation types.

hyperspectral↗
DOE OSTI2023

The significance of biomass densification in biological-based biorefineries: A critical review

Replacing fossil fuels with renewable biofuels derived from lignocellulosic biomass is an important aspect of addressing environmental challenges and developing a sustainable industrial society. Densification overcomes the problems of low bulk density and poor flowability of biomass and has been commercialized for producing solid biofuels, but the overall impact of densification on biological-based biorefineries that primarily target liquid fuels (e.g., ethanol) is still under investigation. Herein, this review provides a thorough summary of the application of densified biomass in biological-based biorefineries. First, the effects of densification parameters as well as variables of biomass materials on the densified products’ quality are reviewed, and different commonly used densification technologies are also discussed and compared. Then, the discussion focuses on the physiochemical modifications of biomass caused by densification that may influence further pretreatment and/or enzymatic hydrolysis in biorefineries. Industrial pelleting has been generally shown to exhibit positive/neutral effects on the enzymatic hydrolysis of multiple biomass feedstocks after pretreatments, indicating the viability of using pellets as starting feedstocks in biorefineries. Densification causes structural disruption of biomass, which may facilitate further biochemical conversions. Integrating biomass densification in the feedstock supply chain is feasible for large-scale biorefineries to overcome the techno-economic barriers and become profitable. Suggestions are presented for the efficiency enhancement and cost reduction in densification-based biorefineries. To enable wider applications, it is now the time to employ more demo and full-scale activities in different regions of the world pushing the research and innovation of densification integrated within biorefining.

09 BIOMASS FUELS↗
DOE OSTI2020

Selenium tolerance, accumulation, localization and speciation in a Cardamine hyperaccumulator and a non-hyperaccumulator

Cardamine violifolia (family Brassicaceae ) is the first discovered selenium hyperaccumulator from the genus Cardamine with unique properties in terms of selenium accumulation, i.e., high abundance of selenolanthionine. In our study, a fully comprehensive experiment was conducted with the comparison of a non-hyperaccumulator Cardamine species, Cardamine pratensis , covering growth characteristics, chlorophyll fluorescence, spatial selenium/sulfur distribution patterns through elemental analyses (synchrotron-based X-Ray Fluorescence and ICP-OES) and speciation data through selenium K-edge micro X-ray absorption near-edge structure analysis (μXANES) and strong cation exchange (SCX)-ICP-MS. The results revealed remarkable differences in contrast to other selenium hyperaccumulators as neither Cardamine pecies showed evidence of growth stimulation by selenium. Also, selenite uptake was not inhibited by phosphate for either of the Cardamine species. Sulfate inhibited selenate uptake, but the two Cardamine species did not show any difference in this respect. However, μXRF derived speciation maps and selenium/sulfur uptake characteristics provided results that are similar to other formerly reported hyperaccumulator and non-hyperaccumulator Brassicaceae species. μXANES showed organic selenium, “C-Se-C”, in seedlings of both species and also in mature C. violifolia plants. In contrast, selenate-supplied mature C. pratensis contained approximately half “C-Se-C” and half selenate. SCX-ICP-MS data showed evidence of the lack of selenocystine in any of the Cardamine plant extracts. Thus, C. violifolia shows clear selenium-related physiological and biochemical differences compared to C. pratensis and other selenium hyperaccumulators.

Environmental Sciences & Ecology↗
DOE OSTI2021

Exposure, health effects, sensing, and remediation of the emerging PFAS contaminants – Scientific challenges and potential research directions

Per- and polyfluoroalkyl substances (PFAS) make up a large group of persistent anthropogenic chemicals which are difficult to degrade and/or destroy. PFAS are an emerging class of contaminants, but little is known about the long-term health effects related to exposure. In addition, technologies to identify levels of contamination in the environment and to remediate contaminated sites are currently inadequate. In this opinion-type discussion paper, a team of researchers from the University of Connecticut and the University at Albany discuss the scientific challenges in their specific but intertwined PFAS research areas, including rapid and low-cost detection, energy-saving remediation, the role of T helper cells in immunotoxicity, and the biochemical and molecular effects of PFAS among community residents with measurable PFAS concentrations. Potential research directions that may be employed to address those challenges and improve the understanding of sensing, remediation, exposure to, and health effects of PFAS are then presented. Here, we hope our account of emerging problems related to PFAS contamination will encourage a broad range of scientific experts to bring these research initiatives addressing PFAS into play on a national scale.

Challenges and research directions↗
DOE OSTI2021

Ecological theory applied to environmental metabolomes reveals compositional divergence despite conserved molecular properties

Stream and river systems transport and process substantial amounts of dissolved organic matter (DOM) from terrestrial and aquatic sources to the ocean, with global biogeochemical implications. However, the underlying mechanisms affecting the spatiotemporal organization of DOM composition are under-investigated. To understand the principles governing DOM composition, we leverage the recently proposed synthesis of metacommunity ecology and metabolomics, termed ‘meta-metabolome ecology.’ Applying this novel approach to a freshwater ecosystem, we demonstrated that despite similar molecular properties across metabolomes, metabolite identity significantly diverged due to environmental filtering and variations in putative biochemical transformations. We refer to this phenomenon as ‘thermodynamic redundancy,’ which is analogous to the ecological concept of functional redundancy. We suggest that under thermodynamic redundancy, divergent metabolomes can support equivalent biogeochemical function just as divergent ecological communities can support equivalent ecosystem function. As these analyses are performed in additional ecosystems, potentially generalizable principles, like thermodynamic redundancy, can be revealed and provide insight into DOM dynamics.

Danczak, Robert E.↗
DOE OSTI2021

Performance characterization of nanofiltration, reverse osmosis, and ion exchange technologies for acetic acid separation

A major obstacle to widespread implementation of bio-based fuels and chemicals is the lack of efficient and costeffective separation methods. To purify acetic acid produced by biochemical conversion of biomass via anaerobic digestion, this work employs two commonly used separation technologies: (1) ion-exchange (IX) resin and (2) pressure-driven membranes. This study tested five commercially available strong- and weak-base anion-exchange resins and five commercially available nanofiltration (NF) and reverse osmosis (RO) membranes. The pH of the feed solution significantly affected adsorption capacity. At pH 6.3, a strong-base IX resin (IRN-78) performed best (95.1% acetate removal). With strong-base IX resins, the Langmuir isotherm model fit well, whereas for weak-base IX resins, the Freundlich isotherm provided a better fit. A pseudo-second-order kinetic model fit well for both IRN-78 and IRA-67. Regarding membrane separation, RO (BW30XFR membrane) achieved the highest rejection (98.6% acetate rejection), whereas an NF membrane (NF*) achieved the best combination of permeate flux (105 L/(h∙m 2 )) and rejection (83.1% acetate rejection). For membrane performance, the experimental data were fit using the solution diffusion model. Increased pH in the feed solution lowered permeate flux but increased acetic acid rejection. Finally, then the acetic acid concentration in the feed solution increased, both permeate flux and acetic acid rejection decreased for membrane NF*.

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

Antecedent conditions determine the biogeochemical response of coastal soils to seawater exposure

Coastal landscapes are increasingly exposed to seawater due to sea level rise and extreme weather events. The biogeochemical responses of these vulnerable ecosystems are poorly understood, limiting our ability to predict how their role in local and global biogeochemical cycles will shift under future conditions. Here we evaluate how antecedent conditions influence the biogeochemical response of soil to seawater inundation events based on a 42-day laboratory incubation experiment with soils collected from a natural salinity gradient across a coastal floodplain. We quantified influences of seawater inundation on intact soil cores through high-frequency carbon dioxide (CO 2 ) and methane (CH 4 ) gas fluxes measurements as well as ultrahigh resolution characterization of organic matter chemistry and metabolites. Mean CO 2 and CH 4 fluxes were higher after inundation compared to control cores for soils that had low in situ salinity. Soils with low in situ salinity also exhibited significant shifts in organic matter profiles after inundation, with surficial soils (0-7.5 cm) becoming more enriched in phenolic compounds, compared to deeper soils (7.5-15 cm). The number of biochemical transformations inferred from mass spectrometry increased significantly after inundation for soils with low in situ salinity, and a trend of decreasing thermodynamic favorability with increasing salinity was observed, suggesting an increase in microbial activity following inundation. Our results suggest that seawater inundation of low-salinity terrestrial environments can lead to increased microbial activity, preferential depletion of more thermodynamically favorable compounds, and increased soil carbon release. Finally, we conclude that the biogeochemical impacts of future seawater exposure will be modulated by antecedent conditions associated with landscape position within coastal watersheds.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Spatial access and resource limitations control carbon mineralization in soils

Core-scale soil carbon fluxes are ultimately regulated by pore-scale dynamics of substrate availability and microbial access. These are constrained by physicochemical and biochemical phenomena (e.g. spatial access and hydrologic connectivity, physical occlusion, adsorption-desorption with mineral surfaces, nutrient and resource limitations). We conducted an experiment to determine how spatial access and resource limitations influence core-scale SOM mineralization, and how these are regulated by antecedent moisture conditions. Intact soil cores were incubated at field-moist vs. drought conditions, after which they were saturated from above (to simulate precipitation) or below (to simulate groundwater recharge). Soluble C (acetate) and N (nitrate) forms were added to some cores during the rewetting process to alleviate potential nutrient limitations. Respiration showed evidence of C-limitation, as C amendments increased the cumulative CO 2 evolved. Drought and rewetting increased soil respiration, and these soils also exhibited increased complex aromatic molecules in porewater. This newly available substrate appeared to alleviate nutrient limitations on respiration, because there were no further respiration increases with subsequent C and N amendments. We expected that respiration would be influenced by wetting direction, as simulated precipitation would mobilize C from the surface. However, as a main effect, this response was seen only in the C-amended soils, indicating that surface-C may not have been bioavailable. At the pore scale (porewater samples), compared to the baseline soils, drought and the C, N amendments caused a net loss of identified molecules when the soils were rewet from below, whereas wetting from above caused a net increase in identified molecules. This indicates that as soils were wet from below, the fresh inputs simulated the C-and N-limited microbial populations present deeper in the soil profile. Furthermore, our experiment highlights the complex and interactive role of antecedent moisture conditions, wetting direction, and resource limitations in driving core-scale C fluxes.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Revising the dynamic energy budget theory with a new reserve mobilization rule and three example applications to bacterial growth

Dynamic energy budget (DEB) theory has been applied to model a wide range of organisms, including microbes. In the standard DEB model, biomass is partitioned into reserve and structural compartments, where reserve biomass is mobilized in a pseudolinear manner (while the reserve biomass density, defined as the ratio between reserve and structural biomass, decays linearly) to drive maintenance and the growth of structural biomass (and extracellular enzyme production if it is considered). However, the linear dynamics of the reserve biomass density makes the standard DEB model incapable of explaining the slowdown of microbial growth at high reserve density that is caused by macromolecular crowding effect which reduces biochemical reaction rates (a typical situation occurs when microbes are experiencing severe moisture stress) and is inconsistent with the observation that intracellular enzymatic reactions generally follow non-linear kinetics. By partitioning biomass into reserve, kinetic, and structural compartments, we show here that the Equilibrium Chemistry Approximation (ECA) kinetics can be used to represent enzymatically catalyzed reserve biomass mobilization that can then drive the kinetic and structural biomass synthesis. This revised DEB model better represents the tradeoff in ribosome allocation for structural growth and internal enzyme production, is structurally compatible with metabolic models of cell individuals, and includes the standard DEB model and the popular compromise model as special cases for representing population growth. We then applied the revised DEB model to interpret components of bacterial respiration, their dependence on substrate availability, and emergent microbial carbon use efficiency dynamics for an exponentially growing population. We found that the revised DEB model enables a better understanding of bacterial substrates use (carbon in our examples) than that can be derived from a few other models in the literature. In particular, the revised DEB model explains why carbon use efficiency may first increase, then plateau, and finally decrease with growth rate (and substrate uptake rate), as a function of proteomics. Additionally, the revised DEB model explains why the kinetic biomass compartment needs to be divided to reasonably incorporate proteomic control of microbial growth.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Artificial extracellular matrix scaffolds of mobile molecules enhance maturation of human stem cell-derived neurons

Human induced pluripotent stem cell (hiPSC) technologies offer a unique resource for modeling neurological diseases. However, iPSC models are fraught with technical limitations including abnormal aggregation and inefficient maturation of differentiated neurons. These problems are in part due to the absence of synergistic cues of the native extracellular matrix (ECM). We report on the use of three artificial ECMs based on peptide amphiphile (PA) supramolecular nanofibers. All nanofibers display the laminin-derived IKVAV signal on their surface but differ in the nature of their non-bioactive domains. Here, we find that nanofibers with greater intensity of internal supramolecular motion have enhanced bioactivity toward hiPSC-derived motor and cortical neurons. Proteomic, biochemical, and functional assays reveal that highly mobile PA scaffolds caused enhanced β1-integrin pathway activation, reduced aggregation, increased arborization, and matured electrophysiological activity of neurons. Our work highlights the importance of designing biomimetic ECMs to study the development, function, and dysfunction of human neurons.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2022

Soft X-ray tomography to map and quantify organelle interactions at the mesoscale

Inter-organelle interactions are a vital part of normal cellular function; however, these have proven difficult to quantify due to the range of scales encountered in cell biology and the throughput limitations of traditional imaging approaches. Here, we demonstrate that soft X-ray tomography (SXT) can be used to rapidly map ultrastructural reorganization and inter-organelle interactions in intact cells. SXT takes advantage of the naturally occurring, differential X-ray absorption of the carbon-rich compounds in each organelle. Specifically, we use SXT to map the spatiotemporal evolution of insulin vesicles and their co-localization and interaction with mitochondria in pancreatic β cells during insulin secretion and in response to different stimuli. We quantify changes in the morphology, biochemical composition, and relative position of mitochondria and insulin vesicles. These findings highlight the importance of a comprehensive and unbiased mapping at the mesoscale to characterize cell reorganization that would be difficult to detect with other existing methodologies.

3D cell mapping↗
DOE OSTI2023

Epitope topography of agonist antibodies to the checkpoint inhibitory receptor BTLA

B and T lymphocyte attenuator (BTLA) is an attractive target for a new class of therapeutics that attempt to rebalance the immune system by agonizing checkpoint inhibitory receptors (CIRs). Herpesvirus entry mediator (HVEM) binds BTLA in both trans- and cis-orientations. We report here the development and structural characterization of three humanized BTLA agonist antibodies, 22B3, 25F7, and 23C8. We determined the crystal structures of the antibody-BTLA complexes, showing that these antibodies bind distinct and non-overlapping epitopes of BTLA. While all three antibodies activate BTLA, 22B3 mimics HVEM binding to BTLA and shows the strongest agonistic activity in functional cell assays and in an imiquimod-induced mouse model of psoriasis. 22B3 is also capable of modulating HVEM signaling through the BTLA-HVEM cis-interaction. The data obtained from crystal structures, biochemical assays, and functional studies provide a mechanistic model of HVEM and BTLA organization on the cell surface and informed the discovery of a highly active BTLA agonist.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2026

Classifying biophysical subpopulations of insulin secretory granules using quantitative whole-cell structure analysis

Pancreatic beta cells contain insulin secretory granules (ISGs), organelles where proinsulin is converted into insulin. As ISGs mature, they undergo extensive biophysical remodeling, producing a spectrum of subpopulations with heterogeneous molecular and spatial characteristics. However, systematic methods to define ISG subpopulations remain underdeveloped. To address this gap in knowledge, we employed soft X-ray tomography (SXT), which can quantitatively measure the biochemical density of ISGs within whole beta cells. Using unsupervised clustering, we classified subpopulations based on molecular density, size, and spatial positioning. Across different insulin secretory stimuli, we observed shifts toward mature and releasable subtypes, demonstrating that exogenous signals can dynamically remodel ISG subpopulation distributions. We extended this methodology to primary beta cells characterized using volume electron microscopy (vEM). Integrating subpopulations from SXT and vEM uncovered insights inaccessible by a single method in isolation. This strategy establishes a framework for defining therapeutic approaches aimed at enriching physiologically beneficial ISG subpopulations.

dense-core granules↗
DOE OSTI2021

Plant biochemistry influences tropospheric ozone formation, destruction, deposition, and response

Tropospheric ozone (O 3 ) is among the most damaging air pollutant to plants. Plants alter the atmospheric O 3 concentration in two distinct ways: (i) by the emission of volatile organic compounds (VOCs) that are precursors of O 3 ; and (ii) by dry deposition, which includes diffusion of O 3 into vegetation through stomata and destruction by nonstomatal pathways. Isoprene, monoterpenes, and higher terpenoids are emitted by plants in quantities that alter tropospheric O 3 . Deposition of O 3 into vegetation is related to stomatal conductance, leaf structural traits, and the detoxification capacity of the apoplast. The biochemical fate of O 3 once it enters leaves and reacts with aqueous surfaces is largely unknown, but new techniques for the tracking and identification of initial products have the potential to open the black box.

59 BASIC BIOLOGICAL SCIENCES↗