Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “biophysical models”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Competition between inside-out unfolding and pathogenic aggregation in an amyloid-forming β-propeller

Studies of folded-to-misfolded transitions using model protein systems reveal a range of unfolding needed for exposure of amyloid-prone regions for subsequent fibrillization. Here, we probe the relationship between unfolding and aggregation for glaucoma-associated myocilin. Mutations within the olfactomedin domain of myocilin (OLF) cause a gain-of-function, namely cytotoxic intracellular aggregation, which hastens disease progression. Aggregation by wild-type OLF (OLF WT ) competes with its chemical unfolding, but only below the threshold where OLF loses tertiary structure. Representative moderate (OLF D380A ) and severe (OLF I499F ) disease variants aggregate differently, with rates comparable to OLF WT in initial stages of unfolding, and variants adopt distinct partially folded structures seen along the OLF WT urea-unfolding pathway. Whether initiated with mutation or chemical perturbation, unfolding propagates outward to the propeller surface. In sum, for this large protein prone to amyloid formation, the requirement for a conformational change to promote amyloid fibrillization leads to direct competition between unfolding and aggregation.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

3D-Scaffold: A Deep Learning Framework to Generate 3D Coordinates of Drug-like Molecules with Desired Scaffolds

The prerequisite of therapeutic drug design is to identify novel molecules with desired biophysical and biochemical properties. Deep generative models have demonstrated their ability to find such molecules by exploring a huge chemical space efficiently. An effective way to obtain molecules with desired target properties is the preservation of critical scaffolds in the generation process. To this end, we propose a domain aware generative framework called 3D-Scaffold that takes 3D coordinates of a desired scaffold as an input and generates 3D coordinates of novel therapeutic candidates as an output while always preserving the desired scaffolds in generated structures. We show that our framework generates predominantly valid, unique, novel, and experimentally synthesizable molecules that have drug-like properties similar to the molecules in the training set. Using domain specific datasets, we generate covalent and non-covalent antiviral inhibitors. Therefore, to measure the success of our framework in generating therapeutic candidates, generated structures were subjected to high throughput virtual screening via docking simulations, which shows favorable interaction against SARS-CoV-2 main protease and non-structural protein endoribonuclease (NSP15) targets. Most importantly, our model performs well with relatively small volumes of training data and generalizes to new scaffolds, making it applicable to other domain.

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

Monatomic ions influence substrate permeation across bacterial microcompartment shells

Abstract Bacterial microcompartments (BMCs) are protein organelles consisting of an inner enzymatic core encased within a selectively permeable shell. BMC shells are modular, tractable architectures that can be repurposed with new interior enzymes for biomanufacturing purposes. The permeability of BMC shells is function-specific and regulated by biophysical properties of the shell subunits, especially its pores. We hypothesized that ions may interact with pore residues in a manner that influences the substrate permeation process. In vitro activity comparisons between native and broken BMCs demonstrated that increasing NaCl negatively affects permeation rates. Molecular dynamics simulations of the dominant shell protein (BMC-H) revealed that chloride ions preferentially occupy the positive pore, hindering substrate permeation, while sodium cations remain excluded. Overall, these results demonstrate that shell properties influence ion permeability and leverages the integration of experimental and computational techniques to improve our understanding of BMC shells towards their repurposing for biotechnological applications.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Toward designing effective exascale scientific computing workflows: experiences and best practices

Many fields within scientific computing have embraced advances in big-data analysis and machine learning, which often requires the deployment of large, distributed and complicated workflows that may combine training neural networks, performing simulations, running inference, and performing database queries and data analysis in asynchronous, parallel and pipelined execution frameworks. Such a shift has brought into focus the need for scalable, efficient workflow management solutions with reproducibility, error and provenance handling, traceability, and checkpoint-restart capabilities, among other needs. Here, we discuss challenges and best-practices for deploying exascale-generation computational science workflows on resources at the Oak Ridge Leadership Computing Facility (OLCF). We present our experiences with large-scale deployment of distributed workflows on the Summit supercomputer, including for bioinformatics and computational biophysics, materials science, and deep learning model optimization. We also present problems and solutions created by working within a Python-centric software base on traditional HPC systems, and discuss steps that will be required before the convergence of HPC, AI, and data science can be fully realized. Our results point to a wealth of exciting new possibilities for harnessing this convergence to tackle new scientific challenges.

Coletti, Mark↗
DOE OSTI2023

Biophysical Impact of Land-Use and Land-Cover Change on Subgrid Temperature in CMIP6 Models

In this study, we investigate the air temperature response to land-use and land-cover change (LULCC; cropland expansion and deforestation) using subgrid land model output generated by a set of CMIP6 model simulations. Our study is motivated by the fact that ongoing land-use activities are occurring at local scales, typically significantly smaller than the resolvable scale of a grid cell in Earth system models. It aims to explore the potential for a multimodel approach to better characterize LULCC local climatic effects. On an annual scale, the CMIP6 models are in general agreement that croplands are warmer than primary and secondary land (psl; mainly forests, grasslands, and bare ground) in the tropics and cooler in the mid–high latitudes, except for one model. The transition from warming to cooling occurs at approximately 40°N. Although the surface heating potential, which combines albedo and latent heat flux effects, can explain reasonably well the zonal mean latitudinal subgrid temperature variations between crop and psl tiles in the historical simulations, it does not provide a good prediction on subgrid temperature for other land tile configurations (crop vs forest; grass vs forest) under Shared Socioeconomic Pathway 5–8.5 (SSP5–8.5) forcing scenarios. A subset of simulations with the CESM2 model reveals that latitudinal subgrid temperature variation is positively related to variation in net surface shortwave radiation and negatively related to variation in the surface energy redistribution factor, with a dominant role from the latter south of 30°N. We suggest that this emergent relationship can be used to benchmark the performance of land surface parameterizations and for prediction of local temperature response to LULCC.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Evaluating post-fire watershed response to varying burn severity and precipitation regimes using fully-distributed and integrated hydrologic models

Wildfires can cause significant changes in vegetation and soil, which may lead to increased surface runoff and soil erosion, thereby affecting water cycling within ecosystems. This study uses the Advanced Terrestrial Simulator (ATS), an integrated and fully distributed hydrologic model at the watershed scale, to examine post-fire hydrologic responses in selected watersheds with varying burn severities in the Pacific Northwest region of the United States. The model integrates surface overland flow, subsurface flow, and canopy biophysical processes. We have developed a new fire module in ATS to account for changes in soil hydraulic properties caused by fire in the topsoil layer. Modeling results indicate that, in the year following a high-severity burn, watershed-averaged evapotranspiration decreases by 25%. Additionally, post-fire peak flows increase by 18-29% in watersheds burned with medium to high severity due to changes in soil properties. Conversely, a low-severity burn results in less than a 1% increase in post-fire peak flow. Furthermore, a high-severity fire causes a 38% reduction in the infiltration rate within the affected watershed during the first post-fire wet season. Hypothetical numerical experiments with varying precipitation regimes after a high-severity fire show that post-fire peak flows can increase by 1-29% due to fire-induced changes in soil hydraulic properties. This study highlights the importance of using fully distributed hydrologic models to quantify disturbance-feedback loops, which are essential for understanding the complexities brought about by spatial heterogeneity in post-fire landscapes.

ATS↗
DOE OSTI2024

Automated AI-driven Molecular Design for Therapeutic Discovery

In recent years, artificial intelligence and machine learning (AI/ML) approaches have revolutionized the process of designing new therapeutics, enabling scientists to rapidly respond to emerging threats from various pathogens. A prime example is the SARS-CoV-2 main protease, a key target for the development of antiviral inhibitors. In this study, we employed a novel, integrated approach that combines AI-driven iterative design of inhibitor candidates, screening based on physio-chemical properties and toxicity, physics-based computational modeling of protein-inhibitor interactions, and AI-assisted analysis of Native MS biophysical assay and characterization of designed candidates. Our deep learning 3D-scaffold model, which uses an input scaffold as a starting point, generated tens of thousands of compounds while preserving the key scaffold. To optimize these candidates, we calculated a comprehensive set of 136 descriptors, including both 2D and 3D molecular features, for compounds targeting the SARS-CoV-2 Main protease (Mpro) and a neurodegenerative disease-associated protein, cyclophilin (Cyp). The generated compounds were initially filtered based on their properties and then ranked according to their predicted binding affinity using our automated modeling and ML methods. Experimental validation of the Mpro candidates showing inhibitory activity demonstrates that our workflow can expedite the therapeutic discovery.

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

Validation of the Community Land Model Version 5 over the Contiguous United States (CONUS) using in situ and remote sensing data sets

The Community Land Model (CLM) is an effective tool to simulate the biophysical and biogeochemical processes and their interactions with the atmosphere. Although CLM Version 5 (CLM5) constitutes various updates in these processes, its performance in simulating energy, water and carbon cycles over the Contiguous United States (CONUS) at scales which land surface changes and hydrometeorological and hydroclimatological applications are more locally relevant is yet to be assessed. In this study, we conducted three simulations at 0.125? during 1979-2018 over the CONUS using different configurations of CLM, namely CLM5-biogeochemistry (CLM5BGC), CLM4.5BGC, and CLM5-satellite phenology (CLM5SP). We validated and compared their simulations against multiple remote-sensed and in-situ datasets. Overall, the parametric and structural updates (e.g., carbon cost for nitrogen uptake, variable soil thickness, dry surface layer) in CLM5 improve its ability in capturing terrestrial biogeochemical dynamics. The low evapotranspiration in CLM5BGC is associated with biases in simulating vegetation phenological characteristics rather than soil water limitations. The mismatch between CLM5BGC-simulated peak leaf area index and reference data can be attributed to CLM5BGC's inability in simulating phenology of trees and grasses. The differences between CLM-simulated irrigation and reference estimates can be attributed to differences between processes represented in models and in reality, and uncertainties in input and validation datasets. Evaluation against observations at small catchments suggest that hydrologic parameters needed to be calibrated to improve simulations of runoff, especially subsurface runoff. Additional efforts are needed to incorporate spatially-distributed plant phenology and physiology parameters and regional-specific agricultural management practices (e.g., planting, harvest).

Cheng, Yanyan↗
DOE OSTI2020

Uncovering a membrane-distal conformation of KRAS available to recruit RAF to the plasma membrane

The small GTPase KRAS is localized at the plasma membrane where it functions as a molecular switch, coupling extracellular growth factor stimulation to intracellular signaling networks. In this process, KRAS recruits effectors, such as RAF kinase, to the plasma membrane where they are activated by a series of complex molecular steps. Defining the membrane-bound state of KRAS is fundamental to understanding the activation of RAF kinase and in evaluating novel therapeutic opportunities for the inhibition of oncogenic KRAS-mediated signaling. We combined multiple biophysical measurements and computational methodologies to generate a consensus model for authentically processed, membrane-anchored KRAS. In contrast to the two membrane-proximal conformations previously reported, we identify a third significantly populated state using a combination of neutron reflectivity, fast photochemical oxidation of proteins (FPOP), and NMR. In this highly populated state, which we refer to as “membrane-distal” and estimate to comprise ~90% of the ensemble, the G-domain does not directly contact the membrane but is tethered via its C-terminal hypervariable region and carboxymethylated farnesyl moiety, as shown by FPOP. Subsequent interaction of the RAF1 RAS binding domain with KRAS does not significantly change G-domain configurations on the membrane but affects their relative populations. Overall, our results are consistent with a directional fly-casting mechanism for KRAS, in which the membrane-distal state of the G-domain can effectively recruit RAF kinase from the cytoplasm for activation at the membrane.

KRAS↗
DOE OSTI2022

Exploring the Role of Cryptic Nitrogen Fixers in Terrestrial Ecosystems: A Frontier in Nitrogen Cycling Research

Biological nitrogen fixation represents the largest natural flux of new nitrogen (N) into terrestrial ecosystems, providing a critical N source to support net primary productivity of both natural and agricultural systems. When they are common, symbiotic associations between plants and bacteria can add more than 100 kg N ha –1 y –1 to ecosystems. Yet, these associations are uncommon in many terrestrial ecosystems. In most cases, N inputs derive from more cryptic sources, including mutualistic and/or free-living microorganisms in soil, plant litter, decomposing roots and wood, lichens, insects, and mosses, among others. Unfortunately, large gaps remain in the understanding of cryptic N fixation. Here, we conducted a literature review to explore rates, patterns, and controls of cryptic N fixation in both unmanaged and agricultural ecosystems. Our analysis indicates that, as is common with N fixation, rates are highly variable across most cryptic niches, with N inputs in any particular cryptic niche ranging from near zero to more than 20 kg ha –1 y –1 . Such large variation underscores the need for more comprehensive measurements of N fixation by organisms not in symbiotic relationships with vascular plants in terrestrial ecosystems, as well as identifying the factors that govern cryptic N fixation rates. We highlight several challenges, opportunities, and priorities in this important research area, and we propose a conceptual model that posits an interacting hierarchy of biophysical and biogeochemical controls over N fixation that should generate valuable new hypotheses and research.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Assembly of JAZ–JAZ and JAZ–NINJA complexes in jasmonate signaling

Jasmonates (JAs) are plant hormones with crucial roles in development and stress resilience. They activate MYC transcription factors by mediating the proteolysis of MYC inhibitors called JAZ proteins. In the absence of JA, JAZ proteins bind and inhibit MYC through the assembly of MYC–JAZ–Novel Interactor of JAZ (NINJA)–TPL repressor complexes. However, JAZ and NINJA are predicted to be largely intrinsically unstructured, which has precluded their experimental structure determination. Through a combination of biochemical, mutational, and biophysical analyses and AlphaFold-derived ColabFold modeling, we characterized JAZ–JAZ and JAZ–NINJA interactions and generated models with detailed, high-confidence domain interfaces. We demonstrate that JAZ, NINJA, and MYC interface domains are dynamic in isolation and become stabilized in a stepwise order upon complex assembly. By contrast, most JAZ and NINJA regions outside of the interfaces remain highly dynamic and cannot be modeled in a single conformation. Our data indicate that the small JAZ Zinc finger expressed in Inflorescence Meristem (ZIM) motif mediates JAZ–JAZ and JAZ–NINJA interactions through separate surfaces, and our data further suggest that NINJA modulates JAZ dimerization. This study advances our understanding of JA signaling by providing insights into the dynamics, interactions, and structure of the JAZ–NINJA core of the JA repressor complex.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Structural basis for cloverleaf RNA-initiated viral genome replication

The genomes of positive-strand RNA viruses serve as a template for both protein translation and genome replication. In enteroviruses, a cloverleaf RNA structure at the 5' end of the genome functions as a switch to transition from viral translation to replication by interacting with host poly(C)-binding protein 2 (PCBP2) and the viral 3CD pro protein. We determined the structures of cloverleaf RNA from coxsackievirus and poliovirus. Cloverleaf RNA folds into an H-type four-way junction and is stabilized by a unique adenosine-cytidine-uridine (A•C-U) base triple involving the conserved pyrimidine mismatch region. The two PCBP2 binding sites are spatially proximal and are located on the opposite end from the 3CD pro binding site on cloverleaf. We determined that the A•C-U base triple restricts the flexibility of the cloverleaf stem–loops resulting in partial occlusion of the PCBP2 binding site, and elimination of the A•C-U base triple increases the binding affinity of PCBP2 to the cloverleaf RNA. Based on the cloverleaf structures and biophysical assays, we propose a new mechanistic model by which enteroviruses use the cloverleaf structure as a molecular switch to transition from viral protein translation to genome replication.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Employing artificial intelligence to steer exascale workflows with colmena

Computational workflows are a common class of application on supercomputers, yet the loosely coupled and heterogeneous nature of workflows often fails to take full advantage of their capabilities. We created Colmena to leverage the massive parallelism of a supercomputer by using Artificial Intelligence (AI) to learn from and adapt a workflow as it executes. Colmena allows scientists to define how their application should respond to events (e.g., task completion) as a series of cooperative agents. In this paper, we describe the design of Colmena, the challenges we overcame while deploying applications on exascale systems, and the science workflows we have enhanced through interweaving AI. The scaling challenges we discuss include developing steering strategies that maximize node utilization, introducing data fabrics that reduce communication overhead of data-intensive tasks, and implementing workflow tasks that cache costly operations between invocations. These innovations coupled with a variety of application patterns accessible through our agent-based steering model have enabled science advances in chemistry, biophysics, and materials science using different types of AI. In conclusion, our vision is that Colmena will spur creative solutions that harness AI across many domains of scientific computing.

Workflows↗
DOE OSTI2024

Global and regional hydrological impacts of global forest expansion

Abstract. Large-scale reforestation, afforestation, and forest restoration schemes have gained global support as climate change mitigation strategies due to their significant carbon dioxide removal (CDR) potential. However, there has been limited research into the unintended consequences of forestation from a biophysical perspective. In the Community Earth System Model version 2 (CESM2), we apply a global forestation scenario, within a Paris Agreement-compatible warming scenario, to investigate the land surface and hydroclimate response. Compared to a control scenario where land use is fixed to present-day levels, the forestation scenario is up to 2 °C cooler at low latitudes by 2100, driven by a 10 % increase in evaporative cooling in forested areas. However, afforested areas where grassland or shrubland are replaced lead to a doubling of plant water demand in some tropical regions, causing significant decreases in soil moisture (∼ 5 % globally, 5 %–10 % regionally) and water availability (∼ 10 % globally, 10 %–15 % regionally) in regions with increased forest cover. While there are some increases in low cloud and seasonal precipitation over the expanded tropical forests, with enhanced negative cloud radiative forcing, the impacts on large-scale precipitation and atmospheric circulation are limited. This contrasts with the precipitation response to simulated large-scale deforestation found in previous studies. The forestation scenario demonstrates local cooling benefits without major disruption to global hydrodynamics beyond those already projected to result from climate change, in addition to the cooling associated with CDR. However, the water demands of extensive forestation, especially afforestation, have implications for its viability, given the uncertainty in future precipitation changes.

Environmental Sciences & Ecology↗
DOE OSTI2023

Data and scripts associated with the manuscript evaluating the hydrologic responses of the Pacific Northwest watersheds to wildfires (v2)

This data package is associated with the publication “Evaluating Post-fire Watershed Response to Varying Burn Severity and Precipitation Regimes Using Fully-distributed and Integrated Hydrologic Models” submitted to Journal of Hydrology (Li et al. 2025). In this study, we employed the Advanced Terrestrial Simulator (ATS), an integrated watershed model that couples surface flow, subsurface flow, and canopy biophysical processes, to investigate post-fire hydrologic responses in a few selected watersheds with varying burn severity.The data package contains the required input data (meteorological forcing, Leaf Area Index, wildfire burn severities, etc.) to run the model, configuration files, the Jupyter notebooks in Python to pre-process and post-process data, the figures in the manuscript, and the modeling output files. The variables include watershed-averaged evapotranspiration, watershed-averaged surface/subsurface/canopy water content, and river discharge at watershed outlet.The data package contains a file-level metadata that lists and describes all the files contained in the data package (ATS_flmd.csv), a data dictionary file that defines columns headers across all csv files contained in the data package (ATS_dd.csv), a data package level readme file (the current file), and four zipped folders.The ‘data’ folder provides data needed to run the model in .h5, .i2s, .xyz, .shp, and .exo formats. The sub-folders are for each data types. The ‘model’ folder provides input files (.xml format) and essential model outputs. Each sub-folder provides the files from each simulated watershed. The ‘notebooks’ folder provides the Jupyter notebooks (.ipynb format) for pre- and post- processing model files, and for producing the figures in the manuscript. The ‘figures’ folder provides the figures associated with manuscript in .pdf and .png formats.The ‘model’ folder and the ‘data’ folder have been split into 5GB-large pieces using the Linux command ‘split -b 5120m model.zip model.zip.’ and ‘split -b 5120m data.zip data.zip.’, respectively. They can be merged back using the Linux command ‘cat model.zip.* > model.zip’ and ‘cat data.zip.* > data.zip’, respectively.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Complex water networks visualized by cryogenic electron microscopy of RNA

The stability and function of biomolecules are directly influenced by their myriad interactions with water. Here we investigated water through cryogenic electron microscopy (cryo-EM) on a highly solvated molecule: the Tetrahymena ribozyme. By using segmentation-guided water and ion modelling (SWIM), an approach combining resolvability and chemical parameters, we automatically modelled and cross-validated water molecules and Mg 2+ ions in the ribozyme core, revealing the extensive involvement of water in mediating RNA non-canonical interactions. Unexpectedly, in regions where SWIM does not model ordered water, we observed highly similar densities in both cryo-EM maps. In many of these regions, the cryo-EM densities superimpose with complex water networks predicted by molecular dynamics, supporting their assignment as water and suggesting a biophysical explanation for their elusiveness to conventional atomic coordinate modelling. Our study demonstrates an approach to unveil both rigid and flexible waters that surround biomolecules through cryo-EM map densities, statistical and chemical metrics, and molecular dynamics simulations.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Empirical estimation of weather-driven yield shocks using biophysical characteristics for U.S. rainfed and irrigated maize, soybeans, and winter wheat

Agricultural yields are highly susceptible to changes in weather system patterns, including both annual and sub-annual changes in temperature and precipitation. Understanding the impacts of future changes in meteorological variables on crop yields have been widely studied in both empirical and process-based models. This work presents a structural econometric approach that combines historical weather data with information about the biophysical growth cycles of maize, winter wheat, and soy to predict the year-to-year yield responses for rainfed crops. Modeled soil moisture and temperature are taken as key predictors, which allows testing the fitted rainfed model’s ability to predict irrigated yields by assuming irrigation produces the numerically optimal level of soil moisture. This approach is grounded in known biophysical processes, which improves confidence in the model predictions. It produces estimates of interannual variability in yields, which enable study of the effects of extreme events on agricultural productivity. Finally, it enables prediction of the potential impacts of changing weather patterns on irrigated crops in areas that are currently primarily rainfed. We present the results of the empirical model, fitted with rainfed data; out-of-sample validation on irrigated crops; and projections of yield shocks under multiple future climate scenarios. Under a bias-corrected GFDL RCP8.5 scenario, this approach predicts average yield change, relative to 2006-2020 average yields, across U.S. counties in the 2040-2060 period of -16.9% and -13.4% for rainfed and irrigated maize, of -19.1% and -18.4% for rainfed and irrigated soy, and of -4.0% and -2.0% for rainfed and irrigated winter wheat. And in the 2070-2090 period of -30.3% and -28.8 % for rainfed and irrigated maize, of -33.1% and -32.8% for rainfed and irrigated soy, and of -6.5% and -5.2%for rainfed and irrigated winter wheat.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Substrate binding-induced conformational transitions in the omega-3 fatty acid transporter MFSD2A

Major Facilitator Superfamily Domain containing 2 A (MFSD2A) is a transporter that is highly enriched at the blood-brain and blood-retinal barriers, where it mediates Na + -dependent uptake of ω-3 fatty acids in the form of lysolipids into the brain and eyes, respectively. Despite recent structural insights, it remains unclear how this process is initiated, and driven by Na + . Here, we perform Molecular Dynamics simulations which demonstrate that substrates enter outward facing MFSD2A from the outer leaflet of the membrane via lateral openings between transmembrane helices 5/8 and 2/11. The substrate headgroup enters first and engages in Na + -bridged interactions with a conserved glutamic acid, while the tail is surrounded by hydrophobic residues. This binding mode is consistent with a “trap-and-flip” mechanism and triggers transition to an occluded conformation. Furthermore, using machine learning analysis, we identify key elements that enable these transitions. These results advance our molecular understanding of the MFSD2A transport cycle.

59 BASIC BIOLOGICAL SCIENCES↗