Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “protein complex structure prediction”

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 109 records · Page 6

Editorial: Structure and mechanism of microbial membrane active transporters

Membrane active transporters play essential roles in microbial physiology. They couple energy transduction to conformational changes that drive translocation of nutrients, substrates and ions, as well as molecular communication. The structure and function of microbial membrane active transporters are highly diverse. Typical examples include the primary active transporters in the ATP-binding cassette (ABC) superfamily (Thomas and Tampé, 2020; Davidson et al., 2008; Locher et al., 2002), the secondary active transporters in the Major Facilitator Superfamily (MFS) (Drew et al., 2021; Kaback and Guan, 2019), and the ligand-gated porins in the TonB-dependent transporter (TBDT) family (Klebba et al., 2021). As structural, proteogenomic, and computational methods advance, active transporters are increasingly recognized as dynamic molecular machines whose mechanisms can now be visualized and modeled with remarkable precision, building on decades of biochemical and biophysical discovery that established the foundations of this field. The transporter studies recruited in this Research Topic provide us with new insights into the field including structure-function of sugar transporters in yeast, structural prediction and classification of ABC complexes in Bacillus subtilis, Type VI Secretion System (T6SS) in Bacteroides fragilis, amino acids uptake in Escherichia coli and bacterial spore germination.

mechanism↗

Structure-Function and Nearest Neighbor Interactions of the Cytochrome b 6 f Complex

Test prediction from our (i) high resolution (2.5 Å) structure of cyanobacterial cytochrome b 6 f complex (Structure, 22, 1008-1015, 2014), and (ii) analysis of quinone/quinol traffic within the b 6 f complex (Biophys. J., 107, 1620-1628, 2014) concerning the rate-limiting step in the linear (non-cyclic) electron transport chain. It was concluded that it arises from access of plastoquinol, the electron/proton donor, to the high potential electron carriers, i. e., the Fe 2 S 2 and cytochrome f proteins. The directed motion of the quinol proceeds via a narrow channel on the lumen or electrochemically positive, p- (lumen) side of the b 6 f complex which was defined by crystal structures of the b 6 f complex with quinol analogue inhibitors (J. Mol. Biol., 370, 39-52, 2007).

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Deuteron off-resonance rotating frame relaxation for the characterization of slow motions in rotating and static solid-state proteins

Here, we demonstrate the feasibility of deuterium solid-state NMR off-resonance rotating frame relaxation measurements for studies of slow motions in biomolecular solids. The pulse sequence, which includes adiabatic pulses for magnetization alignment, is illustrated for static and magic-angle spinning conditions away from rotary resonances. We apply the measurements for three systems with selective deuterium labels at methyl groups: a) a model compound, Fluorenylmethyloxycarbonyl methionine-D 3 amino acid, for which the principles of the measurements and corresponding motional modeling based on rotameric interconversions are demonstrated; b) amyloid-β 1-40 fibrils labeled at a single alanine methyl group located in the disordered N-terminal domain. This system has been extensively studied in prior work and here serves as a test of the method for complex biological systems. The essential features of the dynamics consist of large-scale rearrangements of the disordered N-terminal domain and the conformational exchange between the free and bound forms of the domain, the latter one due to transient interactions with the structured core of the fibrils. and c) a 15-residue helical peptide which belongs to the predicted α-helical domain near the N-terminus of apolipoprotein B. The peptide is solvated with triolein and incorporates a selectively labeled leucine methyl groups. The method permits model refinement, indicating rotameric interconversions with a distribution of rate constants.

59 BASIC BIOLOGICAL SCIENCES↗

Probing a Silent Metal: A Combined X-ray Absorption and Emission Spectroscopic Study of Biologically Relevant Zinc Complexes

As the second most common transition metal in the human body, Zinc is of great interest in research but has few viable routes for its direct structural study in biological systems. Herein, Zn valence-to-core X-ray emission spectroscopy (VtC XES) and Zn K-edge X-ray absorption spectroscopy (XAS) are presented as a means to understand the local structure of zinc in biological systems through the application of these methods to a series of biologically relevant model complexes. Taken together, the Zn K-edge XAS and VtC XES provide a means to establish the ligand identity, local geometry, and metal-ligand bond lengths. Experimental results were supported by strongly correlated density-functional-theory-based calculations, which, coupled with the complementarity of the two spectroscopic approaches, further enables future applications to protein systems to be examined in a predictive fashion. In this regard, test calculations were performed to showcase the sensitivity of these techniques to realistic biological applications.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

An atypical BRCT–BRCT interaction with the XRCC1 scaffold protein compacts human DNA Ligase IIIα within a flexible DNA repair complex

The XRCC1–DNA ligase IIIα complex (XL) is critical for DNA single-strand break repair, a key target for PARP inhibitors in cancer cells deficient in homologous recombination. Here, we combined biophysical approaches to gain insights into the shape and conformational flexibility of the XL as well as XRCC1 and DNA ligase IIIα (LigIIIα) alone. Structurally-guided mutational analyses based on the crystal structure of the human BRCT–BRCT heterodimer identified the network of salt bridges that together with the N-terminal extension of the XRCC1 C-terminal BRCT domain constitute the XL molecular interface. Coupling size exclusion chromatography with small angle X-ray scattering and multiangle light scattering (SEC-SAXS–MALS), we determined that the XL is more compact than either XRCC1 or LigIIIα, both of which form transient homodimers and are highly disordered. The reduced disorder and flexibility allowed us to build models of XL particles visualized by negative stain electron microscopy that predict close spatial organization between the LigIIIα catalytic core and both BRCT domains of XRCC1. Together our results identify an atypical BRCT–BRCT interaction as the stable nucleating core of the XL that links the flexible nick sensing and catalytic domains of LigIIIα to other protein partners of the flexible XRCC1 scaffold.

59 BASIC BIOLOGICAL SCIENCES↗

Purification of the small mechanosensitive channel of Escherichia coli (MscS): the subunit structure, conduction, and gating characteristics in liposomes

The small mechanosensitive channel, MscS, is a part of the turgor-driven solute efflux system that protects bacteria from lysis in the event of osmotic downshift. It has been identified in Escherichia coli as a product of the orphan yggB gene, now called mscS (Levina et al., 1999, EMBO J. 18:1730). Here I show that that the isolated 31-kDa MscS protein is sufficient to form a functional mechanosensitive channel gated directly by tension in the lipid bilayer. MscS-6His complexes purified in the presence of octylglucoside and lipids migrate in a high-resolution gel-filtration column as particles of approximately 200 kDa. Consistent with that, the protein cross-linking patterns predict a hexamer. The channel reconstituted in soybean asolectin liposomes was activated by pressures of 20-60 mm Hg and displayed the same asymmetric I-V curve and slight anionic preference as in situ. At the same time, the single-channel conductance is proportional to the buffer conductivity in a wide range of salt concentrations. The rate of channel activation in response to increasing pressure gradient across the patch was slower than the rate of closure in response to decreasing steps of pressure gradient. Therefore, the open probability curves were recorded with descending series of pressures. Determination of the curvature of patches by video imaging permitted measurements of the channel activity as a function of membrane tension (gamma). Po(gamma) curves had the midpoint at 5.5 +/- 0.1 dyne/cm and gave estimates for the energy of opening DeltaG = 11.4 +/- 0.5 kT, and the transition-related area change DeltaA = 8.4 +/- 0.4 nm(2) when fitted with a two-state Boltzmann model. The correspondence between channel properties in the native and reconstituted systems is discussed.

NASA Discipline Cell Biology↗

Human XPG nuclease structure, assembly, and activities with insights for neurodegeneration and cancer from pathogenic mutations

Xeroderma pigmentosum group G (XPG) protein is both a functional partner in multiple DNA damage responses (DDR) and a pathway coordinator and structure-specific endonuclease in nucleotide excision repair (NER). Different mutations in the XPG gene ERCC5 lead to either of two distinct human diseases: Cancer-prone xeroderma pigmentosum (XP-G) or the fatal neurodevelopmental disorder Cockayne syndrome (XP-G/CS). To address the enigmatic structural mechanism for these differing disease phenotypes and for XPG’s role in multiple DDRs, here we determined the crystal structure of human XPG catalytic domain (XPGcat), revealing XPG-specific features for its activities and regulation. Furthermore, XPG DNA binding elements conserved with FEN1 superfamily members enable insights on DNA interactions. Notably, all but one of the known pathogenic point mutations map to XPGcat, and both XP-G and XP-G/CS mutations destabilize XPG and reduce its cellular protein levels. Mapping the distinct mutation classes provides structure-based predictions for disease phenotypes: Residues mutated in XP-G are positioned to reduce local stability and NER activity, whereas residues mutated in XP-G/CS have implied long-range structural defects that would likely disrupt stability of the whole protein, and thus interfere with its functional interactions. Combined data from crystallography, biochemistry, small angle X-ray scattering, and electron microscopy unveil an XPG homodimer that binds, unstacks, and sculpts duplex DNA at internal unpaired regions (bubbles) into strongly bent structures, and suggest how XPG complexes may bind both NER bubble junctions and replication forks. Collective results support XPG scaffolding and DNA sculpting functions in multiple DDR processes to maintain genome stability.

59 BASIC BIOLOGICAL SCIENCES↗

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↗

Drug discovery efforts at George Mason University

With over 39,000 students, and research expenditures in excess of $200 million, George Mason University (GMU) is the largest R1 (Carnegie Classification of very high research activity) university in Virginia. Mason scientists have been involved in the discovery and development of novel diagnostics and therapeutics in areas as diverse as infectious diseases and cancer. Below are highlights of the efforts being led by Mason researchers in the drug discovery arena. To enable targeted cellular delivery, and non-biomedical applications, Veneziano and colleagues have developed a synthesis strategy that enables the design of self-assembling DNA nanoparticles (DNA origami) with prescribed shape and size in the 10 to 100 nm range. The nanoparticles can be loaded with molecules of interest such as drugs, proteins and peptides, and are a promising new addition to the drug delivery platforms currently in use. The investigators also recently used the DNA origami nanoparticles to fine tune the spatial presentation of immunogens to study the impact on B cell activation. These studies are an important step towards the rational design of vaccines for a variety of infectious agents. To elucidate the parameters for optimizing the delivery efficiency of lipid nanoparticles (LNPs), Buschmann, Paige and colleagues have devised methods for predicting and experimentally validating the pKa of LNPs based on the structure of the ionizable lipids used to formulate the LNPs. These studies may pave the way for the development of new LNP delivery vehicles that have reduced systemic distribution and improved endosomal release of their cargo post administration. To better understand protein-protein interactions and identify potential drug targets that disrupt such interactions, Luchini and colleagues have developed a methodology that identifies contact points between proteins using small molecule dyes. The dye molecules noncovalently bind to the accessible surfaces of a protein complex with very high affinity, but are excluded from contact regions. When the complex is denatured and digested with trypsin, the exposed regions covered by the dye do not get cleaved by the enzyme, whereas the contact points are digested. The resulting fragments can then be identified using mass spectrometry. The data generated can serve as the basis for designing small molecules and peptides that can disrupt the formation of protein complexes involved in disease processes. For example, using peptides based on the interleukin 1 receptor accessory protein (IL-1RAcP), Luchini, Liotta, Paige and colleagues disrupted the formation of IL-1/IL-R/IL-1RAcP complex and demonstrated that the inhibition of complex formation reduced the inflammatory response to IL-1B. Working on the discovery of novel antimicrobial agents, Bishop, van Hoek and colleagues have discovered a number of antimicrobial peptides from reptiles and other species. DRGN-1, is a synthetic peptide based on a histone H1-derived peptide that they had identified from Komodo Dragon plasma. DRGN-1 was shown to disrupt bacterial biofilms and promote wound healing in an animal model. The peptide, along with others, is being developed and tested in preclinical studies. Other research by van Hoek and colleagues focuses on in silico antimicrobial peptide discovery, screening of small molecules for antibacterial properties, as well as assessment of diffusible signal factors (DFS) as future therapeutics. The above examples provide insight into the cutting-edge studies undertaken by GMU scientists to develop novel methodologies and platform technologies important to drug discovery.

59 BASIC BIOLOGICAL SCIENCES↗

Bipartite chromatin recognition by Hop1 from two diverged Holozoa

In meiosis, ploidy reduction is driven by a complex series of DNA breakage and recombination events between homologous chromosomes, orchestrated by meiotic HORMA domain proteins (HORMADs). Meiotic HORMADs possess a central chromatin binding region (CBR) whose architecture varies across eukaryotic groups. Here, we determine high-resolution crystal structures of the meiotic HORMAD CBR from two diverged aquatic Holozoa,Schistosoma mansoniandPatiria miniata, which reveal tightly associated plant homeodomain (PHD) and winged helix-turn-helix (wHTH) domains. We show that PHD–wHTH CBRs bind duplex DNA through their wHTH domains, and identify key residues that disrupt this interaction. Combining experimental and predicted structures, we show that the CBRs’ PHDs likely interact with the tail of histone H3, and may discriminate between unmethylated and trimethylated H3 lysine 4. Finally, we show that Holozoa Hop1 CBRs bind nucleosomes in vitro in a bipartite manner involving both the PHD and wHTH domain. Our data reveal how meiotic HORMADs with PHD–wHTH CBRs can bind chromatin and potentially discriminate between chromatin states to drive meiotic recombination to specific chromosomal regions.

Life Sciences & Biomedicine - Other Topics↗

Cholesterol modulates membrane elasticity via unified biophysical laws

Cholesterol and lipid unsaturation underlie a balance of opposing forces that features prominently in adaptive cell responses to diet and environmental cues. These competing factors have resulted in contradictory observations of membrane elasticity across different measurement scales, requiring chemical specificity to explain incompatible structural and elastic effects. Here, we demonstrate that – unlike macroscopic observations – lipid membranes exhibit a unified elastic behavior in the mesoscopic regime between molecular and macroscopic dimensions. Using nuclear spin techniques and computational analysis, we find that mesoscopic bending moduli follow a universal dependence on the lipid packing density regardless of cholesterol content, lipid unsaturation, or temperature. Our observations reveal that compositional complexity can be explained by simple biophysical laws that directly map membrane elasticity to molecular packing associated with biological function, curvature transformations, and protein interactions. The obtained scaling laws closely align with theoretical predictions based on conformational chain entropy and elastic stress fields. These findings provide unique insights into the membrane design rules optimized by nature and unlock predictive capabilities for guiding the functional performance of lipid-based materials in synthetic biology and real-world applications.

Kumarage, Teshani [Virginia Polytechnic Inst. and ↗

Isolation and characterization of IgG3 glycan-targeting antibodies with exceptional cross-reactivity for diverse viral families

Broadly reactive antibodies that target sequence-diverse antigens are of interest for vaccine design and monoclonal antibody therapeutic development because they can protect against multiple strains of a virus and provide a barrier to evolution of escape mutants. Using LIBRA-seq (linking B cell receptor to antigen specificity through sequencing) data for the B cell repertoire of an individual chronically infected with human immunodeficiency virus type 1 (HIV-1), we identified a lineage of IgG3 antibodies predicted to bind to HIV-1 Envelope (Env) and influenza A Hemagglutinin (HA). Two lineage members, antibodies 2526 and 546, were confirmed to bind to a large panel of diverse antigens, including several strains of HIV-1 Env, influenza HA, coronavirus (CoV) spike, hepatitis C virus (HCV) E protein, Nipah virus (NiV) F protein, and Langya virus (LayV) F protein. We found that both antibodies bind to complex glycans on the antigenic surfaces. Antibody 2526 targets the stem region of influenza HA and the N-terminal domain (NTD) region of SARS-CoV-2 spike. A crystal structure of 2526 Fab bound to mannose revealed the presence of a glycan-binding pocket on the light chain. Antibody 2526 cross-reacted with antigens from multiple pathogens and displayed no signs of autoreactivity. These features distinguish antibody 2526 from previously described glycan-reactive antibodies. Further study of this antibody class may aid in the selection and engineering of broadly reactive antibody therapeutics and can inform the development of effective vaccines with exceptional breadth of pathogen coverage.

Microbiology↗

Selective deuteration of an RNA:RNA complex for structural analysis using small-angle scattering

The structures of RNA:RNA complexes regulate many biological processes. Despite their importance, protein-free RNA:RNA complexes represent a tiny fraction of experimentally determined structures. Here, we describe a joint small-angle X-ray and neutron scattering (SAXS/SANS) approach to structurally interrogate conformational changes in a model RNA:RNA complex. Using SAXS, we measured the solution structures of the individual RNAs and of the overall RNA:RNA complex. With SANS, we demonstrate, as a proof of principle, that isotope labeling and contrast matching (CM) can be combined to probe the bound state structure of an RNA within a selectively deuterated RNA:RNA complex. Furthermore, we show that experimental scattering data can validate and improve predicted AlphaFold 3 RNA:RNA complex structures to reflect its solution structure. In conclusion, our work demonstrates that in silico modeling, SAXS, and CM-SANS can be used in concert to directly analyze conformational changes within RNAs when in complex, enhancing our understanding of RNA structure in functional assemblies.

HIV-1 dimerization initiation site↗

Feeling the Strain: Quantifying Ligand Deformation in Photosynthesis

Structural distortion of protein-bound ligands can play a critical role in enzyme function by tuning the electronic and chemical properties of the ligand molecule. However, quantifying these effects is difficult due to the limited resolution of protein structures and the difficulty of generating accurate structural restrains for non-protein ligands. Here, we seek to quantify these effects through a statistical analysis of ligand distortion in Chlorophyll (Chl) proteins (CP), where ring deformation is thought to play a role in energy and electron transfer. To assess the accuracy of ring-deformation estimates from available structural data, we take advantage of the C 2 symmetry of Photosystem II (PSII), comparing ring-deformation estimates for equivalent sites both within and between 113 distinct X-ray and Cryogenic electron microscopy (CryoEM) PSII structures. Significantly, we find that several deformation modes exhibit considerable variability in predictions, even for equivalent monomers, down to 2 Å resolution, to an extent that probably prevents their utilization in optical calculations. We further find that refinement restrains play a critical role in determining deformation values to resolution as low as 2 Å. However, for those modes that are well-resolved in the structural data, ring deformation in PSII is strongly conserved across all species tested, from cyanobacteria to algae. Furthermore, these results highlight both the opportunities and limitations inherent in the structure-based analysis of the bioenergetic and optical properties of CPs and other protein-ligand complexes.

14 SOLAR ENERGY↗

The regulation of the cell wall by glycosylphosphatidylinositol-anchored proteins in Arabidopsis

A polysaccharides-based cell wall covers the plant cell, shaping it and protecting it from the harsh environment. Cellulose microfibrils constitute the cell wall backbone and are embedded in a matrix of pectic and hemicellulosic polysaccharides and glycoproteins. Various environmental and developmental cues can regulate the plant cell wall, and diverse glycosylphosphatidylinositol (GPI)-anchored proteins participate in these regulations. GPI is a common lipid modification on eukaryotic proteins, which covalently tethers the proteins to the membrane lipid bilayer. Catalyzed by a series of enzymic complexes, protein precursors are post-translationally modified at their hydrophobic carboxyl-terminus in the endomembrane system and anchored to the lipid bilayer through an oligosaccharidic GPI modification. Ultimately, mature proteins reach the plasma membrane via the secretory pathway facing toward the apoplast and cell wall in plants. In Arabidopsis, more than three hundred GPI-anchored proteins (GPI-APs) have been predicted, and many are reported to be involved in diverse regulations of the cell wall. In this review, we summarize GPI-APs involved in cell wall regulation. GPI-APs are proposed to act as structural components of the cell wall, organize cellulose microfibrils at the cell surface, and during cell wall integrity signaling transduction. Besides regulating protein trafficking, the GPI modification is potentially governed by a GPI shedding system that cleaves and releases the GPI-anchored proteins from the plasma membrane into the cell wall.

59 BASIC BIOLOGICAL SCIENCES↗

Atypical Divergence of SARS-CoV-2 Orf8 from Orf7a within the Coronavirus Lineage Suggests Potential Stealthy Viral Strategies in Immune Evasion

Orf8, one of the most puzzling genes in the SARS lineage of coronaviruses, marks a unique and striking difference in genome organization between SARS-CoV-2 and SARS-CoV-1. Here, using sequence comparisons, we unequivocally reveal the distant sequence similarities between SARS-CoV-2 Orf8 with its SARS-CoV-1 counterparts and the X4-like genes of coronaviruses, including its highly divergent “paralog” gene Orf7a, whose product is a potential immune antagonist of known structure. Supervised sequence space walks unravel identity levels that drop below 10% and yet exhibit subtle conservation patterns in this novel superfamily, characterized by an immunoglobulin-like beta sandwich topology. We document the high accuracy of the sequence space walk process in detail and characterize the subgroups of the superfamily in sequence space by systematic annotation of gene and taxon groups. While SARS-CoV-1 Orf7a and Orf8 genes are most similar to bat virus sequences, their SARS-CoV-2 counterparts are closer to pangolin virus homologs, reflecting the fine structure of conservation patterns within the SARS-CoV-2 genomes. The divergence between Orf7a and Orf8 is exceptionally idiosyncratic, since Orf7a is more constrained, whereas Orf8 is subject to rampant change, a peculiar feature that may be related to hitherto-unknown viral infection strategies. Despite their common origin, the Orf7a and Orf8 protein families exhibit different modes of evolutionary trajectories within the coronavirus lineage, which might be partly attributable to their complex interactions with the mammalian host cell, reflected by a multitude of functional associations of Orf8 in SARS-CoV-2 compared to a very small number of interactions discovered for Orf7a.

60 APPLIED LIFE SCIENCES↗

Mutations of Omicron Variant at the Interface of the Receptor Domain Motif and Human Angiotensin-Converting Enzyme-2

The most recent Omicron variant of SARS-CoV-2 has caused global concern and anxiety. The only thing certain about this strain, with a large number of mutations in the spike protein, is that it spreads quickly, seems to evade immune defense, and mitigates the benefits of existing vaccines. Based on the ultra-large-scale ab initio computational modeling of the receptor binding motif (RBM) and the human angiotensin-converting enzyme-2 (ACE2) interface, we provide the details of the effect of Omicron mutations at the fundamental atomic scale level. In-depth analysis anchored in the novel concept of amino acid-amino acid bond pair units (AABPU) indicates that mutations in the Omicron variant are connected with (i) significant changes in the shape and structure of AABPU components, together with (ii) significant increase in the positive partial charge, which facilitates the interaction with ACE2. We have identified changes in bonding due to mutations in the RBM. The calculated bond order, based on AABPU, reveals that the Omicron mutations increase the binding strength of RBM to ACE2. Our findings correlate with and are instrumental to explain the current observations and can contribute to the prediction of next potential new variant of concern.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Autotrophic biofilms sustained by deeply sourced groundwater host diverse bacteria implicated in sulfur and hydrogen metabolism

Abstract Background Biofilms in sulfide-rich springs present intricate microbial communities that play pivotal roles in biogeochemical cycling. We studied chemoautotrophically based biofilms that host diverse CPR bacteria and grow in sulfide-rich springs to investigate microbial controls on biogeochemical cycling. Results Sulfide springs biofilms were investigated using bulk geochemical analysis, genome-resolved metagenomics, and scanning transmission X-ray microscopy (STXM) at room temperature and 87 K. Chemolithotrophic sulfur-oxidizing bacteria, including Thiothrix and Beggiatoa , dominate the biofilms, which also contain CPR Gracilibacteria, Absconditabacteria, Saccharibacteria, Peregrinibacteria, Berkelbacteria, Microgenomates, and Parcubacteria. STXM imaging revealed ultra-small cells near the surfaces of filamentous bacteria that may be CPR bacterial episymbionts. STXM and NEXAFS spectroscopy at carbon K and sulfur L 2,3 edges show that filamentous bacteria contain protein-encapsulated spherical elemental sulfur granules, indicating that they are sulfur oxidizers, likely Thiothrix . Berkelbacteria and Moranbacteria in the same biofilm sample are predicted to have a novel electron bifurcating group 3b [NiFe]-hydrogenase, putatively a sulfhydrogenase, potentially linked to sulfur metabolism via redox cofactors. This complex could potentially contribute to symbioses, for example, with sulfur-oxidizing bacteria such as Thiothrix that is based on cryptic sulfur cycling. One Doudnabacteria genome encodes adjacent sulfur dioxygenase and rhodanese genes that may convert thiosulfate to sulfite. We find similar conserved genomic architecture associated with CPR bacteria from other sulfur-rich subsurface ecosystems. Conclusions Our combined metagenomic, geochemical, spectromicroscopic, and structural bioinformatics analyses of biofilms growing in sulfide-rich springs revealed consortia that contain CPR bacteria and sulfur-oxidizing Proteobacteria, including Thiothrix , and bacteria from a new family within Beggiatoales. We infer roles for CPR bacteria in sulfur and hydrogen cycling.

59 BASIC BIOLOGICAL SCIENCES↗