Virtual Cell Virtual Cell WebsiteSchool of Health Sciences WebsitePurdue University Website
About the Virtual CellDownload VCExamplesTheoryUser ManualHelpExecute VC Simulation
TLK Model Input File #4
Sample Input Files

Tip: Comments are indicated in blue.  For additional details, click on the keywords highlighted in green. Reasonable default values are specified for all model inputs.

/*
TLK Model Calibration Example: Use this sample file to begin analyzing delayed-plating cell survival data (no repopulation effects; GF=0.0).  For small datasets (e.g., only low or high dose rate survival data are available), it may be prudent to perform the data analysis using the LPL or RMR model instead of the TLK.  Nevertheless, the TLK model can still provide some useful insights into the radiobiology -- even for small datasets. For example, the TLK might be used to investigate the impact on parameter estimates of biphasic DSB rejoining kinetics.  Recommendations for TLK model fitting are provided below.  For large datasets (multiple dose rates and/or split-dose survival data), the TLK can provide better fits to the data than the LPL or RMR model.

Small datasets (low or high dose rate data only): Use the default parameters indicated below and treat the ETA and PHI parameters as adjustable (2 or 3 parameter fit).

Large datasets (multiple dose rates): Use the default parameters indicated below.  An adequate fit to most datasets can be obtained by treating 3 to 5 of the TLK model inputs as adjustable.  Adjust ETA, PHI (both kinds of DSB), and RHT for simple DSBs to perform a 3 parameter fit.  For a 4 parameter fit, adjust ETA, PHI, and the RHT for both simple and complex DSBs.  And for a 5 parameter fit, adjust ETA, PHI, the RHT for simple and complex DSBs, and the initial yield of complex DSBs (DSB).  Start with a 3 or 4 parameter fit and increase the number of adjustable parameters until you obtain a good fit to the data.
*/

   MODEL: DRM=TLK CKM=QECK         !Model Specification
   CELL: DNA=5.667D+09 NC=46
   IRAD: DSB=25 FCB=0.2                   !Radiation damage to DNA
   CLRP: RHT=HT4R ETA=2.75E-04 PHI=0.01 A0=PMR  GAM=0.25
   CCKM: TCC=10 TPOT=20 GF=0.0 N0=1 KAP=1.0D+38 VOL=1
   RBM: KX=6.0E+06 FMAX=1   !Analysis of PFGE data

!  TLK model inputs (simple and complex double strand breaks)
   HT4R = 0.25 7 /             !Half-times for DSB repair
   PMR = 0.75 0.25 /         !Probability DSB is rejoined correctly

!  Parameters that impact on solution accuracy
   SOLVER: DBLV=0 ETOL=1.0E-09 MXSS=999999 HMAX=10

!  Parameters that control the time steps and stopping criteria
   SIMCON: FSDX=1 TSAX=12 TCUT=96 FRDL=1.0D+00 RCUT=1.0D-09

   BCONF: MITR=10000 MXFN=100000 SGT=1.0D-04 SST=1.0D-09 FCT=1.0D-09
   MCAT: NPAR=4  !npar=2 or 3 for small datasets, npar=3-5 for large datasets
   AMP1: CID=CLRP PNAM=ETA   PLB=0.0 PUB=1.0D-02 PIG=2.0E-04
   AMP2: CID=CLRP PNAM=PHI    PLB=0.0 PUB=0.25 PIG=0.01
   AMP3: CID=CLRP PNAM=RHT   PLB=0.01 PUB=12 PIG=0.15 LI=1    !Simple DSB
   AMP4: CID=CLRP PNAM=RHT   PLB=0.01 PUB=12 PIG=8 LI=2        !Complex DSB
   AMP5: CID=IRAD  PNAM=DSB   PLB=20 PUB=60.0 PIG=25              !Total DSB


/*
  Figure of Merit (FOM) Weight
  1: cell survival data
  2: DSB rejoining data (FAR)
  3: DSB rejoining data (Avg. DSB/cell)
  4: Transformed cells per survivor
  5: Transformed cells per irradiated cell
 */

  FOM: UW=LFWT
  LFWT = 1 0 0 0 0 /


School of Health Sciences
Purdue University
Disclaimer

Last updated: 10 June, 2011
Webmaster