Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “protein sequence”

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 1,009 records · Page 56
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

The contributions from the progenitor genomes of the mesopolyploid Brassiceae are evolutionarily distinct but functionally compatible

The members of the tribe Brassiceae share a whole-genome triplication (WGT), and one proposed model for its formation is a two-step pair of hybridizations producing hexaploid descendants. However, evidence for this model is incomplete, and the evolutionary and functional constraints that drove evolution after the hexaploidy are even less understood. Here, we report a new genome sequence of Crambe hispanica , a species sister to most sequenced Brassiceae. Using this new genome and three others that share the hexaploidy, we traced the history of gene loss after the WGT using the Polyploidy Orthology Inference Tool (POInT). We confirm the two-step formation model and infer that there was a significant temporal gap between those two allopolyploidizations, with about a third of the gene losses from the first two subgenomes occurring before the arrival of the third. We also, for the 90,000 individual genes in our study, make parental subgenome assignments, inferring, with measured uncertainty, from which of the progenitor genomes of the allohexaploidy each gene derives. We further show that each subgenome has a statistically distinguishable rate of homoeolog losses. There is little indication of functional distinction between the three subgenomes: the individual subgenomes show no patterns of functional enrichment, no excess of shared protein-protein or metabolic interactions between their members, and no biases in their likelihood of having experienced a recent selective sweep. We propose a "mix and match" model of allopolyploidy, in which subgenome origin drives homoeolog loss propensities but where genes from different subgenomes function together without difficulty.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Adaptive language model training for molecular design

Abstract The vast size of chemical space necessitates computational approaches to automate and accelerate the design of molecular sequences to guide experimental efforts for drug discovery. Genetic algorithms provide a useful framework to incrementally generate molecules by applying mutations to known chemical structures. Recently, masked language models have been applied to automate the mutation process by leveraging large compound libraries to learn commonly occurring chemical sequences (i.e., using tokenization) and predict rearrangements (i.e., using mask prediction). Here, we consider how language models can be adapted to improve molecule generation for different optimization tasks. We use two different generation strategies for comparison, fixed and adaptive. The fixed strategy uses a pre-trained model to generate mutations; the adaptive strategy trains the language model on each new generation of molecules selected for target properties during optimization. Our results show that the adaptive strategy allows the language model to more closely fit the distribution of molecules in the population. Therefore, for enhanced fitness optimization, we suggest the use of the fixed strategy during an initial phase followed by the use of the adaptive strategy. We demonstrate the impact of adaptive training by searching for molecules that optimize both heuristic metrics, drug-likeness and synthesizability, as well as predicted protein binding affinity from a surrogate model. Our results show that the adaptive strategy provides a significant improvement in fitness optimization compared to the fixed pre-trained model, empowering the application of language models to molecular design tasks.

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

Decoding Cancer Variants of Unknown Significance for Helicase–Nuclease–RPA Complexes Orchestrating DNA Repair During Transcription and Replication

All tumors have DNA mutations, and a predictive understanding of those mutations could inform clinical treatments. However, 40% of the mutations are variants of unknown significance (VUS), with the challenge being to objectively predict whether a VUS is pathogenic and supports the tumor or whether it is benign. To objectively decode VUS, we mapped cancer sequence data and evolutionary trace (ET) scores onto crystallography and cryo-electron microscopy structures with variant impacts quantitated by evolutionary action (EA) measures. As tumors depend on helicases and nucleases to deal with transcription/replication stress, we targeted helicase–nuclease–RPA complexes: (1) XPB-XPD (within TFIIH), XPF-ERCC1, XPG, and RPA for transcription and nucleotide excision repair pathways and (2) BLM, EXO5, and RPA plus DNA2 for stalled replication fork restart. As validation, EA scoring predicts severe effects for most disease mutations, but disease mutants with low ET scores not only are likely destabilizing but also disrupt sophisticated allosteric mechanisms. For sites of disease mutations and VUS predicted to be severe, we found strong co-localization to ordered regions. Rare discrepancies highlighted the different survival requirements between disease and tumor mutations, as well as the value of examining proteins within complexes. In a genome-wide analysis of 33 cancer types, we found correlation between the number of mutations in each tumor and which pathways or functional processes in which the mutations occur, revealing different mutagenic routes to tumorigenesis. We also found upregulation of ancient genes including BLM, which supports a non-random and concerted cancer process: reversion to a unicellular, proliferation-uncontrolled, status by breaking multicellular constraints on cell division. Together, these genes and global analyses challenge the binary “driver” and “passenger” mutation paradigm, support a gradient impact as revealed by EA scoring from moderate to severe at a single gene level, and indicate reduced regulation as well as activity. The objective quantitative assessment of VUS scoring and gene overexpression in the context of functional interactions and pathways provides insights for biology, oncology, and precision medicine.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Activation of polycystin-1 signaling by binding of stalk-derived peptide agonists

Polycystin-1 (PC1) is the protein product of thePKD1gene whose mutation causes autosomal dominant Polycystic Kidney Disease (ADPKD). PC1 is an atypical G protein-coupled receptor (GPCR) with an autocatalytic GAIN domain that cleaves PC1 into extracellular N-terminal and membrane-embedded C-terminal (CTF) fragments. Recently, activation of PC1 CTF signaling was shown to be regulated by a stalk tethered agonist (TA), resembling the mechanism observed for adhesion GPCRs. Here, synthetic peptides of the first 9- (p9), 17- (p17), and 21-residues (p21) of the PC1 stalk TA were shown to re-activate signaling by a stalkless CTF mutant in human cell culture assays. Novel Peptide Gaussian accelerated molecular dynamics (Pep-GaMD) simulations elucidated binding conformations of p9, p17, and p21 and revealed multiple specific binding regions to the stalkless CTF. Peptide agonists binding to the TOP domain of PC1 induced close TOP-putative pore loop interactions, a characteristic feature of stalk TA-mediated PC1 CTF activation. Additional sequence coevolution analyses showed the peptide binding regions were consistent with covarying residue pairs identified between the TOP domain and the stalk TA. These insights into the structural dynamic mechanism of PC1 activation by TA peptide agonists provide an in-depth understanding that will facilitate the development of therapeutics targeting PC1 for ADPKD treatment.

Life Sciences & Biomedicine - Other Topics↗
DOE OSTI2024

Gcn2 structurally mimics and functionally repurposes the HisRS enzyme for the integrated stress response

Protein kinase Gcn2 attenuates protein synthesis in response to amino acid starvation while stimulating translation of a transcriptional activator of amino acid biosynthesis. Gcn2 activation requires a domain related to histidyl-tRNA synthetase (HisRS), the enzyme that aminoacylates tRNA His . While evidence suggests that deacylated tRNA binds the HisRS domain for kinase activation, ribosomal P-stalk proteins have been implicated as alternative activating ligands on stalled ribosomes. We report crystal structures of the HisRS domain ofChaetomium thermophilumGcn2 that reveal structural mimicry of both catalytic (CD) and anticodon-binding (ABD) domains, which in authentic HisRS bind the acceptor stem and anticodon loop of tRNA His . Elements for forming histidyl adenylate and aminoacylation are lacking, suggesting that Gcn2 HisRS was repurposed for kinase activation, consistent with mutations in the CD that dysregulate yeast Gcn2 function. Substituting conserved ABD residues well positioned to contact the anticodon loop or that form a conserved ABD–CD interface impairs Gcn2 function in starved cells. Mimicry in Gcn2 HisRS of two highly conserved structural domains for binding both ends of tRNA—each crucial for Gcn2 function—supports that deacylated tRNAs activate Gcn2 and exemplifies how a metabolic enzyme is repurposed to host new local structures and sequences that confer a novel regulatory function.

Science & Technology - Other Topics↗
DOE OSTI2024

Transcriptomic Analysis of the CAM Species Kalanchoë fedtschenkoi Under Low- and High-Temperature Regimes

Temperature stress is one of the major limiting environmental factors that negatively impact global crop yields. Kalanchoë fedtschenkoi is an obligate crassulacean acid metabolism (CAM) plant species, exhibiting much higher water-use efficiency and tolerance to drought and heat stresses than C 3 or C 4 plant species. Previous studies on gene expression responses to low- or high-temperature stress have been focused on C 3 and C 4 plants. There is a lack of information about the regulation of gene expression by low and high temperatures in CAM plants. To address this knowledge gap, we performed transcriptome sequencing (RNA-Seq) of leaf and root tissues of K. fedtschenkoi under cold (8 °C), normal (25 °C), and heat (37 °C) conditions at dawn (i.e., 2 h before the light period) and dusk (i.e., 2 h before the dark period). Our analysis revealed differentially expressed genes (DEGs) under cold or heat treatment in comparison to normal conditions in leaf or root tissue at each of the two time points. In particular, DEGs exhibiting either the same or opposite direction of expression change (either up-regulated or down-regulated) under cold and heat treatments were identified. In addition, we analyzed gene co-expression modules regulated by cold or heat treatment, and we performed in-depth analyses of expression regulation by temperature stresses for selected gene categories, including CAM-related genes, genes encoding heat shock factors and heat shock proteins, circadian rhythm genes, and stomatal movement genes. Our study highlights both the common and distinct molecular strategies employed by CAM and C 3 /C 4 plants in adapting to extreme temperatures, providing new insights into the molecular mechanisms underlying temperature stress responses in CAM species.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Kingdom-Wide Analysis of Fungal Protein-Coding and tRNA Genes Reveals Conserved Patterns of Adaptive Evolution

Protein-coding genes evolved codon usage bias due to the combined but uneven effects of adaptive and nonadaptive influences. Studies in model fungi agree on codon usage bias as an adaptation for fine-tuning gene expression levels; however, such knowledge is lacking for most other fungi. Our comparative genomics analysis of over 450 species supports codon usage and transfer RNAs (tRNAs) as coadapted for translation speed and this is most likely a realization of convergent evolution. Rather than drift, phylogenetic reconstruction inferred adaptive radiation as the best explanation for the variation of interspecific codon usage bias. Although the phylogenetic signals for individual codon and tRNAs frequencies are lower than expected by genetic drift, we found remarkable conservation of highly expressed genes being codon optimized for translation by the most abundant tRNAs, especially by inosine-modified tRNAs. As an application, we present a sequence-to-expression neural network that uses codons to reliably predict highly expressed transcripts. The kingdom Fungi, with over a million species, includes many key players in various ecosystems and good targets for biotechnology. Collectively, our results have implications for better understanding the evolutionary success of fungi, as well as informing the biosynthetic manipulation of fungal genes.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Structure and assembly of the diiron cofactor in the heme-oxygenase–like domain of the N-nitrosourea–producing enzyme SznF

In biosynthesis of the pancreatic cancer drug streptozotocin, the tridomain nonheme-iron oxygenase SznF hydroxylates N δ and N ω ' of N ω -methyl- L -arginine before oxidatively rearranging the triply modified guanidine to the N-methyl-N-nitrosourea pharmacophore. A previously published structure visualized the monoiron cofactor in the enzyme’s C-terminal cupin domain, which promotes the final rearrangement, but exhibited disorder and minimal metal occupancy in the site of the proposed diiron cofactor in the N-hydroxylating heme-oxygenase–like (HO-like) central domain. We leveraged our recent observation that the N-oxygenating µ-peroxodiiron(III/III) intermediate can form in the HO-like domain after the apo protein self-assembles its diiron(II/II) cofactor to solve structures of SznF with both of its iron cofactors bound. These structures of a biochemically validated member of the emerging heme-oxygenase–like diiron oxidase and oxygenase (HDO) superfamily with intact diiron cofactor reveal both the large-scale conformational change required to assemble the O 2 -reactive Fe 2 (II/II) complex and the structural basis for cofactor instability—a trait shared by the other validated HDOs. During cofactor (dis)assembly, a ligand-harboring core helix dynamically (un)folds. The diiron cofactor also coordinates an unanticipated Glu ligand contributed by an auxiliary helix implicated in substrate binding by docking and molecular dynamics simulations. The additional carboxylate ligand is conserved in another N-oxygenating HDO but not in two HDOs that cleave carbon–hydrogen and carbon–carbon bonds to install olefins. Among ~9,600 sequences identified bioinformatically as members of the emerging HDO superfamily, ~25% conserve this additional carboxylate residue and are thus tentatively assigned as N-oxygenases.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Molecular mechanism of secreted amyloid-β precursor protein in binding and modulating GABABR1a

A recent phenomenal study discovered that the extension domain of secreted amyloid-β precursor protein (sAPP) can bind to the intrinsically disordered sushi 1 domain of the γ-aminobutyric acid type B receptor subunit 1a (GABA B R1a) and modulate its synaptic transmission. The work provided an important structural foundation for the modulation of GABA B R1a; however, the detailed molecular interaction mechanism, crucial for future drug design, remains elusive. Here, we further investigated the dynamical interactions between sAPP peptides and the natively unstructured sushi 1 domain using all-atom molecular dynamics simulations, for both the 17-residue sAPP peptide (APP 17-mer) and its minimally active 9 residue segment (APP 9-mer). We then explored mutations of the APP 9-mer with rigorous free energy perturbation (FEP) calculations. Our in silico mutagenesis studies revealed key residues (D4, W6, and W7) responsible for the binding with the sushi 1 domain. More importantly, one double mutation based on different vertebrate APP sequences from evolution exhibited a stronger binding (ΔΔG = -1.91 ± 0.66 kcal mol -1 ), indicating a potentially enhanced GABA B R1a modulator. These large-scale simulations may provide new insights into the binding mechanism between sAPP and the sushi 1 domain, which could open new avenues in the development of future GABA B R1a-specific therapeutics.

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

Immune activation during Pseudomonas infection causes local cell wall remodeling and alters AGP accumulation

The plant cell boundary generally comprises constituents of the primary and secondary cell wall (CW) that are deposited sequentially during development. Although it is known that the CW acts as a barrier against phytopathogens and undergoes modifications to limit their invasion, the extent, sequence, and requirements of the pathogen-induced modifications of the CW components are still largely unknown, especially at the level of the polysaccharide fraction. To address this significant knowledge gap, we adopted the compatible Pseudomonas syringae–Arabidopsis thaliana system. We found that, despite systemic signaling actuation, Pseudomonas infection leads only to local CW modifications. Furthermore, by utilizing a combination of CW and immune signaling-deficient mutants infected with virulent or non-virulent bacteria, we demonstrated that the pathogen-induced changes in CW polysaccharides depend on the combination of pathogen virulence and the host's ability to mount an immune response. This results in a pathogen-driven accumulation of CW hexoses, such as galactose, and an immune signaling-dependent increase in CW pentoses, mainly arabinose, and xylose. Our analyses of CW changes during disease progression also revealed a distinct spatiotemporal pattern of arabinogalactan protein (AGP) deposition and significant modifications of rhamnogalacturonan sidechains. Furthermore, genetic analyses demonstrated a critical role of AGPs, specifically of the Arabinoxylan Pectin Arabinogalactan Protein1, in limiting pathogen growth. Collectively, our results provide evidence for the actuation of significant remodeling of CW polysaccharides in a compatible host-pathogen interaction, and, by identifying AGPs as critical elements of the CW in plant defense, they pinpoint opportunities to improve plants against diverse pathogens.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Metagenome-assembled genomes of phytoplankton microbiomes from the Arctic and Atlantic Oceans

Abstract Background Phytoplankton communities significantly contribute to global biogeochemical cycles of elements and underpin marine food webs. Although their uncultured genomic diversity has been estimated by planetary-scale metagenome sequencing and subsequent reconstruction of metagenome-assembled genomes (MAGs), this approach has yet to be applied for complex phytoplankton microbiomes from polar and non-polar oceans consisting of microbial eukaryotes and their associated prokaryotes. Results Here, we have assembled MAGs from chlorophyll a maximum layers in the surface of the Arctic and Atlantic Oceans enriched for species associations (microbiomes) with a focus on pico- and nanophytoplankton and their associated heterotrophic prokaryotes. From 679 Gbp and estimated 50 million genes in total, we recovered 143 MAGs of medium to high quality. Although there was a strict demarcation between Arctic and Atlantic MAGs, adjacent sampling stations in each ocean had 51–88% MAGs in common with most species associations between Prasinophytes and Proteobacteria . Phylogenetic placement revealed eukaryotic MAGs to be more diverse in the Arctic whereas prokaryotic MAGs were more diverse in the Atlantic Ocean. Approximately 70% of protein families were shared between Arctic and Atlantic MAGs for both prokaryotes and eukaryotes. However, eukaryotic MAGs had more protein families unique to the Arctic whereas prokaryotic MAGs had more families unique to the Atlantic. Conclusion Our study provides a genomic context to complex phytoplankton microbiomes to reveal that their community structure was likely driven by significant differences in environmental conditions between the polar Arctic and warm surface waters of the tropical and subtropical Atlantic Ocean.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

A Single Amino Acid Change in Nramp6 from Sedum Alfredii Hance Affects Cadmium Accumulation

SaNramp6 in Sedum alfredii encodes a membrane-localized metal transporter. We isolated the SaNramp6h allele from the hyperaccumulating ecotype (HE) of S. alfredii . When this allele was expressed in transgenic yeast and Arabidopsis thaliana , it enhanced their cadmium (Cd) sensitivity by increased Cd transport and accumulation. We isolated another allele, SaNramp6n , from a nonhyperaccumulating ecotype (NHE) of S. alfredii . Amino acid sequence comparisons revealed three amino acid differences between SaNramp6h and SaNramp6n. We investigated the Cd transport activity of the Nramp6 allele, and determined which residues are essential for the transport activity. We conducted structure-function analyses of SaNramp6 based on site-directed mutagenesis and functional assays of the mutants in yeast and Arabidopsis . The three residues that differed between SaNramp6h and SaNramp6n were mutated. Only the L157P mutation of SaNramp6h impaired Cd transport. The other mutations, S218N and T504A, did not affect the transport activity of SaNramp6h, indicating that these residues are not essential for metal selectivity. Transgenic plants overexpressing SaNramp6h L157P showed altered metal accumulation in shoots and roots. Our results suggest that the conserved site L157 is essential for the high metal transport activity of SaNramp6h. This information may be useful for limiting or increasing Cd transport by other plant natural resistance associated macrophage protein (NRAMP) proteins.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Combinatorial Glycomic Analyses to Direct CAZyme Discovery for the Tailored Degradation of Canola Meal Non-Starch Dietary Polysaccharides

Canola meal (CM), the protein-rich by-product of canola oil extraction, has shown promise as an alternative feedstuff and protein supplement in poultry diets, yet its use has been limited due to the abundance of plant cell wall fibre, specifically non-starch polysaccharides (NSP) and lignin. The addition of exogenous enzymes to promote the digestion of CM NSP in chickens has potential to increase the metabolizable energy of CM. We isolated chicken cecal bacteria from a continuous-flow mini-bioreactor system and selected for those with the ability to metabolize CM NSP. Of 100 isolates identified, Bacteroides spp. and Enterococcus spp. were the most common species with these capabilities. To identify enzymes specifically for the digestion of CM NSP, we used a combination of glycomics techniques, including enzyme-linked immunosorbent assay characterization of the plant cell wall fractions, glycosidic linkage analysis (methylation-GC-MS analysis) of CM NSP and their fractions, bacterial growth profiles using minimal media supplemented with CM NSP, and the sequencing and de novo annotation of bacterial genomes of high-efficiency CM NSP utilizing bacteria. The SACCHARIS pipeline was used to select plant cell wall active enzymes for recombinant production and characterization. This approach represents a multidisciplinary innovation platform to bioprospect endogenous CAZymes from the intestinal microbiota of herbivorous and omnivorous animals which is adaptable to a variety of applications and dietary polysaccharides.

glycome profiling↗
DOE OSTI2023

Eukaryotic RNA-guided endonucleases evolved from a unique clade of bacterial enzymes

Abstract RNA-guided endonucleases form the crux of diverse biological processes and technologies, including adaptive immunity, transposition, and genome editing. Some of these enzymes are components of insertion sequences (IS) in the IS200/IS605 and IS607 transposon families. Both IS families encode a TnpA transposase and a TnpB nuclease, an RNA-guided enzyme ancestral to CRISPR-Cas12s. In eukaryotes, TnpB homologs occur as two distinct types, Fanzor1s and Fanzor2s. We analyzed the evolutionary relationships between prokaryotic TnpBs and eukaryotic Fanzors, which revealed that both Fanzor1s and Fanzor2s stem from a single lineage of IS607 TnpBs with unusual active site arrangement. The widespread nature of Fanzors implies that the properties of this particular lineage of IS607 TnpBs were particularly suited to adaptation in eukaryotes. Biochemical analysis of an IS607 TnpB and Fanzor1s revealed common strategies employed by TnpBs and Fanzors to co-evolve with their cognate transposases. Collectively, our results provide a new model of sequential evolution from IS607 TnpBs to Fanzor2s, and Fanzor2s to Fanzor1s that details how genes of prokaryotic origin evolve to give rise to new protein families in eukaryotes.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

Bulk segregant and QTL analysis of maltose metabolism in Saccharomyces eubayanus

This study investigates natural variation in maltose utilization among Saccharomyces eubayanus strains isolated from different tree hosts. An F2 population was generated from a cross between strain QC18 (Quercus-associated) and strain CL467.1 (Nothofagus-associated). Bulk segregant analysis (BSA) was performed on the F2 offspring, and QTL mapping was conducted by whole-genome sequencing of two pools of 30 segregants each, representing the lowest and highest growth in high maltose concentrations. For QTL analysis, the parental strains QC18 and CL467.1 were also sequenced.

evolutionary plasticity↗
DOE OSTI2020

Metabolic versatility of the nitrite-oxidizing bacterium Nitrospira marina and its proteomic response to oxygen-limited conditions

The genus Nitrospira is the most widespread group of nitrite-oxidizing bacteria and thrives in diverse natural and engineered ecosystems. Nitrospira marina Nb-295 T was isolated from the ocean over 30 years ago; however, its genome has not yet been analyzed. Here, we investigated the metabolic potential of N. marina based on its complete genome sequence and performed physiological experiments to test genome-derived hypotheses. Our data confirm that N. marina benefits from additions of undefined organic carbon substrates, has adaptations to resist oxidative, osmotic, and UV light-induced stress and low dissolved p CO 2 , and requires exogenous vitamin B 12 . In addition, N. marina is able to grow chemoorganotrophically on formate, and is thus not an obligate chemolithoautotroph. We further investigated the proteomic response of N. marina to low (~5.6 µM) O 2 concentrations. The abundance of a potentially more efficient CO 2 -fixing pyruvate:ferredoxin oxidoreductase (POR) complex and a high-affinity cbb 3 -type terminal oxidase increased under O 2 limitation, suggesting a role in sustaining nitrite oxidation-driven autotrophy. This putatively more O 2 -sensitive POR complex might be protected from oxidative damage by Cu/Zn-binding superoxide dismutase, which also increased in abundance under low O 2 conditions. Furthermore, the upregulation of proteins involved in alternative energy metabolisms, including Group 3b [NiFe] hydrogenase and formate dehydrogenase, indicate a high metabolic versatility to survive conditions unfavorable for aerobic nitrite oxidation. In summary, the genome and proteome of the first marine Nitrospira isolate identifies adaptations to life in the oxic ocean and provides insights into the metabolic diversity and niche differentiation of NOB in marine environments.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Large Scale MD to Predict Epitope Regions in HIV Env [Slides]

Highly dense carbohydrates located on the surface of the HIV Env protein play a key role in immune evasion. Such evolutionary adaptation hampers any attempt to obtain a full mechanistic understanding of the role played by the glycans in protecting the virus against an effective immune response. Moreover, and due to their chemical variability, an accurate molecular understanding of the so called “glycan shield” is still limited by the lack of effective resolution of state-of-the-art experimental technics. Here, we have used extensive computational modelling in order to fill this gap, addressing the presence of a large glycan variability as observed experimentally. Based on an automated pipeline, we were able to assemble, set-up and simulate via Molecular dynamics hundreds of different glycosylated Env variants at nearly atomic resolution, recapitulating the glycosylation distributions observed experimentally. Results from these simulations were subjected to machine learning and very accurate prediction of simulation derived glycan shielding areas of each glycan as a function of static sequence features. Such predictive models of per-glycan shielding, incorporating both glycan dynamics and heterogeneity, were used to develop a novel sequence-based glycan shield mapping strategy. Parallel to these studies, we also developed an accurate machine learning approach to predict glycan heterogeneity data using sequence features and found good prediction accuracy.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Comparison of PsbQ and Psb27 in photosystem II provides insight into their roles

Photosystem II (PSII) catalyzes the oxidation of water at its active site that harbors a high-valent inorganic Mn 4 CaO x cluster called the oxygen-evolving complex (OEC). Extrinsic subunits generally serve to protect the OEC from reductants and stabilize the structure, but diversity in the extrinsic subunits exists between phototrophs. Recent cryo-electron microscopy experiments have provided new molecular structures of PSII with varied extrinsic subunits. We focus on the extrinsic subunit PsbQ, that binds to the mature PSII complex, and on Psb27, an extrinsic subunit involved in PSII biogenesis. PsbQ and Psb27 share a similar binding site and have a four-helix bundle tertiary structure, suggesting they are related. Here, we use sequence alignments, structural analyses, and binding simulations to compare PsbQ and Psb27 from different organisms. We find no evidence that PsbQ and Psb27 are related despite their similar structures and binding sites. Evolutionary divergence within PsbQ homologs from different lineages is high, probably due to their interactions with other extrinsic subunits that themselves exhibit vast diversity between lineages. This may result in functional variation as exemplified by large differences in their calculated binding energies. Psb27 homologs generally exhibit less divergence, which may be due to stronger evolutionary selection for certain residues that maintain its function during PSII biogenesis which is consistent with their more similar calculated binding energies between organisms. Previous experimental inconsistencies, low confidence binding simulations, and recent structural data suggest that Psb27 is likely to exhibit flexibility that may be an important characteristic of its activity. Furthermore, the analysis provides insight into the functions and evolution of PsbQ and Psb27, and an unusual example of proteins with similar tertiary structures and binding sites that probably serve different roles.

59 BASIC BIOLOGICAL SCIENCES↗