Measuring the signs of the methyl <Superscript>1</Superscript>H chemical shift differences between major and ‘invisible’ minor protein conformational states using methyl <Superscript>1</Superscript>H multi-quantum spectroscopy

Carr–Purcell–Meiboom–Gill (CPMG) type relaxation dispersion experiments are now routinely used to characterise protein conformational dynamics that occurs on the μs to millisecond (ms) timescale between a visible major state and ‘invisible’ minor states. The exchange rate(s) (\( k_{{{\text{ex}}}} \)), population(s) of the minor state(s) and the absolute value of the chemical shift difference \(|{\Delta \varpi }|\) (ppm) between different exchanging states can be extracted from the CPMG data. However the sign of \({\Delta \varpi }\) that is required to reconstruct the spectrum of the ‘invisible’ minor state(s) cannot be obtained from CPMG data alone. Building upon the recently developed triple quantum (TQ) methyl \( ^{1} {\text{H}} \) CPMG experiment (Yuwen in Angew Chem 55:11490–11494, 2016) we have developed pulse sequences that use carbon detection to generate and evolve single quantum (SQ), double quantum (DQ) and TQ coherences from methyl protons in the indirect dimension to measure the chemical exchange-induced shifts of the SQ, DQ and TQ coherences from which the sign of \({\Delta \varpi }\) is readily obtained for two state exchange. Further a combined analysis of the CPMG data and the difference in exchange induced shifts between the SQ and DQ resonances and between the SQ and TQ resonances improves the estimates of exchange parameters like the population of the minor state. We demonstrate the use of these experiments on two proteins undergoing exchange: (1) the ~ 18 kDa cavity mutant of T4 Lysozyme (\( k_{{{\text{ex}}}} \sim\,3500{\text{ s}}^{{ - 1}} \)) and (2) the \(\sim\,4.7\) kDa Peripheral Sub-unit Binding Domain (PSBD) from the acetyl transferase of Bacillus stearothermophilus (\(k_{ex} \sim\,13,000\hbox { s}^{-1}\)).     

Improved identifiability of myocardial material parameters by an energy-based cost function

Myocardial stiffness is a valuable clinical biomarker for the monitoring and stratification of heart failure (HF). Cardiac finite element models provide a biomechanical framework for the assessment of stiffness through the determination of the myocardial constitutive model parameters. The reported parameter intercorrelations in popular constitutive relations, however, obstruct the unique estimation of material parameters and limit the reliable translation of this stiffness metric to clinical practice. Focusing on the role of the cost function (CF) in parameter identifiability, we investigate the performance of a set of geometric indices (based on displacements, strains, cavity volume, wall thickness and apicobasal dimension of the ventricle) and a novel CF derived from energy conservation. Our results, with a commonly used transversely isotropic material model (proposed by Guccione et al.), demonstrate that a single geometry-based CF is unable to uniquely constrain the parameter space. The energy-based CF, conversely, isolates one of the parameters and in conjunction with one of the geometric metrics provides a unique estimation of the parameter set. This gives rise to a new methodology for estimating myocardial material parameters based on the combination of deformation and energetics analysis. The accuracy of the pipeline is demonstrated in silico, and its robustness in vivo, in a total of 8 clinical data sets (7 HF and one control). The mean identified parameters of the Guccione material law were \(C_1=3000\pm 1700\,\hbox {Pa}\) and \(\alpha =45\pm 25\) (\(b_f=25\pm 14\), \(b_{ft}=11\pm 6\), \(b_{t}=9\pm 5\)) for the HF cases and \(C_1=1700\,\hbox {Pa}\) and \(\alpha =15\) (\(b_f=8\), \(b_{ft}=4\), \(b_{t}=3\)) for the healthy case.     

A Mathematical Model Supports a Key Role for Ae4 (Slc4a9) in Salivary Gland Secretion

We develop a mathematical model of a salivary gland acinar cell with the objective of investigating the role of two \(\mathrm{Cl}^-/\mathrm{HCO}_3^-\) exchangers from the solute carrier family 4 (Slc4), Ae2 (Slc4a2) and Ae4 (Slc4a9), in fluid secretion. Water transport in this type of cell is predominantly driven by \(\mathrm{Cl}^-\) movement. Here, a basolateral \(\mathrm{Na}^+/ \mathrm{K}^+\) adenosine triphosphatase pump (NaK-ATPase) and a \(\mathrm{Na}^+\)–\(\mathrm{K}^+\)–\(2 \mathrm{Cl}^-\) cotransporter (Nkcc1) are primarily responsible for concentrating the intracellular space with \(\mathrm{Cl}^-\) well above its equilibrium potential. Gustatory and olfactory stimuli induce the release of \(\mathrm{Ca}^{2+}\) ions from the internal stores of acinar cells, which triggers saliva secretion. \(\mathrm{Ca}^{2+}\)-dependent \(\mathrm{Cl}^-\) and \(\mathrm{K}^+\) channels promote ion secretion into the luminal space, thus creating an osmotic gradient that promotes water movement in the secretory direction. The current model for saliva secretion proposes that \(\mathrm{Cl}^-/ \mathrm{HCO}_3^-\) anion exchangers (Ae), coupled with a basolateral \(\mathrm{Na}^+/\hbox {proton}\) (\(\hbox {H}^+\)) (Nhe1) antiporter, regulate intracellular pH and act as a secondary \(\mathrm{Cl}^-\) uptake mechanism (Nauntofte in Am J Physiol Gastrointest Liver Physiol 263(6):G823–G837, 1992; Melvin et al. in Annu Rev Physiol 67:445–469, 2005.  https://doi.org/10.1146/annurev.physiol.67.041703.084745). Recent studies demonstrated that Ae4 deficient mice exhibit an approximate \(30\%\) decrease in gland salivation (Peña-Münzenmayer et al. in J Biol Chem 290(17):10677–10688, 2015). Surprisingly, the same study revealed that absence of Ae2 does not impair salivation, as previously suggested. These results seem to indicate that the Ae4 may be responsible for the majority of the secondary \(\mathrm{Cl}^-\) uptake and thus a key mechanism for saliva secretion. Here, by using ‘in-silico’ Ae2 and Ae4 knockout simulations, we produced mathematical support for such controversial findings. Our results suggest that the exchanger’s cotransport of monovalent cations is likely to be important in establishing the osmotic gradient necessary for optimal transepithelial fluid movement.     

Plausibility and parameter sensitivity of micro-finite element-based joint load prediction at the proximal femur

A micro-finite element-based method to estimate the bone loading history based on bone architecture was recently presented in the literature. However, a thorough investigation of the parameter sensitivity and plausibility of this method to predict joint loads is still missing. The goals of this study were (1) to analyse the parameter sensitivity of the joint load predictions at one proximal femur and (2) to assess the plausibility of the results by comparing load predictions of ten proximal femora to in vivo hip joint forces measured with instrumented prostheses (available from www.orthoload.com). Joint loads were predicted by optimally scaling the magnitude of four unit loads (inclined \(-20^{\circ }\) to \(100^{\circ }\) with respect to the vertical axis) applied to micro-finite element models created from high-resolution computed tomography scans (\(30.3~\upmu \)m voxel size). Parameter sensitivity analysis was performed by varying a total of nine parameters and showed that predictions of the peak load directions (range 10\(^{\circ }\)–\(30^{\circ }\)) are more robust than the predicted peak load magnitudes (range 2344.8–4689.5 N). Comparing the results of all ten femora with the in vivo loading data of ten subjects showed that peak loads are plausible both in terms of the load direction (in vivo: \(18.2\pm 2.0^{\circ }\), predicted: \(20.0^{\circ }\)) and magnitude (in vivo: \(2707.6\pm 443.3~\hbox {N}\), predicted: \(3372.2\pm 597.9~\hbox {N}\)). Overall, this study suggests that micro-finite element-based joint load predictions are both plausible and robust in terms of the predicted peak load direction, but predicted load magnitudes should be interpreted with caution.     

Prediction of heterosis in rice based on divergence of morphological and molecular markers

Identifying the best performing hybrid without a field test was essential to save resources and time. In this study, the genetic divergence was estimated using morphological and expressed sequence tag (EST)-derived simple sequence repeats (SSR) markers. Cluster analysis showed that APMS6A and RPHR 1005 belong to groups I and II, respectively, and the hybrid combination recorded the highest mean grain yield of 32.25 g among generated 40 \(\hbox {F}_{1}\hbox {s}\) with standard heterosis of 8.4% over hybrid check, KRH2. The coefficient of marker polymorphism (CMP) value was calculated based on EST-SSR markers; it ranged from 0.40 to 0.80, and a higher CMP value of 0.80 was obtained for the parental combination APMS6A \(\times \) RPHR1005. We predicted heterosis for 40 \(\hbox {F}_{1}\hbox {s}\) based on correlation between CMP and standard heterosis in different traits with standard varietal and hybrid checks indicating positive correlation and significant value for grain yield per plant (\(r=0.58\)**), productivity per day (\(r=0.54\)**), productive tillers (\(r=0.34\)*) and panicle weight (\(r=0.42\)**). This study revealed that the relationship of molecular marker heterozygosity, along with the combining ability, high mean value of different traits, grouping of parental lines based on morphological and molecular characterization is helpful to identify heterotic patterns in rice.     

DFT study of nano zinc/copper voltaic cells

To facilitate the development of new materials for use in batteries, it is necessary to develop ab initio full-electron computational techniques for modeling potential new battery materials. Here, we tested density functional theory procedures that are accurate enough to obtain the energetics of a zinc/copper voltaic cell. We found the magnitude of the zero-point energy correction to be 0.01–0.2 kcal/mol per atom or molecule and the magnitude of the dispersion correction to be 0.1–0.6 kcal/mol per atom or molecule for Zn _n , (H₂O) _n , \( \mathrm{Zn}{\left({\mathrm{H}}_2\mathrm{O}\right)}_n^{2+} \), \( \mathrm{Cu}{\left({\mathrm{H}}_2\mathrm{O}\right)}_n^{2+} \), and Cu _n . Counterpoise correction significantly affected the values of ?\( {E}_n^{\mathrm{abs}} \), ?\( {E}_n^{\mathrm{coh}} \), and ?E^solv by 1.0–3.1 kcal/mol per atom or molecule at the B3PW91/6-31G(d) level of theory, but by only 0.04–0.4 kcal/mol per atom or molecule at the B3PW91/cc-pVTZ level of theory. The application of B3PW91/6-31G(d) yielded results that differed from macroscopic experimental values by 0.1–7.1 kcal/mol per atom or molecule, whereas applying B3PW91/cc-pVTZ produced results that differed from macroscopic experimental values by 0.1–4.8 kcal/mol per atom or molecule, with the smallest differences occurring for reactions with a small macroscopic experimental ?E and the largest differences occurring for reactions with a large macroscopic experimental ?E, implying size consistency.     

Contribution of left ventricular residual stress by myocytes and collagen: existence of inter-constituent mechanical interaction

We quantify the contribution of myocytes, collagen fibers and their interactions to the residual stress field found in the left ventricle (LV) using both experimental and theoretical methods. Ring tissue samples extracted from normal rat, male and female, LV were treated with collagenase and decellularization to isolate myocytes and collagen fibers, respectively. Opening angle tests were then performed on these samples as well as intact tissue samples containing both constituents that served as control. Our results show that the collagen fibers are the main contributor to the residual stress fields found in the LV. Specifically, opening angle measured in collagen-only samples (106.45\(^\circ \) ± 23.02\(^\circ \)) and myocytes-only samples (21.00\(^\circ \) ± 4.37\(^\circ \)) was significantly higher and lower than that of the control (57.88\(^\circ \) ± 12.29\(^\circ \)), respectively. A constrained mixture (CM) modeling framework was then used to infer these experimental results. We show that the framework cannot reproduce the opening angle found in the intact tissue with measurements made on the collagen-only and myocytes-only samples. Given that the CM framework assumes that each constituent contributes to the overall mechanics simply by their mere presence, this result suggests the existence of some myocyte–collagen mechanical interaction that cannot be ignored in the LV. We then propose an extended CM formulation that takes into account of the inter-constituent mechanical interaction in which constituents are deformed additionally when they are physically combined into a mixture. We show that the intact tissue opening angle can be recovered in this framework.     

Genetic variability and structure of an isolated population of <Emphasis Type="Italic">Ambystoma altamirani</Emphasis>, a mole salamander that lives in the mountains of one of the largest urban areas in the world

Amphibians are globally threatened by habitat loss and fragmentation; species within the order Ambystoma are not the exception, as there are 18 species of mole salamanders in México, of which 16 are endemic and all species are under some national or international status of protection. The mole salamander, Ambystoma altamirani is a microendemic species, which is distributed in central México, within the trans-Mexican volcanic belt, and is one of the most threatened species due to habitat destruction and the introduction of exotic species. Nine microsatellite markers were used to determine the genetic structure, genetic variability, effective population size, presence of bottlenecks and inbreeding coefficient of one population of A. altamirani to generate information which might help to protect and conserve this threatened species. We found two genetic subpopulations with significant level of genetic structure (\(F_{\mathrm{ST}}= 0.005\)) and high levels of genetic variability (\(H_{\mathrm{o}}= 0.883\); \(H_{\mathrm{e}}= 0.621\)); we also found a small population size (\(N_{\mathrm{e}} = 8.8\)), the presence of historical (\(M =\) 0.486) and recent bottlenecks under IAM and TPM models, with a low, but significant coefficient of inbreeding (\(F_{\mathrm{IS}} = -\)0.451). This information will help us to raise conservation strategies of this microendemic mole salamander species.     

Critical zone properties control the fate of nitrogen during experimental rainfall in montane forests of the Colorado Front Range

Several decades of research in alpine ecosystems have demonstrated links among the critical zone, hydrologic response, and the fate of elevated atmospheric nitrogen (N) deposition. Less research has occurred in mid-elevation forests, which may be important for retaining atmospheric N deposition. To explore the fate of N in the montane zone, we conducted plot-scale experimental rainfall events across a north–south transect within a catchment of the Boulder Creek Critical Zone Observatory. Rainfall events mimicked relatively common storms (20–50% annual exceedance probability) and were labeled with ¹⁵N-nitrate (\( {\text{NO}}_{3}^{ - } \)) and lithium bromide tracers. For 4 weeks, we measured soil–water and leachate concentrations of Br^?, \( {}^{15}{\text{NO}}_{3}^{ - } , \) and \( {\text{NO}}_{3}^{ - } \) daily, followed by recoveries of ¹⁵N species in bulk soils and microbial biomass. Tracers moved immediately into the subsurface of north-facing slope plots, exhibiting breakthrough at 10 and 30 cm over 22 days. Conversely, little transport of Br^? or \( {}^{15}{\text{NO}}_{3}^{ - } \) occurred in south-facing slope plots; tracers remained in soil or were lost via pathways not measured. Hillslope position was a significant determinant of soil ¹⁵N-\( {\text{NO}}_{3}^{ - } \) recoveries, while soil depth and time were significant determinants of ¹⁵N recovery in microbial biomass. Overall, ¹⁵N recovery in microbial biomass and leachate was greater in upper north-facing slope plots than lower north-facing (toeslope) and both south-facing slope plots in August; by October, ¹⁵N recovery in microbial N biomass within south-facing slope plots had increased substantially. Our results point to the importance of soil properties in controlling the fate of N in mid-elevation forests during the summer season.     

Twist buckling behavior of arteries

Arteries are often subjected to torsion due to body movement and surgical procedures. While it is essential that arteries remain stable and patent under twisting loads, the stability of arteries under torsion is poorly understood. The goal of this work was to experimentally investigate the buckling behavior of arteries under torsion and to determine the critical buckling torque, the critical buckling twist angle, and the buckling shape. Porcine common carotid arteries were slowly twisted in vitro until buckling occurred while subjected to a constant axial stretch ratio (1.1, 1.3, 1.5 (in vivo level) and 1.7) and lumen pressure (20, 40, 70 and 100 mmHg). Upon buckling, the arteries snapped to form a kink. For a group of six arteries, the axial stretch ratio significantly affected the critical buckling torque (\(p<0.002\)) and the critical buckling twist angle (\(p<0.001\)). Lumen pressure also significantly affected the critical buckling torque (\(p<0.001\)) but had no significant effect on the critical twist angle (\(p=0.067\)). Convex material constants for a Fung strain energy function were determined and fit well with the axial force, lumen pressure, and torque data measured pre-buckling. The material constants are valid for axial stretch ratios, lumen pressures, and rotation angles of 1.3–1.5, 20–100 mmHg, and 0–270\(^\circ \), respectively. The current study elucidates the buckling behavior of arteries under torsion and provides new insight into mechanical instability of blood vessels.     

Physical parameter estimation from porcine ex vivo vocal fold dynamics in an inverse problem framework

This study presents a framework for a direct comparison of experimental vocal fold dynamics data to a numerical two-mass-model (2MM) by solving the corresponding inverse problem of which parameters lead to similar model behavior. The introduced 2MM features improvements such as a variable stiffness and a modified collision force. A set of physiologically sensible degrees of freedom is presented, and three optimization algorithms are compared on synthetic vocal fold trajectories. Finally, a total of 288 high-speed video recordings of six excised porcine larynges were optimized to validate the proposed framework. Particular focus lay on the subglottal pressure, as the experimental subglottal pressure is directly comparable to the model subglottal pressure. Fundamental frequency, amplitude and objective function values were also investigated. The employed 2MM is able to replicate the behavior of the porcine vocal folds very well. The model trajectories’ fundamental frequency matches the one of the experimental trajectories in \(98.6\%\) of the recordings. The relative error of the model trajectory amplitudes is on average \(9.5\%\). The experiments feature a mean subglottal pressure of 10.16 (SD \(= 2.31\)) \({\text {cmH}}_2{\text {O}}\); in the model, it was on average 7.61 (SD \(= 2.40\)) \({\text {cmH}}_2{\text {O}}\). A tendency of the model to underestimate the subglottal pressure is found, but the model is capable of inferring trends in the subglottal pressure. The average absolute error between the subglottal pressure in the model and the experiment is 2.90 (SD \(= 1.80\)) \({\text {cmH}}_2{\text {O}}\) or \(27.5\%\). A detailed analysis of the factors affecting the accuracy in matching the subglottal pressure is presented.     

Competition of Multiple Species in Advective Environments

We study the effect of changes in flow speed on competition of an arbitrary number of species living in advective environments, such as streams and rivers. We begin with a spatial Lotka–Volterra model which is described by n reaction–diffusion–advection equations with Danckwerts boundary conditions. Using the dominant eigenvalue \(\lambda \le 0\) of the diffusion–advection operator subject to boundary conditions, we reduce the model to a system of ordinary differential equations. We impose a “transitive arrangement” of the competitors in terms of their interspecific coefficients and growth rates, which means that in the absence of advection, we have the following situation: for all \(1\le i<j\le n\), species i out-competes species j, while species j has higher intrinsic growth rate than species i. Changing advection speed in the original spatial model corresponds to changing the value of \(\lambda \) in the spatially implicit model. Considering the cases of the odd and even n separately, we obtain explicit intervals of the values of \(\lambda \) that allow all n species to be present in the habitat (coexistence interval). Stability of this equilibrium is shown for \(n\le 4\).     

Ultrasound palpation for fast in-situ quantification of articular cartilage stiffness,thickness and relaxation capacity

Most current cartilage testing devices require the preparation of excised samples and therefore do not allow intra-operative application for diagnostic purposes. The gold standard during open or arthroscopic surgery is still the subjective perception of manual palpation. This work presents a new diagnostic method of ultrasound palpation (USP) to acquire applied stress and strain data during manual palpation of articular cartilage. With the proposed method, we obtain cartilage thickness and stiffness. Moreover, repeated palpations allow the quantification of relaxation effects. USP measurements on elastomer phantoms demonstrated very good repeatability for both, stage-guided (97.2%) and handheld (96.0%) applications. The USP measurements were compared with conventional indentation experiments and revealed very good agreement on elastomer phantoms (\(r = 0.98\)) and good agreement on porcine cartilage samples (\(r = 0.76\)). Artificially degenerated cartilage samples showed reduced stiffness, weak capacity to relax after palpation and an increase of stiffness of approximately 50% with each single palpation. Intact cartilage was measured by USP directly at the patella (in situ) and after excision and removal of the subchondral bone (ex situ), leading to stiffness values of \(12.1\pm 5.5\) and \(8.5\pm 5.9\,\hbox {MPa}\) (\(p<0.05\)), respectively. The results demonstrate the potential of the USP system for cartilage testing, its sensitivity to degenerative changes and as a method for quantifying relaxation processes by means of repeated palpations. Furthermore, the differences in the results of in-situ and ex-situ measurements are of general interest, since such comparison has not been reported previously. We point out the limited comparability of ex-situ cartilage with its in-situ biomechanical behavior.     

Simulation of extracellular matrix remodeling by fibroblast cells in soft three-dimensional bioresorbable scaffolds

To culture functional soft tissues and organs in three-dimensional scaffolds, it is essential to elucidate the optimal scaffold mechanical properties. However, mechanoregulated soft tissue remodeling is not well understood. In this study, we hypothesized that individual cells are capable of remodeling extracellular matrix within a short proximity of themselves in order to match the stiffness of the broader surrounding matrix. This theory was implemented in a three-dimensional finite element model to simulate soft tissue remodeling of human fibroblast cells in two collagen–chitosan scaffolds with different mechanical properties. Simulation results closely matched with previously reported experimental data, showing that soft tissue growth in compliant (Scaf-A, 4.30 kPa) and stiff (Scaf-B, 17.03 kPa) scaffolds led to an almost eightfold difference in the resulting stiffnesses after 10 days (8.40 kPa for Scaf-A, 59.25 kPa for Scaf-B). Furthermore, varying the simulated rate for tissue remodeling by \(\pm \)50 % caused unequal changes in the resulting stiffness (+3.6 and \(-\)23 % for Scaf-A, +5 and \(-\)17 % for Scaf-B), and \(\pm \)50 % changes in the assumed upper limit on tissue stiffness only had larger effects on the stiff scaffold (+42 and \(-\)44 % for Scaf-B). These results reinforce the notion that soft tissue remodeling is governed by the stiffness of the surrounding matrix, until meeting an upper limit on tissue stiffness.     

Linear viscoelasticity - bone volume fraction relationships of bovine trabecular bone

Trabecular bone has been previously recognized as time-dependent (viscoelastic) material, but the relationships of its viscoelastic behaviour with bone volume fraction (BV/TV) have not been investigated so far. Therefore, the aim of the present study was to quantify the time-dependent viscoelastic behaviour of trabecular bone and relate it to BV/TV. Uniaxial compressive creep experiments were performed on cylindrical bovine trabecular bone samples (\(\textit{n}\,{=}\,13\)) at loads corresponding to physiological strain level of 2000 \({\upmu }{\upvarepsilon }\). We assumed that the bone behaves in a linear viscoelastic manner at this low strain level and the corresponding linear viscoelastic parameters were estimated by fitting a generalized Kelvin–Voigt rheological model to the experimental creep strain response. Strong and significant power law relationships (\(r^2\,{=}\,0.73,\ p\,{<}\,0.001\)) were found between time-dependent creep compliance function and BV/TV of the bone. These BV/TV-based material properties can be used in finite element models involving trabecular bone to predict time-dependent response. For users’ convenience, the creep compliance functions were also converted to relaxation functions by using numerical interconversion methods and similar power law relationships were reported between time-dependent relaxation modulus function and BV/TV.     

Cell Volume Regulation in the Proximal Tubule of Rat Kidney

We developed a dynamic model of a rat proximal convoluted tubule cell in order to investigate cell volume regulation mechanisms in this nephron segment. We examined whether regulatory volume decrease (RVD), which follows exposure to a hyposmotic peritubular solution, can be achieved solely via stimulation of basolateral K\(^+\) and \(\hbox {Cl}^-\) channels and \(\hbox {Na}^+\)–\(\hbox {HCO}_3^-\) cotransporters. We also determined whether regulatory volume increase (RVI), which follows exposure to a hyperosmotic peritubular solution under certain conditions, may be accomplished by activating basolateral \(\hbox {Na}^+\)/H\(^+\) exchangers. Model predictions were in good agreement with experimental observations in mouse proximal tubule cells assuming that a 10% increase in cell volume induces a fourfold increase in the expression of basolateral K\(^+\) and \(\hbox {Cl}^-\) channels and \(\hbox {Na}^+\)–\(\hbox {HCO}_3^-\) cotransporters. Our results also suggest that in response to a hyposmotic challenge and subsequent cell swelling, \(\hbox {Na}^+\)–\(\hbox {HCO}^-_3\) cotransporters are more efficient than basolateral K\(^+\) and \(\hbox {Cl}^-\) channels at lowering intracellular osmolality and reducing cell volume. Moreover, both RVD and RVI are predicted to stabilize net transcellular \(\hbox {Na}^+\) reabsorption, that is, to limit the net \(\hbox {Na}^+\) flux decrease during a hyposmotic challenge or the net \(\hbox {Na}^+\) flux increase during a hyperosmotic challenge.     

Optimal Therapy Scheduling Based on a Pair of Collaterally Sensitive Drugs

Despite major strides in the treatment of cancer, the development of drug resistance remains a major hurdle. One strategy which has been proposed to address this is the sequential application of drug therapies where resistance to one drug induces sensitivity to another drug, a concept called collateral sensitivity. The optimal timing of drug switching in these situations, however, remains unknown. To study this, we developed a dynamical model of sequential therapy on heterogeneous tumors comprised of resistant and sensitive cells. A pair of drugs (DrugA, DrugB) are utilized and are periodically switched during therapy. Assuming resistant cells to one drug are collaterally sensitive to the opposing drug, we classified cancer cells into two groups, \(A_\mathrm{R}\) and \(B_\mathrm{R}\), each of which is a subpopulation of cells resistant to the indicated drug and concurrently sensitive to the other, and we subsequently explored the resulting population dynamics. Specifically, based on a system of ordinary differential equations for \(A_\mathrm{R}\) and \(B_\mathrm{R}\), we determined that the optimal treatment strategy consists of two stages: an initial stage in which a chosen effective drug is utilized until a specific time point, T, and a second stage in which drugs are switched repeatedly, during which each drug is used for a relative duration (i.e., \(f \Delta t\)-long for DrugA and \((1-f) \Delta t\)-long for DrugB with \(0 \le f \le 1\) and \(\Delta t \ge 0\)). We prove that the optimal duration of the initial stage, in which the first drug is administered, T, is shorter than the period in which it remains effective in decreasing the total population, contrary to current clinical intuition. We further analyzed the relationship between population makeup, \(\mathcal {A/B} = A_\mathrm{R}/B_\mathrm{R}\), and the effect of each drug. We determine a critical ratio, which we term \(\mathcal {(A/B)}^{*}\), at which the two drugs are equally effective. As the first stage of the optimal strategy is applied, \(\mathcal {A/B}\) changes monotonically to \(\mathcal {(A/B)}^{*}\) and then, during the second stage, remains at \(\mathcal {(A/B)}^{*}\) thereafter. Beyond our analytic results, we explored an individual-based stochastic model and presented the distribution of extinction times for the classes of solutions found. Taken together, our results suggest opportunities to improve therapy scheduling in clinical oncology.     

Towards the prediction of flow-induced shear stress distributions experienced by breast cancer cells in the lymphatics

Tumour metastasis in the lymphatics is a crucial step in the progression of breast cancer. The dynamics by which breast cancer cells (BCCs) travel in the lymphatics remains poorly understood. The goal of this work is to develop a model capable of predicting the shear stresses metastasising BCCs experience using numerical and experimental techniques. This paper models the fluidic transport of large particles (\(\eta =d_{\mathrm{p}}/W=0.1-0.4\) where \(d_{\mathrm{p}}\) is the particle diameter and W is the channel width) subjected to lymphatic flow conditions (\({ Re}=0.04\)), in a \(100\times 100\,\upmu \hbox {m}\) microchannel. The feasibility of using the dynamic fluid body interaction (DFBI) method to predict particle motion was assessed, and particle tracking experiments were performed. The experiments found that particle translational velocity decreased from the undisturbed fluid velocity with increasing particle size (5–14% velocity lag for \(\eta =0.1-0.3\)). DFBI simulations were found to better predict particle behaviour than theoretical predictions; however, mesh restrictions in the near-wall region (\(0.2\,\mathrm{W}>y>0.8\,\mathrm{W}\)) result in computationally expensive models. The simulations were in good agreement with the experiments (\(<12\%\) difference) across the channel (\(0.2\,\mathrm{W}\le y\le 0.8\,\mathrm{W}\)), with differences up to 25% in the near-wall region. Particles experience a range of shear stresses (0.002–0.12 Pa) and spatial shear gradients (\(0.004-0.137\,\hbox {Pa}/\upmu \hbox {m}\)) depending on their size and radial position. The predicted shear gradients are far in excess of values associated with BCC apoptosis (\(0.004-0.023\,\hbox {Pa}/\upmu \hbox {m}\)). Increasing our understanding of the shear stress magnitudes and gradients experienced by BCCs could be leveraged to elucidate whether a particular BCC size or location exists that encourages metastasis within the lymphatics.     

Joint non-uniform sampling of all incremented time delays for quicker acquisition in protein relaxation studies

NMR relaxometry plays crucial role in studies of protein dynamics. The measurement of longitudinal and transverse relaxation rates of \(^{15}\)N is the main source of information on backbone motions. However, even the most basic approach exploiting a series of \(^{15}\)N HSQC spectra can require several hours of measurement time. Standard non-uniform sampling (NUS), i.e. random under-sampling of indirect time domain, typically cannot reduce this by more than 2–4\(\times\) due to relatively low “compressibility” of these spectra. In this paper we propose an extension of NUS to relaxation delays. The two-dimensional space of \(t_1\)/\(t_{relax}\) is sampled in a way similar to NUS of \(t_1\)/\(t_2\) domain in 3D spectra. The signal is also processed in a way similar to that known from 3D NUS spectra i.e. using one of the most popular compressed sensing algorithms, iterative soft thresholding. The 2D Fourier transform matrix is replaced with mixed inverse Laplace-Fourier transform matrix. The peak positions in resulting 3D spectrum are characterized by two frequency coordinates and relaxation rate and thus no additional fitting of exponential curves is required. The method is tested on three globular proteins, providing satisfactory results in a time corresponding to acquisition of two conventional \(^{15}\)N HSQC spectra.