Search

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. Author summary: Amyloid fibrils are fibrillar protein structures involved in many neurodegenerative illnesses, such as Parkinson's disease or Alzheimer's disease. To propagate in disease, these misfolded protein aggregates must grow and divide to proliferate. Therefore, the intrinsic characteristics of their division, including the division rate and the pattern of division in terms of whether the fibrils are likely to break in the middle or at the edges, impact the disease aetiology. Here, we discovered mathematical formulae that can be used to directly extract the fibril division characteristics from recent experiments data obtained from time-dependent fibril length distribution measurements. We explain how these formulae can be used, and prove the robustness of the division rate formula where small errors in the measurement leads to small errors in the division rate. We also demonstrate that the mathematical formula is not robust enough to precisely decipher the pattern of division in the data, and suggest instead new future experimental design with short time measurements in experiments starting with fibril suspensions where all fibrils have similar size, which would be suitable to provide an improved estimates. [ABSTRACT FROM AUTHOR]

We consider a scalar conservation law with a flux containing spatial heterogeneities of bounded variation, where the number of discontinuities may be infinite. We prove existence of an adapted entropy solution in the sense of Audusse and Perthame in the BV framework, under assumptions that guarantee uniqueness of the solution, plus one additional technical assumption. In a counterexample, we prove that existence in the BV framework fails without this additional technical assumption. [ABSTRACT FROM AUTHOR]

We developed a mathematical model of calcium (Ca2+) 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 Ca2+ 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 Ca2+ from the vasa recta and loops of Henle generates a significant axial Ca2+ concentration gradient in the medullary interstitium. In the base case, the urinary Ca2+ concentration and fractional excretion are predicted as 2.7 mM and 0.32%, respectively. Urinary Ca2+ excretion is found to be strongly modulated by water and NaCl reabsorption along the nephron. Our simulations also suggest that Ca2+ molar flow and concentration profiles differ significantly between superficial and deep nephrons, such that the latter deliver less Ca2+ to the collecting duct. Finally, our results suggest that the DCT and CNT can act to counteract upstream variations in Ca2+ transport but not always sufficiently to prevent hypercalciuria. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

Abstract We consider the fragmentation equation ∂ f ∂ t (t , x) = − B (x) f (t , x) + ∫ y = x y = ∞ k (y , x) B (y) f (t , y) d y , and address the question of estimating the fragmentation parameters – i.e. the division rate B (x) and the fragmentation kernel k (y , x) – from measurements of the size distribution f (t , ⋅) at various times. This is a natural question for any application where the sizes of the particles are measured experimentally whereas the fragmentation rates are unknown, see for instance Xue and Radford (2013) [26] for amyloid fibril breakage. Under the assumption of a polynomial division rate B (x) = α x γ and a self-similar fragmentation kernel k (y , x) = 1 y k 0 (x y) , we use the asymptotic behavior proved in Escobedo et al. (2004) [11] to obtain uniqueness of the triplet (α , γ , k 0) and a representation formula for k 0. To invert this formula, one of the delicate points is to prove that the Mellin transform of the asymptotic profile never vanishes, what we do through the use of the Cauchy integral. [ABSTRACT FROM AUTHOR]