Search NASA⌕ Search

SEARCH · Search NASA

Results for “Biologically plausible”

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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Electromechanical memcapacitor model offering biologically plausible spiking.

In this article, we introduce a new nanoscale electromechanical device - a leaky memcapacitor - and show that it may be useful for the hardware implementation of spiking neurons. The leaky memcapacitor is a movableplate capacitor that becomes quite conductive when the plates come close to each other. The equivalent circuit of the leaky memcapacitor involves a memcapacitive and memristive system connected in parallel. In the leaky memcapacitor, resistance and capacitance depend on the same internal state variable, which is the displacement of the movable plate. We have performed a comprehensive analysis showing that several types of spiking observed in biological neurons can be implemented with the leaky memcapacitor. Significant attention is paid to the dynamic properties of the model. As in leaky memcapacitors the capacitive, leaking resistive, and reset functionalities are implemented naturally within the same device structure, their use will simplify the creation of spiking neural networks.

Zhang, Zixi↗
DOE OSTI2025

Backpropagation-based learning with local derivative approximation and memory replay in biologically plausible neural systems

When learning, the brain modifies individual synaptic connections to reach a desired behavior. Animal and human brains have been shown to be incredibly capable of learning complex and varied functions across a wide variety of tasks. In recent years, artificial neural networks, inspired by human and animal brains, have shown great capabilities in learning a wide variety of difficult tasks. However, artificial neural networks primarily teach themselves through the use of backpropagation, a learning method which has no clear analogue within the brain. Additionally, Artificial Neural Networks primarily use continuous activation functions, which differ significantly from the spiking neuronal behavior present in the brain. In this paper, we discuss and demonstrate a biologically plausible learning method that approximates backpropagation through two techniques on Spiking Neural Networks. First, we show that the local temporal derivatives that are necessary for backpropagation can be approximately recovered through reconstruction using spike timings. Second, we show that through learning during a sleep phase, inspired by neuroscience research into memory replay, the localized parallel feedback path can learn to approximate the derivative through the forward path weight matrix, thus solving the weight transport problem. Lastly, we demonstrate that the combination of these two methods can approach or exceed the accuracy of backpropagation-based methods for a variety of neuromorphic vision tasks while maintaining biological plausibility.

42 ENGINEERING↗
DOE OSTI2026

A Unidirectional Two-Compartment Neuron Circuit with On-chip STDP learning

Most neuromorphic chips implement the single-compartment point neuron model where synapse circuits connect directly to a leaky integrate and fire (LIF) soma circuit. However, when using a biologically plausible soma circuit (e.g., Hodgkin-Huxley neuron model), an interface circuitry, such as a current conveyor circuit, is needed to transmit synaptic current to the soma circuit. This is especially true for ultra-low power neuron circuits, where membrane capacitance is on the order of 20 fF. This need for an interface circuit arises because the parasitic capacitance and leakage current caused by fabrication mismatch and second-order effects of the output transistors in the synapse circuits can disturb the spiking dynamics of the soma circuit if connected without an interface. Using an interface circuit to isolate the soma’s membrane capacitor from synapses resolves this issue. We propose to use a unidirectional resistor (a transconductance circuit) to connect the synapse and soma circuits instead of conventional current conveyor circuits. Using a biologically plausible spike pattern detection model, we show that the on-chip spike-timing-dependent plasticity (STDP) learning performance of the proposed unidirectional two-compartment neuron circuit is similar to a single-compartment circuit (with a current conveyor as an interface) and additionally, it is more power-efficient and biologically plausible. The chip is fabricated in the Taiwan Semiconductor Manufacturing Company (TSMC) 250 nm technology node and comprises a single neuron circuit.

Gautam, Ashish [ORNL]↗
DOE OSTI2025

Sulfur Species in Zinc-Rich Condylar Zones of a Rat Temporomandibular Joint

We performed synchrotron-based micro-X-ray absorption near-edge spectroscopy (µ-XANES) coupled with micro-X-ray fluorescence (µ-XRF) for the identification of elements that included biometal zinc (Zn) and nonmetal sulfur (S) (and its species) in the condylar zones of a rat temporomandibular joint (TMJ). Zone-specific spatial localization of biometal Zn and nonmetal S from a materials viewpoint when correlated with hypoxia inducible factor-1α (HIF-1α) (a surrogate for tissue oxygenation) can provide insights into Zn-specific redox pathways at the vulnerable subchondral interface. Histologic localization of Zn, HIF-1α, and sulfur-rich proteoglycans (PGs) were mapped using an optical microscope. The µ-XRF maps coupled with site-specific micro-X-ray diffraction (µ-XRD) patterns were used to underline Zn-incorporated biological apatite in the subchondral bone and the bone of a rat TMJ condyle. Results demonstrated an association between Zn, PG, and HIF-1α histologic maps with µ-XRF, µ-XANES, and µ-XRD data and provided insights into plausible biological S species in Zn-enriched zones of a rat TMJ condyle. Spatially localized Zn and S underscore their roles in cell and tissue functions in a zone-specific manner. Elemental Zn with organic and inorganic S species at the cartilage-bone interface and the biomineral phase of Zn-enriched biological apatite from subchondral bone to condylar bone were ascertained using µ-XRF-XANES and µ-XRF-XRD. The coupled µ-XRF-XANES in situ complemented with µ-XRF-XRD in situ and immunohistochemistry provided valuable biological insights into zone-specific biological pathways in rat TMJ condyles. Based on these data, we present a workflow to reliably map and correlate S species within Zn-enriched regions of cartilage, bone, and their interface. We suggest the use of this correlative and complementary microspectroscopic spatial information for zone-specific localization of biometal Zn and nonmetal S to gain insights into plausible microanatomy-specific oxidative stress in the TMJ.

apatites↗
DOE OSTI2024

Chaotic neural dynamics facilitate probabilistic computations through sampling

Cortical neurons exhibit highly variable responses over trials and time. Theoretical works posit that this variability arises potentially from chaotic network dynamics of recurrently connected neurons. Here, we demonstrate that chaotic neural dynamics, formed through synaptic learning, allow networks to perform sensory cue integration in a sampling-based implementation. We show that the emergent chaotic dynamics provide neural substrates for generating samples not only of a static variable but also of a dynamical trajectory, where generic recurrent networks acquire these abilities with a biologically plausible learning rule through trial and error. Furthermore, the networks generalize their experience in the stimulus-evoked samples to the inference without partial or all sensory information, which suggests a computational role of spontaneous activity as a representation of the priors as well as a tractable biological computation for marginal distributions. These findings suggest that chaotic neural dynamics may serve for the brain function as a Bayesian generative model.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2025

CRCNS22 Learning Rules in the Hippocampus and their Mapping to Neuromorphic Systems (Final Technical Report)

Large scale biologically-realistic computational models are key to investigating the interplay between structure and function in nervous systems, thus paving the way to new clinical methods and neuro-inspired computing solutions. This project focuses on the hippocampus, in particular the CA3-CA1 regions, due to their role in associative learning and memory, pattern separation and completion, and spatial navigation. Investigations into the neuronal organization and learning rule(s) of this circuit can shed light into how declarative memories are formed, stored, recalled and forgotten and inform computational, experimental and clinical neuroscience work. Our project aims at developing a novel data-driven methodology supported by a broad heterogeneous base of neuroscience experimental knowledge and inspired from advances in computer science and engineering. Specifically, this work will benchmark existing and new learning rules within a full-scale spiking neural network simulation of the CA3-CA1 region. The model will be based on an open-source repository, called the Hippocampome, which contains neuronal morphologies, firing patterns, synapse probabilities, and most other required parameters for all known neuron types in the rodent hippocampal formation. The model will be first trained in a supervised fashion for associative memory tasks using backpropagation through time traditionally used in computer science, enhanced with a new technique called the surrogate gradient method. This optimization method will be used to obtain a global loss minimization, but it is not biologically inspired as it assumes the use of data not locally available to the synapses. However, we propose its use as a benchmarking tool, to compare the training performance of local biologically plausible and hardware-mappable learning rules at scale. New rules or combinations will be proposed and tested as needed, based on the obtained results. Progress in this area will also drive the development of novel hardware-mappable algorithms for continual lifelong learning and categorization of new events from few presented examples. This project goes beyond the existing state-of-the-art by looking at large scale realistic neuronal circuits as networks trainable via global optimization methods such as surrogate gradient descent. The objective function of the brain that supports learning is largely unknown, but it is likely that it operates through local learning rules. Studying network trajectories around local minima as proposed in this work represents a useful strategy for understanding whether a network is training by using a specific (set of) learning rule(s). Starting from a completely untrained network is a challenging test since it is difficult to determine how the learning rule affects the trajectory of the network. This interdisciplinary project will help understand what rule governs learning in these regions or if multiple learning rules are involved. The work will develop a robust methodology to measure if the network is converging to the target solution, oscillating around it, or diverging away.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Emergence of catalytic function in prebiotic information-coding polymers

Life as we know it relies on the interplay between catalytic activity and information processing carried out by biological polymers. Here we present a plausible pathway by which a pool of prebiotic information-coding oligomers could acquire an early catalytic function, namely sequence-specific cleavage activity. Starting with a system capable of non-enzymatic templated replication, we demonstrate that even non-catalyzed spontaneous cleavage would promote proliferation by generating short fragments that act as primers. Furthermore, we show that catalytic cleavage function can naturally emerge and proliferate in this system. Specifically, a cooperative catalytic network with four subpopulations of oligomers is selected by the evolution in competition with chains lacking catalytic activity. The cooperative system emerges through the functional differentiation of oligomers into catalysts and their substrates. The model is inspired by the structure of the hammerhead RNA enzyme as well as other DNA- and RNA-based enzymes with cleavage activity that readily emerge through natural or artificial selection. We identify the conditions necessary for the emergence of the cooperative catalytic network. In particular, we show that it requires the catalytic rate enhancement over the spontaneous cleavage rate to be at least 10 2 –10 3 , a factor consistent with the existing experiments. The evolutionary pressure leads to a further increase in catalytic efficiency. The presented mechanism provides an escape route from a relatively simple pairwise replication of oligomers toward a more complex behavior involving catalytic function. This provides a bridge between the information-first origin of life scenarios and the paradigm of autocatalytic sets and hypercycles, albeit based on cleavage rather than synthesis of reactants.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Emergence of catalytic function in prebiotic information-coding polymers

Life as we know it relies on the interplay between catalytic activity and information processing carried out by biological polymers. Here we present a plausible pathway by which a pool of prebiotic information-coding oligomers could acquire an early catalytic function, namely sequence-specific cleavage activity. Starting with a system capable of non-enzymatic templated replication, we demonstrate that even non-catalyzed spontaneous cleavage would promote proliferation by generating short fragments that act as primers. Furthermore, we show that catalytic cleavage function can naturally emerge and proliferate in this system. Specifically, a cooperative catalytic network with four subpopulations of oligomers is selected by the evolution in competition with chains lacking catalytic activity. The cooperative system emerges through the functional differentiation of oligomers into catalysts and their substrates. The model is inspired by the structure of the hammerhead RNA enzyme as well as other DNA- and RNA-based enzymes with cleavage activity that readily emerge through natural or artificial selection. We identify the conditions necessary for the emergence of the cooperative catalytic network. In particular, we show that it requires the catalytic rate enhancement over the spontaneous cleavage rate to be at least 10 2 –10 3 , a factor consistent with the existing experiments. The evolutionary pressure leads to a further increase in catalytic efficiency. The presented mechanism provides an escape route from a relatively simple pairwise replication of oligomers toward a more complex behavior involving catalytic function. This provides a bridge between the information-first origin of life scenarios and the paradigm of autocatalytic sets and hypercycles, albeit based on cleavage rather than synthesis of reactants.

Tkachenko, Alexei V. (ORCID:000000031291243X)↗
DOE OSTI2024

Digital Multiphase Composites via Additive Manufacturing

Abstract Mechanical properties of traditional engineering materials are typically coupled to each other, presenting a challenge to practitioners with multi‐dimensional material property requirements. In this work, continuous, independent control over multiple mechanical properties is demonstrated in composite materials realized using additive manufacturing. For the first time, composites additively manufactured from rigid plastic, soft elastomer, and liquid constituents are experimentally characterized, demonstrating materials which span four orders of magnitude in modulus and two orders of magnitude in toughness. By forming analytical mappings between relative concentrations of constituents at the microscale and resulting macroscale material properties, inverse material design is enabled; the method is showcased by printing artifacts with prescribed toughness and elasticity distributions. The properties of these composites are placed in the context of biological tissues, showing they have promise as mechanically plausible tissue mimics.

Smith, Lawrence T.↗
DOE OSTI2025

Gut microbiome changes with micronutrient supplementation in children with attention–deficit/hyperactivity disorder: the MADDY study

Micronutrients have demonstrated promise in managing inattention and emotional dysregulation in children with attention-deficit/hyperactivity disorder (ADHD). The biological mechanism by which micronutrients improve these symptoms remains unclear. One plausible pathway is through the gut-brain axis, the bi-directional communication network that links the gastrointestinal tract with the brain. This study examines changes in gut microbiome composition and diversity after micronutrients supplementation in children with ADHD (N=44) and sheds light on potential mechanisms responsible for the response to micronutrients as measured by clinician-rated global impression. Participants from this investigation represent a sub-group of the Micronutrients for ADHD in Youth (MADDY) study, a double blind randomized controlled study in which participants received either micronutrients or a placebo for 8 weeks, followed by an 8-week open label extension with micronutrients for all participants. Stool samples collected at baseline, week 8, and week 16 were analyzed using 16S rRNA amplicon sequencing targeting the V4 hypervariable region. Pairwise compositional analyses served as the primary means for investigating changes in gut microbiome composition between micronutrients versus placebo groups and responders versus non-responders. A significant change in microbial evenness, as measured by alpha diversity, was observed following micronutrients, and the phylum Actinobacteriota decreased in the micronutrients group compared to placebo. Additionally, two bacterial families: Rikenellaceae and Oscillospiraceae, exhibited a significant increase in change of gut microbiome composition following micronutrients between responders and non-responders. These findings suggest that micronutrients modulated the composition of the gut microbiome and point towards specific bacterial changes associated with response to micronutrients.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2025

Fundamental Research Aimed at Diverting Excess Reducing Power in Photosynthesis to Orthogonal Metabolic Pathways

Photosystems are incredible biological machines that use sunlight to drive the conversion of carbon dioxide to sugar. The amount of sunlight available for photosynthesis sometimes exceeds the amount of energy plants can use. This excess energy has to be safely dissipated through non-productive biological processes. The ultimate goal of this project is to understand whether we can utilize that otherwise unused excess energy. In our previous work, we showed that, in principle, it is possible to attach a catalyst to photosystem I and generate H2 using light. Our current strategy is to genetically fuse parts of the photosystem I complex with a recently discovered oxygen-tolerant [FeFe] hydrogenase. Our rationale is that such chimeric proteins may potentially result in the natural incorporation of the photosystem I-hydrogenase link using the inherent genetic machinery of the cell. In this project, we aim to verify that light-driven hydrogen production in this construction is possible. Throughout the project, we designed nanoconstructs that showcase the plausibility of this technology, at least in vitro. We take advantage of these constructs to investigate details of the coupling between photosystem I and a H2-producing enzyme called [FeFe] hydrogenase. This part of the project reveals details of the electron transfer between photosystem I and the attached hydrogenase, providing information that can lead to new strategies for improved biological photocatalysis. We also researched efficient and robust tethering of the [FeFe] hydrogenase to photosystem I in cyanobacteria. This work will highlight successful design strategies to guide the future development of photosynthetic biohybrids. Uncovering the principles governing the utilization of otherwise unusable energy significantly further our understanding of cyanobacterial photosynthesis. The work proposed establishes the feasibility of diverting excess energy under high light conditions to orthogonal enzymatic pathways and set design rules for efficient utilization of such a strategy for scientific and industrial applications in biosensing, renewable energy, and high-value chemicals production. The work addresses the DOE-BES Photosynthetic Systems program goal to develop a multidimensional understanding of photosystems that would provide specific metrics that instruct strategies for improving biological photosynthesis and for guiding the future development of bioreactors and biomimetic energy systems.

Photosynthetic systems, hydrogenase, cyanobacteria↗
DOE OSTI2024

Major impacts of widespread structural variation on sorghum

Genetic diversity is critical to crop breeding and improvement, and dissection of the genomic variation underlying agronomic traits can both assist breeding and give insight into basic biological mechanisms. Although recent genome analyses in plants reveal many structural variants (SVs), most current studies of crop genetic variation are dominated by single-nucleotide polymorphisms (SNPs). The extent of the impact of SVs on global trait variation, as well as their utility in genome-wide selection, is not yet understood. In this study, we built an SV data set based on whole-genome resequencing of diverse sorghum lines (n = 363), validated the correlation of photoperiod sensitivity and variety type, and identified SV hotspots underlying the divergent evolution of cellulosic and sweet sorghum. In addition, we showed the complementary contribution of SVs for heritability of traits related to sorghum adaptation. Importantly, inclusion of SV polymorphisms in association studies revealed genotype–phenotype associations not observed with SNPs alone. Three-way genome-wide association studies (GWAS) based on whole-genome SNP, SV, and integrated SNP + SV data sets showed substantial associations between SVs and sorghum traits. The addition of SVs to GWAS substantially increased heritability estimates for some traits, indicating their important contribution to functional allelic variation at the genome level. Our discovery of the widespread impacts of SVs on heritable gene expression variation could render a plausible mechanism for their disproportionate impact on phenotypic variation. This study expands our knowledge of SVs and emphasizes the extensive impacts of SVs on sorghum.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2026

Leveraging 13C-Labeling to Assign Molecular Formulas to Unknown Yeast Metabolites

Mass spectrometry analyses have identified tens of thousands of unknown small molecule-associated peaks in different biological specimens. Notably, even the simplest and best studied organisms like Escherichia coli and Saccharomyces cerevisiae yield thousands of unknown peaks. A key question is how many of these reflect actual novel endogenous metabolites. To explore this, Mahieu and Patti used complete 13 C -labeling in E. coli to credential peaks as biological. This reduced the number of unknowns by more than 90%. Here, we carry out similar uniform 13 C-labeling in the Baker’s yeast S. cerevisiae and two less-studied bioenergy-relevant yeasts Rhodotorula toruloides (lipid producer) and Issatchenkia orientalis (organic acid producer). Identification of unknown metabolite peaks and their molecular formulas is facilitated through software tailored for 13 C labeling data and resulting knowledge of carbon atom count. A classification model evaluates the plausibility of each candidate formula, with peaks lacking plausible candidate formulas unlikely to reflect metabolite molecular ions. This approach prioritizes about one hundred candidate abundant unknown metabolites with logical molecular formulas. Most of these are species-specific rather than conserved across yeasts, and more are found in the nonmodel yeasts than S. cerevisiae. Thus, 13 C-labeling data on unknown metabolites highlights the potential for discovering new metabolites and pathways in nonmodel yeasts.

Carbon↗
DOE OSTI2024

Structural and Interactional Analysis of the Flavonoid Pathway Proteins: Chalcone Synthase, Chalcone Isomerase and Chalcone Isomerase-like Protein

Chalcone synthase (CHS) and chalcone isomerase (CHI) catalyze the first two committed steps of the flavonoid pathway that plays a pivotal role in the growth and reproduction of land plants, including UV protection, pigmentation, symbiotic nitrogen fixation, and pathogen resistance. Based on the obtained X-ray crystal structures of CHS, CHI, and chalcone isomerase-like protein (CHIL) from the same monocotyledon, Panicum virgatum, along with the results of the steady-state kinetics, spectroscopic/thermodynamic analyses, intermolecular interactions, and their effect on each catalytic step are proposed. In addition, PvCHI’s unique activity for both naringenin chalcone and isoliquiritigenin was analyzed, and the observed hierarchical activity for those type-I and -II substrates was explained with the intrinsic characteristics of the enzyme and two substrates. The structure of PvCHS complexed with naringenin supports uncompetitive inhibition. PvCHS displays intrinsic catalytic promiscuity, evident from the formation of p-coumaroyltriacetic acid lactone (CTAL) in addition to naringenin chalcone. In the presence of PvCHIL, conversion of p-coumaroyl-CoA to naringenin through PvCHS and PvCHI displayed ~400-fold increased Vmax with reduced formation of CTAL by 70%. Supporting this model, molecular docking, ITC (Isothermal Titration Calorimetry), and FRET (Fluorescence Resonance Energy Transfer) indicated that both PvCHI and PvCHIL interact with PvCHS in a non-competitive manner, indicating the plausible allosteric effect of naringenin on CHS. Significantly, the presence of naringenin increased the affinity between PvCHS and PvCHIL, whereas naringenin chalcone decreased the affinity, indicating a plausible feedback mechanism to minimize spontaneous incorrect stereoisomers. These are the first findings from a three-body system from the same species, indicating the importance of the macromolecular assembly of CHS-CHI-CHIL in determining the amount and type of flavonoids produced in plant cells.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Dirigent isoflavene-forming PsPTS2: 3D structure, stereochemical, and kinetic characterization comparison with pterocarpan-forming PsPTS1 homolog in pea

Pea (legume plant defense) phytoalexin pterocarpans (–)-maackiain and (+)-pisatin have opposite C6a and C11a configurations, but how this occurs biosynthetically is unknown. Herein, pea dirigent-protein (DP) PsPTS2, which generates 7,2'-dihydroxy-4',5'-methylenedioxyisoflav-3-ene (DMDIF), was investigated as regards stereoselectivity towards four possible 7,2'-dihydroxy-4',5'-methylenedioxyisoflavan-4-ol (DMDI) stereoisomers. Each DMDI stereoisomer configuration was determined using NMR spectroscopy, electronic circular dichroism, and molecular orbital analyses. PsPTS2 efficiently converted cis-(3R,4R)-DMDI into DMDIF, this being 20-fold faster than with the trans-(3R,4S)-isomer. The 4R-configured substrate’s near ß-axial OH orientation significantly enhanced its leaving group abilities in creating intermediate A-ring mono-quinone methide (QM), whereas the 4S-isomer’s 4-OH group α-equatorial orientation was a poorer leaving group. Substrate docking simulations also indicated that the 4R-configured ß-axial OH was closest to the Asp catalytic center, whereas the 4S-isomer’s α-equatorial OH was further away. Neither cis-(3S,4S)- nor trans-(3S, 4R)-DMDIs were substrates, even with the former having the same C3/C4 stereochemistry as (+)-pisatin. PsPTS2, assayed with cis-(3R,4R)-7,2'-dihydroxy-4'-methoxyisoflavan-4-ol [cis-(3R,4R)-DMI] and its C3/C4 stereoisomers, gave 2',7-dihydroxy-4'-methoxyisoflav-3-ene (DMIF), with the same substrate stereoselectivity. DP homologs may exist in licorice (Glycyrrhiza pallidiflora) and the tree legume Bolusanthus speciosus, as DMIF is found in both. PsPTS1 preferentially utilized cis-(3R,4R)-DMDI to give (–)-maackiain, this being 2200-fold more efficient than with cis-(3R,4R)-DMI. PsPTS1 also slowly converted trans-(3S,4R)-DMDI into (+)-maackiain, again reflecting the better 4R configured OH leaving group. Pea PsPTS2 and PsPTS1 provisionally provide a facile DP-engendered means to enable differing C6a and C11a configurations in (+)-pisatin and (–)-maackiain. Both DP transformations are considered to occur via DP-engendered generation of an identical mono-QM bound intermediate, which with PsPTS2 either re-aromatizes to give DMDIF or with PsPTS1 undergoes intramolecular cyclization to afford (–)-maackiain. Substrate docking simulations using PsPTS2 (2.5 Å resolution) and PsPTS1 (1.5 Å resolution) plausibly indicate that cis-(3R,4R)-DMDI binds in the anti-configuration in the PsPTS2 active site to afford DMDIF, and the syn-configuration in PsPTS1 to give maackiain.

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

Combining Reactive Quantum-Mechanical Molecular-Dynamics Simulations with Mutagenesis, Crystallography, and Enzyme Kinetics to Reveal Plausible Steps of Isocyanide Hydratase Catalysis

A complete understanding of enzyme mechanisms requires atomistic details of chemical reactions. Quantum-based molecular dynamics simulations (QMD) are a potential source of this information, but trade-offs between accuracy and computational cost have limited their use. We previously developed extended Lagrangian Born–Oppenheimer molecular dynamics (XL-BOMD) methods that leverage a negligible compromise in accuracy to substantially decrease the cost of QMD simulations. Here, we develop a reactive QMD approach using the latest XL-BOMD formulation, which enables efficient simulations of highly reactive systems, and use it to investigate mechanisms of intermediate formation in isocyanide hydratase (ICH) catalysis. In QMD simulations, molecular analogs of ICH active site residues reacted with para-nitrophenyl isocyanide, forming a thioimidate. Analysis of simulated atomic configurational and charge dynamics revealed a pathway where protonation of the isocyanide carbon occurs prior to thioimidate formation and suggested a possible role of Asp17 as a proton donor in the early phase of ICH catalysis. To test whether the pathway seen using the reactive QMD approach might be relevant to ICH catalysis, we performed X-ray crystallography and pre-steady-state enzyme kinetics studies of wild-type and D17N mutant ICH. Both the structure and kinetics are sensitive to the D17N mutation in a manner that is consistent with the order of the reaction steps seen in the simulations. Mobile protons play essential roles in many enzymes, yet they are difficult to observe experimentally, making the ordering of proton-dependent steps ambiguous in many enzyme mechanisms. The ability to directly simulate model reactions for the design of experiments that provide information about enzyme mechanisms involving mobile protons demonstrates the significance of our reactive QMD approach and motivates further biological applications.

36 MATERIALS SCIENCE↗
DOE OSTI2025

The Role of Bedrock Circulation Depth and Porosity in Mountain Streamflow Response to Prolonged Drought

Quantitative understanding is lacking on how the depth of active groundwater circulation in bedrock affects mountain streamflow response to a multi-year drought. We use an integrated hydrological model to explore the sensitivity of a variety of streamflow metrics to bedrock circulation depth and porosity under a plausible extreme drought scenario lasting up to 5 years. Endmember depth versus hydraulic conductivity relationships and porosity values for fractured crystalline rock are simulated. With drought, a deeper circulation system with higher drainable porosity more effectively buffers minimum flow and significantly limits perennial stream loss in comparison to a shallow circulation system. Streamflow buffering is accomplished through extensive groundwater storage loss. However, deeper circulation systems experience prolonged recovery from drought in comparison to storage-limited shallow systems. Research highlights the importance of characterizing the deeper bedrock hydrogeology in mountainous watersheds to better understand and predict drought impacts on stream ecosystem health and water resource sustainability.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Structural insights into the unexpected agonism of tetracyclic antidepressants at serotonin receptors 5-HT 1e R and 5-HT 1F R

Serotonin [5-hydroxytryptamine (5-HT)] acts via 13 different receptors in humans. Of these receptor subtypes, all but 5-HT 1e R have confirmed roles in native tissue and are validated drug targets. Despite 5-HT 1e R’s therapeutic potential and plausible druggability, the mechanisms of its activation remain elusive. To illuminate 5-HT 1e R’s pharmacology in relation to the highly homologous 5-HT 1F R, we screened a library of aminergic receptor ligands at both receptors and observe 5-HT 1e R/5-HT 1F R agonism by multicyclic drugs described as pan-antagonists at 5-HT receptors. Potent agonism by tetracyclic antidepressants mianserin, setiptiline, and mirtazapine suggests a mechanism for their clinically observed antimigraine properties. Using cryo-EM and mutagenesis studies, we uncover and characterize unique agonist-like binding poses of mianserin and setiptiline at 5-HT 1e R distinct from similar drug scaffolds in inactive-state 5-HTR structures. Together with computational studies, our data suggest that these binding poses alongside receptor-specific allosteric coupling in 5-HT 1e R and 5-HT 1F R contribute to the agonist activity of these antidepressants.

59 BASIC BIOLOGICAL SCIENCES↗