Identification of potent EGFR inhibitors from TCM Database@Taiwan

Overexpression of epidermal growth factor receptor (EGFR) has been associated with cancer. Targeted inhibition of the EGFR pathway has been shown to limit proliferation of cancerous cells. Hence, we employed Traditional Chinese Medicine Database (TCM nawiaT@esabataD) (http://tcm.cmu.edu.tw) to identify potential EGFR inhibitor. Multiple Linear Regression (MLR), Support Vector Machine (SVM), Comparative Molecular Field Analysis (CoMFA), and Comparative Molecular Similarities Indices Analysis (CoMSIA) models were generated using a training set of EGFR ligands of known inhibitory activities. The top four TCM candidates based on DockScore were 2-O-caffeoyl tartaric acid, Emitine, Rosmaricine, and 2-O-feruloyl tartaric acid, and all had higher binding affinities than the control Iressa®. The TCM candidates had interactions with Asp855, Lys716, and Lys728, all which are residues of the protein kinase binding site. Validated MLR (r² = 0.7858) and SVM (r² = 0.8754) models predicted good bioactivity for the TCM candidates. In addition, the TCM candidates contoured well to the 3D-Quantitative Structure-Activity Relationship (3D-QSAR) map derived from the CoMFA (q² = 0.721, r² = 0.986) and CoMSIA (q² = 0.662, r² = 0.988) models. The steric field, hydrophobic field, and H-bond of the 3D-QSAR map were well matched by each TCM candidate. Molecular docking indicated that all TCM candidates formed H-bonds within the EGFR protein kinase domain. Based on the different structures, H-bonds were formed at either Asp855 or Lys716/Lys728. The compounds remained stable throughout molecular dynamics (MD) simulation. Based on the results of this study, 2-O-caffeoyl tartaric acid, Emitine, Rosmaricine, and 2-O-feruloyl tartaric acid are suggested to be potential EGFR inhibitors.     

Predicting total, abdominal, visceral and hepatic adiposity with circulating biomarkers in caucasian and Japanese american women

Background

Characterization of abdominal and intra-abdominal fat requires imaging, and thus is not feasible in large epidemiologic studies.

Objective

We investigated whether biomarkers may complement anthropometry (body mass index [BMI], waist circumference [WC], and waist-hip ratio [WHR]) in predicting the size of the body fat compartments by analyzing blood biomarkers, including adipocytokines, insulin resistance markers, sex steroid hormones, lipids, liver enzymes and gastro-neuropeptides.

Methods

Fasting levels of 58 blood markers were analyzed in 60 healthy, Caucasian or Japanese American postmenopausal women who underwent anthropometric measurements, dual energy X-ray absorptiometry (DXA), and abdominal magnetic resonance imaging. Total, abdominal, visceral and hepatic adiposity were predicted based on anthropometry and the biomarkers using Random Forest models.

Results

Total body fat was well predicted by anthropometry alone (R² = 0.85), by the 5 best predictors from the biomarker model alone (leptin, leptin-adiponectin ratio [LAR], free estradiol, plasminogen activator inhibitor-1 [PAI1], alanine transaminase [ALT]; R² = 0.69), or by combining these 5 biomarkers with anthropometry (R² = 0.91). Abdominal adiposity (DXA trunk-to-periphery fat ratio) was better predicted by combining the two types of predictors (R² = 0.58) than by anthropometry alone (R² = 0.53) or the 5 best biomarkers alone (25(OH)-vitamin D₃, insulin-like growth factor binding protein-1 [IGFBP1], uric acid, soluble leptin receptor [sLEPR], Coenzyme Q10; R² = 0.35). Similarly, visceral fat was slightly better predicted by combining the predictors (R² = 0.68) than by anthropometry alone (R² = 0.65) or the 5 best biomarker predictors alone (leptin, C-reactive protein [CRP], LAR, lycopene, vitamin D₃; R² = 0.58). Percent liver fat was predicted better by the 5 best biomarker predictors (insulin, sex hormone binding globulin [SHBG], LAR, alpha-tocopherol, PAI1; R² = 0.42) or by combining the predictors (R² = 0.44) than by anthropometry alone (R² = 0.29).

Conclusion

The predictive ability of anthropometry for body fat distribution may be enhanced by measuring a small number of biomarkers. Studies to replicate these data in men and other ethnic groups are warranted.     

Allelic variants of melanocortin 3 receptor gene (MC3R) and weight loss in obesity: a randomised trial of hypo-energetic high- versus low-fat diets

Introduction

The melanocortin system plays an important role in energy homeostasis. Mice genetically deficient in the melanocortin-3 receptor gene have a normal body weight with increased body fat, mild hypophagia compared to wild-type mice. In humans, Thr6Lys and Val81Ile variants of the melanocortin-3 receptor gene (MC3R) have been associated with childhood obesity, higher BMI Z-score and elevated body fat percentage compared to non-carriers. The aim of this study is to assess the association in adults between allelic variants of MC3R with weight loss induced by energy-restricted diets.

Subjects and Methods

This research is based on the NUGENOB study, a trial conducted to assess weight loss during a 10-week dietary intervention involving two different hypo-energetic (high-fat and low-fat) diets. A total of 760 obese patients were genotyped for 10 single nucleotide polymorphisms covering the single exon of MC3R gene and its flanking regions, including the missense variants Thr6Lys and Val81Ile. Linear mixed models and haplotype-based analysis were carried out to assess the potential association between genetic polymorphisms and differential weight loss, fat mass loss, waist change and resting energy expenditure changes.

Results

No differences in drop-out rate were found by MC3R genotypes. The rs6014646 polymorphism was significantly associated with weight loss using co-dominant (p = 0.04) and dominant models (p = 0.03). These p-values were not statistically significant after strict control for multiple testing. Haplotype-based multivariate analysis using permutations showed that rs3827103–rs1543873 (p = 0.06), rs6014646–rs6024730 (p = 0.05) and rs3746619–rs3827103 (p = 0.10) displayed near-statistical significant results in relation to weight loss. No other significant associations or gene*diet interactions were detected for weight loss, fat mass loss, waist change and resting energy expenditure changes.

Conclusion

The study provided overall sufficient evidence to support that there is no major effect of genetic variants of MC3R and differential weight loss after a 10-week dietary intervention with hypo-energetic diets in obese Europeans.     

Prolactin secretion in healthy adults is determined by gender, age and body mass index

Background

Prolactin (PRL) secretion is quantifiable as mean, peak and nadir PRL concentrations, degree of irregularity (ApEn, approximate entropy) and spikiness (brief staccato-like fluctuations).

Hypothesis

Distinct PRL dynamics reflect relatively distinct (combinations of) subject variables, such as gender, age, and BMI.

Location

Clinical Research Unit.

Subjects

Seventy-four healthy adults aged 22–77 yr (41 women and 33 men), with BMI 18.3–39.4 kg/m².

Measures

Immunofluorometric PRL assay of 10-min samples collected for 24 hours.

Results

Mean 24-h PRL concentration correlated jointly with gender (P<0.0001) and BMI (P = 0.01), but not with age (overall R² = 0.308, P<0.0001). Nadir PRL concentration correlated with gender only (P = 0.017) and peak PRL with gender (P<0.001) and negatively with age (P<0.003), overall R² = 0.325, P<0.0001. Forward-selection multivariate regression of PRL deconvolution results demonstrated that basal (nonpulsatile) PRL secretion tended to be associated with BMI (R² = 0.058, P = 0.03), pulsatile secretion with gender (R² = 0.152, P = 0.003), and total secretion with gender and BMI (R² = 0.204, P<0.0001). Pulse mass was associated with gender (P = 0.001) and with a negative tendency to age (P = 0.038). In male subjects older than 50 yr (but not in women) approximate entropy was increased (0.942±0.301 vs. 1.258±0.267, P = 0.007) compared with younger men, as well as spikiness (0.363±0.122 vs. 0463±2.12, P = 0.031). Cosinor analysis disclosed higher mesor and amplitude in females than in men, but the acrophase was gender-independent. The acrophase was determined by age and BMI (R² = 0.186, P = 0.001).

Conclusion

In healthy adults, selective combinations of gender, age, and BMI specify distinct PRL dynamics, thus requiring balanced representation of these variables in comparative PRL studies.     

Lung function is associated with arterial stiffness in children

Background

In older adults, an independent association exists between impaired lung function and cardiovascular disease. This interaction might be related to the effects of aging and/or smoking. In order to explore possible childhood antecedents to this association, we hypothesized that decreased lung function and vascular stiffness might be related, in early life.

Objective

To determine the relationship between lung function and carotid augmentation index (AIx), a measure of vascular stiffness, in 8-year old children.

Methods

Data on brachial blood pressure, lung function (FEV₁, FVC, FEV₁/FVC, obtained by spirometry) and carotid AIx75 (AIx standardised to an arbitrary heart rate of 75 beats per minute, obtained by applanation tonometry) was available in 249 community-based 8-year old children. These healthy children had been subjects in a randomised controlled trial of two interventions (omega-3 fatty acid supplementation and house-dust mite avoidance) to prevent asthma. Smoking in pregnancy and childhood environmental tobacco smoke (ETS) exposure was prospectively collected by questionnaire. The association between lung function and carotid AIx75 was assessed in multivariate models that included sex, height, smoking status during pregnancy, ETS exposure and randomisation groups (house dust mite avoidance and dietary intervention) as covariates.

Results

In the fully adjusted models, Carotid AIx75 was independently associated with FEV1 (standardised β = −0.17,b = −6.72, partial R² = .02, p = 0.03), FVC (standardised β = −0.29, b = −9.31, partial R² = 0.04, p<0.001) and FEV1/FVC (standardised β = .13, b = 18.4, partial R² = 0.02, p = 0.04).

Conclusion

Lower lung volumes are associated with increased vascular stiffness at an early age. The interaction between lung function and vascular stiffness may thus represent more than just age-related alterations in both the pulmonary and vascular systems.     

Design and development of topoisomerase inhibitors using molecular modelling studies

Topoisomerase inhibitors are used as anticancer and antibacterial agents. A series of novel 2,4,6-tri-substituted pyridine derivatives reported as topoisomerase inhibitors were used for quantitative structure–activity relationship (QSAR) study. In order to understand the structural requirement of these topoisomerase inhibitors, a ligand-based pharmacophore and atom-based 3D-QSAR model have been developed. A five-point pharmacophore with one hydrophobic group (H4), four aromatic rings (R5, R6, R7 and R8) was obtained. The pharmacophore hypothesis yielded a 3D-QSAR model with good partial least-square (PLS) statistic results. The training set correlation is characterized by PLS factors (r² = 0.7892, SD = 0.2948, F = 49.9, P = 1.379). The test set correlation is characterized by PLS factors (q² = 0.7776, root mean squared error = 0.2764, Pearson R = 0.8926). The docking study revealed the binding orientations of these inhibitors at active site amino acid residues of topoisomerases enzyme. The results of pharmacophore hypothesis and 3D-QSAR provided the detail structural insights as well as highlighted the important binding features of novel 2,4,6-tri-substituted pyridine derivatives and can be developed as potent topoisomerase inhibitors.

Figure

Open in a separate windowKey structural requirement for topoisomerase activity     

Enrichment of variations in KIR3DL1/S1 and KIR2DL2/L3 among H1N1/09 ICU patients: an exploratory study

Background

Infection by the pandemic influenza A (H1N1/09) virus resulted in significant pathology among specific ethnic groups worldwide. Natural Killer (NK) cells are important in early innate immune responses to viral infections. Activation of NK cells, in part, depend on killer-cell immunoglobulin-like receptors (KIR) and HLA class I ligand interactions. To study factors involved in NK cell dysfunction in overactive immune responses to H1N1 infection, KIR3DL1/S1 and KIR2DL2/L3 allotypes and cognate HLA ligands of H1N1/09 intensive-care unit (ICU) patients were determined.

Methodology and Findings

KIR3DL1/S1, KIR2DL2/L3, and HLA -B and -C of 51 H1N1/09 ICU patients and 105 H1N1-negative subjects (St. Theresa Point, Manitoba) were characterized. We detected an increase of 3DL1 ligand-negative pairs (3DL1/S1⁺ Bw6⁺ Bw4⁻), and a lack of 2DL1 HLA-C2 ligands, among ICU patients. They were also significantly enriched for 2DL2/L3 ligand-positive pairs (P<0.001, Pc<0.001; Odds Ratio:6.3158, CI95%:2.481–16.078). Relative to St. Theresa aboriginals (STh) and Venezuelan Amerindians (VA), allotypes enriched among aboriginal ICU patients (Ab) were: 2DL3 (Ab>VA, P = 0.024, Pc = 0.047; Odds Ratio:2.563, CI95%:1.109–5.923), 3DL1*00101 (Ab>VA, P<0.001, Pc<0.001), 3DL1*01502 (Ab>STh, P = 0.034, Pc = 0.268), and 3DL1*029 (Ab>STh, P  = 0.039, Pc = 0.301). Aboriginal patients ligand-positive for 3DL1/S1 and 2DL1 had the lowest probabilities of death (R_d) (R_d = 28%), compared to patients that were 3DL1/S1 ligand-negative (R_d = 52%) or carried 3DL1*029 (R_d = 52%). Relative to Caucasoids (CA), two allotypes were enriched among non-aboriginal ICU patients (NAb): 3DL1*00401 (NAb>CA, P<0.001, Pc<0.001) and 3DL1*01502 (CAP = 0.012, Pc = 0.156). Non-aboriginal patients with ligands for all three KIRs (3DL1/S1, 2DL2/L3, and 2DL1) had the lowest probabilities of death (R_d = 36%), compared to subjects with 3DL1*01502 (R_d = 48%) and/or 3DL1*00401 (R_d = 58%).

Conclusions

Specific KIR3DL1/S1 allotypes, 3DL1/S1 and 2DL1 ligand-negative pairs, and 2DL2/L3 ligand-positive pairs were enriched among ICU patients. This suggests a possible association with NK cell dysfunction in patients with overactive immune responses to H1N1/09, leading to severe disease.     

Memory enhancement by traditional Chinese medicine?

Abstract

Cognitive repair by insulin-like growth factor-I (IGF-I) through activation of insulin-like growth factor-I receptor (IGF-IR) is well established, but not used for clinical therapy due to its link to cancer. We hypothesize that IGF-IR activation rather than IGF-I per se may be essential for cognitive repair and attempted to identify ligands from traditional Chinese medicine (TCM) with drug-like potential towards IGF-IR. TCM ligands, 3-(2-carboxyphenyl)-4(3H)-quinazolinone from Isatisin digotica, (+)-N-methyllaurotetanine from Lindera aggregate, and (+)-1(R)-Coclaurine from Nelumbonucifera Gaertn, exhibited high binding affinities and good blood brain barrier (BBB) penetration crucial for accessing IGF-IR. Stable complex formation of the candidates was observed during molecular dynamics (MD) simulation. Interactions with Leu975 and Gly1055 or Asp1056 were important for ligand binding. Amino acid distance analysis revealed residues 974/975, 984–986, 996–1006, 1040–1056, and 1122–1135 as “hotspots” for ligand binding in IGF-IR. Versatile entry pathways for the TCM candidates suggest high accessibility to the binding site. Blockage of the binding site opening by the TCM candidates limits binding site access by other compounds. Multiple linear regression (R²?=?0.9715), support vector machine (R²?=?0.9084), Bayesian network (R²?=?0.8233) comparative molecular field analysis (CoMFA, R²?=?0.9941), and comparative molecular similarity indices analysis (CoMSIA, R²?=?0.9877) models consistently suggest that the TCM candidates might exert bioactivity on IGF-IR. Contour of representative MD conformations to CoMFA and CoMSIA maps exhibits similar results. Properties including BBB passage, evidence of ability to form stable complexes with IGF-IR by MD simulation, and predicted bioactivity suggest that the TCM candidates have drug-like properties and might have potential as cognitive-enhancing drugs.

An animated interactive 3D complement (I3DC) is available in Proteopedia at http://proteopedia.org/w/Journal:JBSD:38     

The obesity-associated polymorphisms FTO rs9939609 and MC4R rs17782313 and endometrial cancer risk in non-Hispanic white women

Overweight and obesity are strongly associated with endometrial cancer. Several independent genome-wide association studies recently identified two common polymorphisms, FTO rs9939609 and MC4R rs17782313, that are linked to increased body weight and obesity. We examined the association of FTO rs9939609 and MC4R rs17782313 with endometrial cancer risk in a pooled analysis of nine case-control studies within the Epidemiology of Endometrial Cancer Consortium (E2C2). This analysis included 3601 non-Hispanic white women with histologically-confirmed endometrial carcinoma and 5275 frequency-matched controls. Unconditional logistic regression models were used to assess the relation of FTO rs9939609 and MC4R rs17782313 genotypes to the risk of endometrial cancer. Among control women, both the FTO rs9939609 A and MC4R rs17782313 C alleles were associated with a 16% increased risk of being overweight (p = 0.001 and p = 0.004, respectively). In case-control analyses, carriers of the FTO rs9939609 AA genotype were at increased risk of endometrial carcinoma compared to women with the TT genotype [odds ratio (OR)  = 1.17; 95% confidence interval (CI): 1.03–1.32, p = 0.01]. However, this association was no longer apparent after adjusting for body mass index (BMI), suggesting mediation of the gene-disease effect through body weight. The MC4R rs17782313 polymorphism was not related to endometrial cancer risk (per allele OR = 0.98; 95% CI: 0.91–1.06; p = 0.68). FTO rs9939609 is a susceptibility marker for white non-Hispanic women at higher risk of endometrial cancer. Although FTO rs9939609 alone might have limited clinical or public health significance for identifying women at high risk for endometrial cancer beyond that of excess body weight, further investigation of obesity-related genetic markers might help to identify the pathways that influence endometrial carcinogenesis.     

Modulating Uranium Binding Affinity in Engineered Calmodulin EF-Hand Peptides: Effect of Phosphorylation

To improve our understanding of uranium toxicity, the determinants of uranyl affinity in proteins must be better characterized. In this work, we analyzed the contribution of a phosphoryl group on uranium binding affinity in a protein binding site, using the site 1 EF-hand motif of calmodulin. The recombinant domain 1 of calmodulin from A. thaliana was engineered to impair metal binding at site 2 and was used as a structured template. Threonine at position 9 of the loop was phosphorylated in vitro, using the recombinant catalytic subunit of protein kinase CK2. Hence, the T₉TKE₁₂ sequence was substituted by the CK2 recognition sequence TAAE. A tyrosine was introduced at position 7, so that uranyl and calcium binding affinities could be determined by following tyrosine fluorescence. Phosphorylation was characterized by ESI-MS spectrometry, and the phosphorylated peptide was purified to homogeneity using ion-exchange chromatography. The binding constants for uranyl were determined by competition experiments with iminodiacetate. At pH 6, phosphorylation increased the affinity for uranyl by a factor of ∼5, from K_d = 25±6 nM to K_d = 5±1 nM. The phosphorylated peptide exhibited a much larger affinity at pH 7, with a dissociation constant in the subnanomolar range (K_d = 0.25±0.06 nM). FTIR analyses showed that the phosphothreonine side chain is partly protonated at pH 6, while it is fully deprotonated at pH 7. Moreover, formation of the uranyl-peptide complex at pH 7 resulted in significant frequency shifts of the ν_as(P-O) and ν_s(P-O) IR modes of phosphothreonine, supporting its direct interaction with uranyl. Accordingly, a bathochromic shift in ν_as(UO₂)²⁺ vibration (from 923 cm⁻¹ to 908 cm⁻¹) was observed upon uranyl coordination to the phosphorylated peptide. Together, our data demonstrate that the phosphoryl group plays a determining role in uranyl binding affinity to proteins at physiological pH.     

No differential regulation of dopamine transporter (DAT) and vesicular monoamine transporter 2 (VMAT2) binding in a primate model of Parkinson disease

Radioligands for DAT and VMAT2 are widely used presynaptic markers for assessing dopamine (DA) nerve terminals in Parkinson disease (PD). Previous in vivo imaging and postmortem studies suggest that these transporter sites may be regulated as the numbers of nigrostriatal neurons change in pathologic conditions. To investigate this issue, we used in vitro quantitative autoradioradiography to measure striatal DAT and VMAT2 specific binding in postmortem brain from 14 monkeys after unilateral internal carotid artery infusion of 1-Methyl-4-Phenyl-1,2,3,6-tetrahydropyridine (MPTP) with doses varying from 0 to 0.31 mg/kg. Quantitative estimates of the number of tyrosine hydroxylase (TH)-immunoreactive (ir) neurons in substantia nigra (SN) were determined with unbiased stereology, and quantitative autoradiography was used to measure DAT and VMAT2 striatal specific binding. Striatal VMAT2 and DAT binding correlated with striatal DA (r_s = 0.83, r_s = 0.80, respectively, both with n = 14, p<0.001) but only with nigra TH-ir cells when nigral cell loss was 50% or less (r = 0.93, n = 8, p = 0.001 and r = 0.91, n = 8, p = 0.002 respectively). Reduction of VMAT2 and DAT striatal specific binding sites strongly correlated with each other (r = 0.93, n = 14, p<0.0005). These similar changes in DAT and VMAT2 binding sites in the striatal terminal fields of the surviving nigrostriatal neurons demonstrate that there is no differential regulation of these two sites at 2 months after MPTP infusion.     

Prediction of promiscuous p-glycoprotein inhibition using a novel machine learning scheme

Background

P-glycoprotein (P-gp) is an ATP-dependent membrane transporter that plays a pivotal role in eliminating xenobiotics by active extrusion of xenobiotics from the cell. Multidrug resistance (MDR) is highly associated with the over-expression of P-gp by cells, resulting in increased efflux of chemotherapeutical agents and reduction of intracellular drug accumulation. It is of clinical importance to develop a P-gp inhibition predictive model in the process of drug discovery and development.

Methodology/Principal Findings

An in silico model was derived to predict the inhibition of P-gp using the newly invented pharmacophore ensemble/support vector machine (PhE/SVM) scheme based on the data compiled from the literature. The predictions by the PhE/SVM model were found to be in good agreement with the observed values for those structurally diverse molecules in the training set (n = 31, r ² = 0.89, q ² = 0.86, RMSE = 0.40, s = 0.28), the test set (n = 88, r ² = 0.87, RMSE = 0.39, s = 0.25) and the outlier set (n = 11, r ² = 0.96, RMSE = 0.10, s = 0.05). The generated PhE/SVM model also showed high accuracy when subjected to those validation criteria generally adopted to gauge the predictivity of a theoretical model.

Conclusions/Significance

This accurate, fast and robust PhE/SVM model that can take into account the promiscuous nature of P-gp can be applied to predict the P-gp inhibition of structurally diverse compounds that otherwise cannot be done by any other methods in a high-throughput fashion to facilitate drug discovery and development by designing drug candidates with better metabolism profile.     

Genetic Modulation of Lipid Profiles following Lifestyle Modification or Metformin Treatment: The Diabetes Prevention Program

Weight-loss interventions generally improve lipid profiles and reduce cardiovascular disease risk, but effects are variable and may depend on genetic factors. We performed a genetic association analysis of data from 2,993 participants in the Diabetes Prevention Program to test the hypotheses that a genetic risk score (GRS) based on deleterious alleles at 32 lipid-associated single-nucleotide polymorphisms modifies the effects of lifestyle and/or metformin interventions on lipid levels and nuclear magnetic resonance (NMR) lipoprotein subfraction size and number. Twenty-three loci previously associated with fasting LDL-C, HDL-C, or triglycerides replicated (P = 0.04–1×10⁻¹⁷). Except for total HDL particles (r = −0.03, P = 0.26), all components of the lipid profile correlated with the GRS (partial |r| = 0.07–0.17, P = 5×10⁻⁵–1×10⁻¹⁹). The GRS was associated with higher baseline-adjusted 1-year LDL cholesterol levels (β = +0.87, SEE±0.22 mg/dl/allele, P = 8×10⁻⁵, P _interaction = 0.02) in the lifestyle intervention group, but not in the placebo (β = +0.20, SEE±0.22 mg/dl/allele, P = 0.35) or metformin (β = −0.03, SEE±0.22 mg/dl/allele, P = 0.90; P _interaction = 0.64) groups. Similarly, a higher GRS predicted a greater number of baseline-adjusted small LDL particles at 1 year in the lifestyle intervention arm (β = +0.30, SEE±0.012 ln nmol/L/allele, P = 0.01, P _interaction = 0.01) but not in the placebo (β = −0.002, SEE±0.008 ln nmol/L/allele, P = 0.74) or metformin (β = +0.013, SEE±0.008 nmol/L/allele, P = 0.12; P _interaction = 0.24) groups. Our findings suggest that a high genetic burden confers an adverse lipid profile and predicts attenuated response in LDL-C levels and small LDL particle number to dietary and physical activity interventions aimed at weight loss.     

Allergic asthmatics show divergent lipid mediator profiles from healthy controls both at baseline and following birch pollen provocation

Background

Asthma is a respiratory tract disorder characterized by airway hyper-reactivity and chronic inflammation. Allergic asthma is associated with the production of allergen-specific IgE and expansion of allergen-specific T-cell populations. Progression of allergic inflammation is driven by T-helper type 2 (Th2) mediators and is associated with alterations in the levels of lipid mediators.

Objectives

Responses of the respiratory system to birch allergen provocation in allergic asthmatics were investigated. Eicosanoids and other oxylipins were quantified in the bronchoalveolar lumen to provide a measure of shifts in lipid mediators associated with allergen challenge in allergic asthmatics.

Methods

Eighty-seven lipid mediators representing the cyclooxygenase (COX), lipoxygenase (LOX) and cytochrome P450 (CYP) metabolic pathways were screened via LC-MS/MS following off-line extraction of bronchoalveolar lavage fluid (BALF). Multivariate statistics using OPLS were employed to interrogate acquired oxylipin data in combination with immunological markers.

Results

Thirty-two oxylipins were quantified, with baseline asthmatics possessing a different oxylipin profile relative to healthy individuals that became more distinct following allergen provocation. The most prominent differences included 15-LOX-derived ω-3 and ω-6 oxylipins. Shared-and-Unique-Structures (SUS)-plot modeling showed a correlation (R² = 0.7) between OPLS models for baseline asthmatics (R²Y[cum] = 0.87, Q²[cum] = 0.51) and allergen-provoked asthmatics (R²Y[cum] = 0.95, Q²[cum] = 0.73), with the majority of quantified lipid mediators and cytokines contributing equally to both groups. Unique structures for allergen provocation included leukotrienes (LTB₄ and 6-trans-LTB₄), CYP-derivatives of linoleic acid (epoxides/diols), and IL-10.

Conclusions

Differences in asthmatic relative to healthy profiles suggest a role for 15-LOX products of both ω-6 and ω-3 origin in allergic inflammation. Prominent differences at baseline levels indicate that non-symptomatic asthmatics are subject to an underlying inflammatory condition not observed with other traditional mediators. Results suggest that oxylipin profiling may provide a sensitive means of characterizing low-level inflammation and that even individuals with mild disease display distinct phenotypic profiles, which may have clinical ramifications for disease.     

A measure of the promiscuity of proteins and characteristics of residues in the vicinity of the catalytic site that regulate promiscuity

Promiscuity, the basis for the evolution of new functions through ‘tinkering’ of residues in the vicinity of the catalytic site, is yet to be quantitatively defined. We present a computational method Promiscuity Indices Estimator (PROMISE) - based on signatures derived from the spatial and electrostatic properties of the catalytic residues, to estimate the promiscuity (PromIndex) of proteins with known active site residues and 3D structure. PromIndex reflects the number of different active site signatures that have congruent matches in close proximity of its native catalytic site, the quality of the matches and difference in the enzymatic activity. Promiscuity in proteins is observed to follow a lognormal distribution (μ = 0.28, σ = 1.1 reduced chi-square = 3.0E-5). The PROMISE predicted promiscuous functions in any protein can serve as the starting point for directed evolution experiments. PROMISE ranks carboxypeptidase A and ribonuclease A amongst the more promiscuous proteins. We have also investigated the properties of the residues in the vicinity of the catalytic site that regulates its promiscuity. Linear regression establishes a weak correlation (R²∼0.1) between certain properties of the residues (charge, polar, etc) in the neighborhood of the catalytic residues and PromIndex. A stronger relationship states that most proteins with high promiscuity have high percentages of charged and polar residues within a radius of 3 Å of the catalytic site, which is validated using one-tailed hypothesis tests (P-values∼0.05). Since it is known that these characteristics are key factors in catalysis, their relationship with the promiscuity index cross validates the methodology of PROMISE.     

Caloric restriction induces changes in insulin and body weight measurements that are inversely associated with subsequent weight regain

Background

Successful weight maintenance following weight loss is challenging for many people. Identifying predictors of longer-term success will help target clinical resources more effectively. To date, focus has been predominantly on the identification of predictors of weight loss. The goal of the current study was to determine if changes in anthropometric and clinical parameters during acute weight loss are associated with subsequent weight regain.

Methodology

The study consisted of an 8-week low calorie diet (LCD) followed by a 6-month weight maintenance phase. Anthropometric and clinical parameters were analyzed before and after the LCD in the 285 participants (112 men, 173 women) who regained weight during the weight maintenance phase. Mixed model ANOVA, Spearman correlation, and linear regression were used to study the relationships between clinical measurements and weight regain.

Principal Findings

Gender differences were observed for body weight and several clinical parameters at both baseline and during the LCD-induced weight loss phase. LCD-induced changes in BMI (Spearman’s ρ = 0.22, p = 0.0002) were inversely associated with weight regain in both men and women. LCD-induced changes in fasting insulin (ρ = 0.18, p = 0.0043) and HOMA-IR (ρ = 0.19, p = 0.0023) were also associated independently with weight regain in both genders. The aforementioned associations remained statistically significant in regression models taking account of variables known to independently influence body weight.

Conclusions/Significance

LCD-induced changes in BMI, fasting insulin, and HOMA-IR are inversely associated with weight regain in the 6-month period following weight loss.     

Association of variants at 1q32 and STAT3 with ankylosing spondylitis suggests genetic overlap with Crohn's disease

Ankylosing spondylitis (AS) is a common inflammatory arthritic condition. Overt inflammatory bowel disease (IBD) occurs in about 10% of AS patients, and in addition 70% of AS cases may have subclinical terminal ileitis. Spondyloarthritis is also common in IBD patients. We therefore tested Crohn''s disease susceptibility genes for association with AS, aiming to identify pleiotropic genetic associations with both diseases. Genotyping was carried out using Sequenom and Applied Biosystems TaqMan and OpenArray technologies on 53 markers selected from 30 Crohn''s disease associated genomic regions. We tested genotypes in a population of unrelated individual cases (n = 2,773) and controls (n = 2,215) of white European ancestry for association with AS. Statistical analysis was carried out using a Cochran-Armitage test for trend in PLINK. Strong association was detected at chr1q32 near KIF21B (rs11584383, P = 1.6×10⁻¹⁰, odds ratio (OR) = 0.74, 95% CI:0.68–0.82). Association with disease was also detected for 2 variants within STAT3 (rs6503695, P = 4.6×10⁻⁴. OR = 0.86 (95% CI:0.79–0.93); rs744166, P = 2.6×10⁻⁵, OR = 0.84 (95% CI:0.77–0.91)). Association was confirmed for IL23R (rs11465804, P = 1.2×10⁻⁵, OR = 0.65 (95% CI:0.54–0.79)), and further associations were detected for IL12B (rs10045431, P = 5.2×10⁻⁵, OR = 0.83 (95% CI:0.76–0.91)), CDKAL1 (rs6908425, P = 1.1×10⁻⁴, OR = 0.82 (95% CI:0.74–0.91)), LRRK2/MUC19 (rs11175593, P = 9.9×10⁻⁵, OR = 1.92 (95% CI: 1.38–2.67)), and chr13q14 (rs3764147, P = 5.9×10⁻⁴, OR = 1.19 (95% CI: 1.08–1.31)). Excluding cases with clinical IBD did not significantly affect these findings. This study identifies chr1q32 and STAT3 as ankylosing spondylitis susceptibility loci. It also further confirms association for IL23R and detects suggestive association with another 4 loci. STAT3 is a key signaling molecule within the Th17 lymphocyte differentiation pathway and further enhances the case for a major role of this T-lymphocyte subset in ankylosing spondylitis. Finally these findings suggest common aetiopathogenic pathways for AS and Crohn''s disease and further highlight the involvement of common risk variants across multiple diseases.     

Functional Characterization Improves Associations between Rare Non-Synonymous Variants in CHRNB4 and Smoking Behavior

Smoking is the leading cause of preventable death worldwide. Accordingly, effort has been devoted to determining the genetic variants that contribute to smoking risk. Genome-wide association studies have identified several variants in nicotinic acetylcholine receptor genes that contribute to nicotine dependence risk. We previously undertook pooled sequencing of the coding regions and flanking sequence of the CHRNA5, CHRNA3, CHRNB4, CHRNA6 and CHRNB3 genes and found that rare missense variants at conserved residues in CHRNB4 are associated with reduced risk of nicotine dependence among African Americans. We identified 10 low frequency (<5%) non-synonymous variants in CHRNB4 and investigated functional effects by co-expression with normal α3 or α4 subunits in human embryonic kidney cells. Voltage-clamp was used to obtain acetylcholine and nicotine concentration–response curves and qRT-PCR, western blots and cell-surface ELISAs were performed to assess expression levels. These results were used to functionally weight genetic variants in a gene-based association test. We find that there is a highly significant correlation between carrier status weighted by either acetylcholine EC₅₀ (β = −0.67, r² = 0.017, P = 2×10⁻⁴) or by response to low nicotine (β = −0.29, r² = 0.02, P = 6×10⁻⁵) when variants are expressed with the α3 subunit. In contrast, there is no significant association when carrier status is unweighted (β = −0.04, r² = 0.0009, P = 0.54). These results highlight the value of functional analysis of variants and the advantages to integrating such data into genetic studies. They also suggest that an increased sensitivity to low concentrations of nicotine is protective from the risk of developing nicotine dependence.