Qui, presentiamo un protocollo per ottenere dati di tre colori smFRET e la sua analisi con un ensemble 3D modello di Markov nascosto. Con questo approccio, gli scienziati possono estrarre informazioni cinetiche da sistemi proteici complessi, tra cui il cooperativity o interazioni correlate.
Trasferimento di energia di risonanza di singola molecola Förster (smFRET) è diventata una tecnica ampiamente utilizzata biofisica per studiare la dinamica delle biomolecole. Per molte macchine molecolari in una proteine delle cellule devono agire insieme ai partner di interazione in un ciclo funzionale per adempiere il loro compito. L’estensione di due colori di multi-colore smFRET permette di sondare contemporaneamente più di una interazione o cambiamento conformazionale. Questo non solo aggiunge una nuova dimensione agli esperimenti di smFRET, ma offre anche la possibilità di studiare direttamente la sequenza di eventi e di rilevare interazioni correlate quando si utilizza un campione immobilizzato ed una fluorescenza di riflessione interna totale microscopio (TIRFM). Di conseguenza, multi-colore smFRET è un versatile strumento per lo studio biomolecolare complessi in maniera quantitativa e dettagliatamente in precedenza irraggiungibili.
Qui, dimostriamo come superare le sfide speciali di multi-colore smFRET esperimenti sulle proteine. Vi presentiamo protocolli dettagliati per ottenere i dati e per l’estrazione di informazioni cinetiche. Questo include i criteri di selezione traccia, separazione di stato e il recupero delle traiettorie di stato dai dati rumorosi utilizzando un ensemble 3D modello di Markov nascosto (HMM). Rispetto ad altri metodi, le informazioni cinetiche non viene recuperate da istogrammi di tempo di permanenza, ma direttamente dall’HMM. Il quadro di massima verosimiglianza ci permette di valutare criticamente il modello cinetico e fornire significative incertezze per le tariffe.
Applicando il nostro metodo per la proteina di scossa di calore 90 (Hsp90), siamo in grado di districarsi tra l’associazione del nucleotide e globali cambiamenti conformazionali della proteina. Questo ci permette di osservare direttamente la cooperatività tra le due tasche di associazione del nucleotide del dimero Hsp90.
Molte proteine assolvere la loro funzione in dinamici complessi con altre molecole, mediati dai cambiamenti conformazionali e transitorie associazioni su una vasta gamma di scale cronologiche1,2,3. Accoppiato ad una fonte di energia esterna (ad es., ATP) queste interazioni dinamiche possono portare alla direzionalità in un ciclo funzionale e, infine, mantenere il non-equilibrio stazionario in una cella, il prerequisito per la vita.
Al fine di comprendere appieno queste macchine molecolari, una descrizione statica guidata da studi strutturali non è sufficiente. Inoltre, è essenziale avere conoscenza del modello cinetico e per determinare le costanti cinetiche tasso. Diversi metodi esistenti consentono ai ricercatori di studiare le dinamiche delle interazioni binarie tra due molecole di interesse, ad esempio, la risonanza plasmonica di superficie, metodi di rilassamento con una lettura spettroscopica (ad es., salto o flusso interrotto tecniche) e risonanza magnetica nucleare. Tuttavia, la loro applicabilità è nella maggior parte dei casi limitati a semplici sistemi di due stati (ad esempio, un limite e uno stato non associato) a causa della media inerente agli esperimenti di massa. Nei casi in cui siano coinvolti più Stati o intermedi, cedono solo una miscela complessa delle costanti di tasso. Metodi di singola molecola come pinzette ottiche o magnetiche o due colori smFRET, cioè, un donatore e un fluoroforo accettore, con un campione di superficie-immobilizzata possono recuperare le costanti di velocità per tutti osservati cambiamenti conformazionali. Tuttavia, quando si tratta di interazioni che interessano più di un binding sito, questi metodi rimangono limitati e le informazioni sulle interazioni possibili correlazioni di due (o più) saranno sempre accessibili via indirette conclusioni da una serie di esperimenti.
Multi-colore smFRET4,5,6,7,8,9 offre l’opportunità di studiare l’interazione tra questi componenti direttamente, in tempo reale e sotto condizioni quasi-fisiologiche10. Questo permette di indagare, ad esempio, l’associazione dipendente dalla conformazione di un ligando o di un’altra proteina8,9,11. L’approccio globale presentato qui è per etichettare la proteina di interesse alle posizioni specifiche, per allegare una proteina alla superficie della camera di misura e per monitorare l’intensità di fluorescenza nel tempo su un TIRFM di tipo a prisma (per dettagli Vedi 9 , 12). la prossimità spaziale delle tinture diverse quindi può essere determinata dal trasferimento di energia tra di loro. Strategie di etichettatura può variare da proteina a proteina (rivisto in 13) ed esistono linee guida per evitare artefatti nelle misure smFRET14.
Poiché un colorante donatore può trasferire energia al ricettore diversi coloranti in un esperimento di smFRET multi-colore, la posizione relativa di tutti i coloranti non è accessibile dall’eccitazione di uno tintura solo15,16. Ma in combinazione con alternanza di eccitazione laser (ALEX17e ha commentato in 18) questo metodo fornisce tutte le informazioni spazio-temporali a frazioni di secondo e risoluzione sub-nanometrica.
In linea di principio, ad alta risoluzione informazioni strutturali possono essere ottenute utilizzando le distanze di Inter-tintura calcolato dalla combinazione di tutte le intensità di fluorescenza in un esperimento di smFRET multi-colore con ALEX. Tuttavia, qui ci concentriamo sulla identificazione dello stato e separazione così come l’estrazione di modelli cinetici, dove è indispensabile smFRET multi-colore. Quando “solo” determinazione della struttura di triangolazione è voluta, una serie di esperimenti di smFRET di due colori più semplice con alto rapporto segnale-rumore può essere eseguita12,19.
Usiamo la fluorescenza parziale () come un proxy per il trasferimento di energia tra due fluorofori7. Il PF è calcolato dall’intensità di fluorescenza analoga per l’efficienza FRET di un esperimento di due colori:
Dove, è l’intensità di emissione canale em dopo l’eccitazione con colore ex, e c è il ricettore con la lunghezza d’onda più lunga. Canali di rilevazione rappresentano la stessa posizione nella camera del campione ma registrare diversi campi di spettro della luce di fluorescenza. Lo stesso identificatore per eccitazione e di emissione sono utilizzati nel presente protocollo (cioè, “blu”, “verde” e “rosso”).
A causa delle carenze sperimentali le intensità di fluorescenza misurata dipendono non solo sul trasferimento di energia, ma anche sulle proprietà fluoroforo e installazione. Per ottenere l’efficienza di trasferimento di vera energia tra due fluorofori, le intensità misurate devono essere corretto. La procedura seguente si basa sul riferimento9. Fattori di correzione per perdita apparente (LC, vale a dire, la rilevazione di fotoni da un fluoroforo in un canale destinata per un altro colorante) e gamma apparente (ag, vale a dire, il rendimento quantico di fluorescenza del colorante e la efficienza di rivelazione del canale) sono ottenuti da tracce di singola molecola che mostrano un accettore candeggio evento.
La dispersione del colorante donatore in ogni canale possibile acceptor viene calcolata da tutti i punti dati dalle tracce di fluorescenza registrata dove l’accettore sbiancato ma il donatore è ancora fluorescente ():
La mediana dell’istogramma perdita viene utilizzata come il fattore di perdita apparente. Dopo la correzione per la perdita, il fattore gamma apparente è determinato dallo stesso insieme di tracce. Viene calcolato dividendo il cambiamento di fluorescenza nel canale ricettore dal cambiamento di fluorescenza nel canale donatore con lo sbiancamento del colorante accettore:
Dove c è ancora una volta il canale di rilevamento per il ricettore con la lunghezza d’onda più lunga. La mediana della distribuzione risultante viene utilizzata come il fattore di correzione apparente.
La corretta intensità in ciascun canale sono stati ottenuti da:
Il PF è quindi calcolato in base alla:
Diverse popolazioni possono essere separati nello spazio multi-dimensionale, attraversato dal s PF. La posizione e la larghezza di ogni stato è determinata inserendo i dati con funzioni gaussiane multi-dimensionale. Successiva ottimizzazione di uno eh globale basato su tutte le tracce PF fornisce una descrizione quantitativa della cinetica osservata. Anche piccoli cambiamenti delle tariffe sono rilevabili.
HMMs consentono di dedurre un modello di stato da una raccolta di tracce tempo rumoroso. Il sistema è considerato di essere in uno di una serie di discreti, stati nascosti in qualsiasi dato momento e l’effettiva osservazione (cioè, l’emissione) è una funzione probabilistica di questo stato nascosto20. Nel caso di dati di smFRET TIRFM, l’emissione probabilità bio a stato possono essere modellati da funzioni di densità di probabilità gaussiana continuo. A intervalli di tempo discreti regolarmente spaziati, le transizioni da uno a un altro stato possono verificarsi secondo la probabilità di transizione che è tempo-invariante e dipende solo dallo stato corrente. La matrice di transizione A contiene queste probabilità di transizione unij tra tutti gli stati nascosti. La distribuzione di stato iniziale dà le probabilità relative allo stato
per il primo punto di tempo di una traccia di tempo. Utilizzando un approccio di massima verosimiglianza, questi parametri possono essere ottimizzati per meglio descrivere i dati con l’avanti-indietro e Baum-Welch algoritmi20,21. Questo produce gli stimatori di massima verosimiglianza (MLE). Infine, la sequenza di stato che più probabilmente prodotta la traiettoria delle osservazioni può essere dedotta con l’algoritmo di Viterbi. A differenza di altre analisi HMM del smFRET dati24,25,26 non usiamo il HMM come una mera “smoothing” dei dati ma il modello di stato cinetico l’Estratto dal set di dati senza la necessità per il montaggio di tempo di permanenza istogrammi27. Analisi HMM sono fatto con in-House scripts utilizzando Igor Pro. Implementazione del codice si basa su riferimento21. Forniamo un kit di software e dati esemplari sulla nostra pagina Web per seguire le sezioni 5 e 6 del presente protocollo (https://www.singlemolecule.uni-freiburg.de/software/3d-fret). Software completo è disponibile su richiesta.
Tempo punti nei dati con PF < -1 o PF > 2 in qualsiasi canale di rilevamento sono assegnate la probabilità di emissione minima per tutti gli Stati (10-200). Questo previene le transizioni artificiale in questi punti di dati.
I parametri per la probabilità di emissione sono ottenuti dal fit dell’istogramma 3D PF con funzioni gaussiane come descritto al punto 5.7. Questi parametri sono tenuti fissi durante l’ottimizzazione dell’eh.
Nell’approccio presentato, il vettore di distribuzione stato iniziale e la matrice di transizione vengono utilizzati globalmente per descrivere l’intero ensemble di tracce. Essi vengono aggiornati basato su tutti i N molecole dal set di dati secondo riferimento27.
Parametri di avvio per la distribuzione di stato iniziale sono determinati dalle proiezioni 2D dell’istogramma PF (punto 5.3) e le probabilità di transizione sono riportate a 0,05 fatta eccezione per le probabilità di rimanere nello stesso stato, che sono scelti tali che la probabilità di lasciare un certo stato viene normalizzata all’unità.
Un metodo di analisi di probabilità viene utilizzato per dare gli intervalli di confidenza (CIs) per tutte le transizione tariffe21,22, che servono come stime significative per loro incertezza. Per calcolare i limiti dell’IC per un tasso specifico, la probabilità di transizione di interesse è fissata su un valore diverso dal MLE. Questo produce la prova modello λ’. Un test di ratio (LR) probabilità del rischio dato il set di dati 0 viene eseguita secondo:
Il 95% di confidenza associato per il parametro è raggiunto quando LR supera 3.841, il quantile 95% di un X2-distribuzione con un grado di libertà22,23.
Il potere del metodo è dimostrato utilizzando la Hsp90. Questa proteina abbondante è trovata nei batteri e negli eucarioti ed è parte della risposta stress cellulare28. È un promettente bersaglio di farmaci nel trattamento di cancro29. Hsp90 è un omodimero con una tasca di legame del nucleotide nel dominio N-terminale di ogni subunità30. Può subire transizioni tra almeno due conformazioni globalmente distinti, uno chiuso e uno N-terminale a forma di V aperta, conformazione19,31,32. La natura dimerica direttamente solleva il problema dell’interazione tra i due siti di legame del nucleotide in Hsp90.
Di seguito, forniamo un protocollo passo-passo per l’acquisizione dati e l’analisi di un esperimento di tre colori smFRET il lievito Hsp90 e del nucleotide. L’associazione dipendente dalla conformazione di fluorescente contrassegnati AMP-PNP (AMP-PNP *, un analogo non hydrolyzable dell’ATP) viene analizzato. L’applicazione della procedura descritta permette lo studio di associazione del nucleotide e allo stesso tempo i cambiamenti conformazionali di Hsp90 e quindi rivela la cooperatività tra le due tasche di associazione del nucleotide di Hsp90.
Vi presentiamo la procedura sperimentale per ottenere dati di smFRET tre colori per un sistema complesso della proteina e una descrizione dettagliata dell’analisi di queste misurazioni. Questo approccio offre la possibilità unica di valutare direttamente la correlazione tra più siti di interazione o cambiamenti conformazionali.
Al fine di ottenere dati di singola molecola multi-colori adatti sulle proteine è importante eseguire misurazioni riproducibili a un livello basso rumore. Ciò può …
The authors have nothing to disclose.
Questo lavoro è finanziato dalla Fondazione di ricerca tedesca (INST 39/969-1) e il Consiglio europeo della ricerca attraverso l’accordo di sovvenzione del CER n. 681891.
Setup | |||
vibration-damped optical table | Newport, Irvine, CA, USA | RS2000 | |
OBIS 473nm LX 75mW LASER | Coherent Inc, Santa Clara, CA, USA | 1185052 | |
OBIS 532nm LS 50mW LASER | Coherent Inc, Santa Clara, CA, USA | 1261779 | |
OBIS 594nm LS 60mW LASER | Coherent Inc, Santa Clara, CA, USA | 1233470 | |
OBIS 637nm LX 140mW LASER | Coherent Inc, Santa Clara, CA, USA | 1196625 | |
laser control unit | Coherent Inc, Santa Clara, CA, USA | 1234465 | Scientific Remote |
aspheric telescope lenses | Thorlabs Inc, Newton, New Jersey, USA | d=25.4mm, f=50mm and f=100mm | |
CF ex1 | AHF analysentechnik AG, Tübingen, Germany | ZET 473/10 | cleanup filter excitation |
CF ex2 | AHF analysentechnik AG, Tübingen, Germany | ZET 532/10 | cleanup filter excitation |
CF ex3 | AHF analysentechnik AG, Tübingen, Germany | ZET 594/10 | cleanup filter excitation |
CF ex4 | Thorlabs Inc, Newton, New Jersey, USA | FL635-10 | cleanup filter excitation |
DM ex1 | AHF analysentechnik AG, Tübingen, Germany | ZQ594RDC | dichroic mirror excitation |
DM ex2 | AHF analysentechnik AG, Tübingen, Germany | 570DCXR | dichroic mirror excitation |
DM ex3 | AHF analysentechnik AG, Tübingen, Germany | ZQ491RDC | dichroic mirror excitation |
AOTFnC-Vis | AA Opto-Electronic, Orsay, France | ||
λ/4 plate | Thorlabs Inc, Newton, New Jersey, USA | AQWP05M-600 | |
CFI Apo TIRF 100x | Nikon Instruments Inc, Melville, NY, USA | high-NA objective | |
piezo focus positioner MIPOS 250 CAP | piezosystem jena GmbH, Jena, Germany | Piezo Controller NV 40/1 CLE | |
piezo stepper | Newport, Irvine, CA, USA | PZA12 | PZC200-KT NanoPZ Actuator Kit |
achromatic aspheric lenses | Qioptiq Photonics GmbH & Co. KG, Göttingen, Germany | G322-304-000 | d=50mm, f=200mm |
adjustable optical slit | Owis GmbH, Staufen i. Br., Germany | 27.160.1212 | max. aperture 12 x 12 mm |
DM det1 | AHF analysentechnik AG, Tübingen, Germany | T 600 LPXR | dichroic mirror detection |
DM det2 | AHF analysentechnik AG, Tübingen, Germany | H 560 LPXR superflat | dichroic mirror detection |
DM det3 | AHF analysentechnik AG, Tübingen, Germany | HC BS R635 | dichroic mirror detection |
BP det1 | AHF analysentechnik AG, Tübingen, Germany | 525/40 BrightLine HC | bandpass filter detection |
BP det2 | AHF analysentechnik AG, Tübingen, Germany | 586/20 BrightLine HC | bandpass filter detection |
BP det3 | AHF analysentechnik AG, Tübingen, Germany | 631/36 BrightLine HC | bandpass filter detection |
BP det4 | AHF analysentechnik AG, Tübingen, Germany | 700/75 ET Bandpass | bandpass filter detection |
optical shutters detection | Vincent Associates, Rochester, NY, USA | Uniblitz VS25S2T0 | |
EMCCD iXon Ultra 897 | Andor Technology Ltd, Belfast, Northern Ireland | ||
digital I/O card, PCIe-6535 | National Instruments, Austin, Texas, USA | ||
syringe pump | Harvard Apparatus, Holliston, MA, USA | PHD22/2000 | |
Name | Company | Catalog Number | Comments |
Flow chamber | |||
quartz slides | G. Finkenbeiner Inc, Waltham, MA, USA | Spectrosil2000, h=3mm | |
TEGADERM film | 3M Deutschland GmbH, Neuss, Germany | 1626W | 10 x 12cm |
spray adhesive | 3M Deutschland GmbH, Neuss, Germany | Photo Mount 050777 | |
glycerol | Carl Zeiss AG, Oberkochen, Germany | Immersol G | |
immersion oil | OLYMPUS EUROPA SE & CO. KG, Hamburg, Germany | IMMOIL-F30CC | |
prism | Vogelsberger Quarzglastechnik GmbH, Hauzenberg, Germany | Suprasil1 | |
aluminium prism holder | custom built | ||
hollow setscrews | Thorlabs Inc, Newton, New Jersey, USA | with custom drilling | |
Tygon S3 E-3603 tubing | neoLab Migge GmbH, Heidelberg, Germany | 2-4450 | ACF00001 |
PTFE tubing | Bohlender GmbH, Grünsfeld, Germany | S1810-08 | |
Name | Company | Catalog Number | Comments |
Sample | |||
yeast Hsp90 D61C, Q385C_biotin | UniProt ID P02829 | ||
Maleimide derivatives of Atto488, Atto550 | ATTO-TEC GmbH, Siegen, Germany | ||
AMP-PNP* | Jena Bioscience, Jena, Germany | γ-[(6-Aminohexyl)-imido]-AMP-PNP-Atto647N | |
Fluospheres | Thermo Fisher Scientific, Waltham, MA, USA | F8764 | amine-modified, 0.2 μm, yellow-green fluorescent |
Name | Company | Catalog Number | Comments |
Software | |||
Andor Solis | Andor Technology Ltd, Belfast, Northern Ireland | version 4.30 | |
LabVIEW | National Instruments, Austin, Texas, USA | version 2012, 32bit; misc. hardware control | |
MDS control software | AA Opto-Electronic, Orsay, France | version 2.03a | |
Coherent Connection | Coherent Inc, Santa Clara, CA, USA | version 3 | |
Igor Pro | WaveMetrics Inc, Portland, OR, USA | version 6.37 |