Search NASA⌕ Search

SEARCH · Search NASA

Results for “Network motifs”

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.

142 records · Page 8

Coordination Stoichiometry Effects on the Binding Hierarchy of Histamine and Imidazole–M 2+ Complexes

Histidine–M 2+ coordination bonds are a recognized bond motif in biogenic materials with high hardness and extensibility, which has led to growing interest in their use in soft materials for mechanical function. However, the effect of different metal ions on the stability of the coordination complex remains poorly understood, complicating their implementation in metal-coordinated polymer materials. Herein, rheology experiments and density functional theory calculations are used to characterize the stability of coordination complexes and establish the binding hierarchy of histamine and imidazole with Ni 2+ , Cu 2+ , and Zn 2+ . It is found that the binding hierarchy is driven by the specific affinity of the metal ions to different coordination states, which can be macroscopically tuned by changing the metal-to-ligand stoichiometry of the metal-coordinated network. These findings facilitate the rational selection of metal ions for optimizing the mechanical properties of metal-coordinated materials.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Electronic Structure Progression across the ACu 2 Q 2 (MQ 2 ) n Semiconductor Series

Moving beyond the engineering of known materials and elemental substitution within common structure types is critical for designing unique properties. Here, we present a new homologous series, ACu 2 Q 2 (MQ 2 ) n (A = Sr, Ba, 2Na; MQ 2 = ZrS2, HfSe 2 ), establishing 11 new members. In β-BaCu 2 Q 2 and Na 2 Cu 2 Se 2 (n = 0), the [Cu 2 Q 2 ] 2– motifs extend in two dimensions, whereas those in α-BaCu 2 Q 2 extend in three dimensions. Hence, there are two structural evolutions within the ACu 2 Q 2 (MQ 2 )n family driven by the polymorphism of the host structures. MQ 2 (n → ∞) displays 2D layers of edge-sharing [MQ 6 ] 8– octahedra that are 1-octahedron-thick and connected via van der Waals bonding. Each insertion of MQ 2 into ACu 2 Q 2 incorporates [MQ 6 ] 8– octahedra extending infinitely in one direction, confined to being 1-octahedron-thick in the second direction, with n controlling the number of [MQ 6 ] 8– octahedra in the third direction. Therefore, increasing n predictively expands the [Cu 2 MnQ 2n+2 ] 2– network in the direction controlled by n and relative to the A + /A 2+ ions. We demonstrate that for a given set of elements, one can enforce a “Host(Insertion)n” formula to systematically evolve both crystal and electronic structures from the host (n = 0) to the insertion (n → ∞) materials through intermediate values of n. Specifically, the energetic misalignment of the electronic band extrema of the parents predictively determines the band extrema and the resulting band gaps of all intermediate n members.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Revealing the Reaction Path of UVC Bond Rupture in Cyclic Disulfides with Ultrafast X-ray Scattering

Disulfide bonds are ubiquitous molecular motifs that influence the tertiary structure and biological functions of many proteins. Yet it is well known that the disulfide bond is photolabile when exposed to UVC radiation. The deep-UV induced S-S bond fragmentation kinetics on very fast time scales are especially pivotal to fully understand the photostability and photodamage repair mechanisms in proteins. In 1,2-dithiane, the smallest saturated cyclic molecule that mimics biologically active species with S-S bonds, we investigate the photochemistry upon 200 nm excitation by femtosecond time-resolved X-ray scattering in the gas phase using an X-ray free electron laser. In the femtosecond time domain, we discover a very fast reaction that generates molecular fragments with one and two sulfur atoms. On picosecond and nanosecond timescales, a complex network of reactions unfolds that, ultimately, completes the sulfur dissociation from the parent molecule.

74 ATOMIC AND MOLECULAR PHYSICS↗

Evolutionary history, potential intermediate animal host, and cross-species analyses of SARS-CoV-2

To investigate the evolutionary history of the recent outbreak of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) in China, a total of 70 genomes of virus strains from China and elsewhere with sampling dates between 24 December 2019 and 3 February 2020 were analyzed. To explore the potential intermediate animal host of the SARS-CoV-2 virus, we reanalyzed virome data sets from pangolins and representative SARS-related coronaviruses isolates from bats, with particular attention paid to the spike glycoprotein gene. We performed phylogenetic, split network, transmission network, likelihood-mapping, and comparative analyses of the genomes. Additionally, based on Bayesian time-scaled phylogenetic analysis using the tip-dating method, we estimated the time to the most recent common ancestor and evolutionary rate of SARS-CoV-2, which ranged from 22 to 24 November 2019 and 1.19 to 1.31 × 10 -3 substitutions per site per year, respectively. Our results also revealed that the BetaCoV/bat/Yunnan/RaTG13/2013 virus was more similar to the SARS-CoV-2 virus than the coronavirus obtained from the two pangolin samples (SRR10168377 and SRR10168378). We also identified a unique peptide (PRRA) insertion in the human SARS-CoV-2 virus, which may be involved in the proteolytic cleavage of the spike protein by cellular proteases, and thus could impact host range and transmissibility. Interestingly, the coronavirus carried by pangolins did not have the RRAR motif. Therefore, we concluded that the human SARS-CoV-2 virus, which is responsible for the recent outbreak of COVID-19, did not come directly from pangolins.

59 BASIC BIOLOGICAL SCIENCES↗

Understanding the Impact of Data Staging for Coupled Scientific Workflows

We report the rate of data generated by cutting-edge experimental science facilities and large-scale simulations enabled by current high-performance computing (HPC) systems has continued to grow at a far greater pace than the development of the network and storage capabilities on which these systems rely. To cope with this challenge, scientist are moving toward the creation of autonomous experiments and HPC simulations using machine learning. However, efficiently moving, storing, and processing large amounts of data away from the point of origin presents an incredible challenge. In-memory computing, in situ analysis, data staging, and data streaming are recognized viable alternatives to traditional file-based methods for transferring data between coupled workflows. However, the performance trade-offs and limitations for these methods are not fully understood when used in HPC applications. This article presents a comprehensive performance assessment of the current solutions for data staging when applied to applications that are not necessary I/O intensive which makes them not ideal candidates for these methods. Our study is based on experiments running at scale on Oak Ridge National Laboratory's Summit supercomputer using applications and simulations that cover typical computational motifs and patterns. We investigated the usability and cost/benefit trade-offs of staging algorithms for HPC applications under different scenarios and highlight opportunities for optimizing the dataflow between coupled simulation workflows.

97 MATHEMATICS AND COMPUTING↗

DeepGRN: prediction of transcription factor binding site across cell-types using attention-based deep neural networks

Abstract Background Due to the complexity of the biological systems, the prediction of the potential DNA binding sites for transcription factors remains a difficult problem in computational biology. Genomic DNA sequences and experimental results from parallel sequencing provide available information about the affinity and accessibility of genome and are commonly used features in binding sites prediction. The attention mechanism in deep learning has shown its capability to learn long-range dependencies from sequential data, such as sentences and voices. Until now, no study has applied this approach in binding site inference from massively parallel sequencing data. The successful applications of attention mechanism in similar input contexts motivate us to build and test new methods that can accurately determine the binding sites of transcription factors. Results In this study, we propose a novel tool (named DeepGRN) for transcription factors binding site prediction based on the combination of two components: single attention module and pairwise attention module. The performance of our methods is evaluated on the ENCODE-DREAM in vivo Transcription Factor Binding Site Prediction Challenge datasets. The results show that DeepGRN achieves higher unified scores in 6 of 13 targets than any of the top four methods in the DREAM challenge. We also demonstrate that the attention weights learned by the model are correlated with potential informative inputs, such as DNase-Seq coverage and motifs, which provide possible explanations for the predictive improvements in DeepGRN. Conclusions DeepGRN can automatically and effectively predict transcription factor binding sites from DNA sequences and DNase-Seq coverage. Furthermore, the visualization techniques we developed for the attention modules help to interpret how critical patterns from different types of input features are recognized by our model.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Gaseous cyclodextrin- closo -dodecaborate complexes χCD·B 12 X 12 2– (χ = α, β, and γ; X = F, Cl, Br, and I): electronic structures and intramolecular interactions

A fundamental understanding of cyclodextrin- closo -dodecaborate inclusion complexes is of great interest in supramolecular chemistry. Herein, we report a systematic investigation on the electronic structures and intramolecular interactions of perhalogenated closo -dodecaborate dianions B 12 X 12 2– (X = F, Cl, Br and I) binding to α-, β-, and γ-cyclodextrins (CDs) in the gas phase using combined negative ion photoelectron spectroscopy (NIPES) and density functional theory (DFT) calculations. The vertical detachment energy (VDE) of each complex and electronic stabilization of each dianion due to the CD binding (ΔVDE, relative to the corresponding isolated B 12 X 12 2– ) are determined from the experiments along α-, β- and γ-CD in the form of VDE (ΔVDE): 4.00 (2.10), 4.33 (2.43), and 4.30 (2.40) eV in X = F; 4.09 (1.14), 4.64 (1.69), and 4.69 (1.74) eV in X = Cl; 4.11 (0.91), 4.58 (1.38), and 4.70 (1.50) eV in X = Br; and 3.54 (0.74), 3.88 (1.08), and 4.05 (1.25) eV in X = I, respectively. All complexes have significantly higher VDEs than the corresponding isolated dodecaborate dianions with ΔVDE spanning from 0.74 eV at (α, I) to 2.43 eV at (β, F), sensitive to both host CD size and guest substituent X. DFT-optimized complex structures indicate that all B 12 X 12 2– prefer binding to the wide openings of CDs with the insertion depth and binding motif strongly dependent on the CD size and halogen X. Dodecaborate anions with heavy halogens, i.e. , X = Cl, Br, and I, are found outside of α-CD, while B 12 F 12 2– is completely wrapped by γ-CD. Partial embedment of B 12 X 12 2– into CDs is observed for the other complexes via multipronged B–X … H–O/C interlocking patterns. The simulated spectra based on the density of states agree well with those of the experiments and the calculated VDEs well reproduce the experimental trends. Molecular orbital analyses suggest that the spectral features at low binding energies originated from electrons detached from the dodecaborate dianion, while those at higher binding energies are derived from electron detachment from CDs. Energy decomposition analyses reveal that the electrostatic interaction plays a dominating role in contributing to the host–guest interactions for the X = F series partially due to the formation of a O/C–H … X–B hydrogen bonding network, and the dispersion forces gradually become important with the increase of halogen size.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Combinatorial transcription factor binding encodes cis -regulatory wiring of mouse forebrain GABAergic neurogenesis

Transcription factors (TFs) bind combinatorially to cis-regulatory elements, orchestrating transcriptional programs. Although studies of chromatin state and chromosomal interactions have demonstrated dynamic neurodevelopmental cis-regulatory landscapes, parallel understanding of TF interactions lags. To elucidate combinatorial TF binding driving mouse basal ganglia development, we integrated chromatin immunoprecipitation sequencing (ChIP-seq) for twelve TFs, H3K4me3-associated enhancer-promoter interactions, chromatin and gene expression data, and functional enhancer assays. We identified sets of putative regulatory elements with shared TF binding (TF-pRE modules) that orchestrate distinct processes of GABAergic neurogenesis and suppress other cell fates. The majority of pREs were bound by one or two TFs; however, a small proportion were extensively bound. These sequences had exceptional evolutionary conservation and motif density, complex chromosomal interactions, and activity as in vivo enhancers. Our results provide insights into the combinatorial TF-pRE interactions that activate and repress expression programs during telencephalon neurogenesis and demonstrate the value of TF binding toward modeling developmental transcriptional wiring.

59 BASIC BIOLOGICAL SCIENCES↗

Theory-Guided Inelastic Neutron Scattering of Crystalline Alkaline Aluminate Salts Bearing Principal Motifs of Solution-State Species

Aluminate salts precipitated from caustic alkaline solutions exhibit a correlation between the anionic speciation and the identity of the alkali cation in the precipitate, with the aluminate ions occurring either in monomeric (Al(OH) 4 – ) or dimeric (Al 2 O(OH) 6 2– ) forms. The origin of this correlation is poorly understood as are the roles that oligomeric aluminate species play in determining the solution structure, prenucleation clusters, and precipitation pathways. Characterization of aluminate solution speciation with vibrational spectroscopy results in spectra that are difficult to interpret because the ions access a diverse and dynamic configurational space. To investigate the Al(OH) 4 – and Al 2 O(OH) 6 2– anions within a well-defined crystal lattice, inelastic neutron scattering (INS) and Raman spectroscopic data were collected and simulated by density functional theory for K 2 [Al 2 O(OH) 6 ], Rb 2 [Al 2 O(OH) 6 ], and Cs[Al(OH) 4 ]·2H 2 O. These structures capture archetypal solution aluminate species: the first two salts contain dimeric Al 2 O(OH) 6 2– anions, while the third contains the monomeric Al(OH) 4 – anion. Here, comparisons were made to the INS and Raman spectra of sodium aluminate solutions frozen in a glassy state. In contrast to solution systems, the crystal lattice of the salts results in well-defined vibrations and associated resolved bands in the INS spectra. The use of a theory-guided analysis of the INS of this solid alkaline aluminate series revealed that differences were related to the nature of the hydrogen-bonding network and showed that INS is a sensitive probe of the degree of completeness and strength of the bond network in hydrogen-bonded materials. Results suggest that the ionic size may explain cation-specific differences in crystallization pathways in alkaline aluminate salts.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Active learning of ternary alloy structures and energies

Abstract Machine learning models with uncertainty quantification have recently emerged as attractive tools to accelerate the navigation of catalyst design spaces in a data-efficient manner. Here, we combine active learning with a dropout graph convolutional network (dGCN) as a surrogate model to explore the complex materials space of high-entropy alloys (HEAs). We train the dGCN on the formation energies of disordered binary alloy structures in the Pd-Pt-Sn ternary alloy system and improve predictions on ternary structures by performing reduced optimization of the formation free energy, the target property that determines HEA stability, over ensembles of ternary structures constructed based on two coordinate systems: (a) a physics-informed ternary composition space, and (b) data-driven coordinates discovered by the Diffusion Maps manifold learning scheme. Both reduced optimization techniques improve predictions of the formation free energy in the ternary alloy space with a significantly reduced number of DFT calculations compared to a high-fidelity model. The physics-based scheme converges to the target property in a manner akin to a depth-first strategy, whereas the data-driven scheme appears more akin to a breadth-first approach. Both sampling schemes, coupled with our acquisition function, successfully exploit a database of DFT-calculated binary alloy structures and energies, augmented with a relatively small number of ternary alloy calculations, to identify stable ternary HEA compositions and structures. This generalized framework can be extended to incorporate more complex bulk and surface structural motifs, and the results demonstrate that significant dimensionality reduction is possible in thermodynamic sampling problems when suitable active learning schemes are employed.

Chemistry↗

Porous Semiconducting K–Sn–Mo–S Aerogel: Synthesis, Local Structure, and Ion-Exchange Properties

Chalcogenide-based aerogels are emerging porous semiconducting nanomaterials that appeal to applications in clean energy and the environment. Here, we report a novel gel, potassium–tin–molybdenum–sulfides (KTMS), that integrates the electrostatically bound K + ions in the covalent network of Sn–Mo–S. Its gelation requires a concurrent reduction of Mo 6+ → Mo 4+/5+ and the oxidation of S 2– → Sn – (n ≈ 1) and Sn 2+ → Sn 4+ . KTMS is an amorphous semiconductor showing quantum confinement effects on band gap energies, 2.1 → 1.4 → 0.9 eV for its wet- → aero- → xerogels. Synchrotron X-ray pair distribution function (PDF) and extended X-ray absorption fine structure (EXAFS) revealed a complex local structure of KTMS consisting of molecular Mo 2 (S 2 ) 6 and Mo 3 S(S 2 ) 6 clusters. In addition, the Sn–S coordination is related to crystalline Na4Sn3S8 and SnS2. KTMS also demonstrated the removal of the radionuclides of Cs + , Sr 2+ , and UO 2 2+ from ppm to ppb levels with distribution constants (Kd) up to ≥104 mL/g. Notably, despite the lack of atomic periodicity in the amorphous KTMS, the K+ ion is ion-exchangeable with chemically diverse Sr 2+ , Cs + , and UO 2 2+ in aqueous solutions; especially the ion-exchange properties of Sr 2+ and UO 2 2+ ≡(O=U=O) 2+ is not known to any chalcogels known to date. The sequestration of Cs + and Sr 2+ was achieved by the exchange of K + in the amorphous KTMS, and the removal of [O=U 6+ =O] 2+ synergistically involves surface sorption via -S····U 6+ =O 2 2+ covalent interactions and ion-exchange via the hard–soft Lewis acid–base paradigm. Overall, cooperative roles played by the diverse bonding motifs, surface-exposed Lewis basic frameworks, and polarizability of the (poly)sulfides make it an exceptional adsorbent for chemically diverse radioactive species. This finding will guide the design of superior sorbents for chemically distinct metal ion separation.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Descriptors of water aggregation

For this work, we rely on a total of 23 (cluster size, 8 structural, and 14 connectivity) descriptors to investigate structural patterns and connectivity motifs associated with water cluster aggregation. In addition to the cluster size n (number of molecules), the 8 structural descriptors can be further categorized into (i) one-body (intramolecular): covalent OH bond length (r OH ) and HOH bond angle (θ HOH ), (ii) two-body: OO distance (r OO ), OHO angle (θ OHO ), and HOOX dihedral angle ($\phi$ HOOX ), where X lies on the bisector of the HOH angle, (iii) three-body: OOO angle (θ OOO ), and (iv) many-body: modified tetrahedral order parameter (q) to account for two-, three-, four-, five-coordinated molecules (q m , m = 2, 3, 4, 5) and radius of gyration (R g ). The 14 connectivity descriptors are all many-body in nature and consist of the AD, AAD, ADD, AADD, AAAD, AAADD adjacencies [number of hydrogen bonds accepted (A) and donated (D) by each water molecule], Wiener index, Average Shortest Path Length, hydrogen bond saturation (% HB), and number of non-short-circuited three-membered cycles, four-membered cycles, five-membered cycles, six-membered cycles, and seven-membered cycles. We mined a previously reported database of 4 948 959 water cluster minima for (H 2 O) n , n = 3–25 to analyze the evolution and correlation of these descriptors for the clusters within 5 kcal/mol of the putative minima. It was found that r OH and % HB correlated strongly with cluster size n, which was identified as the strongest predictor of energetic stability. Marked changes in the adjacencies and cycle count were observed, lending insight into changes in the hydrogen bond network upon aggregation. A Principal Component Analysis (PCA) was employed to identify descriptor dependencies and group clusters into specific structural patterns across different cluster sizes. The results of this study inform our understanding of how water clusters evolve in size and what appropriate descriptors of their structural and connectivity patterns are with respect to system size, stability, and similarity. The approach described in this study is general and can be easily extended to other hydrogen-bonded systems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Small Signaling Peptides in Sorghum bicolor : Integrating Phylogeny and Gene Expression to Characterize Roles in Stem Development

Small signaling peptides (SSPs) are important regulators of plant growth, development, and responses to biotic and abiotic stress, yet their role in the C4 grass Sorghum bicolor is largely uncharacterized. To help fill this knowledge gap, 219 sorghum genes that encode SSPs were identified based on SSP sequences previously identified in Arabidopsis thaliana, Zea mays, Oryza sativa, Triticum aestivum , and Brachypodium distachyon . The 219 sorghum SSP-encoding genes were assigned to 19 gene families, analyzed for the presence of motifs, and aligned with genes that encode SSPs in other plants using phylogenetic analysis. Sorghum genes in 12 of the 19 SSP gene families had not been previously characterized. Expression of the 219 SSP-encoding genes in sorghum organs, during stem development, and in stem tissues and cell types revealed distinct spatial, temporal, and developmental patterns of expression. Genes associated with the SbCEP and SbRGF families were preferentially expressed in roots, whereas SbEPF genes were expressed in stem epidermal and pith parenchyma cells and panicles. The expression of genes during bioenergy sorghum stem growth and development was investigated because stems account for ~80% of harvested biomass and serve as conduits for water and nutrient transport between leaves and roots. During stem development, 28 SSP genes in several families ( CLE, EPF, CEP, GASS, PSY, ES, PSK, CAPE, POE ) were expressed at higher levels in zones of cell proliferation. For example, the TDIF homologs SbCLE41 and SbCLE42 were expressed at high levels in nascent stem nodes where they may regulate vascular bundle cambial activity and cell differentiation. A different set of 15 genes in the CIF, POE, CAPE, PSY, CEP, RALF , and CLE families were expressed at higher levels in zones of stem tissue differentiation highlighted by elevated expression of five SbRALFR s in the stem nodal plexus. Cell type–specific expression of many sorghum genes that encode SSPs was observed in fully elongated internodes indicating gene expression is regulated with high spatial resolution. Overall, the results provide a foundation of information for analysis of SSP function in sorghum that can be integrated with knowledge of sorghum gene regulatory networks to modulate traits important for production of sorghum crops.

bioenergy sorghum↗

Ca X ML: Chemistry‐informed machine learning explains mutual changes between protein conformations and calcium ions in calcium‐binding proteins using structural and topological features

Proteins' flexibility is a feature in communicating changes in cell signaling instigated by binding with secondary messengers, such as calcium ions, associated with the coordination of muscle contraction, neurotransmitter release, and gene expression. When binding with the disordered parts of a protein, calcium ions must balance their charge states with the shape of calcium-binding proteins and their versatile pool of partners depending on the circumstances they transmit. Accurately determining the ionic charges of those ions is essential for understanding their role in such processes. However, it is unclear whether the limited experimental data available can be effectively used to train models to accurately predict the charges of calcium-binding protein variants. Here, we developed a chemistry-informed, machine-learning algorithm that implements a game theoretic approach to explain the output of a machine-learning model without the prerequisite of an excessively large database for high-performance prediction of atomic charges. We used the ab initio electronic structure data representing calcium ions and the structures of the disordered segments of calcium-binding peptides with surrounding water molecules to train several explainable models. Network theory was used to extract the topological features of atomic interactions in the structurally complex data dictated by the coordination chemistry of a calcium ion, a potent indicator of its charge state in protein. Our design created a computational tool of Ca X ML, which provided a framework of explainable machine learning model to annotate ionic charges of calcium ions in calcium-binding proteins in response to the chemical changes in an environment. Our framework will provide new insights into protein design for engineering functionality based on the limited size of scientific data in a genome space.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Mechanisms of Ethylene Epoxidation over Silver from Machine Learning-Accelerated First-Principles Modeling and Microkinetic Simulations

In this work we employed machine learning-augmented density functional theory (DFT) thermodynamic calculations to assess the stability of different AgO x structures under catalytic ethylene epoxidation reaction conditions. We found that there are multiple AgO x surface motifs that could co-exist under the relevant conditions. These included Ag surface oxides (e.g., AgO_p(4 × 4) and Ag 1.83 O) and atomic oxygen-covered Ag(111) surfaces. Furthermore, we employed DFT calculations to evaluate the energetics of different reaction mechanisms by which ethylene and oxygen can react on these surfaces. These studies revealed several energetically viable reaction pathways for ethylene epoxidation. Microkinetic modeling analyses, based on the DFT-calculated reaction pathways, showed that ethylene epoxidation can proceed on all surfaces and that multiple pathways, including those involving Langmuir–Hinshelwood and Eley–Rideal mechanisms, could be involved in selective and unselective reactions. The diversity of mechanisms that we discovered in the context of the relatively simple ethylene epoxidation reaction on Ag suggests that the richness and complexity of surface chemistry are most likely a rule rather than an exception in heterogeneous catalytic chemical transformations on metal surfaces and that the concept of a single or even a dominant mechanism and reaction intermediates might need to be revisited for many reactions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Atlas of Transcription Factor Binding Sites from ENCODE DNase Hypersensitivity Data across 27 Tissue Types

Characterizing the tissue-specific binding sites of transcription factors (TFs) is essential to reconstruct gene regulatory networks and predict functions for non-coding genetic variation. DNase-seq footprinting enables the prediction of genome-wide binding sites for hundreds of TFs simultaneously. Despite the public availability of high-quality DNase-seq data from hundreds of samples, a comprehensive, up-to-date resource for the locations of genomic footprints is lacking. Here, we develop a scalable footprinting workflow using two state-of-the-art algorithms: Wellington and HINT. We apply our workflow to detect footprints in 192 ENCODE DNase-seq experiments and predict the genomic occupancy of 1,515 human TFs in 27 human tissues. We validate that these footprints overlap true-positive TF binding sites from ChIP-seq. We demonstrate that the locations, depth, and tissue specificity of footprints predict effects of genetic variants on gene expression and capture a substantial proportion of genetic risk for complex traits.

59 BASIC BIOLOGICAL SCIENCES↗