Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “protein function”

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,045 records · Page 58
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Recombinant production of a functional SARS-CoV-2 spike receptor binding domain in the green algae Chlamydomonas reinhardtii

Recombinant production of viral proteins can be used to produce vaccine antigens or reagents to identify antibodies in patient serum. Minimally, these proteins must be correctly folded and have appropriate post-translation modifications. Here we report the production of the SARS-CoV-2 spike protein Receptor Binding Domain (RBD) in the green algae Chlamydomonas . RBD fused to a fluorescent reporter protein accumulates as an intact protein when targeted for ER-Golgi retention or secreted from the cell, while a chloroplast localized version is truncated. The ER-retained RBD fusion protein was able to bind the human ACE2 receptor, the host target of SARS-CoV-2, and was specifically out-competed by mammalian cell-produced recombinant RBD, suggesting that the algae produced proteins are sufficiently post-translationally modified to act as authentic SARS-CoV-2 antigens. Because algae can be grown at large scale very inexpensively, this recombinant protein may be a low cost alternative to other expression platforms.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Proton transfer reactions in donor site mutants of ESR, a retinal protein from Exiguobacterium sibiricum

Light-driven proton transport by microbial retinal proteins such as archaeal bacteriorhodopsin involves carboxylic residues as internal proton donors to the catalytic center which is a retinal Schiff base (SB). The proton donor, Asp96 in bacteriorhodopsin, supplies a proton to the transiently deprotonated Schiff base during the photochemical cycle. Subsequent proton uptake resets the protonated state of the donor. Here this two step process became a distinctive signature of retinal based proton pumps. Similar steps are observed also in many natural variants of bacterial proteorhodopsins and xanthorhodopsins where glutamic acid residues serve as a proton donor. Recently, however, an exception to this rule was found. A retinal protein from Exiguobacterium sibiricum, ESR, contains a Lys residue in place of Asp or Glu, which facilitates proton transfer from the bulk to the SB. Lys96 can be functionally replaced with the more common donor residues, Asp or Glu. Proton transfer to the SB in the mutants containing these replacements (K96E and K96D/A47T) is much faster than in the proteins lacking the proton donor (K96A and similar mutants), and in the case of K96D/A47T, comparable with that in the wild type, indicating that carboxylic residues can replace Lys96 as proton donors in ESR. We show here that there are important differences in the functioning of these residues in ESR from the way Asp96 functions in bacteriorhodopsin. Reprotonation of the SB and proton uptake from the bulk occur almost simultaneously during the M to N transition (as in the wild type ESR at neutral pH), whereas in bacteriorhodopsin these two steps are well separated in time and occur during the M to N and N to O transitions, respectively.

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

Structural characterization of high-protein system through ultra-small and small-angle X-ray scattering

High-protein systems exhibit a hierarchical structure consisting of interconnected multiscale length assemblies. Recent investigations are directed toward identifying these structural units under various processing and environmental conditions to establish structure–function correlations. Ultra-small and small-angle X-ray scattering (USAXS/SAXS) has become a crucial tool for characterizing the structures of proteins and their clusters/aggregates, ranging from nanometers to micrometers, with minimal disruption to their original state. Here, this review first describes the facilities and principles of X-ray scattering, followed by discussions on the analysis of scattering data, including the interpretation of fitting models. It then delves into the main applications of USAXS/SAXS in plant and dairy protein-rich systems. Future strategies to enhance the utilization of scattering techniques for elucidating the structure of high-protein systems are also included.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2022

Allosteric control of electron transfer in nitrogenase and nitrogenase-like enzymes (Final Report)

More than a third of all known proteins to date are either oxidoreductases or metalloproteins. The diversity of metal cofactors and redox centers in these proteins make them attractive candidates for engineering protein-based biocatalysts, biosensors, and biofuel cells, among other necessary tools for the bioeconomy. Many of these enzymes function as large multi-subunit complexes arranged as symmetric or asymmetric assemblies, and the purpose of this architecture is unclear. The goal of the project was to address this fundamental question by investigating the mechanism of action of two oxidoreductases: Nitrogenase and DPOR (Dark-Operative Protochlorophyllide Oxidoreductase, a nitrogenase-like oxidoreductase). Nitrogenase and DPOR play key roles in nitrogen fixation and photosynthesis, processes essential to sustaining life. Both proteins are arranged as two symmetrical halves with each half capable of substrate reduction. For more than half a century, substrate reduction by nitrogenase has been presumed to occur independently within the two identical halves. Funded by this DOE grant we have now shown that both nitrogenase and DPOR function asymmetrically. Electron transfer in one half suppresses activity in the other. Thus, we have demonstrated that the two halves do not function in an independent manner; rather, catalytic events in one-half allosterically suppress activity and electron transfer (ET) in the other. This alternating-site mechanism has significant implications on our current understanding of substrate reduction by these enzymes and offers a new mechanistic paradigm for the reactions they catalyze. Our work now provides a mechanistic map of the how the two symmetric catalytic halves work together to establish and maintain the sequential-site mechanism through multiple ET cycles to achieve substrate reduction. There is substantive need for bio-inspired tools that enable reduction of nitrogen and other substrates, as well as photosynthesis, for energy and food production. Understanding the workings of enzymes like nitrogenase and DPOR that mediate multi-electron redox chemistry is vital for designing such synthetic tools. The timing of our research coincides with the biotechnological breakthrough of light-driven artificial N2 reduction by nitrogenase, achieved using the native enzyme with both halves arranged in a mirror-like architecture as catalyst. Together with such discoveries, our research now reveal unknown facets of the reaction mechanism and hold promise for designing more effective catalysts.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Selected humanization of yeast U1 snRNP leads to global suppression of pre-mRNA splicing and mitochondrial dysfunction in the budding yeast

The recognition of the 5′ splice site (5′ ss) is one of the earliest steps of pre-mRNA splicing. To better understand, the mechanism and regulation of 5′ ss recognition, we selectively humanized components of the yeast U1 (yU1) snRNP to reveal the function of these components in 5′ ss recognition and splicing. We targeted U1C and Luc7, two proteins that interact with and stabilize the yU1 snRNA and the 5′ ss RNA duplex. We replaced the zinc-finger (ZnF) domain of yeast U1C (yU1C) with its human counterpart, which resulted in a cold-sensitive growth phenotype and moderate splicing defects. We next added an auxin-inducible degron to yeast Luc7 (yLuc7) protein (to mimic the lack of Luc7Ls in human U1 snRNP). We found that Luc7-depleted yU1 snRNP resulted in the concomitant loss of Prp40 and Snu71 (two other essential yU1 snRNP proteins), and further biochemical analyses suggest a model of how these three proteins interact with each other in the U1 snRNP. The loss of these proteins resulted in a significant growth retardation accompanied by a global suppression of pre-mRNA splicing. The splicing suppression led to mitochondrial dysfunction as revealed by a release of Fe 2+ into the growth medium and an induction of mitochondrial reactive oxygen species. Together, these observations indicate that the human U1C ZnF can substitute that of yeast, Luc7 is essential for the incorporation of the Luc7–Prp40–Snu71 trimer into yU1 snRNP, and splicing plays a major role in the regulation of mitochondrial function in yeast.

Biochemistry & Molecular Biology↗
DOE OSTI2023

A single amino acid change led to structural and functional differentiation of PvHd1 to control flowering in switchgrass

Abstract Switchgrass, a forage and bioenergy crop, occurs as two main ecotypes with different but overlapping ranges of adaptation. The two ecotypes differ in a range of characteristics, including flowering time. Flowering time determines the duration of vegetative development and therefore biomass accumulation, a key trait in bioenergy crops. No causal variants for flowering time differences between switchgrass ecotypes have, as yet, been identified. In this study, we mapped a robust flowering time quantitative trait locus (QTL) on chromosome 4K in a biparental F2 population and characterized the flowering-associated transcription factor gene PvHd1, an ortholog of CONSTANS in Arabidopsis and Heading date 1 in rice, as the underlying causal gene. Protein modeling predicted that a serine to glycine substitution at position 35 (p.S35G) in B-Box domain 1 greatly altered the global structure of the PvHd1 protein. The predicted variation in protein compactness was supported in vitro by a 4 °C shift in denaturation temperature. Overexpressing the PvHd1-p.35S allele in a late-flowering CONSTANS-null Arabidopsis mutant rescued earlier flowering, whereas PvHd1-p.35G had a reduced ability to promote flowering, demonstrating that the structural variation led to functional divergence. Our findings provide us with a tool to manipulate the timing of floral transition in switchgrass cultivars and, potentially, expand their cultivation range.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Structural insights into the assembly and the function of the plant oxidative phosphorylation system

One of the key functions of mitochondria is the production of ATP to support cellular metabolism and growth. The last step of mitochondrial ATP synthesis is performed by the oxidative phosphorylation (OXPHOS) system, an ensemble of protein complexes embedded in the inner mitochondrial membrane. In the last 25 yr, many structures of OXPHOS complexes and supercomplexes have been resolved in yeast, mammals, and bacteria. However, structures of plant OXPHOS enzymes only became available very recently. In this review, we highlight the plant-specific features revealed by the recent structures and discuss how they advance our understanding of the function and assembly of plant OXPHOS complexes. Here we also propose new hypotheses to be tested and discuss older findings to be re-evaluated. Further biochemical and structural work on the plant OXPHOS system will lead to a deeper understanding of plant respiration and its regulation, with significant agricultural, environmental, and societal implications.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Comparative analysis of KRAS4a and KRAS4b splice variants reveals distinctive structural and functional properties

KRAS, the most frequently mutated oncogene in human cancer, produces two isoforms, KRAS4a and KRAS4b, through alternative splicing. These isoforms differ in exon 4, which encodes the final 15 residues of the G-domain and hypervariable regions (HVRs), vital for trafficking and membrane localization. While KRAS4b has been extensively studied, KRAS4a has been largely overlooked. Our multidisciplinary study compared the structural and functional characteristics of KRAS4a and KRAS4b, revealing distinct structural properties and thermal stability. Position 151 influences KRAS4a’s thermal stability, while position 153 affects binding to RAF1 CRD protein. Nuclear magnetic resonance analysis identified localized structural differences near sequence variations and provided a solution-state conformational ensemble. Notably,KRAS4aexhibits substantial transcript abundance in bile ducts, liver, and stomach, with transcript levels approachingKRAS4bin the colon and rectum. Functional disparities were observed in full-length KRAS variants, highlighting the impact of HVR variations on interaction with trafficking proteins and downstream effectors like RAF and PI3K within cells.

Science & Technology - Other Topics↗
DOE OSTI2021

The encapsulin from Thermotoga maritima is a flavoprotein with a symmetry matched ferritin-like cargo protein

Abstract Bacterial nanocompartments, also known as encapsulins, are an emerging class of protein-based ‘organelles’ found in bacteria and archaea. Encapsulins are virus-like icosahedral particles comprising a ~ 25–50 nm shell surrounding a specific cargo enzyme. Compartmentalization is thought to create a unique chemical environment to facilitate catalysis and isolate toxic intermediates. Many questions regarding nanocompartment structure–function remain unanswered, including how shell symmetry dictates cargo loading and to what extent the shell facilitates enzymatic activity. Here, we explore these questions using the model Thermotoga maritima nanocompartment known to encapsulate a redox-active ferritin-like protein. Biochemical analysis revealed the encapsulin shell to possess a flavin binding site located at the interface between capsomere subunits, suggesting the shell may play a direct and active role in the function of the encapsulated cargo. Furthermore, we used cryo-EM to show that cargo proteins use a form of symmetry-matching to facilitate encapsulation and define stoichiometry. In the case of the Thermotoga maritima encapsulin, the decameric cargo protein with fivefold symmetry preferentially binds to the pentameric-axis of the icosahedral shell. Taken together, these observations suggest the shell is not simply a passive barrier—it also plays a significant role in the structure and function of the cargo enzyme.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Accurate Free Energies for Complex Condensed-Phase Reactions Using an Artificial Neural Network Corrected DFTB/MM Methodology

Semiempirical methods like density functional tight-binding (DFTB) allow extensive phase space sampling, making it possible to generate free energy surfaces of complex reactions in condensed-phase environments. Such a high efficiency often comes at the cost of reduced accuracy, which may be improved by developing a specific reaction parametrization (SRP) for the particular molecular system. Thiol–disulfide exchange is a nucleophilic substitution reaction that occurs in a large class of proteins. Its proper description requires a high-level ab initio method, while DFT-GAA and hybrid functionals were shown to be inadequate, and so is DFTB due to its DFT-GGA descent. We develop an SRP for thiol–disulfide exchange based on an artificial neural network (ANN) implementation in the DFTB+ software and compare its performance to that of a standard SRP approach applied to DFTB. Furthermore, as an application, we use both new DFTB-SRP as components of a QM/MM scheme to investigate thiol–disulfide exchange in two molecular complexes: a solvated model system and a blood protein. Demonstrating the strengths of the methodology, highly accurate free energy surfaces are generated at a low cost, as the augmentation of DFTB with an ANN only adds a small computational overhead.

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

Examining polymer‐protein biophysical interactions with small‐angle x‐ray scattering and quartz crystal microbalance with dissipation

Abstract Polymer‐protein hybrids can be deployed to improve protein solubility and stability in denaturing environments. While previous work used robotics and active machine learning to inform new designs, further biophysical information is required to ascertain structure–function behavior. Here, we show the value of tandem small‐angle x‐ray scattering (SAXS) and quartz crystal microbalance with dissipation (QCMD) experiments to reveal detailed polymer‐protein interactions with horseradish peroxidase (HRP) as a test case. Of particular interest was the process of polymer‐protein complex formation under thermal stress whereby SAXS monitors formation in solution while QCMD follows these dynamics at an interface. The radius of gyration ( R g ) of the protein as measured by SAXS does not change significantly in the presence of polymer under denaturing conditions, but thickness and dissipation changes were observed in QCMD data. SAXS data with and without thermal stress were utilized to create bead models of the potential complexes and denatured enzyme, and each model fit provided insight into the degree of interactions. Additionally, QCMD data demonstrated that HRP deforms by spreading upon surface adsorption at low concentration as shown by longer adsorption times and smaller frequency shifts. In contrast, thermally stressed and highly inactive HRP had faster adsorption kinetics. The combination of SAXS and QCMD serves as a framework for biophysical characterization of interactions between proteins and polymers which could be useful in designing polymer‐protein hybrids.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2021

Orthogonal Degron System for Controlled Protein Degradation in Cyanobacteria

Synechococcus elongatus PCC 7942 is a model cyanobacterium for study of the circadian clock, photosynthesis, and bioproduction of chemicals, yet nearly 40% of its gene identities and functions remain unknown, in part due to limitations of the existing genetic toolkit. While classical techniques for the study of genes (e.g., deletion or mutagenesis) can yield valuable information about the absence of a gene and its associated protein, there are limits to these approaches, particularly in the study of essential genes. In this research, we developed a tool for inducible degradation of target proteins in S. elongatus by adapting a method using degron tags from the Mesoplasma florum transfer-messenger RNA (tmRNA) system. We observed that M. florum lon protease can rapidly degrade exogenous and native proteins tagged with the cognate sequence within hours of induction. We used this system to inducibly degrade the essential cell division factor, FtsZ, as well as shell protein components of the carboxysome. Our results have implications for carboxysome biogenesis and the rate of carboxysome turnover during cell growth. Lon protease control of proteins offers an alternative approach for the study of essential proteins and protein dynamics in cyanobacteria.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Preclinical Risk Evaluation of Normal Tissue Injury With Novel Radiosensitizers

Genotoxic damage induced by radiation triggers a highly coordinated DNA damage response, and molecular inhibitors of key nodes within this complex response network can profoundly enhance the antitumor efficacy of radiation. This is especially true for drugs targeting the catalytic subunit of DNA-dependent protein kinase, which is a core component of the nonhomologous end-joining DNA repair pathway, and ataxia telangiectasia mutated, which coordinates cell cycle arrest, apoptosis, and DNA repair functionalities after radiation exposure. Unlike the more modest in vitro radiosensitizing effects seen with classic sensitizing agents such as cisplatin, 5-fluorouracil, or taxanes, DNA-dependent protein kinase or ataxia telangiectasia mutated inhibitors provide much more robust sensitizing effects in vitro, as might be anticipated from targeting these key DNA repair modulators. However, patients with homozygous inactivating mutations of ataxia telangiectasia mutated or mice with homozygous defects in DNA-dependent protein kinase (severe combined immunodeficiency) have profoundly enhanced acute normal tissue radiation reactions. Therefore, there is significant potential that the combination of small molecule inhibitors of these kinases with radiation could cause similar dose-limiting acute normal tissue toxicities. Similarly, although less understood, inhibition of these DNA repair response pathways could markedly increase the risk of late radiation toxicities. Because these potent radiosensitizers could be highly useful to improve local control of otherwise radiation-resistant tumors, understanding the potential for elevated risks of radiation injury is essential for optimizing therapeutic ratio and developing safe and informative clinical trials. In this review, we will discuss 2 straightforward models to assess the potential for enhanced mucosal toxicity in the oral cavity and small intestine established in our laboratories. We also will discuss similar strategies for evaluating potential drug–radiation interactions with regard to increased risks of debilitating late effects.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2023

Mechanistic Insights into Peptide Binding and Deactivation of an Adhesion G Protein-Coupled Receptor

Adhesion G protein-coupled receptors (ADGRGs) play critical roles in the reproductive, neurological, cardiovascular, and endocrine systems. In particular, ADGRG2 plays a significant role in Ewing sarcoma cell proliferation, parathyroid cell function, and male fertility. In 2022, a cryo-EM structure was reported for the active ADGRG2 bound by an optimized peptide agonist IP15 and the Gs protein. The IP15 peptide agonist was also modified to antagonists 4PH-E and 4PH-D with mutations of the 4PH residue to Glu and Asp, respectively. However, experimental structures of inactive antagonist-bound ADGRs remain to be resolved, and the activation mechanism of ADGRs such as ADGRG2 is poorly understood. Here, we applied Gaussian accelerated molecular dynamics (GaMD) simulations to probe conformational dynamics of the agonist- and antagonist-bound ADGRG2. By performing GaMD simulations, we were able to identify important low-energy conformations of ADGRG2 in the active, intermediate, and inactive states, as well as explore the binding conformations of each peptide. Moreover, our simulations revealed critical peptide-receptor residue interactions during the deactivation of ADGRG2. In conclusion, through GaMD simulations, we uncovered mechanistic insights into peptide (agonist and antagonist) binding and deactivation of the ADGRG2. These findings will potentially facilitate rational design of new peptide modulators of ADGRG2 and other ADGRs.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2026

Novel Cell-Type-Specific Drought-Responsive Proteins in Root Tips of Field-Grown Perennial Switchgrass

The root-tip region of plants, including the root cap, forms the most basal terminal of the root and exhibits a high degree of cellular complexity in terms of morphology, cytological function, and interaction with environmental cues in the soil. Cells in this region follow a developmental trajectory, transitioning from stem cells to meristematic cells, and ultimately to fully differentiated cell types. However, our understanding of root-tip cell-type specific proteomic responses to abiotic stresses, such as drought, particularly under field conditions, remains limited. This study aimed to identify spatially resolved, cell type-specific proteomes in switchgrass (Panicum virgatum) root tips under drought stress. Root tips were collected from seven-year-old, field-grown switchgrass ‘Alamo’ plants excavated under both well-watered and long-term drought conditions. Cell type-specific proteins were identified using laser capture microdissection (LCM) coupled with nanoPOTS (Nanodroplet Processing in One Pot for Trace Samples) and nano-LC-MS proteomics analysis. Five distinct cell types were targeted: (1) cells in the quiescent center and stem cell niche (QuC), (2) protodermal epidermal cells (PEC) in the meristematic zone, (3) epidermal cells in the transition and elongation zones above the root cap (Epi), (4) peripheral root cap cells (PRC), forming 2–3 layers below the PEC and 1–2 layers above the root border cells, and (5) columella root cap cells (Col) comprising of the columella initials and a single underlying layer of cells undergoing active growth. Principal component analysis (PCA) revealed clear separation among the five targeted cell types, confirming distinct proteomic profiles. Proteins predominantly enriched in each cell type were linked to distinct cellular functions, with QuC cells showing involvement in chromosomal behavior, DNA replication, and mitosis—key processes for stem cell niche regulation. Drought stress resulted in alterations of proteostasis, as evidenced by significant decreases in ribosomal proteins and increases in protein synthesis inhibitors. Moreover, drought stress induced unique cell-type–specific proteins involved in phytohormone biosynthesis and signaling pathways, including auxin, cytokinin, and jasmonic acid. In particular, QuC cells were more highly enriched in proteins associated with DNA repair and mitotic processes. Metabolic pathways related to amino acids, carbohydrates, and lipids were differentially affected in a cell-type–dependent manner, whereas general stress-responsive proteins exhibited consistent changes across all five cell types. Overall, this study provides unique spatially resolved, cell-type-specific proteomic profiles in root tips, representing a significant advancement in our understanding of the cellular mechanisms underlying plant responses to drought stress in natural field conditions.

perennial grass↗
DOE OSTI2020

Generation of camelina mid-oleic acid seed oil by identification and stacking of fatty acid biosynthetic mutants

Modifying oilseeds to obtain a desired fatty acid composition is often necessary to enable use as feedstocks for specific applications such as food processing, biofuels, or biolubricants. A mutant population of camelina (Camelina sativa), an emerging specialty oilseed crop, was screened by high-throughput gas chromatography for lines with altered seed oil fatty acid composition. By leveraging knowledge of fatty acid synthesis in Arabidopsis thaliana, mutations in orthologs of FATTY ACID ELONGASE1 (FAE1), FATTY ACID DESATURASE2 (FAD2), FATTY ACID DESATURASE3 (FAD3), and β-KETO-ACYL-ACP SYNTHASE II (KASII; FAB1) were identified. The mutations altered conserved amino acid residues in the encoded proteins. The ability of the mutations in FAE1, FAD2 and FAD3 to affect enzyme function was demonstrated by comparing in vivo activities of wild-type and mutant alleles in yeast. In addition, expression of wild-type cDNA in camelina complemented fatty acid phenotypes of these mutants. As camelina has a hexaploid genome, the effect of a mutation in one of the three homeologs for each gene resulted in no or less severe growth phenotypes compared to similar mutations in Arabidopsis. Mid-oleic oils with nearly 40 % oleic acid and reduced very long-chain (≤C20) fatty acid content were obtained by crossing to obtain a fae1c/fad2a/fae1a/fad3a quadruple mutant. Little effect on total seed oil content was observed in the stacked mutant line. Finally, the resulting mid-oleic acid oil had improved oxidative stability due to reductions in polyunsaturated fatty acid content, increasing its utility for biofuels and other applications.

09 BIOMASS FUELS↗
DOE OSTI2022

Characterization of Heparin’s Conformational Ensemble by Molecular Dynamics Simulations and Nuclear Magnetic Resonance Spectroscopy

Heparin is a highly charged, polysulfated polysaccharide and serves as an anticoagulant. Heparin binds to multiple proteins throughout the body, suggesting a large range of potential therapeutic applications. Although its function has been characterized in multiple physiological contexts, heparin’s solution conformational dynamics and structure–function relationships are not fully understood. Molecular dynamics (MD) simulations facilitate the analysis of a molecule’s underlying conformational ensemble, which then provides important information necessary for understanding structure–function relationships. However, for MD simulations to afford meaningful results, they must both provide adequate sampling and accurately represent the energy properties of a molecule. The aim of this study is to compare heparin’s conformational ensemble using two well-developed force fields for carbohydrates, known as GLYCAM06 and CHARMM36, using replica exchange molecular dynamics (REMD) simulations, and to validate these results with NMR experiments. In this work, the anticoagulant sequence, an ultra-low-molecular-weight heparin, known as Arixtra (fondaparinux, sodium), was simulated with both parameter sets. The results suggest that GLYCAM06 matches experimental nuclear magnetic resonance three-bond J-coupling values measured for Arixtra better than CHARMM36. In addition, NOESY and ROESY experiments suggest that Arixtra is very flexible in the sub-millisecond time scale and does not adopt a unique structure at 25 C. Moreover, GLYCAM06 affords a much more dynamic conformational ensemble for Arixtra than CHARMM36.

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

Molecular Modeling and Molecular Dynamics Simulation of a Packed and Intact Bacterial Microcompartment

Bacterial microcompartments (BMCs) are protein-bound organelles found in some bacteria which encapsulate enzymes for enhanced catalytic activity. These compartments spatially sequester enzymes within semipermeable shell proteins and are packed full of enzyme cargoes and metabolites as they fulfill their function. Coupling together recent SAXS and proteomics work, it is possible to develop molecular models for these microcompartments and interrogate enzyme and metabolite dynamics within. Our primary goal of this study is to quantify the permeability of metabolite glyceraldehyde-3-phosphate (G3P) and dihydroxyacetone phosphate (DHAP) across the BMC shell through classical molecular dynamics simulation. The Haliangium ochraceum model of BMC shell (PDB: 6MZX) was used to model an intact BMC of approximately 10 million atoms. Working at this scale presented its own challenges in managing large data sets, with multiple challenges and hardware advances discussed that facilitated this work. Over approximately 750 ns of aggregate simulation, we see multiple permeation events for these metabolites that were added at high concentration through the pores present within BMC shell tiles. When compared to independent permeability estimates for the same metabolites determined through replica exchange umbrella sampling simulations, the permeabilities varied by approximately 3 orders of magnitude. Regardless, the permeability coefficients for both G3P and DHAP are highly similar and very high, such that only very small concentration gradients can be maintained across the BMC shell between the cytosol and BMC interior. The large simulation systems also facilitated comparisons for molecular diffusivity in the crowded environment within the BMC shell. By our estimates, the viscosity within a packed BMC shell is at least 10-fold higher than it would be in neat solution and is the real driver for varying permeability estimates we obtained through simulation. These findings will be used as design inputs for future bioengineering efforts to make products from BMCs, highlighting how permeable BMC shells can be.

Diffusion↗