Molecular dynamics simulations of a multicellular model with cell-cell interactions and Hippo signaling pathway

PLoS Comput Biol. 2024 Nov 11;20(11):e1012536. doi: 10.1371/journal.pcbi.1012536. eCollection 2024 Nov.

Abstract

The transcriptional coactivator Yes-associated protein (YAP)/transcriptional co-activator with PDZ binding motif (TAZ) induces cell proliferation through nuclear localization at low cell density. Conversely, at extremely high cell density, the Hippo pathway, which regulates YAP/TAZ, is activated. This activation leads to the translocation of YAP/TAZ into the cytoplasm, resulting in cell cycle arrest. Various cancer cells have several times more YAP/TAZ than normal cells. However, it is not entirely clear whether this several-fold increase in YAP/TAZ alone is sufficient to overcome proliferation inhibition (contact inhibition) under high-density conditions, thereby allowing continuous proliferation. In this study, we construct a three-dimensional (3D) mathematical model of cell proliferation incorporating the Hippo-YAP/TAZ pathway. Herein, a significant innovation in our approach is the introduction of a novel modeling component that inputs cell density, which reflects cell dynamics, into the Hippo pathway and enables the simulation of cell proliferation as the output response. We assume such 3D model with cell-cell interactions by solving reaction and molecular dynamics (MD) equations by applying adhesion and repulsive forces that act between cells and frictional forces acting on each cell. We assume Lennard-Jones (12-6) potential with a softcore character so that each cell secures its exclusive domain. We set cell cycles composed of mitotic and cell growth phases in which cells divide and grow under the influence of cell kinetics. We perform mathematical simulations at various YAP/TAZ levels to investigate the extent of YAP/TAZ increase required for sustained proliferation at high density. The results show that a twofold increase in YAP/TAZ levels of cancer cells was sufficient to evade cell cycle arrest compared to normal cells, enabling cells to continue proliferating even under high-density conditions. Finally, this mathematical model, which incorporates cell-cell interactions and the Hippo-YAP/TAZ pathway, may be applicable for evaluating cancer malignancy based on YAP/TAZ levels, developing drugs to suppress the abnormal proliferation of cancer cells, and determining appropriate drug dosages. The source codes are freely available.

MeSH terms

  • Adaptor Proteins, Signal Transducing / metabolism
  • Cell Communication* / physiology
  • Cell Proliferation* / physiology
  • Hippo Signaling Pathway*
  • Humans
  • Models, Biological
  • Molecular Dynamics Simulation*
  • Protein Serine-Threonine Kinases* / metabolism
  • Signal Transduction* / physiology
  • Transcription Factors* / metabolism
  • YAP-Signaling Proteins / metabolism

Substances

  • Protein Serine-Threonine Kinases
  • Transcription Factors
  • YAP-Signaling Proteins
  • Adaptor Proteins, Signal Transducing

Grants and funding

The work was supported by the following grants: M.T. received funding from the Japan Science and Technology Agency (JST), CREST (Grant No. JPMJCR2022). The funder’s website: https://projectdb.jst.go.jp/grant/JST-PROJECT-20348631/. H.M. received funding from the Japan Society for the Promotion of Science (JSPS) (Grant No. 24K18108). The funder’s website: https://kaken.nii.ac.jp/grant/KAKENHI-PROJECT-24K18108. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.