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.     

Backbone resonance assignments for the SET domain of the human methyltransferase NSD2

Aberrant NSD2 methyltransferase activity is implicated as the oncogenic driver in multiple myeloma, suggesting opportunities for novel therapeutic intervention. The methyltransferase activity of NSD2 resides in its catalytic SET domain, which is conserved among most lysine methyltransferases. Here we report the backbone \(\hbox {H}^{\mathrm{N}}\), N, C\(^{\prime }\), \(\hbox {C}^\alpha\) and side-chain \(\hbox {C}^\beta\) assignments of a 25 kDa NSD2 SET domain construct, spanning residues 991–1203. A chemical shift analysis of C\(^{\prime }\), \(\hbox {C}^\alpha\) and \(\hbox {C}^\beta\) resonances predicts a secondary structural pattern that is in agreement with homology models.     

Delay decomposition approach to $$\mathcal {H}_{\infty }$$ filtering analysis of genetic oscillator networks with time-varying delays

In this paper, the \(\mathcal {H}_{\infty }\) filtering problem is treated for N coupled genetic oscillator networks with time-varying delays and extrinsic molecular noises. Each individual genetic oscillator is a complex dynamical network that represents the genetic oscillations in terms of complicated biological functions with inner or outer couplings denote the biochemical interactions of mRNAs, proteins and other small molecules. Throughout the paper, first, by constructing appropriate delay decomposition dependent Lyapunov–Krasovskii functional combined with reciprocal convex approach, improved delay-dependent sufficient conditions are obtained to ensure the asymptotic stability of the filtering error system with a prescribed \(\mathcal {H}_{\infty }\) performance. Second, based on the above analysis, the existence of the designed \(\mathcal {H}_{\infty }\) filters are established in terms of linear matrix inequalities with Kronecker product. Finally, numerical examples including a coupled Goodwin oscillator model are inferred to illustrate the effectiveness and less conservatism of the proposed techniques.     

A finite-element model of mechanosensation by a Pacinian corpuscle cluster in human skin

The Pacinian corpuscle (PC) is the cutaneous mechanoreceptor responsible for sensation of high-frequency (20–1000 Hz) vibrations. PCs lie deep within the skin, often in multicorpuscle clusters with overlapping receptive fields. We developed a finite-element mechanical model of one or two PCs embedded within human skin, coupled to a multiphysics PC model to simulate action potentials elicited by each PC. A vibration was applied to the skin surface, and the resulting mechanical signal was analyzed using two metrics: the deformation amplitude ratio (\({\rho }_{\mathrm{1S}} \), \({\rho }_{\mathrm{2S}} )\) and the phase shift of the vibration (\({\delta }_{\mathrm{S}1}^{\mathrm{mech}} \), \({\delta }_{\mathrm{S}2}^{\mathrm{mech}} )\) between the stimulus and the PC. Our results showed that the amplitude attenuation and phase shift at a PC increased with distance from the stimulus to the PC. Differences in amplitude (\(\rho _{12} )\) and phase shift (\({\delta }_{12}^{\mathrm{mech}} )\) between the two PCs in simulated clusters directly affected the interspike interval between the action potentials elicited by each PC (\({\delta }_{12}^{\mathrm{spike}} )\). While \({\delta }_{12}^{\mathrm{mech}} \) had a linear relationship with \({\delta }_{12}^{\mathrm{spike}} \), \(\rho _{12} \)’s effect on \({\delta }_{12}^{\mathrm{spike}} \) was greater for lower values of \(\rho _{12} \). In our simulations, the separation between PCs and the distance of each PC from the stimulus location resulted in differences in amplitude and phase shift at each PC that caused \({\delta }_{12}^{\mathrm{spike}} \) to vary with PC location. Our results suggest that PCs within a cluster receive different mechanical stimuli which may enhance source localization of vibrotactile stimuli, drawing parallels to sound localization in binaural hearing.     

Caspase-1-Mediated Pyroptosis of the Predominance for Driving CD4$$^{}$$ T Cells Death: A Nonlocal Spatial Mathematical Model

Caspase-1-mediated pyroptosis is the predominance for driving CD4\(^{+}\) T cells death. Dying infected CD4\(^{+}\) T cells can release inflammatory signals which attract more uninfected CD4\(^{+}\) T cells to die. This paper is devoted to developing a diffusive mathematical model which can make useful contributions to understanding caspase-1-mediated pyroptosis by inflammatory cytokines IL-1\(\beta \) released from infected cells in the within-host environment. The well-posedness of solutions, basic reproduction number, threshold dynamics are investigated for spatially heterogeneous infection. Travelling wave solutions for spatially homogeneous infection are studied. Numerical computations reveal that the spatially heterogeneous infection can make \(\mathscr {R}_0>1\), that is, it can induce the persistence of virus compared to the spatially homogeneous infection. We also find that the random movements of virus have no effect on basic reproduction number for the spatially homogeneous model, while it may result in less infection risk for the spatially heterogeneous model, under some suitable parameters. Further, the death of infected CD4\(^{+}\) cells which are caused by pyroptosis can make \(\mathscr {R}_0<1\), that is, it can induce the extinction of virus, regardless of whether or not the parameters are spatially dependent.     

Analysis of two susceptibility SNPs in HLA region and evidence of interaction between rs6457617 in <Emphasis Type="BoldItalic">HLA-DQB1</Emphasis> and <Emphasis Type="BoldItalic">HLA-DRB1</Emphasis>*04 locus on Tunisian rheumatoid arthritis

Previous genomewide association studies (GWAS) and meta-analyses have enumerated several genes/loci in major histocompatibility complex region, which are consistently associated with rheumatoid arthritis (RA) in different ethnic populations. Given the genetic heterogeneity of the disease, it is necessary to replicate these susceptibility loci in other populations. In this case, we investigate the analysis of two SNPs, rs13192471 and rs6457617, from the human leukocyte antigen (HLA) region with the risk of RA in Tunisian population. These SNPs were previously identified to have a strong RA association signal in several GWAS studies. A case–control sample composed of 142 RA patients and 123 healthy controls was analysed. Genotyping of rs13192471 and rs6457617 was carried out using real-time PCR methods by TaqMan allelic discrimination assay. A trend of significant association was found in rs6457617 TT genotype with susceptibility to RA (\(P = 0.04\), \(p_{c} = 0.08\), \(\hbox {OR} = 1.73\)). Moreover, using multivariable analysis, the combination of rs6457617*TT–HLA-DRB1*\(04^{+}\) increased risk of RA (\(\hbox {OR} = 2.38\)), which suggest a gene–gene interaction event between rs6457617 located within the HLA-DQB1 and HLA-DRB1. Additionally, haplotypic analysis highlighted a significant association of rs6457617*T–HLA-DRB1*\(04^{+}\) haplotype with susceptibility to RA (\(P = 0.018\), \(p_{c} = 0.036\), \(\hbox {OR} = 1.72\)). An evidence of association was shown subsequently in \(\hbox {antiCCP}^{+}\) subgroup with rs6457617 both in T allele and TT genotype (\(P = 0.01\), \(p_{c} = 0.03\), \(\hbox {OR} = 1.66\) and \(P = 0.008\), \(p_{c} = 0.024\), \(\hbox {OR} = 1.28\), respectively). However, no association was shown for rs13192471 polymorphism with susceptibility and severity to RA. This study suggests the involvement of rs6457617 locus as risk variant for susceptibility/severity to RA in Tunisian population. Secondly, it highlights the gene–gene interaction between HLA-DQB1 and HLA-DRB1.     

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.     

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.     

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.     

Cellular mechanosensitivity to substrate stiffness decreases with increasing dissimilarity to cell stiffness

Computational modelling has received increasing attention to investigate multi-scale coupled problems in micro-heterogeneous biological structures such as cells. In the current study, we investigated for a single cell the effects of (1) different cell-substrate attachment (2) and different substrate modulus \(\textit{E}_\mathrm{s}\) on intracellular deformations. A fibroblast was geometrically reconstructed from confocal micrographs. Finite element models of the cell on a planar substrate were developed. Intracellular deformations due to substrate stretch of \(\lambda =1.1\), were assessed for: (1) cell-substrate attachment implemented as full basal contact (FC) and 124 focal adhesions (FA), respectively, and \(\textit{E}_\mathrm{s}\,=\,\)140 KPa and (2) \(\textit{E}_\mathrm{s}\,=\,10\), 140, 1000, and 10,000 KPa, respectively, and FA attachment. The largest strains in cytosol, nucleus and cell membrane were higher for FC (1.35\(\text {e}^{-2}\), 0.235\(\text {e}^{-2}\) and 0.6\(\text {e}^{-2}\)) than for FA attachment (0.0952\(\text {e}^{-2}\), 0.0472\(\text {e}^{-2}\) and 0.05\(\text {e}^{-2}\)). For increasing \(\textit{E}_\mathrm{s}\), the largest maximum principal strain was 4.4\(\text {e}^{-4}\), 5\(\text {e}^{-4}\), 5.3\(\text {e}^{-4}\) and 5.3\(\text {e}^{-4}\) in the membrane, 9.5\(\text {e}^{-4}\), 1.1\(\text {e}^{-4}\), 1.2\(\text {e}^{-3}\) and 1.2\(\text {e}^{-3}\) in the cytosol, and 4.5\(\text {e}^{-4}\), 5.3\(\text {e}^{-4}\), 5.7\(\text {e}^{-4}\) and 5.7\(\text {e}^{-4}\) in the nucleus. The results show (1) the importance of representing FA in cell models and (2) higher cellular mechanical sensitivity for substrate stiffness changes in the range of cell stiffness. The latter indicates that matching substrate stiffness to cell stiffness, and moderate variation of the former is very effective for controlled variation of cell deformation. The developed methodology is useful for parametric studies on cellular mechanics to obtain quantitative data of subcellular strains and stresses that cannot easily be measured experimentally.     

Stability of a Stochastic Model of an SIR Epidemic with Vaccination

We prove almost sure exponential stability for the disease-free equilibrium of a stochastic differential equations model of an SIR epidemic with vaccination. The model allows for vertical transmission. The stochastic perturbation is associated with the force of infection and is such that the total population size remains constant in time. We prove almost sure positivity of solutions. The main result concerns especially the smaller values of the diffusion parameter, and describes the stability in terms of an analogue \(\mathcal{R}_\sigma\) of the basic reproduction number \(\mathcal{R}_0\) of the underlying deterministic model, with \(\mathcal{R}_\sigma \le \mathcal{R}_0\). We prove that the disease-free equilibrium is almost sure exponentially stable if \(\mathcal{R}_\sigma <1\).     

Structure,electronic properties,and NBO and TD-DFT analyses of nickel(II), zinc(II), and palladium(II) complexes based on Schiff-base ligands

In this work we studied the structural and electronic properties of the metal–Schiff base complexes Ni\( {\mathrm{L}}_2^2 \) (1), Pd\( {\mathrm{L}}_2^1 \) (2), Zn\( {\mathrm{L}}_2^2 \) (3), and Ni\( {\mathrm{L}}_2^1 \)(4), where L¹ and L² are Schiff bases synthesized from salicylaldehyde and 2-hydroxy-5-methylbenzaldehyde, respectively. Natural bond analysis showed that in complexes 1 and 2, the metal ion coordinates to the ligands through electron donation from lone pairs on ligand nitrogen and oxygen atoms to s and d orbitals on the metal ion. In complex 3, metal–N and metal–O bonds are formed through charge transfer from the lone pairs on nitrogen and oxygen atoms to an s orbital of Zn. Dimethylation of the phenolate rings in the ligands decreases the energy gap and redshifts the spectrum of the nickel complex. The main absorptions observed were assigned on the basis of singlet-state transitions. The simulated spectra of the two complexes 1 and 2 are characterized by excited states with ligand-to-ligand charge-transfer (LLCT), metal-to-ligand charge-transfer (MLCT), ligand-to-metal charge-transfer (LMCT), and metal-centered (MC) character.

Open image in new window

Graphical abstract Geometric structure of the palladium complex.

     

Isolation and characterization of microsatellite markers in the spiny lobster, <Emphasis Type="Italic">Panulirus echinatus</Emphasis> Smith, 1869 (Decapoda: Palinuridae) by Illumina MiSeq sequencing

Okra’s (Abelmoschus esculentus (L.) Moench) commercial cultivation is threatened in the tropics due to high incidence of yellow vein mosaic virus (YVMV) disease. Okra geneticists across the world tried to understand the inheritance pattern of YVMV disease tolerance without much success. Therefore, the inheritance pattern of YVMV disease in okra was revisited by employing six generations (\(\hbox {P}_{1}\), \(\hbox {P}_{2}\), \(\hbox {F}_{1}\), \(\hbox {F}_{2}\), \(\hbox {BC}_{1}\) and \(\hbox {BC}_{2}\)) of four selected crosses (one tolerant \(\times \) tolerant, two tolerant \(\times \) susceptible and one susceptible \(\times \) susceptible) using two tolerant (BCO-1 and Lal Bhendi) and two susceptible (Japanese Jhar Bhendi and PAN 2127) genotypes. Qualitative genetic analysis was done on the basis of segregation pattern of tolerant and susceptible plants in \(\hbox {F}_{2}\) and backcross generations of all the four crosses. It revealed that a single dominant gene along with some minor factors governed the disease tolerant trait in both the tolerant parents used. However, it was observed that genes governing disease tolerance identified in both the tolerant variety used was different. It could be concluded that the gene governing YVMV disease tolerance in okra was genotype specific. Further, duplicate gene action as evident from an approximate ratio of 15 : 1 (tolerant : susceptible) in the \(\hbox {F}_{2}\) population in the cross of two tolerant varieties gave a scope of increasing the tolerance level of the hybrid plants when both the tolerant genes are brought together. However, generation mean analysis revealed involvement of both additive and nonadditive effects in the inheritance of disease tolerance. Thus, the present study confirms that a complicated genetic inheritance pattern is involved in the disease tolerance against YVMV trait. The major tolerance genes could be transferred to other okra varieties, but the tolerance breaking virus strains might not allow them to achieve tolerance in stable condition. Therefore, accumulation of additional genes may be needed for a sustainable tolerance phenotype in okra.     

Three years of soil respiration in a mature eucalypt woodland exposed to atmospheric CO<Subscript>2</Subscript> enrichment

Large amounts of atmospheric N deposition cause negative effects on ecosystems. Effective mitigation strategies require the sources of N deposition to be identified and the contributions from individual sources to be quantified. Determination of the isotopic composition represents a useful approach in source apportionment. In this study, the δ¹⁵N-NH_x of wet and dry atmospheric deposition and the main NH₃ emission sources were analyzed at an urban, a suburban and a rural site in the Taihu Lake region of China. The 2-year average δ¹⁵N-\( {\text{NH}}_{4}^{ + } \) of precipitation was ? 3.0 ± 2.3, ? 3.1 ± 2.8 and ? 0.5 ± 2.8‰ for the urban, suburban and rural sites, respectively. These values were much lower than the corresponding values for particulate \( {\text{NH}}_{4}^{ + } \) (15.9, 15.2 and 14.3‰ at the urban, suburban and rural sites, respectively), and much higher than those of gaseous δ¹⁵N-NH₃ (? 16.7, ? 18.2 and ? 17.4‰ at the urban, suburban and rural sites, respectively). The δ¹⁵N-NH₃ of NH₃ from the main emission sources ranged from ? 30.8 to ? 3.3‰ for volatilized fertilizer, from ? 35.1 to ? 10.5‰ for emissions from a pig farm, and ? 24.7 to ? 11.3‰ for emissions from a dairy farm. Temporal variations of deposition δ¹⁵N-NH_x indicated that δ¹⁵N-NH_x values were lower in summer and autumn, but higher in winter and spring for both precipitation \( {\text{NH}}_{4}^{ + } \)-N and gaseous NH₃-N. Weather conditions such as temperature and precipitation significantly influenced the spatial and temporal distribution of isotope values of the deposition. Analysis of δ¹⁵N-NH_x in deposition and emission sources identified volatilized fertilizer and livestock wastes as the origins of both gaseous NH₃-N and precipitation \( {\text{NH}}_{4}^{ + } \)-N over the region. A stable isotope mixing model estimated that volatilized fertilizer and animal excreta contributed more than 65% to precipitation \( {\text{NH}}_{4}^{ + } \)-N, more than 60% to particulate \( {\text{NH}}_{4}^{ + } \)-N, and more than 75% to gaseous NH₃-N.     

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.     

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.     

A new indicator-based many-objective ant colony optimizer for continuous search spaces

In this paper, we propose a novel multi-objective ant colony optimizer (called iMOACO\(_{\mathbb {R}}\)) for continuous search spaces, which is based on ACO\(_{\mathbb {R}}\) and the R2 performance indicator. iMOACO\(_{\mathbb {R}}\) is the first multi-objective ant colony optimizer (MOACO) specifically designed to tackle continuous many-objective optimization problems (i.e., multi-objective optimization problems having four or more objectives). Our proposed iMOACO\(_{\mathbb {R}}\) is compared to three state-of-the-art multi-objective evolutionary algorithms (NSGA-III, MOEA/D and SMS-EMOA) and a MOACO algorithm called MOACO\(_{\mathbb {R}}\) using standard test problems and performance indicators taken from the specialized literature. Our experimental results indicate that iMOACO\(_{\mathbb {R}}\) is very competitive with respect to NSGA-III and MOEA/D and it is able to outperform SMS-EMOA and MOACO\(_{\mathbb {R}}\) in most of the test problems adopted.