Mutation-induced conformational remodeling explains gain- and loss-of-function phenotypes in Euglena gracilis phytoene synthase

Resumo Técnico: Remodelação Conformacional Induzida por Mutação na Fitoeno Sintase de Euglena gracilis

Definição do Problema
A fitoeno sintase (CrtB) catalisa a primeira etapa limitante na biossíntese de carotenoides, tornando-se um alvo crítico para a engenharia metabólica. Embora o rastreamento experimental em Euglena gracilis tenha identificado mutações específicas (Y197A, E121H e L349F) que resultam em fenótipos contrastantes de ganho de função (produção aumentada de carotenoides) e perda de função (redução da produção), os mecanismos estruturais que impulsionam essas diferenças permanecem obscuros. Modelos estruturais estáticos e observações experimentais isoladas não conseguem explicar como substituições de um único resíduo alteram a atividade enzimática. O estudo aborda a lacuna na compreensão da base estrutural dinâmica que vincula mutações específicas à produtividade alterada de carotenoides.

Metodologia
Os autores empregaram simulações de Dinâmica Molecular (MD) da forma apo (livre de substrato) da CrtB de E. gracilis para investigar a dinâmica conformacional da enzima selvagem e de três variantes: Y197A, E121H e L349F. A abordagem computacional incluiu:

• Análise Energética: Cálculos MM-GBSA e MM-PBSA para avaliar as perturbações energéticas dependentes de mutação e a estabilidade ($\Delta\Delta G$).
• Caracterização Eletrostática: Análise PROPKA para determinar mudanças nos pontos isoelétricos (pI) e na carga líquida em pH 7,0.
• Métricas Estruturais Globais: Análise de Área de Superfície Acessível ao Solvente (SASA) e Raio de Giro ($R_g$) para avaliar a compacidade global e a amostragem conformacional.
• Dinâmica Local: Análise de Flutuação de Raiz Quadrada Média (RMSF) e $\Delta$RMSF para mapear mudanças na mobilidade específica de resíduos.
• Arquitetura do Sítio Ativo: Exame detalhado da acessibilidade ao solvente, volumes do bolso de ligação e redes de pontes de hidrogênio dentro dos motivos ricos em Asp (essenciais para a coordenação de Mg²⁺) e da aba hidrofóbica (hydrophobic flap).

Principais Resultados
As simulações revelaram que as mutações induzem uma remodelação distinta do panorama de energia conformacional, afetando as redes de interação local e a arquitetica do sítio ativo sem causar mudanças estruturais globais substanciais:

1. Mudanças Energéticas e Eletrostáticas:

   • E121H exibiu a maior redistribuição da dinâmica conformacional, com valores positivos de $\Delta\Delta G$ indicando desestabilização energética em relação à proteína nativa. Também apresentou o pI previsto mais alto e uma leve carga positiva em pH neutro.
   • Y197A exibiu mudanças energéticas neutras a favoráveis.
   • L349F induziu perturbações energéticas modestas, mas mostrou desestabilização pronunciada nas previsões de energia livre de enovelamento via PROPKA.

2. Amostragem Conformacional e Flexibilidade:

   • Variantes de ganho de função (Y197A, E121H): Ambas mantiveram abas hidrofóbicas relativamente acessíveis. A E121H apresentou a distribuição de $R_g$ mais ampla, enquanto Y197A e E121H exibiram distribuições de SASA mais amplas do que a proteína nativa, sugerindo uma amostragem conformacional aumentada mantendo a compacidade geral.
   • Variante de perda de função (L349F): Assemelhou-se estreitamente à proteína nativa nas distribuições globais de SASA e $R_g$, mas exibiu uma constrição localizada dramática no bolso da aba hidrofóbica (redução de volume de aproximadamente 35%).

3. Redes de Interação Local e Remodelação do Sítio Ativo:

   • Aba Hidrofóbica: A mutação L349F causou uma redução significativa no volume do bolso da aba e na acessibilidade ao solvente, consistente com uma conformação mais fechada. Em contraste, as variantes de ganho de função mantiveram maior exposição da aba.
   • Motivos Ricos em Asp: Todos os mutantes mostraram aumento na ocupação de pontes de hidrogênio na região Asp-rica 1. Notavelmente, a L349F exibiu a maior ocupação, apesar da redução no desempenho catalítico, sugerindo que a eficiência catalítica depende de um equilíbrio entre estabilidade e flexibilidade, em vez de uma estabilização local máxima.
   • Rigidez em L349F: A análise estrutural do mutante L349F revelou que o resíduo Tyr127 adotou uma conformação mais rígida, formando uma rede de interação persistente envolvendo interações $\pi$-cátion e $\pi$-ânion com Arg154 e Glu158. Isso contrasta com a proteína nativa e as variantes de ganho de função, onde a Tyr127 permaneceu conformacionalmente flexível.

Significância e Alegações
O artigo afirma fornecer um arcabouço estrutural que vincula a remodelação conformacional induzida por mutação ao redor do centro catalítico às diferenças observadas experimentalmente na produção de carotenoides. Os autores sustentam que os fenótipos contrastantes surgem não meramente de perturbações estruturais locais nos locais de mutação, mas da redistribuição da amostragem conformacional ao redor de motivos funcionalmente importantes, especificamente a aba hidrofóbica e o motivo Asp-rico 1 adjacente.

O estudo destaca que o fenótipo de perda de função em L349F é caracterizado por uma "constrição pronunciada" da aba hidrofóbica e uma "rede de interação rígida", enquanto as variantes de ganho de função mantêm uma "aba hidrofóbica relativamente acessível" e maior flexibilidade. Essas descobertas sugerem que a organização estrutural dinâmica é um determinante crítico da função da CrtB. Os autores concluem que esses insights oferecem uma base para a engenharia racional de fitoeno sintases, embora notem explicitamente que simulações futuras incorporando substratos (GGPP) e íons Mg²⁺ são necessárias para elucidar totalmente como essas mudanças conformacionais intrínsecas influenciam o reconhecimento do substrato e a catálise.