Similarity of decay-associated spectra for tryptophan fluorescence of proteins with different structures

Tryptophan fluorescence lifetimes were analyzed for three proteins: human serum albumin, bovine serum albumin, and bacterial luciferase, which contain one, two, and seven tryptophan residues, respectively. For all of the proteins, the fluorescence decays were fitted by three lifetimes: τ₁ = 6–7 ns, τ₂ = 2.0–2.3 ns, and τ₃ ≤ 0.1 ns (the native state), and τ₁ = 4.4–4.6 ns, τ₂ = 1.7–1.8 ns, and τ₃ ≤ 0.1 ns (the denatured state). Corresponding decay-associated spectra had similar peak wavelengths and spectrum half-widths both in the native state (\(\lambda _{\max }^{{\tau _1}} = 324nm\), \(\lambda _{\max }^{{\tau _2}} = 328nm\), and \(\lambda _{\max }^{{\tau _3}} = 315nm\)), and in the denatured state (\(\lambda _{\max }^{{\tau _1}} = 350nm\), \(\lambda _{\max }^{{\tau _2}} = 343nm\), and \(\lambda _{\max }^{{\tau _3}} = 317nm\)). The differences in the steady-state spectra of the studied proteins were accounted for the individual ratio of the lifetime component contributions. The lifetime components were compared with a classification of tryptophan residues in the structure of these proteins within the discrete states model.     

The gene pool of the Belgorod oblast population: II. “Family name portraits” in groups of districts with different degrees of subdivision and the role of migrations in their formation

The frequencies and spectra of surnames have been analyzed in groups of raions (districts) of the Belgorod oblast (region) with different degrees of population subdivision. The “family name portraits” of districts with low (0.00003 < < f*_r < 0.00022, \(\overline {f_r^ * } \) = 0.00015) and moderate (0.00023 < f*_r < 0.00042, \(\overline {f_r^ * } \) = 0.00029) inbreeding levels are similar both to each other and to the “family name portrait” of the Belgorod oblast as a whole. Districts with high subdivision levels (0.00043 < f*_r < 0.00125, \(\overline {f_r^ * } \) = 0.00072) had very distinctive surname spectra and the highest surname frequencies. Intense immigration to the Belgorod oblast significantly affects its population genetic structure, decreasing the population subdivision.     

Efficient generation of negative hydrogen ions in a low-voltage cesium-hydrogen discharge

The possibility of achieving the high density of negative hydrogen ions \(N_{H^ - } \) in a low-voltage cesium-hydrogen discharge is investigated. The \(N_{H^ - } \) density is determined experimentally from the absorption of laser radiation due to the photodetachment of electrons from H^? ions. The discharge plasma is investigated by the probe technique. The populations of the excited states of Cs atoms are determined from their emission intensities. With an input power of W≈(15–25) W/cm² in the discharge, densities of \(N_{H^ - } \sim (10^{12} - 10^{13} )cm^{ - 3} \) are achieved. The self-consistent calculations of the plasma parameters in the discharge gap agree well with the experimental results. The absorption of laser radiation due to the photoionization of Cs atoms is investigated. It is shown that the role of this absorption mechanism is negligible.     

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.     

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

The basic flow equations of electrophysiology in the presence of chemical reactions: I. Development of the equations

Starting from the basic flux equation, it is possible to obtain an integral form relating the current componentsI _i at an arbitrary pointr ₂ to the distribution of mobilities and concentrationsc _i, potential forces\(\bar \mu \), and chemical productivityp _i without any restrictive assumptions such as constant mobilities, constant field, steady state, or electrical neutrality. The equation is

$$\begin{gathered} I_i (r_2 ) = G_i (r_2 )\left[ {\Delta \bar \mu _i - \int_{r_1 }^{r_2 } {z_i } FA\left( {p_i - dc_i /dt} \right)\left( {\frac{1}{{G_i (r)}}} \right)dr} \right]; \hfill \\ G_i (r) = 1/\int_{r_1 }^r {\frac{{dr}}{{z_i^2 F^2 c_i u_i }}.} \hfill \\ \end{gathered} $$     

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.     

Measurements of the plasma density in the FTU tokamak by a pulsed time-of-flight X-wave refractometer

On-line control over the plasma density in tokamaks (especially, in long-term discharges) requires reliable measurements of the averaged plasma density. For this purpose, a new method of density measurements—a pulsed time-of-flight plasma refractometry—was developed and tested in the T-11M tokamak. This method allows one to determine the averaged density from the measured time delay of nanosecond microwave pulses propagating through the plasma. For an O-wave, the measured time delay is proportional to the line-averaged density and is independent of the density profile (f?f _p ) τ_o ≈ k _o \(\tfrac{1}{{f^2 }}\mathop \smallint \limits_l \) N(x dx. A similar formula is valid for an X-wave: τ_X = ≈ k _x \(\tfrac{{f^2 + f_c^2 }}{{(f^2 - f_c^2 )^2 }}\mathop \smallint \limits_l \) N(x)dx. Here, f is the frequency of the probing wave, f _p is the plasma frequency, l= 4 a is the path length for two-pass probing in the equatorial plane, a is the plasma minor radius, k _O and k _X are numerical factors, f _c is the electron-cyclotron frequency at the axis of the plasma column, and f _p ?f _c , f. Measurements of the time delay provide the same information as plasma interferometry, though they do no employ the effect of interference. When the conditions f _p ?f _c , f are not satisfied, the measured time delay depends on the shape of the density profile. In this case, in order to determine the average density regardless of the density profile, it is necessary to perform simultaneous measurements at several probing frequencies in order to determine the average density. In ITER (Bt ～ 5T), a spectral window between the lower and upper cutoff frequencies in the range of 50–100 GHz can be used for pulsed time-of-flight X-wave refractometry. This appreciably simplifies the diagnostics and eliminates the problem of the first mirror. In this paper, the first results obtained in the FTU tokamak with a prototype of the ITER pulsed time-of-flight refractometer are presented. The geometry and layout of experiments similar to the planned ITER experiments are described. The density measured by pulsed time-of-flight refractometry is shown to agree well with the results obtained in FTU with a two-frequency scanning IR interferometer. The results obtained are analyzed, and the future experiments are discussed.     

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.     

Generation of ion-acoustic and magnetoacoustic waves in an RF helicon discharge

A study is made of the generation of ion-acoustic and magnetoacoustic waves in a discharge excited in an external magnetic field by an electromagnetic wave in the whistler frequency range (ω_LH ? ω ? ω_He, where ω_LH = $\sqrt {\omega _{He} \omega _{Hi} } $ and ω_He and ω_Hi are the electron and ion gyrofrequencies, respectively). The excitation of acoustic waves is attributed to the decay of a high-frequency hybrid mode forming a plasma waveguide into low-frequency acoustic waves and new high-frequency waves that satisfy both the decay conditions and the waveguide dispersion relations. The excitation of acoustic waves is resonant in character because the conditions for the generation of waveguide modes and for the occurrence of the corresponding nonlinear wave processes should be satisfied simultaneously. An unexpected effect is the generation of magnetoacoustic waves by whistlers. A diagnostic technique is proposed that allows one to determine the thermal electron velocity by analyzing decay conditions and dispersion relations for waves in the discharge channel.     

Theoretical Study of Sensitivity and Localized Surface Plasmon Resonance of Ag-Dielectric Core-Shell Multi-layered Nanosphere

Localized surface plasmon resonances (LSPRs) of Ag-dielectric-Ag multi-layered nanoshell are studied by quasi-static approximation and plasmon hybridization theory. Absorption properties of multi-layered nanoshell with the silver core and nanoshell separated by a dielectric layer exhibit strong coupling between the core and nanoshell. The result shows absorption spectrum of LSPRS is influenced by the refractive index of surrounding medium, the dielectric constant of middle dielectric layer, the thickness of inner core radius and outer shell layer. LSPR shift of the longest wavelength \(\left |\omega _{-}^{-}\right >\) is red-shifted with increasing the inner core radius. It is interesting to find that longer wavelength \(\left |\omega _{-}^{+}\right >\) mode is mainly effected by the ratio constant of the surrounding medium refractive index ε ₄ to the middle layer dielectric constant ε ₂. \(\left |\omega _{-}^{+}\right >\) mode takes place a blue-shift with increasing inner core radius when ε ₂ > ε ₄, a red-shift when ε ₂ < ε ₄, and no-shifting when ε ₂ = ε ₄. However, the influence of dielectric layer radius to \(\left |\omega _{-}^{+}\right >\) mode shows the different property as that of increasing the inner core radius. The underlying mechanisms are analyzed with the plasmon hybridization theory and the distribution of induced charge interaction between the inner core and outer shell. In addition, the influence of core radius, middle dielectric layer radius and outer shell radius to sensitivity of Ag-dielectric-Ag multi-layered nanoshell are also reported, a higher sensitivity could be gotten by adjusting geometrical parameters. Our theoretical study could give an easy way to analyze properties of the core-shell nanosphere based on plasmon hybridization theory and the induced charge interaction, and usefully broaden the applications in nano-optics.     

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

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.     

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.     

Theoretical Studies of Tunable Localized Surface Plasmon Resonance of Gold-Dielectric Multilayered Nanoshells

Tunable properties of localized surface plasmon resonances (LSPR) of gold-dielectric multilayered nanoshells are studied by quasi-static theory and plasmon hybridization theory. Multilayered nanoshells with the gold core and nanoshell separated by a spacer layer exhibit strong coupling between the core and nanoshell plasmon resonance modes. It is found that the absorption spectra characteristics of LSPR are sensitive to multiple parameters including the surrounding medium refractive index, the dielectric constant of spacer layer, the radius of inner core gold sphere, outer shell layer thickness, and their coupling strength. The results show that LSPR is mainly influenced by the ratio of spacer layer dielectric constant ε ₂ to surrounding medium dielectric constant ε ₄. Absorption spectrum of \(\left |\omega _{-}^{+}\right \rangle \) mode is red-shifted with increasing core radius when ε ₂ > ε ₄. It is surprising to find that LSPR is blue-shifted with increasing core radius when ε ₂ < ε ₄, and no shift when ε ₂ = ε ₄. These interesting contrary shifts of \(\left |\omega _{-}^{+}\right \rangle \) mode with different ratios ε ₂/ε ₄ are well analysed with plasmon hybridization theory and the distributions of induced charges interaction between the inner core and outer shell. In addition, for the sake of clarity, the distributions of electric filed intensity at their plasmon resonance wavelengths are also calculated. This work may provide an alternative approach to analyse property of the core-shell nanoshell particles based on plasmon hybridization theory and the induced charge interaction.     

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.