Search NASA⌕ Search

SEARCH · Search NASA

Results for “biochemical simulations”

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

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

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

Trapped-ion quantum simulation of electron transfer models with tunable dissipation

Electron transfer is at the heart of many fundamental physical, chemical, and biochemical processes essential for life. The exact simulation of these reactions is often hindered by the large number of degrees of freedom and by the essential role of quantum effects. Here, we experimentally simulate a paradigmatic model of molecular electron transfer using a multispecies trapped-ion crystal, where the donor-acceptor gap, the electronic and vibronic couplings, and the bath relaxation dynamics can all be controlled independently. By manipulating both the ground-state and optical qubits, we observe the real-time dynamics of the spin excitation, measuring the transfer rate in several regimes of adiabaticity and relaxation dynamics. Our results provide a testing ground for increasingly rich models of molecular excitation transfer processes that are relevant for molecular electronics and light-harvesting systems.

Science & Technology - Other Topics↗
DOE OSTI2022

Structure of maize BZR1-type β-amylase BAM8 provides new insights into its noncatalytic adaptation

Plant β-amylase (BAM) proteins play an essential role in growth, development, stress response, and hormone regulation. Despite their typical (β/α) 8 barrel structure as active catalysts in starch breakdown, catalytically inactive BAMs are implicated in diverse yet elusive functions in plants. Further, the noncatalytic BAM7/8 contain N-terminal BZR1 domains and were shown to be involved in the regulation of brassinosteroid signaling and possibly serve as sensors of yet an uncharacterized metabolic signal. While the structures of several catalytically active BAMs have been reported, structural characterization of the catalytically inactive BZR1-type BAMs remain unknown. Here, we determine the crystal structure of β-amylase domain of Zea mays BAM8/BES1/BZR1-5 and provide comprehensive insights into its noncatalytic adaptation. Using structural-guided comparison combined with biochemical analysis and molecular dynamics simulations, we revealed conformational changes in multiple distinct highly conserved regions resulting in rearrangement of the binding pocket. Altogether, this study adds a new layer of understanding to starch breakdown mechanism and elucidates the acquired adjustments of noncatalytic BZR1-type BAMs as putative regulatory domains and/or metabolic sensors in plants.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

Structural determinants for pH-dependent activation of a plant metacaspase

Arabidopsis thaliana metacaspase 9 (AtMC9) plays roles in clearing dead cells, forming xylem vessels, and regulating immunity and programmed cell death in plants. The protease's activation is controlled by pH levels, but the exact structural mechanism behind this has not been elucidated. In this work, high-resolution crystal structures for AtMC9 at active (pH 5.5 and pH 4.2) and inactive (pH 7.5) conditions are reported. The three structures are similar except for local conformations where their hydrogen bonding interactions with solvents are mediated through the protonation of specific titratable amino acid residues' side chains. By combining structural analysis, molecular dynamics simulations under constant pHs, and biochemical assays coupled with site-directed mutagenesis, we show that the regulation of AtMC9 activation involves multiple titratable glutamate and histidine residues across the three domains of p20, linker, and p10. Specifically, deprotonated Glu112, His193, and His208 can suppress AtMC9 proteolytic activity, while protonation of Glu255 and His307 at acidic pH may promote it. This study provides valuable insights into the pH-dependent activation of AtMC9 and could potentially lead to improving crops with enhanced immunity and controlled cell death, ultimately increasing agricultural productivity.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Ratchet, swivel, tilt and roll: a complete description of subunit rotation in the ribosome

Protein synthesis by the ribosome requires large-scale rearrangements of the ‘small’ subunit (SSU; ~1 MDa), including inter- and intra-subunit rotational motions. However, with nearly 2000 structures of ribosomes and ribosomal subunits now publicly available, it is exceedingly difficult to design experiments based on analysis of all known rotation states. To overcome this, we developed an approach where the orientation of each SSU head and body is described in terms of three angular coordinates (rotation, tilt and tilt direction) and a single translation. By considering the entire RCSB PDB database, we describe 1208 fully-assembled ribosome complexes and 334 isolated small subunits, which span >50 species. This reveals aspects of subunit rearrangements that are universal, and others that are organism/domain-specific. For example, we show that tilt-like rearrangements of the SSU body (i.e. ‘rolling’) are pervasive in both prokaryotic and eukaryotic (cytosolic and mitochondrial) ribosomes. As another example, domain orientations associated with frameshifting in bacteria are similar to those found in eukaryotic ribosomes. Together, this study establishes a common foundation with which structural, simulation, single-molecule and biochemical efforts can more precisely interrogate the dynamics of this prototypical molecular machine.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

GTP binding by Arabidopsis extra-large G protein 2 is not essential for its functions

The extra-large guanosine-5'-triphosphate (GTP)-binding protein 2, XLG2, is an unconventional Gα subunit of the Arabidopsis (Arabidopsis thaliana) heterotrimeric GTP-binding protein complex with a major role in plant defense. In vitro biochemical analyses and molecular dynamic simulations show that affinity of XLG2 for GTP is two orders of magnitude lower than that of the conventional Gα, AtGPA1. Here we tested the physiological relevance of GTP binding by XLG2. We generated an XLG2(T476N) variant with abolished GTP binding, as confirmed by in vitro GTPγS binding assay. Yeast three-hybrid, bimolecular fluorescence complementation, and split firefly-luciferase complementation assays revealed that the nucleotide-depleted XLG2(T476N) retained wild-type XLG2-like interactions with the Gβγ dimer and defense-related receptor-like kinases. Both wild-type and nucleotide-depleted XLG2(T476N) restored the defense responses against Fusarium oxysporum and Pseudomonas syringae compromised in the xlg2 xlg3 double mutant. Additionally, XLG2(T476N) was fully functional restoring stomatal density, root growth, and sensitivity to NaCl, but failed to complement impaired germination and vernalization-induced flowering. We conclude that XLG2 is able to function in a GTP-independent manner and discuss its possible mechanisms of action.

Plant Sciences↗
DOE OSTI2020

Structural adaptation of oxygen tolerance in 4-hydroxybutyrl-CoA dehydratase, a key enzyme of archaeal carbon fixation

Abstract Autotrophic microorganisms that convert inorganic carbon into organic matter were key players in the evolution of life on Earth. As the early atmosphere became oxygenated, microorganisms needed to develop mechanisms for oxygen protection, especially those relying on enzymes containing oxygen-sensitive metal clusters (e.g., Fe-S). Here we investigated how 4-hydroxybutyryl-CoA dehydratase (4HBD) - the key enzyme of the 3-hydroxypropionate/4-hydroxybutyrate (HP/HB) cycle for CO 2 -fixation - adapted as conditions shifted from anoxic to oxic. 4HBD is found in both anaerobic bacteria and aerobic ammonia-oxidizing archaea (AOA). The oxygen-sensitive bacterial 4HBD and oxygen-tolerant archaeal 4HBD share 59 % amino acid identity. To examine the structural basis of oxygen tolerance in archaeal 4HBD, we determined the atomic resolution structure of the enzyme. Two tunnels providing access to the canonical [4Fe-4S] cluster in oxygen-sensitive bacterial 4HBD were closed with four conserved mutations found in all aerobic AOA and other archaea. Further biochemical experiments and molecular dynamics simulations support our findings that restricting access to the active site is the key to oxygen tolerance, explaining how active site evolution drove a major evolutionary transition. Significance statement Autotrophy (primary production) was the first life strategy on Earth. Before photosynthesis (using solar energy to fix carbon dioxide), life relied on chemical reactions for energy. These chemosynthetic reactions are present in all domains of life, including archaea possessing the most energy-efficient carbon fixation pathway - the 3-hydroxypropionate/4-hydroxybutyrate cycle. This efficiency results from enzyme modifications, including enhanced enzyme stability and catalysis of multiple reactions. We reveal the first structure of aerobic 4-hydroxybutyryl-CoA dehydratase (4HBD) from ammonia-oxidizing archaea. These archaea are among the most abundant organisms on the planet, and their 4HBD active site evolved oxygen tolerance to support aerobic metabolism. This modification can provide further insight into enzyme evolution on early earth, as photosynthesis developed and began oxygenating the atmosphere.

DeMirci, Hasan (ORCID:0000000291355397)↗
DOE OSTI2020

Molecular Mechanism of Polysaccharide Acetylation by the Arabidopsis Xylan O-acetyltransferase XOAT1

Xylans are a major component of plant cell walls. O-Acetyl moieties are the dominant backbone substituents of glucuronoxylan in dicots and play a major role in the polymer-polymer interactions that are crucial for wall architecture and normal plant development. Here, we describe the biochemical, structural, and mechanistic characterization of Arabidopsis (Arabidopsis thaliana) xylan O-acetyltransferase 1 (XOAT1), a member of the plant-specific Trichome Birefringence Like (TBL) family. Detailed characterization of XOAT1-catalyzed reactions by real-time NMR confirms that it exclusively catalyzes the 2-O-acetylation of xylan, followed by nonenzymatic acetyl migration to the O-3 position, resulting in products that are monoacetylated at both O-2 and O-3 positions. In addition, we report the crystal structure of the catalytic domain of XOAT1, which adopts a unique conformation that bears some similarities to the a/ß/a topology of members of the GDSL-like lipase/acylhydrolase family. Finally, we use a combination of biochemical analyses, mutagenesis, and molecular simulations to show that XOAT1 catalyzes xylan acetylation through formation of an acyl-enzyme intermediate, Ac–Ser-216, by a double displacement bi-bi mechanism involving a Ser-His-Asp catalytic triad and unconventionally uses an Arg residue in the formation of an oxyanion hole.

acetylation↗
DOE OSTI2022

Lytic polysaccharide monooxygenase increases cellobiohydrolases activity by promoting decrystallization of cellulose surface

Efficient depolymerization of crystalline cellulose requires cooperation between multiple cellulolytic enzymes. Through biochemical approaches, molecular dynamics (MD) simulation, and single-molecule observations using high-speed atomic force microscopy (HS-AFM), we quantify and track synergistic activity for cellobiohydrolases (CBHs) with a lytic polysaccharide monooxygenase (LPMO) from Phanerochaete chrysosporium. Increasing concentrations of LPMO (AA9D) increased the activity of a glycoside hydrolase family 6 CBH, Cel6A, whereas the activity of a family 7 CBH (Cel7D) was enhanced only at lower concentrations of AA9D. MD simulation suggests that the result of AA9D action to produce chain breaks in crystalline cellulose can oxidatively disturb the crystalline surface by disrupting hydrogen bonds. HS-AFM observations showed that AA9D increased the number of Cel7D molecules moving on the substrate surface and increased the processivity of Cel7D, thereby increasing the depolymerization performance, suggesting that AA9D not only generates chain ends but also amorphizes the crystalline surface, thereby increasing the activity of CBHs.

09 BIOMASS FUELS↗
DOE OSTI2020

Thermosteam: BioSTEAM's Premier Thermodynamic Engine

Thermosteam is a thermodynamic engine capable of solving mass and energy balances, estimating mixture properties, solving thermodynamic phase equilibria, and modeling stoichiometric reactions. All chemical data in Thermosteam is imported from the chemicals library, an open-source compilation of data and functions for the estimation of pure component chemical and mixture properties. Thermosteam’s fast and flexible platform has enabled the evaluation of conceptual and emerging biochemical production processes. The Biorefinery Simulation and Techno-Economic Analysis Modules (BioSTEAM) — capable of modeling reactors, distillation columns, heat exchangers, and other unit operations — has adopted Thermosteam as its premier thermodynamic engine. Published biorefinery designs modeled in BioSTEAM implement thermodynamic property packages created with Thermosteam, including a cornstover biorefinery for the production of cellulosic ethanol, a lipid-cane biorefinery for the co-production of ethanol and biodiesel, and a wheatstraw biorefinery for the production of cellulosic ethanol.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

The Roles of Tidal Marshes in the Estuarine Biochemical Processes: A Numerical Modeling Study

Observations suggest that the existence of tidal marsh can alter the oxygen and nutrient dynamics in adjacent water bodies, but assessing the impacts of large tidal marshes on an estuary is challenging. In this study, we use a modeling approach to investigate the roles of tidal marshes on the estuarine biochemical processes. The marsh model, which simulates the ecological functions of marshes at seasonal and annual time-scales, is embedded inside an unstructured-grid three-dimensional hydrodynamic and eutrophication model (SCHISM-ICM). This modeling system simulates the growth and metabolism of the tidal marshes and links biological processes to nutrient dynamics in the water column and sediment. This model dynamically simulates nutrient recycling and physical transport of the materials between marshes and open water through wetting-drying processes. This coupled model system is validated and successfully applied to the York River Estuary. Model results suggest that tidal marshes influence the local diurnal dissolved oxygen (DO) cycle by exporting dissolved organic carbon and high sediment oxygen demand in the marsh system through the tidal exchange. The high deposition rates of organics and diurnal DO cycle enhance the sediment release of phosphorus. On the other hand, marshes tend to decrease dissolved inorganic nitrogen in the water column by settling particulate nutrients and enhancing the denitrification process. The study demonstrates that tidal marshes exert substantial impacts on the estuarine biochemical processes. The developed tidal marsh model enhances eutrophication modeling and advances the understanding of the feedback effects between marsh biogeochemistry and estuarine eutrophication processes on a systemic scale.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Machine Learning‐Augmented Molecular Dynamics Simulations (MD) Reveal Insights Into the Disconnect Between Affinity and Activation of ZTP Riboswitch Ligands

Abstract The challenge of targeting RNA with small molecules necessitates a better understanding of RNA–ligand interaction mechanisms. However, the dynamic nature of nucleic acids, their ligand‐induced stabilization, and how conformational changes influence gene expression pose significant difficulties for experimental investigation. This work employs a combination of computational and experimental methods to address these challenges. By integrating structure‐informed design, crystallography, and machine learning‐augmented all‐atom molecular dynamics simulations (MD), we synthesized, biophysically and biochemically characterized, and studied the dissociation of a library of small molecule activators of the 5‐aminoimidazole–4–carboxamide ribonucleotide triphosphate (ZTP) riboswitch, a ligand‐binding RNA motif that regulates bacterial gene expression. We uncovered key interaction mechanisms, revealing valuable insights into the role of ligand binding kinetics on riboswitch activation. Further, we established that ligand on‐rates determine activation potency as opposed to binding affinity and elucidated RNA structural differences, which provide mechanistic insights into the interplay of RNA structure on riboswitch activation.

Chemistry↗
DOE OSTI2026

Unraveling design principles of protein landscapes in photosynthetic membranes in plant chloroplasts

The supramolecular organization of proteins within photosynthetic membranes is crucial for energy conversion in plants. Here, we introduce an analytical and computational pipeline that integrates high-resolution cryo–scanning electron microscopy, biochemical quantification, advanced Monte Carlo computer simulations, and statistical methods to elucidate the elusive protein landscapes of grana membranes in intact Arabidopsis leaves. Our integrated analysis challenges the prevailing view that particles on the exoplasmic fracture faces in freeze-fracture samples represent photosystem II exclusively. Instead, these particles also include cytochrome b 6 f complexes. Furthermore, our steric clash analysis demonstrates that stacked membranes contain a mixture of larger PSII supercomplexes (C 2 S 2 M 2 and C 2 S 2 ) in addition to a smaller complex (C 2 ). This suggests that in vivo PSII supercomplexes exist in an equilibrium distribution of differing sizes. Furthermore, we discovered that, although size exclusion effects govern the global protein arrangement, local packing exhibits orientational order indicative of lateral attractive protein-protein interactions.

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

Insights into the dual cleavage activity of the GH16 laminarinase enzyme class on β-1,3 and β-1,4 glycosidic bonds

Glycoside hydrolases (GH) are involved in the degradation of a wide diversity of carbohydrates and present several biotechnological applications. Many GH families are composed of enzymes with a single well-defined specificity. In contrast, enzymes from the GH16 family can act on a range of different polysaccharides, including β-glucans and galactans. SCLam, a GH16 member derived from a soil metagenome, an endo-β-1,3(4)-glucanase (EC 3.2.1.6), can cleave both β-1,3 and β-1,4 glycosidic bonds in glucans, such as laminarin, barley β-glucan, and cello-oligosaccharides. A similar cleavage pattern was previously reported for other GH16 family members. However, the molecular mechanisms for this dual cleavage activity on (1,3)- and (1,4)-β-D-glycosidic bonds by laminarinases have not been elucidated. In this sense, we determined the X-ray structure of a presumably inactive form of SCLam co-crystallized with different oligosaccharides. The solved structures revealed general bound products that are formed due to residual activities of hydrolysis and transglycosylation. Biochemical and biophysical analyses and molecular dynamics simulations help to rationalize differences in activity towards different substrates. Our results depicted a bulky aromatic residue near the catalytic site critical to select the preferable configuration of glycosidic bonds in the binding cleft. Altogether, these data contribute to understanding the structural basis of recognition and hydrolysis of β-1,3 and β-1,4 glycosidic linkages of the laminarinase enzyme class, which is valuable for future studies on the GH16 family members and applications related to biomass conversion into feedstocks and bioproducts.

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

Machine learning for surrogate process models of bioproduction pathways

Technoeconomic analysis and life-cycle assessment are critical to guiding and prioritizing bench-scale experiments and to evaluating economic and environmental performance of biofuel or biochemical production processes at scale. Traditionally, commercial process simulation tools have been used to develop detailed models for these purposes. However, developing and running such models can be costly and computationally intensive, which limits the degree to which they can be shared and reproduced in the broader research community. This study evaluates the potential of an automated machine learning approach to develop surrogate models based on conventional process simulation models. The analysis focuses on several high-value biofuels and bioproducts for which pathways of production from biomass feedstocks have been well-established. The results demonstrate that surrogate models can be an accurate and effective tool for approximating the cost, mass and energy balance outputs of more complex process simulations at a fraction of the computational expense.

09 BIOMASS FUELS↗
DOE OSTI2020

Polymerization and editing modes of a high-fidelity DNA polymerase are linked by a well-defined path

Proofreading by replicative DNA polymerases is a fundamental mechanism ensuring DNA replication fidelity. In proofreading, mis-incorporated nucleotides are excised through the 3'-5' exonuclease activity of the DNA polymerase holoenzyme. The exonuclease site is distal from the polymerization site, imposing stringent structural and kinetic requirements for efficient primer strand transfer. Yet, the molecular mechanism of this transfer is not known. Here we employ molecular simulations using recent cryo-EM structures and biochemical analyses to delineate an optimal free energy path connecting the polymerization and exonuclease states of E. coli replicative DNA polymerase Pol III. We identify structures for all intermediates, in which the transitioning primer strand is stabilized by conserved Pol III residues along the fingers, thumb and exonuclease domains. We demonstrate switching kinetics on a tens of milliseconds timescale and unveil a complete pol-to-exo switching mechanism, validated by targeted mutational experiments.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Mechanism of tethered agonist-mediated signaling by polycystin-1

Mutations of polycystin-1 (PC1) are the major cause (85% of cases) of autosomal dominant polycystic kidney disease (ADPKD), which is the fourth leading cause of kidney failure. PC1 is thought to function as an atypical G protein-coupled receptor, yet the mechanism by which PC1 regulates G-protein signaling remains poorly understood. A significant portion of ADPKD mutations of PC1 encode a protein with defects in maturation or reduced function that may be amenable to functional rescue. In this work, we have combined complementary biochemical and cellular assay experiments and accelerated molecular simulations, which revealed an allosteric transduction pathway in activation of the PC1 C-terminal fragment. Our findings will facilitate future rational drug design efforts targeting the PC1 signaling function. Polycystin-1 (PC1) is an important unusual G protein-coupled receptor (GPCR) with 11 transmembrane domains, and its mutations account for 85% of cases of autosomal dominant polycystic kidney disease (ADPKD). PC1 shares multiple characteristics with Adhesion GPCRs. These include a GPCR proteolysis site that autocatalytically divides these proteins into extracellular, N-terminal, and membrane-embedded, C-terminal fragments (CTF), and a tethered agonist (TA) within the N-terminal stalk of the CTF that is suggested to activate signaling. However, the mechanism by which a TA can activate PC1 is not known. Here, we have combined functional cellular signaling experiments of PC1 CTF expression constructs encoding wild type, stalkless, and three different ADPKD stalk variants with all-atom Gaussian accelerated molecular dynamics (GaMD) simulations to investigate TA-mediated signaling activation. Correlations of residue motions and free-energy profiles calculated from the GaMD simulations correlated with the differential signaling abilities of wild type and stalk variants of PC1 CTF. They suggested an allosteric mechanism involving residue interactions connecting the stalk, Tetragonal Opening for Polycystins (TOP) domain, and putative pore loop in TA-mediated activation of PC1 CTF. Key interacting residues such as N3074–S3585 and R3848–E4078 predicted from the GaMD simulations were validated by mutagenesis experiments. Together, these complementary analyses have provided insights into a TA-mediated activation mechanism of PC1 CTF signaling, which will be important for future rational drug design targeting PC1.

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

Tight Coupling of Surface and In-Plant Biochemistry and Convection Governs Key Fine Particulate Components over the Amazon Rainforest

We report that throughout the troposphere over the Amazon rainforest, isoprene epoxydiol-secondary organic aerosols (IEPOX-SOA) are key components of submicrometer particle mass. IEPOX-SOA is thought to form by multiphase chemical pathways in the atmosphere, but using detailed model simulations, we show that these pathways are strongly inhibited by the solid particle phase state prevailing in the upper troposphere. We find that 2-methyltetrol gases formed by in-plant biochemical oxidation pathways, their partitioning to fine particles, and dissolution in cloud droplets can explain over 90% of IEPOX-SOA mass concentrations in the upper troposphere. In agreement with previous PM10 chiral measurements, our simulations indicate that plant biochemical processes are the dominant source of particulate 2-methyltetrols in the lower troposphere as well (>70%). Our results imply that in-plant biochemistry plays a central but previously unidentified role in aerosol processes and atmosphere-biosphere-climate interactions over the rainforest.

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

Whole cell simulator

A whole cell model (WCM) is a comprehensive multi-scale computational model representing all the known biochemical processes in a cell. It relies on a variety of intracellular pathway models and omics data. This computational framework enables seamless integration of diverse simulation methods used in WCM such as stochastic simulation algorithm (SSA), ordinary differential equations (ODE), flux balance analysis (FBA), and logic-based approaches. These methods run simultaneously, not only for whole pathways but also for subsets of reactions. Furthermore, we aim to enable dynamic switching between the methods when beneficial. Is

GEORGAKOUDIS, GIORGIS↗