Search NASA⌕ Search

SEARCH · Search NASA

Results for “in silico”

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

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

At least 109 records · Page 6
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

CATMoS: Collaborative Acute Toxicity Modeling Suite

Background: Humans are exposed to tens of thousands of chemical substances that need to be assessed for their potential toxicity. Acute systemic toxicity testing serves as the basis for regulatory hazard classification, labeling, and risk management. However, it is cost- and time-prohibitive to evaluate all new and existing chemicals using traditional rodent acute toxicity tests. In silico models built using existing data facilitate rapid acute toxicity predictions without using animals. Objectives: The U.S. Interagency Coordinating Committee on the Validation of Alternative Methods Acute Toxicity Workgroup organized an international collaboration to develop in silico models for predicting acute oral toxicity based on five different endpoints: LD50 value, U.S. Environmental Protection Agency hazard categories, Globally Harmonized System for Classification and Labelling hazard categories, very toxic chemicals (LD50 =50 mg/kg), and non-toxic chemicals (LD50 >2000 mg/kg). Methods: An acute oral toxicity data inventory for 11,992 chemicals was compiled, split into training and evaluation sets, and made available to 35 participating international research groups that submitted a total of 139 predictive models. Predictions that fell within the applicability domains of the submitted models were evaluated using external validation sets. These were then combined into consensus models to leverage strengths of individual approaches. Results: The resulting consensus predictions, which leverage the collective strengths of each individual model, form the Collaborative Acute Toxicity Modeling Suite (CATMoS). CATMoS demonstrated high performance in terms of accuracy and robustness when compared to in vivo results. Discussion: CATMoS is being evaluated by regulatory agencies for its utility and applicability as a potential replacement for in vivo rat acute oral toxicity studies. CATMoS predictions for over 800,000 chemicals have been made available via the NTP’s Integrated Chemical Environment. The models are also implemented in a free, standalone open-source tool, OPERA, which allows predictions of new and untested chemicals to be made.

63 RADIATION, THERMAL, AND OTHER ENVIRON. POLLUTAN↗
DOE OSTI2020

Evaluating kratom alkaloids using PHASE

Kratom is a botanical substance that is marketed and promoted in the US for pharmaceutical opioid indications despite having no US Food and Drug Administration approved uses. Kratom contains over forty alkaloids including two partial agonists at the mu opioid receptor, mitragynine and 7-hydroxymitragynine, that have been subjected to the FDA’s scientific and medical evaluation. However, pharmacological and toxicological data for the remaining alkaloids are limited. Therefore, we applied the Public Health Assessment via Structural Evaluation (PHASE) protocol to generate in silico binding profiles for 25 kratom alkaloids to facilitate the risk evaluation of kratom. PHASE demonstrates that kratom alkaloids share structural features with controlled opioids, indicates that several alkaloids bind to the opioid, adrenergic, and serotonin receptors, and suggests that mitragynine and 7-hydroxymitragynine are the strongest binders at the mu opioid receptor. Subsequently, the in silico binding profiles of a subset of the alkaloids were experimentally verified at the opioid, adrenergic, and serotonin receptors using radioligand binding assays. The verified binding profiles demonstrate the ability of PHASE to identify potential safety signals and provide a tool for prioritizing experimental evaluation of high-risk compounds.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2021

Robust Molecular Predictive Methods for Novel Polymer Discovery and Applications

Polymeric materials are ubiquitous in modern society and they play an instrumental role in almost all industries, undoubtedly including the energy and environment sectors. Increased demand of energy and awareness to sustainability both necessitates the development of novel polymers with enhanced properties. Unfortunately, their structural and behavioral complexity render such discovery challenging and impeded. To address this problem, scientists are developing various computational modeling techniques and leveraging their power to depict the relationship between structural characteristics of polymers and their properties (such as rheological behaviors), and use such prediction to guide the design and syntheses of novel polymeric materials with enhanced performances. Unfortunately, predicting the relationships between polymer structure and composition with rheological properties via atomistic modeling is still a major challenge because of the extended time and length scales involved. Studying dynamic shear viscosity and linear viscoelasticity using molecular models requires capabilities that have been elusive, including representation of large molecular weight chains with an effective internal scale capable of describing entanglement, shear-rates that are in the s-1 scale with accurate quantitative stresses, and chemically-realistic combinations of both homogeneous and heterogeneous systems. Motivated by these unmet challenges, the overall technical objective of this DOE-STTR Phase II project is to develop robust molecular predictive methods for advanced polymer discovery and applications and especially for designing and demonstrating the “smart” polymer-based waterflooding enhanced oil recovery (EOR) process. In particular, we apply state-of-the-art molecular modeling methods developed by our academic partner, Materials Stimulation Center (MSC) at California Institute of Technology (Caltech), to facilitate and accelerate the experimental discovery processes. During the Phase I of this project, we had focused on development and demonstration of the molecular modeling methods to describe rheological properties of non-Newtonian polymer fluids, and to improve our fundamental understandings of shear-thickening mechanism and kinetics. In Phase II, we further apply the theoretical models to guide our experimental programs to improve our design of smart rheology modifier (SRM) polymers and their optimization for EOR. Specifically, we have three objectives in the Phase II study: (1) to further improve out computational modeling methods, coupling with the advanced machine learning algorithms; (2) to develop cost-effective and efficient SRM-flooding process suitable for EOR applications under typical reservoir conditions; and (3) to further explore the application of our molecular predictive models for innovative material discovery in other industrial applications. The recent development of our multiscale predictive framework allows the successful prediction of rheological properties from the chemical structure for polymers of experimentally relevant molecular weights, and provides an in-silico machine learning engine for screening novel compositions and structures with optimized non-Newtonian response, required for both shear-thinning and shear-thickening applications. Our framework provides: (1) procedures and tools for systematic coarsening from atomistic models and reverse mapping of coarse-grain models to atomistic, (2) unique ab initio methods to characterize the atomistic origin of colloidal and interfacial interactions and phenomena, (3) systematic structure and composition builders based on practical descriptors that drive rheological changes in polymer melts and diluted polymer mixtures, (4) a rheological properties engine capable of predicting viscosity in the zero-shear limit and under realistic dynamic conditions (for shear-rates commensurate with experiments) for large heterogeneous systems, (5) coarse-grain force fields with improved non-bond descriptions based on accurate quantum mechanics, (6) an in-silico screening machine learning engine that feeds from the systematic model builders to cover the descriptors search space, computes the rheological properties from converged trajectories spanning sub-milliseconds and ranks them for each structure/composition using an automated viscosity-vs-shear rate fitness function that can be tuned for shear-thickening, shear-thinning and other rheological responses.

02 PETROLEUM↗
DOE OSTI2023

gRNA-SeqRET: a universal tool for targeted and genome-scale gRNA design and sequence extraction for prokaryotes and eukaryotes

High-throughput genetic screening is frequently employed to rapidly associate gene with phenotype and establish sequence-function relationships. With the advent of CRISPR technology, and the ability to functionally interrogate previously genetically recalcitrant organisms, non-model organisms can be investigated using pooled guide RNA (gRNA) libraries and sequencing-based assays to quantitatively assess fitness of every targeted locus in parallel. To aid the construction of pooled gRNA assemblies, we have developed an in silico design workflow for gRNA selection using the gRNA Sequence Region Extraction Tool (gRNA-SeqRET). Built upon the previously developed CCTop, gRNA-SeqRET enables automated, scalable design of gRNA libraries that target user-specified regions or whole genomes of any prokaryote or eukaryote. Additionally, gRNA-SeqRET automates the bulk extraction of any regions of sequence relative to genes or other features, aiding in the design of homology arms for insertion or deletion constructs. We also assess in silico the application of a designed gRNA library to other closely related genomes and demonstrate that for very closely related organisms Average Nucleotide Identity (ANI) > 95% a large fraction of the library may be of relevance. The gRNA-SeqRET web application pipeline can be accessed at https://grna.jgi.doe.gov. The source code is comprised of freely available software tools and customized Python scripts, and is available at https://bitbucket.org/berkeleylab/grnadesigner/src/master/ under a modified BSD open-source license (https://bitbucket.org/berkeleylab/grnadesigner).

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Structural analysis and functional evaluation of the disordered ß–hexosyltransferase region from Hamamotoa (Sporobolomyces) singularis

Hamamotoa (Sporobolomyces) singularis codes for an industrially important membrane bound ß-hexosyltransferase (BHT), (BglA, UniprotKB: Q564N5) that has applications in the production of natural fibers such as galacto-oligosaccharides (GOS) and natural sugars found in human milk. When heterologously expressed by Komagataella phaffii GS115, BHT is found both membrane bound and soluble secreted into the culture medium. In silico structural predictions and crystal structures support a glycosylated homodimeric enzyme and the presence of an intrinsically disordered region (IDR) with membrane binding potential within its novel N-terminal region (1–110 amino acids). Additional in silico analysis showed that the IDR may not be essential for stable homodimerization. Thus, we performed progressive deletion analyses targeting segments within the suspected disordered region, to determine the N-terminal disorder region’s impact on the ratio of membrane-bound to secreted soluble enzyme and its contribution to enzyme activity. The ratio of the soluble secreted to membrane-bound enzyme shifted from 40% to 53% after the disordered N-terminal region was completely removed, while the specific activity was unaffected. Furthermore, functional analysis of each glycosylation site found within the C-terminal domain revealed reduced total secreted protein activity by 58%–97% in both the presence and absence of the IDR, indicating that glycosylation at all four locations is required by the host for the secretion of active enzyme and independent of the removed disordered N-terminal region. Overall, the data provides evidence that the disordered region only partially influences the secretion and membrane localization of BHT.

59 BASIC BIOLOGICAL SCIENCES↗
NASA NTRS2019

Circulating miRNA Signature Predicts Health Risks Associated with Radiation and Microgravity

The many known health risks currently associated with space travel include increased risk of cardiovascular disease, cancer, central nervous system related diseases, muscle degeneration, and changes with host-gut microbiome interactions that can have profound impact with these and other health risks. The majority of the risk from space travel stem of the two components of the space environment which are microgravity and radiation. From our earlier work (Beheshti et al, PLOS One, 2018), we predicted that there is a systemic component of the host that causes general increased health risks due to spaceflight driven by a circulating microRNA (miRNA) signature consisting of 13 miRNAs that directly regulates both p53 and TGF1. MiRNAs are small non-coding RNA molecules with a negative and post-transcriptional regulation on gene expression) are increasingly recognized as major systemic regulators of responses to stressors, including microgravity, oxidative stress, and DNA damage. In addition, due to the size and stability of miRNAs, it is known that miRNAs can circulate throughout the body and have been found in the majority of the bodily fluids including blood, urine, saliva, and tears. Here, we start to dissect the actual impact of this miRNA signature on both the radiation and microgravity components and prove that this miRNA signature actually exists in the circulation of a host. To achieve this, we obtained multiple tissues including, serum, liver, and spleen and utilizing droplet digital PCR (ddPCR), we start to show how this circulating miRNA signature impacts which component of the spaceflight. The tissue was obtained from experiments performed on C57BL/6 male mice (N=10 for each condition) that were hindlimb unloaded (HU) to simulated microgravity, irradiated with 2Gy gamma (IR), HU plus IR, and control mice under normal conditions. It was shown that these miRNAs were present in the serum as predicted by the in silico prediction from our earlier predictions. The HU vs Controls show significant increases of the predicted miRNAs in the serum for more than half of the miRNA signature, with remaining miRNAs increasing comparing to the controls close to statistical significance. IR vs control mice showed increases for the miRNAs, but not has pronounced as the HU conditions. Finally, the combination of the HU+IR vs controls showed increases for the majority of the miRNA signature. The data indicates that the miRNA signature originally predicted through in silico methods is mainly associated with the microgravity component and is circulating throughout the host resulting in a systemic impact of the miRNAs on the host. These miRNAs are shown in the literature to potentially increase health risks associated with several diseases. In addition, we have begun testing the potential of utilizing antagonists to this miRNA signature to act as a potential countermeasure to mitigate radiation impact on the organism. This work demonstrates for the first time the potential of a minimally invasive novel biomarker and countermeasure that can be used to mitigate both radiation and microgravity effects.

Beheshti, Afshin↗
NASA NTRSDate not supplied

Evaluation of Human Spaceflight-Related Tissue Weight Relief Using Whole Body Finite Element Model Simulations

Tissue Weight Relief (TWR) is a physiological condition observed in human spaceflight. It not only impacts the injury biomechanics of soft tissue but also the physiological responses of the cardiovascular system both due to fluid redistribution and the effect of tissue-related transmural pressure on the large venous blood vessels. Understanding the effects of tissue weight relief is especially important because of the role it may play in understanding the cause of Space Associated Neuro-Ocular Syndrome (SANS). SANS can be characterized by a number of ocular changes which reduce visual acuity and SANS related symptoms occur in up to 51% of astronauts. A prevailing theory for the causation of SANS is that of headward (cephalad) fluid shift and a prolonged increase of Intracranial Pressure (ICP) similar to intracranial hypertension, which is not fully supported by the experimental data or astronaut symptom reporting. However, it is still believed that SANS is caused by a pressure change in the eye and the surrounding tissues. It has been proposed that TWR plays a substantial role in affecting internal pressures and fluid shifts in microgravity. In this effort, two whole-body Finite Element (FE) models – Elemance and THUMS – are used to ascertain the microgravity-associated TWR of the musculature surrounding the lower body veins. Elemance and THUMS are physics-based computational models that have been validated and verified for several automotive and domestic applications, and as such, can simulate the relief of soft tissue weight due to changes in the gravitation vector. Specifically, the current effort modeled the transition of the gravitational vector from 1G to 0G, applied across the whole-body model in a supine position. For each 1G to 0G transition simulation, the lower body vein’s transmural pressure-time profile was extracted and averaged around the anterior portion of the thigh muscle. The ascertained transmural pressure changes from 1G to 0G transition are given in Figure 1 for the Elemance and the THUMS FE models. The transmural pressure changes of 10 mmHg and 21 mmHg are in the same order of magnitude as Lu’s value of 44 mmHg. It is to be noted that Lu implemented a 0D to 1D lumped parameter model and the Elemance and THUMS are 3D higher order computational models. This proof-of-concept approach demonstrates that TWR pressure can be adequately estimated with in silico techniques however, further in silico investigations need to be conducted to address the unique contributions to the transmural pressure from each of the computational models.

Finite element modeling↗
DOE OSTI2021

Elucidation of Aromatic Catabolic Pathways in White-Rot Fungi

This project aims to investigate the hypothesis that white-rot fungi can simultaneously depolymerize lignin extracellularly and catabolize depolymerization products intracellularly as carbon and energy sources. Evaluating this hypothesis will provide deeper understanding of the role of white-rot fungi in facilitating carbon sequestration in Nature. Additionally, identifying the most promising fungal strains for lignin turnover and catabolism will catalyze future efforts in genetic tool development to enable metabolic engineering in white-rot fungi for lignin bioconversion to bioproducts. Lignin is the second most abundant plant-based biopolymer on Earth and represents up to 40% of the energy density of lignocellulosic biomass. Even though lignin is a massive natural carbon and energy reservoir, only a small group of basidiomycete fungi, namely white-rot fungi (WRF), have evolved the ability to efficiently depolymerize and mineralize lignin to CO2 and H2O. Considerable research efforts have been undertaken to understand how WRF depolymerize lignin but the biochemical reactions that convert lignin into CO2 have been largely neglected. In fact, it is unclear if WRF intracellularly catabolize lignin-derived aromatic compounds to utilize them as a carbon and energy source, or rather if lignin is depolymerized and mineralized extracellularly merely to facilitate access to cellulose and hemicellulose for use as a primary carbon source. To date, we have employed 13C-isotope labeling, systems biology approaches, and in vitro enzyme assays to definitively demonstrate that two WRF, Trametes versicolor and Gelatoporia (Ceriporiopsis) subvermispora, funnel carbon from lignin-derived aromatic compounds into central carbon metabolism via intracellular catabolic pathways [1]. Specifically, 13C-isotopic labeling approaches showed that these WRF utilize poplar-derived aromatic compounds (e.g. 4-hydroxybenzoic acid (4-HBA)) as a carbon source. In silico genome analysis led us to hypothesize a complete catabolic pathway for 4-HBA and identify multiple homologous sequences for enzymes with putative oxidative decarboxylase, hydroxylase, and ring-opening dioxygenase activities, which are among the main biochemical reactions acting on aromatic compounds. Spatial and differential proteomic and metabolomic analyses supported the proposed catabolic pathways and showed alternative catabolic steps in T. versicolor that were not present in G. subvermispora. Based on the in silico, proteomics, and transcriptomics results, we down-selected enzymes for further in vitro characterization, and we have assigned a function to six fungal enzymes (including oxidative decarboxylases, hydroxylases, and ring-opening dioxygenases). Interestingly, even though we selected homologous enzyme pairs from both WRF with similar -omics trends, in a few cases only one of the studied fungi showed activity for the proposed substrate. Based on all the observations from this study [1], we hypothesized that 4-HBA preferentially undergoes oxidative decarboxylation to hydroquinone and subsequent hydroxylation to 1,2,4-benzenetriol in G. subvermispora before ring cleavage, whereas 4-HBA would preferentially undergo hydroxylation to protocatechuate and further oxidative decarboxylation to 1,2,4-benzenetriol in T. versicolor. Examining additional decarboxylases and hydroxylases as well as enzymes from other protein families that can perform the same or similar oxidative reactions, such as cytochromes P450 with aromatic hydroxylation activity, will also be key for elucidating enzyme preferences for specific substrates. Overall, this work forms the foundation of a new research area based on lignin catabolism by WRF, which could be further exploited to convert the undervalued biopolymer lignin into value-added compounds.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2025

Paired Neural Network for Matching Experimental and Predicted Infrared Spectra

Here, we present a novel machine learning (ML)-based scoring technique for determining the similarity between experimental and predicted infrared (IR) spectra for identification purposes. IR spectroscopy is a powerful technique used to identify the molecular structure and composition of a sample by measuring the unique vibrational frequency pattern of the molecule’s functional groups. Molecular identifications are often made by comparing experimental and reference spectra. However, the limited number of reference spectra available in spectral libraries can confound the identification process. Alternative identification procedures rely on in silico techniques to simulate spectra for a wide range of molecules. However, scoring spectral similarity between an experimental query and computationally predicted reference remains a significant challenge. Our proposed ML-based scoring technique overcomes these barriers by accurately and efficiently determining spectral similarity.

Neural Network↗
DOE OSTI2022

Interleukin-2 superkines by computational design

Significance While computational engineering of therapeutic proteins is a desirable goal, in practice the optimization of protein–protein interactions requires substantial experimental intervention. We present here a computational approach that focuses on stabilizing core protein structures rather than engineering the protein–protein interface. Using this approach, we designed thermostabilized interleukin-2 (IL-2) variants that bind tightly to their receptor without experimental optimization, mimicking the properties of the yeast-display engineered IL-2 variant “super-2.” Our results suggest that structure-guided stabilization may be a general method for in silico affinity maturation of protein–protein interactions.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Heracles: Predictive Tools for Opioid Crisis Intervention - m/q Initiative Project Report

The opioid crisis in the United States is being fueled primarily by fentanyl and its molecular analogs, which can be anywhere from 50 to 1,000 times more potent than morphine. Fentanyl itself is straightforward to synthesize; furthermore, the structure is such that fentanyl’s flexible, rotatable side chains are easy to modify to create new analogs. Reference-free computational techniques to predict and identify new fentanyls have the potential to provide a desperately needed preemptive advantage to regulatory stakeholders and toxicologists. The computational pipeline Heracles was developed with this preemptive advantage in mind. Heracles has two primary components: 1) the creation of an in silico library of putative fentanyl analogs, and 2) a downselection pipeline to prioritize generated fentanyl analogs predicted to be potent and easy to synthesize. Experimental observables were also predicted for prioritized analogs, with validation of the observables begun. Heracles has demonstrated potential to aid in the advancement of reference-free paradigms while providing new tools to first responders and other stakeholders attempting to mitigate the opioid crisis.

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

A goldilocks computational protocol for inhibitor discovery targeting DNA damage responses including replication-repair functions

While many researchers can design knockdown and knockout methodologies to remove a gene product, this is mainly untrue for new chemical inhibitor designs that empower multifunctional DNA Damage Response (DDR) networks. Here, we present a robust Goldilocks (GL) computational discovery protocol to efficiently innovate inhibitor tools and preclinical drug candidates for cellular and structural biologists without requiring extensive virtual screen (VS) and chemical synthesis expertise. By computationally targeting DDR replication and repair proteins, we exemplify the identification of DDR target sites and compounds to probe cancer biology. Our GL pipeline integrates experimental and predicted structures to efficiently discover leads, allowing early-structure and early-testing (ESET) experiments by many laboratories. By employing an efficient VS protocol to examine protein-protein interfaces (PPIs) and allosteric interactions, we identify ligand binding sites beyond active sites, leveraging in silico advances for molecular docking and modeling to screen PPIs and multiple targets. A diverse 3,174 compound ESET library combines Diamond Light Source DSI-poised, Protein Data Bank fragments, and FDA-approved drugs to span relevant chemotypes and facilitate downstream hit evaluation efficiency for academic laboratories. Two VS per library and multiple ranked ligand binding poses enable target testing for several DDR targets. This GL library and protocol can thus strategically probe multiple DDR network targets and identify readily available compounds for early structural and activity testing to overcome bottlenecks that can limit timely breakthrough drug discoveries. By testing accessible compounds to dissect multi-functional DDRs and suggesting inhibitor mechanisms from initial docking, the GL approach may enable more groups to help accelerate discovery, suggest new sites and compounds for challenging targets including emerging biothreats and advance cancer biology for future precision medicine clinical trials.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Investigation and Expression Analysis of R2R3-MYBs and Anthocyanin Biosynthesis-Related Genes during Seed Color Development of Common Bean ( Phaseolus vulgaris )

Anthocyanins are responsible for the coloration of common bean seeds, and their accumulation is positively correlated with the expression level of anthocyanin biosynthetic genes. The MBW (MYB-bHLH-WD40) complex is thought to regulate the expression of these genes, and MYB proteins, which are a key factor in activating anthocyanin pathway genes, have been identified in several plants. This study demonstrated gene structures, chromosomal placements, gene duplications of R2R3-MYBs, miRNAs associated with R2R3-MYBs, and the interaction of these genes with other flavonoid regulatory genes. qRT-PCR was used to investigate the role of specific R2R3-MYBs and flavonoid genes in common bean seed color development. As a result of a comprehensive analysis with the help of in silico tools, we identified 160 R2R3-MYB genes in the common bean genome. We divided these genes into 16 classes on the basis of their intron-exon and motif structures. Except for three, the rest of the common bean R2R3-MYB members were distributed to all chromosomes with different densities, primarily located on chromosomes 3 and 8. We identified a total of 44 duplicated gene pairs dispersed across 11 chromosomes and evolved under purifying selection (Ka/Ks < 1), 19 of which were derived from a whole-genome duplication. Our research uncovered 25 putative repressor PvMYB proteins that contain the EAR motif. Additionally, fifty different cis-regulatory elements regulated by light, stress, and hormone were identified. Within the genome of the common bean, we discovered a total of 36 microRNAs that target a total of 72 R2R3-MYB transcripts. The effect of 16 R2R3-MYB genes and 16 phenylpropanoid pathway genes, selected on the basis of their interaction in the protein-protein interaction map, playing role in the regulation of seed coat color development was evaluated using qRT-PCR in 5 different tissues at different developmental stages. The results revealed that these specific genes have different expression levels during different developmental periods, with higher levels in the pod filling and early pod stages than in the rest of the developmental periods. Furthermore, it was shown that PvTT8 (bHLH), PvTT2 (PvMYB42), PvMYB113, PvTTG1, and PvWD68 genes have effects on the regulation of seed coat color. The findings of this study, which is the first to use whole-genome analysis to identify and characterize the R2R3-MYB genes in common bean, may serve as a reference for future functional research in the legume.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Lipocalin Blc is a potential heme-binding protein

Lipocalins are a superfamily of functionally diverse proteins defined by a well-conserved tertiary structure despite variation in sequence. Lipocalins bind and transport small hydrophobic molecules in organisms of all kingdoms. However, there is still uncertainty regarding the function of some members of the family, including bacterial lipocalin Blc from Escherichia coli. Here, we present evidence that lipocalin Blc may be involved in heme binding, trans-periplasmic transport, or heme storage. This conclusion is supported by a cocrystal structure, mass-spectrometric data, absorption titration, and in silico analysis. Binding of heme is observed at low micromolar range with one-to-one ligand-to-protein stoichiometry. However, the absence of classical coordination to the iron atom leaves the possibility that the primary ligand of Blc is another tetrapyrrole.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Structure- and Interaction-Based Design of Anti-SARS-CoV-2 Aptamers

Aptamer selection against novel infections is a complicated and time-consuming approach. Synergy can be achieved by using computational methods together with experimental procedures. In this study, we aim to develop a reliable methodology for a rational aptamer in silico et vitro design. The new approach combines multiple steps: (1) Molecular design, based on screening in a DNA aptamer library and directed mutagenesis to fit the protein tertiary structure; (2) 3D molecular modeling of the target; (3) Molecular docking of an aptamer with the protein; (4) Molecular dynamics (MD) simulations of the complexes; (5) Quantum-mechanical (QM) evaluation of the interactions between aptamer and target with further analysis; (6) Experimental verification at each cycle for structure and binding affinity by using small-angle X-ray scattering, cytometry, and fluorescence polarization. By using a new iterative design procedure, structure- and interaction-based drug design (SIBDD), a highly specific aptamer to the receptor-binding domain of the SARS-CoV-2 spike protein, was developed and validated. The SIBDD approach enhances speed of the high-affinity aptamers development from scratch, using a target protein structure. The method could be used to improve existing aptamers for stronger binding. This approach brings to an advanced level the development of novel affinity probes, functional nucleic acids. It offers a blueprint for the straightforward design of targeting molecules for new pathogen agents and emerging variants.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2022

Calibration and Validation of a Mechanistic COVID-19 Model for Translational Quantitative Systems Pharmacology – A Proof-of-Concept Model Development for Remdesivir

With the ongoing global pandemic of coronavirus disease 2019 (COVID-19), there is an urgent need to accelerate the traditional drug development process. Many studies identified potential COVID-19 therapies based on promising nonclinical data. However, the poor translatability from nonclinical to clinical settings has led to failures of many of these drug candidates in the clinical phase. In this study, we propose a mechanism-based, quantitative framework to translate nonclinical findings to clinical outcome. Adopting a modularized approach, this framework includes an in silico disease model for COVID-19 (virus infection and human immune responses) and a pharmacological component for COVID-19 therapies. Here, the disease model was able to reproduce important longitudinal clinical data for patients with mild and severe COVID-19, including viral titer, key immunological cytokines, antibody responses, and time courses of lymphopenia. Using remdesivir as a proof-of-concept example of model development for the pharmacological component, we developed a pharmacological model that describes the conversion of intravenously administered remdesivir as a prodrug to its active metabolite nucleoside triphosphate through intracellular metabolism and connected it to the COVID-19 disease model. After being calibrated with the placebo arm data, our model was independently and quantitatively able to predict the primary endpoint (time to recovery) of the remdesivir clinical study, Adaptive Covid-19 Clinical Trial (ACTT). Our work demonstrates the possibility of quantitatively predicting clinical outcome based on nonclinical data and mechanistic understanding of the disease and provides a modularized framework to aid in candidate drug selection and clinical trial design for COVID-19 therapeutics.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2023

Thermal insulation versus capacitance: A simulation experiment comparing effects of shade and hyporheic exchange on daily and seasonal stream temperature cycles

In streams where water temperatures stress native biota, management of riparian shade or hyporheic exchange are both considered viable management strategies for reducing the peaks of daily and seasonal stream channel temperature cycles. Although shade and hyporheic exchange may have similar effects on stream temperatures, their mechanisms differ. Improved understanding of the heat-exchange mechanisms influenced by shade and hyporheic exchange will aid in the appropriate application of either stream temperature management strategy. To illustrate a conceptual model highlighting shade as ‘thermal insulation’ and hyporheic exchange imparting ‘thermal capacitance’ to a stream reach, we conducted an in-silico simulation modelling experiment increasing shade or hyporheic exchange parameters on an idealized, hypothetical stream. We assessed the potential effects of increasing shade or hyporheic exchange on a stream reach using an established process-based heat-energy budget model of stream-atmosphere heat exchange and incorporated an advection-driven hyporheic heat exchange routine. The model tracked heat transport through the hyporheic zone and exchange with the stream channel, while including the effects of hyporheic water age distribution on upwelling hyporheic temperatures. Results showed that shade and hyporheic exchange similarly damped diurnal temperature cycles and differentially altered seasonal cycles of our theoretical stream. In winter, hyporheic exchange warmed simulated channel temperatures whereas shade had little effect. In summer, both shade and hyporheic exchange cooled channel temperatures, though the effects of shade were more pronounced. Our simple-to-grasp analogies of ‘thermal insulation’ for shade effects and ‘thermal capacitance’ for hyporheic exchange effects on stream temperature encourage more accurate conceptualization of complex, dynamic heat exchange processes among the atmosphere, stream channel, and alluvial aquifer.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Phaseolus vulgaris SUT1.1 is a high affinity sucrose–proton co–transporter

Plant sucrose transporters are required for phloem loading, and therefore are essential for plant growth and development. In common beans (Phaseolus vulgaris) there are only two sucrose transporters functionally characterized. Through a previous RNA-seq study, we identified a putative sucrose transporter in common bean, which we hypothesize to function in import of sucrose into plant cells. In silico analysis revealed that PvSUT1.1 is a putative sucrose-proton co-transporter distinct from other characterized sucrose transporters in common bean indicating that this is a previously undescribed transporter protein in beans. Further analysis revealed that PvSUT1.1 shares high protein sequence homology to the phloem loader Arabidopsis SUC2; both have 12 transmembrane domains, a typical characteristic of plant sucrose transporters. Heterologous expression in yeast further showed PvSUT1.1 to be functional and it imported sucrose into yeast cells with a K m of 0.7 mM sucrose. Import of sucrose through PvSUT1.1 is also pH-dependent with highest uptake at pH 4.0, and activity is lost in the presence of the uncoupler carbonyl cyanide 3-chlorophenylhydrazone. Consistent with identification of PvSUT1.1 as a Type I transporter, PvSUT1.1 also transports esculin. Finally, PvSUT1.1 showed expression in multiple tissues and the protein was localized to the plasma membrane. The results show that PvSUT1.1 is a sucrose transporter that is probably involved in the uptake of sucrose into source and sink cells. The potential role of PvSUT1.1 in leaf phloem loading of sucrose in common beans and its importance in heat tolerance of reproductive tissues are further discussed.

59 BASIC BIOLOGICAL SCIENCES↗