Uso de equações lineares na determinação dos parâmetros de Michaelis-Menten

Use of linear equations to obtain Michaelis-Menten parameters

Michaelis-Menten; Lineweaver-Burk; nonlinear least-square fit

a33v33n7

EDUCAÇÃO

Uso de equações lineares na determinação dos parâmetros de Michaelis-Menten

Use of linear equations to obtain Michaelis-Menten parameters

Nakédia M. F. Carvalho^I; Bianca M. Pires^I; Octavio A. C. Antunes^I,†; e Roberto B. Faria^I,*; Renata E. H. M. B. Osório^II; Clovis Piovezan^II; Ademir Neves^II

^IInstituto de Química, Universidade Federal do Rio de Janeiro, Av. Athos da Silveira Ramos, 149, CT Bl. A, 21941-909 Rio de Janeiro - RJ, Brasil

^IIUniversidade Federal de Santa Catarina, Departamento de Química, Campus Trindade, CP 476, 88040-900 Florianópolis - SC, Brasil

ABSTRACT

The Michaelis-Menten equation is used in many biochemical and bioinorganic kinetic studies involving homogeneous catalysis. Otherwise, it is known that determination of Michaelis-Menten parameters K_M, V_max, and k_cat by the well-known Lineweaver-Burk double reciprocal linear equation does not produce the best values for these parameters. In this paper we present a discussion on different linear equations which can be used to calculate these parameters and we compare their results with the values obtained by the more reliable nonlinear least-square fit.

Keywords: Michaelis-Menten; Lineweaver-Burk; nonlinear least-square fit.

INTRODUÇÃO

No estudo de reações enzimáticas e biomiméticas utiliza-se comumente o esquema de Briggs-Haldane,¹ que considera as reações de formação e dissociação do complexo enzima-substrato, ES, e sua reação de decomposição formando os produtos e regenerando a enzima.

Aplicando-se a este esquema cinético a aproximação do estado estacionário (aproximação de Bodenstein) para a concentração do complexo enzima-substrato obtemos a famosa equação de hipérbole retangular de Michaelis-Menten para a velocidade inicial da reação (Equação 1).^2-5

onde, K_M = (k + k₂)/k₁ é a constante de Michaelis e k₂[E]₀ = V_max é a velocidade máxima da reação para uma dada concentração inicial da enzima, [E]₀, na presença de um excesso de substrato. A constante de velocidade k₂ é também chamada de k_cat ou número de renovação (turnover number, turnover frequency).

Para se obter K_M, V_max e k_cat, o melhor procedimento é o uso de um ajuste não linear por mínimos quadrados (nonlinear least-square method) dos dados experimentais de v₀ contra [S]₀.⁶ Signorella et al.⁷ mostraram recentemente que k_cat é especialmente sensível ao método de ajuste empregado, enfatizando que o ajuste não linear por mínimos quadrados é o mais confiável. Por outro lado, um ajuste linear menos confiável do ponto de vista estatístico pode ser feito usando-se a equação de duplo recíproco de Lineweaver-Burk (Equação 2),⁸ que pode ser obtida após manipulação da Equação 1.

Assim, um gráfico de 1/v₀ contra 1/[S]₀ permite a determinação de 1/V_max como a intercessão da reta com o eixo das ordenadas (coeficiente linear) e K_M/V_max pela inclinação da reta (coeficiente angular), a partir dos quais K_M e V_max podem ser obtidos.

Apesar do fato de que diversos pacotes computacionais atuais oferecem a possibilidade de um ajuste não linear por mínimos quadrados, um ajuste linear é mais simples de ser apresentado e compreendido pelos estudantes, bem como, geralmente, permite uma percepção mais imediata dos pontos experimentais que apresentem maiores desvios de um comportamento do tipo Michaelis-Menten (quando a maioria dos pontos não segue um ajuste linear, isto indica que podemos estar diante de outro tipo de mecanismo onde ocorram, como por exemplo, etapas de inibição enzimática, um aspecto que não será abordado no presente artigo). Entretanto, como é sabido que o uso da Equação 2 não produz valores confiáveis para os parâmetros de Michaelis-Menten,^2,4,7 apresentamos neste trabalho uma discussão sobre quatro diferentes equações lineares que podem ser obtidas por simples manipulação da Equação 1, com o objetivo de escolher qual delas produz os melhores parâmetros. A abordagem que escolhemos para este estudo foi a de comparar os parâmetros de Michaelis-Menten obtidos por cada uma das equações lineares com aqueles obtidos pelo ajuste não linear por mínimos quadrados, que é o mais confiável. Os resultados obtidos mostram que duas dessas equações lineares, uma delas aqui apresentada pela primeira vez, podem ser utilizadas para o cálculo dos parâmetros de Michaelis-Menten, uma vez que os valores obtidos são muito próximos daqueles obtidos por ajuste não linear por mínimos quadrados e significativamente mais confiáveis do que os obtidos pela Equação de Lineweaver-Burk.

MÉTODOS

Os ajustes não lineares por mínimos quadrados dos dados experimentais foram feitos empregando-se o pacote computacional Origin 6.1.⁹ Os ajuste lineares (regressão linear) foram obtidos empregando-se a planilha eletrônica Lotus 1-2-3.¹⁰

Os dados experimentais utilizados são provenientes de estudos cinéticos atualmente sendo feitos em nossos laboratórios envolvendo a síntese de compostos bioinorgânicos miméticos de enzimas naturais. O conjunto de dados 1 (ver Tabela 1) refere-se ao acompanhamento cinético da atividade de catalase (reação de decomposição do peróxido de hidrogênio) de um complexo tetranuclear de ferro(III) com o ligante 1,3-bis[((2-aminoetil)amino]-2-propanol). Os conjuntos de dados 2 e 4 (ver Tabela 3) referem-se à hidrólise do 2,4-bis-dinitro-fenilfostato e do 3,5-di-terc-butilcatecol por um complexo dinuclear de Cu e o conjunto de dados 3 refere-se à oxidação do 3,5-di-terc-butilcatecol catalisada por um complexo dinuclear de Ni.

RESULTADOS E DISCUSSÃO

Diferentes equações lineares podem ser obtidas por meio de simples manipulações da Equação 1, dentre as quais temos as equações de Lineweaver-Burk (Equação 2),⁸ de Hanes (Equação 3)¹¹ e de Eadie-Hofstee (Equação 4).^12,13 A estas podemos acrescentar a Equação 5 que, no melhor do nosso conhecimento, ainda não foi investigada. A partir destas equações os parâmetros de Michaelis-Menten podem ser facilmente obtidos.

A Equação 2 permite a determinação de V_max e K_M a partir de um gráfico de 1/v₀ contra 1/[S]₀, a Equação 3 a partir de um gráfico de [S]₀/v₀ contra [S]₀, a Equação 4 a partir de um gráfico de v₀ contra v₀/[S]₀, e a Equação 5 a partir de um gráfico de [S]₀ contra [S]₀/v₀. O uso das Equações 4 e 5 parece mais simples que as Equações 2 e 3, uma vez que os coeficientes linear e angular das retas fornecem diretamente os valores de V_max e K_M.

A Equação de Hanes (3) foi investigada por Lineweaver e Burk⁸ que concluiram que a Equação 2 produzia melhores resultados para baixas concentrações de substrato. Entretanto, Wilkinson¹⁴ ao fazer uma análise estatística detalhada das Equações 2 e 3 concluiu que esta última era superior. A comparação estatística feita por Dowd e Riggs¹⁵ das Equações 2, 3 e 4 confirmou a inferioridade da Equação 2 e indicou a Equação 3 como uma melhor opção quando o erro em v₀ era pequeno, mas recomendou a Equação 4 quando o erro em v₀ fosse grande.

Embora uma análise estatística da confiabilidade das diferentes equações lineares empregadas para a determinação dos parâmetros de Michaelis-Menten seja relevante,¹⁶ a abordagem que apresentamos neste trabalho foi a de comparar os valores dos parâmetros obtidos pelas diferentes equações lineares com os obtidos pelo método mais confiável do ajuste não linear por mínimos quadrados.^6,7

A Tabela 1 mostra os parâmetros de Michaelis-Menten calculados a partir dos dados experimentais de velocidade inicial (Conjunto 1) para a atividade de catalase (reação de decomposição do peróxido de hidrogênio) de um complexo de ferro(III). Os valores experimentais de velocidade inicial empregados são mostrados na Tabela 2, juntamente com os valores de v₀ calculados pela Equação 1 usando-se os parâmetros de Michaelis-Menten obtidos tanto pelo método de ajuste não linear por mínimos quadrados quanto pelas Equações 2 a 5. A Figura 1 apresenta uma comparação gráfica dos valores de v₀ obtidos.

Como pode ser visto na Tabela 1, é boa a concordância entre os parâmetros obtidos pelo método de ajuste não linear por mínimos quadrados e as Equações 3 e 5, diferentemente dos valores obtidos pelas Equações de Lineweaver-Burk e Eadie-Hofstee (2 e 4, respectivamente). Os parâmetros obtidos pelo uso das Equações 2 e 4, quando substituídos na Equação 1, produzem valores de v₀ que se desviam significativamente dos valores experimentais, especialmente para as concentrações mais altas de substrato, conforme pode ser visto na Figura 1 e Tabela 2.

A fim de obtermos outros parâmetros de Michaelis-Menten que permitissem comparar o desempenho das Equações 2 a 5 em outras condições, aplicamos este mesmo procedimento para mais três conjuntos de dados: conjunto 2 – hidrólise de um fosfato de diarila por um complexo dinuclear de Cu; conjunto 3 – hidrólise de um fosfato de diarila por um complexo dinuclear de Ni; conjunto 4 – oxidação de um catecol substituído catalisada por um complexo dinuclear de Cu.

Como pode ser visto na Tabela 3, as Equações 3 e 5 produzem, novamente, valores significativamente mais próximos daqueles obtidos por ajuste não linear por mínimos quadrados, do que os valores obtidos pelas Equações 2 e 4. Vale notar que para os conjuntos de dados 2 e 3 os valores de v₀ calculados usando-se a Equação de Lineweaver-Burk (2) são maiores do que os valores calculados por ajuste não linear por mínimos quadrados (Tabelas 4 a 6 e Figuras 2 a 4), o que é uma situação diferente daquela mostrada na Tabela 2 e Figura 1. Por outro lado, os parâmetros de Michaelis-Menten obtidos pela Equação 2 para os conjuntos de dados 1 e 4 (Tabelas 2 e 3, respectivamente) produzem valores de v₀ que são menores dos que os calculados por ajuste não linear por mínimos quadrados (Tabelas 2 e ⁶ e Figuras 1 e ⁴). Além disso, os dados experimentais que escolhemos são especialmente difícieis de serem ajustados, em especial o conjunto de dados 3, para o qual mesmo o ajuste não linear por mínimos quadrados se afasta um pouco dos valores experimentais, conforme se pode ver na Figura 3 (ver também Tabelas 5 e 3).

CONCLUSÕES

Quatro equações lineares foram investigadas quanto a sua capacidade de produzir parâmetros de Michaelis-Menten confiáveis. Como resultado foi observado que duas delas, a Equação de Hanes (3) e a Equação apresentada aqui pela primeira vez (5), podem ser empregadas para determinar os parâmetros V_max e K_M através de um ajuste linear (regressão linear), uma vez que os resultados obtidos são bastante similares aos obtidos por ajuste não linear por mínimos quadrados. As equações escolhidas permitem uma abordagem didática simples e direta uma vez que os parâmetros de Michaelis-Menten são obtidos facilmente a partir dos coeficientes angular e linear das respectivas retas representadas por estas equações. A Equação 5 apresenta uma pequena vantagem de manipulação sobre a Equação 3 pelo fato dos valores de V_max e - K_M se apresentarem diretamente como os coeficientes angular e linear da reta. Por outro lado, embora mais confiável, o uso de um ajuste não linear por mínimos quadrados, requer pacotes de programas mais sofisticados e de manipulação mais complexa, exigindo a estimativa de valores iniciais e a realização de diversas iterações até que os valores dos parâmetros a serem obtidos não variem mais após novas iterações.

A comparação com os parâmetros obtidos por ajuste não linear por mínimos quadrados confirma também que a bem conhecida Equação de Lineweaver-Burk (2) e a Equação de Eadie-Hofstee (4) não produzem bons valores para os parâmetros de Michaelis-Menten, conforme já apontado por outros autores.^2,4,7

AGRADECIMENTOS

Ao apoio financeiro do CNPq, CAPES, INCT-Catalysis e FAPERJ.

6. Em ref. 2, p. 91.

Recebido em 4/11/09; aceito em 26/3/10; publicado na web em 20/7/10

• 1. Briggs, G. E.; Haldane, J. B. S.; Biochem. J. 1925, 19, 338.
• 2. Espenson, J. H.; Chemical Kinetics and Reaction Mechanisms, 2^nd ed.; McGraw-Hill: Nova Iorque, 2002.
• 3. Connors, K. A.; Chemical Kinetics. The Study of Reaction Rates in Solution, John Wiley: Nova Iorque, 1990.
• 4. Cornish-Bowden, A.; Fundamentals of Enzyme Kinetics, Portland Press: Londres, 1995.
• 5. Michaelis, L.; Menten, M. L.; Biochem. Z. 1913, 49, 333.
• 7. Signorella, S.; Rompel, A.; Büldt-Karentzopoulos, K.; Krebs, B.; Pecoraro, V. L.; Tuchagues, J.-P.; Inorg. Chem. 2007, 46, 10864.
• 8. Lineweaver, H.; Burk, D.; J. Am. Chem. Soc. 1934, 56, 658.
• 9. OriginLab Corporation, Origin 6.1; Northamptom, 2000.
• 10. Lotus Development Corporation; Lotus 1-2-3 Release 9 for Windows; 1998.
• 11. Hanes, C. S.; Biochem. J. 1932, 26, 1406.
• 12. Eadie, G. S.; J. Biol. Chem. 1942, 146, 85.
• 13. Hofstee, B. H. J.; J. Biol. Chem., 1952, 199, 357.
• 14. Wilkinson, G. N.; Biochem. J. 1961, 80, 324.
• 15. Dowd, J. E.; Riggs, D. S.; J. Biol. Chem. 1965, 240, 863.
• 16. Cornish-Bowden, A.; Eisenthal, R.; Biochem. J. 1974, 139, 721.

†

In Memoriam

Datas de Publicação

• Publicação nesta coleção
  27 Set 2010
• Data do Fascículo
  2010

Histórico

• Recebido
  04 Nov 2009
• Aceito
  20 Jul 2010
• Revisado
  26 Mar 2010