Proteomics evaluation of five economical commercial abundant protein depletion kits for enrichment of diseases-specific biomarkers from blood serum

Blood serum is arguably the most analyzed biofluid for disease prediction and diagnosis. Herein, we benchmarked five different serum abundant protein depletion (SAPD) kits with regard to the identification of disease-specific biomarkers in human serum using bottom-up proteomics. As expected, the IgG removal efficiency among the SAPD kits is highly variable, ranging from 70% to 93%. A pairwise comparison of database search results showed a 10%–19% variation in protein identification among the kits. Immunocapturing-based SAPD kits against IgG and albumin outperformed the others in the removal of these two abundant proteins. Conversely, non-antibody-based methods (i.e., kits using ion exchange resins) and kits leveraging a multi-antibody approach were proven to be less efficient in depleting IgG/albumin from samples but led to the highest number of identified peptides. Notably, our results indicate that different cancer biomarkers could be enriched up to 10% depending on the utilized SAPD kit compared with the undepleted sample. Additionally, functional analysis of the bottom-up proteomic results revealed that different SAPD kits enrich distinct disease- and pathway-specific protein sets. Overall, our study emphasizes that a careful selection of the appropriate commercial SAPD kit is crucial for the analysis of disease biomarkers in serum by shotgun proteomics.     

Breast Cancer Biology and Ethnic Disparities in Breast Cancer Mortality in New Zealand: A Cohort Study

Introduction

Indigenous Māori women have a 60% higher breast cancer mortality rate compared with European women in New Zealand. We investigated differences in cancer biological characteristics and their impact on breast cancer mortality disparity between Māori and NZ European women.

Materials and Methods

Data on 2849 women with primary invasive breast cancers diagnosed between 1999 and 2012 were extracted from the Waikato Breast Cancer Register. Differences in distribution of cancer biological characteristics between Māori and NZ European women were explored adjusting for age and socioeconomic deprivation in logistic regression models. Impacts of socioeconomic deprivation, stage and cancer biological characteristics on breast cancer mortality disparity between Māori and NZ European women were explored in Cox regression models.

Results

Compared with NZ European women (n=2304), Māori women (n=429) had significantly higher rates of advanced and higher grade cancers. Māori women also had non-significantly higher rates of ER/PR negative and HER-2 positive breast cancers. Higher odds of advanced stage and higher grade remained significant for Māori after adjusting for age and deprivation. Māori women had almost a 100% higher age and deprivation adjusted breast cancer mortality hazard compared with NZ European women (HR=1.98, 1.55-2.54). Advanced stage and lower proportion of screen detected cancer in Māori explained a greater portion of the excess breast cancer mortality (HR reduction from 1.98 to 1.38), while the additional contribution through biological differences were minimal (HR reduction from 1.38 to 1.35).

Conclusions

More advanced cancer stage at diagnosis has the greatest impact while differences in biological characteristics appear to be a minor contributor for inequities in breast cancer mortality between Māori and NZ European women. Strategies aimed at reducing breast cancer mortality in Māori should focus on earlier diagnosis, which will likely have a greater impact on reducing breast cancer mortality inequity between Māori and NZ European women.     

A substrate-specific cytochrome P450 monooxygenase, CYP6AB11, from the polyphagous navel orangeworm (Amyelois transitella)

The navel orangeworm Amyelois transitella (Walker) (Lepidoptera: Pyralidae) is a serious pest of many tree crops in California orchards, including almonds, pistachios, walnuts and figs. To understand the molecular mechanisms underlying detoxification of phytochemicals, insecticides and mycotoxins by this species, full-length CYP6AB11 cDNA was isolated from larval midguts using RACE PCR. Phylogenetic analysis of this insect cytochrome P450 monooxygenase established its evolutionary relationship to a P450 that selectively metabolizes imperatorin (a linear furanocoumarin) and myristicin (a natural methylenedioxyphenyl compound) in another lepidopteran species. Metabolic assays conducted with baculovirus-expressed P450 protein, P450 reductase and cytochrome b₅ on 16 compounds, including phytochemicals, mycotoxins, and synthetic pesticides, indicated that CYP6AB11 efficiently metabolizes imperatorin (0.88 pmol/min/pmol P450) and slowly metabolizes piperonyl butoxide (0.11 pmol/min/pmol P450). LC-MS analysis indicated that the imperatorin metabolite is an epoxide generated by oxidation of the double bond in its extended isoprenyl side chain. Predictive structures for CYP6AB11 suggested that its catalytic site contains a doughnut-like constriction over the heme that excludes aromatic rings on substrates and allows only their extended side chains to access the catalytic site. CYP6AB11 can also metabolize the principal insecticide synergist piperonyl butoxide (PBO), a synthetic methylenedioxyphenyl compound, albeit slowly, which raises the possibility that resistance may evolve in this species after exposure to synergists under field conditions.     

Homology modeling of plant cytochrome P450s

Plant P450 monooxygenases represent the largest family of plant proteins and the largest collection of P450s available for comparative studies and biotechnological applications. They have been shown to catalyze a diverse array of difficult chemical reactions and have demonstrated potential to be used in pharmacological, agronomic and phytoremediative applications. Central to our use of these catalytically competent enzymes is the need to understand their interactions with substrates. Because most characterized plant P450s are membrane-bound proteins that are resistant to standard X-ray and NMR structure determinations, homology modeling represents a reliable and relatively rapid alternative method for analyzing structure-function relationships and predicting substrates for many P450s that are only now being characterized. These methods, which are being widely used in mammalian P450 structure-function studies, can allow plant biologists to define critical residues interacting with substrates and, in a directed fashion, alter the reactivities of individual monooxygenases. The homology modelings that have been done on a limited number of plant P450s and the site-directed mutations that validate them indicate that current modeling and substrate docking procedures are capable of providing structural explanations for sequence variants as well as for predicting functional characteristics of undefined P450s.     

Epac and phospholipase Cepsilon regulate Ca2+ release in the heart by activation of protein kinase Cepsilon and calcium-calmodulin kinase II

Recently, we identified a novel signaling pathway involving Epac, Rap, and phospholipase C (PLC)epsilon that plays a critical role in maximal beta-adrenergic receptor (betaAR) stimulation of Ca2+-induced Ca2+ release (CICR) in cardiac myocytes. Here we demonstrate that PLCepsilon phosphatidylinositol 4,5-bisphosphate hydrolytic activity and PLCepsilon-stimulated Rap1 GEF activity are both required for PLCepsilon-mediated enhancement of sarcoplasmic reticulum Ca2+ release and that PLCepsilon significantly enhances Rap activation in response to betaAR stimulation in the heart. Downstream of PLCepsilon hydrolytic activity, pharmacological inhibition of PKC significantly inhibited both betaAR- and Epac-stimulated increases in CICR in PLCepsilon+/+ myocytes but had no effect in PLCepsilon-/- myocytes. betaAR and Epac activation caused membrane translocation of PKCepsilon in PLCepsilon+/+ but not PLCepsilon-/- myocytes and small interfering RNA-mediated PKCepsilon knockdown significantly inhibited both betaAR and Epac-mediated CICR enhancement. Further downstream, the Ca2+/calmodulin-dependent protein kinase II (CamKII) inhibitor, KN93, inhibited betaAR- and Epac-mediated CICR in PLCepsilon+/+ but not PLCepsilon-/- myocytes. Epac activation increased CamKII Thr286 phosphorylation and enhanced phosphorylation at CamKII phosphorylation sites on the ryanodine receptor (RyR2) (Ser2815) and phospholamban (Thr17) in a PKC-dependent manner. Perforated patch clamp experiments revealed that basal and betaAR-stimulated peak L-type current density are similar in PLCepsilon+/+ and PLCepsilon-/- myocytes suggesting that control of sarcoplasmic reticulum Ca2+ release, rather than Ca2+ influx through L-type Ca2+ channels, is the target of regulation of a novel signal transduction pathway involving sequential activation of Epac, PLCepsilon, PKCepsilon, and CamKII downstream of betaAR activation.     

EQ-5D-3L Derived Population Norms for Health Related Quality of Life in Sri Lanka

Background

Health Related Quality of Life (HRQoL) is an important outcome measure in health economic evaluation that guides health resource allocations. Population norms for HRQoL are an essential ingredient in health economics and in the evaluation of population health. The aim of this study was to produce EQ-5D-3L-derived population norms for Sri Lanka.

Method

A population sample (n =  780) was selected from four districts of Sri Lanka. A stratified cluster sampling approach with probability proportionate to size was employed. Twenty six clusters of 30 participants each were selected; each participant completed the EQ-5D-3L in a face-to-face interview. Utility weights for their EQ-5D-3L health states were assigned using the Sri Lankan EQ-5D-3L algorithm. The population norms are reported by age and socio-economic variables.

Results

The EQ-5D-3L was completed by 736 people, representing a 94% response rate. Sixty per cent of the sample reported being in full health. The percentage of people responding to any problems in the five EQ-5D-3L dimensions increased with age. The mean EQ-5D-3L weight was 0.85 (SD 0.008; 95%CI 0.84-0.87). The mean EQ-5D-3L weight was significantly associated with age, housing type, disease experience and religiosity. People above 70 years of age were 7.5 times more likely to report mobility problems and 3.7 times more likely to report pain/discomfort than those aged 18-29 years. Those with a tertiary education were five times less likely to report any HRQoL problems than those without a tertiary education. A person living in a shanty was 4.3 more likely to have problems in usual activities than a person living in a single house.

Conclusion

The population norms in Sri Lanka vary with socio-demographic characteristics. The socioeconomically disadvantaged have a lower HRQoL. The trends of population norms observed in this lower middle income country were generally similar to those previously reported in high income countries.     

VITAL NMR: using chemical shift derived secondary structure information for a limited set of amino acids to assess homology model accuracy

Homology modeling is a powerful tool for predicting protein structures, whose success depends on obtaining a reasonable alignment between a given structural template and the protein sequence being analyzed. In order to leverage greater predictive power for proteins with few structural templates, we have developed a method to rank homology models based upon their compliance to secondary structure derived from experimental solid-state NMR (SSNMR) data. Such data is obtainable in a rapid manner by simple SSNMR experiments (e.g., ¹³C–¹³C 2D correlation spectra). To test our homology model scoring procedure for various amino acid labeling schemes, we generated a library of 7,474 homology models for 22 protein targets culled from the TALOS+/SPARTA+ training set of protein structures. Using subsets of amino acids that are plausibly assigned by SSNMR, we discovered that pairs of the residues Val, Ile, Thr, Ala and Leu (VITAL) emulate an ideal dataset where all residues are site specifically assigned. Scoring the models with a predicted VITAL site-specific dataset and calculating secondary structure with the Chemical Shift Index resulted in a Pearson correlation coefficient (−0.75) commensurate to the control (−0.77), where secondary structure was scored site specifically for all amino acids (ALL 20) using STRIDE. This method promises to accelerate structure procurement by SSNMR for proteins with unknown folds through guiding the selection of remotely homologous protein templates and assessing model quality.     

Identification of cysteines involved in S-nitrosylation, S-glutathionylation, and oxidation to disulfides in ryanodine receptor type 1

The skeletal muscle Ca(2+)-release channel (ryanodine receptor type 1 (RyR1)) is a redox sensor, susceptible to reversible S-nitrosylation, S-glutathionylation, and disulfide oxidation. So far, Cys-3635 remains the only cysteine residue identified as functionally relevant to the redox sensing properties of the channel. We demonstrate that expression of the C3635A-RyR1 mutant in RyR1-null myotubes alters the sensitivity of the ryanodine receptor to activation by voltage, indicating that Cys-3635 is involved in voltage-gated excitation-contraction coupling. However, H(2)O(2) treatment of C3635A-RyR1 channels or wild-type RyR1, following their expression in human embryonic kidney cells, enhances [(3)H]ryanodine binding to the same extent, suggesting that cysteines other than Cys-3635 are responsible for the oxidative enhancement of channel activity. Using a combination of Western blotting and sulfhydryl-directed fluorescent labeling, we found that two large regions of RyR1 (amino acids 1-2401 and 3120-4475), previously shown to be involved in disulfide bond formation, are also major sites of both S-nitrosylation and S-glutathionylation. Using selective isotopecoded affinity tag labeling of RyR1 and matrix-assisted laser desorption/ionization time-of-flight mass spectroscopy, we identified, out of the 100 cysteines in each RyR1 subunit, 9 that are endogenously modified (Cys-36, Cys-315, Cys-811, Cys-906, Cys-1591, Cys-2326, Cys-2363, Cys-3193, and Cys-3635) and another 3 residues that were only modified with exogenous redox agents (Cys-253, Cys-1040, and Cys-1303). We also identified the types of redox modification each of these cysteines can undergo. In summary, we have identified a discrete subset of cysteines that are likely to be involved in the functional response of RyR1 to different redox modifications (S-nitrosylation, S-glutathionylation, and oxidation to disulfides).     

Common active site architecture and binding strategy of four phenylpropanoid P450s from Arabidopsis thaliana as revealed by molecular modeling

Despite extensive primary sequence diversity, crystal structures of several bacterial cytochrome P450 monooxygenases (P450s) and a single eukaryotic P450 indicate that these enzymes share a structural core of alpha-helices and beta-sheets and vary in the loop regions contacting individual substrates. To determine the extent to which individual structural features are conserved among divergent P450s existing in a single biosynthetic pathway, we have modeled the structures of four highly divergent P450s (CYP73A5, CYP84A1, CYP75B1, CYP98A3) in the Arabidopsis phenylpropanoid pathway synthesizing lignins, flavonoids and anthocyanins. Analysis of these models has indicated that, despite primary sequence identities as low as 13%, the structural cores and several loop regions of these P450s are highly conserved. Substrate docking indicated that all four enzymes employ a common strategy to identify their substrates in that their cinnamate-derived substrates align along helix I with their aromatic ring positioned towards the C-terminus of this helix and their aliphatic tails positioned towards the N-terminus. Further similarity was observed in the way the substrates contact the consensus P450 substrate recognition sites (SRS). Residues predicted to contact the aromatic ring region exist in SRS5, SRS6 and the C-terminal portion of SRS4 and residues contacting the distal end of each substrate exist in SRS1, SRS2 and the N-terminal portion of SRS4. Alignments of the regions contacting the aromatic ring region indicate that SRS4, SRS5 and SRS6 share higher degrees of sequence conservation than found in SRS1, SRS2 or the full-length protein.