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.     

WNT10A variants in relation to nonsyndromic hypodontia in eastern Slovak population

Nonsyndromic hypodontia is a congenital absence of less than six permanent teeth, with a most common subtype maxillary lateral incisor agenesis (MLIA). Mutations in several genes have been described in severe tooth agenesis. The aim of this study was to search for the variants in wingless-type MMTV-integration site family member (WNT10A), paired box 9 (PAX9) and axis inhibitor 2 (AXIN2) genes, and investigate their potential role in the pathogenesis of non-syndromic hypodontia. Clinical examination and panoramic radiograph were performed in the cohort of 60 unrelated Slovak patients of Caucasian origin with nonsyndromic hypodontia including 37 MLIA cases and 48 healthy controls. Genomic DNA was isolated from buccal swabs and Sanger sequencing of WNT10A, PAX9 and AXIN2 was performed. Altogether, we identified 23 single-nucleotide variants, of which five were novel. We have found three rare nonsynonymous variants in WNT10A (p.Gly165Arg; p.Gly213Ser and p.Phe228Ile) in eight (13.33%) of 60 patients. Analysis showed potentially damaged WNT10A variant p.Phe228Ile predominantly occurred only in MLIA patients, and with a dominant form of tooth agenesis (odds ratio \(({\hbox {OR}}_{\mathrm{dom}}) = 9.841\); \(P=0.045\); 95% confidence interval (CI) 0.492–196.701; \({\hbox {OR}}_{\mathrm{rec}} = 0.773\); \(P =1.000\); 95% CI 0.015–39.877). In addition, the WNT10A variant p.Phe228Ile showed a trend associated with familial nonsyndromic hypodontia (\(P =0.024\); OR = 1.20; 95% CI 0.97–1.48). After Bonferroni correction, these effects remained with borderline tendencies. Using a 3D WNT10A protein model, we demonstrated that the variant Phe228Ile changes the protein secondary structure. In PAX9 and AXIN2, common variants were detected. Our findings suggest that the identified WNT10A variant p.Phe228Ile could represent risk for the inherited nonsyndromic hypodontia underlying MLIA. However, further study in different populations is required.     

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.     

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.     

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.     

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.

     

Forecasting next season’s <Emphasis Type="Italic">Ixodes ricinus</Emphasis> nymphal density: the example of southern Germany 2018

The castor bean tick, Ixodes ricinus (L.) (Ixodida: Ixodidae), is the principal vector of pathogens causing tick-borne encephalitis or Lyme borreliosis in Europe. It is therefore of general interest to make an estimate of the density of I. ricinus for the whole year at the beginning of the tick season. There are two necessary conditions for making a successful prediction: a long homogeneous time series of observed tick density and a clear biological relationship between environmental predictors and tick density. A 9-year time series covering the period 2009–2017 of nymphal I. ricinus flagged at monthly intervals in southern Germany has been used. With the hypothesis that I. ricinus density is triggered by the fructification of the European beech 2 years before, the mean annual temperature of the previous year, and the current mean winter temperature (December–February), a forecast of the annual nymphal tick density has been made. Therefore, a Poisson regression model was generated resulting in an explained variance of 93.4% and an error of \(\hbox {RMSE} = 21\) ticks per \(100\,\hbox {m}^2\) (annual \(\hbox {MEAN} = 260\) collected ticks/\(100\,\hbox {m}^2\)). An independent verification of the forecast for the year 2017 resulted in 187 predicted versus 180 observed nymphs per \(100\,\hbox {m}^2\). For the year 2018 a relatively high number of 443 questing I. ricinus nymphs per \(100\,\hbox {m}^2\) is forecasted, i.e., a “good” tick year.     

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.     

Understanding the effects of isolation on seed and pollen flow,spatial genetic structure and effective population size of the dioecious tropical tree species <Emphasis Type="Italic">Myracrodruon urundeuva</Emphasis>

This study examines the levels of gene flow, the distance and the patterns of pollen and seed dispersal, the intra-population spatial genetic structure (SGS) and the effective population size of a spatially isolated Myracrodruon urundeuva population using five microsatellite loci. The study was carried out in the Paulo de Faria Ecological Station, São Paulo State, Brazil and included the sampling and mapping of 467 adult-trees and 149 juveniles. Open-pollinated seeds (514) from 29 seed-trees were also sampled and genotyped. Significant SGS was detected in both adult (S _p  = 0.0269) and juveniles trees (S _p  = 0.0246), indicating short-distance seed dispersal. Using maternity analysis, all juveniles had the mother-tree assigned within the stand. A father-tree within the stand was also assigned for 97.3% of the juveniles and 98.4% of offspring. The average pollen dispersal distance measured in juveniles \( \left( {\hat{\delta } = 1 3 8\pm 1 6 9 {\text{ m}},{\text{ mean}} \pm {\text{SD}}} \right) \) and offspring \( \left( {\hat{\delta } = 2 5 2\pm 20 4 {\text{ m}}} \right) \) were higher than the average seed dispersal distance measured in juveniles \( \left( {\hat{\delta } = 1 2 4\pm 1 50{\text{ m}}} \right) \). About 70% of the pollen from juveniles and 51% from offspring traveled less than 200 m and, 72% of the seeds traveled less than 50 m. The effective population size of the studied sample indicates that the 467 adult-trees and 145 juveniles correspond respectively to 335 and 63 individuals that are neither inbred nor relatives. The results are discussed in relation to their impact on seed collection practices and genetic conservation.     

Association of <Emphasis Type="BoldItalic">lactase</Emphasis> 13910 C/T polymorphism with bone mineral density and fracture risk: a meta-analysis

A number of studies have investigated the association of lactase (LCT, C/T-13910) gene polymorphism with bone mineral density (BMD) and fracture risk, but previous results were inconclusive. In this study, a meta-analysis was performed to quantify the association of LCT (C/T-13910) polymorphism with BMD and fracture risk. Eligible publications were searched in the PubMed, Web of Science, Embase databases, Google Scholar, Yahoo and Baidu. Pooled weighed mean difference (WMD) or odds ratio (OR) with their 95% confidence interval (CI) were calculated using a fixed-effects or random-effects model. A total of nine articles with 8871 subjects were investigated in the present meta-analysis. Overall, the TT/TC genotypes of LCT 13910 C/T polymorphism showed significantly higher BMD than those with the CC genotype at femur neck (FN) (\(\hbox {WMD} = 0.011\,\hbox {g/cm}^{2}\), 95% CI \(=\) 0.004–0.018, \(P = 0.003\)). Besides, LCT 13910 C/T polymorphism may decrease the risk of any site fractures (for TT versus TC \(+\) CC, OR \(=\) 0.813, 95% CI \(=\) 0.704–0.938, \(P = 0.005\); for T allele versus C allele, OR \(=\) 0.885, 95% CI \(=\) 0.792–0.989, \(P = 0.032\)). However, there was no significant association of LCT 13910 C/T polymorphism with BMD at lumbar spine and risk of vertebral fractures under all genetic contrast models (all P values were \({>}0.05\)). The meta-analysis suggests that there are significant effects of LCT 13910 C/T polymorphism on BMD and fracture risk. Large-scale studies with different ethnic populations will be needed to further investigate the possible race-specific effect of LCT 13910 C/T polymorphism on BMD and fracture risk.     

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 computational investigation on the antioxidant potential of myricetin 3,4′-di-<Emphasis Type="Italic">O</Emphasis>-<Emphasis Type="Italic">α</Emphasis>-<Emphasis Type="Italic">L</Emphasis>-rhamnopyranoside

In this work, we present a computational study on the antioxidant potential of myricetin 3,4\(^{\prime }\)-di-O-α-L-rhamnopyranoside (Compound M). A density functional theory (DFT) approach with the B3LYP and LC-ωPBE functionals and with both the 6-311G(d,p) and 6-311+G(d,p) basis sets was used. The focus of the investigation was on the structural and energetic parameters including both bond dissociation enthalpies (BDEs) and ionization potentials (IPs), which provide information on the potential antioxidant activity. The properties computed were compared with BDEs and IPs available in the literature for myricetin, a compound well known for presenting antioxidant activity (and the parent molecule of the compound of interest in the present work). Myricetin 3,4\(^{\prime }\)-di-O-α-L-rhamnopyranoside presented the lowest BDE to be 79.13 kcal/mol (as determined using B3LYP/6-311G(d,p) in water) while myricetin has a quite similar value (within 3.4 kcal/mol). IPs computed in the gas phase [B3LYP/6-311G(d,p)] are 157.18 and 161.4 kcal/mol for myricetin 3,4\(^{\prime }\)-di-O-α-L-rhamnopyranoside and myricetin, respectively. As the values of BDEs are considerably lower than the ones probed for IPs (in the gas phase or in any given solvent environment), the hydrogen atom transfer mechanism is preferred over the single electron transfer mechanism. The BDEs obtained suggest that myricetin 3,4\(^{\prime }\)-di-O-α-L-rhamnopyranoside can present antioxidant potential as good as the parent molecule myricetin (a well-known antioxidant). Therefore, experimental tests on the antioxidant activity of Compound M are encouraged.     

Invasions Slow Down or Collapse in the Presence of Reactive Boundaries

Motivated by the propagation of thin bacterial films around planar obstacles, this paper considers the dynamics of travelling wave solutions to the Fisher–KPP equation \(u_t = u(1-u) + u_{xx} + u_{yy}\) in a planar strip \(-\infty< x < \infty \), \(0 \le y \le L\). We examine the propagation of fronts in the presence of a mixed boundary condition (also referred to as a ‘partially absorbing’ or ‘reactive’ boundary) \(u_y = \alpha u\), with \(\alpha >0\), at \(y=0\). The presence of boundary conditions of this kind leads to the development of front solutions that propagate in x but contain transverse structure in y. Motivated by the observation that the speed of propagation in the Fisher–KPP equation is determined (for exponentially decaying initial conditions) by the behaviour at the leading edge, we analyse the linearised Fisher–KPP equation in order to estimate the speed of the stable travelling front, a function of the width L and the imposed boundary conditions. For wide strips the speed estimate based on the linearised equation agrees well with the results of numerical simulations. For narrow channels numerical simulations indicate that the stable front propagates more slowly, and for sufficiently small L or sufficiently large \(\alpha \) the front speed falls to zero and the front collapses. The reason for the collapse is the non-existence, far behind the front, of a stable positive equilibrium solution u(x, y). While existence of these equilibrium states can be demonstrated via phase plane arguments, the investigation of stability is similar to calculations of critical patch sizes carried out in similar ecological models.     

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.     

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

Oestrogen regulates the expression of cathepsin E-A-like gene through ER$$\upbeta $$ in liver of chicken ( Gallus gallus)

The cathepsin E-A-like, also known as ‘similar to nothepsin’, is a new member of the aspartic protease family, which may take part in processing of egg yolk macromolecules, due to it was identified in the chicken egg-yolk. Previously, studies have suggested that the expression of cathepsin E-A-like increased gradually during sexual maturation of pullets, but the exact regulation mechanism is poorly understood. In this study, to gain insight into the function and regulation mechanism of the gene in egg-laying hen, we cloned the cathepsin E-A-like gene and evaluated its evolutionary origin by using both phylogenetic and syntenic methods. The mode of the gene expression regulation was analysed through stimulating juvenile hens with \(17\upbeta \)-estradiol and chicken embryo hepatocytes with \(17\upbeta \)-estradiol combined with oestrogen receptor antagonists including MPP, ICI 182,780 and tamoxifen. Our results showed that cathepsin E-A-like was an orthologoues gene with nothepsin, which is present in birds but not in mammals. The expression of cathepsin E-A-like significantly increased in a dose-dependent manner after the juvenile hens were treated with \(17\upbeta \)-estradiol (\(P~<~0.05\)). Compared with the \(17\upbeta \)-estradiol treatment group, the expression of cathepsin E-A-like was not significantly changed when the hepatocytes were treated with \(17\upbeta \)-estradiol combined with MPP (\(P~<~0.05\)). In contrast, compared with the \(17\upbeta \)-estradiol combined with MPP treatment group, the expression of cathepsin E-A-like was significantly downregulated when the hepatocytes were treated with \(17\upbeta \)-estradiol combined with tamoxifen or ICI 182,780 (\(P~<~0.05\)). These results demonstrated that cathepsin E-A-like shared the same evolutionary origin with nothepsin. The expression of cathepsin E-A-like was regulated by oestrogen, and the regulative effect was predominantly mediated through ER-\(\upbeta \) in liver of chicken.     

Energetics and reactivity of small beryllium deuterides

Enthalpies and free energies of reaction for small neutral and charged beryllium deuterides BeD, BeD₂, and BeD₃ that have been calculated are reported for a temperature range of 0 K to 1000 K. We discuss probable dissociation channels and possible ways of producing BeD by localizing the relevant transition states and by calculating corresponding rate constants. BeD and BeD⁺ are found to be the most stable ones among the considered compounds. BeD₂ and \( {\mathrm{BeD}}_2^{+} \) are more likely to decompose into Be^0,+ + D₂ than into BeD^0,+ + D. The metastable BeD₃ and \( {\mathrm{BeD}}_3^{+} \) predominantly decompose into BeD^0,+ + D₂. In light of our results on the reaction energetics, we can interpret the pathways for production of BeD via BeD₂ and BeD₃ intermediates observed in molecular dynamics simulations.     

Interaction-induced electric (hyper)polarizability in the dihydrogen–neon pair: basis set and electron correlation effects

We investigated the interaction (hyper)polarizability of neon–dihydrogen pairs by performing high-level ab initio calculations with atom/molecule-specific, purpose-oriented Gaussian basis sets. We obtained interaction-induced electric properties at the SCF, MP2, and CCSD levels of theory. At the CCSD level, for the T-shaped configuration, around the respective potential minimum of 6.437 a₀, the interaction-induced mean first hyperpolarizability varies for 5?<? R/a₀?<?10 as

$$ \left[{\overline{\beta}}_{\mathrm{int}}(R)\hbox{-} {\overline{\beta}}_{\mathrm{int}}\left({R}_{\mathrm{e}}\right)\right]/{e}^3{a_0}^3{E_{\mathrm{h}}}^{-2}=-0.91\left(R\hbox{-} {R}_{\mathrm{e}}\right)+0.50{\left(R\hbox{-} {R}_{\mathrm{e}}\right)}^2\hbox{--} 0.13{\left(R\hbox{-} {R}_{\mathrm{e}}\right)}^3+0.01{\left(R\hbox{-} {R}_{\mathrm{e}}\right)}^4. $$

Again, at the CCSD level, but for the L-shaped configuration around the respective potential minimum of 6.572 a₀, this property varies for 5?<? R/a₀?<?10 as

$$ \left[{\overline{\beta}}_{\mathrm{int}}(R)\hbox{-} {\overline{\beta}}_{\mathrm{int}}\left({R}_{\mathrm{e}}\right)\right]/{e}^3{a_0}^3{E_{\mathrm{h}}}^{-2}=-1.33\left(R\hbox{-} {R}_{\mathrm{e}}\right)+0.75{\left(R\hbox{-} {R}_{\mathrm{e}}\right)}^2-0.20{\left(R\hbox{-} {R}_{\mathrm{e}}\right)}^3+0.02{\left(R\hbox{-} {R}_{\mathrm{e}}\right)}^4. $$

Open image in new window

Graphical Abstract Interaction-induced mean dipole polarizability (\( \overline{a} \)) for the T-shaped configuration of H₂–Ne calculated at the SCF, MP2, and CCSD levels of theory