Usando três cores único-molécula FRET para estudar a correlação de interações da proteína

Resumo

Aqui, apresentamos um protocolo para obter dados de três cores de smFRET e sua análise com um ensemble 3D Hidden Markov modelo. Com esta abordagem, cientistas podem extrair informações cinéticas de sistemas complexos proteína, incluindo cooperatividade ou interações correlacionadas.

Resumo

Transferência de energia de ressonância do único-molécula Förster (smFRET) tornou-se uma técnica amplamente utilizada biofísica para estudar a dinâmica de biomoléculas. Para muitas máquinas moleculares em proteínas uma célula tem de agir em conjunto com parceiros de interação em um ciclo funcional a cumprir sua tarefa. A extensão de duas cores para cor multi smFRET possibilita a sonda simultaneamente mais de uma interação ou mudança conformacional. Isto não só acrescenta uma nova dimensão às experiências de smFRET, mas também oferece a possibilidade única de estudar diretamente a sequência de eventos e para detectar interações correlacionadas ao usar uma amostra de imobilizado e uma fluorescência de reflexão interna total microscópio (TIRFM). Portanto, multi cor smFRET é uma ferramenta versátil para estudar biomolecular complexos de forma quantitativa e em um detalhe anteriormente inatingível.

Aqui, demonstramos como superar os desafios especiais de multi-cor smFRET experimentos sobre proteínas. Apresentamos a protocolos detalhados para obtenção dos dados e para extrair as informações cinéticas. Isso inclui os critérios de seleção de rastreamento, separação do estado e a recuperação de trajetórias de estado a partir dos dados ruidosos usando um ensemble 3D Hidden modelo de Markov (HMM). Em comparação com outros métodos, as informações cinéticas não são recuperadas de histogramas de tempo de interrupção mas diretamente do HMM. O quadro de máxima verossimilhança nos permite avaliar criticamente o modelo cinético e prever as taxas significativas incertezas.

Através da aplicação de nosso método para a proteína de choque do calor 90 (Hsp90), somos capazes de separar a ligação de nucleotídeos e as globais mudanças conformacionais da proteína. Isso nos permite observar diretamente a cooperatividade entre os dois bolsos de vinculação do nucleotide do dímero Hsp90.

Introdução

Muitas proteínas cumprem sua função em dinâmicos complexos com outras moléculas, mediadas por mudanças conformacionais e associações transientes em uma ampla gama de escalas temporais^1,2,3. Acoplado a uma fonte de energia externa (por exemplo, ATP) essas interações dinâmicas podem levar a direcionalidade de um ciclo funcional e, finalmente, manter o estado de desequilíbrio constante em uma célula, o pré-requisito para a vida.

A fim de compreender plenamente essas máquinas moleculares, uma descrição estática guiada por estudos estruturais não é suficiente. Além disso, é essencial ter conhecimento do modelo cinético subjacente e determinar as constantes de velocidade cinética. Vários métodos existentes permitem aos pesquisadores estudar a dinâmica de interações de binárias entre duas moléculas de interesse, por exemplo, ressonância de plasmon de superfície, métodos de relaxamento com uma leitura espectroscópico (por exemplo, o salto ou fluxo parou técnicas) e ressonância magnética nuclear. No entanto, sua aplicabilidade é na maioria dos casos limitados a simples sistemas de dois Estados (por exemplo, um limite e um estado desacoplado) devido a média de inerente às experiências em massa. Em casos onde a mais Estados ou intermediários estão envolvidos, rendem apenas uma mistura complexa das constantes de velocidade. Único-molécula métodos tais como Pinça óptica ou magnética ou smFRET de duas cores, ou seja, um doador e um aceitador fluoróforo, com uma amostra de superfície-imobilizado podem recuperar as constantes de taxa para todos observadas mudanças conformacionais. No entanto, quando se trata de interações que afetam o sítio de ligação a mais de um, esses métodos permanecem limitados e as informações sobre as interações de possível correlação dos dois (ou mais) só será acessíveis via indirectas conclusões de um conjunto de experiências.

Multi cor smFRET^4,5,6,7,8,9 oferece a oportunidade de estudar a interação entre esses componentes directamente, em tempo real e sob condições fisiológicas perto de¹⁰. Isso permite que um para investigar, por exemplo, a vinculação dependente da conformação de um ligante ou outra proteína^8,9,11. A abordagem global apresentada aqui é para rotular as somáticas de interesse em posições específicas, para anexar uma proteína para a superfície da câmara de medição e para controlar a intensidade da fluorescência ao longo do tempo em um prisma-tipo TIRFM (para detalhes veja ^{9 , 12}). a proximidade espacial das tinturas diferentes então pode ser determinada a partir a transferência de energia entre eles. Estratégias de rotulagem pode variar de proteína a proteína (revista em ¹³) e orientações para evitar artefatos nas medições smFRET existem¹⁴.

Desde que um corante doador pode transferir energia para aceitador diferentes corantes em um experimento smFRET multi cor, a posição relativa de todos os corantes não é acessível a partir da excitação de um corante sozinho^15,16. Mas em combinação com a alternância de excitação do laser (ALEX¹⁷e revisto em ¹⁸) este método fornece todas as informações de espaço-temporal em segundo e resolução de sub nanômetros.

No principal, de alta resolução, informações estruturais podem ser conseguidas usando-se as distâncias de tintura inter calculadas a partir da combinação de todas as intensidades de fluorescência em um smFRET de cor multi experiência com ALEX. No entanto, aqui focamos na identificação de estado e de separação, bem como a extração de modelos cinéticos, onde smFRET multi cor é indispensável. Quando "só" determinação de estrutura por triangulação é desejada, um conjunto de experiências mais simples de duas cores smFRET com alta relação sinal-ruído pode ser realizada de^12,19.

Nós usamos a fluorescência parcial () como um proxy para a transferência de energia entre dois fluorophores⁷. A PF é calculado a partir da intensidade de fluorescência análoga à eficiência de traste de um experimento de duas cores:

Onde, é a intensidade de emissão canal em após excitação com cor ex, e c é o aceitador com o comprimento de onda mais longo. Canais de deteção representam a mesma posição na câmara de amostra, mas registros diferentes intervalos espectrais da luz da fluorescência. O mesmo identificador para excitação e emissão são utilizados neste protocolo (ou seja, "azul", "verde" e "vermelho").

Por causa de deficiências experimentais as intensidades de fluorescência medido dependem não só sobre a transferência de energia mas também Propriedades fluoróforo e instalação. Para obter a eficiência de transferência de energia verdadeira entre dois fluorophores, as intensidades medidas tem que ser corrigido. O procedimento a seguir baseia-se na referência⁹. Fatores de correção para vazamento aparente (LC, ou seja, a detecção de fótons de um fluoróforo num canal designada para outro corante) e aparente gama (ag, ou seja, o rendimento quântico de fluorescência do corante e a eficiência de deteção do canal) são obtidos a partir de vestígios de único-molécula que mostram um aceitador evento de branqueamento.

O escapamento do doador corante em todos os canais possíveis aceitador é calculado a partir de todos os pontos de dados nos vestígios de fluorescência gravada onde o corante aceitador branqueada, mas o doador é ainda fluorescente ():

A mediana do histograma escapamento é usada como o fator de vazamento aparente. Após a correção para detecção de fugas, o fator gama aparente é determinado a partir do mesmo conjunto de traços. É calculado dividindo-se a mudança da fluorescência no canal aceitador pela mudança de fluorescência no canal ao descoramento do corante aceitador de doador:

Onde c é novamente o canal de deteção para o aceitador com o comprimento de onda mais longo. A mediana da distribuição resultante é usada como o fator de correção aparente.

As intensidades corrigidas em cada canal são obtidas por:

A PF é então calculado de acordo com:

Diferentes populações podem ser separadas no espaço multidimensional ocupado pela s da PF. A posição e a largura de cada Estado é determinado pelo ajuste os dados com funções Gaussian multidimensionais. Otimização subsequente de um HMM global com base em todos os vestígios PF fornece uma descrição quantitativa de cinética observada. Mesmo as pequenas alterações das taxas são detectáveis.

HMMs fornecem uma maneira de se inferir um modelo de estado de uma coleção de traços de tempo barulhento. O sistema é considerado em parte de um conjunto de discreta, Estados ocultos em um determinado momento e a observação real (ou seja, a emissão) é uma função probabilística deste estado oculto²⁰. No caso dos dados de smFRET TIRFM, a emissão probabilidades b_eu por Estado, pode ser modelado por funções de densidade de probabilidade gaussiana contínua. Em pontos regularmente espaçados de tempo discreto, transições de um para outro Estado podem ocorrer de acordo com a probabilidade de transição que é invariante no tempo e só depende do estado atual. A matriz de transição A contém estas probabilidades de transição um_ij entre todos os Estados-Membros escondidos. A distribuição de estado inicial dá as probabilidades específicas de estado para o primeiro ponto de tempo de um rastreamento de tempo. Usando uma abordagem de máxima probabilidade, estes parâmetros podem ser otimizados para melhor descrevem os dados com o para diante-para trás e Baum-Welch algoritmos^20,21. Isso produz os estimadores de máxima verossimilhança (MLE). Finalmente, a sequência de estado que provavelmente produzido a trajetória de observações pode ser inferida com o algoritmo de Viterbi. Em contraste com outras analises HMM smFRET dados^24,25,26 não usamos o HMM como um mero "suavização" da dados, mas o modelo de estado cinético extraído do conjunto de dados sem a necessidade para caber o tempo de permanência histogramas de²⁷. HMM análise é feita com in-house scripts usando Igor Pro. Aplicação do código baseia-se na referência²¹. Nós fornecemos um kit de software e dados exemplares na nossa página da Web para seguir seções 5 e 6 do presente protocolo (https://www.singlemolecule.uni-freiburg.de/software/3d-fret). Software completo está disponível mediante pedido.

Pontos nos dados com PF < -1 ou PF > 2 em qualquer canal de deteção são atribuídos a probabilidade de emissão mínima para todos os Estados (10^-200). Isso impede que transições artificiais nestes pontos de dados.

Os parâmetros para as probabilidades de emissão são obtidos a partir do ajuste do histograma PF 3D com funções Gaussian conforme descrito na etapa 5.7. Esses parâmetros são mantidos fixos durante a otimização do HMM.

Na abordagem apresentada, o vetor de distribuição do estado inicial e a matriz de transição são utilizados mundialmente para descrever todo o conjunto de traços. Eles são atualizados com base em todas as moléculas de N do conjunto de dados de acordo com a referência²⁷.

Parâmetros para a distribuição de estado inicial são determinados de projeções 2D do histograma PF (etapa 5.3) e as probabilidades de transição são definidas como 0,05 com exceção as probabilidades de permanecer no mesmo estado, que são escolhidos tais que a probabilidade de sair de um determinado estado é normalizada à unidade.

Um método de criação de perfil de probabilidade é usado para dar intervalos de confiança (CIs) para todos os transição taxas^21,22, que servem como estimativas significativas para sua incerteza. Para calcular os limites da IC para uma taxa específica, a probabilidade de transição de juros é fixa para um valor diferente do MLE. Isto rende o teste modelo λ'. Um teste de ratio (LR) probabilidade da probabilidade dado o conjunto de dados 0 é realizada de acordo com:

A confiança de 95% vinculada para o parâmetro é alcançado quando LR excede 3,841, o Quantil de 95% de um ^X2-distribuição com um grau de liberdade^22,23.

O poder do método é demonstrado usando o Hsp90. Esta proteína abundante é encontrada em bactérias e eucariotas e é parte do estresse celular resposta²⁸. É um alvo de droga promissora no tratamento de câncer²⁹. Hsp90 é um homodímero com um bolso de ligação de nucleotídeos do domínio N-terminal de cada subunidade³⁰. Ele pode sofrer transições pelo menos duas conformações distintas globalmente, um fechado e um N-terminal conformação aberta, em forma de V^19,31,32. A natureza dimérica diretamente a questão da interação entre os dois locais de ligação de nucleotídeos em Hsp90.

A seguir, fornecemos um protocolo passo a passo para a aquisição de dados e análise de um experimento de smFRET três cores na levedura Hsp90 e nucleotídeos. A vinculação dependente da conformação do AMP-PNP fluorescente etiquetadas (AMP-PNP *, um analog não hidrolisável do ATP) é analisada. A aplicação do procedimento descrito permite o estudo da ligação de nucleótidos e ao mesmo tempo as mudanças conformacionais de Hsp90 e, assim, revela a cooperatividade entre os dois bolsos de vinculação de nucleotídeo de Hsp90.

Protocolo

1. instalação e pré-requisitos

1. Realizar as medições de cor multi smFRET em um TIRFM de prisma-tipo. Uma descrição de uma configuração de duas cores como uma publicação de JoVE é dado referência a¹².
   1. Construa um TIRFM de várias cor. Um layout geral é detalhado em ⁹.
   2. Uso comutável, diodo bombeou onda contínua lasers de estado sólido, que tornam a utilização de obturadores mecânicos nos caminhos de excitação desnecessárias.
   3. Emprega um prisma assimétrico, alongado que impede a reflexão traseira do feixe da excitação da parte traseira para entrar o objetivo.
   4. Use lentes de silicone acromática de esféricos fundido de 2 polegadas nos caminhos de deteção que coletar tanta luz quanto possível e evitar autofluorescência e aberrações, por exemplo, as distorções nas regiões fora do centro da imagem.
   5. Focar cada caminho de deteção do chip do EMCCD com uma lente separada. Isso permite que o foco ideal de cada canal de deteção.
      Atenção: Lasers de classe 3B são usados na TIRFM. Isto significa que eles são perigosos, se o olho é exposto diretamente, mas reflexões difusas não são prejudiciais. Garantir a conformidade com laser precauções de segurança, de acordo com as regulamentações governamentais locais antes que o sistema é operado.
2. Determine os factores de correcção para propriedades de configuração e fluoróforo previamente usando amostras de dsDNA.
   1. Usar um alto-FRET dsDNA amostra para cada corante em combinação com o aceitador tendo o mais longo comprimento de onda de excitação (para a configuração apresentada: Atto488-Atto647N, Atto550-Atto647N, Atto594-Atto647N). Certifique-se de que o DNA é modificado adicionalmente com uma biotina.
   2. Diluir a amostra a 5 nM com tampão TNM (5 mM Tris pH 7,5, 5 mM de NaCl, 20 mM MgCl₂) e 2 mM Trolox (use também esse buffer para medição).
   3. Imobilize o dsDNA conforme descrito para Hsp90 em etapas 2.5 e 2.7.
   4. Calcule os fatores de correção para vazamento aparente (LC) e aparente gama (ag) de traços de único-molécula que mostram um evento de branqueamento do aceitador.
3. Construa uma câmara de fluxo que é um sanduíche de uma lâmina de quartzo passivadas PEG/biotina-PEG, uma película fina que é adesivo em ambos os lados e uma lamínula. Para protocolos de detalhe para a limpeza do quartzo slides e passivação consulte referência⁹.
   1. Uso de espessura (3mm) quartzo slides para evitar geometricamente a coleção do laser espalhados na interface entre o prisma e o slide funcionalizados quartzo quartzo-glicerol-quartzo.
   2. Use uma fina (40 µm) filme que é pulverizado com adesivo no lado não-adesivo de vedação. A película fina reduz a distância entre as moléculas de superfície-anexado e objectivo. Aquecer a 80 ° C e pressione.
   3. Coloca uma lamínula por cima. Aquecer a 80 ° C e pressione. Use uma gota de glicerol, ao colocar a câmara de fluxo sobre o prisma.
      Nota: Os materiais de prisma e os slides de quartzo, bem como o glicerol são índice de correspondência.
4. Express Hsp90 de Saccharomyces cerevisiae na forma de dois mutantes de ponto único de cisteína nas posições D61 ou Q385. Adicione um motivo de bobina coiled do C-terminal para evitar a dissociação do dímero em concentrações picomolar. Rotular as proteínas mutantes separadamente e trocar os monómeros para obter heterodímeros rotulados com Atto488 na posição do ácido aminado 61 e Atto550 no ácido aminado posição 385⁹.

2. medição

1. Inicie o software da câmera e definir os parâmetros de imagem conforme especificado abaixo:
   1. Regule a temperatura para o sensor de resfriamento mais baixo possível (-95 ° C com refrigeração a água externa) para diminuir o ruído de corrente escuro.
   2. Usar as configurações da câmera que são otimizadas para a gravação do único-molécula: velocidade de deslocamento vertical 3.3 µs, tensão de pulso de disparo vertical normal, 17 horizontal MHz 16 bits lidos, pré-amplificação ganho 3, ganho de elétron multiplicador 1.000.
   3. Defina o acionamento do alarme a aquisição "Externo" e o tempo de exposição para 70 filmes de registro MS., com um comprimento de 750 ciclos de aquisição.
      Nota: Desliga a luz do quarto quando a câmera está adquirindo para evitar a saturação do sensor EMCCD.
2. Crie uma pasta no local SSD para a medição. No software configurações vão para o piloto < Auto-Save > habilitar < Auto-save > em escolher o formato de arquivo "Tiff" para aquisição de filme. Selecione a pasta sobre o SSD como local de salvamento automático.
3. Inicie o software que controles o acústico-Óptico Sintonizável filtro (AOTF), o software que controla a operação dos lasers e o software de gatilho que sincroniza os lasers, AOTF, persianas no caminho de deteção e câmeras. Ajustar a potência do laser com o AOTF (ca. 3 mW antes de entrar o prisma) e carregar o padrão correto de disparo.
4. Montar o porta-amostras com o prisma e a câmara de fluxo, conecte a tubulação, coloque o tubo de entrada em um copo de microcentrifuga e conecte o tubo de saída de uma bomba de seringa. Lave a câmara com cerca de 150 buffer µ l, alinhe a excitação feixe e foco. Bleach todos os contaminadores na superfície fluorescentes por movendo-se lentamente ao longo do intervalo de detecção completa do slide com um poder do laser de cerca de 10 mW para todos os lasers. Isso leva cerca de 1h.
   Nota: se não for indicado o contrário, o buffer usado contém 40 mM HEPES pH 7,5, KCl, de 150 mM e 10 mM MgCl₂.
5. Descarga de aproximadamente 300 µ l de uma solução de NeutrAvidin (0,25 mg/mL em buffer) para a câmara e incubar durante 1 min. Flush fora NeutrAvidin desacoplado com tampão e flush aproximadamente 300 ml de uma solução de BSA (0,5 mg/mL em buffer) através da câmara.
   Nota: Este bloco restante superfície functionalization defeitos por adsorção inespecíficas da BSA para a superfície.
6. Imobilizar a amostra por carregar a câmara de fluxo com aproximadamente 150 µ l do biotinilado e rotulado Hsp90 em concentrações crescentes (diluídas em tampão + 0,5 mg/mL BSA) até uma densidade de superfície suficiente é atingida, que é geralmente o caso numa concentração de 5-22:00. Lave desvinculado da proteína com buffer de aproximadamente 300 µ l + 0,5 mg/mL BSA.
7. Descarga em 150 µ l de 25 nM AMP-PNP * no tampão + 0,5 mg/mL BSA. Deixe incubar durante 5 min e repita este passo uma vez para garantir a concentração correta do nucleotide.
   Nota: Para experiências na presença de AMP-PNP adicionais, sem rótulo, também adicionar 250 µM AMP-PNP.
8. Começam com a aquisição de dados. Para uma quantidade adequada de dados adquira cerca de 20 filmes, o que leva cerca de 1,5 h.
   1. Mova a posição da câmara de amostra perpendicular ao feixe de excitação com um seletor de piezo para alterar o campo de visão.
   2. Ajuste o foco com o z-piezo que controla a altura do objectivo, se necessário. Isso não deve ser necessário muitas vezes, quando a câmara de medição é montada sem inclinação.
   3. Preparar a gravação das câmeras pressionando < leve sinal > no software da câmera e iniciar os ciclos de excitação/aquisição usando o botão < Iniciar > no software do gatilho. Isso inicia a aquisição da intensidade da fluorescência.
9. Execute um registo de canal pela primeira gravação de um filme com grânulos fluorescentes que mostram a emissão de fluorescência na faixa espectral de todos os canais de deteção da instalação. Em seguida, detectar as posições do grânulo do filme de calibração por procurar os pontos mais brilhantes e determinar que a posição central de uma gaussiana se encaixam no perfil de intensidade. Salve as coordenadas dos grânulos que são encontrados em todos os canais e ajuste os dois o deslocamento mapeamento em x - e y-direção com uma 2D polinomial de grau três⁹.

3. selecção dos traços de único-molécula

1. A análise dos dados é feita com in-house scripts usando Igor Pro. Carregar todos os scripts necessários através da abertura de "iniTIRF.ipf" e selecione o tipo correto de experiência.
   Nota: no seguinte, < botão > especifica elementos clicáveis no menu ou a interface do usuário. Chamadas de função são indicadas entre aspas, por exemplo, "Print"Hello World"". Esses comandos podem ser colados na linha de comando do Igor Pro (sem as aspas de fechamento).
2. Verifique se que os parâmetros para o registo do canal de deteção são carregados.
3. Inicie a GUI clicando < smFRET novo | Análise GUI >.
4. Carrega os filmes (ou seja, a sequência de quadros com 512 x 512 pixels armazenados como pilhas TIFF de 16 bits) que seguram a intensidade nos respectivos canais. Fazer isso pressionando o botão < filme de carga > e selecionando os arquivos das câmeras "mestre" e "escravo" so-called um após o outro.
5. Identifica as posições dos potenciais moléculas simples, procurando os pontos mais brilhantes na soma dos cinco primeiros quadros em um determinado canal de deteção. Calcule as posições correspondentes em todos os outros canais de deteção de mapeamento do canal. Para obter o rastreamento de intensidade de fluorescência, soma-se a intensidade de um quadrado de pixel em torno da posição central para cada quadro. Fazer isso pressionando o botão < encontrar vestígios > no GUI.
   Nota: O comprimento do lado do quadrado (em pixel) é dada por: 2 * < soma Pxs > + 1.
6. Para cada molécula, calcule um traço comum intensidade crua como a soma de todos os vestígios deste local com a mesma cor de excitação. Avalie o perfil de intensidade da molécula em todos os canais para os seguintes critérios:
   1. Um planalto aproximadamente plano em intensidade crua conjunta e um branqueamento único passo para todas as cores de excitação, comportamento anticorrelated nos canais apropriados de deteção, detecção de fluorescência vermelha (indicando FRET para um limite AMP-PNP *) pelo menos uma vez dentro do rastreamento e não há várias etapas na fluorescência vermelha, que indicaria a presença de dois AMP-PNP * acoplado a um dímero de Hsp90.
7. Salve os traços de fluorescência do local para posterior análise, se estes critérios forem satisfeitos. Fazer isso selecionando o rastreamento com o cursor e pressionando o botão < salvar > no gráfico de "prazos". Inspecione manualmente os traços de intensidade em todos os três canais de deteção após excitação azul para cerca de 200 moléculas por filme.

4. cálculo dos traços fluorescência parcial

1. Exiba todos os vestígios de intensidade de fluorescência para uma das moléculas salvas. Use os cursores no gráfico para selecionar intervalos de tempo.
2. Selecione um intervalo de tempo onde todos os fluorophores são clareados já. A intensidade média do fundo é calculada a partir desse intervalo e subtraída o rastreamento de intensidade em cada canal. Fazer isso pressionando o botão < Background >.
3. Selecione o intervalo de eficiência FRET, onde pelo menos ambos os corantes anexados a Hsp90 (Atto488 e Atto550) estão presentes. Certifique-se de excluir rastreamentos que contêm um evento piscando (ver Figura 3B). Esses eventos são caracterizados por uma intensidade de fluorescência de drop-in em um canal sem um aumento que acompanha qualquer outro canal.
4. Calcule os vestígios PF . Fazer isso pressionando o botão < PF Calc >. Os factores de correcção predefinidos para vazamento aparente (LC) e aparente gama (ag) são aplicados a intensidade crua para corrigir para foto-física e propriedades de configuração.

5. população seleção e ajuste de histograma 3D

1. Remova moléculas que mostram uma baixa relação sinal-ruído nos vestígios PF . Moléculas que excedem o intervalo [-1; 2] em qualquer vestígio PF para mais de 10% dos quadros são removidas do conjunto de dados. Fazer isso, executar o "RemoveTracesLowSNR()" na janela de comando.
2. Calcule binned projeções 2D de dados da PF . Lote sobre e sobre no intervalo [-0.5; 1.5] com uma resolução de 100 x 100 caixas. Para fazer isso, execute:
   1. "HistFret2D ("r_b","r_g", binHist = 100)"
   2. "HistFret2D ("r_b","g_b", binHist = 100); MoveWindow 553.5, 42,5, 1055.25, 508,25"
3. Determine a população relativa de cada Estado distinguível nas projecções 2D.
   1. Trazer o gráfico apropriado para a frente e executar o "panelHist2DCount()".
   2. Pressione a tecla < Init > e desenhar um polígono de mão livre em volta do pico.
   3. Clique no botão < contagem >. O número de pontos de dados no polígono e o número total de pontos de dados na projeção são impressos na janela de comando.
4. Preparar um histograma 3D dos dados PF executando "HistFret3D ("g_b","r_b","r_g")".
5. Normalize o histograma 3D para uma integral de 1. Execute o seguinte:
   1. "Fit0 de NewDataFolder/S"
   2. "Duplicata/O:: FRET: Hist3D, Hist3D"
   3. "Div variável/g = sum(Hist3D)*(DimDelta(Hist3D,0)) ^ 3"
   4. "Hist3D = div; Impressão div"
6. Fornece parâmetros iniciais para o Gaussian 3D se encaixa e preparem estruturas de dados necessárias.
   1. Execute "Gauss3D_initParam(); editar W_coef_old."
   2. Adicione as populações do estado até o fim do vetor parâmetro.
   3. Executar "Gauss3D_prepareFit()."
      Nota: W_coef_old é um vetor que contém os parâmetros iniciais para o ajuste. Por Estado, isso significa e a população do estado, que é concatenada ao final do vetor. Certifique-se que a matriz de covariância é simétrica.
7. Cabe a soma de S funções Gaussian 3D para o histograma PF 3D, com S , sendo o número de Estados distinguíveis.
   1. Executar "do3D()." Isto pode demorar uma hora ou mais em um escritório normal PC, dependendo da qualidade dos parâmetros iniciais.
   2. Execute "postprocessFitMultiGauss3D(); evalFitMultiGauss3D(); editar W_coef."
8. Exiba o resultado de apto. Para cada uma das duas projeções 2D, use os seguintes comandos:
   1. "contourPF3D_new(0); contourPF3D_new(1); contourPF3D_new(2); contourPF3D_new(3); contourPF3D_new(4)"
   2. "contourPF3D_colorize()"

6. cinética análise 3D Ensemble HMM

1. Prepare um ensemble HMM executar para extrair as informações cinéticas. Um HMM é otimizado de todas as moléculas no conjunto de dados. Use as informações obtidas na etapa anterior para definir a posição e a largura de cada Estado no espaço 3D PF .
   1. Inicializar a interface do usuário HMM (< HMM | Init HMM >) e escolha o número apropriado de Estados-Membros (no caso dos dados Hsp90 isso significa < NumStates > = 5), o número de dimensões do sinal de entrada (< NumDims > = 3) e o tipo da entrada (< Input Type > = "FRET bgr 3D").
2. Otimizar os parâmetros do HMM, deixando o software convergem a probabilidade do HMM (execute "prepENS_CONVERGE_gB(GetDataFolder(1), -14)") até que a mudança na matriz de transição, em comparação com a iteração anterior cai abaixo de um limite (10^-14 para a soma da mudança absoluta para cada probabilidade de transição). Isto rende o MLE para as probabilidades de transição em cerca de uma hora em um PC normal do escritório.
3. Repita a seleção da população, encaixe gaussiana e a otimização de HMM para subconjuntos de dados (por exemplo, 75% do conjunto de dados completo). Se fundiu-se o conjunto de dados completo de experiências diferentes, repita a otimização também para cada um dos experimentos simples. Análise de subconjuntos permite estimar a incerteza da seleção manual da população e a variabilidade dentro do conjunto de dados.
4. Calcule o CI para as probabilidades de transição, que o relatório sobre a heterogeneidade do conjunto de dados e a precisão do HMM.
   1. Obter uma estimativa aproximada dos limites CI executando "$cd (raiz: path3Dimport +"HMM"); loop_getCI_estimate_limits()."
   2. Calcule os limites exatos da IC, executando:
   3. "loop_getCI_HMM_converge(1)"
   4. "CIresults_conv_new()"
   5. "cd:: HMM_CIresult; reportCI_conv()"
   6. "cd:: cmp_CI_conv; CI_plot2 ("HMM", doAppend = 0) "
5. Condense as informações cinéticas disponíveis para simplificar a interpretação.
   1. Colete informações sobre o tempo um nucleotídeo rotulado permanece ligado à Hsp90. Para fazer isso, o colapso dos Estados que são obrigados a AMP-PNP * (ou seja, S₁e S₂ , S₃, ver também a Figura 2) e compilar o histograma de tempo de permanência. Execute o "cd $(raiz: path3Dimport +"HMM"); collapse_states_get_DT({0,1,1,1,0}) ", que combina habita de Estados S₀ e S₄ bem como S₁e S₂S₃.
   2. Extrair os tempos de permanência do caminho de Viterbi para cada Estado de interesse e comparar o histograma resultante de tempo Habitai para diferentes condições experimentais. Execute "plot_collapsed_DT_Hist (wDTo_01110_record1)."
   3. Para uma imagem mais detalhada, colapso dos Estados que têm distintas PF mas são funcionalmente idênticos para facilitar a análise de dados adicional, por exemplo, no caso de um experimento de três cores com Hsp90 etiquetados e nucleotídeos, Estados S₂ e S ₃ pode ser recolhido. Executar "collapse_states_get_DT({0,1,2,2,4})."

Resultados

Medições de cor multi smFRET permitem a detecção directa de correlação entre dois ou mais sites de interação distintos. Isso torna a técnica exclusiva para investigar sistemas multicomponentes, tais como os complexos da proteína. Focalizamos a apresentação de um experimento de três cores de smFRET aqui, que serve como um exemplo ilustrativo.

O fluxo de trabalho geral do método é mostrado na Figura 1...

Discussão

Apresentamos o procedimento experimental para obter dados de três cores smFRET para um sistema complexo de proteínas e uma descrição passo a passo da análise dessas medições. Essa abordagem oferece a possibilidade única de avaliar diretamente a correlação entre vários sites de interação ou mudanças conformacionais.

A fim de obter dados de único-molécula multi cores adequados em proteínas é importante realizar as medições podem ser reproduzidas em um baixo nível de ruído. Isto pode ser consegui...

Materiais

Lista de materiais utilizados neste artigo

Nome	Empresa	Número de catálogo	Comentários
Configuração
mesa óptica amortecida por vibração	Newport, Irvine, CA, EUA	RS2000
OBIS 473nm LX 75mW LASER	Coherent Inc, Santa Clara, CA, EUA	1185052
OBIS 532nm LS 50mW LASER	Coherent Inc, Santa Clara, CA, EUA	1261779
OBIS 594nm LS 60mW LASER	Coherent Inc, Santa Clara, CA, EUA	1233470
OBIS 637nm LX 140mW LASER	Coherent Inc, Santa Clara, CA, EUA	1196625
unidade de controle a laser	Coherent Inc, Santa Clara, CA, EUA	1234465	Scientific Remote
lentes de telescópio asférico	Thorlabs Inc, Newton, New Jersey, EUA		d=25.4mm, f=50mm e f=100mm
CF ex1	AHF analysentechnik AG, Tü bingen, Alemanha	ZET 473/10	limpeza excitação do filtro
CF ex2	AHF analysentechnik AG, Tü bingen, Alemanha	ZET 532/10	limpeza excitação do filtro
CF ex3	AHF analysentechnik AG, Tü bingen, Alemanha	ZET 594/10	excitação do filtro de limpeza
CF ex4	Thorlabs Inc, Newton, New Jersey, EUA	FL635-10	excitação do filtro de limpeza
DM ex1	AHF analysentechnik AG, Tü bingen, Alemanha	ZQ594RDC	excitação dicróica do espelho
DM ex2	AHF analysentechnik AG, Tü bingen, Alemanha	570DCXR	excitação de espelho dicróico
DM ex3	AHF analysentechnik AG, Tü bingen, Alemanha	ZQ491RDC	excitação de espelho dicróico
AOTFnC-Vis	AA Opto-Electronic, Orsay, França
λ/4 placa	Thorlabs Inc, Newton, New Jersey, EUA	AQWP05M-600
CFI Apo TIRF 100x	Nikon Instruments Inc, Melville, NY, EUA		de alta objetiva
NA MIPOS 250 CAP	piezosystem jena GmbH, Jena, Controlador		NV 40/1 CLE
passo piezo	Newport, Irvine, CA, EUA	PZA12	PZC200-KT Kit de atuador NanoPZ
lentes asféricas acromáticas	Qioptiq Photonics GmbH & Co. KG, Gö ttingen, Alemanha	G322-304-000	d=50mm, f=200mm
fenda óptica ajustável	Owis GmbH, Staufen i. Br., Alemanha	27.160.1212	abertura máxima 12 x 12 mm
DM det1	AHF analysentechnik AG, Tü bingen, Alemanha	T 600 LPXR	detecção de espelho dicróico
DM det2	AHF analysentechnik AG, Tü bingen, Alemanha	H 560 LPXR detecção	espelho dicróico
DM det3	AHF analysentechnik AG, Tü bingen, Alemanha	HC BS R635	detecção de espelho dicróico
BP det1	AHF analysentechnik AG, Tü bingen, Alemanha	525/40 BrightLine HC	detecção de filtro passa-banda
BP det2	AHF analysentechnik AG, Tü bingen, Alemanha	586/20 BrightLine HC	detecção de filtro passa-banda
BP det3	AHF analysentechnik AG, Tü bingen, Alemanha	631/36 BrightLine HC	detecção de filtro passa-banda
BP det4	AHF analysentechnik AG, Tü bingen, Alemanha	700/75 ET Detecção	de filtro passa-banda de passagem
de banda detecção de persianas ópticas	Vincent Associates, Rochester, NY, EUA	Uniblitz VS25S2T0
EMCCD iXon Ultra 897	Andor Technology Ltd, Belfast, Irlanda
placa de E/S digital, PCIe-6535	National Instruments, Austin, Texas, EUA
seringa	Harvard Apparatus, Holliston, MA, EUA	PHD22/2000
Nome	Empresa<	strong>Número de catálogo	Comentários
Lâminas de câmara de fluxo
quartz	G. Finkenbeiner Inc, Waltham, MA, EUA		Spectrosil2000, h=3mm
Filme TEGADERM	3M Deutschland GmbH, Neuss, Alemanha	1626W 10	x 12cm
adesivo spray	3M Deutschland GmbH, Neuss, Alemanha	Montagem para fotos 050777
glicerol	Carl Zeiss AG, Oberkochen, Alemanha	Immersol G
	OLYMPUS EUROPA SE & CO. KG, Hamburgo, Alemanha	Prisma IMMOIL-F30CC
	Vogelsberger Quarzglastechnik GmbH, Hauzenberg, Alemanha		Suprasil1
			de
de fixação ocos	Thorlabs Inc, Newton, New Jersey, EUA		com perfuração personalizada
Tygon S3 E-3603 tubulação	neoLab Migge GmbH, Heidelberg, Alemanha	2-4450	ACF00001
tubulação PTFE	Bohlender GmbH, Grü nsfeld, Alemanha	S1810-08
Nome	Empresa	Número de catálogo	Comentários<
strong>Amostra
yeast Hsp90 D61C, Q385C_biotin			UniProt ID P02829
Maleimide derivados de Atto488, Atto550	ATTO-TEC GmbH, Siegen, Alemanha
AMP-PNP*	Jena Bioscience, Jena, Alemanha		γ-[(6-Aminohexyl)-imido]-AMP-PNP-Atto647N
Fluosferas	Thermo Fisher Scientific, Waltham, MA, EUA	F8764	modificado por amina, 0,2 μ m, amarelo-esverdeado fluorescente
Nome<	strong>Empresa<	strong>Número de catálogo<	strong>Comentários
Software
Andor Solis	Andor Technology Ltd, Belfast, Irlanda		versão 4.30
LabVIEW	National Instruments, Austin, Texas, EUA		versão 2012, 32 bits; controle de hardware diverso
Software de controle MDS	AA Opto-Electronic, Orsay, França		versão 2.03a
Coherent Connection	Coherent Inc, Santa Clara, CA, EUA		versão 3
Igor Pro	WaveMetrics Inc, Portland, OR, EUA		versão 6.37