A blind benchmark of analysis tools to infer kinetic rate constants from single-molecule FRET trajectories

Single-molecule FRET (smFRET) is a versatile technique to study the dynamics and function of biomolecules since it makes nanoscale movements detectable as fluorescence signals. The powerful ability to infer quantitative kinetic information from smFRET data is, however, complicated by experimental limitations. Diverse analysis tools have been developed to overcome these hurdles but a systematic comparison is lacking. Here, we report the results of a blind benchmark study assessing eleven analysis tools used to infer kinetic rate constants from smFRET trajectories. We test them against simulated and experimental data containing the most prominent difficulties encountered in analyzing smFRET experiments: different noise levels, varied model complexity, non-equilibrium dynamics, and kinetic heterogeneity. Our results highlight the current strengths and limitations in inferring kinetic information from smFRET trajectories. In addition, we formulate concrete recommendations and identify key targets for future developments, aimed to advance our understanding of biomolecular dynamics through quantitative experiment-derived models.

A blind benchmark of analysis tools to infer kinetic rate constants from single-molecule FRET trajectories

Single-molecule FRET (smFRET) is a versatile technique to study the dynamics and function of biomolecules since it makes nanoscale movements detectable as fluorescence signals. The powerful ability to infer quantitative kinetic information from smFRET data is, however, complicated by experimental limitations. Diverse analysis tools have been developed to overcome these hurdles but a systematic comparison is lacking. Here, we report the results of a blind benchmark study assessing eleven analysis tools used to infer kinetic rate constants from smFRET trajectories. We test them against simulated and experimental data containing the most prominent difficulties encountered in analyzing smFRET experiments: different noise levels, varied model complexity, non-equilibrium dynamics, and kinetic heterogeneity. Our results highlight the current strengths and limitations in inferring kinetic information from smFRET trajectories. In addition, we formulate concrete recommendations and identify key targets for future developments, aimed to advance our understanding of biomolecular dynamics through quantitative experiment-derived models.

Non-ergodicity of a globular protein extending beyond its functional timescale

Internal motions of folded proteins have been assumed to be ergodic, i.e., that the dynamics of a single protein molecule averaged over a very long time resembles that of an ensemble. Here, by performing single-molecule fluorescence resonance energy transfer (smFRET) experiments and molecular dynamics (MD) simulations of a multi-domain globular protein, cytoplasmic protein-tyrosine phosphatase (SHP2), we demonstrate that the functional inter-domain motion is observationally non-ergodic over the time spans 10⁻¹² to 10⁻⁷ s and 10⁻¹ to 10² s. The difference between observational non-ergodicity and simple non-convergence is discussed. In comparison, a single-strand DNA of similar size behaves ergodically with an energy landscape resembling a one-dimensional linear chain. The observed non-ergodicity results from the hierarchical connectivity of the high-dimensional energy landscape of the protein molecule. As the characteristic time for the protein to conduct its dephosphorylation function is ∼10 s, our findings suggest that, due to the non-ergodicity, individual, seemingly identical protein molecules can be dynamically and functionally different.

Internal motions of folded proteins have been assumed to be ergodic, i.e., that the dynamics of a single protein molecule averaged over a very long time resembles that of an ensemble.

Integrative structural dynamics probing of the conformational heterogeneity in synaptosomal-associated protein 25

SNAP-25 (synaptosomal-associated protein of 25 kDa) is a prototypical intrinsically disordered protein (IDP) that is unstructured by itself but forms coiled-coil helices in the SNARE complex. With high conformational heterogeneity, detailed structural dynamics of unbound SNAP-25 remain elusive. Here, we report an integrative method to probe the structural dynamics of SNAP-25 by combining replica-exchange discrete molecular dynamics (rxDMD) simulations and label-based experiments at ensemble and single-molecule levels. The rxDMD simulations systematically characterize the coil-to-molten globular transition and reconstruct structural ensemble consistent with prior ensemble experiments. Label-based experiments using Förster resonance energy transfer and double electron-electron resonance further probe the conformational dynamics of SNAP-25. Agreements between simulations and experiments under both ensemble and single-molecule conditions allow us to assign specific helix-coil transitions in SNAP-25 that occur in submillisecond timescales and potentially play a vital role in forming the SNARE complex. We expect that this integrative approach may help further our understanding of IDPs.

How Membrane Geometry Regulates Protein Sorting Independently of Mean Curvature

Biological membranes have distinct geometries that confer specific functions. However, the molecular mechanisms underlying the phenomenological geometry/function correlations remain elusive. We studied the effect of membrane geometry on the localization of membrane-bound proteins. Quantitative comparative experiments between the two most abundant cellular membrane geometries, spherical and cylindrical, revealed that geometry regulates the spatial segregation of proteins. The measured geometry-driven segregation reached 50-fold for membranes of the same mean curvature, demonstrating a crucial and hitherto unaccounted contribution by Gaussian curvature. Molecular-field theory calculations elucidated the underlying physical and molecular mechanisms. Our results reveal that distinct membrane geometries have specific physicochemical properties and thus establish a ubiquitous mechanistic foundation for unravelling the conserved correlations between biological function and membrane polymorphism.

Coordinated protein and DNA conformational changes govern mismatch repair initiation by MutS

MutS homologs identify base-pairing errors made in DNA during replication and initiate their repair. In the presence of adenosine triphosphate, MutS induces DNA bending upon mismatch recognition and subsequently undergoes conformational transitions that promote its interaction with MutL to signal repair. In the absence of MutL, these transitions lead to formation of a MutS mobile clamp that can move along the DNA. Previous single-molecule FRET (smFRET) studies characterized the dynamics of MutS DNA-binding domains during these transitions. Here, we use protein–DNA and DNA–DNA smFRET to monitor DNA conformational changes, and we use kinetic analyses to correlate DNA and protein conformational changes to one another and to the steps on the pathway to mobile clamp formation. The results reveal multiple sequential structural changes in both MutS and DNA, and they suggest that DNA dynamics play a critical role in the formation of the MutS mobile clamp. Taking these findings together with data from our previous studies, we propose a unified model of coordinated MutS and DNA conformational changes wherein initiation of mismatch repair is governed by a balance of DNA bending/unbending energetics and MutS conformational changes coupled to its nucleotide binding properties.

Spontaneous Switching among Conformational Ensembles in Intrinsically Disordered Proteins

The common conception of intrinsically disordered proteins (IDPs) is that they stochastically sample all possible configurations driven by thermal fluctuations. This is certainly true for many IDPs, which behave as swollen random coils that can be described using polymer models developed for homopolymers. However, the variability in interaction energy between different amino acid sequences provides the possibility that some configurations may be strongly preferred while others are forbidden. In compact globular IDPs, core hydration and packing density can vary between segments of the polypeptide chain leading to complex conformational dynamics. Here, we describe a growing number of proteins that appear intrinsically disordered by biochemical and bioinformatic characterization but switch between restricted regions of conformational space. In some cases, spontaneous switching between conformational ensembles was directly observed, but few methods can identify when an IDP is acting as a restricted chain. Such switching between disparate corners of conformational space could bias ligand binding and regulate the volume of IDPs acting as structural or entropic elements. Thus, mapping the accessible energy landscape and capturing dynamics across a wide range of timescales are essential to recognize when an IDP is acting as such a switch.

Single Molecule FRET: A Powerful Tool to Study Intrinsically Disordered Proteins

Intrinsically disordered proteins (IDPs) are often modeled using ideas from polymer physics that suggest they smoothly explore all corners of configuration space. Experimental verification of this random, dynamic behavior is difficult as random fluctuations of IDPs cannot be synchronized across an ensemble. Single molecule fluorescence (or Förster) resonance energy transfer (smFRET) is one of the few approaches that are sensitive to transient populations of sub-states within molecular ensembles. In some implementations, smFRET has sufficient time resolution to resolve transitions in IDP behaviors. Here we present experimental issues to consider when applying smFRET to study IDP configuration. We illustrate the power of applying smFRET to IDPs by discussing two cases in the literature of protein systems for which smFRET has successfully reported phosphorylation-induced modification (but not elimination) of the disordered properties that have been connected to impacts on the related biological function. The examples we discuss, PAGE4 and a disordered segment of the GluN2B subunit of the NMDA receptor, illustrate the great potential of smFRET to inform how IDP function can be regulated by controlling the detailed ensemble of disordered states within biological networks.

Motor-like DNA motion due to an ATP-hydrolyzing protein under nanoconfinement

We report that long double-stranded DNA confined to quasi-1D nanochannels undergoes superdiffusive motion under the action of the enzyme T4 DNA ligase in the presence of necessary co-factors. Inside the confined environment of the nanochannel, double-stranded DNA molecules stretch out due to self-avoiding interactions. In absence of a catalytically active enzyme, we see classical diffusion of the center of mass. However, cooperative interactions of proteins with the DNA can lead to directed motion of DNA molecules inside the nanochannel. Here we show directed motion in this configuration for three different proteins (T4 DNA ligase, MutS, E. coli DNA ligase) in the presence of their energetic co-factors (ATP, NAD⁺).

Single molecule studies of DNA mismatch repair

DNA mismatch repair, which involves is a widely conserved set of proteins, is essential to limit genetic drift in all organisms. The same system of proteins plays key roles in many cancer related cellular transactions in humans. Although the basic process has been reconstituted in vitro using purified components, many fundamental aspects of DNA mismatch repair remain hidden due in part to the complexity and transient nature of the interactions between the mismatch repair proteins and DNA substrates. Single molecule methods offer the capability to uncover these transient but complex interactions and allow novel insights into mechanisms that underlie DNA mismatch repair. In this review, we discuss applications of single molecule methodology including electron microscopy, atomic force microscopy, particle tracking, FRET, and optical trapping to studies of DNA mismatch repair. These studies have led to formulation of mechanistic models of how proteins identify single base mismatches in the vast background of matched DNA and signal for their repair.

Dynamics of MutS-mismatched DNA complexes are predictive of their repair phenotypes

MutS recognizes base–base mismatches and base insertions/deletions (IDLs) in newly replicated DNA. Specific interactions between MutS and these errors trigger a cascade of protein–protein interactions that ultimately lead to their repair. The inability to explain why different DNA errors are repaired with widely varying efficiencies in vivo remains an outstanding example of our limited knowledge of this process. Here, we present single-molecule Förster resonance energy transfer measurements of the DNA bending dynamics induced by Thermus aquaticus MutS and the E41A mutant of MutS, which is known to have error specific deficiencies in signaling repair. We compared three DNA mismatches/IDLs (T-bulge, GT, and CC) with repair efficiencies ranging from high to low. We identify three dominant DNA bending states [slightly bent/unbent (U), intermediately bent (I), and significantly bent (B)] and find that the kinetics of interconverting among states varies widely for different complexes. The increased stability of MutS–mismatch/IDL complexes is associated with stabilization of U and lowering of the B to U transition barrier. Destabilization of U is always accompanied by a destabilization of B, supporting the suggestion that B is a “required” precursor to U. Comparison of MutS and MutS-E41A dynamics on GT and the T-bulge suggests that hydrogen bonding to MutS facilitates the changes in base–base hydrogen bonding that are required to achieve the U state, which has been implicated in repair signaling. Taken together with repair propensities, our data suggest that the bending kinetics of MutS–mismatched DNA complexes may control the entry into functional pathways for downstream signaling of repair.

Cancer/testis antigen PAGE4, a regulator of c-Jun transactivation, is phosphorylated by homeodomain-interacting protein kinase 1, a component of the stress-response pathway

Prostate-associated gene 4 (PAGE4) is a cancer/testis antigen that is typically restricted to the testicular germ cells but is aberrantly expressed in cancer. Furthermore, PAGE4 is developmentally regulated with dynamic expression patterns in the developing prostate and is also a stress-response protein that is upregulated in response to cellular stress. PAGE4 interacts with c-Jun, which is activated by the stress-response kinase JNK1, and plays an important role in the development and pathology of the prostate gland. Here, we have identified homeodomain-interacting protein kinase 1 (HIPK1), also a component of the stress-response pathway, as a kinase that phosphorylates PAGE4 at T51. We show that phosphorylation of PAGE4 is critical for its transcriptional activity since mutating this T residue abolishes its ability to potentiate c-Jun transactivation. In vitro single molecule FRET indicates phosphorylation results in compaction of (still) intrinsically disordered PAGE4. Interestingly, however, while our previous observations indicated that the wild-type nonphosphorylated PAGE4 protein interacted with c-Jun [RajagopalanK. et al. (2014) Biochim, Biophys. Acta1842, 154−16324263171], here we show that phosphorylation of PAGE4 weakens its interaction with c-Jun in vitro. These data suggest that phosphorylation induces conformational changes in natively disordered PAGE4 resulting in its decreased affinity for c-Jun to promote interaction of c-Jun with another, unidentified, partner. Alternatively, phosphorylated PAGE4 may induce transcription of a novel partner, which then potentiates c-Jun transactivation. Regardless, the present results clearly implicate PAGE4 as a component of the stress-response pathway and uncover a novel link between components of this pathway and prostatic development and disease.

The Stress-response protein prostate-associated gene 4, interacts with c-Jun and potentiates its transactivation

The Cancer/Testis Antigen (CTA), Prostate-associated Gene 4 (PAGE4), is a stress-response protein that is upregulated in prostate cancer (PCa) especially in precursor lesions that result from inflammatory stress. In cells under stress, translocation of PAGE4 to mitochondria increases while production of reactive oxygen species decreases. Furthermore, PAGE4 is also upregulated in human fetal prostate, underscoring its potential role in development. However, the proteins that interact with PAGE4 and the mechanisms underlying its pleiotropic functions in prostatic development and disease remain unknown. Here, we identified c-Jun as a PAGE4 interacting partner. We show that both PAGE4 and c-Jun are overexpressed in the human fetal prostate; and in cell-based assays, PAGE4 robustly potentiates c-Jun transactivation. Single-molecule Förster resonance energy transfer experiments indicate that upon binding to c-Jun, PAGE4 undergoes conformational changes. However, no interaction is observed in presence of BSA or unilamellar vesicles containing the mitochondrial inner membrane diphosphatidylglycerol lipid marker cardiolipin. Together, our data indicate that PAGE4 specifically interacts with c-Jun and that, conformational dynamics may account for its observed pleiotropic functions. To our knowledge, this is the first report demonstrating crosstalk between a CTA and a proto-oncogene. Disrupting PAGE4/c-Jun interactions using small molecules may represent a novel therapeutic strategy for PCa.

Beyond the random coil: stochastic conformational switching in intrinsically disordered proteins

Intrinsically disordered proteins (IDPs) participate in critical cellular functions that exploit the flexibility and rapid conformational fluctuations of their native state. Limited information about the native state of IDPs can be gained by the averaging over many heterogeneous molecules that is unavoidable in ensemble approaches. We used single molecule fluorescence to characterize native state conformational dynamics in five synaptic proteins confirmed to be disordered by other techniques. For three of the proteins, SNAP-25, synaptobrevin and complexin, their conformational dynamics could be described with a simple semiflexible polymer model. Surprisingly, two proteins, neuroligin and the NMDAR-2B glutamate receptor, were observed to stochastically switch among distinct conformational states despite the fact that they appeared intrinsically disordered by other measures. The hop-like intramolecular diffusion found in these proteins is suggested to define a class of functionality previously unrecognized for IDPs.

Three-dimensional molecular modeling with single molecule FRET

Single molecule fluorescence energy transfer experiments enable investigations of macromolecular conformation and folding by the introduction of fluorescent dyes at specific sites in the macromolecule. Multiple such experiments can be performed with different labeling site combinations in order to map complex conformational changes or interactions between multiple molecules. Distances that are derived from such experiments can be used for determination of the fluorophore positions by triangulation. When combined with a known structure of the macromolecule(s) to which the fluorophores are attached, a three-dimensional model of the system can be determined. However, care has to be taken to properly derive distance from fluorescence energy transfer efficiency and to recognize the systematic or random errors for this relationship. Here we review the experimental and computational methods used for three-dimensional modeling based on single molecule fluorescence resonance transfer, and describe recent progress in pushing the limits of this approach to macromolecular complexes.

Single-molecule FRET-derived model of the synaptotagmin 1-SNARE fusion complex

Synchronous neurotransmission is triggered when Ca(2+) binds to synaptotagmin 1 (Syt1), a synaptic-vesicle protein that interacts with SNAREs and membranes. We used single-molecule fluorescence resonance energy transfer (FRET) between synaptotagmin's two C2 domains to determine that their conformation consists of multiple states with occasional transitions, consistent with domains in random relative motion. SNARE binding results in narrower intrasynaptotagmin FRET distributions and less frequent transitions between states. We obtained an experimentally determined model of the elusive Syt1-SNARE complex using a multibody docking approach with 34 FRET-derived distances as restraints. The Ca(2+)-binding loops point away from the SNARE complex, so they may interact with the same membrane. The loop arrangement is similar to that of the crystal structure of SNARE-induced Ca(2+)-bound Syt3, suggesting a common mechanism by which the interaction between synaptotagmins and SNAREs aids in Ca(2+)-triggered fusion.

Detecting the conformation of individual proteins in live cells

We combined single molecule fluorescence resonance energy transfer (smFRET) with single particle tracking in live cells to detect the in vivo conformation of individual proteins. We site-specifically labeled recombinant SNARE proteins with a FRET donor and acceptor before microinjecting them into cultured cells. We observed that individual proteins rapidly incorporated into folded complexes at the cell membrane, demonstrating the potential of this method to reveal dynamic interactions within cells.

Comment on "Mie scattering from a sonoluminescing bubble with high spatial and temporal resolution"

A key parameter underlying the existence of sonoluminescence (or SL) is the time dependence of the radius R(t) of the collapsing bubble from which SL originates. With regard to the use of light scattering to measure this quantity, we wish to note that we disagree with the statement of Gompf and Pecha-highly compressed water causes the minimum in scattered light to occur 700 ps before SL-and that this effect leads to an overestimate of the bubble wall velocity. We discuss potential artifacts in their experimental arrangement and reply to their criticisms of our experiments on Mie scattering.

Observation of bubble dynamics within luminescent cavitation clouds: Sonoluminescence at the nano-scale

Measurements of acoustically driven cavitation luminescence indicate that this phenomenon is robust over a huge parameter space ranging from 10 kHz to >10 MHz. The minimum bubble radius achieved is an upper bound for the size of the light-emitting region and ranges from about 1 microm at 15 kHz to tens of nm at 11 MHz. Although lines can be discerned in the spectra of some cavitation clouds, they sit on top of a broadband continuum which can have greater spectral density in the ultraviolet than is observed for resonantly driven sonoluminescence from a single bubble.

Time correlated single photon mie scattering from a sonoluminescing bubble

Application of time correlated single photon counting to pulsed Mie scattering enables one to resolve changes in light scattering to better than 50 ps. This technique is applied to the highly nonlinear motion of a sonoluminescing bubble. Physical processes, such as outgoing shock wave emission, that limit the interpretation of the data are measured with a streak camera and microscopy. Shock speeds about 6 km/s have been observed.