Quantifying internal friction in unfolded and intrinsically disordered proteins with single-molecule spectroscopy

Memory Effects Explain the Fractional Viscosity Dependence of Rates Associated with Internal Friction: Simple Models and Applications to Butane Dihedral Rotation

Barrier-crossing rates of biophysical processes, ranging from simple conformational changes to protein folding, often deviate from the Kramers prediction of an inverse viscosity dependence. In many recent studies, this has been attributed to the presence of internal friction within the system. Previously, we showed that memory-dependent friction arising from the nonequilibrium solvation of a single particle crossing a smooth one-dimensional barrier can also cause such a deviation and be misinterpreted as internal friction. Here we introduce a simple diatom model and show that even in the absence of explicit solvent, internal memory effects arise due to coupling of the reaction coordinate motion with frictionally orthogonal degrees of freedom. This results in a fractional viscosity dependence and a deviation from Kramers' theory, typically attributed to the presence of internal friction. This model therefore mimics several biological processes where a local conformational change of a biomolecule is often influenced by its surroundings. This gives rise to an apparent "internal friction" commonly measured in terms of empirical fitting parameters α and σ. We propose a microscopic measure of this internal friction using Grote-Hynes theory which employs memory-dependent friction. We use butane to demonstrate the effect of coupling strength on the internal friction in realistic systems. This model therefore can serve the purpose of understanding internal friction in biological systems in terms of such coupling.

Structural dynamics of the intrinsically disordered linker region of cardiac troponin T

The cardiac troponin complex, composed of troponins I, T, and C, plays a central role in regulating the calcium-dependent interactions between myosin and the thin filament. Mutations in troponin can cause cardiomyopathies; however, it is still a major challenge to connect how changes in sequence affect troponin's function. Recent high-resolution structures of the thin filament revealed critical insights into the structure-function relationship of troponin, but there remain large, unresolved segments of troponin, including the troponin-T linker region that is a hotspot for cardiomyopathy mutations. This linker region is predicted to be intrinsically disordered, with behaviors that are not well described by traditional structural approaches; however, this proposal has not been experimentally verified. Here, we used a combination of single-molecule Förster resonance energy transfer (FRET), molecular dynamics simulations, and functional reconstitution assays to investigate the troponin-T linker region. We show that in the context of both isolated troponin and the fully regulated troponin complex, the linker behaves as a dynamic, intrinsically disordered region. This region undergoes polyampholyte expansion in the presence of high salt and distinct conformational changes during the assembly of the troponin complex. We also examine the ΔE160 hypertrophic cardiomyopathy mutation in the linker and demonstrate that it does not affect the conformational dynamics of the linker, rather it allosterically affects interactions with other troponin complex subunits, leading to increased molecular contractility. Taken together, our data clearly demonstrate the importance of disorder within the troponin-T linker and provide new insights into the molecular mechanisms driving the pathogenesis of cardiomyopathies.

Disordered regions of human eIF4B orchestrate a dynamic self-association landscape

Eukaryotic translation initiation factor eIF4B is required for efficient cap-dependent translation, it is overexpressed in cancer cells, and may influence stress granule formation. Due to the high degree of intrinsic disorder, eIF4B is rarely observed in cryo-EM structures of translation complexes and only ever by its single structured RNA recognition motif domain, leaving the molecular details of its large intrinsically disordered region (IDR) unknown. By integrating experiments and simulations we demonstrate that eIF4B IDR orchestrates and fine-tunes an intricate transition from monomers to a condensed phase, in which large-size dynamic oligomers form before mesoscopic phase separation. Single-molecule spectroscopy combined with molecular simulations enabled us to characterize the conformational ensembles and underlying intra- and intermolecular dynamics across the oligomerization transition. The observed sensitivity to ionic strength and molecular crowding in the self-association landscape suggests potential regulation of eIF4B nanoscopic and mesoscopic behaviors such as driven by protein modifications, binding partners or changes to the cellular environment.

Direct computations of viscoelastic moduli of biomolecular condensates

Biomolecular condensates are viscoelastic materials defined by time-dependent, sequence-specific complex shear moduli. Here, we show that viscoelastic moduli can be computed directly using a generalization of the Rouse model that leverages information regarding intra- and inter-chain contacts, which we extract from equilibrium configurations of lattice-based Metropolis Monte Carlo (MMC) simulations of phase separation. The key ingredient of the generalized Rouse model is a graph Laplacian that we compute from equilibrium MMC simulations. We compute two flavors of graph Laplacians, one based on a single-chain graph that accounts only for intra-chain contacts, and the other referred to as a collective graph that accounts for inter-chain interactions. Calculations based on the single-chain graph systematically overestimate the storage and loss moduli, whereas calculations based on the collective graph reproduce the measured moduli with greater fidelity. However, in the long time, low-frequency domain, a mixture of the two graphs proves to be most accurate. In line with the theory of Rouse and contrary to recent assertions, we find that a continuous distribution of relaxation times exists in condensates. The single crossover frequency between dominantly elastic vs dominantly viscous behaviors does not imply a single relaxation time. Instead, it is influenced by the totality of the relaxation modes. Hence, our analysis affirms that viscoelastic fluid-like condensates are best described as generalized Maxwell fluids. Finally, we show that the complex shear moduli can be used to solve an inverse problem to obtain the relaxation time spectra that underlie the dynamics within condensates. This is of practical importance given advancements in passive and active microrheology measurements of condensate viscoelasticity.

Room-Temperature Solution Fluorescence Excitation Correlation Spectroscopy

Single-molecule fluorescence spectroscopy is a powerful tool for investigating the physical properties of individual molecules, yet elucidating the fast fluctuation dynamics of freely diffusing single molecules in solution at room temperature, where a variety of chemical and biological processes occur, remains challenging. In this study, we report on fluorescence excitation correlation spectroscopy of room-temperature solutions, which enables the study of spontaneous fluctuation of the excitation spectrum with microsecond time resolution. By employing Fourier transform spectroscopy with broadband femtosecond pulses and time-correlated single-photon counting, we achieved fluorescence excitation spectroscopy of a room-temperature solution at the single-molecule level. Building upon this single-molecule measurement, we obtained an excitation wavelength-resolved fluorescence autocorrelation function in the microsecond to millisecond range, demonstrating the potential of this method to elucidate fast, spontaneous, time-dependent changes of excitation spectra in statistically equilibrated systems. With further development, this method will allow the study of spectral exchange associated with transitions between sub-ensembles of solution-phase molecules with unprecedented time resolution.

Mesoscale properties of biomolecular condensates emerging from protein chain dynamics

Biomolecular condensates form by phase separation of biological polymers. The cellular functions of the resulting membraneless organelles are closely linked to their physical properties over a wide range of length- and timescales: From the nanosecond dynamics of individual molecules and their interactions, to the microsecond translational diffusion of molecules in the condensates, to their viscoelastic properties at the mesoscopic scale. However, it has remained unclear how to quantitatively link these properties across scales. Here we address this question by combining single-molecule fluorescence, correlation spectroscopy, microrheology, and large-scale molecular dynamics simulations on different condensates that are formed by complex coacervation and span about two orders of magnitude in viscosity and their dynamics at the molecular scale. Remarkably, we find that the absolute timescale of protein chain dynamics in the dense phases can be quantitatively and accurately related to translational diffusion and condensate viscosities by Rouse theory of polymer solutions including entanglement. The simulations indicate that the observed wide range of dynamics arises from different contact lifetimes between amino acid residues, which in the mean-field description of the polymer model cause differences in the friction acting on the chains. These results suggest that remarkably simple physical principles can relate the mesoscale properties of biomolecular condensates to their dynamics at the nanoscale.

Nanosecond chain dynamics of single-stranded nucleic acids

The conformational dynamics of single-stranded nucleic acids are fundamental for nucleic acid folding and function. However, their elementary chain dynamics have been difficult to resolve experimentally. Here we employ a combination of single-molecule Förster resonance energy transfer, nanosecond fluorescence correlation spectroscopy, and nanophotonic enhancement to determine the conformational ensembles and rapid chain dynamics of short single-stranded nucleic acids in solution. To interpret the experimental results in terms of end-to-end distance dynamics, we utilize the hierarchical chain growth approach, simple polymer models, and refinement with Bayesian inference to generate structural ensembles that closely align with the experimental data. The resulting chain reconfiguration times are exceedingly rapid, in the 10-ns range. Solvent viscosity-dependent measurements indicate that these dynamics of single-stranded nucleic acids exhibit negligible internal friction and are thus dominated by solvent friction. Our results provide a detailed view of the conformational distributions and rapid dynamics of single-stranded nucleic acids.

Direct computations of viscoelastic moduli of biomolecular condensates

In vitro facsimiles of biomolecular condensates are formed by different types of intrinsically disordered proteins including prion-like low complexity domains (PLCDs). PLCD condensates are viscoelastic materials defined by time-dependent, sequence-specific complex shear moduli. Here, we show that viscoelastic moduli can be computed directly using a generalization of the Rouse model and information regarding intra- and inter-chain contacts that is extracted from equilibrium configurations of lattice-based Metropolis Monte Carlo (MMC) simulations. The key ingredient of the generalized Rouse model is the Zimm matrix that we compute from equilibrium MMC simulations. We compute two flavors of Zimm matrices, one referred to as the single-chain model that accounts only for intra-chain contacts, and the other referred to as a collective model, that accounts for inter-chain interactions. The single-chain model systematically overestimates the storage and loss moduli, whereas the collective model reproduces the measured moduli with greater fidelity. However, in the long time, low-frequency domain, a mixture of the two models proves to be most accurate. In line with the theory of Rouse, we find that a continuous distribution of relaxation times exists in condensates. The single crossover frequency between dominantly elastic versus dominantly viscous behaviors is influenced by the totality of the relaxation modes. Hence, our analysis suggests that viscoelastic fluid-like condensates are best described as generalized Maxwell fluids. Finally, we show that the complex shear moduli can be used to solve an inverse problem to obtain distributions of relaxation times that underlie the dynamics within condensates.

The role of memory-dependent friction and solvent viscosity in isomerization kinetics in viscogenic media

Molecular isomerization kinetics in liquid solvent depends on a complex interplay between the solvent friction acting on the molecule, internal dissipation effects (also known as internal friction), the viscosity of the solvent, and the dihedral free energy profile. Due to the absence of accurate techniques to directly evaluate isomerization friction, it has not been possible to explore these relationships in full. By combining extensive molecular dynamics simulations with friction memory-kernel extraction techniques we consider a variety of small, isomerising molecules under a range of different viscogenic conditions and directly evaluate the viscosity dependence of the friction acting on a rotating dihedral. We reveal that the influence of different viscogenic media on isomerization kinetics can be dramatically different, even when measured at the same viscosity. This is due to the dynamic solute-solvent coupling, mediated by time-dependent friction memory kernels. We also show that deviations from the linear dependence of isomerization rates on solvent viscosity, which are often simply attributed to internal friction effects, are due to the simultaneous violation of two fundamental relationships: the Stokes-Einstein relation and the overdamped Kramers prediction for the barrier-crossing rate, both of which require explicit knowledge of friction.

Individual transcription factors modulate both the micromovement of chromatin and its long-range structure

The control of eukaryotic gene expression is intimately connected to highly dynamic chromatin structures. Gene regulation relies on activator and repressor transcription factors (TFs) that induce local chromatin opening and closing. However, it is unclear how nucleus-wide chromatin organization responds dynamically to the activity of specific TFs. Here, we examined how two TFs with opposite effects on local chromatin accessibility modulate chromatin dynamics nucleus-wide. We combine high-resolution diffusion mapping and dense flow reconstruction and correlation in living cells to obtain an imaging-based, nanometer-scale analysis of local diffusion processes and long-range coordinated movements of both chromatin and TFs. We show that the expression of either an individual transcriptional activator (CDX2) or repressor (SIX6) with large numbers of binding sites increases chromatin mobility nucleus-wide, yet they induce opposite coherent chromatin motions at the micron scale. Hi-C analysis of higher-order chromatin structures shows that induction of the pioneer factor CDX2 leads both to changes in local chromatin interactions and the distribution of A and B compartments, thus relating the micromovement of chromatin with changes in compartmental structures. Given that inhibition of transcription initiation and elongation by RNA Pol II has a partial impact on the global chromatin dynamics induced by CDX2, we suggest that CDX2 overexpression alters chromatin structure dynamics both dependently and independently of transcription. Our biophysical analysis shows that sequence-specific TFs can influence chromatin structure on multiple architectural levels, arguing that local chromatin changes brought by TFs alter long-range chromatin mobility and its organization.

Effects of Proline on Internal Friction in Simulated Folding Dynamics of Several Alanine-Based α-Helical Peptides

We have studied in silico the effect of proline, a model cosolvent, on local and global friction coefficients in (un)folding of several typical alanine-based α-helical peptides. Local friction is related to dwell times of a single, ensemble-averaged hydrogen bond (HB) within each peptide. Global friction is related to energy dissipated in a series of configurational changes of each peptide experienced by increasing the number of HBs during folding. Both of these approaches are important in relation to future atomic force microscopic-based measurements of internal friction via force-clamp single-molecule force spectroscopy. Molecular dynamics (MD) simulations for six peptides, namely, ALA5, ALA8, ALA15, ALA21, (AAQAA)3, and H2N-GN(AAQAA)2G-COONH2, have been conducted at 2 and 5 M proline solutions in water. Using previously obtained MD data for these peptides in pure water as well as upgraded theoretical models, we obtained variations of local and global internal friction coefficients as a function of solution viscosity. The results showed the substantial role of proline in stabilizing the folded state and slowing the overall folding dynamics. Consequently, larger friction coefficients were obtained at larger viscosities. The local and global internal friction, i.e., respective, friction coefficients approximated to zero viscosity, was also obtained. The evolution of friction coefficients with viscosity was weakly dependent on the number of concurrent folding pathways but was rather dominated by a stabilizing effect of proline on the folded states. Obtained values of local and global internal friction showed qualitatively similar results and a clear dependency on the structure of the studied peptide.

Unveiling the photoluminescence dynamics of gold nanoclusters with fluorescence correlation spectroscopy

Gold nanoclusters (AuNCs) have captured significant interest for their photoluminescent properties; however, their rapid photodynamics remain elusive while probed by ensemble-averaging spectroscopy techniques. To address this challenge, we use fluorescence correlation spectroscopy (FCS) to uncover the photoluminescence dynamics of colloidal Au18(SG)14 nanoclusters. Our FCS analysis reveals the photoluminescence (PL) brightness per nanocluster, elucidating the impact of photoexcitation saturation and ligand interactions. Unlike DNA-encapsulated silver nanoclusters, their gold counterparts notably exhibit minimal blinking, with moderate amplitudes and 200 μs characteristic times. Our data also clearly reveal the occurrence of photon antibunching in the PL emission, showcasing the quantum nature of the PL process, with each AuNC acting as an individual quantum source. Using zero-mode waveguide nanoapertures, we achieve a 16-fold enhancement of the PL brightness of individual AuNCs. This constitutes an important enabling proof-of-concept for tailoring emission properties through nanophotonics. Overall, our study bridges the gap between ensemble-averaged techniques and single-molecule spectroscopy, offering new insights into AuNC photodynamics for biosensing and imaging applications.

Quantifying the energy landscape in weakly and strongly disordered frictional media

We investigate the "roughness" of the energy landscape of a system that diffuses in a heterogeneous medium with a random position-dependent friction coefficient α(x). This random friction acting on the system stems from spatial inhomogeneity in the surrounding medium and is modeled using the generalized Caldira-Leggett model. For a weakly disordered medium exhibiting a Gaussian random diffusivity D(x) = kBT/α(x) characterized by its average value ⟨D(x)⟩ and a pair-correlation function ⟨D(x1)D(x2)⟩, we find that the renormalized intrinsic diffusion coefficient is lower than the average one due to the fluctuations in diffusivity. The induced weak internal friction leads to increased roughness in the energy landscape. When applying this idea to diffusive motion in liquid water, the dissociation energy for a hydrogen bond gradually approaches experimental findings as fluctuation parameters increase. Conversely, for a strongly disordered medium (i.e., ultrafast-folding proteins), the energy landscape ranges from a few to a few kcal/mol, depending on the strength of the disorder. By fitting protein folding dynamics to the escape process from a metastable potential, the decreased escape rate conceptualizes the role of strong internal friction. Studying the energy landscape in complex systems is helpful because it has implications for the dynamics of biological, soft, and active matter systems.

Electrostatics and hydrophobicity in the dynamics of intrinsically disordered proteins

Internal friction is a major contribution to the dynamics of intrinsically disordered proteins (IDPs). Yet, the molecular origin of internal friction has so far been elusive. Here, we investigate whether attractive electrostatic interactions in IDPs modulate internal friction differently than the hydrophobic effect. To this end, we used nanosecond fluorescence correlation spectroscopy (nsFCS) and single-molecule Förster resonance energy transfer (FRET) to quantify the conformation and dynamics of the disordered DNA-binding domains Myc, Max and Mad at different salt concentrations. We find that internal friction effects are stronger when the chain is compacted by electrostatic attractions compared to the hydrophobic effect. Although the effect is moderate, the results show that the heteropolymeric nature of IDPs is reflected in their dynamics.

Investigation of Rare Protein Conformational Transitions via Dissipation-Corrected Targeted Molecular Dynamics

To sample rare events, dissipation-corrected targeted molecular dynamics (dcTMD) applies a constant velocity constraint along a one-dimensional reaction coordinate s, which drives an atomistic system from an initial state into a target state. Employing a cumulant approximation of Jarzynski's identity, the free energy ΔG(s) is calculated from the mean external work and dissipated work of the process. By calculating the friction coefficient Γ(s) from the dissipated work, in a second step, the equilibrium dynamics of the process can be studied by propagating a Langevin equation. While so far dcTMD has been mostly applied to study the unbinding of protein-ligand complexes, here its applicability to rare conformational transitions within a protein and the prediction of their kinetics are investigated. As this typically requires the introduction of multiple collective variables {xj} = x, a theoretical framework is outlined to calculate the associated free energy ΔG(x) and friction Γ(x) from dcTMD simulations along coordinate s. Adopting the α-β transition of alanine dipeptide as well as the open-closed transition of T4 lysozyme as representative examples, the virtues and shortcomings of dcTMD to predict protein conformational transitions and the related kinetics are studied.

An automated single-molecule FRET platform for high-content, multiwell plate screening of biomolecular conformations and dynamics

Single-molecule FRET (smFRET) has become a versatile tool for probing the structure and functional dynamics of biomolecular systems, and is extensively used to address questions ranging from biomolecular folding to drug discovery. Confocal smFRET measurements are amongst the widely used smFRET assays and are typically performed in a single-well format. Thus, sampling of many experimental parameters is laborious and time consuming. To address this challenge, we extend here the capabilities of confocal smFRET beyond single-well measurements by integrating a multiwell plate functionality to allow for continuous and automated smFRET measurements. We demonstrate the broad applicability of the multiwell plate assay towards DNA hairpin dynamics, protein folding, competitive and cooperative protein-DNA interactions, and drug-discovery, revealing insights that would be very difficult to achieve with conventional single-well format measurements. For the adaptation into existing instrumentations, we provide a detailed guide and open-source acquisition and analysis software.

Local and Large-Scale Conformational Dynamics in Unfolded Proteins and IDPs. I. Effect of Solvent Viscosity and Macromolecular Crowding

Protein/solvent interactions largely influence protein dynamics, particularly motions in unfolded and intrinsically disordered proteins (IDPs). Here, we apply triplet-triplet energy transfer (TTET) to investigate the coupling of internal protein motions to solvent motions by determining the effect of solvent viscosity (η) and macromolecular crowding on the rate constants of loop formation (kc) in several unfolded polypeptide chains including IDPs. The results show that the viscosity dependence of loop formation depends on amino acid sequence, loop length, and co-solute size. Below a critical size (rc), co-solutes exert a maximum effect, indicating that under these conditions microviscosity experienced by chain motions matches macroviscosity of the solvent. rc depends on chain stiffness and reflects the length scale of the chain motions, i.e., it is related to the persistence length. Above rc, the effect of solvent viscosity decreases with increasing co-solute size. For co-solutes typically used to mimic cellular environments, a scaling of kc ∝ η-0.1 is observed, suggesting that dynamics in unfolded proteins are only marginally modulated in cells. The effect of solvent viscosity on kc in the small co-solute limit (below rc) increases with increasing chain length and chain flexibility. Formation of long and very flexible loops exhibits a kc ∝ η-1 viscosity dependence, indicating full solvent coupling. Shorter and less flexible loops show weaker solvent coupling with values as low as kc ∝ η-0.75 ± 0.02. Coupling of formation of short loops to solvent motions is very little affected by amino acid sequence, but solvent coupling of long-range loop formation is decreased by side chain sterics.

Exploring Intra- and Inter-Regional Interactions in the IDP α-Synuclein Using smFRET and MD Simulations

Theoretical concepts from polymer physics are often used to describe intrinsically disordered proteins (IDPs). However, amino acid interactions within and between regions of the protein can lead to deviations from typical polymer scaling behavior and even to short-lived secondary structures. To investigate the key interactions in the dynamic IDP α-synuclein (αS) at the amino acid level, we conducted single-molecule fluorescence resonance energy transfer (smFRET) experiments and coarse-grained molecular dynamics (CG-MD) simulations. We find excellent agreement between experiments and simulations. Our results show that a physiological salt solution is a good solvent for αS and that the protein is highly dynamic throughout its entire chain, with local intra- and inter-regional interactions leading to deviations from global scaling. Specifically, we observe expansion in the C-terminal region, compaction in the NAC region, and a slightly smaller distance between the C- and N-termini than expected. Our simulations indicate that the compaction in the NAC region results from hydrophobic aliphatic contacts, mostly between valine and alanine residues, and cation-π interactions between lysine and tyrosine. In addition, hydrogen bonds also seem to contribute to the compaction of the NAC region. The expansion of the C-terminal region is due to intraregional electrostatic repulsion and increased chain stiffness from several prolines. Overall, our study demonstrates the effectiveness of combining smFRET experiments with CG-MD simulations to investigate the key interactions in highly dynamic IDPs at the amino acid level.

Extreme dynamics in a biomolecular condensate

Proteins and nucleic acids can phase-separate in the cell to form concentrated biomolecular condensates1-4. The functions of condensates span many length scales: they modulate interactions and chemical reactions at the molecular scale5, organize biochemical processes at the mesoscale6 and compartmentalize cells4. Understanding the underlying mechanisms of these processes will require detailed knowledge of the rich dynamics across these scales7. The mesoscopic dynamics of biomolecular condensates have been extensively characterized8, but their behaviour at the molecular scale has remained more elusive. Here, as an example of biomolecular phase separation, we study complex coacervates of two highly and oppositely charged disordered human proteins9. Their dense phase is 1,000 times more concentrated than the dilute phase, and the resulting percolated interaction network10 leads to a bulk viscosity 300 times greater than that of water. However, single-molecule spectroscopy optimized for measurements within individual droplets reveals that at the molecular scale, the disordered proteins remain exceedingly dynamic, with their chain configurations interconverting on submicrosecond timescales. Massive all-atom molecular dynamics simulations reproduce the experimental observations and explain this apparent discrepancy: the underlying interactions between individual charged side chains are short-lived and exchange on a pico- to nanosecond timescale. Our results indicate that, despite the high macroscopic viscosity of phase-separated systems, local biomolecular rearrangements required for efficient reactions at the molecular scale can remain rapid.

Assessment of models for calculating the hydrodynamic radius of intrinsically disordered proteins

Diffusion measurements by pulsed-field gradient NMR and fluorescence correlation spectroscopy can be used to probe the hydrodynamic radius of proteins, which contains information about the overall dimension of a protein in solution. The comparison of this value with structural models of intrinsically disordered proteins is nonetheless impaired by the uncertainty of the accuracy of the methods for computing the hydrodynamic radius from atomic coordinates. To tackle this issue, we here build conformational ensembles of 11 intrinsically disordered proteins that we ensure are in agreement with measurements of compaction by small-angle x-ray scattering. We then use these ensembles to identify the forward model that more closely fits the radii derived from pulsed-field gradient NMR diffusion experiments. Of the models we examined, we find that the Kirkwood-Riseman equation provides the best description of the hydrodynamic radius probed by pulsed-field gradient NMR experiments. While some minor discrepancies remain, our results enable better use of measurements of the hydrodynamic radius in integrative modeling and for force field benchmarking and parameterization.

Effects of Sequence Composition, Patterning and Hydrodynamics on the Conformation and Dynamics of Intrinsically Disordered Proteins

Intrinsically disordered proteins (IDPs) and intrinsically disordered regions (IDRs) perform diverse functions in cellular organization, transport and signaling. Unlike the well-defined structures of the classical natively folded proteins, IDPs and IDRs dynamically span large conformational and structural ensembles. This dynamic disorder impedes the study of the relationship between the amino acid sequences of the IDPs and their spatial structures and dynamics, with different experimental techniques often offering seemingly contradictory results. Although experimental and theoretical evidence indicates that some IDP properties can be understood based on their average biophysical properties and amino acid composition, other aspects of IDP function are dictated by the specifics of the amino acid sequence. We investigate the effects of several key variables on the dimensions and the dynamics of IDPs using coarse-grained polymer models. We focus on the sequence "patchiness" informed by the sequence and biophysical properties of different classes of IDPs-and in particular FG nucleoporins of the nuclear pore complex (NPC). We show that the sequence composition and patterning are well reflected in the global conformational variables such as the radius of gyration and hydrodynamic radius, while the end-to-end distance and dynamics are highly sequence-specific. We find that in good solvent conditions highly heterogeneous sequences of IDPs can be well mapped onto averaged minimal polymer models for the purpose of prediction of the IDPs dimensions and dynamic relaxation times. The coarse-grained simulations are in a good agreement with the results of atomistic MD. We discuss the implications of these results for the interpretation of the recent experimental measurements, and for the further applications of mesoscopic models of FG nucleoporins and IDPs more broadly.

Illuminating Intrinsically Disordered Proteins with Integrative Structural Biology

Intense study of intrinsically disordered proteins (IDPs) did not begin in earnest until the late 1990s when a few groups, working independently, convinced the community that these 'weird' proteins could have important functions. Over the past two decades, it has become clear that IDPs play critical roles in a multitude of biological phenomena with prominent examples including coordination in signaling hubs, enabling gene regulation, and regulating ion channels, just to name a few. One contributing factor that delayed appreciation of IDP functional significance is the experimental difficulty in characterizing their dynamic conformations. The combined application of multiple methods, termed integrative structural biology, has emerged as an essential approach to understanding IDP phenomena. Here, we review some of the recent applications of the integrative structural biology philosophy to study IDPs.

Molecular Origin of Internal Friction in Intrinsically Disordered Proteins

Protein folding and dynamics are controlled by an interplay of thermal and viscosity effects. The effect of viscous drag through the solvent molecules is described by the classic Kramers theory in the high friction limit, which considers the dampening of the reactant molecules in the solution and quantifies the dependence of the reaction rate on the frictional drag. In addition to the external energy dissipation originating from the surrounding solvent molecules, there is an additional mode of internal energy dissipative force operative within the polypeptide chain reflecting the internal resistance of the chain to its conformational alterations. This dry, solvent-independent intrinsic frictional drag is termed internal friction. In the case of natively folded proteins, the physical origin of internal friction is primarily attributed to the intrachain interactions or other nonnative interactions in their compact states. However, the molecular origin of internal friction in intrinsically disordered proteins (IDPs) remains elusive.In this Account, we address this fundamental issue: what are the principal drivers of viscosity-independent (dry) friction in highly solvated, expanded, conformationally flexible, rapidly fluctuating IDPs that do not possess persistent intrachain interactions? IDPs exhibit diffusive conformational dynamics that is predominantly dominated by the segmental motion of the backbone arising due to the dihedral rotations in the Ramachandran Φ-Ψ space. The physical origin of friction in a complex biopolymeric system such as IDPs can be described by classic polymer models, namely, Rouse/Zimm models with internal friction. These one-dimensional models do not invoke torsional fluctuation components. Kuhn's classic description includes the connection between internal friction and microscopic dihedral hopping. Based on our time-resolved fluorescence anisotropy results, we describe that the sequence-dependent, collective, short-range backbone dihedral rotations govern localized internal friction in an archetypal IDP, namely, α-synuclein. The highly sensitive, residue-specific fluorescence depolarization kinetics offers a potent methodology to characterize and quantify the directional decorrelation engendered due to the short-range dihedral relaxation of the polypeptide backbone in the dihedral space. We utilized this characteristic relaxation time scale as our dynamic readout to quantify the site-specific frictional component. Our linear viscosity-dependent model of torsional relaxation time scale furnished a finite nonzero time constant at the zero solvent viscosity representing the solvent-independent internal friction. These results unveil the effect of the degree of dihedral restraining parameter on the internal friction component by showing that a restrained proline residue imparts higher torsional stiffness in the chain segments and, therefore, exhibits higher internal friction. This Account sheds light on the molecular underpinning of the sequence-specific internal friction in IDPs and will be of interest to unmask the role of internal friction in a diverse range of biomolecular processes involving binding-induced folding, allosteric interaction, protein misfolding and aggregation, and biomolecular condensation via phase separation.

The biophysics of disordered proteins from the point of view of single-molecule fluorescence spectroscopy

Intrinsically disordered proteins (IDPs) and regions (IDRs) have emerged as key players across many biological functions and diseases. Differently from structured proteins, disordered proteins lack stable structure and are particularly sensitive to changes in the surrounding environment. Investigation of disordered ensembles requires new approaches and concepts for quantifying conformations, dynamics, and interactions. Here, we provide a short description of the fundamental biophysical properties of disordered proteins as understood through the lens of single-molecule fluorescence observations. Single-molecule Förster resonance energy transfer (FRET) and fluorescence correlation spectroscopy (FCS) provides an extensive and versatile toolbox for quantifying the characteristics of conformational distributions and the dynamics of disordered proteins across many different solution conditions, both in vitro and in living cells.

Analytical Approaches for Deriving Friction Coefficients for Selected α-Helical Peptides Based Entirely on Molecular Dynamics Simulations

In this paper we derive analytically from molecular dynamics (MD) simulations the friction coefficients related to conformational transitions within several model peptides with α-helical structures. We study a series of alanine peptides with various length from ALA₅ to ALA₂₁ as well as their two derivatives, the (AAQAA)₃ peptide and a 13-residue KR1 peptide that is a derivative of the (AAQAA)₂ peptide with the formula GN(AAQAA)₂G. We use two kinds of approaches to derive their friction coefficients. In the local approach, friction associated with fluctuations of single hydrogen bonds are studied. In the second approach, friction coefficients associated with a folding transitions within the studied peptides are obtained. In both cases, the respective friction coefficients differentiated very well the subtle structural changes between studied peptides and compared favorably to experimentally available data.

Threading single proteins through pores to compare their energy landscapes

Protein function correlates with its structural dynamics. While theoretical approaches to studying protein energy landscapes are well developed, experimental methods that enable probing these landscapes of proteins remain challenging. We used solid-state nanopores to study the translocation behavior of three mutants of a helix bundle protein and quantified the number of energetically accessible conformational states for each mutant. We found that a slower-folding mutant with access to more conformational states translocates faster than a faster-folding mutant with a smaller number of accessible states, suggesting that ease of folding and ease of translocation are at odds in this case.

Translocation of proteins is correlated with structural fluctuations that access conformational states higher in free energy than the folded state. We use electric fields at the solid-state nanopore to control the relative free energy and occupancy of different protein conformational states at the single-molecule level. The change in occupancy of different protein conformations as a function of electric field gives rise to shifts in the measured distributions of ionic current blockades and residence times. We probe the statistics of the ionic current blockades and residence times for three mutants of the λ-repressor family in order to determine the number of accessible conformational states of each mutant and evaluate the ruggedness of their free energy landscapes. Translocation becomes faster at higher electric fields when additional flexible conformations are available for threading through the pore. At the same time, folding rates are not correlated with ease of translocation; a slow-folding mutant with a low-lying intermediate state translocates faster than a faster-folding two-state mutant. Such behavior allows us to distinguish among protein mutants by selecting for the degree of current blockade and residence time at the pore. Based on these findings, we present a simple free energy model that explains the complementary relationship between folding equilibrium constants and translocation rates.

Percolation transition prescribes protein size-specific barrier to passive transport through the nuclear pore complex

Nuclear pore complexes (NPCs) control biomolecular transport in and out of the nucleus. Disordered nucleoporins in the complex's pore form a permeation barrier, preventing unassisted transport of large biomolecules. Here, we combine coarse-grained simulations of experimentally derived NPC structures with a theoretical model to determine the microscopic mechanism of passive transport. Brute-force simulations of protein transport reveal telegraph-like behavior, where prolonged diffusion on one side of the NPC is interrupted by rapid crossings to the other. We rationalize this behavior using a theoretical model that reproduces the energetics and kinetics of permeation solely from statistics of transient voids within the disordered mesh. As the protein size increases, the mesh transforms from a soft to a hard barrier, enabling orders-of-magnitude reduction in permeation rate for proteins beyond the percolation size threshold. Our model enables exploration of alternative NPC architectures and sets the stage for uncovering molecular mechanisms of facilitated nuclear transport.

Toward protein NMR at physiological concentrations by hyperpolarized water-Finding and mapping uncharted conformational spaces

Nuclear magnetic resonance (NMR) spectroscopy is a key method for determining the structural dynamics of proteins in their native solution state. However, the low sensitivity of NMR typically necessitates nonphysiologically high sample concentrations, which often limit the relevance of the recorded data. We show how to use hyperpolarized water by dissolution dynamic nuclear polarization (DDNP) to acquire protein spectra at concentrations of 1 μM within seconds and with a high signal-to-noise ratio. The importance of approaching physiological concentrations is demonstrated for the vital MYC-associated factor X, which we show to switch conformations when diluted. While in vitro conditions lead to a population of the well-documented dimer, concentrations lowered by more than two orders of magnitude entail dimer dissociation and formation of a globularly folded monomer. We identified this structure by integrating DDNP with computational techniques to overcome the often-encountered constraint of DDNP of limited structural information provided by the typically detected one-dimensional spectra.

Variation of Structural and Dynamical Flexibility of Myelin Basic Protein in Response to Guanidinium Chloride

Myelin basic protein (MBP) is intrinsically disordered in solution and is considered as a conformationally flexible biomacromolecule. Here, we present a study on perturbation of MBP structure and dynamics by the denaturant guanidinium chloride (GndCl) using small-angle scattering and neutron spin-echo spectroscopy (NSE). A concentration of 0.2 M GndCl causes charge screening in MBP resulting in a compact, but still disordered protein conformation, while GndCl concentrations above 1 M lead to structural expansion and swelling of MBP. NSE data of MBP were analyzed using the Zimm model with internal friction (ZIF) and normal mode (NM) analysis. A significant contribution of internal friction was found in compact states of MBP that approaches a non-vanishing internal friction relaxation time of approximately 40 ns at high GndCl concentrations. NM analysis demonstrates that the relaxation rates of internal modes of MBP remain unaffected by GndCl, while structural expansion due to GndCl results in increased amplitudes of internal motions. Within the model of the Brownian oscillator our observations can be rationalized by a loss of friction within the protein due to structural expansion. Our study highlights the intimate coupling of structural and dynamical plasticity of MBP, and its fundamental difference to the behavior of ideal polymers in solution.

Measuring Photophysical Transition Rates with Fluorescence Correlation Spectroscopy and Antibunching

We present a new method that combines fluorescence correlation spectroscopy (FCS) on the microsecond time scale with fluorescence antibunching measurements on the nanosecond time scale for measuring photophysical rate constants of fluorescent molecules. The antibunching measurements allow us to quantify the average excitation rate of fluorescent molecules within the confocal detection volume of the FCS measurement setup. Knowledge of this value allows us then to quantify, in an absolute manner, the intersystem crossing rate and triplet state lifetime from the microsecond temporal decay of the FCS curves. We present a theoretical analysis of the method and estimate the maximum bias caused by the averaging of all quantities (excitation rate and photophysical rates) over the confocal detection volume, and we show that this bias is smaller than 5% in most cases. We apply the method for measuring the photophysical rate constants of the widely used dyes Rhodamine 110 and ATTO 655.

NMR Provides Unique Insight into the Functional Dynamics and Interactions of Intrinsically Disordered Proteins

Intrinsically disordered proteins are ubiquitous throughout all known proteomes, playing essential roles in all aspects of cellular and extracellular biochemistry. To understand their function, it is necessary to determine their structural and dynamic behavior and to describe the physical chemistry of their interaction trajectories. Nuclear magnetic resonance is perfectly adapted to this task, providing ensemble averaged structural and dynamic parameters that report on each assigned resonance in the molecule, unveiling otherwise inaccessible insight into the reaction kinetics and thermodynamics that are essential for function. In this review, we describe recent applications of NMR-based approaches to understanding the conformational energy landscape, the nature and time scales of local and long-range dynamics and how they depend on the environment, even in the cell. Finally, we illustrate the ability of NMR to uncover the mechanistic basis of functional disordered molecular assemblies that are important for human health.

Unraveling multi-state molecular dynamics in single-molecule FRET experiments. I. Theory of FRET-lines

Conformational dynamics of biomolecules are of fundamental importance for their function. Single-molecule studies of Förster Resonance Energy Transfer (smFRET) between a tethered donor and acceptor dye pair are a powerful tool to investigate the structure and dynamics of labeled molecules. However, capturing and quantifying conformational dynamics in intensity-based smFRET experiments remains challenging when the dynamics occur on the sub-millisecond timescale. The method of multiparameter fluorescence detection addresses this challenge by simultaneously registering fluorescence intensities and lifetimes of the donor and acceptor. Together, two FRET observables, the donor fluorescence lifetime τD and the intensity-based FRET efficiency E, inform on the width of the FRET efficiency distribution as a characteristic fingerprint for conformational dynamics. We present a general framework for analyzing dynamics that relates average fluorescence lifetimes and intensities in two-dimensional burst frequency histograms. We present parametric relations of these observables for interpreting the location of FRET populations in E-τD diagrams, called FRET-lines. To facilitate the analysis of complex exchange equilibria, FRET-lines serve as reference curves for a graphical interpretation of experimental data to (i) identify conformational states, (ii) resolve their dynamic connectivity, (iii) compare different kinetic models, and (iv) infer polymer properties of unfolded or intrinsically disordered proteins. For a simplified graphical analysis of complex kinetic networks, we derive a moment-based representation of the experimental data that decouples the motion of the fluorescence labels from the conformational dynamics of the biomolecule. Importantly, FRET-lines facilitate exploring complex dynamic models via easily computed experimental observables. We provide extensive computational tools to facilitate applying FRET-lines.

Tug-of-War between Internal and External Frictions and Viscosity Dependence of Rate in Biological Reactions

The role of water in biological processes is studied in three reactions, namely, the Fe-CO bond rupture in myoglobin, GB1 unfolding, and insulin dimer dissociation. We compute both internal and external components of friction on relevant reaction coordinates. In all of the three cases, the cross-correlation between forces from protein and water is found to be large and negative that serves to reduce the total friction significantly, increase the calculated reaction rate, and weaken solvent viscosity dependence. The computed force spectrum reveals bimodal 1/f noise, suggesting the use of a non-Markovian rate theory.

Release of linker histone from the nucleosome driven by polyelectrolyte competition with a disordered protein

Highly charged intrinsically disordered proteins are essential regulators of chromatin structure and transcriptional activity. Here we identify a surprising mechanism of molecular competition that relies on the pronounced dynamical disorder present in these polyelectrolytes and their complexes. The highly positively charged human linker histone H1.0 (H1) binds to nucleosomes with ultrahigh affinity, implying residence times incompatible with efficient biological regulation. However, we show that the disordered regions of H1 retain their large-amplitude dynamics when bound to the nucleosome, which enables the highly negatively charged and disordered histone chaperone prothymosin α to efficiently invade the H1-nucleosome complex and displace H1 via a competitive substitution mechanism, vastly accelerating H1 dissociation. By integrating experiments and simulations, we establish a molecular model that rationalizes the remarkable kinetics of this process structurally and dynamically. Given the abundance of polyelectrolyte sequences in the nuclear proteome, this mechanism is likely to be widespread in cellular regulation.

Spatial organization of chromosomes leads to heterogeneous chromatin motion and drives the liquid- or gel-like dynamical behavior of chromatin

Chromosome organization and dynamics are involved in regulating many fundamental processes such as gene transcription and DNA repair. Experiments unveiled that chromatin motion is highly heterogeneous inside cell nuclei, ranging from a liquid-like, mobile state to a gel-like, rigid regime. Using polymer modeling, we investigate how these different physical states and dynamical heterogeneities may emerge from the same structural mechanisms. We found that the formation of topologically associating domains (TADs) is a key driver of chromatin motion heterogeneity. In particular, we showed that the local degree of compaction of the TAD regulates the transition from a weakly compact, fluid state of chromatin to a more compact, gel state exhibiting anomalous diffusion and coherent motion. Our work provides a comprehensive study of chromosome dynamics and a unified view of chromatin motion enabling interpretation of the wide variety of dynamical behaviors observed experimentally across different biological conditions, suggesting that the "liquid" or "solid" state of chromatin are in fact two sides of the same coin.

Synergies of Single Molecule Fluorescence and NMR for the Study of Intrinsically Disordered Proteins

Single molecule fluorescence and nuclear magnetic resonance spectroscopy (NMR) are two very powerful techniques for the analysis of intrinsically disordered proteins (IDPs). Both techniques have individually made major contributions to deciphering the complex properties of IDPs and their interactions, and it has become evident that they can provide very complementary views on the distance-dynamics relationships of IDP systems. We now review the first approaches using both NMR and single molecule fluorescence to decipher the molecular properties of IDPs and their interactions. We shed light on how these two techniques were employed synergistically for multidomain proteins harboring intrinsically disordered linkers, for veritable IDPs, but also for liquid–liquid phase separated systems. Additionally, we provide insights into the first approaches to use single molecule Förster resonance energy transfer (FRET) and NMR for the description of multiconformational models of IDPs.

On distributions of barrier crossing times as observed in single-molecule studies of biomolecules

Single-molecule experiments that monitor time evolution of molecular observables in real time have expanded beyond measuring transition rates toward measuring distributions of times of various molecular events. Of particular interest is the first-passage time for making a transition from one molecular configuration ( a ) to another ( b ) and conditional first-passage times such as the transition path time, which is the first-passage time from a to b conditional upon not leaving the transition region intervening between a and b . Another experimentally accessible (but not yet studied experimentally) observable is the conditional exit time, i.e., the time to leave the transition region through a specified boundary. The distributions of such times contain a wealth of mechanistic information about the transitions in question. Here, we use the first and the second (and, if desired, higher) moments of these distributions to characterize their relative width for the model in which the experimental observable undergoes Brownian motion in a potential of mean force. We show that although the distributions of transition path times are always narrower than exponential (in that the ratio of the standard deviation to the distribution's mean is always less than 1), distributions of first-passage times and of conditional exit times can be either narrow or broad, in some cases displaying long power-law tails. The conditional exit time studied here provides a generalization of the transition path time that also allows one to characterize the temporal scales of failed barrier crossing attempts.

Quantitative Description of Intrinsically Disordered Proteins Using Single-Molecule FRET, NMR, and SAXS

Studying the conformational landscape of intrinsically disordered and partially folded proteins is challenging and only accessible to a few solution state techniques, such as nuclear magnetic resonance (NMR), small-angle scattering techniques, and single-molecule Förster resonance energy transfer (smFRET). While each of the techniques is sensitive to different properties of the disordered chain, such as local structural propensities, overall dimension, or intermediate- and long-range contacts, conformational ensembles describing intrinsically disordered proteins (IDPs) accurately should ideally respect all of these properties. Here we develop an integrated approach using a large set of FRET efficiencies and fluorescence lifetimes, NMR chemical shifts, and paramagnetic relaxation enhancements (PREs), as well as small-angle X-ray scattering (SAXS) to derive quantitative conformational ensembles in agreement with all parameters. Our approach is tested using simulated data (five sets of PREs and 15 FRET efficiencies) and validated experimentally on the example of the disordered domain of measles virus phosphoprotein, providing new insights into the conformational landscape of this viral protein that comprises transient structural elements and is more compact than an unfolded chain throughout its length. Rigorous cross-validation using FRET efficiencies, fluorescence lifetimes, and SAXS demonstrates the predictive nature of the calculated conformational ensembles and underlines the potential of this strategy in integrative dynamic structural biology.

Protein-Bath Coupling of an Internal Reaction Coordinate at Intermediate Time Scales

Thermally activated barrier-crossing processes are central to protein reaction kinetics. A determining factor for such kinetics is the extent to which the protein's motions are coupled to the surrounding bath. It is understood that slow large-scale conformational motions are strongly coupled to the environment, while fast librational motions are uncoupled. However, less is known about protein-bath coupling of reaction coordinates located on the interior of a protein and with dynamics on intermediate time scales. In this work, we use single molecule 2D fluorescence lifetime correlation spectroscopy to study the microsecond chemical reaction occurring in the chromophore pocket of eGFP. The equilibrium reaction involves a dihedral rotation of a glutamic acid residue and a rearrangement of the local hydrogen-bonding network surrounding the endogenous chromophore, with no accompanying large-scale conformational changes. We observe that the internal chemical reaction is coupled to the solvent viscosity, though the scaling deviates from Kramers' behavior. We attribute this deviation to the internal friction of the protein, which weakens the protein-solvent coupling at high viscosity and intermediate time scales.

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.

Probing the Internal Dynamics and Shape of Simple Peptides in Urea, Guanidinium Hydrochloride, and Proline Solutions with Time-Resolved Fluorescence Anisotropy and Atomistic Cosolvent Simulations

Picosecond time-resolved fluorescence anisotropy was used to measure the effect of denaturants and osmolytes on the reorientation dynamics of the simplest dipeptide. The solvent denaturants guanidinium hydrochloride (gdm), urea, and the osmolyte proline were used at several concentrations. Analysis of the concentration dependence of denaturants at a fixed temperature showed faster and slower reorientation time in two different denaturants at a nearly identical solvent viscosity (η). The reorientation time τ significantly deviates from Kramers' theory (τ ∝ η¹) in the high friction limit for guanidinium and urea with r ≈ 0.4 and r ≈ 0.6 at pH 7.2, respectively. In proline, τ is nearly proportional to η. Atomistic molecular dynamics simulations of the dipeptide in identical cosolvents showed excellent agreement with the measured rotational orientation time. The dipeptide dihedral (ϕ, ψ) isomerization times in water and 6 M urea are almost identical and significantly slower in guanidinium. If a faster and slower reorientation time can be associated with the compact and expanded shapes, the fractional viscosity dependence for guanidinium and urea may result from the fact that internal dynamics of peptides in these cosolvents involve higher and lower internal friction within the dynamic elements.

Diffusion of a disordered protein on its folded ligand

Flexibility in complexes between intrinsically disordered proteins and folded ligands is widespread in nature. However, timescales and spatial amplitudes of such dynamics remained unexplored for most systems. Our results show that the disordered cytoplasmic tail of the cell adhesion protein E-cadherin diffuses across the entire surface of its folded binding partner β-catenin at fast submillisecond timescales. The nanometer amplitude of these motions could allow kinases to access their recognition motifs without requiring a dissociation of the complex. We expect that the rugged energy landscape found in the E-cadherin/β-catenin complex is a defining feature of dynamic and partially disordered protein complexes.

Intrinsically disordered proteins often form dynamic complexes with their ligands. Yet, the speed and amplitude of these motions are hidden in classical binding kinetics. Here, we directly measure the dynamics in an exceptionally mobile, high-affinity complex. We show that the disordered tail of the cell adhesion protein E-cadherin dynamically samples a large surface area of the protooncogene β-catenin. Single-molecule experiments and molecular simulations resolve these motions with high resolution in space and time. Contacts break and form within hundreds of microseconds without a dissociation of the complex. The energy landscape of this complex is rugged with many small barriers (3 to 4 k_BT) and reconciles specificity, high affinity, and extreme disorder. A few persistent contacts provide specificity, whereas unspecific interactions boost affinity.

Integrating single-molecule spectroscopy and simulations for the study of intrinsically disordered proteins

Over the last two decades, intrinsically disordered proteins and protein regions (IDRs) have emerged from a niche corner of biophysics to be recognized as essential drivers of cellular function. Various techniques have provided fundamental insight into the function and dysfunction of IDRs. Among these techniques, single-molecule fluorescence spectroscopy and molecular simulations have played a major role in shaping our modern understanding of the sequence-encoded conformational behavior of disordered proteins. While both techniques are frequently used in isolation, when combined they offer synergistic and complementary information that can help uncover complex molecular details. Here we offer an overview of single-molecule fluorescence spectroscopy and molecular simulations in the context of studying disordered proteins. We discuss the various means in which simulations and single-molecule spectroscopy can be integrated, and consider a number of studies in which this integration has uncovered biological and biophysical mechanisms.

Slow Escape from a Helical Misfolded State of the Pore-Forming Toxin Cytolysin A

The pore-forming toxin cytolysin A (ClyA) is expressed as a large α-helical monomer that, upon interaction with membranes, undergoes a major conformational rearrangement into the protomer conformation, which then assembles into a cytolytic pore. Here, we investigate the folding kinetics of the ClyA monomer with single-molecule Förster resonance energy transfer spectroscopy in combination with microfluidic mixing, stopped-flow circular dichroism experiments, and molecular simulations. The complex folding process occurs over a broad range of time scales, from hundreds of nanoseconds to minutes. The very slow formation of the native state occurs from a rapidly formed and highly collapsed intermediate with large helical content and nonnative topology. Molecular dynamics simulations suggest pronounced non-native interactions as the origin of the slow escape from this deep trap in the free-energy surface, and a variational enhanced path-sampling approach enables a glimpse of the folding process that is supported by the experimental data.

How important are fluctuations in the treatment of internal friction in polymers?

The Rouse model with internal friction (RIF), a widely used theoretical framework to interpret the effects of internal friction on conformational transitions in biomolecules, is shown to be an approximate treatment that is based on preaveraging internal friction. By comparison with Brownian dynamics simulations of an exact coarse-grained model that incorporates fluctuations in internal friction, the accuracy of the preaveraged model predictions is examined both at and away from equilibrium. While the two models predict intrachain autocorrelations that approach each other for long enough chain segments, they differ in their predictions for shorter segments. Furthermore, the two models differ qualitatively in their predictions for the chain extension and viscosity in shear flow, which is taken to represent a prototypical out-of-equilibrium condition.

A multi-step nucleation process determines the kinetics of prion-like domain phase separation

Compartmentalization by liquid-liquid phase separation (LLPS) has emerged as a ubiquitous mechanism underlying the organization of biomolecules in space and time. Here, we combine rapid-mixing time-resolved small-angle X-ray scattering (SAXS) approaches to characterize the assembly kinetics of a prototypical prion-like domain with equilibrium techniques that characterize its phase boundaries and the size distribution of clusters prior to phase separation. We find two kinetic regimes on the micro- to millisecond timescale that are distinguished by the size distribution of clusters. At the nanoscale, small complexes are formed with low affinity. After initial unfavorable complex assembly, additional monomers are added with higher affinity. At the mesoscale, assembly resembles classical homogeneous nucleation. Careful multi-pronged characterization is required for the understanding of condensate assembly mechanisms and will promote understanding of how the kinetics of biological phase separation is encoded in biomolecules.

Does Electric Friction Matter in Living Cells?

The thermal motion of charged proteins causes randomly fluctuating electric fields inside cells. According to the fluctuation-dissipation theorem, there is an additional friction force associated with such fluctuations. However, the impact of these fluctuations on the diffusion and dynamics of proteins in the cytoplasm is unclear. Here, we provide an order-of-magnitude estimate of this effect by treating electric field fluctuations within a generalized Langevin equation model with a time-dependent friction memory kernel. We find that electric friction is generally negligible compared to solvent friction. However, a significant slowdown of protein diffusion and dynamics is expected for biomolecules with high net charges such as intrinsically disordered proteins and RNA. The results show that direct contacts between biomolecules in a cell are not necessarily required to alter their dynamics.