Computational study on neuronal activities arising in the pre-Bötzinger complex

Experimental investigations have shown that the pre-Bötzinger complex (pre-BötC) within the mammalian brainstem generates the inspiratory phase of respiratory rhythm. Based on a single-compartment model of a pre-BötC inspiratory neuron, we, in this paper, use semi-analytical, numerical as well as fast-slow dynamical methods to investigate the effects of sodium conductance (\(g_{\text{Na}}\)) and potassium conductance (\(g_{{\text{K}}}\)) on the firing activities of pre-BötC and try to reveal the dynamical mechanisms behind them. We show how \(g_{{\text{Na}}}\) and \(g_{\text{K}}\) affect the bifurcations of the fast-subsystem and how the the firing patterns of pre-BötC transit according to the bifurcations.     

Influence of muscle groups’ activation on proximal femoral growth tendency

Muscle and joint contact force influence stresses at the proximal growth plate of the femur and thus bone growth, affecting the neck shaft angle (NSA) and femoral anteversion (FA). This study aims to illustrate how different muscle groups’ activation during gait affects NSA and FA development in able-bodied children. Subject-specific femur models were developed for three able-bodied children (ages 6, 7, and 11 years) using magnetic resonance images. Contributions of different muscle groups—hip flexors, hip extensors, hip adductors, hip abductors, and knee extensors—to overall hip contact force were computed. Specific growth rate for the growth plate was computed, and the growth was simulated in the principal stress direction at each element in the growth front. The predicted growth indicated decreased NSA and FA (of about \(0.1 {^{\circ }}\) over a four-month period) for able-bodied children. Hip abductors contributed the most, and hip adductors, the least, to growth rate. All muscles groups contributed to a decrease in predicted NSA (\(\sim \)0.01\({^{\circ }}\)–0.04\({^{\circ }})\) and FA (\(\sim \)0.004\({^{\circ }}\)–\(0.2{^{\circ }}\)), except hip extensors and hip adductors, which showed a tendency to increase the FA (\(\sim \)0.004\({^{\circ }}\)–\(0.02{^{\circ }}\)). Understanding influences of different muscle groups on long bone growth tendency can help in treatment planning for growing children with affected gait.     

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.     

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.     

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}\)).     

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.     

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.     

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.     

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.     

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.     

Trend analysis of variations in carbon stock using stock big data

Changes in land use affect the terrestrial carbon stock through changes in the land cover. Research on land use and analysis of variations in carbon stock have practical applications in the optimization of land use and the mitigation of climate change effects. This study was conducted in Baixiang and Julu counties in the Taihang Piedmont by employing the trend analysis method to characterize the variation in county land use and carbon stock. The findings show that in both counties, agricultural and unused land areas decreased while built-up land area increased, and the reduction in cropland was the main reason behind the agricultural land reduction. An inflection point appeared on the cropland curves of Julu, because the cropland area decreased by 1576.97 hm\(^{2}\) from 2004 to 2006. Cropland area in Baixiang decreased from 1996 to 1998 by a total of 129.89 hm\(^{2}\) and then remained relatively stable after 1998. The total carbon storage and variation in land use in the two counties displayed similar trends. Total carbon reserves in Julu increased by 2.76 \(\times \) 10\(^{4}\) tC (carbon equivalent), while those in Baixiang decreased by 0.63 \(\times \) 10\(^{4}\) tC. Carbon stock of built-up land in Julu and Baixiang increased by 2.44 \(\times \) 10\(^{4}\) and 1.22 \(\times \) 10\(^{4}\) tC, respectively.     

Glutathione <Emphasis Type="BoldItalic">S</Emphasis>-transferase P1 gene polymorphisms and susceptibility to coronary artery disease in a subgroup of north Indian population

The present study aimed to investigate the association of \(\hbox {g}.313\hbox {A}{>}\hbox {G}\) and \(\hbox {g}.341\hbox {C}{>}\hbox {T}\) polymorphisms of GSTP1 with coronary artery disease (CAD) in a subgroup of north Indian population. In the present case–control study, CAD patients (\(n = 200\)) and age-matched, sex-matched and ethnicity-matched healthy controls (\(n = 200\)) were genotyped for polymorphisms in GSTP1 using polymerase chain reaction-restriction fragment length polymorphism (PCR-RFLP) method. Genotype distribution of \(\hbox {g}.313\hbox {A}{>}\hbox {G}\) and \(\hbox {g}.341\hbox {C}{>}\hbox {T}\) polymorphisms of GSTP1 gene was significantly different between cases and controls (\(P = 0.005\) and 0.024, respectively). Binary logistic regression analysis showed significant association of A/G (odds ratio (OR): 1.6, 95% CI: 1.08–2.49, \(P = 0.020\)) and G/G (OR: 3.1, 95% CI: 1.41–6.71, P \(=\) 0.005) genotypes of GSTP1 \(\hbox {g}.313\hbox {A}{\!>\!}\hbox {G}\), and C/T (OR: 5.8, 95% CI: 1.26–26.34, \(P = 0.024\)) genotype of GSTP1 \(\hbox {g}.341\hbox {C}{>}\hbox {T}\) with CAD. The A/G and G/G genotypes of \(\hbox {g}.313\hbox {A}{>}\hbox {G}\) and C/T genotype of \(\hbox {g}.341\hbox {C}{>}\hbox {T}\) conferred 6.5-fold increased risk for CAD (OR: 6.5, 95% CI: 1.37–31.27, \(P = 0.018\)). Moreover, the recessive model of GSTP1 \(\hbox {g}.313\hbox {A}{>}\hbox {G}\) is the best fit inheritance model to predict the susceptible gene effect (OR: 2.3, 95% CI: 1.11–4.92, \(P = 0.020\)). In conclusion, statistically significant associations of GSTP1 \(\hbox {g}.313\hbox {A}{>}\hbox {G}\) (A/G, G/G) and \(\hbox {g}.341\hbox {C}{>}\hbox {T}\) (C/T) genotypes with CAD were observed.     

Acclimation to hypercarbia protects cardiac contractility and alters tissue carbohydrate metabolism in the Amazonian armored catfish <Emphasis Type="Italic">Pterygoplichthys pardalis</Emphasis>

The armored catfish Pterygoplichthys pardalis tolerates environmental hypercarbia, high partial pressures of CO₂ (\(P_{{{\text{CO}}_{ 2} }}\)), by preferentially protecting intracellular pH (pHi) in the face of extracellular acidosis. This response is associated with ionic changes which may disrupt contractility in cardiac muscle, and it is not known whether acclimation to hypercarbia provides protection against these changes. We studied the influence of different \(P_{{{\text{CO}}_{ 2} }}\) acclimation histories on cardiac muscle function using isometrically contracting ventricular strip preparations. Fish were held for >4 months at 21 mmHg \(P_{{{\text{CO}}_{ 2} }}\) and then exposed to normocarbia (6 mmHg \(P_{{{\text{CO}}_{ 2} }}\)) for either 15 h or 5–6 days. Acclimation to chronic hypercarbia eliminated the negative inotropic effects of in vitro hypercarbia, decreased extracellular Ca²⁺ sensitivity, and reduced maximum pacing frequency in ventricular strip preparations. Fish acclimated to chronic hypercarbia also exhibited hepatic glycogen and plasma glucose accumulation, and lower plasma lactate levels compared to fish acclimated to normocarbia for 5–6 days. We suggest chronic hypercarbia may induce cardiac remodeling to protect contractility and reduce the energetic demands of pHi regulation. The activation of HCO₃ ^? synthesis pathways may decrease glucose utilization and enhance carbohydrate stores, potentially providing protection against hypoxia, a stressor frequently encountered in conjunction with hypercarbia in the Amazon.     

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.     

Potential use of sulfite as a supplemental electron donor for wastewater denitrification

Biological denitrification typically requires the addition of a supplemental electron donor, which can add a significant operating expense to wastewater treatment facilities. Most common electron donors are organic, but reduced inorganic sulfur compounds (RISCs), such as sulfide (HS^?) and elemental sulfur (S⁰), may be more cost-effective. S⁰ is an inexpensive and well characterized electron donor, but it provides slow denitrification rates due to its low solubility. A lesser-known RISC is sulfite (\({\text{SO}}_{3}^{2 - }\)), which can be easily produced from S⁰ by a simple combustion process. Unlike S⁰, \({\text{SO}}_{3}^{2 - }\) is highly soluble, and therefore may provide higher denitrification rates. However, very little is known about microbial denitrification with \({\text{SO}}_{3}^{2 - }\). Also, \({\text{SO}}_{3}^{2 - }\) is a strong reductant that reacts abiotically with oxygen and has toxic effects on microorganisms. This paper reviews \({\text{SO}}_{3}^{2 - }\) in the environment, \({\text{SO}}_{3}^{2 - }\) chemistry, microbiology, toxicity, and its potential use for denitrification. Since \({\text{SO}}_{3}^{2 - }\) is an intermediate in the sulfur oxidation pathway of most sulfur-oxidizing microorganisms, it is an energetic electron donor and it should select for a \({\text{SO}}_{3}^{2 - }\)-oxidizing community. Our review of the literature, as well as our own lab experience, suggests that \({\text{SO}}_{3}^{2 - }\) can effectively serve as an electron donor for denitrification. Further research is needed to determine the kinetics of \({\text{SO}}_{3}^{2 - }\)-based denitrification, its toxic threshold for sulfur-oxidizing microorganisms, and its potential inhibition of sensitive species such as nitrifying microorganisms and potential formation of nitrous oxide. Its effect on sludge settling efficiency also should be explored.     

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.     

<Superscript>1</Superscript>H, <Superscript>13</Superscript>C, <Superscript>15</Superscript>N resonance assignments of the extracellular loop 1 domain (ECL1) of <Emphasis Type="Italic">Streptococcus pneumoniae</Emphasis> D39 FtsX,an essential cell division protein

FtsX is an integral membrane protein from Streptococcus pneumoniae (pneumococcus) that harbors an extracellular loop 1 domain (\({\text{FtsX}}^{\text{ECL1}}_{Spn}\)) that interacts with PcsB, an peptidoglycan hydrolase that is essential for cell growth and division. Here, we report nearly complete backbone and side chain resonance assignments and a secondary structural analysis of \({\text{FtsX}}^{\text{ECL1}}_{Spn}\) (residues 47–168 of FtsX) as first steps toward structure determination of \({\text{FtsX}}^{\text{ECL1}}_{Spn}\).     

Phylogenetics beyond biology

Evolutionary processes have been described not only in biology but also for a wide range of human cultural activities including languages and law. In contrast to the evolution of DNA or protein sequences, the detailed mechanisms giving rise to the observed evolution-like processes are not or only partially known. The absence of a mechanistic model of evolution implies that it remains unknown how the distances between different taxa have to be quantified. Considering distortions of metric distances, we first show that poor choices of the distance measure can lead to incorrect phylogenetic trees. Based on the well-known fact that phylogenetic inference requires additive metrics, we then show that the correct phylogeny can be computed from a distance matrix \({\mathbf {D}}\) if there is a monotonic, subadditive function \(\zeta\) such that \(\zeta ^{-1}({\mathbf {D}})\) is additive. The required metric-preserving transformation \(\zeta\) can be computed as the solution of an optimization problem. This result shows that the problem of phylogeny reconstruction is well defined even if a detailed mechanistic model of the evolutionary process remains elusive.