Arrangement Of Binding Sites Research Articles

End Group Effects on Anion Binding in Tetraglycine Peptide: A Computational Study.

The importance of anions in various processes has led to a search for molecules that can effectively recognize and interact with these anions. This study explores how the tetraglycine (Gly)4 peptide in its zwitterionic, neutral, and terminally capped forms acts as a receptor for H2PO4- and HSO4- anions within the framework of supramolecular host-guest chemistry. Using molecular dynamics (MD) simulations, we obtained the conformations of the receptor-anion complexes. Density functional theory (DFT), quantifies the complexes' interaction energies in both gas and solvent phases. Proton transfer within the zwitterionic complex with H2PO4- anion alters peptide charge distribution, affecting its conformation and binding site arrangement, as analysed by quantum mechanics/molecular mechanics (QM/MM) methods. Symmetry-adapted perturbation theory (SAPT) and noncovalent interactions analysis highlight the role of electrostatic interactions in these receptor-anion complexes. It emphasizes the key interactions such as N-H···O and O-H···O=C between the peptide backbone and anions and elucidates the molecular recognition mechanism driven by crucial noncovalent interactions. The termination of the peptide's end groups modulates anion binding sites from the backbone to the charged N-terminal, resulting in distinct binding sites. Our findings provide insights for designing peptides tailored to function as anion receptors in diverse supramolecular chemistry applications.

Generating information-dense promoter sequences with optimal string packing.

Dense arrangements of binding sites within nucleotide sequences can collectively influence downstream transcription rates or initiate biomolecular interactions. For example, natural promoter regions can harbor many overlapping transcription factor binding sites that influence the rate of transcription initiation. Despite the prevalence of overlapping binding sites in nature, rapid design of nucleotide sequences with many overlapping sites remains a challenge. Here, we show that this is an NP-hard problem, coined here as the nucleotide String Packing Problem (SPP). We then introduce a computational technique that efficiently assembles sets of DNA-protein binding sites into dense, contiguous stretches of double-stranded DNA. For the efficient design of nucleotide sequences spanning hundreds of base pairs, we reduce the SPP to an Orienteering Problem with integer distances, and then leverage modern integer linear programming solvers. Our method optimally packs sets of 20-100 binding sites into dense nucleotide arrays of 50-300 base pairs in 0.05-10 seconds. Unlike approximation algorithms or meta-heuristics, our approach finds provably optimal solutions. We demonstrate how our method can generate large sets of diverse sequences suitable for library generation, where the frequency of binding site usage across the returned sequences can be controlled by modulating the objective function. As an example, we then show how adding additional constraints, like the inclusion of sequence elements with fixed positions, allows for the design of bacterial promoters. The nucleotide string packing approach we present can accelerate the design of sequences with complex DNA-protein interactions. When used in combination with synthesis and high-throughput screening, this design strategy could help interrogate how complex binding site arrangements impact either gene expression or biomolecular mechanisms in varied cellular contexts.

Targeted mutagenesis of specific genomic DNA sequences in animals for the in vivo generation of variant libraries.

Understanding how the number, placement and affinity of transcription factor binding sites dictates gene regulatory programs remains a major unsolved challenge in biology, particularly in the context of multicellular organisms. To uncover these rules, it is first necessary to find the binding sites within a regulatory region with high precision, and then to systematically modulate this binding site arrangement while simultaneously measuring the effect of this modulation on output gene expression. Massively parallel reporter assays (MPRAs), where the gene expression stemming from 10,000s of in vitro-generated regulatory sequences is measured, have made this feat possible in high-throughput in single cells in culture. However, because of lack of technologies to incorporate DNA libraries, MPRAs are limited in whole organisms. To enable MPRAs in multicellular organisms, we generated tools to create a high degree of mutagenesis in specific genomic loci in vivo using base editing. Targeting GFP integrated in genome of Drosophila cell culture and whole animals as a case study, we show that the base editor AIDevoCDA1 stemming from sea lamprey fused to nCas9 is highly mutagenic. Surprisingly, longer gRNAs increase mutation efficiency and expand the mutating window, which can allow the introduction of mutations in previously untargetable sequences. Finally, we demonstrate arrays of >20 gRNAs that can efficiently introduce mutations along a 200bp sequence, making it a promising tool to test enhancer function in vivo in a high throughput manner.

TFSyntax: a database of transcription factors binding syntax in mammalian genomes.

In mammals, transcriptional factors (TFs) drive gene expression by binding to regulatory elements in a cooperative manner. Deciphering the rules of such cooperation is crucial to obtain a full understanding of cellular homeostasis and development. Although this is a long-standing topic, there is no comprehensive database for biologists to access the syntax of TF binding sites. Here we present TFSyntax (https://tfsyntax.zhaopage.com), a database focusing on the arrangement of TF binding sites. TFSyntax maps the binding motif of 1299 human TFs and 890 mouse TFs across 382 cells and tissues, representing the most comprehensive TF binding map to date. In addition to location, TFSyntax defines motif positional preference, density and colocalization within accessible elements. Powered by a series of functional modules based on web interface, users can freely search, browse, analyze, and download data of interest. With comprehensive characterization of TF binding syntax across distinct tissues and cell types, TFSyntax represents a valuable resource and platform for studying the mechanism of transcriptional regulation and exploring how regulatory DNA variants cause disease.

The relative binding position of Nck and Grb2 adaptors impacts actin-based motility of Vaccinia virus.

Phosphotyrosine (pTyr) motifs in unstructured polypeptides orchestrate important cellular processes by engaging SH2-containing adaptors to assemble complex signalling networks. The concept of phase separation has recently changed our appreciation of multivalent networks, however, the role of pTyr motif positioning in their function remains to be explored. We have now investigated this parameter in the operation of the signalling cascade driving actin-based motility and spread of Vaccinia virus. This network involves two pTyr motifs in the viral protein A36 that recruit the adaptors Nck and Grb2 upstream of N-WASP and Arp2/3 complex-mediated actin polymerisation. Manipulating the position of pTyr motifs in A36 and the unrelated p14 from Orthoreovirus, we find that only specific spatial arrangements of Nck and Grb2 binding sites result in robust N-WASP recruitment, Arp2/3 complex driven actin polymerisation and viral spread. This suggests that the relative position of pTyr adaptor binding sites is optimised for signal output. This finding may explain why the relative positions of pTyr motifs are frequently conserved in proteins from widely different species. It also has important implications for regulation of physiological networks, including those undergoing phase transitions.

AdpA Positively Regulates Morphological Differentiation and Chloramphenicol Biosynthesis in Streptomyces venezuelae.

ABSTRACTIn members of genus Streptomyces, AdpA is a master transcriptional regulator that controls the expression of hundreds of genes involved in morphological differentiation, secondary metabolite biosynthesis, chromosome replication, etc. However, the function of AdpASv, an AdpA ortholog of Streptomyces venezuelae, is unknown. This bacterial species is a natural producer of chloramphenicol and has recently become a model organism for studies on Streptomyces. Here, we demonstrate that AdpASv is essential for differentiation and antibiotic biosynthesis in S. venezuelae and provide evidence suggesting that AdpASv positively regulates its own gene expression. We speculate that the different modes of AdpA-dependent transcriptional autoregulation observed in S. venezuelae and other Streptomyces species reflect the arrangement of AdpA binding sites in relation to the transcription start site. Lastly, we present preliminary data suggesting that AdpA may undergo a proteolytic processing and we speculate that this may potentially constitute a novel regulatory mechanism controlling cellular abundance of AdpA in Streptomyces.IMPORTANCE Streptomyces are well-known producers of valuable secondary metabolites which include a large variety of antibiotics and important model organisms for developmental studies in multicellular bacteria. The conserved transcriptional regulator AdpA of Streptomyces exerts a pleiotropic effect on cellular processes, including the morphological differentiation and biosynthesis of secondary metabolites. Despite extensive studies, the function of AdpA in these processes remains elusive. This work provides insights into the role of a yet unstudied AdpA ortholog of Streptomyces venezuelae, now considered a novel model organism. We found that AdpA plays essential role in morphological differentiation and biosynthesis of chloramphenicol, a broad-spectrum antibiotic. We also propose that AdpA may undergo a proteolytic processing that presumably constitutes a novel mechanism regulating cellular abundance of this master regulator.

Characterizing microRNA-mediated modulation of gene expression noise and its effect on synthetic gene circuits.

MicroRNAs (miRNAs) have been shown to modulate gene expression noise, but less is known about how miRNAs with different properties may regulate noise differently. Here, we investigate the role of competing RNAs and the composition of miRNA response elements (MREs) in modulating noise. We find that weak competing RNAs could introduce lower noise than strong competing RNAs. In comparison with a single MRE, both repetitive and composite MREs can reduce the noise at low expression, but repetitive MREs can elevate the noise remarkably at high expression. We further observed the behavior of a synthetic cell-type classifier with miRNAs as inputs and find that miRNAs and MREs that could introduce higher noise tend to enhance cell state transition. These results provide a systematic and quantitative understanding of the function of miRNAs in controlling gene expression noise and the utilization of miRNAs to modulate the behavior of synthetic gene circuits.

DNA sequence-directed cooperation between nucleoid-associated proteins

SummaryNucleoid-associated proteins (NAPs) are a class of highly abundant DNA-binding proteins in bacteria and archaea. While both the composition and relative abundance of the NAPs change during the bacterial growth cycle, surprisingly little is known about their crosstalk in mutually binding and stabilizing higher-order nucleoprotein complexes in the bacterial chromosome. Here, we use atomic force microscopy and solid-state nanopores to investigate long-range nucleoprotein structures formed by the binding of two major NAPs, FIS and H-NS, to DNA molecules with distinct binding site arrangements. We find that spatial organization of the protein binding sites can govern the higher-order architecture of the nucleoprotein complexes. Based on sequence arrangement the complexes differed in their global shape and compaction as well as the extent of FIS and H-NS binding. Our observations highlight the important role the DNA sequence plays in driving structural differentiation within the bacterial chromosome.

Focus on Human Monoamine Transporter Selectivity. New Human DAT and NET Models, Experimental Validation, and SERT Affinity Exploration.

The most commonly used antidepressant drugs are the serotonin transporter inhibitors. Their effects depend strongly on the selectivity for a single monoamine transporter compared to other amine transporters or receptors, and the selectivity is roughly influenced by the spatial protein structure. Here, we provide a computational study on three human monoamine transporters, i.e., DAT, NET, and SERT. Starting from the construction of hDAT and hNET models, whose three-dimensional structure is unknown, and the prediction of the binding pose for 19 known inhibitors, 3D-QSAR models of three human transporters were built. The training set variability, which was high in structure and activity profile, was validated using a set of in-house compounds. Results concern more than one aspect. First of all, hDAT and hNET three-dimensional structures were built, validated, and compared to the hSERT one; second, the computational study highlighted the differences in binding site arrangement statistically correlated to inhibitor selectivity; third, the profiling of new inhibitors pointed out a conservation of the inhibitory activity trend between rabbit and human SERT with a difference of about 1 order of magnitude; fourth, binding and functional studies confirmed 4-(benzyloxy)-4-phenylpiperidine 20a–d and 21a–d as potent SERT inhibitors. In particular, one of the compounds (compound 20b) revealed a higher affinity for SERT than paroxetine in human platelets.

Regulation of spatiotemporal limits of developmental gene expression via enhancer grammar

The regulatory specificity of a gene is determined by the structure of its enhancers, which contain multiple transcription factor binding sites. A unique combination of transcription factor binding sites in an enhancer determines the boundary of target gene expression, and their disruption often leads to developmental defects. Despite extensive characterization of binding motifs in an enhancer, it is still unclear how each binding site contributes to overall transcriptional activity. Using live imaging, quantitative analysis, and mathematical modeling, we measured the contribution of individual binding sites in transcriptional regulation. We show that binding site arrangement within the Rho-GTPase component t48 enhancer mediates the expression boundary by mainly regulating the timing of transcriptional activation along the dorsoventral axis of Drosophila embryos. By tuning the binding affinity of the Dorsal (Dl) and Zelda (Zld) sites, we show that single site modulations are sufficient to induce significant changes in transcription. Yet, no one site seems to have a dominant role; rather, multiple sites synergistically drive increases in transcriptional activity. Interestingly, Dl and Zld demonstrate distinct roles in transcriptional regulation. Dl site modulations change spatial boundaries of t48, mostly by affecting the timing of activation and bursting frequency rather than transcriptional amplitude or bursting duration. However, modulating the binding site for the pioneer factor Zld affects both the timing of activation and amplitude, suggesting that Zld may potentiate higher Dl recruitment to target DNAs. We propose that such fine-tuning of dynamic gene control via enhancer structure may play an important role in ensuring normal development.

Drug ReposER: a web server for predicting similar amino acid arrangements to known drug binding interfaces for potential drug repositioning.

A common drug repositioning strategy is the re-application of an existing drug to address alternative targets. A crucial aspect to enable such repurposing is that the drug's binding site on the original target is similar to that on the alternative target. Based on the assumption that proteins with similar binding sites may bind to similar drugs, the 3D substructure similarity data can be used to identify similar sites in other proteins that are not known targets. The Drug ReposER (DRug REPOSitioning Exploration Resource) web server is designed to identify potential targets for drug repurposing based on sub-structural similarity to the binding interfaces of known drug binding sites. The application has pre-computed amino acid arrangements from protein structures in the Protein Data Bank that are similar to the 3D arrangements of known drug binding sites thus allowing users to explore them as alternative targets. Users can annotate new structures for sites that are similarly arranged to the residues found in known drug binding interfaces. The search results are presented as mappings of matched sidechain superpositions. The results of the searches can be visualized using an integrated NGL viewer. The Drug ReposER server has no access restrictions and is available at http://mfrlab.org/drugreposer/.

Capturing Auxin Response Factors Syntax Using DNA Binding Models

Auxin is a key hormone performing a wealth of functions throughout the life cycle of plants. It acts largely by regulating genes at the transcriptional level through a family of transcription factors called auxin response factors (ARFs). Even though all ARF monomers analyzed so far bind a similar DNA sequence, there is evidence that ARFs differ in their target genomic regions and regulated genes. Here, we report the use of position weight matrices (PWMs) to model ARF DNA binding specificity based on published DNA affinity purification sequencing (DAP-seq) data. We found that the genome binding of two ARFs (ARF2 and ARF5/Monopteros [MP]) differ largely because these two factors have different preferred ARF binding site (ARFbs) arrangements (orientation and spacing). We illustrated why PWMs are more versatile to reliably identify ARFbs than the widely used consensus sequences and demonstrated their power with biochemical experiments in the identification of the regulatory regions of IAA19, an well-characterized auxin-responsive gene. Finally, we combined gene regulation by auxin with ARF-bound regions and identified specific ARFbs configurations that are over-represented in auxin-upregulated genes, thus deciphering the ARFbs syntax functional for regulation. Our study provides a general method to exploit the potential of genome-wide DNA binding assays and to decode gene regulation.

Characterization of an RNase with two catalytic centers. Human RNase6 catalytic and phosphate-binding site arrangement favors the endonuclease cleavage of polymeric substrates

BackgroundHuman RNase6 is a small cationic antimicrobial protein that belongs to the vertebrate RNaseA superfamily. All members share a common catalytic mechanism, which involves a conserved catalytic triad, constituted by two histidines and a lysine (His15/His122/Lys38 in RNase6 corresponding to His12/His119/Lys41 in RNaseA). Recently, our first crystal structure of human RNase6 identified an additional His pair (His36/His39) and suggested the presence of a secondary active site. MethodsIn this work we have explored RNase6 and RNaseA subsite architecture by X-ray crystallography, site-directed mutagenesis and kinetic characterization. ResultsThe analysis of two novel crystal structures of RNase6 in complex with phosphate anions at atomic resolution locates a total of nine binding sites and reveals the contribution of Lys87 to phosphate-binding at the secondary active center. Contribution of the second catalytic triad residues to the enzyme activity is confirmed by mutagenesis. RNase6 catalytic site architecture has been compared with an RNaseA engineered variant where a phosphate-binding subsite is converted into a secondary catalytic center (RNaseA-K7H/R10H). ConclusionsWe have identified the residues that participate in RNase6 second catalytic triad (His36/His39/Lys87) and secondary phosphate-binding sites. To note, residues His39 and Lys87 are unique within higher primates. The RNaseA/RNase6 side-by-side comparison correlates the presence of a dual active site in RNase6 with a favored endonuclease-type cleavage pattern. General significanceAn RNase dual catalytic and extended binding site arrangement facilitates the cleavage of polymeric substrates. This is the first report of the presence of two catalytic centers in a single monomer within the RNaseA superfamily.

MBNL splicing activity depends on RNA binding site structural context.

Muscleblind-like (MBNL) proteins are conserved RNA-binding factors involved in alternative splicing (AS) regulation during development. While AS is controlled by distribution of MBNL paralogs and isoforms, the affinity of these proteins for specific RNA-binding regions and their location within transcripts, it is currently unclear how RNA structure impacts MBNL-mediated AS regulation. Here, we defined the RNA structural determinants affecting MBNL-dependent AS activity using both cellular and biochemical assays. While enhanced inclusion of MBNL-regulated alternative exons is controlled by the arrangement and number of MBNL binding sites within unstructured RNA, when these sites are embedded in a RNA hairpin MBNL binds preferentially to one side of stem region. Surprisingly, binding of MBNL proteins to RNA targets did not entirely correlate with AS efficiency. Moreover, comparison of MBNL proteins revealed structure-dependent competitive behavior between the paralogs. Our results showed that the structure of targeted RNAs is a prevalent component of the mechanism of alternative splicing regulation by MBNLs.

Sialyllactose-Modified Three-Way Junction DNA as Binding Inhibitor of Influenza Virus Hemagglutinin

Sialic acid present on the cell surface is recognized by hemagglutinin (HA) on the influenza virus in the first step of infection. Therefore, a compound that can efficiently interfere with the interaction between sialic acid and HA might inhibit infection and allow detection of the influenza virus. We focused on the spatial arrangement of sialic acid binding sites on HA and developed 2,3-sialyllactose (2,3-SL)-modified three-way junction (3WJ) DNA molecules with a topology similar to that of sialic acid binding sites. 3WJ DNA with three 2,3-SL residues on each DNA strand showed (8.0 × 104)-fold higher binding affinity for influenza virus A/Puerto Rico/08/34 (H1N1) compared to the 2,3-SL. This result indicated that the glycocluster effect due to clustering on one DNA arm and optimal spatial arrangement of the 3WJ DNA improved the weak interactions between a sialic acid and its binding site on HA. This 3WJ DNA compound has possible application as an inhibitor of influenza infection and for virus sensing.

Binding maps for the study and prediction of bimetallic catalyst surface reactions: The case of methanol oxidation

AbstractFocusing on the competing pathways of methanol oxidation on platinum and platinum/gold bimetallic catalysts, we explore a novel density functional theory (DFT)‐based approach to the study of reactions on catalyst surfaces. Traditionally, DFT has been used to compute binding energies of products and intermediates as proxies for catalytic activity, and to compute full reaction pathways and their activation energy barriers. Merging the computational simplicity and intuitive clarity of binding energy calculations with the site sensitivity of transition state calculations, we construct maps of the binding energies of relevant atoms and molecules at all sites on a surface. We show that knowledge of the arrangement of strong and weak binding sites on a surface is powerful in rationalizing the ease with which a reaction step proceeds on a given local motif of surface atoms. We highlight the prospects and challenges of this approach toward catalyst screening and prediction.

Abstract 5232: The molecular basis for specific recognition of the biologically relevant hybrid-2 type human telomeric G-quadruplex by epiberberine

Abstract G-quadruplex structures have been shown to form in human telomeres, and the formation of G-quadruplexes can inhibit telomerase, which plays a key role in cancers by stabilizing telomere length and integrity thereby granting limitless replicative potential. The human telomeric G-quadruplex is thus considered to be a potential target for cancer therapeutics. In physiologically relevant potassium solution, human telomeric DNA sequences form two equilibrating hybrid-type G-quadruplex structures, with the hybrid-2 structure being the predominant form in extended sequences. We discovered that epiberberine (EPI), a naturally occurring protoberberine alkaloid, can specifically bind the hybrid-2 human telomeric G-quadruplex and can induce the conversion to the hybrid-2 structure. Using nuclear magnetic resonance (NMR) spectroscopy, we determined the solution structure of the 1:1 complex of EPI and hybrid-2 human telomeric G-quadruplex. Our NMR structure shows EPI bound at the 5’ end of the hybrid-2 telomeric quadruplex with an unexpectedly large drug-induced conformational change in the flanking and loop regions, creating a very well-defined drug binding pocket with extensive capping structures. The EPI molecule and 5’ flanking adenine form an induced quasi-triad plane which is intercalated between the external 5’ tetrad and capping structures. Two layers of new capping structures are formed with the 5’ flanking segment and the second TTA lateral loop to cover the “intercalated quasi-triad”. Notably, the well-defined EPI-induced multi-layer binding-site arrangement is only possible in the hybrid-2 folding topology of telomeric DNA. Our results provide important insights into specific targeting of the physiologically relevant hybrid-2 human telomeric G-quadruplex by a small molecule compound; EPI’s asymmetric crescent-shape and the positioning of its dioxolane moiety, as well as the hybrid-2 folding and the loop and flanking bases in the human telomeric DNA sequence, all contribute toward the specific recognition of EPI. In addition, we have conducted polymerase stop assays, which confirmed that EPI can stabilize the human telomeric G-quadruplex to inhibit polymerase activity. Citation Format: Clement Lin, Guanhui Wu, Yong Shao, Danzhou Yang. The molecular basis for specific recognition of the biologically relevant hybrid-2 type human telomeric G-quadruplex by epiberberine [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2017; 2017 Apr 1-5; Washington, DC. Philadelphia (PA): AACR; Cancer Res 2017;77(13 Suppl):Abstract nr 5232. doi:10.1158/1538-7445.AM2017-5232

Model-guided combinatorial optimization of complex synthetic gene networks.

Constructing gene circuits that satisfy quantitative performance criteria has been a long‐standing challenge in synthetic biology. Here, we show a strategy for optimizing a complex three‐gene circuit, a novel proportional miRNA biosensor, using predictive modeling to initiate a search in the phase space of sensor genetic composition. We generate a library of sensor circuits using diverse genetic building blocks in order to access favorable parameter combinations and uncover specific genetic compositions with greatly improved dynamic range. The combination of high‐throughput screening data and the data obtained from detailed mechanistic interrogation of a small number of sensors was used to validate the model. The validated model facilitated further experimentation, including biosensor reprogramming and biosensor integration into larger networks, enabling in principle arbitrary logic with miRNA inputs using normal form circuits. The study reveals how model‐guided generation of genetic diversity followed by screening and model validation can be successfully applied to optimize performance of complex gene networks without extensive prior knowledge.

Transcriptomics-Guided Design of Synthetic Promoters for a Mammalian System.

Despite recent advances in improving titers for therapeutic proteins such as antibodies to the 10 g/L scale, these high yields can only be achieved in select mammalian hosts. Regardless of the host or product, strong promoters are required to obtain high levels of transgene expression. However, the promoters employed to drive this expression are rather limited in variety and are usually either viral-derived or screened empirically during vector design. To begin to move away from viral parts, we employed a more systematic approach to identify and design new synthetic promoters using endogenous elements. To do so, we established a workflow to design these elements by (1) analyzing the transcriptomics profile of a specific cell line under a desired, representative cell culture condition, (2) identifying key genetic motifs using bioinformatics that can be used to rationally construct synthetic promoters, (3) building synthetic promoters using conventional DNA synthesis and molecular biology techniques, and (4) evaluating the performance of these synthetic promoters using model proteins. The resulting promoters perform comparably to the hCMV IE promoter variants tested, but with endogenous components. During this design-build-test cycle we also investigated the underlying design rules for transcription factor binding site arrangement in synthetic promoters. Overall, this approach of using an "omics-guided" workflow for designing synthetic promoters facilitates the construction of high expression vectors for immediate use in current production hosts.

SiteOut: An Online Tool to Design Binding Site-Free DNA Sequences.

DNA-binding proteins control many fundamental biological processes such as transcription, recombination and replication. A major goal is to decipher the role that DNA sequence plays in orchestrating the binding and activity of such regulatory proteins. To address this goal, it is useful to rationally design DNA sequences with desired numbers, affinities and arrangements of protein binding sites. However, removing binding sites from DNA is computationally non-trivial since one risks creating new sites in the process of deleting or moving others. Here we present an online binding site removal tool, SiteOut, that enables users to design arbitrary DNA sequences that entirely lack binding sites for factors of interest. SiteOut can also be used to delete sites from a specific sequence, or to introduce site-free spacers between functional sequences without creating new sites at the junctions. In combination with commercial DNA synthesis services, SiteOut provides a powerful and flexible platform for synthetic projects that interrogate regulatory DNA. Here we describe the algorithm and illustrate the ways in which SiteOut can be used; it is publicly available at https://depace.med.harvard.edu/siteout/.