Sturrock Marc, Hellander Andreas, Aldakheel Sahar, Petzold Linda, Chaplain Mark A J
University of Dundee, Dundee, DD1 4HN, UK.
Bull Math Biol. 2014 Apr;76(4):766-98. doi: 10.1007/s11538-013-9842-5. Epub 2013 May 18.
Hes1 is a member of the family of basic helix-loop-helix transcription factors and the Hes1 gene regulatory network (GRN) may be described as the canonical example of transcriptional control in eukaryotic cells, since it involves only the Hes1 protein and its own mRNA. Recently, the Hes1 protein has been established as an excellent target for an anti-cancer drug treatment, with the design of a small molecule Hes1 dimerisation inhibitor representing a promising if challenging approach to therapy. In this paper, we extend a previous spatial stochastic model of the Hes1 GRN to include nuclear transport and dimerisation of Hes1 monomers. Initially, we assume that dimerisation occurs only in the cytoplasm, with only dimers being imported into the nucleus. Stochastic simulations of this novel model using the URDME software show that oscillatory dynamics in agreement with experimental studies are retained. Furthermore, we find that our model is robust to changes in the nuclear transport and dimerisation parameters. However, since the precise dynamics of the nuclear import of Hes1 and the localisation of the dimerisation reaction are not known, we consider a second modelling scenario in which we allow for both Hes1 monomers and dimers to be imported into the nucleus, and we allow dimerisation of Hes1 to occur everywhere in the cell. Once again, computational solutions of this second model produce oscillatory dynamics in agreement with experimental studies. We also explore sensitivity of the numerical solutions to nuclear transport and dimerisation parameters. Finally, we compare and contrast the two different modelling scenarios using numerical experiments that simulate dimer disruption, and suggest a biological experiment that could distinguish which model more faithfully captures the Hes1 GRN.
Hes1是碱性螺旋-环-螺旋转录因子家族的成员,Hes1基因调控网络(GRN)可被描述为真核细胞中转录控制的典型例子,因为它仅涉及Hes1蛋白及其自身的mRNA。最近,Hes1蛋白已被确立为抗癌药物治疗的一个极佳靶点,设计一种小分子Hes1二聚化抑制剂代表了一种虽具挑战性但很有前景的治疗方法。在本文中,我们扩展了先前的Hes1 GRN空间随机模型,以纳入Hes1单体的核转运和二聚化。最初,我们假设二聚化仅发生在细胞质中,只有二聚体被导入细胞核。使用URDME软件对这个新模型进行的随机模拟表明,与实验研究一致的振荡动力学得以保留。此外,我们发现我们的模型对核转运和二聚化参数的变化具有鲁棒性。然而,由于Hes1核输入的精确动力学以及二聚化反应的定位尚不清楚,我们考虑第二种建模方案,即允许Hes1单体和二聚体都被导入细胞核,并允许Hes1在细胞内各处发生二聚化。同样,这个第二种模型的计算解产生了与实验研究一致的振荡动力学。我们还探讨了数值解对核转运和二聚化参数的敏感性。最后,我们使用模拟二聚体破坏的数值实验比较和对比了这两种不同的建模方案,并提出了一个生物学实验,该实验可以区分哪种模型更忠实地捕捉了Hes1 GRN。