Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Sequence”

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 361 records · Page 20
All sourcesNASA NTRSDOE OSTI
DOE OSTI2016

Supramolecular polymerization of a prebiotic nucleoside provides insights into the creation of sequence-controlled polymers

The self-assembly of a nucleoside on Au(111) was studied to ascertain whether polymerization on well-defined substrates constitutes a promising approach for making sequence-controlled polymers. Scanning tunneling microscopy and density functional theory were used to investigate the self-assembly on Au(111) of (RS)-N9-(2,3-dihydroxypropyl)adenine (DHPA), a plausibly prebiotic nucleoside analog of adenosine. It is found that DHPA molecules self-assemble into a hydrogen-bonded polymer that grows almost exclusively along the herringbone reconstruction pattern, has a two component sequence that is repeated over hundreds of nanometers, and is erasable with electron-induced excitation. Although the sequence is simple, more complicated ones are envisioned if two or more nucleoside types are combined. Because polymerization occurs on a substrate in a dry environment, the success of each combination can be gauged with high-resolution imaging and accurate modeling techniques. The resulting characteristics make nucleoside self-assembly on a substrate an attractive approach for designing sequence-controlled polymers. Moreover, by choosing plausibly prebiotic nucleosides, insights may be provided into how nature created the first sequence-controlled polymers capable of storing information. Such insights, in turn, can inspire new ways of synthesizing sequence-controlled polymers.

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

Aging of galaxies along the morphological sequence, marked by bulge growth and disk quenching

Aims: We revisit the color bimodality of galaxies using the extensive EFIGI morphological classification of nearby galaxies. Methods: The galaxy profiles from the Sloan Digital Sky Survey (SDSS) gri images were decomposed as a bulge and a disk by controlled profile modeling with the Euclid SourceXtractor++ software. The spectral energy distributions from our resulting gri SDSS photometry complemented with Galaxy Evolution Explorer (GALEX) NUV photometry were fitted with the ZPEG software and PEGASE.2 templates in order to estimate the stellar masses and specific star formation rates (sSFR) of whole galaxies as well as their bulge and disk components. Results: The absolute NUV-r color versus stellar mass diagram shows a continuous relationship between the present sSFR of galaxies and their stellar mass, which spans all morphological types of the Hubble sequence monotonously. Irregular galaxies to intermediate-type Sab spirals make up the “Blue Cloud” across 4 orders of magnitude in stellar mass but a narrow range of sSFR. This mass build-up of spiral galaxies requires major mergers, in agreement with their frequently perturbed isophotes. At high mass, the Blue Cloud leads to the “Green Plain”, dominated by S0a and Sa early-type spirals. It was formerly called the “Green Valley”, due to its low density, but we rename it because of its wide stretch and nearly flat density over ~2 mag in NUV-r color (hence sSFR), despite a limited range of stellar mass (1 order of magnitude). The Green Plain links up the “Red Sequence”, containing all lenticular and elliptical galaxies with a 2 order of magnitude mass interval, and systematically higher masses for the ellipticals. We confirm that the Green Plain cannot be studied using u - r optical colors because it is overlayed by the Red Sequence, hence NUV data are necessary. Galaxies across the Green Plain undergo a marked growth by a factor 2 to 3 in their bulge-to-total mass ratio and a systematic profile change from pseudo to classical bulges, as well as a significant reddening due to star formation fading in their disks. The Green Plain is also characterized by a maximum stellar mass of 10 11.7 M ⊙ beyond which only elliptical galaxies exist, hence supporting the scenario of ellipticals partly forming by major mergers of massive disk galaxies. Conclusions: The EFIGI attributes indicate that dynamical processes (spiral arms and isophote distortions) contribute to the scatter of the Main Sequence of star-forming galaxies (Blue Cloud), via the enhancement of star formation (flocculence, HII regions). The significant bulge growth across the Green Plain confirms that it is a transition region, and excludes a predominantly quick transit due to rapid quenching. The high frequency of bars for all spirals as well as the stronger spiral arms and flocculence in the knee of the Green Plain suggest that internal dynamics, likely triggered by flybys or (mainly minor) mergers, may be the key to the bulge growth of massive disk galaxies, which is a marker of the aging of galaxies from star forming to quiescence. The Hubble sequence can then be considered as an inverse sequence of galaxy physical evolution.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2019

Steric complementarity directs sequence promiscuous leader binding in RiPP biosynthesis

Enzymes that generate ribosomally synthesized and posttranslationally modified peptide (RiPP) natural products have garnered significant interest, given their ability to produce large libraries of chemically diverse scaffolds. Such RiPP biosynthetic enzymes are predicted to bind their corresponding peptide substrates through sequence-specific recognition of the leader sequence, which is removed after the installation of posttranslational modifications on the core sequence. The conservation of the leader sequence within a given RiPP class, in otherwise disparate precursor peptides, further supports the notion that strict sequence specificity is necessary for leader peptide engagement. Here, we demonstrate that leader binding by a biosynthetic enzyme in the lasso peptide class of RiPPs is directed by a minimal number of hydrophobic interactions. Biochemical and structural data illustrate how a single leader-binding domain can engage sequence-divergent leader peptides using a conserved motif that facilitates hydrophobic packing. We report the presence of this simple motif in noncognate peptides results in low micromolar affinity binding by binding domains from several different lasso biosynthetic systems. We also demonstrate that these observations likely extend to other RiPP biosynthetic classes. The portability of the binding motif opens avenues for the engineering of semisynthetic hybrid RiPP products.

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

Neural networks to learn protein sequence–function relationships from deep mutational scanning data

Understanding the relationship between protein sequence and function is necessary to design new and useful proteins with applications in bioenergy, medicine, and agriculture. The mapping from sequence to function is tremendously complex because it involves thousands of molecular interactions that are coupled over multiple lengths and timescales. We show that neural networks can learn the sequence–function mapping from large protein datasets. Neural networks are appealing for this task because they can learn complicated relationships from data, make few assumptions about the nature of the sequence–function relationship, and can learn general rules that apply across the length of the protein sequence. We demonstrate that learned models can be applied to design new proteins with properties that exceed natural sequences.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Prevention efficacy of the broadly neutralizing antibody VRC01 depends on HIV-1 envelope sequence features

In the Antibody Mediated Prevention (AMP) trials (HVTN 704/HPTN 085 and HVTN 703/HPTN 081), prevention efficacy (PE) of the monoclonal broadly neutralizing antibody (bnAb) VRC01 (vs. placebo) against HIV-1 acquisition diagnosis varied according to the HIV-1 Envelope (Env) neutralization sensitivity to VRC01, as measured by 80% inhibitory concentration (IC80). Here, we performed a genotypic sieve analysis, a complementary approach to gaining insight into correlates of protection that assesses how PE varies with HIV-1 sequence features. We analyzed HIV-1 Env amino acid (AA) sequences from the earliest available HIV-1 RNA-positive plasma samples from AMP participants diagnosed with HIV-1 and identified Env sequence features that associated with PE. The strongest Env AA sequence correlate in both trials was VRC01 epitope distance that quantifies the divergence of the VRC01 epitope in an acquired HIV-1 isolate from the VRC01 epitope of reference HIV-1 strains that were most sensitive to VRC01-mediated neutralization. In HVTN 704/HPTN 085, the Env sequence-based predicted probability that VRC01 IC80 against the acquired isolate exceeded 1 µg/mL also significantly associated with PE. In HVTN 703/HPTN 081, a physicochemical-weighted Hamming distance across 50 VRC01 binding-associated Env AA positions of the acquired isolate from the most VRC01-sensitive HIV-1 strain significantly associated with PE. These results suggest that incorporating mutation scoring by BLOSUM62 and weighting by the strength of interactions at AA positions in the epitope:VRC01 interface can optimize performance of an Env sequence-based biomarker of VRC01 prevention efficacy. Future work could determine whether these results extend to other bnAbs and bnAb combinations.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2013

DETECTION OF HIGH-DEGREE PROGRADE SECTORAL MODE SEQUENCES IN THE A-STAR KIC 8054146?

This paper examines the 46 frequencies found in the {delta} Sct star KIC 8054146 involving a frequency spacing of exactly 2.814 cycles day{sup -1} (32.57 {mu}Hz), which is also a dominant low-frequency peak near or equal to the rotational frequency. These 46 frequencies range up to 146 cycles day{sup -1}. Three years of Kepler data reveal distinct sequences of these equidistantly spaced frequencies, including the basic sequence and side lobes associated with other dominant modes (i.e., small amplitude modulations). The amplitudes of the basic sequence show a high-low pattern. The basic sequence follows the equation f{sub m} = 2.8519 + m*2.81421 cycles day{sup -1} with m ranging from 25 to 35. The zero-point offset and the lack of low-order harmonics eliminate an interpretation in terms of a Fourier series of a non-sinusoidal light curve. The exactness of the spacing eliminates high-order asymptotic pulsation. The frequency pattern is not compatible with simple hypotheses involving single or multiple spots, even with differential rotation. The basic high-frequency sequence is interpreted in terms of prograde sectoral modes. These can be marginally unstable, while their corresponding low-degree counterparts are stable due to stronger damping. The measured projected rotation velocity (300 km s{sup -1}) indicates that the star rotates with {approx}>70% of the Keplerian break-up velocity. This suggests a near equator-on view. We qualitatively examine the visibility of prograde sectoral high-degree g-modes in integrated photometric light in such a geometrical configuration and find that prograde sectoral modes can reproduce the frequencies and the odd-even amplitude pattern of the high-frequency sequence.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

wavess 1.2: presenting an HLA-aware within-host virus sequence simulation framework

Motivation Understanding how virus sequences are shaped by selection can inform vaccine design and transmission inference. Modeling within-host evolution to interrogate these questions requires a detailed mechanistic framework that accurately captures sequence diversification. The CD8 + cytotoxic T-lymphocyte (CTL) response plays an important role in immune-mediated selection and can leave strong signatures in virus sequences; however, existing sequence-based within-host virus modeling frameworks do not explicitly include a human leukocyte antigen (HLA)-aware CTL response. Results We extended our previously published within-host sequence evolution simulator, wavess, to include an explicit CTL response, and share a method for identifying HLA-specific CTL epitopes given a founder virus sequence. We also updated the model to permit a variable recombination rate, which allows for modeling non-adjacent genes, segmented genomes, and recombination hotspots. These extensions to wavess allow for more accurate simulation of viruses and virus genes, particularly in regions of the genome where the immune response is dominated by CTLs (rather than antibodies). It also provides the foundation for investigations of how these newly-added biological mechanisms influence within-host evolution. Availability and implementation The core of wavess is written in Python 3, with helper functions written in R. It is available at https://github.com/MolEvolEpid/wavess.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2018

MicroPheno: predicting environments and host phenotypes from 16S rRNA gene sequencing using a k-mer based representation of shallow sub-samples

Microbial communities play important roles in the function and maintenance of various biosystems, ranging from the human body to the environment. A major challenge in microbiome research is the classification of microbial communities of different environments or host phenotypes. The most common and cost-effective approach for such studies to date is 16S rRNA gene sequencing. Recent falls in sequencing costs have increased the demand for simple, efficient and accurate methods for rapid detection or diagnosis with proved applications in medicine, agriculture and forensic science. We describe a reference- and alignment-free approach for predicting environments and host phenotypes from 16S rRNA gene sequencing based on k-mer representations that benefits from a bootstrapping framework for investigating the sufficiency of shallow sub-samples. Deep learning methods as well as classical approaches were explored for predicting environments and host phenotypes. A k-mer distribution of shallow sub-samples outperformed Operational Taxonomic Unit (OTU) features in the tasks of body-site identification and Crohn’s disease prediction. Aside from being more accurate, using k-mer features in shallow sub-samples allows (i) skipping computationally costly sequence alignments required in OTU-picking and (ii) provided a proof of concept for the sufficiency of shallow and short-length 16S rRNA sequencing for phenotype prediction. In addition, k-mer features predicted representative 16S rRNA gene sequences of 18 ecological environments, and 5 organismal environments with high macro-F1 scores of 0.88 and 0.87. For large datasets, deep learning outperformed classical methods such as Random Forest and Support Vector Machine.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Z-Sequence: photometric redshift predictions for galaxy clusters with sequential random k-nearest neighbours

ABSTRACT We introduce Z-Sequence, a novel empirical model that utilizes photometric measurements of observed galaxies within a specified search radius to estimate the photometric redshift of galaxy clusters. Z-Sequence itself is composed of a machine learning ensemble based on the k-nearest neighbours algorithm. We implement an automated feature selection strategy that iteratively determines appropriate combinations of filters and colours to minimize photometric redshift prediction error. We intend for Z-Sequence to be a standalone technique but it can be combined with cluster finders that do not intrinsically predict redshift, such as our own DEEP-CEE. In this proof-of-concept study, we train, fine-tune, and test Z-Sequence on publicly available cluster catalogues derived from the Sloan Digital Sky Survey. We determine the photometric redshift prediction error of Z-Sequence via the median value of |Δ$z$|/(1 + $z$) (across a photometric redshift range of 0.05 ≤ $z$ ≤ 0.6) to be ∼0.01 when applying a small search radius. The photometric redshift prediction error for test samples increases by 30–50 per cent when the search radius is enlarged, likely due to line-of-sight interloping galaxies. Eventually, we aim to apply Z-Sequence to upcoming imaging surveys such as the Legacy Survey of Space and Time to provide photometric redshift estimates for large samples of as yet undiscovered and distant clusters.

Chan, Matthew C.↗
DOE OSTI2014

Impact of sequencing depth in ChIP-seq experiments

In a chromatin immunoprecipitation followed by high-throughput sequencing (ChIP-seq) experiment, an important consideration in experimental design is the minimum number of sequenced reads required to obtain statistically significant results. We present an extensive evaluation of the impact of sequencing depth on identification of enriched regions for key histone modifications (H3K4me3, H3K36me3, H3K27me3 and H3K9me2/me3) using deep-sequenced datasets in human and fly. We propose to define sufficient sequencing depth as the number of reads at which detected enrichment regions increase <1% for an additional million reads. Although the required depth depends on the nature of the mark and the state of the cell in each experiment, we observe that sufficient depth is often reached at <20 million reads for fly. For human, there are no clear saturation points for the examined datasets, but our analysis suggests 40–50 million reads as a practical minimum for most marks. We also devise a mathematical model to estimate the sufficient depth and total genomic coverage of a mark. Lastly, we find that the five algorithms tested do not agree well for broad enrichment profiles, especially at lower depths. Our findings suggest that sufficient sequencing depth and an appropriate peak-calling algorithm are essential for ensuring robustness of conclusions derived from ChIP-seq data.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2016

Mapping DNA polymerase errors by single-molecule sequencing

Genomic integrity is compromised by DNA polymerase replication errors, which occur in a sequence-dependent manner across the genome. Accurate and complete quantification of a DNA polymerase's error spectrum is challenging because errors are rare and difficult to detect. We report a high-throughput sequencing assay to map in vitro DNA replication errors at the single-molecule level. Unlike previous methods, our assay is able to rapidly detect a large number of polymerase errors at base resolution over any template substrate without quantification bias. To overcome the high error rate of high-throughput sequencing, our assay uses a barcoding strategy in which each replication product is tagged with a unique nucleotide sequence before amplification. Here, this allows multiple sequencing reads of the same product to be compared so that sequencing errors can be found and removed. We demonstrate the ability of our assay to characterize the average error rate, error hotspots and lesion bypass fidelity of several DNA polymerases.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2018

RNAcentral: a hub of information for non-coding RNA sequences

RNAcentral is a comprehensive database of non-coding RNA (ncRNA) sequences, collating information on ncRNA sequences of all types from a broad range of organisms. We have recently added a new genome mapping pipeline that identifies genomic locations for ncRNA sequences in 296 species. We have also added several new types of functional annotations, such as tRNA secondary structures, Gene Ontology annotations, and miRNA-target interactions. A new quality control mechanism based on Rfam family assignments identifies potential contamination, incomplete sequences, and more. The RNAcentral database has become a vital component of many workflows in the RNA community, serving as both the primary source of sequence data for academic and commercial groups, as well as a source of stable accessions for the annotation of genomic and functional features. These examples are facilitated by an improved RNAcentral web interface, which features an updated genome browser, a new sequence feature viewer, and improved text search functionality.

Sweeney, Blake A.↗
DOE OSTI2020

Adaptive Anomaly Detection for Dynamic Clinical Event Sequences

Over the past decade, health information technology (IT) has enabled the amount of digital information stored in electronic health records (EHRs) to expand greatly. However, according to some studies, hazards in health IT can lead to changes in clinical decisions, care processes, and care outcomes, as well as other issues. Thus, the effects of health IT hazards on patient safety have been at the forefront of recent patient safety research. Nonetheless, hazard detection in health IT remains a challenge. In this paper, the authors assume that safety-related issues in health IT would exhibit anomalous characteristics in EHR data. Although all hazards will exhibit some anomalous characteristics, not all anomalies can be regarded as hazards. The authors hypothesize that errors in health IT could lead to interruptions in the sequence of clinical actions. To this end, the problem of detecting anomalous sequences in big EHR data is considered. This paper focuses on dynamic event sequences, which are a series of clinical actions in motion. The authors propose an adaptive anomaly detection approach that uses higher-order network representation to detect anomalous sequences. Furthermore, the authors propose a contiguous subsequence anomaly detection approach that identifies abnormal subsequences in the detected anomalous sequences. The proposed approaches are tested by using synthetic and real-world EHR data. The proposed methods outperform existing state of the art anomaly detection techniques. To reduce the computational complexity associated with the operational implementation of the proposed approaches, the Apache Spark environment was leveraged, and a much shorter run time together with improved performance were achieved, especially for data with more than 60,000 sequences.

Niu, Haoran↗
DOE OSTI2012

pGraph: Efficient Parallel Construction of Large-Scale Protein Sequence Homology Graphs

Detecting sequence homology between protein sequences is a fundamental problem in computational molecular biology, with a pervasive application in nearly all analyses that aim to structurally and functionally characterize protein molecules. While detecting the homology between two protein sequences is relatively inexpensive, detecting pairwise homology for a large number of protein sequences can become computationally prohibitive for modern inputs, often requiring millions of CPU hours. Yet, there is currently no robust support to parallelize this kernel. In this paper, we identify the key characteristics that make this problemparticularly hard to parallelize, and then propose a new parallel algorithm that is suited for detecting homology on large data sets using distributed memory parallel computers. Our method, called pGraph, is a novel hybrid between the hierarchical multiple-master/worker model and producer-consumer model, and is designed to break the irregularities imposed by alignment computation and work generation. Experimental results show that pGraph achieves linear scaling on a 2,048 processor distributed memory cluster for a wide range of inputs ranging from as small as 20,000 sequences to 2,560,000 sequences. In addition to demonstrating strong scaling, we present an extensive report on the performance of the various system components and related parametric studies.

Environmental Molecular Sciences Laboratory↗
DOE OSTI2014

TH-C-BRD-07: Minimizing Dose Uncertainty for Spot Scanning Beam Proton Therapy of Moving Tumor with Optimization of Delivery Sequence

Purpose: Intensity modulated proton therapy (IMPT) has been shown to be able to reduce dose to normal tissue compared to intensity modulated photon radio-therapy (IMRT), and has been implemented for selected lung cancer patients. However, respiratory motion-induced dose uncertainty remain one of the major concerns for the radiotherapy of lung cancer, and the utility of IMPT for lung patients was limited because of the proton dose uncertainty induced by motion. Strategies such as repainting and tumor tracking have been proposed and studied but repainting could result in unacceptable long delivery time and tracking is not yet clinically available. We propose a novel delivery strategy for spot scanning proton beam therapy. Method: The effective number of delivery (END) for each spot position in a treatment plan was calculated based on the parameters of the delivery system, including time required for each spot, spot size and energy. The dose uncertainty was then calculated with an analytical formula. The spot delivery sequence was optimized to maximize END and minimize the dose uncertainty. 2D Measurements with a detector array on a 1D moving platform were performed to validate the calculated results. Results: 143 2D measurements on a moving platform were performed for different delivery sequences of a single layer uniform pattern. The measured dose uncertainty is a strong function of the delivery sequence, the worst delivery sequence results in dose error up to 70% while the optimized delivery sequence results in dose error of <5%. END vs. measured dose uncertainty follows the analytical formula. Conclusion: With optimized delivery sequence, it is feasible to minimize the dose uncertainty due to motion in spot scanning proton therapy.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

SU-F-T-350: Continuous Leaf Optimization (CLO) for IMRT Leaf Sequencing

Purpose: To study a new step-and-shoot IMRT leaf sequencing model that avoids the two main pitfalls of conventional leaf sequencing: (1) target fluence being stratified into a fixed number of discrete levels and/or (2) aperture leaf positions being restricted to a discrete set of locations. These assumptions induce error into the sequence or reduce the feasible region of potential plans, respectively. Methods: We develop a one-dimensional (single leaf pair) methodology that does not make assumptions (1) or (2) that can be easily extended to a multi-row model. The proposed continuous leaf optimization (CLO) methodology takes in an existing set of apertures and associated intensities, or solution “seed,” and improves the plan without the restrictiveness of 1or (2). It then uses a first-order descent algorithm to converge onto a locally optimal solution. A seed solution can come from models that assume (1) and (2), thus allowing the CLO model to improve upon existing leaf sequencing methodologies. Results: The CLO model was applied to 208 generated target fluence maps in one dimension. In all cases for all tested sequencing strategies, the CLO model made improvements on the starting seed objective function. The CLO model also was able to keep MUs low. Conclusion: The CLO model can improve upon existing leaf sequencing methods by avoiding the restrictions of (1) and (2). By allowing for more flexible leaf positioning, error can be reduced when matching some target fluence. This study lays the foundation for future models and solution methodologies that can incorporate continuous leaf positions explicitly into the IMRT treatment planning model. Supported by Cancer Prevention & Research Institute of Texas (CPRIT) - ID RP150485.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

TU-H-CAMPUS-JeP3-05: Adaptive Determination of Needle Sequence HDR Prostate Brachytherapy with Divergent Needle-By-Needle Delivery

Purpose: To develop a new method which adaptively determines the optimal needle insertion sequence for HDR prostate brachytherapy involving divergent needle-by-needle dose delivery by e.g. a robotic device. A needle insertion sequence is calculated at the beginning of the intervention and updated after each needle insertion with feedback on needle positioning errors. Methods: Needle positioning errors and anatomy changes may occur during HDR brachytherapy which can lead to errors in the delivered dose. A novel strategy was developed to calculate and update the needle sequence and the dose plan after each needle insertion with feedback on needle positioning errors. The dose plan optimization was performed by numerical simulations. The proposed needle sequence determination optimizes the final dose distribution based on the dose coverage impact of each needle. This impact is predicted stochastically by needle insertion simulations. HDR procedures were simulated with varying number of needle insertions (4 to 12) using 11 patient MR data-sets with PTV, prostate, urethra, bladder and rectum delineated. Needle positioning errors were modeled by random normally distributed angulation errors (standard deviation of 3 mm at the needle’s tip). The final dose parameters were compared in the situations where the needle with the largest vs. the smallest dose coverage impact was selected at each insertion. Results: Over all scenarios, the percentage of clinically acceptable final dose distribution improved when the needle selected had the largest dose coverage impact (91%) compared to the smallest (88%). The differences were larger for few (4 to 6) needle insertions (maximum difference scenario: 79% vs. 60%). The computation time of the needle sequence optimization was below 60s. Conclusion: A new adaptive needle sequence determination for HDR prostate brachytherapy was developed. Coupled to adaptive planning, the selection of the needle with the largest dose coverage impact increases chances of reaching the clinical constraints. M. Borot de Battisti is funded by Philips Medical Systems Nederland B.V.; M. Moerland is principal investigator on a contract funded by Philips Medical Systems Nederland B.V.; G. Hautvast and D. Binnekamp are fulltime employees of Philips Medical Systems Nederland B.V.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2021

Quantification of Cas9 binding and cleavage across diverse guide sequences maps landscapes of target engagement

The RNA-guided nuclease Cas9 has unlocked powerful methods for perturbing both the genome through targeted DNA cleavage and the regulome through targeted DNA binding, but limited biochemical data have hampered efforts to quantitatively model sequence perturbation of target binding and cleavage across diverse guide sequences. We present scalable, sequencing-based platforms for high-throughput filter binding and cleavage and then perform 62,444 quantitative binding and cleavage assays on 35,047 on- and off-target DNA sequences across 90 Cas9 ribonucleoproteins (RNPs) loaded with distinct guide RNAs. We observe that binding and cleavage efficacy, as well as specificity, vary substantially across RNPs; canonically studied guides often have atypically high specificity; sequence context surrounding the target modulates Cas9 on-rate; and Cas9 RNPs may sequester targets in nonproductive states that contribute to “proofreading” capability. Lastly, we distill our findings into an interpretable biophysical model that predicts changes in binding and cleavage for diverse target sequence perturbations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗