Learning rates and states from biophysical time series: a Bayesian approach to model selection and single-molecule FRET data

Simple methods to determine the dissociation constant, Kd

The determination of the dissociation constant (Kd) is pivotal in biochemistry and pharmacology for understanding binding affinities in chemical reactions, which is crucial for drug development and comprehending biological systems. Here, we introduce a single-molecule fluorescence resonance energy transfer-based method for determining Kd, alongside the conventional electrophoretic mobility shift assay method of Kd, offering insights into thermodynamic interactions between proteins and substrates. The single-molecule fluorescence resonance energy transfer approach is highlighted for its ability to accurately measure binding and dissociation kinetics through fluorescence labeling and the intrinsic nature of protein-DNA interactions, representing a significant advancement in the fields of molecular biology and pharmacology.

Directional change during active diffusion of viral ribonucleoprotein particles through cytoplasm

A mesh of cytoskeletal fibers, consisting of microtubules, intermediate filaments, and fibrous actin, prevents the Brownian diffusion of particles with a diameter larger than 0.10 μm, such as vesicular stomatitis virus ribonucleoprotein (RNP) particles, in mammalian cells. Nevertheless, RNP particles do move in random directions but at a lower rate than Brownian diffusion, which is thermally driven. This nonthermal biological transport process is called "active diffusion" because it is driven by ATP. The ATP powers motor proteins such as myosin II. The motor proteins bend and cross-link actin fibers, causing the mesh to jiggle. Until recently, little was known about how RNP particles get through the mesh. It has been customary to analyze the tracks of particles like RNPs by computing the slope of the ensemble-averaged mean-squared displacement of the particles as a signature of mechanism. Although widely used, this approach "loses information" about the timing of the switches between physical mechanisms. It has been recently shown that machine learning composed of variational Bayesian analysis, Gaussian mixture models, and hidden Markov models can use "all the information" in a single track to reveal that that the positions of RNP particles are spatially clustered. Machine learning assigns a number, called a state, to each cluster. RNP particles remain in one state for 0.2-1.0 s before switching (hopping) to a different state. This earlier work is here extended to analyze the movements of a particle within a state and to determine particle directionality within and between states.

Increasing the accuracy of single-molecule data analysis using tMAVEN

Time-dependent single-molecule experiments contain rich kinetic information about the functional dynamics of biomolecules. A key step in extracting this information is the application of kinetic models, such as hidden Markov models (HMMs), which characterize the molecular mechanism governing the experimental system. Unfortunately, researchers rarely know the physicochemical details of this molecular mechanism a priori, which raises questions about how to select the most appropriate kinetic model for a given single-molecule data set and what consequences arise if the wrong model is chosen. To address these questions, we have developed and used time-series modeling, analysis, and visualization environment (tMAVEN), a comprehensive, open-source, and extensible software platform. tMAVEN can perform each step of the single-molecule analysis pipeline, from preprocessing to kinetic modeling to plotting, and has been designed to enable the analysis of a single-molecule data set with multiple types of kinetic models. Using tMAVEN, we have systematically investigated mismatches between kinetic models and molecular mechanisms by analyzing simulated examples of prototypical single-molecule data sets exhibiting common experimental complications, such as molecular heterogeneity, with a series of different types of HMMs. Our results show that no single kinetic modeling strategy is mathematically appropriate for all experimental contexts. Indeed, HMMs only correctly capture the underlying molecular mechanism in the simplest of cases. As such, researchers must modify HMMs using physicochemical principles to avoid the risk of missing the significant biological and biophysical insights into molecular heterogeneity that their experiments provide. By enabling the facile, side-by-side application of multiple types of kinetic models to individual single-molecule data sets, tMAVEN allows researchers to carefully tailor their modeling approach to match the complexity of the underlying biomolecular dynamics and increase the accuracy of their single-molecule data analyses.

Stabilization of the hexasome intermediate during histone exchange by yeast SWR1 complex

The yeast SWR1 complex catalyzes the exchange of histone H2A/H2B dimers in nucleosomes with Htz1/H2B dimers. We use cryoelectron microscopy to determine the structure of an enzyme-bound hexasome intermediate in the reaction pathway of histone exchange, in which an H2A/H2B dimer has been extracted from a nucleosome prior to the insertion of a dimer comprising Htz1/H2B. The structure reveals a key role for the Swc5 subunit in stabilizing the unwrapping of DNA from the histone core of the hexasome. By engineering a crosslink between an Htz1/H2B dimer and its chaperone protein Chz1, we show that this blocks histone exchange by SWR1 but allows the incoming chaperone-dimer complex to insert into the hexasome. We use this reagent to trap an SWR1/hexasome complex with an incoming Htz1/H2B dimer that shows how the reaction progresses to the next step. Taken together the structures reveal insights into the mechanism of histone exchange by SWR1 complex.

Foundation model for efficient biological discovery in single-molecule data

Modern data-intensive techniques offer ever deeper insights into biology, but render the process of discovery increasingly complex. For example, exploiting the unique ability of single-molecule fluorescence microscopy (SMFM)1-5. to uncover rare but critical intermediates often demands manual inspection of time traces and iterative ad hoc approaches that are difficult to systematize. To facilitate systematic and efficient discovery from SMFM data, we introduce META-SiM, a transformer-based foundation model pre-trained on diverse SMFM analysis tasks. META-SiM achieves high performance-rivaling best-in-class algorithms-on a broad range of analysis tasks including trace selection, classification, segmentation, idealization, and stepwise photobleaching analysis. Additionally, the model produces high-dimensional embedding vectors that encapsulate detailed information about each trace, which the web-based META-SiM Projector (https://www.simol-projector.org) casts into lower-dimensional space for efficient whole-dataset visualization, labeling, comparison, and sharing. Combining this Projector with the objective metric of Local Shannon Entropy enables rapid identification of condition-specific behaviors, even if rare or subtle. As a result, by applying META-SiM to an existing single-molecule Förster resonance energy transfer (smFRET) dataset6, we discover a previously unobserved intermediate state in pre-mRNA splicing. META-SiM thus removes bottlenecks, improves objectivity, and both systematizes and accelerates biological discovery in complex single-molecule data.

eIF1 and eIF5 dynamically control translation start site fidelity

Translation initiation defines the identity of a synthesized protein through selection of a translation start site on a messenger RNA. This process is essential to well-controlled protein synthesis, modulated by stress responses, and dysregulated in many human diseases. The eukaryotic initiation factors eIF1 and eIF5 interact with the initiator methionyl-tRNAi Met on the 40S ribosomal subunit to coordinate start site selection. Here, using single-molecule analysis of in vitro reconstituted human initiation combined with translation assays in cells, we examine eIF1 and eIF5 function. During translation initiation on a panel of RNAs, we monitored both proteins directly and in real time using single-molecule fluorescence. As expected, eIF1 loaded onto mRNAs as a component of the 43S initiation complex. Rapid (~ 2 s) eIF1 departure required a translation start site and was delayed by alternative start sites and a longer 5' untranslated region (5'UTR). After its initial departure, eIF1 rapidly and transiently sampled initiation complexes, with more prolonged sampling events on alternative start sites. By contrast, eIF5 only transiently bound initiation complexes late in initiation immediately prior to association of eIF5B, which allowed joining of the 60S ribosomal subunit. eIF5 association required the presence of a translation start site and was inhibited and destabilized by alternative start sites. Using both knockdown and overexpression experiments in human cells, we validated that eIF1 and eIF5 have opposing roles during initiation. Collectively, our findings demonstrate how multiple eIF1 and eIF5 binding events control start-site selection fidelity throughout initiation, which is tuned in response to changes in the levels of both proteins.

Sequence-specific dynamic DNA bending explains mitochondrial TFAM's dual role in DNA packaging and transcription initiation

Mitochondrial transcription factor A (TFAM) employs DNA bending to package mitochondrial DNA (mtDNA) into nucleoids and recruit mitochondrial RNA polymerase (POLRMT) at specific promoter sites, light strand promoter (LSP) and heavy strand promoter (HSP). Herein, we characterize the conformational dynamics of TFAM on promoter and non-promoter sequences using single-molecule fluorescence resonance energy transfer (smFRET) and single-molecule protein-induced fluorescence enhancement (smPIFE) methods. The DNA-TFAM complexes dynamically transition between partially and fully bent DNA conformational states. The bending/unbending transition rates and bending stability are DNA sequence-dependent-LSP forms the most stable fully bent complex and the non-specific sequence the least, which correlates with the lifetimes and affinities of TFAM with these DNA sequences. By quantifying the dynamic nature of the DNA-TFAM complexes, our study provides insights into how TFAM acts as a multifunctional protein through the DNA bending states to achieve sequence specificity and fidelity in mitochondrial transcription while performing mtDNA packaging.

RNase E searches for cleavage sites in RNA by linear diffusion: direct evidence from single-molecule FRET

The ability of obstacles in cellular transcripts to protect downstream but not upstream sites en masse from attack by RNase E has prompted the hypothesis that this mRNA-degrading endonuclease may scan 5'-monophosphorylated RNA linearly for cleavage sites, starting at the 5' end. However, despite its proposed regulatory importance, the migration of RNase E on RNA has never been directly observed. We have now used single-molecule FRET to monitor the dynamics of this homotetrameric enzyme on RNA. Our findings reveal that RNase E slides along unpaired regions of RNA without consuming a molecular source of energy such as ATP and that its forward progress can be impeded when it encounters a large structural discontinuity. This movement, which is bidirectional, occurs in discrete steps of variable length and requires an RNA ligand much longer than needed to occupy a single RNase E subunit. These results indicate that RNase E scans for cleavage sites by one-dimensional diffusion and suggest a possible molecular mechanism.

Short-range translocation by a restriction enzyme motor triggers diffusion along DNA

Cleavage of bacteriophage DNA by the Type III restriction-modification enzymes requires long-range interaction between DNA sites. This is facilitated by one-dimensional diffusion ('DNA sliding') initiated by ATP hydrolysis catalyzed by a superfamily 2 helicase-like ATPase. Here we combined ultrafast twist measurements based on plasmonic DNA origami nano-rotors with stopped-flow fluorescence and gel-based assays to examine the role(s) of ATP hydrolysis. Our data show that the helicase-like domain has multiple roles. First, this domain stabilizes initial DNA interactions alongside the methyltransferase subunits. Second, it causes environmental changes in the flipped adenine base following hydrolysis of the first ATP. Finally, it remodels nucleoprotein interactions via constrained translocation of a ∼ 5 to 22-bp double stranded DNA loop. Initiation of DNA sliding requires 8-15 bp of DNA downstream of the motor, corresponding to the site of nuclease domain binding. Our data unify previous contradictory communication models for Type III enzymes.

Stepwise activation of a metabotropic glutamate receptor

Metabotropic glutamate receptors belong to a family of G protein-coupled receptors that are obligate dimers and possess a large extracellular ligand-binding domain that is linked via a cysteine-rich domain to their 7-transmembrane domain1. Upon activation, these receptors undergo a large conformational change to transmit the ligand binding signal from the extracellular ligand-binding domain to the G protein-coupling 7-transmembrane domain2. In this manuscript, we propose a model for a sequential, multistep activation mechanism of metabotropic glutamate receptor subtype 5. We present a series of structures in lipid nanodiscs, from inactive to fully active, including agonist-bound intermediate states. Further, using bulk and single-molecule fluorescence imaging, we reveal distinct receptor conformations upon allosteric modulator and G protein binding.

Bayesian Mechanistic Inference, Statistical Mechanics, and a New Era for Monte Carlo

On the one hand, much of computational chemistry is concerned with "bottom-up" calculations which elucidate observable behavior starting from exact or approximated physical laws, a paradigm exemplified by typical quantum mechanical calculations and molecular dynamics simulations. On the other hand, "top down" computations aiming to formulate mathematical models consistent with observed data, e.g., parametrizing force fields, binding or kinetic models, have been of interest for decades but recently have grown in sophistication with the use of Bayesian inference (BI). Standard BI provides an estimation of parameter values, uncertainties, and correlations among parameters. Used for "model selection," BI can also distinguish between model structures such as the presence or absence of individual states and transitions. Fortunately for physical scientists, BI can be formulated within a statistical mechanics framework, and indeed, BI has led to a resurgence of interest in Monte Carlo (MC) algorithms, many of which have been directly adapted from or inspired by physical strategies. Certain MC algorithms─notably procedures using an "infinite temperature" reference state─can be successful in a 5-20 parameter BI context which would be unworkable in molecular spaces of 103 coordinates and more. This Review provides a pedagogical introduction to BI and reviews key aspects of BI through a physical lens, setting the computations in terms of energy landscapes and free energy calculations and describing promising sampling algorithms. Statistical mechanics and basic probability theory also provide a reference for understanding intrinsic limitations of Bayesian inference with regard to model selection and the choice of priors.

A Sequential Binding Mechanism for 5' Splice Site Recognition and Modulation for the Human U1 snRNP

Splice site recognition is essential for defining the transcriptome. Drugs like risdiplam and branaplam change how U1 snRNP recognizes particular 5' splice sites (5'SS) and promote U1 snRNP binding and splicing at these locations. Despite the therapeutic potential of 5'SS modulators, the complexity of their interactions and snRNP substrates have precluded defining a mechanism for 5'SS modulation. We have determined a sequential binding mechanism for modulation of -1A bulged 5'SS by branaplam using a combination of ensemble kinetic measurements and colocalization single molecule spectroscopy (CoSMoS). Our mechanism establishes that U1-C protein binds reversibly to U1 snRNP, and branaplam binds to the U1 snRNP/U1-C complex only after it has engaged a -1A bulged 5'SS. Obligate orders of binding and unbinding explain how reversible branaplam interactions cause formation of long-lived U1 snRNP/5'SS complexes. Branaplam is a ribonucleoprotein, not RNA duplex alone, targeting drug whose action depends on fundamental properties of 5'SS recognition.

Mechanisms of active diffusion of vesicular stomatitis virus inclusion bodies and cellular early endosomes in the cytoplasm of mammalian cells

Viral and cellular particles too large to freely diffuse have two different types of mobility in the eukaryotic cell cytoplasm: directed motion mediated by motor proteins moving along cytoskeletal elements with the particle as its load, and motion in random directions mediated by motor proteins interconnecting cytoskeletal elements. The latter motion is referred to as "active diffusion." Mechanisms of directed motion have been extensively studied compared to mechanisms of active diffusion, despite the observation that active diffusion is more common for many viral and cellular particles. Our previous research showed that active diffusion of vesicular stomatitis virus (VSV) ribonucleoproteins (RNPs) in the cytoplasm consists of hopping between traps and that actin filaments and myosin II motors are components of the hop-trap mechanism. This raises the question whether similar mechanisms mediate random motion of larger particles with different physical and biological properties. Live-cell fluorescence imaging and a variational Bayesian analysis used in pattern recognition and machine learning were used to determine the molecular mechanisms of random motion of VSV inclusion bodies and cellular early endosomes. VSV inclusion bodies are membraneless cellular compartments that are the major sites of viral RNA synthesis, and early endosomes are representative of cellular membrane-bound organelles. Like VSV RNPs, inclusion bodies and early endosomes moved from one trapped state to another, but the distance between states was inconsistent with hopping between traps, indicating that the apparent state-to-state movement is mediated by trap movement. Like VSV RNPs, treatment with the actin filament depolymerizing inhibitor latrunculin A increased VSV inclusion body mobility by increasing the size of the traps. In contrast neither treatment with latrunculin A nor depolymerization of microtubules by nocodazole treatment affected the size of traps that confine early endosome mobility, indicating that intermediate filaments are likely major trap components for these cellular organelles.

Structural analysis and conformational dynamics of a holo-adhesion GPCR reveal interplay between extracellular and transmembrane domains

Adhesion G Protein-Coupled Receptors (aGPCRs) are key cell-adhesion molecules involved in numerous physiological functions. aGPCRs have large multi-domain extracellular regions (ECR) containing a conserved GAIN domain that precedes their seven-pass transmembrane domain (7TM). Ligand binding and mechanical force applied on the ECR regulate receptor function. However, how the ECR communicates with the 7TM remains elusive, because the relative orientation and dynamics of the ECR and 7TM within a holoreceptor is unclear. Here, we describe the cryo-EM reconstruction of an aGPCR, Latrophilin3/ADGRL3, and reveal that the GAIN domain adopts a parallel orientation to the membrane and has constrained movement. Single-molecule FRET experiments unveil three slow-exchanging FRET states of the ECR relative to the 7TM within the holoreceptor. GAIN-targeted antibodies, and cancer-associated mutations at the GAIN-7TM interface, alter FRET states, cryo-EM conformations, and receptor signaling. Altogether, this data demonstrates conformational and functional coupling between the ECR and 7TM, suggesting an ECR-mediated mechanism of aGPCR activation.

Integrating magnetic tweezers and single-molecule FRET: A comprehensive approach to studying local and global conformational changes simultaneously

This chapter presents the integration of magnetic tweezers with single-molecule FRET technology, a significant advancement in the study of nucleic acids and other biological systems. We detail the technical aspects, challenges, and current status of this hybrid technique, which combines the global manipulation and observation capabilities of magnetic tweezers with the local conformational detection of smFRET. This innovative approach enhances our ability to analyze and understand the molecular mechanics of biological systems. The chapter serves as our first formal documentation of this method, offering insights and methodologies developed in our laboratory over the past decade.

Increasing the accuracy of single-molecule data analysis using tMAVEN

Time-dependent single-molecule experiments contain rich kinetic information about the functional dynamics of biomolecules. A key step in extracting this information is the application of kinetic models, such as hidden Markov models (HMMs), which characterize the molecular mechanism governing the experimental system. Unfortunately, researchers rarely know the physico-chemical details of this molecular mechanism a priori, which raises questions about how to select the most appropriate kinetic model for a given single-molecule dataset and what consequences arise if the wrong model is chosen. To address these questions, we have developed and used time-series Modeling, Analysis, and Visualization ENvironment (tMAVEN), a comprehensive, open-source, and extensible software platform. tMAVEN can perform each step of the single-molecule analysis pipeline, from pre-processing to kinetic modeling to plotting, and has been designed to enable the analysis of a single-molecule dataset with multiple types of kinetic models. Using tMAVEN, we have systematically investigated mismatches between kinetic models and molecular mechanisms by analyzing simulated examples of prototypical single-molecule datasets exhibiting common experimental complications, such as molecular heterogeneity, with a series of different types of HMMs. Our results show that no single kinetic modeling strategy is mathematically appropriate for all experimental contexts. Indeed, HMMs only correctly capture the underlying molecular mechanism in the simplest of cases. As such, researchers must modify HMMs using physico-chemical principles to avoid the risk of missing the significant biological and biophysical insights into molecular heterogeneity that their experiments provide. By enabling the facile, side-by-side application of multiple types of kinetic models to individual single-molecule datasets, tMAVEN allows researchers to carefully tailor their modeling approach to match the complexity of the underlying biomolecular dynamics and increase the accuracy of their single-molecule data analyses.

Inferring pointwise diffusion properties of single trajectories with deep learning

To characterize the mechanisms governing the diffusion of particles in biological scenarios, it is essential to accurately determine their diffusive properties. To do so, we propose a machine-learning method to characterize diffusion processes with time-dependent properties at the experimental time resolution. Our approach operates at the single-trajectory level predicting the properties of interest, such as the diffusion coefficient or the anomalous diffusion exponent, at every time step of the trajectory. In this way, changes in the diffusive properties occurring along the trajectory emerge naturally in the prediction and thus allow the characterization without any prior knowledge or assumption about the system. We first benchmark the method on synthetic trajectories simulated under several conditions. We show that our approach can successfully characterize both abrupt and continuous changes in the diffusion coefficient or the anomalous diffusion exponent. Finally, we leverage the method to analyze experiments of single-molecule diffusion of two membrane proteins in living cells: the pathogen-recognition receptor DC-SIGN and the integrin α5β1. The analysis allows us to characterize physical parameters and diffusive states with unprecedented accuracy, shedding new light on the underlying mechanisms.

The mechanism of mRNA activation

During translation initiation, messenger RNA molecules must be identified and activated for loading into a ribosome. In this rate-limiting step, the heterotrimeric protein eukaryotic initiation factor eIF4F must recognize and productively interact with the 7-methylguanosine cap at the 5' end of the messenger RNA and subsequently activate the message. Despite its fundamental, regulatory role in gene expression, the molecular events underlying cap recognition and messenger RNA activation remain mysterious. Here, we generate a unique, single-molecule fluorescence imaging system to interrogate the dynamics with which eIF4F discriminates productive and non-productive locations on full-length, native messenger RNA molecules. At the single-molecule level, we observe stochastic sampling of eIF4F along the length of the messenger RNA and identify allosteric communication between the eIF4F subunits which ultimately drive cap-recognition and subsequent activation of the message. Our experiments uncover novel functions for each subunit of eIF4F and we conclude by presenting a model for messenger RNA activation which precisely defines the composition of the activated message. This model provides a general framework for understanding how messenger RNA molecules may be discriminated from one another, and how other RNA-binding proteins may control the efficiency of translation initiation.

Modulation of protein-DNA binding reveals mechanisms of spatiotemporal gene control in early Drosophila embryos

It is well known that enhancers regulate the spatiotemporal expression of their target genes by recruiting transcription factors (TFs) to the cognate binding sites in the region. However, the role of multiple binding sites for the same TFs and their specific spatial arrangement in determining the overall competency of the enhancer has yet to be fully understood. In this study, we utilized the MS2-MCP live imaging technique to quantitatively analyze the regulatory logic of the snail distal enhancer in early Drosophila embryos. Through systematic modulation of Dorsal and Twist binding motifs in this enhancer, we found that a mutation in any one of these binding sites causes a drastic reduction in transcriptional amplitude, resulting in a reduction in mRNA production of the target gene. We provide evidence of synergy, such that multiple binding sites with moderate affinities cooperatively recruit more TFs to drive stronger transcriptional activity than a single site. Moreover, a Hidden Markov-based stochastic model of transcription reveals that embryos with mutated binding sites have a higher probability of returning to the inactive promoter state. We propose that TF-DNA binding regulates spatial and temporal gene expression and drives robust pattern formation by modulating transcriptional kinetics and tuning bursting rates.

DNA strand breaks and gaps target retroviral intasome binding and integration

Retrovirus integration into a host genome is essential for productive infections. The integration strand transfer reaction is catalyzed by a nucleoprotein complex (Intasome) containing the viral integrase (IN) and the reverse transcribed (RT) copy DNA (cDNA). Previous studies suggested that DNA target-site recognition limits intasome integration. Using single molecule Förster resonance energy transfer (smFRET), we show prototype foamy virus (PFV) intasomes specifically bind to DNA strand breaks and gaps. These break and gap DNA discontinuities mimic oxidative base excision repair (BER) lesion-processing intermediates that have been shown to affect retrovirus integration in vivo. The increased DNA binding events targeted strand transfer to the break/gap site without inducing substantial intasome conformational changes. The major oxidative BER substrate 8-oxo-guanine as well as a G/T mismatch or +T nucleotide insertion that typically introduce a bend or localized flexibility into the DNA, did not increase intasome binding or targeted integration. These results identify DNA breaks or gaps as modulators of dynamic intasome-target DNA interactions that encourage site-directed integration.

Single-molecule dynamics show a transient lipopolysaccharide transport bridge

Gram-negative bacteria are surrounded by two membranes. A special feature of the outer membrane is its asymmetry. It contains lipopolysaccharide (LPS) in the outer leaflet and phospholipids in the inner leaflet1-3. The proper assembly of LPS in the outer membrane is required for cell viability and provides Gram-negative bacteria intrinsic resistance to many classes of antibiotics. LPS biosynthesis is completed in the inner membrane, so the LPS must be extracted, moved across the aqueous periplasm that separates the two membranes and translocated through the outer membrane where it assembles on the cell surface4. LPS transport and assembly requires seven conserved and essential LPS transport components5 (LptA-G). This system has been proposed to form a continuous protein bridge that provides a path for LPS to reach the cell surface6,7, but this model has not been validated in living cells. Here, using single-molecule tracking, we show that Lpt protein dynamics are consistent with the bridge model. Half of the inner membrane Lpt proteins exist in a bridge state, and bridges persist for 5-10 s, showing that their organization is highly dynamic. LPS facilitates Lpt bridge formation, suggesting a mechanism by which the production of LPS can be directly coupled to its transport. Finally, the bridge decay kinetics suggest that there may be two different types of bridges, whose stability differs according to the presence (long-lived) or absence (short-lived) of LPS. Together, our data support a model in which LPS is both a substrate and a structural component of dynamic Lpt bridges that promote outer membrane assembly.

Zero-Mode Waveguide Nanowells for Single-Molecule Detection in Living Cells

Single-molecule fluorescence imaging experiments generally require sub-nanomolar protein concentrations to isolate single protein molecules, which makes such experiments challenging in live cells due to high intracellular protein concentrations. Here, we show that single-molecule observations can be achieved in live cells through a drastic reduction in the observation volume using overmilled zero-mode waveguides (ZMWs- subwavelength-size holes in a metal film). Overmilling of the ZMW in a palladium film creates a nanowell of tunable size in the glass layer below the aperture, which cells can penetrate. We present a thorough theoretical and experimental characterization of the optical properties of these nanowells over a wide range of ZMW diameters and overmilling depths, showing an excellent signal confinement and a 5-fold fluorescence enhancement of fluorescent molecules inside nanowells. ZMW nanowells facilitate live-cell imaging as cells form stable protrusions into the nanowells. Importantly, the nanowells greatly reduce the cytoplasmic background fluorescence, enabling the detection of individual membrane-bound fluorophores in the presence of high cytoplasmic expression levels, which could not be achieved with TIRF microscopy. Zero-mode waveguide nanowells thus provide great potential to study individual proteins in living cells.

Fluorescence Microscopy: a statistics-optics perspective

Fundamental properties of light unavoidably impose features on images collected using fluorescence microscopes. Modeling these features is ever more important in quantitatively interpreting microscopy images collected at scales on par or smaller than light's wavelength. Here we review the optics responsible for generating fluorescent images, fluorophore properties, microscopy modalities leveraging properties of both light and fluorophores, in addition to the necessarily probabilistic modeling tools imposed by the stochastic nature of light and measurement.

Universal features of Nsp1-mediated translational shutdown by coronaviruses

Nonstructural protein 1 (Nsp1) produced by coronaviruses inhibits host protein synthesis. The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) Nsp1 C-terminal domain was shown to bind the ribosomal mRNA channel to inhibit translation, but it is unclear whether this mechanism is broadly used by coronaviruses, whether the Nsp1 N-terminal domain binds the ribosome, or how Nsp1 allows viral RNAs to be translated. Here, we investigated Nsp1 from SARS-CoV-2, Middle East respiratory syndrome coronavirus (MERS-CoV), and Bat-Hp-CoV coronaviruses using structural, biophysical, and biochemical experiments, revealing a conserved role for the C-terminal domain. Additionally, the N-terminal domain of Bat-Hp-CoV Nsp1 binds to the decoding center of the 40S subunit, where it would prevent mRNA and eIF1A accommodation. Structure-based experiments demonstrated the importance of decoding center interactions in all three coronaviruses and showed that the same regions of Nsp1 are necessary for the selective translation of viral RNAs. Our results provide a mechanistic framework to understand how Nsp1 controls preferential translation of viral RNAs.

Yeast Rad52 is a homodecamer and possesses BRCA2-like bipartite Rad51 binding modes

Homologous recombination (HR) is an essential double-stranded DNA break repair pathway. In HR, Rad52 facilitates the formation of Rad51 nucleoprotein filaments on RPA-coated ssDNA. Here, we decipher how Rad52 functions using single-particle cryo-electron microscopy and biophysical approaches. We report that Rad52 is a homodecameric ring and each subunit possesses an ordered N-terminal and disordered C-terminal half. An intrinsic structural asymmetry is observed where a few of the C-terminal halves interact with the ordered ring. We describe two conserved charged patches in the C-terminal half that harbor Rad51 and RPA interacting motifs. Interactions between these patches regulate ssDNA binding. Surprisingly, Rad51 interacts with Rad52 at two different bindings sites: one within the positive patch in the disordered C-terminus and the other in the ordered ring. We propose that these features drive Rad51 nucleation onto a single position on the DNA to promote formation of uniform pre-synaptic Rad51 filaments in HR.

Hop2-Mnd1 and Swi5-Sfr1 stimulate Dmc1 filament assembly using distinct mechanisms

In meiosis, Dmc1 recombinase and the general recombinase Rad51 are responsible for pairing homologous chromosomes and exchanging strands. Fission yeast (Schizosaccharomyces pombe) Swi5-Sfr1 and Hop2-Mnd1 stimulate Dmc1-driven recombination, but the stimulation mechanism is unclear. Using single-molecule fluorescence resonance energy transfer (smFRET) and tethered particle motion (TPM) experiments, we showed that Hop2-Mnd1 and Swi5-Sfr1 individually enhance Dmc1 filament assembly on single-stranded DNA (ssDNA) and adding both proteins together allows further stimulation. FRET analysis showed that Hop2-Mnd1 enhances the binding rate of Dmc1 while Swi5-Sfr1 specifically reduces the dissociation rate during the nucleation, about 2-fold. In the presence of Hop2-Mnd1, the nucleation time of Dmc1 filaments shortens, and doubling the ss/double-stranded DNA (ss/dsDNA) junctions of DNA substrates reduces the nucleation times in half. Order of addition experiments confirmed that Hop2-Mnd1 binds on DNA to recruit and stimulate Dmc1 nucleation at the ss/dsDNA junction. Our studies directly support the molecular basis of how Hop2-Mnd1 and Swi5-Sfr1 act on different steps during the Dmc1 filament assembly. DNA binding of these accessory proteins and nucleation preferences of recombinases thus dictate how their regulation can take place.

Single-Molecule Reaction-Diffusion

We propose to capture reaction-diffusion on a molecule-by-molecule basis from the fastest acquirable timescale, namely individual photon arrivals. We illustrate our method on intrinsically disordered human proteins, the linker histone H1.0 as well as its chaperone prothymosin α , as these diffuse through an illuminated confocal spot and interact forming larger ternary complexes on millisecond timescales. Most importantly, single-molecule reaction-diffusion, smRD, reveals single molecule properties without trapping or otherwise confining molecules to surfaces. We achieve smRD within a Bayesian paradigm and term our method Bayes-smRD. Bayes-smRD is further free of the average, bulk, results inherent to the analysis of long photon arrival traces by fluorescence correlation spectroscopy. In learning from thousands of photon arrivals continuous spatial positions and discrete conformational and photophysical state changes, Bayes-smRD estimates kinetic parameters on a molecule-by-molecule basis with two to three orders of magnitude less data than tools such as fluorescence correlation spectroscopy thereby also dramatically reducing sample photodamage.

Real-time detection of human telomerase DNA synthesis by multiplexed single-molecule FRET

Genomic stability in proliferating cells critically depends on telomere maintenance by telomerase reverse transcriptase. Here we report the development and proof-of-concept results of a single-molecule approach to monitor the catalytic activity of human telomerase in real time and with single-nucleotide resolution. Using zero-mode waveguides and multicolor FRET, we recorded the processive addition of multiple telomeric repeats to individual DNA primers. Unlike existing biophysical and biochemical tools, the novel approach enables the quantification of nucleotide-binding kinetics before nucleotide incorporation. Moreover, it provides a means to dissect the unique translocation dynamics that telomerase must undergo after synthesis of each hexameric DNA repeat. We observed an unexpectedly prolonged binding dwell time of dGTP in the enzyme active site at the start of each repeat synthesis cycle, suggesting that telomerase translocation is composed of multiple rate-contributing sub-steps that evade classical biochemical analysis.

Modulation of translational decoding by m6A modification of mRNA

N6-methyladenosine (m6A) is an abundant, dynamic mRNA modification that regulates key steps of cellular mRNA metabolism. m6A in the mRNA coding regions inhibits translation elongation. Here, we show how m6A modulates decoding in the bacterial translation system using a combination of rapid kinetics, smFRET and single-particle cryo-EM. We show that, while the modification does not impair the initial binding of aminoacyl-tRNA to the ribosome, in the presence of m6A fewer ribosomes complete the decoding process due to the lower stability of the complexes and enhanced tRNA drop-off. The mRNA codon adopts a π-stacked codon conformation that is remodeled upon aminoacyl-tRNA binding. m6A does not exclude canonical codon-anticodon geometry, but favors alternative more dynamic conformations that are rejected by the ribosome. These results highlight how modifications outside the Watson-Crick edge can still interfere with codon-anticodon base pairing and complex recognition by the ribosome, thereby modulating the translational efficiency of modified mRNAs.

Different conformational dynamics of SNARE protein Ykt6 among yeast and mammals

Ykt6 is one of the most conserved SNARE (N-ethylmaleimide-sensitive factor attachment protein receptor) proteins involved in multiple intracellular membrane trafficking processes. The membrane-anchoring function of Ykt6 has been elucidated to result from its conformational transition from a closed state to an open state. Two ways of regulating the conformational transition were proposed: the C-terminal lipidation and the phosphorylation at the SNARE core. Despite many aspects of common properties, Ykt6 displays differential cellular localizations and functional behaviors in different species, such as yeast, mammals, and worms. The structure-function relationship underlying these differences remains elusive. Here, we combined biochemical characterization, single-molecule FRET measurement, and molecular dynamics simulation to compare the conformational dynamics of yeast and rat Ykt6. Compared to rat Ykt6 (rYkt6), yeast Ykt6 (yYkt6) has more open conformations and could not bind dodecylphosphocholine that inhibits rYkt6 in the closed state. A point mutation T46L/Q57A was shown to be able to convert yYkt6 to a more closed and dodecylphosphocholine-bound state, where Leu46 contributes key hydrophobic interactions for the closed state. We also demonstrated that the phospho-mutation S174D could shift the conformation of rYkt6 to a more open state, but the corresponding mutation S176D in yYkt6 leads to a slightly more closed conformation. These observations shed light on the regulatory mechanism underlying the variations of Ykt6 functions across species.

Conformational dynamics of the RNA binding channel regulates loading and translocation of the DEAH-box helicase Prp43

The DEAH-box helicase Prp43 has essential functions in pre-mRNA splicing and ribosome biogenesis, remodeling structured RNAs. To initiate unwinding, Prp43 must first accommodate a single-stranded RNA segment into its RNA binding channel. This allows translocation of the helicase on the RNA. G-patch (gp) factors activate Prp43 in its cellular context enhancing the intrinsically low ATPase and RNA unwinding activity. It is unclear how the RNA loading process is accomplished by Prp43 and how it is regulated by its substrates, ATP and RNA, and the G-patch partners. We developed single-molecule (sm) FRET reporters on Prp43 from Chaetomium thermophilum to monitor the conformational dynamics of the RNA binding channel in Prp43 in real-time. We show that the channel can alternate between open and closed conformations. Binding of Pfa1(gp) and ATP shifts the distribution of states towards channel opening, facilitating the accommodation of RNA. After completion of the loading process, the channel remains firmly closed during successive cycles of ATP hydrolysis, ensuring stable interaction with the RNA and processive translocation. Without Pfa1(gp), it remains predominantly closed preventing efficient RNA loading. Our data reveal how the ligands of Prp43 regulate the structural dynamics of the RNA binding channel controlling the initial binding of RNA.

Dynamic interplay between target search and recognition for a Type I CRISPR-Cas system

CRISPR-Cas effector complexes enable the defense against foreign nucleic acids and have recently been exploited as molecular tools for precise genome editing at a target locus. To bind and cleave their target, the CRISPR-Cas effectors have to interrogate the entire genome for the presence of a matching sequence. Here we dissect the target search and recognition process of the Type I CRISPR-Cas complex Cascade by simultaneously monitoring DNA binding and R-loop formation by the complex. We directly quantify the effect of DNA supercoiling on the target recognition probability and demonstrate that Cascade uses facilitated diffusion for its target search. We show that target search and target recognition are tightly linked and that DNA supercoiling and limited 1D diffusion need to be considered when understanding target recognition and target search by CRISPR-Cas enzymes and engineering more efficient and precise variants.

Recent Advances in Single-Molecule Tracking and Imaging Techniques

Since the early 1990s, single-molecule detection in solution at room temperature has enabled direct observation of single biomolecules at work in real time and under physiological conditions, providing insights into complex biological systems that the traditional ensemble methods cannot offer. In particular, recent advances in single-molecule tracking techniques allow researchers to follow individual biomolecules in their native environments for a timescale of seconds to minutes, revealing not only the distinct pathways these biomolecules take for downstream signaling but also their roles in supporting life. In this review, we discuss various single-molecule tracking and imaging techniques developed to date, with an emphasis on advanced three-dimensional (3D) tracking systems that not only achieve ultrahigh spatiotemporal resolution but also provide sufficient working depths suitable for tracking single molecules in 3D tissue models. We then summarize the observables that can be extracted from the trajectory data. Methods to perform single-molecule clustering analysis and future directions are also discussed.

Universal features of Nsp1-mediated translational shutdown by coronaviruses

Nonstructural protein 1 (Nsp1) produced by coronaviruses shuts down host protein synthesis in infected cells. The C-terminal domain of SARS-CoV-2 Nsp1 was shown to bind to the small ribosomal subunit to inhibit translation, but it is not clear whether this mechanism is broadly used by coronaviruses, whether the N-terminal domain of Nsp1 binds the ribosome, or how Nsp1 specifically permits translation of viral mRNAs. Here, we investigated Nsp1 from three representative Betacoronaviruses - SARS-CoV-2, MERS-CoV, and Bat-Hp-CoV - using structural, biophysical, and biochemical assays. We revealed a conserved mechanism of host translational shutdown across the three coronaviruses. We further demonstrated that the N-terminal domain of Bat-Hp-CoV Nsp1 binds to the decoding center of the 40S subunit, where it would prevent mRNA and eIF1A binding. Structure-based biochemical experiments identified a conserved role of these inhibitory interactions in all three coronaviruses and showed that the same regions of Nsp1 are responsible for the preferential translation of viral mRNAs. Our results provide a mechanistic framework to understand how Betacoronaviruses overcome translational inhibition to produce viral proteins.

Observation of structural switch in nascent SAM-VI riboswitch during transcription at single-nucleotide and single-molecule resolution

Growing RNAs fold differently as they are transcribed, which modulates their finally adopted structures. Riboswitches regulate gene expression by structural change, which are sensitive to co-transcriptionally structural biology. Here we develop a strategy to track the structural change of RNAs during transcription at single-nucleotide and single-molecule resolution and use it to monitor individual transcripts of the SAM-VI riboswitch (riboSAM) as transcription proceeds, observing co-existence of five states in riboSAM. We report a bifurcated helix in one newly identified state from NMR and single-molecule FRET (smFRET) results, and its presence directs the translation inhibition in our cellular translation experiments. A model is proposed to illustrate the distinct switch patterns and gene-regulatory outcome of riboSAM when SAM is present or absent. Our strategy enables the precise mapping of RNAs' conformational landscape during transcription, and may combine with detection methods other than smFRET for structural studies of RNAs in general.

Basic helix-loop-helix pioneer factors interact with the histone octamer to invade nucleosomes and generate nucleosome-depleted regions

Nucleosomes drastically limit transcription factor (TF) occupancy, while pioneer transcription factors (PFs) somehow circumvent this nucleosome barrier. In this study, we compare nucleosome binding of two conserved S. cerevisiae basic helix-loop-helix (bHLH) TFs, Cbf1 and Pho4. A cryo-EM structure of Cbf1 in complex with the nucleosome reveals that the Cbf1 HLH region can electrostatically interact with exposed histone residues within a partially unwrapped nucleosome. Single-molecule fluorescence studies show that the Cbf1 HLH region facilitates efficient nucleosome invasion by slowing its dissociation rate relative to DNA through interactions with histones, whereas the Pho4 HLH region does not. In vivo studies show that this enhanced binding provided by the Cbf1 HLH region enables nucleosome invasion and ensuing repositioning. These structural, single-molecule, and in vivo studies reveal the mechanistic basis of dissociation rate compensation by PFs and how this translates to facilitating chromatin opening inside cells.

Optimal inference of molecular interaction dynamics in FRET microscopy

Intensity-based time-lapse fluorescence resonance energy transfer (FRET) microscopy has been a major tool for investigating cellular processes, converting otherwise unobservable molecular interactions into fluorescence time series. However, inferring the molecular interaction dynamics from the observables remains a challenging inverse problem, particularly when measurement noise and photobleaching are nonnegligible-a common situation in single-cell analysis. The conventional approach is to process the time-series data algebraically, but such methods inevitably accumulate the measurement noise and reduce the signal-to-noise ratio (SNR), limiting the scope of FRET microscopy. Here, we introduce an alternative probabilistic approach, B-FRET, generally applicable to standard 3-cube FRET-imaging data. Based on Bayesian filtering theory, B-FRET implements a statistically optimal way to infer molecular interactions and thus drastically improves the SNR. We validate B-FRET using simulated data and then apply it to real data, including the notoriously noisy in vivo FRET time series from individual bacterial cells to reveal signaling dynamics otherwise hidden in the noise.

Reliability and accuracy of single-molecule FRET studies for characterization of structural dynamics and distances in proteins

Single-molecule Förster-resonance energy transfer (smFRET) experiments allow the study of biomolecular structure and dynamics in vitro and in vivo. We performed an international blind study involving 19 laboratories to assess the uncertainty of FRET experiments for proteins with respect to the measured FRET efficiency histograms, determination of distances, and the detection and quantification of structural dynamics. Using two protein systems with distinct conformational changes and dynamics, we obtained an uncertainty of the FRET efficiency ≤0.06, corresponding to an interdye distance precision of ≤2 Å and accuracy of ≤5 Å. We further discuss the limits for detecting fluctuations in this distance range and how to identify dye perturbations. Our work demonstrates the ability of smFRET experiments to simultaneously measure distances and avoid the averaging of conformational dynamics for realistic protein systems, highlighting its importance in the expanding toolbox of integrative structural biology.

Single-photon smFRET. I: Theory and conceptual basis

We present a unified conceptual framework and the associated software package for single-molecule Förster resonance energy transfer (smFRET) analysis from single-photon arrivals leveraging Bayesian nonparametrics, BNP-FRET. This unified framework addresses the following key physical complexities of a single-photon smFRET experiment, including: 1) fluorophore photophysics; 2) continuous time kinetics of the labeled system with large timescale separations between photophysical phenomena such as excited photophysical state lifetimes and events such as transition between system states; 3) unavoidable detector artefacts; 4) background emissions; 5) unknown number of system states; and 6) both continuous and pulsed illumination. These physical features necessarily demand a novel framework that extends beyond existing tools. In particular, the theory naturally brings us to a hidden Markov model with a second-order structure and Bayesian nonparametrics on account of items 1, 2, and 5 on the list. In the second and third companion articles, we discuss the direct effects of these key complexities on the inference of parameters for continuous and pulsed illumination, respectively.

Single-photon smFRET: II. Application to continuous illumination

Here we adapt the Bayesian nonparametrics (BNP) framework presented in the first companion article to analyze kinetics from single-photon, single-molecule Förster resonance energy transfer (smFRET) traces generated under continuous illumination. Using our sampler, BNP-FRET, we learn the escape rates and the number of system states given a photon trace. We benchmark our method by analyzing a range of synthetic and experimental data. Particularly, we apply our method to simultaneously learn the number of system states and the corresponding kinetics for intrinsically disordered proteins using two-color FRET under varying chemical conditions. Moreover, using synthetic data, we show that our method can deduce the number of system states even when kinetics occur at timescales of interphoton intervals.

CLASP2 facilitates dynamic actin filament organization along the microtubule lattice

Coordination between the microtubule and actin networks is essential for cell motility, neuronal growth cone guidance, and wound healing. Members of the CLASP (cytoplasmic linker-associated protein) family of proteins have been implicated in the cytoskeletal cross-talk between microtubules and actin networks; however, the molecular mechanisms underlying the role of CLASP in cytoskeletal coordination are unclear. Here, we investigate CLASP2α's cross-linking function with microtubules and F-actin. Our results demonstrate that CLASP2α cross-links F-actin to the microtubule lattice in vitro. We find that the cross-linking ability is retained by L-TOG2-S, a minimal construct containing the TOG2 domain and serine-arginine-rich region of CLASP2α. Furthermore, CLASP2α promotes the accumulation of multiple actin filaments along the microtubule, supporting up to 11 F-actin landing events on a single microtubule lattice region. CLASP2α also facilitates the dynamic organization of polymerizing actin filaments templated by the microtubule network, with F-actin forming bridges between individual microtubules. Finally, we find that depletion of CLASPs in vascular smooth muscle cells results in disorganized actin fibers and reduced coalignment of actin fibers with microtubules, suggesting that CLASP and microtubules contribute to higher-order actin structures. Taken together, our results indicate that CLASP2α can directly cross-link F-actin to microtubules and that this microtubule-CLASP-actin interaction may influence overall cytoskeletal organization in cells.

The nucleosome unwrapping free energy landscape defines distinct regions of transcription factor accessibility and kinetics

Transcription factors (TF) require access to target sites within nucleosomes to initiate transcription. The target site position within the nucleosome significantly influences TF occupancy, but how is not quantitatively understood. Using ensemble and single-molecule fluorescence measurements, we investigated the targeting and occupancy of the transcription factor, Gal4, at different positions within the nucleosome. We observe a dramatic decrease in TF occupancy to sites extending past 30 base pairs (bp) into the nucleosome which cannot be explained by changes in the TF dissociation rate or binding site orientation. Instead, the nucleosome unwrapping free energy landscape is the primary determinant of Gal4 occupancy by reducing the Gal4 binding rate. The unwrapping free energy landscape defines two distinct regions of accessibility and kinetics with a boundary at 30 bp into the nucleosome where the inner region is over 100-fold less accessible. The Gal4 binding rate in the inner region no longer depends on its concentration because it is limited by the nucleosome unwrapping rate, while the frequency of nucleosome rewrapping decreases because Gal4 exchanges multiple times before the nucleosome rewraps. Our findings highlight the importance of the nucleosome unwrapping free energy landscape on TF occupancy and dynamics that ultimately influences transcription initiation.

Inhibited complete folding of consecutive human telomeric G-quadruplexes

Noncanonical DNA structures, termed G-quadruplexes, are present in human genomic DNA and are important elements in many DNA metabolic processes. Multiple sites in the human genome have G-rich DNA stretches able to support formation of several consecutive G-quadruplexes. One of those sites is the telomeric overhang region that has multiple repeats of TTAGGG and is tightly associated with both cancer and aging. We investigated the folding of consecutive G-quadruplexes in both potassium- and sodium-containing solutions using single-molecule FRET spectroscopy, circular dichroism, thermal melting and molecular dynamics simulations. Our observations show coexistence of partially and fully folded DNA, the latter consisting of consecutive G-quadruplexes. Following the folding process over hours in sodium-containing buffers revealed fast G-quadruplex folding but slow establishment of thermodynamic equilibrium. We find that full consecutive G-quadruplex formation is inhibited by the many DNA structures randomly nucleating on the DNA, some of which are off-path conformations that need to unfold to allow full folding. Our study allows describing consecutive G-quadruplex formation in both nonequilibrium and equilibrium conditions by a unified picture, where, due to the many possible DNA conformations, full folding with consecutive G-quadruplexes as beads on a string is not necessarily achieved.

Modulation of protein-DNA binding reveals mechanisms of spatiotemporal gene control in early Drosophila embryos

It is well known that enhancers regulate the spatiotemporal expression of their target genes by recruiting transcription factors (TFs) to the cognate binding sites in the region. However, the role of multiple binding sites for the same TFs and their specific spatial arrangement in determining the overall competency of the enhancer has yet to be fully understood. In this study, we utilized the MS2-MCP live imaging technique to quantitatively analyze the regulatory logic of the snail distal enhancer in early Drosophila embryos. Through systematic modulation of Dorsal and Twist binding motifs in this enhancer, we found that a mutation in any one of these binding sites causes a drastic reduction in transcriptional amplitude, resulting in a reduction in total mRNA production of the target gene. We provide evidence of synergy, such that multiple binding sites with moderate affinities cooperatively recruit more TFs to drive stronger transcriptional activity than a single site. Moreover, a Hidden Markov-based stochastic model of transcription reveals that embryos with mutated binding sites have a higher probability of returning to the inactive promoter state. We propose that TF-DNA binding regulates spatial and temporal gene expression and drives robust pattern formation by modulating transcriptional kinetics and tuning bursting rates.

Probing Transient Riboswitch Structures via Single Molecule Accessibility Analysis

Riboswitches are a class of RNA motifs in the untranslated regions of bacterial messenger RNAs (mRNAs) that can adopt different conformations to regulate gene expression. The binding of specific small molecule or ion ligands, or other RNAs, influences the conformation the riboswitch adopts. Single Molecule Kinetic Analysis of RNA Transient Structure (SiM-KARTS) offers an approach for probing this structural isomerization, or conformational switching, at the level of single mRNA molecules. SiM-KARTS utilizes fluorescently labeled, short, sequence-complementary DNA or RNA oligonucleotide probes that transiently access a specific RNA conformation over another. Binding and dissociation to a surface-immobilized target RNA of arbitrary length are monitored by Total Internal Reflection Fluorescence Microscopy (TIRFM) and quantitatively analyzed, via spike train and burst detection, to elucidate the rate constants of isomerization, revealing mechanistic insights into riboswitching.

Integration Host Factor Binds DNA Holliday Junctions

Integration host factor (IHF) is a nucleoid-associated protein involved in DNA packaging, integration of viral DNA and recombination. IHF binds with nanomolar affinity to duplex DNA containing a 13 bp consensus sequence, inducing a bend of ~160° upon binding. We determined that IHF binds to DNA Four-way or Holliday junctions (HJ) with high affinity regardless of the presence of the consensus sequence, signifying a structure-based mechanism of recognition. Junctions, important intermediates in DNA repair and homologous recombination, are dynamic and can adopt either an open or stacked conformation, where the open conformation facilitates branch migration and strand exchange. Using ensemble and single molecule Förster resonance energy transfer (FRET) methods, we investigated IHF-induced changes in the population distribution of junction conformations and determined that IHF binding shifts the population to the open conformation. Further analysis of smFRET dynamics revealed that even in the presence of protein, the junctions remain dynamic as fast transitions are observed for the protein-bound open state. Protein binding alters junction conformational dynamics, as cross correlation analyses reveal the protein slows the transition rate at 1 mM Mg²⁺ but accelerates the transition rate at 10 mM Mg²⁺. Stopped flow kinetic experiments provide evidence for two binding steps, a rapid, initial binding step followed by a slower step potentially associated with a conformational change. These measurements also confirm that the protein remains bound to the junction during the conformer transitions and further suggest that the protein forms a partially dissociated state that allows junction arms to be dynamic. These findings, which demonstrate that IHF binds HJs with high affinity and stabilizes junctions in the open conformation, suggest that IHF may play multiple roles in the processes of integration and recombination in addition to stabilizing bacterial biofilms.

A scanning-to-incision switch in TFIIH-XPG induced by DNA damage licenses nucleotide excision repair

Nucleotide excision repair (NER) is critical for removing bulky DNA base lesions and avoiding diseases. NER couples lesion recognition by XPC to strand separation by XPB and XPD ATPases, followed by lesion excision by XPF and XPG nucleases. Here, we describe key regulatory mechanisms and roles of XPG for and beyond its cleavage activity. Strikingly, by combing single-molecule imaging and bulk cleavage assays, we found that XPG binding to the 7-subunit TFIIH core (coreTFIIH) stimulates coreTFIIH-dependent double-strand (ds)DNA unwinding 10-fold, and XPG-dependent DNA cleavage by up to 700-fold. Simultaneous monitoring of rates for coreTFIIH single-stranded (ss)DNA translocation and dsDNA unwinding showed XPG acts by switching ssDNA translocation to dsDNA unwinding as a likely committed step. Pertinent to the NER pathway regulation, XPG incision activity is suppressed during coreTFIIH translocation on DNA but is licensed when coreTFIIH stalls at the lesion or when ATP hydrolysis is blocked. Moreover, ≥15 nucleotides of 5'-ssDNA is a prerequisite for efficient translocation and incision. Our results unveil a paired coordination mechanism in which key lesion scanning and DNA incision steps are sequentially coordinated, and damaged patch removal is only licensed after generation of ≥15 nucleotides of 5'-ssDNA, ensuring the correct ssDNA bubble size before cleavage.

A quantitative model for the dynamics of target recognition and off-target rejection by the CRISPR-Cas Cascade complex

CRISPR-Cas effector complexes recognise nucleic acid targets by base pairing with their crRNA which enables easy re-programming of the target specificity in rapidly emerging genome engineering applications. However, undesired recognition of off-targets, that are only partially complementary to the crRNA, occurs frequently and represents a severe limitation of the technique. Off-targeting lacks comprehensive quantitative understanding and prediction. Here, we present a detailed analysis of the target recognition dynamics by the Cascade surveillance complex on a set of mismatched DNA targets using single-molecule supercoiling experiments. We demonstrate that the observed dynamics can be quantitatively modelled as a random walk over the length of the crRNA-DNA hybrid using a minimal set of parameters. The model accurately describes the recognition of targets with single and double mutations providing an important basis for quantitative off-target predictions. Importantly the model intrinsically accounts for observed bias regarding the position and the proximity between mutations and reveals that the seed length for the initiation of target recognition is controlled by DNA supercoiling rather than the Cascade structure.

Altered tRNA dynamics during translocation on slippery mRNA as determinant of spontaneous ribosome frameshifting

When reading consecutive mRNA codons, ribosomes move by exactly one triplet at a time to synthesize a correct protein. Some mRNA tracks, called slippery sequences, are prone to ribosomal frameshifting, because the same tRNA can read both 0- and –1-frame codon. Using smFRET we show that during EF-G-catalyzed translocation on slippery sequences a fraction of ribosomes spontaneously switches from rapid, accurate translation to a slow, frameshifting-prone translocation mode where the movements of peptidyl- and deacylated tRNA become uncoupled. While deacylated tRNA translocates rapidly, pept-tRNA continues to fluctuate between chimeric and posttranslocation states, which slows down the re-locking of the small ribosomal subunit head domain. After rapid release of deacylated tRNA, pept-tRNA gains unconstrained access to the –1-frame triplet, resulting in slippage followed by recruitment of the –1-frame aa-tRNA into the A site. Our data show how altered choreography of tRNA and ribosome movements reduces the translation fidelity of ribosomes translocating in a slow mode.

Entropy Driving the Mg2+-Induced Folding of TPP Riboswitch RNA

Mg²⁺ is well known to facilitate the structural folding of RNA. However, the thermodynamic and dynamic roles of Mg²⁺ in RNA folding remain elusive. Here, we exploit single-molecule fluorescence resonance energy transfer (smFRET) and isothermal titration calorimetry (ITC) to study the mechanism of Mg²⁺ in facilitating the folding of thiamine pyrophosphate (TPP) riboswitch RNA. The results of smFRET identify that the presence of Mg²⁺ compacts the RNA and enlarges the conformational dispersity among individual RNA molecules, resulting in a large gain of entropy. The compact yet flexible conformations triggered by Mg²⁺ may help the riboswitch recognize its specific ligand and further fold. This is supported by the ITC experiments, in which the Mg²⁺-induced RNA folding is driven by entropy (ΔS) instead of enthalpy (ΔH). Our results complement the understanding of the Mg²⁺-induced RNA folding. The strategy developed in this work can be used to model other RNAs' folding under different conditions.