Cellerator: a system for simulating biochemical reaction networks
Explore the source record for details and available documents.
signal transduction↗
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.
Biochemical Process Modeling and Simulation (BPMS)
The Biochemical Process Modeling and Simulation project aims to reduce the cost and time of research by applying theory, modeling, and simulation to the most relevant bottlenecks in the biochemical process. We use molecular modeling, quantum mechanics, metabolic modeling, fluid dynamics, and reaction-diffusion methods in close collaboration with pretreatment, hydrolysis, upgrading, and TEA. The project's outcomes are increased yields and efficiency of the biochemical process, added value to products, and reduced price of fuels by specifically targeting catalytic efficiency, reactor design, enzyme efficiency, and microbial design. We work closely with experimental projects to identify problems and iterate with experiments to find and refine solutions. By working with experimentalists, we decide on problems that can be solved with simulation that could otherwise not be solved or would take too long with experiment alone to reach BETO's targets. Over the years, we have produced solutions that have resulted in determining the most likely fatty-acid derivative for passive transport out of bacteria that upgrade biomass, and we have also designed enzyme mutations for enhanced lignin upgrading. Metabolic models have been developed to tune the activity of 2,3 butanediol production for the 2030 target. A computational method to deliver understanding of how complex omics data can be interpreted in the metabolic pathways of organisms used in the Agile Biofoundry. We have found methods to overcome specific barriers and continue to develop those methods. Our reactor studies have guided the design of both the microbes and reactors for aerobic and micro-aerobic production at all scales and have been instrumental in improving the accuracy of techno-economic analysis models. This project is essential in the process of selecting the final processes for 2030 SAF production targets. More specifically, recently, we have: 1) Predicted the strength of the basic structural interactions in commodity plastics to provide guidance for plastics upcycling strategies. 2) Developed computational tool to improve the characterization of lignin-derived compounds 3) Developed new methodologies to enable Machine Learning-based Directed Evolution for protein engineering. 4) Developed Machine Learning methods to predict protein promiscuity and mutations to further improve microbial and enzymatic driven processes and demonstrated the utility of ML approaches to engineering proteins from sparse experimental datasets. 5) Developed new methods to enable high-fidelity simulation of aerobic fermentation at industrial scale and resolving mismatch of time scales through subcycling/operator splitting 7) Identified the difficulty in preventing local high-oxygen conditions in industrial bubble columns, which leads to less-desirable acetoin production, suggesting future research directions in alternative reactor configurations (e.g loop reactors, shallow-channel reactors).
BIOMASS FUELS↗
BETO 2021 Peer Review - Biochemical Process Modeling and Simulation (BPMS)
The Biochemical Process Modeling and Simulation project aims to reduce the cost and time of research by applying theory, modeling, and simulation to the most relevant bottlenecks in the biochemical process. We use molecular modeling, quantum mechanics, metabolic modeling, fluid dynamics, and reaction-diffusion methods in close collaboration with pretreatment, hydrolysis, upgrading, and TEA. The project's outcomes are increased yields and efficiency of the biochemical process, added value to products, and reduced price of fuels by specifically targeting catalytic efficiency, reactor design, enzyme efficiency, and microbial design.
29 ENERGY PLANNING, POLICY, AND ECONOMY↗
A Green's Function Approach to Simulate DNA Damage by the Indirect Effect
The DNA damage is of fundamental importance in the understanding of the effects of ionizing radiation. DNA is damaged by the direct effect of radiation (e.g. direct ionization) and by indirect effect (e.g. damage by.OH radicals created by the radiolysis of water). Despite years of research, many questions on the DNA damage by ionizing radiation remains. In the recent years, the Green's functions of the diffusion equation (GFDE) have been used extensively in biochemistry [1], notably to simulate biochemical networks in time and space [2]. In our future work on DNA damage, we wish to use an approach based on the GFDE to refine existing models on the indirect effect of ionizing radiation on DNA. To do so, we will use the code RITRACKS [3] developed at the NASA Johnson Space Center to simulate the radiation track structure and calculate the position of radiolytic species after irradiation. We have also recently developed an efficient Monte‐Carlo sampling algorithm for the GFDE of reversible reactions with an intermediate state [4], which can be modified and adapted to simulate DNA damage by free radicals. To do so, we will use the known reaction rate constants between radicals (OH, eaq, H,...) and the DNA bases, sugars and phosphates and use the sampling algorithms to simulate the diffusion of free radicals and chemical reactions with DNA. These techniques should help the understanding of the contribution of the indirect effect in the formation of DNA damage and double‐strand breaks.
Plante, Ianik↗
A Radiation Chemistry Code Based on the Greens Functions of the Diffusion Equation
Ionizing radiation produces several radiolytic species such as.OH, e-aq, and H. when interacting with biological matter. Following their creation, radiolytic species diffuse and chemically react with biological molecules such as DNA. Despite years of research, many questions on the DNA damage by ionizing radiation remains, notably on the indirect effect, i.e. the damage resulting from the reactions of the radiolytic species with DNA. To simulate DNA damage by ionizing radiation, we are developing a step-by-step radiation chemistry code that is based on the Green's functions of the diffusion equation (GFDE), which is able to follow the trajectories of all particles and their reactions with time. In the recent years, simulations based on the GFDE have been used extensively in biochemistry, notably to simulate biochemical networks in time and space and are often used as the "gold standard" to validate diffusion-reaction theories. The exact GFDE for partially diffusion-controlled reactions is difficult to use because of its complex form. Therefore, the radial Green's function, which is much simpler, is often used. Hence, much effort has been devoted to the sampling of the radial Green's functions, for which we have developed a sampling algorithm This algorithm only yields the inter-particle distance vector length after a time step; the sampling of the deviation angle of the inter-particle vector is not taken into consideration. In this work, we show that the radial distribution is predicted by the exact radial Green's function. We also use a technique developed by Clifford et al. to generate the inter-particle vector deviation angles, knowing the inter-particle vector length before and after a time step. The results are compared with those predicted by the exact GFDE and by the analytical angular functions for free diffusion. This first step in the creation of the radiation chemistry code should help the understanding of the contribution of the indirect effect in the formation of DNA damage and double-strand breaks.
Plante, Ianik↗
BIOCHEMICAL ACTIVITIES OF TERRESTRIAL MICROORGANISMS IN SIMULATED PLANETARY ENVIRONMENTS
Biochemical activities of terrestrial micro. organisms in simulated planetary environments
PLANETARY ATMOSPHERE↗
Blood biochemical and cellular changes during decompression and simulated extravehicular activity
Blood biochemical and cellular parameters were measured in human subjects before and after exposure to a decompression schedule involving 6 h of oxygen prebreathing. The exposure was designed to simulate extravehicular activity for 6 h (subjects performed exercise while exposed to 29.6 kPa). There were no significant differences between blood samples from subjects who were susceptible (n = 11) versus those who were resistant (n = 27) to formation of venous gas emboli. Although several statistically significant (P less than 0.05) changes in blood parameters were observed following the exposure (increases in white blood cell count, prothrombin time, and total bilirubin, and decreases in triglycerides, very-low-density lipoprotein cholesterol, and blood urea nitrogen), the changes were small in magnitude and blood factor levels remained within normal clinical ranges. Thus, the decompression schedule used in this study is not likely to result in blood changes that would pose a threat to astronauts during extravehicular activity.
Decompression Sickness/blood/etiology/prevention &↗
Biochemical activities of terrestrial microorganisms in simulated planetary environments Final report
Effects of simulated anaerobic planetary environment on biochemical activities of terrestrial microorganisms
PLANETARY ENVIRONMENT↗
SBbadger: biochemical reaction networks with definable degree distributions
Abstract Motivation An essential step in developing computational tools for the inference, optimization and simulation of biochemical reaction networks is gauging tool performance against earlier efforts using an appropriate set of benchmarks. General strategies for the assembly of benchmark models include collection from the literature, creation via subnetwork extraction and de novo generation. However, with respect to biochemical reaction networks, these approaches and their associated tools are either poorly suited to generate models that reflect the wide range of properties found in natural biochemical networks or to do so in numbers that enable rigorous statistical analysis. Results In this work, we present SBbadger, a python-based software tool for the generation of synthetic biochemical reaction or metabolic networks with user-defined degree distributions, multiple available kinetic formalisms and a host of other definable properties. SBbadger thus enables the creation of benchmark model sets that reflect properties of biological systems and generate the kinetics and model structures typically targeted by computational analysis and inference software. Here, we detail the computational and algorithmic workflow of SBbadger, demonstrate its performance under various settings, provide sample outputs and compare it to currently available biochemical reaction network generation software. Availability and implementation SBbadger is implemented in Python and is freely available at https://github.com/sys-bio/SBbadger and via PyPI at https://pypi.org/project/SBbadger/. Documentation can be found at https://SBbadger.readthedocs.io. Supplementary information Supplementary data are available at Bioinformatics online.
59 BASIC BIOLOGICAL SCIENCES↗
Cellerator: extending a computer algebra system to include biochemical arrows for signal transduction simulations
Cellerator describes single and multi-cellular signal transduction networks (STN) with a compact, optionally palette-driven, arrow-based notation to represent biochemical reactions and transcriptional activation. Multi-compartment systems are represented as graphs with STNs embedded in each node. Interactions include mass-action, enzymatic, allosteric and connectionist models. Reactions are translated into differential equations and can be solved numerically to generate predictive time courses or output as systems of equations that can be read by other programs. Cellerator simulations are fully extensible and portable to any operating system that supports Mathematica, and can be indefinitely nested within larger data structures to produce highly scaleable models.
Non-NASA Center↗
Thin membrane sensor with biochemical switch
A modular biosensor system for chemical or biological agent detection utilizes electrochemical measurement of an ion current across a gate membrane triggered by the reaction of the target agent with a recognition protein conjugated to a channel blocker. The sensor system includes a bioresponse simulator or biochemical switch module which contains the recognition protein-channel blocker conjugate, and in which the detection reactions occur, and a transducer module which contains a gate membrane and a measuring electrode, and in which the presence of agent is sensed electrically. In the poised state, ion channels in the gate membrane are blocked by the recognition protein-channel blocker conjugate. Detection reactions remove the recognition protein-channel blocker conjugate from the ion channels, thus eliciting an ion current surge in the gate membrane which subsequently triggers an output alarm. Sufficiently large currents are generated that simple direct current electronics are adequate for the measurements. The biosensor has applications for environmental, medical, and industrial use.
Case, George D.↗
Ceramide-1-phosphate transfer protein enhances lipid transport by disrupting hydrophobic lipid–membrane contacts
Cellular distributions of the sphingolipid ceramide-1-phosphate (C1P) impact essential biological processes. C1P levels are spatiotemporally regulated by ceramide-1-phosphate transfer protein (CPTP), which efficiently shuttles C1P between organelle membranes. Yet, how CPTP rapidly extracts and inserts C1P into a membrane remains unknown. Here, we devise a multiscale simulation approach to elucidate biophysical details of CPTP-mediated C1P transport. We find that CPTP binds a membrane poised to extract and insert C1P and that membrane binding promotes conformational changes in CPTP that facilitate C1P uptake and release. By significantly disrupting a lipid’s local hydrophobic environment in the membrane, CPTP lowers the activation free energy barrier for passive C1P desorption and enhances C1P extraction from the membrane. Upon uptake of C1P, further conformational changes may aid membrane unbinding in a manner reminiscent of the electrostatic switching mechanism used by other lipid transfer proteins. Insertion of C1P into an acceptor membrane, eased by a decrease in membrane order by CPTP, restarts the transfer cycle. Most notably, we provide molecular evidence for CPTP’s ability to catalyze C1P extraction by breaking hydrophobic C1P–membrane contacts with compensatory hydrophobic lipid–protein contacts. Our work, thus, provides biophysical insights into how CPTP efficiently traffics C1P between membranes to maintain sphingolipid homeostasis and, additionally, presents a simulation method aptly suited for uncovering the catalytic mechanisms of other lipid transfer proteins.
59 BASIC BIOLOGICAL SCIENCES↗
OpenAWSEM with Open3SPN2: A fast, flexible, and accessible framework for large-scale coarse-grained biomolecular simulations
We present OpenAWSEM and Open3SPN2, new cross-compatible implementations of coarse-grained models for protein (AWSEM) and DNA (3SPN2) molecular dynamics simulations within the OpenMM framework. These new implementations retain the chemical accuracy and intrinsic efficiency of the original models while adding GPU acceleration and the ease of forcefield modification provided by OpenMM’s Custom Forces software framework. By utilizing GPUs, we achieve around a 30-fold speedup in protein and protein-DNA simulations over the existing LAMMPS-based implementations running on a single CPU core. We showcase the benefits of OpenMM’s Custom Forces framework by devising and implementing two new potentials that allow us to address important aspects of protein folding and structure prediction and by testing the ability of the combined OpenAWSEM and Open3SPN2 to model protein-DNA binding. The first potential is used to describe the changes in effective interactions that occur as a protein becomes partially buried in a membrane. We also introduced an interaction to describe proteins with multiple disulfide bonds. Using simple pairwise disulfide bonding terms results in unphysical clustering of cysteine residues, posing a problem when simulating the folding of proteins with many cysteines. We now can computationally reproduce Anfinsen’s early Nobel prize winning experiments by using OpenMM’s Custom Forces framework to introduce a multi-body disulfide bonding term that prevents unphysical clustering. Our protein-DNA simulations show that the binding landscape is funneled towards structures that are quite similar to those found using experiments. In summary, this paper provides a simulation tool for the molecular biophysics community that is both easy to use and sufficiently efficient to simulate large proteins and large protein-DNA systems that are central to many cellular processes. These codes should facilitate the interplay between molecular simulations and cellular studies, which have been hampered by the large mismatch between the time and length scales accessible to molecular simulations and those relevant to cell biology.
59 BASIC BIOLOGICAL SCIENCES↗
Crystal structure of dihydrofolate reductase from the filarial nematode W. bancrofti in complex with NADPH and folate
Lymphatic filariasis is a debilitating illness with an estimated 50 million cases as of 2018. The majority of cases are caused by the parasitic worm W. bancrofti and additional cases by the worms B. malayi and B. timori. Dihydrofolate reductase (DHFR) is an established target in the treatment of cancer, bacterial, and protozoal infections and may be a potential target for drugs targeting parasitic worm infections, including filariasis. Recent studies have shown that known antifolate compounds, including methotrexate, inhibit the activity of W. bancrofti DHFR (WbDHFR). However, the absence of structural information for filarial DHFRs has limited the study of more in-depth structure-function relationships. We report the structure of WbDHFR complexed with NADPH and folate using X-ray diffraction data measured to 2.47 Å resolution. The structure of WbDHFR reveals the usual DHFR fold and is currently only the second nematode DHFR structure in the Protein Data Bank. The equilibrium dissociation constants for NADPH (90 ± 29 nM) and folate (23 ± 4 nM) were determined by equilibrium titrations. The interactions of known antifolates with WbDHFR were analyzed using molecular docking programs and molecular dynamics simulations. Antifolates with a hydrophobic core and extended linker formed favorable interactions with WbDHFR. These combined data should now facilitate the rational design of filarial DHFR inhibitors, which in turn can be used to determine whether DHFR is a viable drug target for filariasis and whether existing antifolates may be repurposed for its treatment.
60 APPLIED LIFE SCIENCES↗
Structural and biochemical insight into a modular β-1,4-galactan synthase in plants
Rhamnogalacturonan I (RGI) is a structurally complex pectic polysaccharide with a backbone of alternating rhamnose and galacturonic acid residues substituted with arabinan and galactan side chains. Galactan synthase 1 (GalS1) transfers galactose and arabinose to either extend or cap the ..beta..-1,4-galactan side chains of RGI, respectively. Here we report the structure of GalS1 from Populus trichocarpa, showing a modular protein consisting of an N-terminal domain that represents the founding member of a new family of carbohydrate-binding module, CBM95, and a C-terminal glycosyltransferase family 92 (GT92) catalytic domain that adopts a GT-A fold. GalS1 exists as a dimer in vitro, with stem domains interacting across the chains in a 'handshake' orientation that is essential for maintaining stability and activity. In addition to understanding the enzymatic mechanism of GalS1, here we gained insight into the donor and acceptor substrate binding sites using deep evolutionary analysis, molecular simulations and biochemical studies. Combining all the results, a mechanism for GalS1 catalysis and a new model for pectic galactan side-chain addition are proposed.
09 BIOMASS FUELS↗
Structural analyses of the PKA RIIβ holoenzyme containing the oncogenic DnaJB1-PKAc fusion protein reveal protomer asymmetry and fusion-induced allosteric perturbations in fibrolamellar hepatocellular carcinoma
When the J-domain of the heat shock protein DnaJB1 is fused to the catalytic (C) subunit of cAMP-dependent protein kinase (PKA), replacing exon 1, this fusion protein, J-C subunit (J-C), becomes the driver of fibrolamellar hepatocellular carcinoma (FL-HCC). Here, we use cryo-electron microscopy (cryo-EM) to characterize J-C bound to RIIβ, the major PKA regulatory (R) subunit in liver, thus reporting the first cryo-EM structure of any PKA holoenzyme. We report several differences in both structure and dynamics that could not be captured by the conventional crystallography approaches used to obtain prior structures. Most striking is the asymmetry caused by the absence of the second cyclic nucleotide binding (CNB) domain and the J-domain in one of the RIIβ:J-C protomers. Using molecular dynamics (MD) simulations, we discovered that this asymmetry is already present in the wild-type (WT) RIIβ 2 C 2 but had been masked in the previous crystal structure. This asymmetry may link to the intrinsic allosteric regulation of all PKA holoenzymes and could also explain why most disease mutations in PKA regulatory subunits are dominant negative. The cryo-EM structure, combined with small-angle X-ray scattering (SAXS), also allowed us to predict the general position of the Dimerization/Docking (D/D) domain, which is essential for localization and interacting with membrane-anchored A-Kinase-Anchoring Proteins (AKAPs). This position provides a multivalent mechanism for interaction of the RIIβ holoenzyme with membranes and would be perturbed in the oncogenic fusion protein. The J-domain also alters several biochemical properties of the RIIβ holoenzyme: It is easier to activate with cAMP, and the cooperativity is reduced. These results provide new insights into how the finely tuned allosteric PKA signaling network is disrupted by the oncogenic J-C subunit, ultimately leading to the development of FL-HCC.
59 BASIC BIOLOGICAL SCIENCES↗
The TPR domain of PgaA is a multifunctional scaffold that binds PNAG and modulates PgaB-dependent polymer processing
The synthesis of exopolysaccharides as biofilm matrix components by pathogens is a crucial factor for chronic infections and antibiotic resistance. Many periplasmic proteins involved in polymer processing and secretion in Gram-negative synthase dependent exopolysaccharide biosynthetic systems have been individually characterized. The operons responsible for the production of PNAG, alginate, cellulose and the Pel polysaccharide each contain a gene that encodes an outer membrane associated tetratricopeptide repeat (TPR) domain containing protein. While the TPR domain has been shown to bind other periplasmic proteins, the functional consequences of these interactions for the polymer remain poorly understood. Herein, we show that the C-terminal TPR region of PgaA interacts with the de-N-acetylase domain of PgaB, and increases its deacetylase activity. Additionally, we found that when the two proteins form a complex, the glycoside hydrolase activity of PgaB is also increased. To better understand structure-function relationships we determined the crystal structure of a stable TPR module, which has a conserved groove formed by three repeat motifs. Tryptophan quenching, mass spectrometry analysis and molecular dynamics simulation studies suggest that the crystallized TPR module can bind PNAG/dPNAG via its electronegative groove on the concave surface, and potentially guide the polymer through the periplasm towards the porin for export. Our results suggest a scaffolding role for the TPR domain that combines PNAG/dPNAG translocation with the modulation of its chemical structure by PgaB.
59 BASIC BIOLOGICAL SCIENCES↗
Functional Characterization and Structural Insights Into Stereoselectivity of Pulegone Reductase in Menthol Biosynthesis
Monoterpenoids are the main components of plant essential oils and the active components of some traditional Chinese medicinal herbs like Mentha haplocalyx Briq ., Nepeta tenuifolia Briq ., Perilla frutescens (L.) Britt and Pogostemin cablin (Blanco) Benth. Pulegone reductase is the key enzyme in the biosynthesis of menthol and is required for the stereoselective reduction of the Δ 2,8 double bond of pulegone to produce the major intermediate menthone, thus determining the stereochemistry of menthol. However, the structural basis and mechanism underlying the stereoselectivity of pulegone reductase remain poorly understood. In this study, we characterized a novel (-)-pulegone reductase from Nepeta tenuifolia ( Nt PR), which can catalyze (-)-pulegone to (+)-menthone and (-)-isomenthone through our RNA-seq, bioinformatic analysis in combination with in vitro enzyme activity assay, and determined the structure of (+)-pulegone reductase from ( Mp PR) by using X-ray crystallography, molecular modeling and docking, site-directed mutagenesis, molecular dynamics simulations, and biochemical analysis. We identified and validated the critical residues in the crystal structure of Mp PR involved in the binding of the substrate pulegone. We also further identified that residues Leu56, Val282, and Val284 determine the stereoselectivity of the substrate pulegone, and mainly contributes to the product stereoselectivity. This work not only provides a starting point for the understanding of stereoselectivity of pulegone reductases, but also offers a basis for the engineering of menthone/menthol biosynthetic enzymes to achieve high-titer, industrial-scale production of enantiomerically pure products.
32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗