If the 13Cdelta2 chemical shift of neutral ("high pH") histidine is >122 ppm, primarily Ndelta1-H tautomer (2) is indicated; if it is <122 ppm, primarily Nepsilon2-H tautomer (1) is indicated. His resonances from the catalytic triad of active serine proteases, for example, are readily distinguished from those of denatured enzyme. The 13Cdelta2 chemical shifts increased by 6.2 ppm for the catalytic histidines in both alpha-lytic protease and subtilisin BPN' in raising the pH from that of imidazolium cation to that of tautomer 2. This tautomer identification method is easy to implement, requiring only bioincorporation of [U-13C] (or the more readily available [U-13C,15N])-histidine. Standard 1H/13C correlation HMQC or HSQC NMR pulse programs then yield the 13Cdelta2 chemical shifts with the benefit of high 1H sensitivity. Because of large one-bond spin-couplings (1JCH approximately 200 Hz), the method should extend to proteins having large 1H and 13C line widths, including very high molecular weights.

There have been several reports on the varying rates of progression among Alzheimer's Disease (AD) patients; however, there has been no quantitative study of the amount of heterogeneity in AD. Obtaining a reliable quantitative measure of AD progression rates and their variances among the patients for each stage of AD is essential for evaluating results of any clinical study. The Global Deterioration Scale (GDS) and Functional Assessment Staging procedure (FAST) characterize seven stages in the course of AD from normal aging to severe dementia. Each GDS/FAST stage has a published mean duration, but the variance is unknown. We use statistical analysis to reconstruct GDS/FAST stage durations in a cohort of 648 AD patients with an average follow-up time of 4.78 years. Calculations for GDS/FAST stages 4–6 reveal that the standard deviations for stage durations are comparable with their mean values, indicating the presence of large variations in the AD progression among patients. Such amount of heterogeneity in the course of progression of AD is consistent with the existence of several sub-groups of AD patients, which differ by their patterns of decline.

Malignant melanoma is a cancer of the skin arising in the melanocytes. We present a mathematical model of melanoma invasion into healthy tissue with an immune response. We use this model as a framework with which to investigate primary tumor invasion and treatment by surgical excision. We observe that the presence of immune cells can destroy tumors, hold them to minimal expansion, or, through the production of angiogenic factors, induce tumorigenic expansion. We also find that the tumor–immune system dynamic is critically important in determining the likelihood and extent of tumor regrowth following resection. We find that small metastatic lesions distal to the primary tumor mass can be held to a minimal size via the immune interaction with the larger primary tumor. Numerical experiments further suggest that metastatic disease is optimally suppressed by immune activation when the primary tumor is moderately, rather than minimally, metastatic. Furthermore, satellite lesions can become aggressively tumorigenic upon removal of the primary tumor and its associated immune tissue. This can lead to recurrence where total cancer mass increases more quickly than in primary tumor invasion, representing a clinically more dangerous disease state. These results are in line with clinical case studies involving resection of a primary melanoma followed by recurrence in local metastases.

Abstract. Quantitative systems pharmacology (QSP) modeling has become increasingly important in pharmaceutical research and development, and is a powerful tool to gain mechanistic insights into the complex dynamics of biological systems in response to drug treatment. However, even once a suitable mathematical framework to describe the pathophysiology and mechanisms of interest is established, final model calibration and the exploration of variability can be challenging and time consuming. QSP models are often formulated as multi-scale, multi-compartment nonlinear systems of ordinary differential equations. Commonly accepted modeling strategies, workflows, and tools have promise to greatly improve the efficiency of QSP methods and improve productivity. In this paper, we present the QSP Toolbox, a set of functions, structure array conventions, and class definitions that computationally implement critical elements of QSP workflows including data integration, model calibration, and variability exploration. We present the application of the toolbox to an ordinary differential equations-based model for antibody drug conjugates. As opposed to a single stepwise reference model calibration, the toolbox also facilitates simultaneous parameter optimization and variation across multiple in vitro, in vivo, and clinical assays to more comprehensively generate alternate mechanistic hypotheses that are in quantitative agreement with available data. The toolbox also includes scripts for developing and applying virtual populations to mechanistic exploration of biomarkers and efficacy. We anticipate that the QSP Toolbox will be a useful resource that will facilitate implementation, evaluation, and sharing of new methodologies in a common framework that will greatly benefit the community.

BackgroundWe study the selection dynamics in a heterogeneous spatial colony of cells. We use two spatial generalizations of the Moran process, which include cell divisions, death and migration. In the first model, migration is included explicitly as movement to a proximal location. In the second, migration is implicit, through the varied ability of cell types to place their offspring a distance away, in response to another cell's death.ResultsIn both models, we find that migration has a direct positive impact on the ability of a single mutant cell to invade a pre-existing colony. Thus, a decrease in the growth potential can be compensated by an increase in cell migration. We further find that the neutral ridges (the set of all types with the invasion probability equal to that of the host cells) remain invariant under the increase of system size (for large system sizes), thus making the invasion probability a universal characteristic of the cells selection status. We find that repeated instances of large scale cell-death, such as might arise during therapeutic intervention or host response, strongly select for the migratory phenotype.ConclusionsThese models can help explain the many examples in the biological literature, where genes involved in cell's migratory and invasive machinery are also associated with increased cellular fitness, even though there is no known direct effect of these genes on the cellular reproduction. The models can also help to explain how chemotherapy may provide a selection mechanism for highly invasive phenotypes.ReviewersThis article was reviewed by Marek Kimmel and Glenn Webb.