Resonance Energy Transfer Data Research Articles

Molecular Dynamics-Assisted Interpretation of Experimentally Determined Intrinsically Disordered Protein Conformational Components: The Case of Human α-Synuclein.

Mass spectrometry and single molecule force microscopy are two experimental approaches able to provide structural information on intrinsically disordered proteins (IDPs). These techniques allow the dissection of conformational ensembles in their main components, although at a low-resolution level. In this work, we interpret the results emerging from these experimental approaches on human alpha synuclein (AS) by analyzing a previously published 73 μs-long molecular-dynamics (MD) simulation of the protein in explicit solvent. We further compare MD-based predictions of single molecule Förster resonance energy transfer (smFRET) data of AS in solution with experimental data. The combined theoretical and experimental data provide a description of AS main conformational ensemble, shedding light into its intramolecular interactions and overall structural compactness. This analysis could be easily transferred to other IDPs.

Millisecond molecular dynamics simulations of KRas-dimer formation and interfaces

Ras dimers have been proposed as building blocks for initiating the extracellular signal-regulated kinase (ERK)/mitogen-activated protein kinase (MAPK) cellular signaling pathway. To better examine the structure of possible dimer interfaces, the dynamics of Ras dimerization, and its potential signaling consequences, we performed molecular dynamics simulations totaling 1 ms of sampling, using an all-atom model of two full-length, farnesylated, guanosine triphosphate (GTP)-bound, wild-type KRas4b proteins diffusing on 29%POPS (1-palmitoyl-2-oleoyl-sn-glycero-3-phospho-L-serine)-mixed POPC (1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine) membranes. Our simulations unveil an ensemble of thermodynamically weak KRas dimers spanning multiple conformations. The most stable conformations, having the largest interface areas, involve helix α2 and a hypervariable region (HVR). Among the dimer conformations, we found that the HVR of each KRas has frequent interactions with various parts of the dimer, thus potentially mediating the dimerization. Some dimer configurations have one KRas G-domain elevated above the lipid bilayer surface by residing on top of the other G-domain, thus likely contributing to the recruitment of cytosolic Raf kinases in the context of a stably formed multi-protein complex. We identified a variant of the α4-α5 KRas-dimer interface that is similar to the interfaces obtained with fluorescence resonance energy transfer (FRET) data of HRas on lipid bilayers. Interestingly, we found two arginine fingers, R68 and R149, that directly interact with the beta-phosphate of the GTP bound in KRas, in a manner similar to what is observed in a crystal structure of GAP-HRas complex, which can facilitate the GTP hydrolysis via the arginine finger of GTPase-activating protein (GAP).

Salt-dependent conformational changes of intrinsically disordered proteins

It is well known that the structure of an intrinsically disordered protein (IDP) can be affected by changes in salt concentration as salt ions screen the electrostatic interactions between charged amino acids. Many IPDs, however, contain more uncharged amino acids than charged ones. These acids are not subjected to electrostatic screening, but are influenced by the salting-out effect. Using recent Förster resonance energy transfer data, we have parameterized the salting-out effect into a previously established coarse-grained model.

Folding and misfolding of potassium channel monomers during assembly and tetramerization

The dynamics and folding of potassium channel pore domain monomers are connected to the kinetics of tetramer assembly. In all-atom molecular dynamics simulations of Kv1.2 and KcsA channels, monomers adopt multiple nonnative conformations while the three helices remain folded. Consistent with this picture, NMR studies also find the monomers to be dynamic and structurally heterogeneous. However, a KcsA construct with a disulfide bridge engineered between the two transmembrane helices has an NMR spectrum with well-dispersed peaks, suggesting that the monomer can be locked into a native-like conformation that is similar to that observed in the folded tetramer. During tetramerization, fluoresence resonance energy transfer (FRET) data indicate that monomers rapidly oligomerize upon insertion into liposomes, likely forming a protein-dense region. Folding within this region occurs along separate fast and slow routes, with τfold ∼40 and 1,500 s, respectively. In contrast, constructs bearing the disulfide bond mainly fold via the faster pathway, suggesting that maintaining the transmembrane helices in their native orientation reduces misfolding. Interestingly, folding is concentration independent despite the tetrameric nature of the channel, indicating that the rate-limiting step is unimolecular and occurs after monomer association in the protein-dense region. We propose that the rapid formation of protein-dense regions may help with the assembly of multimeric membrane proteins by bringing together the nascent components prior to assembly. Finally, despite its name, the addition of KcsA's C-terminal "tetramerization" domain does not hasten the kinetics of tetramerization.

Salt-Dependent Conformational Changes of Intrinsically Disordered Proteins.

The flexible structure of an intrinsically disordered protein (IDP) is known to be perturbed by salt concentrations, which can be understood by electrostatic screening on charged amino acids. However, an IDP usually contains more uncharged residues that are influenced by the salting-out effect. Here we have parametrized the salting-out effect into a coarse-grained model using a set of Förster resonance energy transfer data and verified with experimental salt-dependent liquid-liquid phase separation (LLPS) of 17 proteins. The new model can correctly capture the behavior of 6 more sequences, resulting in a total of 13 when varying salt concentrations. Together with a survey of more than 500 IDP sequences, we conclude that the salting-out effect, which was considered to be secondary to electrostatic screening, is important for IDP sequences with moderately charged residues at physiological salt concentrations. The presented scheme is generally applicable to other computational models for capturing salt-dependent IDP conformations.

High-resolution crystal structure of the Borreliella burgdorferi PlzA protein in complex with c-di-GMP: new insights into the interaction of c-di-GMP with the novel xPilZ domain.

ABSTRACTIn the tick-borne pathogens, Borreliella burgdorferi and Borrelia hermsii, c-di-GMP is produced by a single diguanylate cyclase (Rrp1). In these pathogens, the Plz proteins (PlzA, B and C) are the only c-di-GMP receptors identified to date and PlzA is the sole c-di-GMP receptor found in all Borreliella isolates. Bioinformatic analyses suggest that PlzA has a unique PilZN3-PilZ architecture with the relatively uncommon xPilZ domain. Here, we present the crystal structure of PlzA in complex with c-di-GMP (1.6 Å resolution). This is the first structure of a xPilz domain in complex with c-di-GMP to be determined. PlzA has a two-domain structure, where each domain comprises topologically equivalent PilZ domains with minimal sequence identity but remarkable structural similarity. The c-di-GMP binding site is formed by the linker connecting the two domains. While the structure of apo PlzA could not be determined, previous fluorescence resonance energy transfer data suggest that apo and holo forms of the protein are structurally distinct. The information obtained from this study will facilitate ongoing efforts to identify the molecular mechanisms of PlzA-mediated regulation in ticks and mammals.

Extraction of rapid kinetics from smFRET measurements using integrative detectors

SUMMARYHidden Markov models (HMMs) are used to learn single-molecule kinetics across a range of experimental techniques. By their construction, HMMs assume that single-molecule events occur on slower timescales than those of data acquisition. To move beyond that HMM limitation and allow for single-molecule events to occur on any timescale, we must treat single-molecule events in continuous time as they occur in nature. We propose a method to learn kinetic rates from single-molecule Förster resonance energy transfer (smFRET) data collected by integrative detectors, even if those rates exceed data acquisition rates. To achieve that, we exploit our recently proposed “hidden Markov jump process” (HMJP), with which we learn transition kinetics from parallel measurements in donor and acceptor channels. HMJPs generalize the HMM paradigm in two critical ways: (1) they deal with physical smFRET systems as they switch between conformational states in continuous time, and (2) they estimate transition rates between conformational states directly without having recourse to transition probabilities or assuming slow dynamics. Our continuous-time treatment learns the transition kinetics and photon emission rates for dynamic regimes that are inaccessible to HMMs, which treat system kinetics in discrete time. We validate our framework’s robustness on simulated data and demonstrate its performance on experimental data from FRET-labeled Holliday junctions.

Automatic classification and segmentation of single-molecule fluorescence time traces with deep learning

Traces from single-molecule fluorescence microscopy (SMFM) experiments exhibit photophysical artifacts that typically necessitate human expert screening, which is time-consuming and introduces potential for user-dependent expectation bias. Here, we use deep learning to develop a rapid, automatic SMFM trace selector, termed AutoSiM, that improves the sensitivity and specificity of an assay for a DNA point mutation based on single-molecule recognition through equilibrium Poisson sampling (SiMREPS). The improved performance of AutoSiM is based on accepting both more true positives and fewer false positives than the conventional approach of hidden Markov modeling (HMM) followed by hard thresholding. As a second application, the selector is used for automated screening of single-molecule Förster resonance energy transfer (smFRET) data to identify high-quality traces for further analysis, and achieves ~90% concordance with manual selection while requiring less processing time. Finally, we show that AutoSiM can be adapted readily to novel datasets, requiring only modest Transfer Learning.

Workflow Description to Dynamically Model β-Arrestin Signaling Networks.

Dynamic models of signaling networks allow the formulation of hypotheses on the topology and kinetic rate laws characterizing a given molecular network, in-depth exploration, and confrontation with kinetic biological data. Despite its standardization, dynamic modeling of signaling networks still requires successive technical steps that need to be carefully performed. Here, we detail these steps by going through the mathematical and statistical framework. We explain how it can be applied to the understanding of β-arrestin-dependent signaling networks. We illustrate our methodology through the modeling of β-arrestin recruitment kinetics at the follicle-stimulating hormone (FSH) receptor supported by in-house bioluminescence resonance energy transfer (BRET) data.

Sequence-Dependent Effects of Monovalent Cations on the Structural Dynamics of Trinucleotide-Repeat DNA Hairpins.

Repetitive trinucleotide DNA sequences at specific genetic loci are associated with numerous hereditary, neurodegenerative diseases. The propensity of single-stranded domains containing these sequences to form secondary structure via extensive self-complementarity disrupts normal DNA processing to create genetic instabilities. To investigate these intrastrand structural dynamics, a DNA hairpin system was devised for single-molecule fluorescence study of the folding kinetics and energetics for secondary structure formation between two interacting, repetitive domains with specific numbers of the same trinucleotide motif (CXG), where X = T or A. Single-molecule fluorescence resonance energy transfer (smFRET) data show discrete conformational transitions between unstructured and closed hairpin states. The lifetimes of the closed hairpin states correlate with the number of repeats, with (CTG) N/(CTG) N domains maintaining longer-lived, closed states than equivalent-sized (CAG) N/(CAG) N domains. NaCl promotes similar degree of stabilization for the closed hairpin states of both repeat sequences. Temperature-based, smFRET experiments reveal that NaCl favors hairpin closing for (CAG) N/(CAG) N by preordering single-stranded repeat domains to accelerate the closing transition. In contrast, NaCl slows the opening transition of CTG hairpins; however, it promotes misfolded conformations that require unfolding. Energy diagrams illustrate the distinct folding pathways of (CTG) N and (CAG) N repeat domains and identify features that may contribute to their gene-destabilizing effects.

Imaging Protein-Protein Interactions by Förster Resonance Energy Transfer (FRET) Microscopy in Live Cells.

This updated unit compares three methods to acquire Förster Resonance Energy Transfer (FRET) data in living cells using a confocal microscope: Acceptor photobleaching, Acceptor-sensitized emission FRET, and Donor fluorescence lifetime imaging. Detailed protocols for live cell husbandry, image acquisition, and data analysis are provided. In addition to providing instructions for manufacturer's analysis tool sets, we provide an easy-to-use, MATLAB-based code to calculate FRET efficiency from data obtained using the Acceptor photobleaching or Acceptor-sensitized emission method, which can be freely downloaded. © 2018 by John Wiley & Sons, Inc.

Linking time-series of single-molecule experiments with molecular dynamics simulations by machine learning.

Single-molecule experiments and molecular dynamics (MD) simulations are indispensable tools for investigating protein conformational dynamics. The former provide time-series data, such as donor-acceptor distances, whereas the latter give atomistic information, although this information is often biased by model parameters. Here, we devise a machine-learning method to combine the complementary information from the two approaches and construct a consistent model of conformational dynamics. It is applied to the folding dynamics of the formin-binding protein WW domain. MD simulations over 400 μs led to an initial Markov state model (MSM), which was then "refined" using single-molecule Förster resonance energy transfer (FRET) data through hidden Markov modeling. The refined or data-assimilated MSM reproduces the FRET data and features hairpin one in the transition-state ensemble, consistent with mutation experiments. The folding pathway in the data-assimilated MSM suggests interplay between hydrophobic contacts and turn formation. Our method provides a general framework for investigating conformational transitions in other proteins.

Nuclear Transport and Accumulation of Smad Proteins Studied by Single-Molecule Microscopy

Nuclear translocation of stimulated Smad heterocomplexes is a critical step in the signal transduction of transforming growth factor β (TGF-β) from transmembrane receptors into the nucleus. Specifically, normal nuclear accumulation of Smad2/Smad4 heterocomplexes induced by TGF-β1 is involved in carcinogenesis. However, the relationship between nuclear accumulation and the nucleocytoplasmic transport kinetics of Smad proteins in the presence of TGF-β1 remains obscure. By combining a high-speed single-molecule tracking microscopy and Förster resonance energy transfer technique, we tracked the entire TGF-β1-induced process of Smad2/Smad4 heterocomplex formation, as well as their transport through nuclear pore complexes in live cells, with a high single-molecule localization precision of 2 ms and <20 nm. Our single-molecule Förster resonance energy transfer data have revealed that in TGF-β1-treated cells, Smad2/Smad4 heterocomplexes formed in the cytoplasm, imported through the nuclear pore complexes as entireties, and finally dissociated in the nucleus. Moreover, we found that basal-state Smad2 or Smad4 cannot accumulate in the nucleus without the presence of TGF-β1, mainly because both of them have an approximately twofold higher nuclear export efficiency compared to their nuclear import. Remarkably and reversely, heterocomplexes of Smad2/Smad4 induced by TGF-β1 can rapidly concentrate in the nucleus because of their almost fourfold higher nuclear import rate in comparison with their nuclear export rate. Thus, we believe that the determined TGF-β1-dependent transport configurations and efficiencies for the basal-state Smad or stimulated Smad heterocomplexes elucidate the basic molecular mechanism to understand their nuclear transport and accumulation.

FRET Detects the Size of Nanodomains for Coexisting Liquid-Disordered and Liquid-Ordered Phases

Biomembranes with as few as three lipid components can form coexisting liquid-disordered (Ld) and liquid-ordered (Lo) phases. In the coexistence region of Ld and Lo phases, the lipid mixtures 1,2-distearoyl-sn-glycero-3-phosphocholine (DSPC)/1,2-dioleoyl-sn-glycero-3-phosphocholine (DOPC)/chol or brain sphingomyelin (bSM)/DOPC/chol form micron-scale domains that are easily visualized with light microscopy. Although large domains are not observed in the mixtures DSPC/1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC)/chol and bSM/POPC/chol, lateral heterogeneity is nevertheless detected using techniques with nanometer-scale spatial resolution. We propose a simple and accessible method to measure domain sizes below optical resolution (∼200 nm). We measured nanodomain size for the latter two mixtures by combining experimental Förster resonance energy transfer data with a Monte-Carlo-based analysis. We found a domain radius of 7.5−10 nm for DSPC/POPC/chol, similar to values obtained previously by neutron scattering, and ∼5 nm for bSM/POPC/chol, slightly smaller than measurable by neutron scattering. These analyses also detect the domain-size transition that is observed by fluorescence microscopy in the four-component lipid mixture bSM/DOPC/POPC/chol. Accurate measurements of fluorescent-probe partition coefficients are especially important for the analysis; therefore, we exploit three different methods to measure the partition coefficient of fluorescent molecules between Ld and Lo phases.

Determination of the fluorescence labels location in lipid bilayer based on fluorescence quenching

A comparison of Förster Resonance Energy Transfer (FRET) data of vesicles labeled with fluorescence labels on both inner and outer layers of the bilayer, with ones labeled only on the inner layer of bilayer to permit a determination of the location of fluorescence labels in relation to the bilayer center. The theoretical description of the FRET effect is based on the FRET model generalized and extended by Wolber and Hudson. A useful version of the FRET model implies location of acceptors with respect to a trigonal or tetragonal lattice of acceptors centered on donors located deeper than the acceptor layer. Several fluorescence labels were tested in order to find FRET pairs of labels for which the same quantity of quencher reduced label fluorescence to one half of the initial value before quencher treatment. 2,2,6,6-tetramethylpiperidine-1-oxyl (TEMPO) was used as a fluorescence quencher of both labels located in outer layer of bilayer.

Crystal Structures of Cyanine Fluorophores Stacked onto the End of Double-Stranded RNA

The indodicarbocyanine fluorophores Cy3 and Cy5 are extensively used as donor-acceptor pairs in fluorescence resonance energy transfer experiments, especially those involving single molecules. When terminally attached to double-stranded nucleic acids via the 5′ phosphate group these fluorophores stack onto the ends of the molecule. Knowledge of the positions of the fluorophores is critical to the interpretation of fluorescence resonance energy transfer data. The positions have been demonstrated for double-stranded (ds) DNA using NMR spectroscopy. Here, we have used x-ray crystallography to analyze the location of Cy3 and Cy5 on dsRNA, using complexes of an RNA stem-loop bound to L5 protein determined at 2.4 Å resolution. This confirms the tendency of both fluorophores to stack on the free end of RNA, with the long axis of the fluorophores approximately parallel to that of the terminal basepair. However, the manner of interaction of both Cy3 and Cy5 with the terminus of the dsRNA is significantly different from that deduced for dsDNA using NMR. The fluorophores are stacked on the terminal basepair such that their indole nitrogen atoms lie on the major groove side, and thus their pendant methyl groups are on the minor groove side.

Characterizing Protein Dynamics with Integrative Use of Bulk and Single-Molecule Techniques.

A protein dynamically samples multiple conformations, and the conformational dynamics enables protein function. Most biophysical measurements are ensemble-based, with the observables averaged over all members of the ensemble. Though attainable, the decomposition of the observables to the constituent conformational states can be computationally expensive and ambiguous. Here we show that the incorporation of single-molecule fluorescence resonance energy transfer (smFRET) data resolves the ambiguity and affords protein ensemble structures that are more precise and accurate. Using K63-linked diubiquitin, we characterize the dynamic domain arrangements of the model system, with the use of chemical cross-linking coupled with mass spectrometry (CXMS), small-angle X-ray scattering (SAXS), and smFRET techniques. CXMS allows the modeling of protein conformational states that are alternatives to the crystal structure. SAXS provides ensemble-averaged low-resolution shape information. Importantly, smFRET affords state-specific populations, and the FRET distances validate the ensemble structures obtained by refining against CXMS and SAXS restraints. Together, the integrative use of bulk and single-molecule techniques affords better insight into protein dynamics and shall be widely implemented in structural biology.

TOWARD BAYESIAN INFERENCE OF THE SPATIAL DISTRIBUTION OF PROTEINS FROM THREE-CUBE FÖRSTER RESONANCE ENERGY TRANSFER DATA.

Förster resonance energy transfer (FRET) is a quantum-physical phenomenon where energy may be transferred from one molecule to a neighbor molecule if the molecules are close enough. Using fluorophore molecule marking of proteins in a cell, it is possible to measure in microscopic images to what extent FRET takes place between the fluorophores. This provides indirect information of the spatial distribution of the proteins. Questions of particular interest are whether (and if so to which extent) proteins of possibly different types interact or whether they appear independently of each other. In this paper we propose a new likelihood-based approach to statistical inference for FRET microscopic data. The likelihood function is obtained from a detailed modeling of the FRET data-generating mechanism conditional on a protein configuration. We next follow a Bayesian approach and introduce a spatial point process prior model for the protein configurations depending on hyperparameters quantifying the intensity of the point process. Posterior distributions are evaluated using Markov chain Monte Carlo. We propose to infer microscope-related parameters in an initial step from reference data without interaction between the proteins. The new methodology is applied to simulated and real datasets.

The Arginine Pairs and C-Termini of the Sso7c4 from Sulfolobus solfataricus Participate in Binding and Bending DNA.

The Sso7c4 from Sulfolobus solfataricus forms a dimer, which is believed to function as a chromosomal protein involved in genomic DNA compaction and gene regulation. Here, we present the crystal structure of wild-type Sso7c4 at a high resolution of 1.63 Å, showing that the two basic C-termini are disordered. Based on the fluorescence polarization (FP) binding assay, two arginine pairs, R11/R22′ and R11′/R22, on the top surface participate in binding DNA. As shown in electron microscopy (EM) images, wild-type Sso7c4 compacts DNA through bridging and bending interactions, whereas the binding of C-terminally truncated proteins rigidifies and opens DNA molecules, and no compaction of the DNA occurs. Moreover, the FP, EM and fluorescence resonance energy transfer (FRET) data indicated that the two basic and flexible C-terminal arms of the Sso7c4 dimer play a crucial role in binding and bending DNA. Sso7c4 has been classified as a repressor-like protein because of its similarity to Escherichia coli Ecrep 6.8 and Ecrep 7.3 as well as Agrobacterium tumefaciens ACCR in amino acid sequence. Based on these data, we proposed a model of the Sso7c4-DNA complex using a curved DNA molecule in the catabolite activator protein-DNA complex. The DNA end-to-end distance measured with FRET upon wild-type Sso7c4 binding is almost equal to the distance measured in the model, which supports the fidelity of the proposed model. The FRET data also confirm the EM observation showing that the binding of wild-type Sso7c4 reduces the DNA length while the C-terminal truncation does not. A functional role for Sso7c4 in the organization of chromosomal DNA and/or the regulation of gene expression through bridging and bending interactions is suggested.

Significant Heterogeneity and Slow Dynamics of the Unfolded Ubiquitin Detected by the Line Confocal Method of Single-Molecule Fluorescence Spectroscopy.

The conformation and dynamics of the unfolded state of ubiquitin doubly labeled regiospecifically with Alexa488 and Alexa647 were investigated using single-molecule fluorescence spectroscopy. The line confocal fluorescence detection system combined with the rapid sample flow enabled the characterization of unfolded proteins at the improved structural and temporal resolutions compared to the conventional single-molecule methods. In the initial stage of the current investigation, however, the single-molecule Förster resonance energy transfer (sm-FRET) data of the labeled ubiquitin were flawed by artifacts caused by the adsorption of samples to the surfaces of the fused-silica flow chip and the sample delivery system. The covalent coating of 2-methacryloyloxyethyl phosphorylcholine polymer to the flow chip surface was found to suppress the artifacts. The sm-FRET measurements based on the coated flow chip demonstrated that the histogram of the sm-FRET efficiencies of ubiquitin at the native condition were narrowly distributed, which is comparable to the probability density function (PDF) expected from the shot noise, demonstrating the structural homogeneity of the native state. In contrast, the histogram of the sm-FRET efficiencies of the unfolded ubiquitin obtained at a time resolution of 100 μs was distributed significantly more broadly than the PDF expected from the shot noise, demonstrating the heterogeneity of the unfolded state conformation. The variety of the sm-FRET efficiencies of the unfolded state remained even after evaluating the moving average of traces with a window size of 1 ms, suggesting that conformational averaging of the heterogeneous conformations mostly occurs in the time domain slower than 1 ms. Local structural heterogeneity around the labeled fluorophores was inferred as the cause of the structural heterogeneity. The heterogeneity and slow dynamics revealed by the line confocal tracking of sm-FRET might be common properties of the unfolded proteins.