Objective: A molecular dynamics approach to understanding fundamental mechanisms of combined thermal and osmotic stress induced by thermochemical ablation (TCA) is presented.
Methods: Structural models of fibronectin and fibronectin bound to its integrin receptor provide idealized models for the effects of thermal and osmotic stress in the extracellular matrix. Fibronectin binding to integrin is known to facilitate cell survival. The extracellular environment produced by TCA at the lesion boundary was modelled at 37 °C and 43 °C with added sodium chloride (NaCl) concentrations (0, 40, 80, 160, and 320 mM). Atomistic simulations of solvated proteins were performed using the GROMOS96 force field and TIP3P water model. Computational results were compared with the results of viability studies of human hepatocellular carcinoma (HCC) cell lines HepG2 and Hep3B under matching thermal and osmotic experimental conditions.
Results: Cell viability was inversely correlated with hyperthermal and hyperosmotic stresses. Added NaCl concentrations were correlated with a root mean square fluctuation increase of the fibronectin arginylglycylaspartic acid (RGD) binding domain. Computed interaction coefficients demonstrate preferential hydration of the protein model and are correlated with salt-induced strengthening of hydrophobic interactions. Under the combined hyperthermal and hyperosmotic stress conditions (43 °C and 320 mM added NaCl), the free energy change required for fibronectin binding to integrin was less favorable than that for binding under control conditions (37 °C and 0 mM added NaCl).
Conclusion: Results quantify multiple measures of structural changes as a function of temperature increase and addition of NaCl to the solution. Correlations between cell viability and stability measures suggest that protein aggregates, non-functional proteins, and less favorable cell attachment conditions have a role in TCA-induced cell stress.
Keywords: Hyperthermia; molecular dynamics; thermochemical ablation.