A Quantitative and Predictive Model for RNA Binding by Human Pumilio Proteins.

High-throughput methodologies have enabled routine generation of RNA target sets and sequence motifs for RNA-binding proteins (RBPs). Nevertheless, quantitative approaches are needed to capture the landscape of RNA-RBP interactions responsible for cellular regulation. We have used the RNA-MaP platform to directly measure equilibrium binding for thousands of designed RNAs and to construct a predictive model for RNA recognition by the human Pumilio proteins PUM1 and PUM2. Despite prior findings of linear sequence motifs, our measurements revealed widespread residue flipping and instances of positional coupling. Application of our thermodynamic model to published in vivo crosslinking data reveals quantitative agreement between predicted affinities and in vivo occupancies. Our analyses suggest a thermodynamically driven, continuous Pumilio-binding landscape that is negligibly affected by RNA structure or kinetic factors, such as displacement by ribosomes. This work provides a quantitative foundation for dissecting the cellular behavior of RBPs and cellular features that impact their occupancies.

Crowdsourced RNA design discovers diverse, reversible, efficient, self-contained molecular switches.

Internet-based scientific communities promise a means to apply distributed, diverse human intelligence toward previously intractable scientific problems. However, current implementations have not allowed communities to propose experiments to test all emerging hypotheses at scale or to modify hypotheses in response to experiments. We report high-throughput methods for molecular characterization of nucleic acids that enable the large-scale video game­based crowdsourcing of RNA sensor design, followed by high-throughput functional characterization. Iterative design testing of thousands of crowdsourced RNA sensor designs produced near­thermodynamically optimal and reversible RNA switches that act as self-contained molecular sensors and couple five distinct small molecule inputs to three distinct protein binding and fluorogenic outputs. This work suggests a paradigm for widely distributed experimental bioscience.

Intraflagellar transport velocity is governed by the number of active KIF17 and KIF3AB motors and their motility properties under load.

Homodimeric KIF17 and heterotrimeric KIF3AB are processive, kinesin-2 family motors that act jointly to carry out anterograde intraflagellar transport (IFT), ferrying cargo along microtubules (MTs) toward the tips of cilia. How IFT trains attain speeds that exceed the unloaded rate of the slower, KIF3AB motor remains unknown. By characterizing the motility properties of kinesin-2 motors as a function of load we find that the increase in KIF3AB velocity, elicited by forward loads from KIF17 motors, cannot alone account for the speed of IFT trains in vivo. Instead, higher IFT velocities arise from an increased likelihood that KIF3AB motors dissociate from the MT, resulting in transport by KIF17 motors alone, unencumbered by opposition from KIF3AB. The rate of transport is therefore set by an equilibrium between a faster state, where only KIF17 motors move the train, and a slower state, where at least one KIF3AB motor on the train remains active in transport. The more frequently the faster state is accessed, the higher the overall velocity of the IFT train. We conclude that IFT velocity is governed by (i) the absolute numbers of each motor type on a given train, (ii) how prone KIF3AB is to dissociation from MTs relative to KIF17, and (iii) how prone both motors are to dissociation relative to binding MTs.

Comprehensive and quantitative mapping of RNA-protein interactions across a transcribed eukaryotic genome.

RNA-binding proteins (RBPs) control the fate of nearly every transcript in a cell. However, no existing approach for studying these posttranscriptional gene regulators combines transcriptome-wide throughput and biophysical precision. Here, we describe an assay that accomplishes this. Using commonly available hardware, we built a customizable, open-source platform that leverages the inherent throughput of Illumina technology for direct biophysical measurements. We used the platform to quantitatively measure the binding affinity of the prototypical RBP Vts1 for every transcript in the Saccharomyces cerevisiae genome. The scale and precision of these measurements revealed many previously unknown features of this well-studied RBP. Our transcribed genome array (TGA) assayed both rare and abundant transcripts with equivalent proficiency, revealing hundreds of low-abundance targets missed by previous approaches. These targets regulated diverse biological processes including nutrient sensing and the DNA damage response, and implicated Vts1 in de novo gene "birth." TGA provided single-nucleotide resolution for each binding site and delineated a highly specific sequence and structure motif for Vts1 binding. Changes in transcript levels in vts1Δ cells established the regulatory function of these binding sites. The impact of Vts1 on transcript abundance was largely independent of where it bound within an mRNA, challenging prevailing assumptions about how this RBP drives RNA degradation. TGA thus enables a quantitative description of the relationship between variant RNA structures, affinity, and in vivo phenotype on a transcriptome-wide scale. We anticipate that TGA will provide similarly comprehensive and quantitative insights into the function of virtually any RBP.

High-throughput biochemical profiling reveals sequence determinants of dCas9 off-target binding and unbinding.

The bacterial adaptive immune system CRISPR-Cas9 has been appropriated as a versatile tool for editing genomes, controlling gene expression, and visualizing genetic loci. To analyze Cas9's ability to bind DNA rapidly and specifically, we generated multiple libraries of potential binding partners for measuring the kinetics of nuclease-dead Cas9 (dCas9) interactions. Using a massively parallel method to quantify protein-DNA interactions on a high-throughput sequencing flow cell, we comprehensively assess the effects of combinatorial mismatches between guide RNA (gRNA) and target nucleotides, both in the seed and in more distal nucleotides, plus disruption of the protospacer adjacent motif (PAM). We report two consequences of PAM-distal mismatches: reversal of dCas9 binding at long time scales, and synergistic changes in association kinetics when other gRNA-target mismatches are present. Together, these observations support a model for Cas9 specificity wherein gRNA-DNA mismatches at PAM-distal bases modulate different biophysical parameters that determine association and dissociation rates. The methods we present decouple aspects of kinetic and thermodynamic properties of the Cas9-DNA interaction and broaden the toolkit for investigating off-target binding behavior.

Kinesin processivity is gated by phosphate release.

Kinesin-1 is a dimeric motor protein, central to intracellular transport, that steps hand-over-hand toward the microtubule (MT) plus-end, hydrolyzing one ATP molecule per step. Its remarkable processivity is critical for ferrying cargo within the cell: over 100 successive steps are taken, on average, before dissociation from the MT. Despite considerable work, it is not understood which features coordinate, or "gate," the mechanochemical cycles of the two motor heads. Here, we show that kinesin dissociation occurs subsequent to, or concomitant with, phosphate (P(i)) release following ATP hydrolysis. In optical trapping experiments, we found that increasing the steady-state population of the posthydrolysis ADP · P(i) state (by adding free P(i)) nearly doubled the kinesin run length, whereas reducing either the ATP binding rate or hydrolysis rate had no effect. The data suggest that, during processive movement, tethered-head binding occurs subsequent to hydrolysis, rather than immediately after ATP binding, as commonly suggested. The structural change driving motility, thought to be neck linker docking, is therefore completed only upon hydrolysis, and not ATP binding. Our results offer additional insights into gating mechanisms and suggest revisions to prevailing models of the kinesin reaction cycle.

Compact RNA sensors for increasingly complex functions of multiple inputs.

Designing single molecules that compute general functions of input molecular partners represents a major unsolved challenge in molecular design. Here, we demonstrate that high-throughput, iterative experimental testing of diverse RNA designs crowdsourced from Eterna yields sensors of increasingly complex functions of input oligonucleotide concentrations. After designing single-input RNA sensors with activation ratios beyond our detection limits, we created logic gates, including challenging XOR and XNOR gates, and sensors that respond to the ratio of two inputs. Finally, we describe the OpenTB challenge, which elicited 85-nucleotide sensors that compute a score for diagnosing active tuberculosis, based on the ratio of products of three gene segments. Building on OpenTB design strategies, we created an algorithm Nucleologic that produces similarly compact sensors for the three-gene score based on RNA and DNA. These results open new avenues for diverse applications of compact, single molecule sensors previously limited by design complexity.

Comprehensive sequence-to-function mapping of cofactor-dependent RNA catalysis in the glmS ribozyme.

Massively parallel, quantitative measurements of biomolecular activity across sequence space can greatly expand our understanding of RNA sequence-function relationships. We report the development of an RNA-array assay to perform such measurements and its application to a model RNA: the core glmS ribozyme riboswitch, which performs a ligand-dependent self-cleavage reaction. We measure the cleavage rates for all possible single and double mutants of this ribozyme across a series of ligand concentrations, determining kcat and KM values for active variants. These systematic measurements suggest that evolutionary conservation in the consensus sequence is driven by maintenance of the cleavage rate. Analysis of double-mutant rates and associated mutational interactions produces a structural and functional mapping of the ribozyme sequence, revealing the catalytic consequences of specific tertiary interactions, and allowing us to infer structural rearrangements that permit certain sequence variants to maintain activity.

Automated Design of Diverse Stand-Alone Riboswitches.

Riboswitches that couple binding of ligands to conformational changes offer sensors and control elements for RNA synthetic biology and medical biotechnology. However, design of these riboswitches has required expert intuition or software specialized to transcription or translation outputs; design has been particularly challenging for applications in which the riboswitch output cannot be amplified by other molecular machinery. We present a fully automated design method called RiboLogic for such "stand-alone" riboswitches and test it via high-throughput experiments on 2875 molecules using RNA-MaP (RNA on a massively parallel array) technology. These molecules consistently modulate their affinity to the MS2 bacteriophage coat protein upon binding of flavin mononucleotide, tryptophan, theophylline, and microRNA miR-208a, achieving activation ratios of up to 20 and significantly better performance than control designs. By encompassing a wide diversity of stand-alone switches and highly quantitative data, the resulting ribologic-solves experimental data set provides a rich resource for further improvement of riboswitch models and design methods.

The Mechanochemical Cycle of Mammalian Kinesin-2 KIF3A/B under Load.

The response of motor proteins to external loads underlies their ability to work in teams and determines the net speed and directionality of cargo transport. The mammalian kinesin-2, KIF3A/B, is a heterotrimeric motor involved in intraflagellar transport and vesicle motility in neurons. Bidirectional cargo transport is known to result from the opposing activities of KIF3A/B and dynein bound to the same cargo, but the load-dependent properties of kinesin-2 are poorly understood. We used a feedback-controlled optical trap to probe the velocity, run length, and unbinding kinetics of mouse KIF3A/B under various loads and nucleotide conditions. The kinesin-2 motor velocity is less sensitive than kinesin-1 to external forces, but its processivity diminishes steeply with load, and the motor was observed occasionally to slip and reattach. Each motor domain was characterized by studying homodimeric constructs, and a global fit to the data resulted in a comprehensive pathway that quantifies the principal force-dependent kinetic transitions. The properties of the KIF3A/B heterodimer are intermediate between the two homodimers, and the distinct load-dependent behavior is attributable to the properties of the motor domains and not to the neck linkers or the coiled-coil stalk. We conclude that the force-dependent movement of KIF3A/B differs significantly from conventional kinesin-1. Against opposing dynein forces, KIF3A/B motors are predicted to rapidly unbind and rebind, resulting in qualitatively different transport behavior from kinesin-1.

Examining kinesin processivity within a general gating framework.

Kinesin-1 is a dimeric motor that transports cargo along microtubules, taking 8.2-nm steps in a hand-over-hand fashion. The ATP hydrolysis cycles of its two heads are maintained out of phase by a series of gating mechanisms, which lead to processive runs averaging ~1 µm. A key structural element for inter-head coordination is the neck linker (NL), which connects the heads to the stalk. To examine the role of the NL in regulating stepping, we investigated NL mutants of various lengths using single-molecule optical trapping and bulk fluorescence approaches in the context of a general framework for gating. Our results show that, although inter-head tension enhances motor velocity, it is crucial neither for inter-head coordination nor for rapid rear-head release. Furthermore, cysteine-light mutants do not produce wild-type motility under load. We conclude that kinesin-1 is primarily front-head gated, and that NL length is tuned to enhance unidirectional processivity and velocity.

A universal pathway for kinesin stepping.

Kinesin-1 is an ATP-driven, processive motor that transports cargo along microtubules in a tightly regulated stepping cycle. Efficient gating mechanisms ensure that the sequence of kinetic events proceeds in the proper order, generating a large number of successive reaction cycles. To study gating, we created two mutant constructs with extended neck-linkers and measured their properties using single-molecule optical trapping and ensemble fluorescence techniques. Owing to a reduction in the inter-head tension, the constructs access an otherwise rarely populated conformational state in which both motor heads remain bound to the microtubule. ATP-dependent, processive backstepping and futile hydrolysis were observed under moderate hindering loads. On the basis of measurements, we formulated a comprehensive model for kinesin motion that incorporates reaction pathways for both forward and backward stepping. In addition to inter-head tension, we found that neck-linker orientation is also responsible for ensuring gating in kinesin.

An optical apparatus for rotation and trapping.

We present details of the design, construction, and testing of a single-beam optical tweezers apparatus capable of measuring and exerting torque, as well as force, on microfabricated, optically anisotropic particles (an "optical torque wrench"). The control of angular orientation is achieved by rotating the linear polarization of a trapping laser with an electro-optic modulator (EOM), which affords improved performance over previous designs. The torque imparted to the trapped particle is assessed by measuring the difference between left- and right-circular components of the transmitted light, and constant torque is maintained by feeding this difference signal back into a custom-designed electronic servo loop. The limited angular range of the EOM (+/-180 degrees ) is extended by rapidly reversing the polarization once a threshold angle is reached, enabling the torque clamp to function over unlimited, continuous rotations at high bandwidth. In addition, we developed particles suitable for rotation in this apparatus using microfabrication techniques. Altogether, the system allows for the simultaneous application of forces (approximately 0.1-100 pN) and torques (approximately 1-10,000 pN nm) in the study of biomolecules. As a proof of principle, we demonstrate how our instrument can be used to study the supercoiling of single DNA molecules.