Search NASA⌕ Search

SEARCH · Search NASA

Results for “Sequence motif analysis”

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

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

At least 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Genome-Wide Discovery of G-Quadruplexes in Wheat: Distribution and Putative Functional Roles

G-quadruplexes are nucleic acid secondary structures formed by a stack of square planar G-quartets. G-quadruplexes were implicated in many biological functions including telomere maintenance, replication, transcription, and translation, in many species including humans and plants. For wheat, however, though it is one of the world’s most important staple food, no G-quadruplex studies have been reported to date. Here, we computationally identify putative G4 structures (G4s) in wheat genome for the first time and compare its distribution across the genome against five other genomes (human, maize, Arabidopsis, rice, and sorghum). We identified close to 1 million G4 motifs with a density of 76 G4s/Mb across the whole genome and 93 G4s/Mb over genic regions. Remarkably, G4s were enriched around three regions, two located on the antisense and one on the sense strand at the following positions: 1) the transcription start site (TSS) (antisense), 2) the first coding domain sequence (CDS) (antisense), and 3) the start codon (sense). Functional enrichment analysis revealed that the gene models containing G4 motifs within these peaks were associated with specific gene ontology (GO) terms, such as developmental process, localization, and cellular component organization or biogenesis. We investigated genes encoding MADS-box transcription factors and showed examples of G4 motifs within critical regulatory regions in the VRN-1 genes in wheat. Furthermore, comparison with other plants showed that monocots share a similar distribution of G4s, but Arabidopsis shows a unique G4 distribution. Our study shows for the first time the prevalence and possible functional roles of G4s in wheat.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

A broad specificity β-propeller enzyme from Rhodopseudomonas palustris that hydrolyzes many lactones including γ-valerolactone

Lactones are prevalent in biological and industrial settings, yet there is a lack of information regarding enzymes used to metabolize these compounds. One compound, γ-valerolactone (GVL), is used as a solvent to dissolve plant cell walls into sugars and aromatic molecules for subsequent microbial conversion to fuels and chemicals. Despite the promise of GVL as a renewable solvent for biomass deconstruction, residual GVL can be toxic to microbial fermentation. Here, we identified a Ca 2+ -dependent enzyme from Rhodopseudomonas palustris (Rpa3624) and showed that it can hydrolyze aliphatic and aromatic lactones and esters, including GVL. Maximum-likelihood phylogenetic analysis of other related lactonases with experimentally determined substrate preferences shows that Rpa3624 separates by sequence motifs into a subclade with preference for hydrophobic substrates. Additionally, we solved crystal structures of this β-propeller enzyme separately with either phosphate, an inhibitor, or a mixture of GVL and products to define an active site where calcium-bound water and calcium-bound aspartic and glutamic acid residues make close contact with substrate and product. Our kinetic characterization of WT and mutant enzymes combined with structural insights inform a reaction mechanism that centers around activation of a calcium-bound water molecule promoted by general base catalysis and close contacts with substrate and a potential intermediate. Similarity of Rpa3624 with other β-propeller lactonases suggests this mechanism may be relevant for other members of this emerging class of versatile catalysts.

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

Unveiling a pervasive DNA adenine methylation regulatory network in the early-diverging fungus Rhizopus microsporus

Development of the DNA affinity purification and sequencing (DAP-seq) technique has allowed genome-scale studies of transcription factor (TF)-binding sites with high reproducibility. Here, we apply this technique to the human opportunistic pathogen Rhizopus microsporus, a mucoralean fungus belonging to the understudied group of early-diverging fungi. We characterize genome-wide binding sites of 58 TFs encoded by genes regulated through adenine methylation and representing major TF families. This analysis reveals their binding profiles and recognized sequences, expanding and diversifying the catalog of known fungal motifs. By integrating this data with DNA 6-methyladenine profiling, we uncover the extensive direct and indirect impact of this epigenetic modification on the regulation of gene expression. Furthermore, we use the generated data to identify TFs involved in biologically relevant processes such as zinc metabolism and light response. Our work enhances our understanding of regulatory mechanisms in R. microsporus and provides broader insights into gene regulation across the fungal kingdom.

Lax, Carlos [Universidad de Murcia (Spain)] (ORCID↗
DOE OSTI2024

Library Screening, In Vivo Confirmation, and Structural and Bioinformatic Analysis of Pentapeptide Sequences as Substrates for Protein Farnesyltransferase

Protein farnesylation is a post-translational modification where a 15-carbon farnesyl isoprenoid is appended to the C-terminal end of a protein by farnesyltransferase (FTase). This process often causes proteins to associate with the membrane and participate in signal transduction pathways. The most common substrates of FTase are proteins that have C-terminal tetrapeptide CaaX box sequences where the cysteine is the site of modification. However, recent work has shown that five amino acid sequences can also be recognized, including the pentapeptides CMIIM and CSLMQ. In this work, peptide libraries were initially used to systematically vary the residues in those two parental sequences using an assay based on Matrix Assisted Laser Desorption Ionization–Mass Spectrometry (MALDI-MS). In addition, 192 pentapeptide sequences from the human proteome were screened using that assay to discover additional extended CaaaX-box motifs. Selected hits from that screening effort were rescreened using an in vivo yeast reporter protein assay. The X-ray crystal structure of CMIIM bound to FTase was also solved, showing that the C-terminal tripeptide of that sequence interacted with the enzyme in a similar manner as the C-terminal tripeptide of CVVM, suggesting that the tripeptide comprises a common structural element for substrate recognition in both tetrapeptide and pentapeptide sequences. Molecular dynamics simulation of CMIIM bound to FTase further shed light on the molecular interactions involved, showing that a putative catalytically competent Zn(II)-thiolate species was able to form. Bioinformatic predictions of tetrapeptide (CaaX-box) reactivity correlated well with the reactivity of pentapeptides obtained from in vivo analysis, reinforcing the importance of the C-terminal tripeptide motif. This analysis provides a structural framework for understanding the reactivity of extended CaaaX-box motifs and a method that may be useful for predicting the reactivity of additional FTase substrates bearing CaaaX-box sequences.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Evolutionary Diversification of Host-Targeted Bartonella Effectors Proteins Derived from a Conserved FicTA Toxin-Antitoxin Module

Proteins containing a FIC domain catalyze AMPylation and other post-translational modifications (PTMs). In bacteria, they are typically part of FicTA toxin-antitoxin modules that control conserved biochemical processes such as topoisomerase activity, but they have also repeatedly diversified into host-targeted virulence factors. Among these, Bartonella effector proteins (Beps) comprise a particularly diverse ensemble of FIC domains that subvert various host cellular functions. However, no comprehensive comparative analysis has been performed to infer molecular mechanisms underlying the biochemical and functional diversification of FIC domains in the vast Bep family. Here, we used X-ray crystallography, structural modelling, and phylogenetic analyses to unravel the expansion and diversification of Bep repertoires that evolved in parallel in three Bartonella lineages from a single ancestral FicTA toxin-antitoxin module. Our analysis is based on 99 non-redundant Bep sequences and nine crystal structures. Inferred from the conservation of the FIC signature motif that comprises the catalytic histidine and residues involved in substrate binding, about half of them represent AMP transferases. A quarter of Beps show a glutamate in a strategic position in the putative substrate binding pocket that would interfere with triphosphate-nucleotide binding but may allow binding of an AMPylated target for deAMPylation or another substrate to catalyze a distinct PTM. The β-hairpin flap that registers the modifiable target segment to the active site exhibits remarkable structural variability. The corresponding sequences form few well-defined groups that may recognize distinct target proteins. The binding of Beps to promiscuous FicA antitoxins is well conserved, indicating a role of the antitoxin to inhibit enzymatic activity or to serve as a chaperone for the FIC domain before translocation of the Bep into host cells. Taken together, our analysis indicates a remarkable functional plasticity of Beps that is mostly brought about by structural changes in the substrate pocket and the target dock. These findings may guide future structure–function analyses of the highly versatile FIC domains.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

Resolving local structural motifs across the phase evolution of zinc titanates with computational x-ray absorption spectroscopy

Resolving the local structure motifs that characterize phase evolution as a function of composition is a key challenge in structure characterization of complex materials. Here, in this study, we combine first-principles simulations and x-ray absorption near-edge structures (XANES) analysis to gain insights into the structure evolution revealed by measurements across a combinatorial zinc titanate thin film, which was grown with smoothly varying composition over a wide range of the Ti:Zn ratio. Specifically, we propose a cluster blind-signal-separation (cBSS) method for XANES spectral analysis based on a library of the structures and spectra of representative local motifs. In addition to motifs from zinc titanate crystals, two types of Ti-defect models constructed in this study are key to the understanding of the structure characteristics in the Zn-rich region. The cBSS method makes use of both spectral clustering of the simulated site-XANES spectra library and the BSS procedure to construct high-fidelity spectral basis functions from an experimental spectral sequence. The method provides a rigorous measure of the spectral sensitivity and basis completeness. The results of the XANES analysis are corroborated with other experimental modalities, including x-ray diffraction and spectroscopic ellipsometry, to validate the cBSS method. The calculated motif weights resulting from fitting the XANES spectra with the cBSS basis probe the atomic structure characteristics of both crystalline and amorphous phases as a function of the Ti/Zn composition. The insights of the local structure motif evolution are pivotal to the understanding of the nonmonotonic trend in the optical gap, which may lead to potential applications through tuning the optical properties of zinc titanate. The workflow of the XANES spectral analysis developed in this work can be generalized to construct the structure-property relationship in a broad material space.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Quantifying the Long‐Range Coupling of Electronic Properties in Proteins with ab initio Molecular Dynamics**

Abstract The delicate interplay of covalent and non‐covalent interactions in proteins is inherently quantum mechanical and highly dynamic in nature. To directly interrogate the evolving nature of the electronic structure of proteins, we carry out 100‐ps‐scale ab initio molecular dynamics simulations of three representative small proteins with range‐separated hybrid density functional theory. We quantify the nature and length‐scale of the coupling of residue‐specific charge probability distributions in these proteins. While some nonpolar residues exhibit expectedly narrow charge distributions, most polar and charged residues exhibit broad, multimodal distributions. Even for nonpolar residues, we observe sequence‐specific deviations corresponding to charge accumulation or depletion that would be challenging to capture in a fixed charge force field. We quantify the effect of residue‐residue interactions on charge distributions first with linear cross‐correlations. We then show how additional insight can be gained from evaluating the mutual information of charge distributions. We show that a significant number of residues couple most strongly with residues that are distant in both sequence and space over a range of secondary structures including α‐helical, β‐sheet, disulfide bridging, and lasso motifs. The mutual information analysis is necessary to capture coupling between some polar and charged residues that would be otherwise missed.

Yang, Zhongyue↗
DOE OSTI2023

Plant defensin MtDef4–derived antifungal peptide with multiple modes of action and potential as a bio–inspired fungicide

Chemical fungicides have been instrumental in protecting crops from fungal diseases. However, increasing fungal resistance to many of the single-site chemical fungicides calls for the development of new antifungal agents with novel modes of action (MoA). The sequence-divergent cysteine-rich antifungal defensins with multisite MoA are promising starting templates for design of novel peptide-based fungicides. Here, we experimentally tested such a set of 17-amino-acid peptides containing the γ-core motif of the antifungal plant defensin MtDef4. These designed peptides exhibited antifungal properties different from those of MtDef4. Focused analysis of a lead peptide, GMA4CG_V6, showed that it was a random coil in solution with little or no secondary structure elements. Additionally, it exhibited potent cation-tolerant antifungal activity against the plant fungal pathogen Botrytis cinerea, the causal agent of grey mould disease in fruits and vegetables. Its multisite MoA involved localization predominantly to the plasma membrane, permeabilization of the plasma membrane, rapid internalization into the vacuole and cytoplasm, and affinity for the bioactive phosphoinositides phosphatidylinositol 3-phosphate (PI3P), PI4P, and PI5P. The sequence motif RRRW was identified as a major determinant of the antifungal activity of this peptide. While topical spray application of GMA4CG_V6 on Nicotiana benthamiana and tomato plants provided preventive and curative suppression of grey mould disease symptoms, the peptide was not internalized into plant cells. Our findings open the possibility that truncated and modified defensin-derived peptides containing the γ-core sequence could serve as promising candidates for further development of bio-inspired fungicides.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Hydrated Metal and Metal-Nitrate Complexes in Water: Full Lanthanide(III) Series plus Miscellaneous Metal Ions

This is a dataset of hydrated metal complexes and metal–nitrate hydrated complexes intended for public use, reproducibility, and downstream structural analysis. A key feature is coverage across the full lanthanide(III) series (La–Lu), enabling systematic comparisons of coordination motifs and bonding trends across the entire lanthanide sequence. In addition to the lanthanides, the dataset also includes other metal ions such as UO2(VI), Fe(II), and Fe(III). The dataset provides optimized geometries for hydrated and nitrate-containing hydrated complexes, together with representative ab initio molecular dynamics (AIMD) trajectories saved in standard XYZ formats. The accompanying NWChem input decks enable reproduction of the reported calculations and provide a starting point for extending the simulations to related coordination environments. Computationally, DFT calculations employ the B3LYP functional with DFT-D3BJ dispersion corrections and a COSMO continuum solvent model (dielectric constant 78.4) to represent solvation beyond the explicitly treated first hydration shell. AIMD simulations are performed with the NWChem qmd module at 298 K, integrating nuclear motion with the velocity-Verlet algorithm and controlling temperature using a Nosé–Hoover thermostat. Trajectories are approximately 4.8 ps in length and are used primarily to assess short-time stability of candidate coordination motifs, including (for lanthanides) differences between 8- versus 9-water coordination and comparisons between nitrate-bound and nitrate-free hydrated complexes.

Dinpajooh, Mohammadhasan [Pacific Northwest Nation↗
DOE OSTI2023

Rise and fall of Mott insulating gaps in YNiO 3 paramagnets as a reflection of symmetry breaking and remaking

The YNiO 3 nickelate is a paradigm d-electron oxide that manifests the intriguing temperature-mediated sequence of three phases transitions from (i) magnetically ordered insulator to (ii) paramagnetic (PM) insulator and then to (iii) PM metal. Such phenomena raised the question of the nature of the association of magnetism and structural symmetry breaking with the appearance in (i) and (ii) and disappearance in (iii) of insulating band gaps. It is demonstrated here that first-principles mean-field–like density-functional theory (DFT), driven by molecular dynamics temperature evolution, can describe not only the origin of the magnetically long-range ordered insulating phase (i), but also the creation of an insulating paramagnet (ii) that lacks spin- long-range order, and of a metallic paramagnet (iii) as temperature rises. This approach provides the patterns of structural and magnetic symmetry breaking at different temperatures, in parallel with band gaps obtained when the evolving geometries are used as input to DFT electronic band-structure calculations. This disentangles the complex interplay among spin, charge, and orbital degrees of freedom. Analysis shows that the success in describing the rise and fall of the insulating band gaps along the phase transition sequence is enabled by allowing sufficient flexibility in describing diverse local structural and magnetic motifs as input to DFT. This entails the use of sufficiently large supercells that allow expressing structural disproportionation of octahedra, as well as a description of PM phases as a distribution of local magnetic moments (rather than using a single averaged moment). It appears that the historic dismissal of mean-field–like DFT as being unable to describe such Mott-like transitions was premature, as it was based on consideration of averaged crystallographic unit cells, a description that washes out local symmetry-breaking motifs. The magnetically ordered insulating YNiO 3 phase (i) and the PM insulating phase (ii) result in DFT from allowing symmetry breaking, evident already by considering the athermal internal energy. In contrast, the PM metallic phase (iii) is formed thermally by smearing out thus weakening symmetry breaking. Analysis of snapshots of the different forms of structural vs magnetic symmetry breaking shows that only the loss of the polymorphous distribution of magnetic moments existing in (ii) causes the fall of the band gap, resulting in the metallic state in (iii). Finally, the interesting conclusion is that such a description of the rise [in phases (i) and (ii)] and fall [in phase (iii)] of the insulating gap does not rely on the traditional Mott-like strong correlation understanding, but on breaking and remaking of magnetic and structural symmetries reflected in energy lowering.

36 MATERIALS SCIENCE↗
DOE OSTI2025

On the distribution of stacking faults at dissociated medium-angle grain boundaries: Crystallographic geometry and metastability

Grain boundaries in FCC metals with low stacking-fault energy can form in dissociated configurations of stacking faults. Perhaps the most studied example is 9R stacking at boundaries near Σ3{112}, where the distribution of stacking faults is related to the emission of Shockley partial dislocations. Here, we combine atomic-resolution electron microscopy, atomistic simulations, and dislocation theory to demonstrate that boundaries vicinal to Σ33a support the stabilization of dissociated 9R stacking within a narrow range of inclinations. This boundary is interesting since its misorientation (20.05°) lies in the medium-angle regime, just past the upper misorientation limit for low-angle boundaries, motivating questions for how best to describe it in terms of dislocations. Our HAADF-STEM observations of thin film bicrystals, supported by atomistic modeling, reveal that this inclination dependence arises from specific geometric constraints on the arrangement of Shockley partial dislocations at the interface. Quantification of stacking-fault distributions across multiple boundaries indicates that the density and spacing of faults closely follow the ideal 9R motif, with subtle variations reflecting the complex energy landscape of these boundaries. Through energetic analysis, we establish the presence of competing metastable states enabled by variations in stacking sequences, emphasizing the significant role of crystallographic geometry. We generalize our analysis as a function of misorientation, showing how 9R at the Σ33a boundary is related to previous observations and calculations of HCP at a 29.7° boundary. This study provides a crystallographically grounded framework connecting dislocation structures, stacking-fault distributions, and metastability at grain boundaries in FCC metals.

Atomistic modeling↗
DOE OSTI2020

Analysis of a preQ1-I riboswitch in effector-free and bound states reveals a metabolite-programmed nucleobase-stacking spine that controls gene regulation

Abstract Riboswitches are structured RNA motifs that recognize metabolites to alter the conformations of downstream sequences, leading to gene regulation. To investigate this molecular framework, we determined crystal structures of a preQ1-I riboswitch in effector-free and bound states at 2.00 Å and 2.65 Å-resolution. Both pseudoknots exhibited the elusive L2 loop, which displayed distinct conformations. Conversely, the Shine-Dalgarno sequence (SDS) in the S2 helix of each structure remained unbroken. The expectation that the effector-free state should expose the SDS prompted us to conduct solution experiments to delineate environmental changes to specific nucleobases in response to preQ1. We then used nudged elastic band computational methods to derive conformational-change pathways linking the crystallographically-determined effector-free and bound-state structures. Pathways featured: (i) unstacking and unpairing of L2 and S2 nucleobases without preQ1—exposing the SDS for translation and (ii) stacking and pairing L2 and S2 nucleobases with preQ1—sequestering the SDS. Our results reveal how preQ1 binding reorganizes L2 into a nucleobase-stacking spine that sequesters the SDS, linking effector recognition to biological function. The generality of stacking spines as conduits for effector-dependent, interdomain communication is discussed in light of their existence in adenine riboswitches, as well as the turnip yellow mosaic virus ribosome sensor.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Two Alternative Splicing Variants of AtERF73/HRE1, HRE1α and HRE1β, Have Differential Transactivation Activities in Arabidopsis

AtERF73/HRE1 is an AP2/ERF transcription factor in Arabidopsis and has two distinct alternative splicing variants, HRE1α and HRE1β. In this study, we examined the differences between the molecular functions of HRE1α and HRE1β. We found that HRE1α and HRE1β are both involved in hypoxia response and root development and have transactivation activity. Two conserved motifs in the C-terminal region of HRE1α and HRE1β, EELL and LWSY-like, contributed to their transactivation activity, specifically the four E residues in the EELL motif and the MGLWS amino acid sequence at the end of the LWSY-like motif. The N-terminal region of HRE1β also showed transactivation activity, mediated by the VDDG motif, whereas that of HRE1α did not. The transactivation activity of HRE1β was stronger than that of HRE1α in Arabidopsis protoplasts. Both transcription factors transactivated downstream genes via the GCC box. RNA-sequencing analysis further supported that both HRE1α and HRE1β might regulate gene expression associated with the hypoxia stress response, although they may transactivate different subsets of genes in downstream pathways. Our results, together with previous studies, suggested that HRE1α and HRE1β differentially transactivate downstream genes in hypoxia response and root development in Arabidopsis.

36 MATERIALS SCIENCE↗
DOE OSTI2024

A compact stem-loop DNA aptamer targets a uracil-binding pocket in the SARS-CoV-2 nucleocapsid RNA-binding domain

Abstract SARS-CoV-2 nucleocapsid (N) protein is a structural component of the virus with essential roles in the replication and packaging of the viral RNA genome. The N protein is also an important target of COVID-19 antigen tests and a promising vaccine candidate along with the spike protein. Here, we report a compact stem-loop DNA aptamer that binds tightly to the N-terminal RNA-binding domain of SARS-CoV-2 N protein. Crystallographic analysis shows that a hexanucleotide DNA motif (5′-TCGGAT-3′) of the aptamer fits into a positively charged concave surface of N-NTD and engages essential RNA-binding residues including Tyr109, which mediates a sequence-specific interaction in a uracil-binding pocket. Avid binding of the DNA aptamer allows isolation and sensitive detection of full-length N protein from crude cell lysates, demonstrating its selectivity and utility in biochemical applications. We further designed a chemically modified DNA aptamer and used it as a probe to examine the interaction of N-NTD with various RNA motifs, which revealed a strong preference for uridine-rich sequences. Our studies provide a high-affinity chemical probe for the SARS-CoV-2 N protein RNA-binding domain, which may be useful for diagnostic applications and investigating novel antiviral agents.

Biochemistry & Molecular Biology↗
DOE OSTI2022

Marine DNA methylation patterns are associated with microbial community composition and inform virus-host dynamics

Background: DNA methylation in prokaryotes is involved in many different cellular processes including cell cycle regulation and defense against viruses. To date, most prokaryotic methylation systems have been studied in culturable microorganisms, resulting in a limited understanding of DNA methylation from a microbial ecology perspective. Here, we analyze the distribution patterns of several microbial epigenetics marks in the ocean microbiome through genome-centric metagenomics across all domains of life. Results: We reconstructed 15,056 viral, 252 prokaryotic, 56 giant viral, and 6 eukaryotic metagenome-assembled genomes from northwest Pacific Ocean seawater samples using short- and long-read sequencing approaches. These metagenome-derived genomes mostly represented novel taxa, and recruited a majority of reads. Thanks to single-molecule real-time (SMRT) sequencing technology, base modification could also be detected for these genomes. This showed that DNA methylation can readily be detected across dominant oceanic bacterial, archaeal, and viral populations, and microbial epigenetic changes correlate with population differentiation. Furthermore, our genome-wide epigenetic analysis of Pelagibacter suggests that GANTC, a DNA methyltransferase target motif, is related to the cell cycle and is affected by environmental conditions. Yet, the presence of this motif also partitions the phylogeny of the Pelagibacter phages, possibly hinting at a competitive co-evolutionary history and multiple effects of a single methylation mark. Conclusions: Overall, this study elucidates that DNA methylation patterns are associated with ecological changes and virus-host dynamics in the ocean microbiome.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Genome-Wide Analysis of Terpene Synthase Gene Family in Mentha longifolia and Catalytic Activity Analysis of a Single Terpene Synthase

Terpenoids are a wide variety of natural products and terpene synthase (TPS) plays a key role in the biosynthesis of terpenoids. Mentha plants are rich in essential oils, whose main components are terpenoids, and their biosynthetic pathways have been basically elucidated. However, there is a lack of systematic identification and study of TPS in Mentha plants. In this work, we genome-widely identified and analyzed the TPS gene family in Mentha longifolia, a model plant for functional genomic research in the genus Mentha. A total of 63 TPS genes were identified in the M. longifolia genome sequence assembly, which could be divided into six subfamilies. The TPS-b subfamily had the largest number of genes, which might be related to the abundant monoterpenoids in Mentha plants. The TPS-e subfamily had 18 members and showed a significant species-specific expansion compared with other sequenced Lamiaceae plant species. The 63 TPS genes could be mapped to nine scaffolds of the M. longifolia genome sequence assembly and the distribution of these genes is uneven. Tandem duplicates and fragment duplicates contributed greatly to the increase in the number of TPS genes in M. longifolia. The conserved motifs (RR(X)8W, NSE/DTE, RXR, and DDXXD) were analyzed in M. longifolia TPSs, and significant differentiation was found between different subfamilies. Adaptive evolution analysis showed that M. longifolia TPSs were subjected to purifying selection after the species-specific expansion, and some amino acid residues under positive selection were identified. Furthermore, we also cloned and analyzed the catalytic activity of a single terpene synthase, MlongTPS29, which belongs to the TPS-b subfamily. MlongTPS29 could encode a limonene synthase and catalyze the biosynthesis of limonene, an important precursor of essential oils from the genus Mentha. This study provides useful information for the biosynthesis of terpenoids in the genus Mentha.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Cation and Anion Channelrhodopsins: Sequence Motifs and Taxonomic Distribution

Cation and anion channelrhodopsins (CCRs and ACRs, respectively) primarily from two algal species, Chlamydomonas reinhardtii and Guillardia theta, have become widely used as optogenetic tools to control cell membrane potential with light. We mined algal and other protist polynucleotide sequencing projects and metagenomic samples to identify 75 channelrhodopsin homologs from four channelrhodopsin families, including one revealed in dinoflagellates in this study. We carried out electrophysiological analysis of 33 natural channelrhodopsin variants from different phylogenetic lineages and 10 metagenomic homologs in search of sequence determinants of ion selectivity, photocurrent desensitization, and spectral tuning in channelrhodopsins. Our results show that association of a reduced number of glutamates near the conductance path with anion selectivity depends on a wider protein context, because prasinophyte homologs with a glutamate pattern identical to that in cryptophyte ACRs are cation selective. Desensitization is also broadly context dependent, as in one branch of stramenopile ACRs and their metagenomic homologs, its extent roughly correlates with phylogenetic relationship of their sequences. Regarding spectral tuning, we identified two prasinophyte CCRs with red-shifted spectra to 585 nm. They exhibit a third residue pattern in their retinal-binding pockets distinctly different from those of the only two types of red-shifted channelrhodopsins known (i.e., the CCR Chrimson and RubyACRs). In cryptophyte ACRs we identified three specific residue positions in the retinal-binding pocket that define the wavelength of their spectral maxima. Lastly, we found that dinoflagellate rhodopsins with a TCP motif in the third transmembrane helix and a metagenomic homolog exhibit channel activity. Channelrhodopsins are widely used in neuroscience and cardiology as research tools and are considered prospective therapeutics, but their natural diversity and mechanisms remain poorly characterized. Genomic and metagenomic sequencing projects are producing an ever-increasing wealth of data, whereas biophysical characterization of the encoded proteins lags behind. In this study, we used manual and automated patch clamp recording of representative members of four channelrhodopsin families, including a family in dinoflagellates that we report in this study. Our results contribute to a better understanding of molecular determinants of ionic selectivity, photocurrent desensitization, and spectral tuning in channelrhodopsins.

59 BASIC BIOLOGICAL SCIENCES↗