Homeostatic Noise Buffering in Biomolecular Condensates Hinges on Phase Multiplicity Modulated by Interfacial and Droplet Size Effects

Resumo Técnico: O Amortecimento de Ruído Homeostático em Condensados Biomoleculares Depende da Multiplicidade de Fases Modulada por Efeitos de Interface e de Tamanho de Gotícula

Definição do Problema
Condensados biomoleculares, formados via separação de fases líquido-líquido (LLPS), servem como compartimentos sem membrana que regulam atividades celulares ao enriquecer ou excluir moléculas específicas. Uma função física crítica desses condensados é o amortecimento passivo de flutuações estocásticas de concentração (ruído) na dinâmica celular. Embora a regra de fases de Gibbs sugira que sistemas multifásicos (ex: LLPS ternária) poderiam, teoricamente, oferecer um amortecimento de ruído superior em comparação a sistemas binários ao manter composições de fase constantes sobre um regime de concentração 2D (em oposição a uma linha de ligação 1D), esse potencial não foi substancialmente abordado. Além disso, modelos teóricos existentes frequentemente assumem o limite termodinâmico (volume infinito, $V \to \infty$), negligenciando efeitos de tamanho finito. Isso cria uma desconexão entre teorias analíticas e simulações de dinâmica molecular (MD), que tipicamente modelam "nanocondensados" (dimensões lineares $\sim 15\text{--}50$ nm) onde flutuações interfaciais e térmicas são significativas. O problema central abordado é como o tamanho do condensado e a multiplicidade de fases (binária vs. ternária) influenciam conjuntamente a eficácia do amortecimento de ruído, e se simulações de MD de pequena escala podem modelar de forma confiável o comportamento de condensados maiores e biologicamente realistas.

Metodologia
O estudo emprega uma abordagem multifacetada combinando teoria analítica, simulações de grão grosso e modelagem fenomenológica:

1. Teoria da Aproximação de Fase Aleatória (RPA): Utilizada para modelar os comportamentos de fase em volume infinito ($V \to \infty$) de pares de sequências poliamfolíticas globalmente neutras (modelo de proteínas intrinsecamente desordenadas, IDPs). A teoria considera a conectividade das cadeias, interações Coulomb e volume excluído. Diagramas de fase foram construídos para identificar regiões de LLPS binária e ternária como função da temperatura e do padrão de carga da sequência.
2. Dinâmica Molecular (MD) de Grão Grosso: Simulações foram realizadas usando um potencial "1/3 Lennard-Jones" para modelar sistemas de tamanho finito (tipicamente 500 cadeias, dimensões lineares $\sim 15\text{--}50$ nm). Essas simulações testaram se os comportamentos de fase previstos pela RPA de volume infinito (especificamente a demixagem ternária) são recapitulados em pequenas gotículas e avaliaram o grau de dispersão de concentração.
3. Teoria de Campo de Fase: Um modelo fenomenológico foi desenvolvido para fazer a ponte entre MD e RPA. Esta abordagem aumenta a densidade de energia livre da RPA com um termo de energia interfacial explícito (energia de gradiente) para considerar efeitos de tamanho finito e tensão superficial. Ela permite a solução numérica de distribuições de densidade em equilíbrio em volumes finitos.
4. Monte Carlo Semelhante ao Ensemble de Gibbs: Uma formulação de MC baseada em voxels foi usada para isolar e quantificar a contribuição das flutuações térmicas para a dispersão de concentração em volumes finitos, independentemente de efeitos interfaciais.

Principais Contribuições e Resultados

• Previsão de LLPS Ternária: A teoria RPA prevê que pares específicos de poliamfolitos com padrões de carga de sequência díspares sofrem LLPS ternária (três fases coexistentes) em certas temperaturas, além da LLPS binária. Este regime ternário é caracterizado por uma região triangular no plano de concentração 2D onde qualquer composição inicial dentro do triângulo se separa em três fases fixas.
• Desvios de Tamanho Finito em MD: Simulações de MD de pequenas gotículas ($\sim 15\text{--}50$ nm) revelam que, embora as tendências dependentes da sequência (demixagem vs. mistura) se alinhem com a RPA, a separação de fase estrita prevista pela teoria de volume infinito não é observada. Em vez disso, pequenas gotículas exibem regiões interfaciais substanciais e sobreposições de populações, levando a uma distribuição de saída de concentrações "difusa" mesmo para um único input. Consequentemente, pequenas gotículas falham em alcançar o amortecimento de ruído ideal da LLPS ternária, performando de forma semelhante a sistemas binários em termos de dispersão de concentração.
• Papel dos Efeitos Interfaciais e Térmicos: O estudo demonstra que o desvio do comportamento ideal da RPA em MD é impulsionado por dois efeitos de tamanho finito:
  1. Efeitos Interfaciais: Modelados com sucesso pela teoria de campo de fase com termos de tensão superficial explícitos. À medida que o tamanho do sistema ($L$) aumenta, a contribuição interfacial diminui e o sistema recupera as composições de fase nítidas e fixas previstas pela RPA.
  2. Flutuações Térmicas: A análise de MC confirma que flutuações térmicas contribuem para a dispersão de concentração em volumes pequenos. No entanto, para gotículas com dimensões lineares $\gtrsim 40$ nm, os efeitos térmicos tornam-se negligenciáveis em comparação aos efeitos interfaciais.
• Amortecimento de Ruído Dependente do Tamanho: A eficácia do amortecimento de ruído é fortemente dependente do tamanho.
  • Para pequenas gotículas ($L \lesssim 200b$, onde $b \approx 0,38$ nm, correspondendo a $L \lesssim 80$ nm), sistemas de LLPS ternária e binária exibem capacidades de amortecimento de ruído comparáveis e fracas devido à significativa dispersão de concentração.
  • Para gotas maiores ($L \gtrsim 200b$), sistemas de LLPS ternária demonstram um amortecimento de ruído significativamente superior aos sistemas binários. Neste regime, as composições das três fases coexistentes permanecem invariantes sobre um intervalo de concentração 2D, efetivamente suprimindo o ruído de entrada.
• Reconciliação de Escalas: A análise fornece uma base para extrapolar resultados de simulações de MD de pequena escala para condensados biológicos maiores. Ao incorporar correções de tamanho finito (tensão interfacial e flutuações térmicas) em teorias analíticas, o estudo sugere que o comportamento "ruidoso" observado em simulações típicas de MD é um artefato de tamanho finito que diminui à medida que o tamanho do sistema se aproxima da escala de nanocondensados biológicos ($\sim 500$ nm) e organelas de escala micrométrica.

Significância
O artigo estabelece que a vantagem teórica da LLPS ternária para o amortecimento homeostático de ruído é contingente ao fato de o tamanho do condensado ser suficientemente grande para minimizar os efeitos de tamanho finito. Destaca que, embora os condensados biomoleculares naturais abranjam uma ampla faixa de tamanhos (desde "nanocondensados" até organelas micrométricas), o mecanismo físico de amortecimento de ruído transita de ser dominado por flutuações interfaciais e térmicas em sistemas pequenos para ser governado pelas regras de fase termodinâmica de volume em sistemas maiores.

O estudo ressalta a utilidade das teorias analíticas (RPA e campo de fase) como ferramentas para interpretar e corrigir as limitações de simulações de MD de pequena escala. Sugere que, para sistemas poliamfolíticos, condensados com dimensões lineares superiores a $\sim 80$ nm já podem exibir propriedades de amortecimento de ruído aproximando-se do limite ideal de volume infinito, tornando a LLPS ternária um mecanismo potencialmente robusto para a homeostase em muitos contextos biológicos. Além disso, o trabalho oferece um arcabouço para entender como os padrões de carga específicos de sequência ditam se um sistema sofrerá separação de fase binária ou ternária, oferecendo um caminho para o design racional de condensados multifásicos.