Search NASA⌕ Search

SEARCH · Search NASA

Results for “Amino acid sequence 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 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Molecular evolution of human coronavirus-NL63, -229E, -HKU1 and -OC43 in hospitalized children in China

Human coronaviruses (HCoVs) HCoV-NL63, HCoV-229E, HCoV-HKU1 and HCoV-OC43 have been circulated in the human population worldwide, and they are associated with a broad range of respiratory diseases with varying severity. However, there are neither effective therapeutic drugs nor licensed vaccines available for the treatment and prevention of infections by the four HCoVs. In this study, we collected nasopharyngeal aspirates of children hospitalized for respiratory tract infection in China during 2014–2018 and conducted next-generation sequencing. Sequences of four HCoVs were then selected for an in-depth analysis. Genome sequences of 2 HCoV-NL63, 8 HCoV-229E, 2 HCoV-HKU1, and 6 HCoV-OC43 were obtained. Based on the full-length S gene, a strong temporal signal was found in HCoV-229E and the molecular evolutionary rate was 6 × 10 −4 substitutions/site/year. Based on the maximum-likelihood (ML) phylogenetic tree of complete S gene, we designated H78 as a new sub-genotype C2 of HCoV-HKU1, and the obtained P43 sequence was grouped into the reported novel genotype K of HCoV-OC43 circulating in Guangzhou, China. Based on the complete genome, potential recombination events were found to occur as two phenomena, namely intraspecies and interspecies. Moreover, we observed two amino acid substitutions in the S1 subunit of obtained HCoV-NL63 (G534V) and HCoV-HKU1 (H512R), while residues 534 and 512 are important for the binding of angiotensin-converting enzyme 2 and neutralizing antibodies, respectively. Our findings might provide a clue for the molecular evolution of the four HCoVs and help in the early diagnosis, treatment and prevention of broad-spectrum HCoV infection.

Shao, Nan↗
DOE OSTI2021

Toho-1 β-lactamase: backbone chemical shift assignments and changes in dynamics upon binding with avibactam

Backbone chemical shift assignments for the Toho-1 β-lactamase (263 amino acids, 28.9 kDa) are reported based on triple resonance solution-state NMR experiments performed on a uniformly 2 H, 13 C, 15 N-labeled sample. These assignments allow for subsequent site-specific characterization at the chemical, structural, and dynamical levels. At the chemical level, titration with the non-β-lactam β-lactamase inhibitor avibactam is found to give chemical shift perturbations indicative of tight covalent binding that allow for mapping of the inhibitor binding site. At the structural level, protein secondary structure is predicted based on the backbone chemical shifts and protein residue sequence using TALOS-N and found to agree well with structural characterization from X-ray crystallography. At the dynamical level, model-free analysis of 15 N relaxation data at a single field of 16.4 T reveals well-ordered structures for the ligand-free and avibactam-bound enzymes with generalized order parameters of ~0.85. Complementary relaxation dispersion experiments indicate that there is an escalation in motions on the millisecond timescale in the vicinity of the active site upon substrate binding. The combination of high rigidity on short timescales and active site flexibility on longer timescales is consistent with hypotheses for achieving both high catalytic efficiency and broad substrate specificity: the induced active site dynamics allows variously sized substrates to be accommodated and increases the probability that the optimal conformation for catalysis will be sampled.

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

classLog: Logistic regression for the classification of genetic sequences

Introduction Sequencing and phylogenetic classification have become a common task in human and animal diagnostic laboratories. It is routine to sequence pathogens to identify genetic variations of diagnostic significance and to use these data in realtime genomic contact tracing and surveillance. Under this paradigm, unprecedented volumes of data are generated that require rapid analysis to provide meaningful inference. Methods We present a machine learning logistic regression pipeline that can assign classifications to genetic sequence data. The pipeline implements an intuitive and customizable approach to developing a trained prediction model that runs in linear time complexity, generating accurate output rapidly, even with incomplete data. Our approach was benchmarked against porcine respiratory and reproductive syndrome virus (PRRSv) and swine H1 influenza A virus (IAV) datasets. Trained classifiers were tested against sequences and simulated datasets that artificially degraded sequence quality at 0, 10, 20, 30, and 40%. Results When applied to a poor-quality sequence data, the classifier achieved between >85% to 95% accuracy for the PRRSv and the swine H1 IAV HA dataset and this increased to near perfect accuracy when using the full dataset. The model also identifies amino acid positions used to determine genetic clade identity through a feature selection ranking within the model. These positions can be mapped onto a maximum-likelihood phylogenetic tree, allowing for the inference of clade defining mutations. Discussion Our approach is implemented as a python package with code available at https://github.com/flu-crew/classLog .

Zeller, Michael A.↗
DOE OSTI2022

Natural variation and improved genome annotation of the emerging biofuel crop field pennycress ( Thlaspi arvense )

The Brassicaceae family comprises more than 3,700 species with a diversity of phenotypic characteristics, including seed oil content and composition. Recently, the global interest in Thlaspi arvense L. (pennycress) has grown as the seed oil composition makes it a suitable source for biodiesel and aviation fuel production. However, many wild traits of this species need to be domesticated to make pennycress ideal for cultivation. Molecular breeding and engineering efforts require the availability of an accurate genome sequence of the species. Here, we describe pennycress genome annotation improvements, using a combination of long- and short-read transcriptome data obtained from RNA derived from embryos of 22 accessions, in addition to public genome and gene expression information. Our analysis identified 27,213 protein-coding genes, as well as on average 6,188 biallelic SNPs. In addition, we used the identified SNPs to evaluate the population structure of our accessions. The data from this analysis support that the accession Ames 32872, originally from Armenia, is highly divergent from the other accessions, while the accessions originating from Canada and the United States cluster together. When we evaluated the likely signatures of natural selection from alternative SNPs, we found 7 candidate genes under likely recent positive selection. These genes are enriched with functions related to amino acid metabolism and lipid biosynthesis and highlight possible future targets for crop improvement efforts in pennycress.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Graphene oxide exposure alters gut microbial community composition and metabolism in an in vitro human model

Graphene oxide (GO) nanomaterials have unique physicochemical properties that make them highly promising for biomedical, environmental, and agricultural applications. Despite the increasing interest and the use of GO, assessments of its nanotoxicity have largely not interrogated its potential impact on the gut microbiome. This study addresses an important knowledge gap by investigating the impact of GO exposure- both at low (25 ppm) and high (250ppm) doses and fed (nutrient rich) and fasted (nutrient deplete) conditions- on the gut microbial community structure and function, using an in vitro human colon bioreactor model. 16S rRNA amplicon sequencing revealed that GO exposure resulted in a restructuring of community composition. 25 ppm GO induced a marked decrease in the Bacteroidota phylum and increased the ratio of Firmicutes to Bacteroidota (F/B). Untargeted metabolomics on the supernatants indicated that 25 ppm GO impaired microbial utilization and metabolism of substrates (amino acids, carbohydrate metabolites) and reduced production of beneficial microbial metabolites such as 5-hydroxyindole-3-acetic acid and GABA. Exposure to 250 ppm GO resulted in community composition and metabolome profiles that were very similar to the controls that lacked both GO and digestive enzymes, suggesting that high concentrations of GO may interact with digestive enzymes to form protein coronas, causing their depletion in the gut environment. Differential abundance analyses revealed that 3 genera from the phylum Bacteroidota (Bacteroides, Dysgonomonas, and Parabacteroides) were more abundant after 250 ppm GO exposure, irrespective of feed state. Integrative correlation network analysis indicated that the phylum Bacteroidota showed strong positive correlations to multiple microbial metabolites including GABA and 3-indoleacetic acid, are much larger number of correlations compared to other phyla. These results show that GO exposure has a significant impact on gut microbial community composition and metabolism and different mechanisms are at play for low and high GO concentrations.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Ameloblastin binding to biomimetic models of cell membranes – A continuum of intrinsic disorder

A 37-residue amino acid sequence corresponding to the segment encoded by exon-5 of murine ameloblastin (Ambn), AB2 (Y67-Q103), has been implicated with membrane association, ameloblastin self-assembly, and amelogenin-binding. Here, our aim was to characterize, at the residue level, the structural behavior of AB2 bound to chemical mimics of biological membranes using NMR spectroscopy. To better define the structure of AB2 using NMR-based methods, recombinant 13 C- and 15 N-labelled AB2 (*AB2) was prepared and data collected free in solution and with deuterated dodecylphosphocholine (dPC) micelles, deuterated bicelles, and both small and large unilamellar vesicles. Amide chemical shift and intensity perturbations observed in 1 H- 15 N HSQC spectra of *AB2 in the presence of bicelles and dPC micelles suggest that a region of *AB2, S6-E36 (murine Ambn S68 – E98), associates with the membrane biomimetics. A CSI-3 analysis of the NMR chemical shift assignments for *AB2 free in solution and bound to dPC micelles indicated the peptide remains disordered except for the adoption of a short, 12-residue α-helix, F10-G21 (murine Ambn F72-G83). In dPC micelles, the NOE NMR data was void of patterns characteristic of long-lived helical structure indicating this helix was transient in nature. A continuum of intrinsic disorder in the membrane-bound state may be responsible for ameloblastin’s ability to dynamically interact with multiple partners at the same site during amelogenesis.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Reekeekee- and roodoodooviruses, two different Microviridae clades constituted by the smallest DNA phages

Small circular single-stranded DNA viruses of the Microviridae family are both prevalent and diverse in all ecosystems. They usually harbor a genome between 4.3 and 6.3 kb, with a microvirus recently isolated from a marine Alphaproteobacteria being the smallest known genome of a DNA phage (4.248 kb). A subfamily, Amoyvirinae, has been proposed to classify this virus and other related small Alphaproteobacteria-infecting phages. Here, we report the discovery, in meta-omics data sets from various aquatic ecosystems, of sixteen complete microvirus genomes significantly smaller (2.991–3.692 kb) than known ones. Phylogenetic analysis reveals that these sixteen genomes represent two related, yet distinct and diverse, novel groups of microviruses—amoyviruses being their closest known relatives. We propose that these small microviruses are members of two tentatively named subfamilies Reekeekeevirinae and Roodoodoovirinae. As known microvirus genomes encode many overlapping and overprinted genes that are not identified by gene prediction software, we developed a new methodology to identify all genes based on protein conservation, amino acid composition, and selection pressure estimations. Surprisingly, only four to five genes could be identified per genome, with the number of overprinted genes lower than that in phiX174. These small genomes thus tend to have both a lower number of genes and a shorter length for each gene, leaving no place for variable gene regions that could harbor overprinted genes. Even more surprisingly, these two Microviridae groups had specific and different gene content, and major differences in their conserved protein sequences, highlighting that these two related groups of small genome microviruses use very different strategies to fulfill their lifecycle with such a small number of genes. The discovery of these genomes and the detailed prediction and annotation of their genome content expand our understanding of ssDNA phages in nature and are further evidence that these viruses have explored a wide range of possibilities during their long evolution.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Population genomics of the pathogenic yeast Candida tropicalis identifies hybrid isolates in environmental samples

Candida tropicalis is a human pathogen that primarily infects the immunocompromised. Whereas the genome of one isolate, C. tropicalis MYA-3404, was originally sequenced in 2009, there have been no large-scale, multi-isolate studies of the genetic and phenotypic diversity of this species. Here, we used whole genome sequencing and phenotyping to characterize 77 isolates of C. tropicalis from clinical and environmental sources from a variety of locations. We show that most C. tropicalis isolates are diploids with approximately 2–6 heterozygous variants per kilobase. The genomes are relatively stable, with few aneuploidies. However, we identified one highly homozygous isolate and six isolates of C. tropicalis with much higher heterozygosity levels ranging from 36–49 heterozygous variants per kilobase. Our analyses show that the heterozygous isolates represent two different hybrid lineages, where the hybrids share one parent (A) with most other C. tropicalis isolates, but the second parent (B or C) differs by at least 4% at the genome level. Four of the sequenced isolates descend from an AB hybridization, and two from an AC hybridization. The hybrids are MTLa/α heterozygotes. Hybridization, or mating, between different parents is therefore common in the evolutionary history of C. tropicalis. The new hybrids were predominantly found in environmental niches, including from soil. Hybridization is therefore unlikely to be associated with virulence. In addition, we used genotype-phenotype correlation and CRISPR-Cas9 editing to identify a genome variant that results in the inability of one isolate to utilize certain branched-chain amino acids as a sole nitrogen source.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

NEAR: Neural Embeddings for Amino acid Relationships

Protein language models (PLMs) have recently demonstrated potential to supplant classical protein database search methods based on sequence alignment, but are slower than common alignment-based tools and appear to be prone to a high rate of false labeling. Here, we present NEAR, a method based on neural representation learning that is designed to improve both speed and accuracy of search for likely homologs in a large protein sequence database. NEAR’s ResNet embedding model is trained using contrastive learning guided by trusted sequence alignments. It computes per-residue embeddings for target and query protein sequences, and identifies alignment candidates with a pipeline consisting of residue-level k-NN search and a simple neighbor aggregation scheme. Tests on a benchmark consisting of trusted remote homologs and randomly shuffled decoy sequences reveal that NEAR substantially improves accuracy relative to state-of-the-art PLMs, with lower memory requirements and faster embedding and search speed. While these results suggest that the NEAR model may be useful for standalone homology detection with increased sensitivity over standard alignment-based methods, in this manuscript we focus on a more straightforward analysis of the model’s value as a high-speed pre-filter for sensitive annotation. In that context, NEAR is at least 5x faster than the pre-filter currently used in the widely-used profile hidden Markov model (pHMM) search tool HMMER3, and also outperforms the pre-filter used in our fast pHMM tool, nail.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Metabolic activity and community structure of prokaryotes associated with particles in the twilight zone of the South China Sea

The twilight zone is an important depth of the ocean where particulate organic matter (POM) remineralization takes place, and prokaryotes contribute to more than 70% of the estimated remineralization. However, little is known about the microbial community and metabolic activity associated with different particles in the twilight zone. The composition and distribution of particle-attached prokaryotes in the twilight zone of the South China Sea (SCS) were investigated using high-throughput sequencing and quantitative PCR, together with the Biolog Ecoplate™ microplates culture to analyze the microbial metabolic activity. We found that α- and γ-Proteobacteria dominating at the lower and upper boundary of the twilight zone, respectively; Methanosarcinales and Halobacteriales of the Euyarchaeota occupied in the larger particles at the upper boundary. Similar microbial community existed between euphotic layer and the upper boundary. Higher amount of shared Operational Taxonomic Units (OTUs) in the larger particles along the water depths, might be due to the fast sinking and major contribution of carbon flux of the larger particles from the euphotic layer. In addition to polymers as the major carbon source, carbohydrates and amino acids were preferentially used by microbial community at the upper and lower boundary, respectively. This could potentially be attributed to the metabolic capabilities of attached microbial groups in different particles, and reflected the initial preference of the carbon source by the natural microbes in the twilight zone as well. The microbial structure and carbon metabolic profiles could be complemented with metatranscriptomic analysis in future studies to augment the understanding of the complex carbon cycling pathways in the twilight zone.

Liu, Hao↗
DOE OSTI2021

Towards understanding the formation of internal fragments generated by collisionally activated dissociation for top-down mass spectrometry

Top-down mass spectrometry (TD-MS) generates fragment ions that returns information on the polypeptide amino acid sequence. In addition to terminal fragments, internal fragments that result from multiple cleavage events can also be formed. Traditionally, internal fragments are largely ignored due to a lack of available software to reliably assign them, mainly caused by a poor understanding of their formation mechanism. To accurately assign internal fragments, their formation process needs to be better understood. Here we applied a statistical method to compare fragmentation patterns of internal and terminal fragments of peptides and proteins generated by collisionally activated dissociation (CAD). Internal fragments share similar fragmentation propensities with terminal fragments (e.g., enhanced cleavages N-terminal to proline and C-terminal to acidic residues), suggesting that their formation follows conventional CAD pathways. Internal fragments should be generated by subsequent cleavages of terminal fragments and their formation can be explained by the well-known mobile proton model. Additionally, internal fragments can be coupled with terminal fragments to form complementary product ions that span the entire protein sequence. These enhance our understanding of internal fragment formation and can help improve sequencing algorithms to accurately assign internal fragments, which will ultimately lead to more efficient and comprehensive TD-MS analysis of proteins and proteoforms.

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

Neotropical termite microbiomes as sources of novel plant cell wall degrading enzymes

In this study, we used shotgun metagenomic sequencing to characterise the microbial metabolic potential for lignocellulose transformation in the gut of two colonies of Argentine higher termite species with different feeding habits, Cortaritermes fulviceps and Nasutitermes aquilinus. Our goal was to assess the microbial community compositions and metabolic capacity, and to identify genes involved in lignocellulose degradation. Individuals from both termite species contained the same five dominant bacterial phyla (Spirochaetes, Firmicutes, Proteobacteria, Fibrobacteres and Bacteroidetes) although with different relative abundances. However, detected functional capacity varied, with C. fulviceps (a grass-wood-feeder) gut microbiome samples containing more genes related to amino acid metabolism, whereas N. aquilinus (a wood-feeder) gut microbiome samples were enriched in genes involved in carbohydrate metabolism and cellulose degradation. The C. fulviceps gut microbiome was enriched specifically in genes coding for debranching- and oligosaccharide-degrading enzymes. These findings suggest an association between the primary food source and the predicted categories of the enzymes present in the gut microbiomes of each species. To further investigate the termite microbiomes as sources of biotechnologically relevant glycosyl hydrolases, a putative GH10 endo-β-1,4-xylanase, Xyl10E, was cloned and expressed in Escherichia coli. Functional analysis of the recombinant metagenome-derived enzyme showed high specificity towards beechwood xylan (288.1 IU/mg), with the optimum activity at 50 °C and a pH-activity range from 5 to 10. These characteristics suggest that Xy110E may be a promising candidate for further development in lignocellulose deconstruction applications.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Defining upstream enhancing and inhibiting sequence patterns for plant peroxisome targeting signal type 1 using large–scale in silico and in vivo analyses

Peroxisomes are universal eukaryotic organelles essential to plants and animals. Most peroxisomal matrix proteins carry peroxisome targeting signal type 1 (PTS1), a C-terminal tripeptide. Studies from various kingdoms have revealed influences from sequence upstream of the tripeptide on peroxisome targeting, supporting the view that positive charges in the upstream region are the major enhancing elements. However, a systematic approach to better define the upstream elements influencing PTS1 targeting capability is needed. Here, we used protein sequences from 177 plant genomes to perform large-scale and in-depth analysis of the PTS1 domain, which includes the PTS1 tripeptide and upstream sequence elements. We identified and verified 12 low-frequency PTS1 tripeptides and revealed upstream enhancing and inhibiting sequence patterns for peroxisome targeting, which were subsequently validated in vivo. Follow-up analysis revealed that nonpolar and acidic residues have relatively strong enhancing and inhibiting effects, respectively, on peroxisome targeting. However, in contrast to the previous understanding, positive charges alone do not show the anticipated enhancing effect and that both the position and property of the residues within these patterns are important for peroxisome targeting. We further demonstrated that the three residues immediately upstream of the tripeptide are the core influencers, with a ‘basic-nonpolar-basic’ pattern serving as a strong and universal enhancing pattern for peroxisome targeting. These findings have significantly advanced our knowledge of the PTS1 domain in plants and likely other eukaryotic species as well. The principles and strategies employed in the present study may also be applied to deciphering auxiliary targeting signals for other organelles.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Real-time structural motif searching in proteins using an inverted index strategy

Biochemical and biological functions of proteins are the product of both the overall fold of the polypeptide chain, and, typically, structural motifs made up of smaller numbers of amino acids constituting a catalytic center or a binding site that may be remote from one another in amino acid sequence. Detection of such structural motifs can provide valuable insights into the function(s) of previously uncharacterized proteins. Technically, this remains an extremely challenging problem because of the size of the Protein Data Bank (PDB) archive. Existing methods depend on a clustering by sequence similarity and can be computationally slow. We have developed a new approach that uses an inverted index strategy capable of analyzing >170,000 PDB structures with unmatched speed. The efficiency of the inverted index method depends critically on identifying the small number of structures containing the query motif and ignoring most of the structures that are irrelevant. Our approach (implemented at motif.rcsb.org ) enables real-time retrieval and superposition of structural motifs, either extracted from a reference structure or uploaded by the user. Herein, we describe the method and present five case studies that exemplify its efficacy and speed for analyzing 3D structures of both proteins and nucleic acids.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Snekmer: a scalable pipeline for protein sequence fingerprinting based on amino acid recoding

Abstract Motivation The vast expansion of sequence data generated from single organisms and microbiomes has precipitated the need for faster and more sensitive methods to assess evolutionary and functional relationships between proteins. Representing proteins as sets of short peptide sequences (kmers) has been used for rapid, accurate classification of proteins into functional categories; however, this approach employs an exact-match methodology and thus may be limited in terms of sensitivity and coverage. We have previously used similarity groupings, based on the chemical properties of amino acids, to form reduced character sets and recode proteins. This amino acid recoding (AAR) approach simplifies the construction of protein representations in the form of kmer vectors, which can link sequences with distant sequence similarity and provide accurate classification of problematic protein families. Results Here, we describe Snekmer, a software tool for recoding proteins into AAR kmer vectors and performing either (i) construction of supervised classification models trained on input protein families or (ii) clustering for de novo determination of protein families. We provide examples of the operation of the tool against a set of nitrogen cycling families originally collected using both standard hidden Markov models and a larger set of proteins from Uniprot and demonstrate that our method accurately differentiates these sequences in both operation modes. Availability and implementation Snekmer is written in Python using Snakemake. Code and data used in this article, along with tutorial notebooks, are available at http://github.com/PNNL-CompBio/Snekmer under an open-source BSD-3 license. Supplementary information Supplementary data are available at Bioinformatics Advances online.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Multicellular magnetotactic bacteria are genetically heterogeneous consortia with metabolically differentiated cells

Consortia of multicellular magnetotactic bacteria (MMB) are currently the only known example of bacteria without a unicellular stage in their life cycle. Because of their recalcitrance to cultivation, most previous studies of MMB have been limited to microscopic observations. To study the biology of these unique organisms in more detail, we use multiple culture-independent approaches to analyze the genomics and physiology of MMB consortia at single-cell resolution. We separately sequenced the metagenomes of 22 individual MMB consortia, representing 8 new species, and quantified the genetic diversity within each MMB consortium. This revealed that, counter to conventional views, cells within MMB consortia are not clonal. Single consortia metagenomes were then used to reconstruct the species-specific metabolic potential and infer the physiological capabilities of MMB. To validate genomic predictions, we performed stable isotope probing (SIP) experiments and interrogated MMB consortia using fluorescence in situ hybridization (FISH) combined with nanoscale secondary ion mass spectrometry (NanoSIMS). By coupling FISH with bioorthogonal noncanonical amino acid tagging (BONCAT), we explored their in situ activity as well as variation of protein synthesis within cells. We demonstrate that MMB consortia are mixotrophic sulfate reducers and that they exhibit metabolic differentiation between individual cells, suggesting that MMB consortia are more complex than previously thought. These findings expand our understanding of MMB diversity, ecology, genomics, and physiology, as well as offer insights into the mechanisms underpinning the multicellular nature of their unique lifestyle.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Essential amino acids in the Plant-Conserved and Class-Specific Regions of cellulose synthases

The Plant-Conserved Region (P-CR) and the Class-Specific Region (CSR) are two plant-unique sequences in the catalytic core of cellulose synthases (CESAs) for which specific functions have not been established. Here, we used site-directed mutagenesis to replace amino acids and motifs within these sequences predicted to be essential for assembly and function of CESAs. We developed an in vivo method to determine the ability of mutated CesA1 transgenes to complement an Arabidopsis (Arabidopsis thaliana) temperature-sensitive root-swelling1 (rsw1) mutant. Replacement of a Cys residue in the CSR, which blocks dimerization in vitro, rendered the AtCesA1 transgene unable to complement the rsw1 mutation. Examination of the CSR sequences from 33 diverse angiosperm species showed domains of high-sequence conservation in a class-specific manner but with variation in the degrees of disorder, indicating a nonredundant role of the CSR structures in different CESA isoform classes. The Cys residue essential for dimerization was not always located in domains of intrinsic disorder. Expression of AtCesA1 transgene constructs, in which Pro417 and Arg453 were substituted for Ala or Lys in the coiled-coil of the P-CR, were also unable to complement the rsw1 mutation. Despite an expected role for Arg457 in trimerization of CESA proteins, AtCesA1 transgenes with Arg457Ala mutations were able to fully restore the wild-type phenotype in rsw1. Our data support that Cys662 within the CSR and Pro417 and Arg453 within the P-CR of Arabidopsis CESA1 are essential residues for functional synthase complex formation, but our data do not support a specific role for Arg457 in trimerization in native CESA complexes.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

Energy metric prediction for double insertion mutants via the RoseNet deep learning framework

Studying the structural and functional implications of protein mutations is an important task in computational biology and bioinformatics. We leverage our previously proposed RoseNet neural network architecture to predict energy metrics of proteins with double amino acid insertions or deletions (InDels). We train models on previously generated benchmark datasets containing the exhaustive double InDel mutations for three proteins, as well as an additional three proteins for which ∼145k random mutants, each with two InDels, have been generated. We expand on our previous work by evaluating three additional proteins and analyzing domain features that impact the prediction capabilities of RoseNet. These features include InDels into secondary structures and the solvent accessible surface area (SASA) scores of the residues. We uncover further evidence to support that RoseNet has a higher proficiency of generalizing to unseen residue combinations than unseen insertion positions. We also observe that RoseNet produces higher-quality predictions when inserting into a β-sheet over an α-helix. Additionally, when the insertions fall in an area of high SASA, RoseNet often displays better performance than inserting into areas of low SASA.

59 BASIC BIOLOGICAL SCIENCES↗