Search NASA⌕ Search

SEARCH · Search NASA

Results for “Protein Family”

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 55 records · Page 3

The Legionella Effector SdjA Is a Bifunctional Enzyme That Distinctly Regulates Phosphoribosyl Ubiquitination

Legionella pneumophila promotes its survival and replication in phagocytes by actively modulating cellular processes using effectors injected into host cells by its Dot/Icm type IV secretion system. Many of these effectors function to manipulate the ubiquitin network of infected cells, thus contributing to the biogenesis of the Legionella-containing vacuole (LCV), which is permissive for bacterial replication. Among these, members of the SidE effector family (SidEs) catalyze ubiquitination of functionally diverse host proteins by a mechanism that is chemically distinct from the canonical three-enzyme cascade. The activity of SidEs is regulated by two mechanisms: reversal of the phosphoribosyl ubiquitination by DupA and DupB and direct inactivation by SidJ, which is a calmodulin-dependent glutamylase. In many L. pneumophila strains, SidJ belongs to a two-member protein family. Its homolog SdjA appears to function differently from SidJ despite the high-level similarity in their primary sequences. Here, we found that SdjA is a bifunctional enzyme that exhibits distinct activities toward members of the SidE family. It inhibits the activity of SdeB and SdeC by glutamylation. Unexpectedly, it also functions as a deglutamylase that reverses SidJ-induced glutamylation on SdeA. Our results reveal that an enzyme can catalyze two completely opposite biochemical reactions, which highlights the distinct regulation of phosphoribosyl ubiquitination by the SidJ effector family.

59 BASIC BIOLOGICAL SCIENCES↗

Specialization Restricts the Evolutionary Paths Available to Yeast Sugar Transporters

Functional innovation at the protein level is a key source of evolutionary novelties. The constraints on functional innovations are likely to be highly specific in different proteins, which are shaped by their unique histories and the extent of global epistasis that arises from their structures and biochemistries. These contextual nuances in the sequence–function relationship have implications both for a basic understanding of the evolutionary process and for engineering proteins with desirable properties. Here, we have investigated the molecular basis of novel function in a model member of an ancient, conserved, and biotechnologically relevant protein family. These Major Facilitator Superfamily sugar porters are a functionally diverse group of proteins that are thought to be highly plastic and evolvable. By dissecting a recent evolutionary innovation in an α-glucoside transporter from the yeast Saccharomyces eubayanus, we show that the ability to transport a novel substrate requires high-order interactions between many protein regions and numerous specific residues proximal to the transport channel. To reconcile the functional diversity of this family with the constrained evolution of this model protein, we generated new, state-of-the-art genome annotations for 332 Saccharomycotina yeast species spanning ~400 My of evolution. By integrating phylogenetic and phenotypic analyses across these species, we show that the model yeast α-glucoside transporters likely evolved from a multifunctional ancestor and became subfunctionalized. The accumulation of additive and epistatic substitutions likely entrenched this subfunction, which made the simultaneous acquisition of multiple interacting substitutions the only reasonably accessible path to novelty.

59 BASIC BIOLOGICAL SCIENCES↗

A universal metabolite repair enzyme removes a strong inhibitor of the TCA cycle

Abstract A prevalent side-reaction of succinate dehydrogenase oxidizes malate to enol-oxaloacetate (OAA), a metabolically inactive form of OAA that is a strong inhibitor of succinate dehydrogenase. We purified from cow heart mitochondria an enzyme (OAT1) with OAA tautomerase (OAT) activity that converts enol-OAA to the physiological keto-OAA form, and determined that it belongs to the highly conserved and previously uncharacterized Fumarylacetoacetate_hydrolase_domain-containing protein family. From all three domains of life, heterologously expressed proteins were shown to have strong OAT activity, and ablating the OAT1 homolog caused significant growth defects. In Escherichia coli , expression of succinate dehydrogenase was necessary for OAT1-associated growth defects to occur, and ablating OAT1 caused a significant increase in acetate and other metabolites associated with anaerobic respiration. OAT1 increased the succinate dehydrogenase reaction rate by 35% in in vitro assays with physiological concentrations of both succinate and malate. Our results suggest that OAT1 is a universal metabolite repair enzyme that is required to maximize aerobic respiration efficiency by preventing succinate dehydrogenase inhibition.

59 BASIC BIOLOGICAL SCIENCES↗

Allosterically coupled conformational dynamics in solution prepare the sterol transfer protein StarD4 to release its cargo upon interaction with target membranes

Complex mechanisms regulate the cellular distribution of cholesterol, a critical component of eukaryote membranes involved in regulation of membrane protein functions directly and through the physiochemical properties of membranes. StarD4, a member of the steroidogenic acute regulator-related lipid-transfer (StART) domain (StARD)-containing protein family, is a highly efficient sterol-specific transfer protein involved in cholesterol homeostasis. Its mechanism of cargo loading and release remains unknown despite recent insights into the key role of phosphatidylinositol phosphates in modulating its interactions with target membranes. We have used large-scale atomistic Molecular dynamics (MD) simulations to study how the dynamics of cholesterol bound to the StarD4 protein can affect interaction with target membranes, and cargo delivery. We identify the two major cholesterol (CHL) binding modes in the hydrophobic pocket of StarD4, one near S136&S147 (the Ser-mode), and another closer to the putative release gate located near W171, R92&Y117 (the Trp-mode). We show that conformational changes of StarD4 associated directly with the transition between these binding modes facilitate the opening of the gate. To understand the dynamics of this connection we apply a machine-learning algorithm for the detection of rare events in MD trajectories (RED), which reveals the structural motifs involved in the opening of a front gate and a back corridor in the StarD4 structure occurring together with the spontaneous transition of CHL from the Ser-mode of binding to the Trp-mode. Further analysis of MD trajectory data with the information-theory based NbIT method reveals the allosteric network connecting the CHL binding site to the functionally important structural components of the gate and corridor. Mutations of residues in the allosteric network are shown to affect the performance of the allosteric connection. These findings outline an allosteric mechanism which prepares the CHL-bound StarD4 to release and deliver the cargo when it is bound to the target membrane.

59 BASIC BIOLOGICAL SCIENCES↗

Large language models generate functional protein sequences across diverse families

Deep-learning language models have shown promise in various biotechnological applications, including protein design and engineering. Here, in this paper, we describe ProGen, a language model that can generate protein sequences with a predictable function across large protein families, akin to generating grammatically and semantically correct natural language sentences on diverse topics. The model was trained on 280 million protein sequences from >19,000 families and is augmented with control tags specifying protein properties. ProGen can be further fine-tuned to curated sequences and tags to improve controllable generation performance of proteins from families with sufficient homologous samples. Artificial proteins fine-tuned to five distinct lysozyme families showed similar catalytic efficiencies as natural lysozymes, with sequence identity to natural proteins as low as 31.4%. ProGen is readily adapted to diverse protein families, as we demonstrate with chorismate mutase and malate dehydrogenase.

59 BASIC BIOLOGICAL SCIENCES↗

A widespread family of molecular chaperones promotes the intracellular stability of type VIIb secretion system– exported toxins

To survive in highly competitive environments, bacteria use specialized secretion systems to deliver antibacterial toxins into neighboring cells, thereby inhibiting their growth. In many Gram-positive bacteria, the export of such toxins requires a membrane-bound molecular apparatus known as the type VIIb secretion system (T7SSb). Recently, it was shown that toxin recruitment to the T7SSb requires a physical interaction between a toxin and two or more so-called targeting factors, which harbor key residues required for T7SS-dependent protein export. However, in addition to these targeting factors, some toxins additionally require a protein belonging to the DUF4176 protein family. Here, by examining two toxin–DUF4176 protein pairs, we demonstrate that DUF4176 constitutes a family of toxin-specific molecular chaperones. In addition to being required for toxin stability in producing cells, we find that DUF4176 proteins facilitate toxin export by specifically interacting with a previously uncharacterized intrinsically disordered region found in many T7SS toxins. Using X-ray crystallography, we determine structures of several DUF4176 chaperones in their unbound state, and of a DUF4176 chaperone in complex with the binding site of its cognate toxin. These structures reveal that this binding site consists of a disordered amphipathic α-helix that requires interaction with its cognate chaperone for proper folding. Overall, we have identified a family of secretion system associated molecular chaperones found throughout T7SSb-containing Gram-positive bacteria.

Gkragkopoulou, Polyniki↗

Proteomic and metabolic disturbances in lignin-modified Brachypodium distachyon

Abstract Lignin biosynthesis begins with the deamination of phenylalanine and tyrosine (Tyr) as a key branch point between primary and secondary metabolism in land plants. Here, we used a systems biology approach to investigate the global metabolic responses to lignin pathway perturbations in the model grass Brachypodium distachyon. We identified the lignin biosynthetic protein families and found that ammonia-lyases (ALs) are among the most abundant proteins in lignifying tissues in grasses. Integrated metabolomic and proteomic data support a link between lignin biosynthesis and primary metabolism mediated by the ammonia released from ALs that is recycled for the synthesis of amino acids via glutamine. RNA interference knockdown of lignin genes confirmed that the route of the canonical pathway using shikimate ester intermediates is not essential for lignin formation in Brachypodium, and there is an alternative pathway from Tyr via sinapic acid for the synthesis of syringyl lignin involving yet uncharacterized enzymatic steps. Our findings support a model in which plant ALs play a central role in coordinating the allocation of carbon for lignin synthesis and the nitrogen available for plant growth. Collectively, these data also emphasize the value of integrative multiomic analyses to advance our understanding of plant metabolism.

59 BASIC BIOLOGICAL SCIENCES↗

Structural insights into differences in G protein activation by family A and family B GPCRs

Family B heterotrimeric guanine nucleotide–binding protein (G protein)–coupled receptors (GPCRs) play important roles in carbohydrate metabolism. Recent structures of family B GPCR-G s protein complexes reveal a disruption in the α-helix of transmembrane segment 6 (TM6) not observed in family A GPCRs. To investigate the functional impact of this structural difference, we compared the structure and function of the glucagon receptor (GCGR; family B) with the β 2 adrenergic receptor (β 2 AR; family A). We determined the structure of the GCGR-G s complex by means of cryo–electron microscopy at 3.1-angstrom resolution. This structure shows the distinct break in TM6. Guanosine triphosphate (GTP) turnover, guanosine diphosphate release, GTP binding, and G protein dissociation studies revealed much slower rates for G protein activation by the GCGR compared with the β 2 AR. Fluorescence and double electron-electron resonance studies suggest that this difference is due to the inability of agonist alone to induce a detectable outward movement of the cytoplasmic end of TM6.

59 BASIC BIOLOGICAL SCIENCES↗

Comparative genomics of Aspergillus nidulans and section Nidulantes

Aspergillus nidulans is an important model organism for eukaryotic biology and the reference for the section Nidulantes in comparative studies. In this study, we de novo sequenced the genomes of 25 species of this section. Whole-genome phylogeny of 34 Aspergillus species and Penicillium chrysogenum clarifies the position of clades inside section Nidulantes. Comparative genomics reveals a high genetic diversity between species with 684 up to 2433 unique protein families. Furthermore, we categorized 2118 secondary metabolite gene clusters (SMGC) into 603 families across Aspergilli, with at least 40 % of the families shared between Nidulantes species. Genetic dereplication of SMGC and subsequent synteny analysis provides evidence for horizontal gene transfer of a SMGC. Proteins that have been investigated in A. nidulans as well as its SMGC families are generally present in the section Nidulantes, supporting its role as model organism. The set of genes encoding plant biomass-related CAZymes is highly conserved in section Nidulantes, while there is remarkable diversity of organization of MAT-loci both within and between the different clades. This study provides a deeper understanding of the genomic conservation and diversity of this section and supports the position of A. nidulans as a reference species for cell biology.

Theobald, Sebastian [Technical University of Denma↗

PLAT domain protein 1 (PLAT1/PLAFP) binds to the Arabidopsis thaliana plasma membrane and inserts a lipid

Robust agricultural yields depend on the plant's ability to fix carbon amid variable environmental conditions. Over seasonal and diurnal cycles, the plant must constantly adjust its metabolism according to available resources or external stressors. The metabolic changes that a plant undergoes in response to stress are well understood, but the long-distance signaling mechanisms that facilitate communication throughout the plant are less studied. The phloem is considered the predominant conduit for the bidirectional transport of these signals in the form of metabolites, nucleic acids, proteins, and lipids. Lipid trafficking through the phloem in particular attracted our attention due to its reliance on soluble lipid-binding proteins (LBP) that generate and solubilize otherwise membrane-associated lipids. The Phloem Lipid-Associated Family Protein (PLAFP) from Arabidopsis thaliana is generated in response to abiotic stress as is its lipid-ligand phosphatidic acid (PA). PLAFP is proposed to transport PA through the phloem in response to drought stress. To understand the interactions between PLAFP and PA, nearly 100 independent systems comprised of the protein and one PA, or a plasma membrane containing varying amounts of PA, were simulated using atomistic classical molecular dynamics methods. In these simulations, PLAFP is found to bind to plant plasma membrane models independent of the PA concentration. When bound to the membrane, PLAFP adopts a binding pose where W41 and R82 penetrate the membrane surface and anchor PLAFP. This triggers a separation of the two loop regions containing W41 and R82. Subsequent simulations indicate that PA insert into the β -sandwich of PLAFP, driven by interactions with multiple amino acids besides the W41 and R82 identified during the insertion process. Finally, fine-tuning the protein-membrane and protein-PA interface by mutating a selection of these amino acids may facilitate engineering plant signaling processes by modulating the binding response.

59 BASIC BIOLOGICAL SCIENCES↗

PYK-SubstitutionOME: an integrated database containing allosteric coupling, ligand affinity and mutational, structural, pathological, bioinformatic and computational information about pyruvate kinase isozymes

Interpreting changes in patient genomes, understanding how viruses evolve and engineering novel protein function all depend on accurately predicting the functional outcomes that arise from amino acid substitutions. To that end, the development of first-generation prediction algorithms was guided by historic experimental datasets. However, these datasets were heavily biased toward substitutions at positions that have not changed much throughout evolution (i.e. conserved). Although newer datasets include substitutions at positions that span a range of evolutionary conservation scores, these data are largely derived from assays that agglomerate multiple aspects of function. To facilitate predictions from the foundational chemical properties of proteins, large substitution databases with biochemical characterizations of function are needed. We report here a database derived from mutational, biochemical, bioinformatic, structural, pathological and computational studies of a highly studied protein family—pyruvate kinase (PYK). A centerpiece of this database is the biochemical characterization—including quantitative evaluation of allosteric regulation—of the changes that accompany substitutions at positions that sample the full conservation range observed in the PYK family. We have used these data to facilitate critical advances in the foundational studies of allosteric regulation and protein evolution and as rigorous benchmarks for testing protein predictions. We trust that the collected dataset will be useful for the broader scientific community in the further development of prediction algorithms.

59 BASIC BIOLOGICAL SCIENCES↗

Expression and characterization of spore coat Cot$\mathrm{H}$ kinases from the cellulosomes of anaerobic fungi ( Neocallimastigomycetes )

Anaerobic fungi (Neocallimastigomycetes) found in the guts of herbivores are biomass deconstruction specialists with a remarkable ability to extract sugars from recalcitrant plant material. Anaerobic fungi, as well as many species of anaerobic bacteria, deploy multi-enzyme complexes called cellulosomes, which modularly tether together hydrolytic enzymes, to accelerate biomass hydrolysis. While the majority of genomically encoded cellulosomal genes in Neocallimastigomycetes are biomass degrading enzymes, the second largest family of cellulosomal genes encode spore coat CotH domains, whose contribution to fungal cellulosome and/or cellular function is unknown. Structural bioinformatics of CotH proteins from the anaerobic fungus Piromyces finnis shows anaerobic fungal CotH domains conserve key ATP and Mg 2+ binding motifs from bacterial Bacillus CotH proteins known to act as protein kinases. Experimental characterization further demonstrates ATP hydrolysis activity in the presence and absence of substrate from two cellulosomal P. finnis CotH proteins when recombinantly produced in E. coli. These results present foundational evidence for CotH activity in anaerobic fungi and provide a path towards elucidating the functional contribution of this protein family to fungal cellulosome assembly and activity.

59 BASIC BIOLOGICAL SCIENCES↗

The N‐terminal intrinsically disordered region of Ncb5or docks with the cytochrome b 5 core to form a helical motif that is of ancient origin

Abstract NADH cytochromeb 5 oxidoreductase (Ncb5or) is a cytosolic ferric reductase implicated in diabetes and neurological conditions. Ncb5or comprises cytochromeb 5 (b 5 ) and cytochromeb 5 reductase (b 5 R) domains separated by a CHORD‐Sgt1 (CS) linker domain. Ncb5or redox activity depends on proper inter‐domain interactions to mediate electron transfer from NADH or NADPH via FAD to heme. While full‐length human Ncb5or has proven resistant to crystallization, we have succeeded in obtaining high‐resolution atomic structures of theb 5 domain and a construct containing the CS andb 5 R domains (CS/b 5 R). Ncb5or also contains an N‐terminal intrinsically disordered region of 50 residues that has no homologs in other protein families in animals but features a distinctive, conserved L 34 MDWIRL 40 motif also present in reduced lateral root formation (RLF) protein in rice and increased recombination center 21 in baker's yeast, all attaching to ab 5 domain. After unsuccessful attempts at crystallizing a human Ncb5or construct comprising the N‐terminal region naturally fused to theb 5 domain, we were able to obtain a high‐resolution atomic structure of a recombinant rice RLF construct corresponding to residues 25–129 of human Ncb5or (52% sequence identity; 74% similarity). The structure reveals Trp 120 (corresponding to invariant Trp 37 in Ncb5or) to be part of an 11‐residue α‐helix (S 116 QMDWLKLTRT 126 ) packing against two of the four helices in theb 5 domain that surround heme (α2 and α5). The Trp 120 side chain forms a network of interactions with the side chains of four highly conserved residues corresponding to Tyr 85 and Tyr 88 (α2), Cys 124 (α5), and Leu 47 in Ncb5or. Circular dichroism measurements of human Ncb5or fragments further support a key role of Trp 37 in nucleating the formation of the N‐terminal helix, whose location in the N/b 5 module suggests a role in regulating the function of this multi‐domain redox enzyme. This study revealed for the first time an ancient origin of a helical motif in the N/b 5 module as reflected by its existence in a class of cytochromeb 5 proteins from three kingdoms among eukaryotes.

Biochemistry & Molecular Biology↗

Conformation-specific inhibitors of activated Ras GTPases reveal limited Ras dependency of patient-derived cancer organoids

The small GTPases H, K, and NRAS are molecular switches indispensable for proper regulation of cellular proliferation and growth. Several mutations in the genes encoding members of this protein family are associated with cancer and result in aberrant activation of signaling processes caused by a deregulated recruitment of downstream effector proteins. In this study, we engineered variants of the Ras-binding domain (RBD) of the C-Raf proto-oncogene, Ser/Thr kinase (CRAF). These variants bound with high affinity with the effector-binding site of Ras in an active conformation. Structural characterization disclosed how the newly identified RBD mutations cooperate and thereby enhance affinity with the effector-binding site in Ras compared with WT RBD. The engineered RBD variants closely mimicked the interaction mode of naturally occurring Ras effectors and acted as dominant-negative affinity reagents that block Ras signal transduction. We report experiments with cancer cells showed that expression of these RBD variants inhibits Ras signaling, reducing cell growth and inducing apoptosis. Using these optimized RBD variants, we stratified patient-derived colorectal cancer organoids with known Ras mutational status according to their response to Ras inhibition. These results revealed that the presence of Ras mutations was insufficient to predict sensitivity to Ras inhibition, suggesting that not all of these tumors required Ras signaling for proliferation. In summary, by engineering the Ras/Raf interface of the CRAF-RBD, we identified potent and selective inhibitors of Ras in its active conformation that outcompete binding of Ras-signaling effectors.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Structural- and Functional-Informed Machine Learning for Protein Function Prediction

In this project we aimed to extend methods for protein function prediction to include structural prediction data, and benchmark methods against existing tools. We proposed to apply the method to large metagenome datasets, and develop approaches to examine activity-based protein profiling results for protein function-structure patterns. Nitrogen cycle protein families were previously identified and are used here to provide a proof-of-principle for use of structure prediction in protein function classification.

59 BASIC BIOLOGICAL SCIENCES↗

efam: an e xpanded, metaproteome-supported HMM profile database of viral protein fam ilies

Viruses infect, reprogram and kill microbes, leading to profound ecosystem consequences, from elemental cycling in oceans and soils to microbiome-modulated diseases in plants and animals. Although metagenomic datasets are increasingly available, identifying viruses in them is challenging due to poor representation and annotation of viral sequences in databases. Here, we establish efam, an expanded collection of Hidden Markov Model (HMM) profiles that represent viral protein families conservatively identified from the Global Ocean Virome 2.0 dataset. This resulted in 240 311 HMM profiles, each with at least 2 protein sequences, making efam >7-fold larger than the next largest, pan-ecosystem viral HMM profile database. Adjusting the criteria for viral contig confidence from ‘conservative’ to ‘eXtremely Conservative’ resulted in 37 841 HMM profiles in our efam-XC database. To assess the value of this resource, we integrated efam-XC into VirSorter viral discovery software to discover viruses from less-studied, ecologically distinct oxygen minimum zone (OMZ) marine habitats. This expanded database led to an increase in viruses recovered from every tested OMZ virome by ~24% on average (up to ~42%) and especially improved the recovery of often-missed shorter contigs (<5 kb). Additionally, to help elucidate lesser-known viral protein functions, we annotated the profiles using multiple databases from the DRAM pipeline and virion-associated metaproteomic data, which doubled the number of annotations obtainable by standard, single-database annotation approaches. Together, these marine resources (efam and efam-XC) are provided as searchable, compressed HMM databases that will be updated bi-annually to help maximize viral sequence discovery and study from any ecosystem.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Conserved unique peptide patterns (CUPP) online platform 2.0: implementation of +1000 JGI fungal genomes

Carbohydrate-processing enzymes, CAZymes, are classified into families based on sequence and three-dimensional fold. Because many CAZyme families contain members of diverse molecular function (different EC-numbers), sophisticated tools are required to further delineate these enzymes. Such delineation is provided by the peptide-based clustering method CUPP, Conserved Unique Peptide Patterns. CUPP operates synergistically with the CAZy family/subfamily categorizations to allow systematic exploration of CAZymes by defining small protein groups with shared sequence motifs. The updated CUPP library contains 21,930 of such motif groups including 3,842,628 proteins. The new implementation of the CUPP-webserver, https://cupp.info/, now includes all published fungal and algal genomes from the Joint Genome Institute (JGI), genome resources MycoCosm and PhycoCosm, dynamically subdivided into motif groups of CAZymes. This allows users to browse the JGI portals for specific predicted functions or specific protein families from genome sequences. Thus, a genome can be searched for proteins having specific characteristics. All JGI proteins have a hyperlink to a summary page which links to the predicted gene splicing including which regions have RNA support. The new CUPP implementation also includes an update of the annotation algorithm that uses only a fourth of the RAM while enabling multi-threading, providing an annotation speed below 1 ms/protein.

59 BASIC BIOLOGICAL SCIENCES↗

Cytochrome b 5 : A versatile electron carrier and regulator for plant metabolism

Cytochrome b 5 (CB5) is a small heme-binding protein, known as an electron donor delivering reducing power to the terminal enzymes involved in oxidative reactions. In plants, the CB5 protein family is substantially expanded both in its isoform numbers and cellular functions, compared to its yeast and mammalian counterparts. As an electron carrier, plant CB5 proteins function not only in fatty acid desaturation, hydroxylation and elongation, but also in the formation of specialized metabolites such as flavonoids, phenolic esters, and heteropolymer lignin. Furthermore, plant CB5s are found to interact with different non-catalytic proteins such as ethylene signaling regulator, cell death inhibitor, and sugar transporters, implicating their versatile regulatory roles in coordinating different metabolic and cellular processes, presumably in respect to the cellular redox status and/or carbon availability. Compared to the plentiful studies on biochemistry and cellular functions of mammalian CB5 proteins, the cellular and metabolic roles of plant CB5 proteins have received far less attention. This article summarizes the fragmentary information pertaining to the discovery of plant CB5 proteins, and discusses the conventional and peculiar functions that plant CB5s might play in different metabolic and cellular processes. Gaining comprehensive insight into the biological functions of CB5 proteins could offer effective biotechnological solutions to tailor plant chemodiversity and cellular responses to environment stimuli.

59 BASIC BIOLOGICAL SCIENCES↗