Insights into the dynamic trajectories of protein filament division revealed by numerical investigation into the mathematical model of pure fragmentation.

The dynamics by which polymeric protein filaments divide in the presence of negligible growth, for example due to the depletion of free monomeric precursors, can be described by the universal mathematical equations of 'pure fragmentation'. The rates of fragmentation reactions reflect the stability of the protein filaments towards breakage, which is of importance in biology and biomedicine for instance in governing the creation of amyloid seeds and the propagation of prions. Here, we devised from mathematical theory inversion formulae to recover the division rates and division kernel information from time-dependent experimental measurements of filament size distribution. The numerical approach to systematically analyze the behaviour of pure fragmentation trajectories was also developed. We illustrate how these formulae can be used, provide some insights on their robustness, and show how they inform the design of experiments to measure fibril fragmentation dynamics. These advances are made possible by our central theoretical result on how the length distribution profile of the solution to the pure fragmentation equation aligns with a steady distribution profile for large times.

The Division of Amyloid Fibrils: Systematic Comparison of Fibril Fragmentation Stability by Linking Theory with Experiments.

The division of amyloid protein fibrils is required for the propagation of the amyloid state and is an important contributor to their stability, pathogenicity, and normal function. Here, we combine kinetic nanoscale imaging experiments with analysis of a mathematical model to resolve and compare the division stability of amyloid fibrils. Our theoretical results show that the division of any type of filament results in self-similar length distributions distinct to each fibril type and the conditions applied. By applying these theoretical results to profile the dynamical stability toward breakage for four different amyloid types, we reveal particular differences in the division properties of disease-related amyloid formed from α-synuclein when compared with non-disease associated model amyloid, the former showing lowered intrinsic stability toward breakage and increased likelihood of shedding smaller particles. Our results enable the comparison of protein filaments' intrinsic dynamic stabilities, which are key to unraveling their toxic and infectious potentials.

A Growth-Fragmentation Approach for Modeling Microtubule Dynamic Instability.

Microtubules (MTs) are protein filaments found in all eukaryotic cells which are crucial for many cellular processes including cell movement, cell differentiation, and cell division. Due to their role in cell division, they are often used as targets for chemotherapy drugs used in cancer treatment. Experimental studies of MT dynamics have played an important role in the development and administration of many novel cancer drugs; however, a complete description of MT dynamics is lacking. Here, we propose a new mathematical model for MT dynamics, that can be used to study the effects of chemotherapy drugs on MT dynamics. Our model consists of a growth-fragmentation equation describing the dynamics of a length distribution of MTs, coupled with two ODEs that describe the dynamics of free GTP- and GDP-tubulin concentrations (the individual dimers that comprise of MTs). Here, we prove the well-posedness of our system and perform a numerical exploration of the influence of certain model parameters on the systems dynamics. In particular, we focus on a qualitative description for how a certain class of destabilizing drugs, the vinca alkaloids, alter MT dynamics. Through variation of certain model parameters which we know are altered by these drugs, we make comparisons between simulation results and what is observed in in vitro studies.

A model of calcium transport along the rat nephron.

We developed a mathematical model of calcium (Ca(2+)) transport along the rat nephron to investigate the factors that promote hypercalciuria. The model is an extension of the flat medullary model of Hervy and Thomas (Am J Physiol Renal Physiol 284: F65-F81, 2003). It explicitly represents all the nephron segments beyond the proximal tubules and distinguishes between superficial and deep nephrons. It solves dynamic conservation equations to determine NaCl, urea, and Ca(2+) concentration profiles in tubules, vasa recta, and the interstitium. Calcium is known to be reabsorbed passively in the thick ascending limbs and actively in the distal convoluted (DCT) and connecting (CNT) tubules. Our model predicts that the passive diffusion of Ca(2+) from the vasa recta and loops of Henle generates a significant axial Ca(2+) concentration gradient in the medullary interstitium. In the base case, the urinary Ca(2+) concentration and fractional excretion are predicted as 2.7 mM and 0.32%, respectively. Urinary Ca(2+) excretion is found to be strongly modulated by water and NaCl reabsorption along the nephron. Our simulations also suggest that Ca(2+) molar flow and concentration profiles differ significantly between superficial and deep nephrons, such that the latter deliver less Ca(2+) to the collecting duct. Finally, our results suggest that the DCT and CNT can act to counteract upstream variations in Ca(2+) transport but not always sufficiently to prevent hypercalciuria.