Noncanonical function of MEKK2 and MEK5 PB1 domains for coordinated extracellular signal-regulated kinase 5 and c-Jun N-terminal kinase signaling

MEKK2 and MEK5 encode Phox/Bem1p (PB1) domains that heterodimerize with one another. MEKK2, MEK5, and extracellular signal-related kinase 5 (ERK5) form a ternary complex through interactions involving the MEKK2 and MEK5 PB1 domains and a 34-amino-acid C-terminal extension of the MEK5 PB1 domain. This C-terminal extension encodes an ERK5 docking site required for MEK5 activation of ERK5. The PB1 domains bind in a front-to-back arrangement, with a cluster of basic amino acids in the front of the MEKK2 PB1 domain binding to the back-end acidic clusters of the MEK5 PB1 domain. The C-terminal moiety, including the acidic cluster of the MEKK2 PB1 domain, is not required for MEK5 binding and binds MKK7. Quiescent MEKK2 preferentially binds MEK5, and MEKK2 activation results in ERK5 activation. Activated MEKK2 binds and activates MKK7, leading to JNK activation. The findings define how the MEKK2 and MEK5 PB1 domains are uniquely used for differential binding of two mitogen-activated protein kinase kinases, MEK5 and MKK7, for the coordinated control of ERK5 and c-Jun N-terminal kinase activation.     

Overexpression of proline oxidase induces proline-dependent and mitochondria-mediated apoptosis

The neuropathogenesis of influenza-associated encephalopathy in children and Reye's syndrome remains unclear. A surveillance effort conducted during 2000-2003 in South-West Japan reveals that almost all fatal and handicapped influenza-associated encephalopathy patients exhibit a disorder of mitochondrial β-oxidation with elevated serum acylcarnitine ratios (C_16:0+C_18:1)/C₂. Here we show invasion by a non-neurotropic epidemic influenza A H3N2 virus in cerebral capillaries with progressive brain edema after intranasal infection of mice having impaired mitochondrial β-oxidation congenitally or posteriorly in the newborn/ suckling periods. Mice genetically lacking of carnitine transporter OCTN2, resulting in carnitine deficiency and impaired β-oxidation, exhibited significant higher virus-genome numbers in the brain, accumulation of virus antigen exclusively in the cerebral capillaries and increased brain vascular permeability compared to in wild type mice. Mini-plasmin, which proteolytically potentiates influenza virus multiplication in vivo and destroys the blood-brain barrier, accumulated with virus antigen in the brain capillaries of OCTN2-deficient mice but only a little in wild-type mice. These results suggest that the impaired mitochondrial β-oxidation changes the susceptibility to a non-neurotropic influenza A virus as to multiplication in the brain capillaries and to cause brain edema. These pathological findings in the brain of mice having impaired mitochondrial β-oxidation after influenza virus infection may have implications for human influenza-associated encephalopathy.     

On the NMR analysis of pKa values in the unfolded state of proteins by extrapolation to zero denaturant

Detailed knowledge of the pH-dependence in both folded and unfolded states of proteins is essential to understand the role of electrostatics in protein stability. The increasing number of natively disordered proteins constitutes an excellent source for the NMR analysis of pKa values in the unfolded state of proteins. However, the tendency of many natively disordered proteins to aggregate via intermolecular hydrophobic clusters limits their NMR analysis over a wide pH range. To assess whether the pKa values in natively disordered polypeptides can be extrapolated from NMR measurements in the presence of denaturants, the natively disordered backbone of the C-terminal fragment 75 to 105 of Human Thioredoxin was studied. First, assignments using triple resonance experiments were performed to confirm lack of secondary structure. Then the pH-dependence of the amides and carboxylate side chains of Glu residues (Glu88, Glu95, Glu98, and Glu103) in the pH range from 2.0 to 7.0 was monitored using 2D 1H15N HSQC and 3D C(CO)NH experiments, and the behavior of their amides and corresponding carboxyl groups was compared to confirm the absence of nonlocal interactions. Lastly, the effect of increasing dimethyl urea concentration on the pKa values of these Glu residues was monitored. The results indicate that: (i) the dispersion in the pKa of carboxyl groups and the pH midpoints of amides in Glu residues is about 0.5 pH units and 0.6 pH units, respectively; (ii) the backbone amides of the Glu residues exhibit pH midpoints which are within 0.2 pH units from those of their carboxylates; (iii) the addition of denaturant produces upshifts in the pKa values of Glu residues that are nearly independent of their position in the sequence; and (iv) these upshifts show a nonlinear behavior in denaturant concentration, complicating the extrapolation to zero denaturant. Nevertheless, the relative ordering of the pKa values of Glu residues is preserved over the whole range of denaturant concentrations indicating that measurements at high denaturant concentration (e.g. 4 M dimethyl urea) can yield a qualitatively correct ranking of the pKa of these residues in natively disordered proteins whose pH-dependence cannot be monitored directly by NMR.     

DNA bar coding and pyrosequencing to identify rare HIV drug resistance mutations

Treatment of HIV-infected individuals with antiretroviral agents selects for drug-resistant mutants, resulting in frequent treatment failures. Although the major antiretroviral resistance mutations are routinely characterized by DNA sequencing, treatment failures are still common, probably in part because undetected rare resistance mutations facilitate viral escape. Here we combined DNA bar coding and massively parallel pyrosequencing to quantify rare drug resistance mutations. Using DNA bar coding, we were able to analyze seven viral populations in parallel, overall characterizing 118093 sequence reads of average length 103bp. Analysis of a control HIV mixture showed that resistance mutations present as 5% of the population could be readily detected without false positive calls. In three samples of multidrug-resistant HIV populations from patients, all the drug-resistant mutations called by conventional analysis were identified, as well as four additional low abundance drug resistance mutations, some of which would be expected to influence the response to antiretroviral therapy. Methods for sensitive characterization of HIV resistance alleles have been reported, but only the pyrosequencing method allows all the positions at risk for drug resistance mutations to be interrogated deeply for many HIV populations in a single experiment.     

A unified approach for quantifying, testing and correcting population stratification in case-control association studies

The HapMap project has given case-control association studies a unique opportunity to uncover the genetic basis of complex diseases. However, persistent issues in such studies remain the proper quantification of, testing for, and correction for population stratification (PS). In this paper, we present the first unified paradigm that addresses all three fundamental issues within one statistical framework. Our unified approach makes use of an omnibus quantity (delta), which can be estimated in a case-control study from suitable null loci. We show how this estimated value can be used to quantify PS, to statistically test for PS, and to correct for PS, all in the context of case-control studies. Moreover, we provide guidelines for interpreting values of delta in association studies (e.g., at alpha = 0.05, a delta of size 0.416 is small, a delta of size 0.653 is medium, and a delta of size 1.115 is large). A novel feature of our testing procedure is its ability to test for either strictly any PS or only 'practically important' PS. We also performed simulations to compare our correction procedure with Genomic Control (GC). Our results show that, unlike GC, it maintains good Type I error rates and power across all levels of PS.     

Predicting Secchi disk depth from average beam attenuation in a deep,ultra-clear lake

We addressed potential sources of error in estimating the water clarity of mountain lakes by investigating the use of beam transmissometer measurements to estimate Secchi disk depth. The optical properties Secchi disk depth (SD) and beam transmissometer attenuation (BA) were measured in Crater Lake (Crater Lake National Park, Oregon, USA) at a designated sampling station near the maximum depth of the lake. A standard 20 cm black and white disk was used to measure SD. The transmissometer light source had a nearly monochromatic wavelength of 660 nm and a path length of 25 cm. We created a SD prediction model by regression of the inverse SD of 13 measurements recorded on days when environmental conditions were acceptable for disk deployment with BA averaged over the same depth range as the measured SD. The relationship between inverse SD and averaged BA was significant and the average 95% confidence interval for predicted SD relative to the measured SD was ±1.6 m (range = −4.6 to 5.5 m) or ±5.0%. Eleven additional sample dates tested the accuracy of the predictive model. The average 95% confidence interval for these sample dates was ±0.7 m (range = −3.5 to 3.8 m) or ±2.2%. The 1996–2000 time-series means for measured and predicted SD varied by 0.1 m, and the medians varied by 0.5 m. The time-series mean annual measured and predicted SD’s also varied little, with intra-annual differences between measured and predicted mean annual SD ranging from −2.1 to 0.1 m. The results demonstrated that this prediction model reliably estimated Secchi disk depths and can be used to significantly expand optical observations in an environment where the conditions for standardized SD deployments are limited.     

Seasonal nutrient and plankton dynamics in a physical-biological model of Crater Lake

A coupled 1D physical-biological model of Crater Lake is presented. The model simulates the seasonal evolution of two functional phytoplankton groups, total chlorophyll, and zooplankton in good quantitative agreement with observations from a 10-year monitoring study. During the stratified period in summer and early fall the model displays a marked vertical structure: the phytoplankton biomass of the functional group 1, which represents diatoms and dinoflagellates, has its highest concentration in the upper 40 m; the phytoplankton biomass of group 2, which represents chlorophyta, chrysophyta, cryptomonads and cyanobacteria, has its highest concentrations between 50 and 80 m, and phytoplankton chlorophyll has its maximum at 120 m depth. A similar vertical structure is a reoccurring feature in the available data. In the model the key process allowing a vertical separation between biomass and chlorophyll is photoacclimation. Vertical light attenuation (i.e., water clarity) and the physiological ability of phytoplankton to increase their cellular chlorophyll-to-biomass ratio are ultimately determining the location of the chlorophyll maximum. The location of the particle maxima on the other hand is determined by the balance between growth and losses and occurs where growth and losses equal. The vertical particle flux simulated by our model agrees well with flux measurements from a sediment trap. This motivated us to revisit a previously published study by Dymond et al. (1996). Dymond et al. used a box model to estimate the vertical particle flux and found a discrepancy by a factor 2.5–10 between their model-derived flux and measured fluxes from a sediment trap. Their box model neglected the exchange flux of dissolved and suspended organic matter, which, as our model and available data suggests is significant for the vertical exchange of nitrogen. Adjustment of Dymond et al.’s assumptions to account for dissolved and suspended nitrogen yields a flux estimate that is consistent with sediment trap measurements and our model.     

Germ-line DNA copy number variation frequencies in a large North American population

Genomic copy number variation (CNV) is a recently identified form of global genetic variation in the human genome. The Affymetrix GeneChip 100 and 500 K SNP genotyping platforms were used to perform a large-scale population-based study of CNV frequency. We constructed a genomic map of 578 CNV regions, covering approximately 220 Mb (7.3%) of the human genome, identifying 183 previously unknown intervals. Copy number changes were observed to occur infrequently (<1%) in the majority (>93%) of these genomic regions, but encompass hundreds of genes and disease loci. This North American population-based map will be a useful resource for future genetic studies. Electronic supplementary material The online version of this article (doi:) contains supplementary material, which is available to authorized users.