Modest Influence of FRET Chromophores on the Properties of Unfolded Proteins

Single-molecule Förster resonance energy transfer (FRET) experiments are often used to study the properties of unfolded and intrinsically disordered proteins. Because of their large extinction coefficients and quantum yields, synthetic heteroaromatic chromophores covalently linked to the protein are often used as donor and acceptor fluorophores. A key issue in the interpretation of such experiments is the extent to which the properties of the unfolded chain may be affected by the presence of these chromophores. In this article, we investigate this question using all-atom explicit solvent replica exchange molecular dynamics simulations of three different unfolded or intrinsically disordered proteins. We find that the secondary structure and long-range contacts are largely the same in the presence or absence of the fluorophores, and that the dimensions of the chain with and without chromophores are similar. This suggests that, at least in the cases studied, extrinsic fluorophores have little effect on the structural properties of unfolded or disordered proteins. We also find that the critical FRET orientational factor κ², has an average value and equilibrium distribution very close to that expected for isotropic orientations, which supports one of the assumptions frequently made when interpreting FRET efficiency in terms of distances.     

Using a FRET Library with Multiple Probe Pairs To Drive Monte Carlo Simulations of α-Synuclein

We describe a strategy for experimentally-constraining computational simulations of intrinsically disordered proteins (IDPs), using α-synuclein, an IDP with a central role in Parkinson’s disease pathology, as an example. Previously, data from single-molecule Förster Resonance Energy Transfer (FRET) experiments have been effectively utilized to generate experimentally constrained computational models of IDPs. However, the fluorophores required for single-molecule FRET experiments are not amenable to the study of short-range (<30 Å) interactions. Using ensemble FRET measurements allows one to acquire data from probes with multiple distance ranges, which can be used to constrain Monte Carlo simulations in PyRosetta. To appropriately employ ensemble FRET data as constraints, we optimized the shape and weight of constraining potentials to afford ensembles of structures that are consistent with experimental data. We also used this approach to examine the structure of α-synuclein in the presence of the compacting osmolyte trimethylamine-N-oxide. Despite significant compaction imparted by 2 M trimethylamine-N-oxide, the underlying ensemble of α-synuclein remains largely disordered and capable of aggregation, also in agreement with experimental data. These proof-of-concept experiments demonstrate that our modeling protocol enables one to efficiently generate experimentally constrained models of IDPs that incorporate atomic-scale detail, allowing one to study an IDP under a variety of conditions.     

Comprehensive determination of <Superscript>3</Superscript>J<Subscript>HNHα</Subscript> for unfolded proteins using <Superscript>13</Superscript>C′-resolved spin-echo difference spectroscopy

An experiment is presented to determine ³J_HNHα coupling constants, with significant advantages for applications to unfolded proteins. The determination of coupling constants for the peptide chain using 1D ¹H, or 2D and 3D ¹H-¹⁵N correlation spectroscopy is often hampered by extensive resonance overlap when dealing with flexible, disordered proteins. In the experiment detailed here, the overlap problem is largely circumvented by recording ¹H-¹³C′ correlation spectra, which demonstrate superior resolution for unfolded proteins. J-coupling constants are extracted from the peak intensities in a pair of 2D spin-echo difference experiments, affording rapid acquisition of the coupling data. In an application to the cytoplasmic domain of human neuroligin-3 (hNlg3cyt) data were obtained for 78 residues, compared to 54 coupling constants obtained from a 3D HNHA experiment. The coupling constants suggest that hNlg3cyt is intrinsically disordered, with little propensity for structure.     

Highly Disordered Amyloid-β Monomer Probed by Single-Molecule FRET and MD Simulation

Monomers of amyloid-β (Aβ) protein are known to be disordered, but there is considerable controversy over the existence of residual or transient conformations that can potentially promote oligomerization and fibril formation. We employed single-molecule Förster resonance energy transfer (FRET) spectroscopy with site-specific dye labeling using an unnatural amino acid and molecular dynamics simulations to investigate conformations and dynamics of Aβ isoforms with 40 (Aβ40) and 42 residues (Aβ42). The FRET efficiency distributions of both proteins measured in phosphate-buffered saline at room temperature show a single peak with very similar FRET efficiencies, indicating there is apparently only one state. 2D FRET efficiency-donor lifetime analysis reveals, however, that there is a broad distribution of rapidly interconverting conformations. Using nanosecond fluorescence correlation spectroscopy, we measured the timescale of the fluctuations between these conformations to be ～35 ns, similar to that of disordered proteins. These results suggest that both Aβ40 and Aβ42 populate an ensemble of rapidly reconfiguring unfolded states, with no long-lived conformational state distinguishable from that of the disordered ensemble. To gain molecular-level insights into these observations, we performed molecular dynamics simulations with a force field optimized to describe disordered proteins. We find, as in experiments, that both peptides populate configurations consistent with random polymer chains, with the vast majority of conformations lacking significant secondary structure, giving rise to very similar ensemble-averaged FRET efficiencies.     

Atomic-level characterization of disordered protein ensembles

The roles of unfolded states of proteins in normal folding and in diseases involving aggregation, as well as the prevalence and regulatory functions of intrinsically disordered proteins, have become increasingly recognized. The structural representation of these disordered states as ensembles of interconverting conformers can therefore provide critical insights. Experimental methods can be used to probe ensemble-averaged structural properties of disordered states and computational approaches generate representative ensembles of conformers using experimental restraints. In particular, NMR and small-angle X-ray scattering provide quantitative data that can readily be incorporated into calculations. These techniques have gleaned structural information about denatured, unfolded and intrinsically disordered proteins. The use of experimental data in different computational approaches, including ensemble molecular dynamics simulations and algorithms that assign populations to pregenerated conformers, has highlighted the presence of both local and long-range structure, and the occurrence of native-like and non-native interactions in unfolded and denatured states. Analysis of the resulting ensembles has suggested important implications of this fluctuating structure for folding, aggregation and binding.     

Minimal Effects of Macromolecular Crowding on an Intrinsically Disordered Protein: A Small-Angle Neutron Scattering Study

Small-angle neutron scattering was used to study the effects of macromolecular crowding by two globular proteins, i.e., bovine pancreatic trypsin inhibitor and equine metmyoglobin, on the conformational ensemble of an intrinsically disordered protein, the N protein of bacteriophage λ. The λ N protein was uniformly labeled with ²H, and the concentrations of D₂O in the samples were adjusted to match the neutron scattering contrast of the unlabeled crowding proteins, thereby masking their contribution to the scattering profiles. Scattering from the deuterated λ N was recorded for samples containing up to 0.12 g/mL bovine pancreatic trypsin inhibitor or 0.2 g/mL metmyoglobin. The radius of gyration of the uncrowded protein was estimated to be 30 Å and was found to be remarkably insensitive to the presence of crowders, varying by <2 Å for the highest crowder concentrations. The scattering profiles were also used to estimate the fractal dimension of λ N, which was found to be ∼1.8 in the absence or presence of crowders, indicative of a well-solvated and expanded random coil under all of the conditions examined. These results are contrary to the predictions of theoretical treatments and previous experimental studies demonstrating compaction of unfolded proteins by crowding with polymers such as dextran and Ficoll. A computational simulation suggests that some previous treatments may have overestimated the effective volumes of disordered proteins and the variation of these volumes within an ensemble. The apparent insensitivity of λ N to crowding may also be due in part to weak attractive interactions with the crowding proteins, which may compensate for the effects of steric exclusion.     

Structural Characterization of the Natively Unfolded N-Terminal Domain of Human c-Src Kinase: Insights into the Role of Phosphorylation of the Unique Domain

The N-terminal regions of the members of Src family of non-receptor protein tyrosine kinases are intrinsically unfolded and contain the maximum sequence divergence among them. In this study, we have addressed the structural characterization by nuclear magnetic resonance of this region of 84 residues that encompasses the SH4 and the unique domains (USrc) of the human c-Src. With this aim, the backbone assignment was performed using ¹³C-detected experiments that overcome the spectral resolution problems and the large number of prolines that are typical for intrinsically unfolded proteins. The analysis of the residual dipolar couplings measured for the USrc indicates the presence of a low populated helical structure in the 60-75 region. No long-range contacts between remote fragments of the chain were detected with paramagnetic relaxation enhancement experiments. The structural characterization was extended to two different phosphorylation states of USrc that encompassed three different phosphorylated sites, Ser17, Thr37, and Ser75. The structural and conformational changes upon phosphorylation were monitored through chemical shift perturbations and residual dipolar couplings, indicating that modifications occur at local level and no global rearrangements were apparent. These results suggest a scenario where phosphorylation induces a global electrostatic perturbation that could be involved in the membrane unbinding of c-Src and that could be related with the localization of the enzyme. These observations suggest the unique domain of Src kinases as a source of selectivity and reinforce the relevant role of intrinsically disordered proteins in biological processes.     

Targeting intrinsically disordered proteins in neurodegenerative and protein dysfunction diseases: another illustration of the D2 concept

Many biologically active proteins, which are usually called intrinsically disordered or natively unfolded proteins, lack stable tertiary and/or secondary structure under physiological conditions in vitro. Their functions complement the functional repertoire of ordered proteins, with intrinsically disordered proteins (IDPs) often being involved in regulation, signaling and control. Their amino acid sequences and compositions are very different from those of ordered proteins, making reliable identification of IDPs possible at the proteome level. IDPs are highly abundant in various human diseases, including neurodegeneration and other protein dysfunction maladies and, therefore, represent attractive novel drug targets. Some of the aspects of IDPs, as well as their roles in neurodegeneration and protein dysfunction diseases, are discussed in this article, together with the peculiarities of IDPs as potential drug targets.     

Backbone-driven collapse in unfolded protein chains

Collapse of unfolded protein chains is an early event in folding. It affects structural properties of intrinsically disordered proteins, which take a considerable fraction of the human proteome. Collapse is generally believed to be driven by hydrophobic forces imposed by the presence of nonpolar amino acid side chains. Contributions from backbone hydrogen bonds to protein folding and stability, however, are controversial. To date, the experimental dissection of side-chain and backbone contributions has not yet been achieved because both types of interactions are integral parts of protein structure. Here, we realized this goal by applying mutagenesis and chemical modification on a set of disordered peptides and proteins. We measured the protein dimensions and kinetics of intra-chain diffusion of modified polypeptides at the level of individual molecules using fluorescence correlation spectroscopy, thereby avoiding artifacts commonly caused by aggregation of unfolded protein material in bulk. We found no contributions from side chains to collapse but, instead, identified backbone interactions as a source sufficient to form globules of native-like dimensions. The presence of backbone hydrogen bonds decreased polypeptide water solubility dramatically and accelerated the nanosecond kinetics of loop closure, in agreement with recent predictions from computer simulation. The presence of side chains, instead, slowed loop closure and modulated the dimensions of intrinsically disordered domains. It appeared that the transient formation of backbone interactions facilitates the diffusive search for productive conformations at the early stage of folding and within intrinsically disordered proteins.     

Fluorescence resonance energy transfer-based intracellular assay for the conformation of hepatitis C virus drug target NS5A

Nonstructural protein 5A (NS5A) is essential for hepatitis C virus (HCV) replication and assembly and is a critical drug target. Biochemical data suggest large parts of NS5A are unfolded as an isolated protein, but little is known about its folded state in the cell. We used fluorescence resonance energy transfer (FRET) to probe whether or not different regions of NS5A are in close proximity within the cell. Twenty-three separate reporter constructs were created by inserting one or more fluorophores into different positions throughout the three domains of NS5A. FRET efficiency was maximal when donor and acceptor fluorophores were positioned next to each other but also could be observed when the two fluorophores flanked NS5A domain 1 or domain 3. Informatic and biochemical analysis suggests that large portions of the carboxy terminus of NS5A are in an unfolded and disordered state. Quercetin, a natural product known to disrupt NS5A function in cells, specifically disrupted a conformationally specific domain 3 FRET signal. Intermolecular FRET indicated that the NS5A amino termini, but not other regions, are in close proximity in multimeric complexes. Overall, this assay provides a new window on the intracellular conformation(s) of NS5A and how the conformation changes in response to cellular and viral components of the replication and assembly complex as well as antiviral drugs.     

Determination of tumor necrosis factor receptor-associated factor trimerization in living cells by CFP→YFP→mRFP FRET detected by flow cytometry

The availability of protein fluorophores with appropriate spectral properties has made it possible to employ fluorescence resonance energy transfer (FRET) to assess interactions between three proteins microscopically. Flow cytometry offers excellent sensitivity, effective signal separation and the capacity to assess a large number of events, and, therefore, should be an ideal means to explore protein interactions in living cells. Here, we report a flow-cytometric FRET technique that employed both direct energy transfer from CFP→YFP→mRFP and donor quenching to assess TRAF2 trimerization in living cells. Initially, a series of fusion proteins incorporating CFP, YFP and mRFP with spacers that did or did not permit FRET were employed to document the magnitude of CFP→YFP and YFP→mRFP FRET and to calculate the efficiency of CFP→YFP→mRFP two-step FRET. Based upon this, TRAF2 homotrimerization could be detected. This method should have great utility in studying the dynamics of interactions between three specific proteins in vivo.     

Internalization and Observation of Fluorescent Biomolecules in Living Microorganisms via Electroporation

The ability to study biomolecules in vivo is crucial for understanding their function in a biological context. One powerful approach involves fusing molecules of interest to fluorescent proteins such as GFP to study their expression, localization and function. However, GFP and its derivatives are significantly larger and less photostable than organic fluorophores generally used for in vitro experiments, and this can limit the scope of investigation. We recently introduced a straightforward, versatile and high-throughput method based on electroporation, allowing the internalization of biomolecules labeled with organic fluorophores into living microorganisms. Here we describe how to use electroporation to internalize labeled DNA fragments or proteins into Escherichia coli and Saccharomyces cerevisiæ, how to quantify the number of internalized molecules using fluorescence microscopy, and how to quantify the viability of electroporated cells. Data can be acquired at the single-cell or single-molecule level using fluorescence or FRET. The possibility of internalizing non-labeled molecules that trigger a physiological observable response in vivo is also presented. Finally, strategies of optimization of the protocol for specific biological systems are discussed.     

How multisite phosphorylation impacts the conformations of intrinsically disordered proteins

Phosphorylation of intrinsically disordered proteins (IDPs) can produce changes in structural and dynamical properties and thereby mediate critical biological functions. How phosphorylation effects intrinsically disordered proteins has been studied for an increasing number of IDPs, but a systematic understanding is still lacking. Here, we compare the collapse propensity of four disordered proteins, Ash1, the C-terminal domain of RNA polymerase (CTD2’), the cytosolic domain of E-Cadherin, and a fragment of the p130Cas, in unphosphorylated and phosphorylated forms using extensive all-atom molecular dynamics (MD) simulations. We find all proteins to show V-shape changes in their collapse propensity upon multi-site phosphorylation according to their initial net charge: phosphorylation expands neutral or overall negatively charged IDPs and shrinks positively charged IDPs. However, force fields including those tailored towards and commonly used for IDPs overestimate these changes. We find quantitative agreement of MD results with SAXS and NMR data for Ash1 and CTD2’ only when attenuating protein electrostatic interactions by using a higher salt concentration (e.g. 350 mM), highlighting the overstabilization of salt bridges in current force fields. We show that phosphorylation of IDPs also has a strong impact on the solvation of the protein, a factor that in addition to the actual collapse or expansion of the IDP should be considered when analyzing SAXS data. Compared to the overall mild change in global IDP dimension, the exposure of active sites can change significantly upon phosphorylation, underlining the large susceptibility of IDP ensembles to regulation through post-translational modifications.     

Contribution of Fluorophore Dynamics and Solvation to Resonant Energy Transfer in Protein-DNA Complexes: A Molecular-Dynamics Study

In Förster resonance energy transfer (FRET) experiments, extracting accurate structural information about macromolecules depends on knowing the positions and orientations of donor and acceptor fluorophores. Several approaches have been employed to reduce uncertainties in quantitative FRET distance measurements. Fluorophore-position distributions can be estimated by surface accessibility (SA) calculations, which compute the region of space explored by the fluorophore within a static macromolecular structure. However, SA models generally do not take fluorophore shape, dye transition-moment orientation, or dye-specific chemical interactions into account. We present a detailed molecular-dynamics (MD) treatment of fluorophore dynamics for an ATTO donor/acceptor dye pair and specifically consider as case studies dye-labeled protein-DNA intermediates in Cre site-specific recombination. We carried out MD simulations in both an aqueous solution and glycerol/water mixtures to assess the effects of experimental solvent systems on dye dynamics. Our results unequivocally show that MD simulations capture solvent effects and dye-dye interactions that can dramatically affect energy transfer efficiency. We also show that results from SA models and MD simulations strongly diverge in cases where donor and acceptor fluorophores are in close proximity. Although atomistic simulations are computationally more expensive than SA models, explicit MD studies are likely to give more realistic results in both homogeneous and mixed solvents. Our study underscores the model-dependent nature of FRET analyses, but also provides a starting point to develop more realistic in silico approaches for obtaining experimental ensemble and single-molecule FRET data.     

Protein folding at single-molecule resolution

The protein folding reaction carries great significance for cellular function and hence continues to be the research focus of a large interdisciplinary protein science community. Single-molecule methods are providing new and powerful tools for dissecting the mechanisms of this complex process by virtue of their ability to provide views of protein structure and dynamics without associated ensemble averaging. This review briefly introduces common FRET and force methods, and then explores several areas of protein folding where single-molecule experiments have yielded insights. These include exciting new information about folding landscapes, dynamics, intermediates, unfolded ensembles, intrinsically disordered proteins, assisted folding and biomechanical unfolding. Emerging and future work is expected to include advances in single-molecule techniques aimed at such investigations, and increasing work on more complex systems from both the physics and biology standpoints, including folding and dynamics of systems of interacting proteins and of proteins in cells and organisms. This article is part of a Special Issue entitled: Protein Dynamics: Experimental and Computational Approaches.