Contrast-FEL—A Test for Differences in Selective Pressures at Individual Sites among Clades and Sets of Branches

A number of evolutionary hypotheses can be tested by comparing selective pressures among sets of branches in a phylogenetic tree. When the question of interest is to identify specific sites within genes that may be evolving differently, a common approach is to perform separate analyses on subsets of sequences and compare parameter estimates in a post hoc fashion. This approach is statistically suboptimal and not always applicable. Here, we develop a simple extension of a popular fixed effects likelihood method in the context of codon-based evolutionary phylogenetic maximum likelihood testing, Contrast-FEL. It is suitable for identifying individual alignment sites where any among the K≥2 sets of branches in a phylogenetic tree have detectably different ω ratios, indicative of different selective regimes. Using extensive simulations, we show that Contrast-FEL delivers good power, exceeding 90% for sufficiently large differences, while maintaining tight control over false positive rates, when the model is correctly specified. We conclude by applying Contrast-FEL to data from five previously published studies spanning a diverse range of organisms and focusing on different evolutionary questions.     

Neto2 Interacts with the Scaffolding Protein GRIP and Regulates Synaptic Abundance of Kainate Receptors

Kainate receptors (KARs) are a class of ionotropic glutamate receptors that are expressed throughout the central nervous system. The function and subcellular localization of KARs are tightly regulated by accessory proteins. We have previously identified the single-pass transmembrane proteins, Neto1 and Neto2, to be associated with native KARs. In the hippocampus, Neto1, but not Neto2, controls the abundance and modulates the kinetics of postsynaptic KARs. Here we evaluated whether Neto2 regulates synaptic KAR levels in the cerebellum where Neto1 expression is limited to the deep cerebellar nuclei. In the cerebellum, where Neto2 is present abundantly, we found a ∼40% decrease in GluK2-KARs at the postsynaptic density (PSD) of Neto2-null mice. No change, however, was observed in total level of GluK2-KARs, thereby suggesting a critical role of Neto2 on the synaptic localization of cerebellar KARs. The presence of a putative class II PDZ binding motif on Neto2 led us to also investigate whether it interacts with PDZ domain-containing proteins previously implicated in regulating synaptic abundance of KARs. We identified a PDZ-dependent interaction between Neto2 and the scaffolding protein GRIP. Furthermore, coexpression of Neto2 significantly increased the amount of GRIP associated with GluK2, suggesting that Neto2 may promote and/or stabilize GluK2:GRIP interactions. Our results demonstrate that Neto2, like Neto1, is an important auxiliary protein for modulating the synaptic levels of KARs. Moreover, we propose that the interactions of Neto1/2 with various scaffolding proteins is a critical mechanism by which KARs are stabilized at diverse synapses.     

Solution Structure of the HIV-1 Intron Splicing Silencer and Its Interactions with the UP1 Domain of Heterogeneous Nuclear Ribonucleoprotein (hnRNP) A1

Splicing patterns in human immunodeficiency virus type 1 (HIV-1) are maintained through cis regulatory elements that recruit antagonistic host RNA-binding proteins. The activity of the 3′ acceptor site A7 is tightly regulated through a complex network of an intronic splicing silencer (ISS), a bipartite exonic splicing silencer (ESS3a/b), and an exonic splicing enhancer (ESE3). Because HIV-1 splicing depends on protein-RNA interactions, it is important to know the tertiary structures surrounding the splice sites. Herein, we present the NMR solution structure of the phylogenetically conserved ISS stem loop. ISS adopts a stable structure consisting of conserved UG wobble pairs, a folded 2X2 (GU/UA) internal loop, a UU bulge, and a flexible AGUGA apical loop. Calorimetric and biochemical titrations indicate that the UP1 domain of heterogeneous nuclear ribonucleoprotein A1 binds the ISS apical loop site-specifically and with nanomolar affinity. Collectively, this work provides additional insights into how HIV-1 uses a conserved RNA structure to commandeer a host RNA-binding protein.     

β-Amyloid (Aβ) Oligomers Impair Brain-derived Neurotrophic Factor Retrograde Trafficking by Down-regulating Ubiquitin C-terminal Hydrolase,UCH-L1

We previously found that BDNF-dependent retrograde trafficking is impaired in AD transgenic mouse neurons. Utilizing a novel microfluidic culture chamber, we demonstrate that Aβ oligomers compromise BDNF-mediated retrograde transport by impairing endosomal vesicle velocities, resulting in impaired downstream signaling driven by BDNF/TrkB, including ERK5 activation, and CREB-dependent gene regulation. Our data suggest that a key mechanism mediating the deficit involves ubiquitin C-terminal hydrolase L1 (UCH-L1), a deubiquitinating enzyme that functions to regulate cellular ubiquitin. Aβ-induced deficits in BDNF trafficking and signaling are mimicked by LDN (an inhibitor of UCH-L1) and can be reversed by increasing cellular UCH-L1 levels, demonstrated here using a transducible TAT-UCH-L1 strategy. Finally, our data reveal that UCH-L1 mRNA levels are decreased in the hippocampi of AD brains. Taken together, our data implicate that UCH-L1 is important for regulating neurotrophin receptor sorting to signaling endosomes and supporting retrograde transport. Further, our results support the idea that in AD, Aβ may down-regulate UCH-L1 in the AD brain, which in turn impairs BDNF/TrkB-mediated retrograde signaling, compromising synaptic plasticity and neuronal survival.     

Quantification of the importance of individual steps in the control of aromatic amino acid metabolism.

The quantitative importance of the individual steps of aromatic amino acid metabolism in rat liver was determined by calculation of the respective Control Coefficients (Strengths). The Control Coefficient of tryptophan 2,3-dioxygenase for tryptophan degradation was determined in a variety of physiological conditions and with a range of activities of tryptophan 2,3-dioxygenase. The Control Coefficient varied from 0.75 with basal enzyme activity to 0.25 after maximal induction of the enzyme by dexamethasone. The remainder of the control for tryptophan degradation was associated with the transport of the amino acid across the plasma membrane, with only very small contributions from kynureninase and kynurenine hydroxylase. The Control Coefficients of tyrosine aminotransferase for tyrosine degradation were approx. 0.70 and 0.20 with basal and dexamethasone-induced tyrosine aminotransferase activities respectively; the Control Coefficients of the transport of the amino acid into the cell were 0.22 and 0.58 respectively. Phenylalanine hydroxylase was found to have a Control Coefficient for the degradation of phenylalanine of approx. 0.50 under conditions of basal enzyme activity; after maximal activation by glucagon, the Control Coefficient decreased to 0.12. The transport of phenylalanine was responsible for the remaining control in the pathway. These results have important implications, directly for the regulation of aromatic amino acid metabolism in the liver, and indirectly for the regulation of neuroamine synthesis in the brain.     

Binding of a calcium sensitizer, bepridil, to cardiac troponin C. A fluorescence stopped-flow kinetic, circular dichroism, and proton nuclear magnetic resonance study

Stopped-flow fluorescence kinetic measurements, circular dichroism (CD), and 1H nuclear magnetic resonance (NMR) spectroscopy at 360 MHz have been used to study the interaction of the calcium-channel blocker and calmodulin antagonist bepridil with cardiac troponin C (cTnC) in the presence of calcium. The kinetic data show that bepridil reduces the rate of calcium release only from the low affinity, calcium-specific site and not from the two high affinity calcium/magnesium sites. CD measurements indicate that drug binding leads to a small increase in the alpha-helical content of the complex. 1H NMR shows that the protein binds one equivalent of bepridil, with a dissociation constant of approximately 20 microM, only when the low affinity calcium site is occupied. Exchange is fast or intermediate on the chemical shift time scale. Drug binding is shown to be largely localized in the N-terminal domain, containing the low affinity calcium site, by observing the shifting and broadening of several resonances associated with that domain. These include assigned aromatic signals together with methionyl and other methyl signals. Observation of intermolecular nuclear Overhauser effects was precluded by extensive spectral overlap. Consideration of the data from the three techniques permitted a model of the bepridil-cTnC complex to be constructed, using the model of cTnC derived from the x-ray structure of calmodulin (MacLachlan L. K., Reid, D. G., and Carter, N. (1990) J. Biol. Chem. 265, 9754-9763). Binding of bepridil to a prominent hydrophobic depression in the N-terminal domain can be invoked to explain many of the induced changes in the spectral and kinetic properties of the protein. The implications of the model for the calcium sensitizing action of bepridil are discussed.     

HLA-B51 and HLA-Bw52 differ by only two amino acids which are in the helical region of the alpha 1 domain

Genes encoding the serologically cross-reactive HLA-B51 and HLA-Bw52 molecules were isolated and the exons sequenced. HLA-B51 genes obtained from Caucasian and Oriental individuals were identical. HLA-Bw52 differs from HLA-B51 by four nucleotide substitutions in exon 2 encoding the alpha 1 domain. These comprise one isolated silent substitution in codon 23 and a cluster of three coding substitutions in codons 63 and 67. Amino acid substitutions of N----E at position 63 and F----S at position 67 are the only differences between HLA-B51 and HLA-Bw52 and these residues are postulated to form HLA-B51 specific epitopes. HLA-B51 could have been formed from HLA-Bw52 by the combination of a genetic exchange with HLA-B8 and a point mutation. Similarity of HLA-B51 and HLA-Bw52 with HLA-Bw58 suggest they also share a common ancestor.     

Intracellular transport of class I HLA molecules is affected by polymorphic residues in the binding groove

Intracellular transport of class I MHC complexes is dependent on assembly of class I heavy chains with ₂-microglobulin (₂m) and peptides. This suggests that amino acid residues of individual class I molecules which are important for their stability and transport are likely to include those which contribute to binding of a majority of the cleft-associated peptides. To identify such critical residues, substitutions at polymorphic positions within the peptide binding cleft were introduced into a mutant HLA-A*0201 molecule bearing an additional gly>lys substitution at position 242 (242K). The 242K mutation weakens association of the HLA-A*201 heavy chain with ₂m and was used to enhance potential effects of substitutions in the peptide binding groove on class I stability. Critical in choosing which binding cleft positions to mutate was the observation that HLA-A*6801 was less sensitive to the effects of 242K mutation than HLA-A*0201 and A*6901. This suggested that one or more of the six residues in the 2 domain differing between HLA-A*6901 and A*6801 were likely to affect class I complex stability. Positions 95, 97, 107, 114, 116, and 156 in either 242K or wild-type HLA-A*0201 molecules were therefore each converted to those residues found in HLA-A*6801. One of the second-site substitutions, arg>met at position 97, increased stability and restored surface expression of the 242K molecule. Five other substitutions either had no additional effect or further impaired 242K stability. Substitution of his>arg at position 114 blocked surface expression of both 242K and wild-type HLA-A*0201 molecules. These results demonstrate that polymorphic residues in the binding cleft influence the stability of class I complexes, and suggest that position 97 plays a critical role in stabilizing class I molecules for transport.