Rapid protein domain assignment from amino acid sequence using predicted secondary structure

The elucidation of the domain content of a given protein sequence in the absence of determined structure or significant sequence homology to known domains is an important problem in structural biology. Here we address how successfully the delineation of continuous domains can be accomplished in the absence of sequence homology using simple baseline methods, an existing prediction algorithm (Domain Guess by Size), and a newly developed method (DomSSEA). The study was undertaken with a view to measuring the usefulness of these prediction methods in terms of their application to fully automatic domain assignment. Thus, the sensitivity of each domain assignment method was measured by calculating the number of correctly assigned top scoring predictions. We have implemented a new continuous domain identification method using the alignment of predicted secondary structures of target sequences against observed secondary structures of chains with known domain boundaries as assigned by Class Architecture Topology Homology (CATH). Taking top predictions only, the success rate of the method in correctly assigning domain number to the representative chain set is 73.3%. The top prediction for domain number and location of domain boundaries was correct for 24% of the multidomain set (±20 residues). These results have been put into context in relation to the results obtained from the other prediction methods assessed

The four major N- and C-terminal splice variants of the excitatory amino acid transporter GLT-1 form cell surface homomeric and heteromeric assemblies

The L-glutamate transporter GLT-1 is an abundant CNS membrane protein of the excitatory amino acid transporter (EAAT) family which controls extracellular L-glutamate levels and is important in limiting excitotoxic neuronal death. Using RT-PCR, we have determined that four mRNAs encoding GLT-1 exist in mouse brain, with the potential to encode four GLT-1 isoforms that differ in their N- and C-termini. We expressed all four isoforms (termed MAST-KREK, MPK-KREK, MAST-DIETCI and MPK-DIETCI according to amino acid sequence) in a range of cell lines and primary astrocytes and show that each isoform can reach the cell surface. In transfected HEK-293 or COS-7 cells, all four isoforms support high-affinity sodium-dependent L-glutamate uptake with identical pharmacological and kinetic properties. Inserting a viral epitope (V5, HA or FLAG) into the second extracellular domain of each isoform allowed co-immunoprecipitation and tr-FRET studies using transfected HEK-293 cells. Here we show for the first time that each of the four isoforms are able to combine to form homomeric and heteromeric assemblies, each of which are expressed at the cell surface of primary astrocytes. After activation of protein kinase C by phorbol ester, V5-tagged GLT-1 is rapidly removed from the cell surface of HEK-293 cells and degraded. This study provides direct biochemical evidence for oligomeric assembly of GLT-1 and reports the development of novel tools to provide insight into the trafficking of GLT-1.

IntFOLD: an integrated server for modelling protein structures and functions from amino acid sequences

IntFOLD is an independent web server that integrates our leading methods for structure and function prediction. The server provides a simple unified interface that aims to make complex protein modelling data more accessible to life scientists. The server web interface is designed to be intuitive and integrates a complex set of quantitative data, so that 3D modelling results can be viewed on a single page and interpreted by non-expert modellers at a glance. The only required input to the server is an amino acid sequence for the target protein. Here we describe major performance and user interface updates to the server, which comprises an integrated pipeline of methods for: tertiary structure prediction, global and local 3D model quality assessment, disorder prediction, structural domain prediction, function prediction and modelling of protein-ligand interactions. The server has been independently validated during numerous CASP (Critical Assessment of Techniques for Protein Structure Prediction) experiments, as well as being continuously evaluated by the CAMEO (Continuous Automated Model Evaluation) project. The IntFOLD server is available at: http://www.reading.ac.uk/bioinf/IntFOLD/

Critical role for an acidic amino acid region in platelet signaling by the HemITAM (hemi-immunoreceptor tyrosine-based activation motif) containing receptor CLEC-2 (C-type lectin receptor-2)

CLEC-2 is a member of new family of C-type lectin receptors characterized by a cytosolic YXXL downstream of three acidic amino acids in a sequence known as a hemITAM (hemi-immunoreceptor tyrosine-based activation motif). Dimerization of two phosphorylated CLEC-2 molecules leads to recruitment of the tyrosine kinase Syk via its tandem SH2 domains and initiation of a downstream signaling cascade. Using Syk-deficient and Zap-70-deficient cell lines we show that hemITAM signaling is restricted to Syk and that the upstream triacidic amino acid sequence is required for signaling. Using surface plasmon resonance and phosphorylation studies, we demonstrate that the triacidic amino acids are required for phosphorylation of the YXXL. These results further emphasize the distinct nature of the proximal events in signaling by hemITAM relative to ITAM receptors.

A hierarchical Bayesian model for predicting the functional consequences of amino-acid polymorphisms

Genetic polymorphisms in deoxyribonucleic acid coding regions may have a phenotypic effect on the carrier, e.g. by influencing susceptibility to disease. Detection of deleterious mutations via association studies is hampered by the large number of candidate sites; therefore methods are needed to narrow down the search to the most promising sites. For this, a possible approach is to use structural and sequence-based information of the encoded protein to predict whether a mutation at a particular site is likely to disrupt the functionality of the protein itself. We propose a hierarchical Bayesian multivariate adaptive regression spline (BMARS) model for supervised learning in this context and assess its predictive performance by using data from mutagenesis experiments on lac repressor and lysozyme proteins. In these experiments, about 12 amino-acid substitutions were performed at each native amino-acid position and the effect on protein functionality was assessed. The training data thus consist of repeated observations at each position, which the hierarchical framework is needed to account for. The model is trained on the lac repressor data and tested on the lysozyme mutations and vice versa. In particular, we show that the hierarchical BMARS model, by allowing for the clustered nature of the data, yields lower out-of-sample misclassification rates compared with both a BMARS and a frequen-tist MARS model, a support vector machine classifier and an optimally pruned classification tree.

Single amino acid substitutions in puroindoline-b mutants influence lipid binding properties

External reflectance Fourier transform infrared (ER-FTIR) spectroscopy and surface pressure measurements have been used to characterize the interaction of wild-type puroindoline-b (Pin-b) and two mutant forms featuring single residue substitutions-namely, Gly-46 to Ser-46 (Pin-bH) and Trp-44 to Arg-44 (Pin-bS)-with condensed-phase monolayers of zwitterionic (L-alpha-dipalmitoylphosphatidylcholine, DPPC) and anionic (L-alpha-dipalmitoylphosphatidyl-dl-glycerol, DPPG) phospholipids. The interaction with anionic DPPG monolayers, monitored by surface pressure isotherms, was influenced significantly by mutations in Pin-b (p < 0.05); wild-type Pin-b showed the highest surface pressure change of 10.6 +/- 1.0 mN m(-1), followed by Pin-bH (7.9 +/- 1.6 mN m(-1)) and Pin-bS (6.3 +/- 1.0 mN m(-1)), and the surface pressure isotherm kinetics were also different in each case. Integrated Amide I peak areas from corresponding ER-FTIR spectra confirmed the differences in adsorption kinetics, but also showed that differences in adsorbed amount were less significant, suggesting that mutations influence the degree of penetration into DPPG films. All Pin-b types showed evidence of interaction with DPPC films, detected as changes in surface pressure (5.6 +/- 1.1 mN m(-1)); however, no protein peaks were detected in the ER-FTIR spectra, which indicated that the interaction was via penetration with limited adsorption at the lipid/water interface. The expression of Pin-b mutants is linked to wheat endosperm hardness; therefore, the data presented here suggest that the lipid binding properties may be pivotal within the mechanism for this quality trait. In addition, the data suggest antimicrobial activities of Pin-b mutants would be lower than those of the wild-type Pin-b, because of decreased selectivity toward anionic phospholipids.

A model of net amino acid absorption and utilization by the portal-drained viscera of the lactating dairy cow

A more complete understanding of amino acid ( AA) metabolism by the various tissues of the body is required to improve upon current systems for predicting the use of absorbed AA. The objective of this work was to construct and parameterize a model of net removal of AA by the portal-drained viscera (PDV). Six cows were prepared with arterial, portal, and hepatic catheters and infused abomasally with 0, 200, 400, or 600 g of casein daily. Casein infusion increased milk yield quadratically and tended to increase milk protein yield quadratically. Arterial concentrations of a number of essential AA increased linearly with respect to infusion amount. When infused casein was assumed to have a true digestion coefficient of 0.95, the minimum likely true digestion coefficient for noninfused duodenal protein was found to be 0.80. Net PDV use of AA appeared to be linearly related to total supply (arterial plus absorption), and extraction percentages ranged from 0.5 to 7.25% for essential AA. Prediction errors for portal vein AA concentrations ranged from 4 to 9% of the observed mean concentrations. Removal of AA by PDV represented approximately 33% of total postabsorptive catabolic use, including use during absorption but excluding use for milk protein synthesis, and was apparently adequate to support endogenous N losses in feces of 18.4 g/d. As 69% of this use was from arterial blood, increased PDV catabolism of AA in part represents increased absorption of AA in excess of amounts required by other body tissues. Based on the present model, increased anabolic use of AA in the mammary and other tissues would reduce the catabolic use of AA by the PDV.

An integrative model of amino acid metabolism in the liver of the lactating dairy cow

The objective of this work was to construct a dynamic model of hepatic amino acid metabolism in the lactating dairy cow that could be parameterized using net flow data from in vivo experiments. The model considers 22 amino acids, ammonia, urea, and 13 energetic metabolites, and was parameterized using a steady-state balance model and two in vivo, net flow experiments conducted with mid-lactation dairy cows. Extracellular flows were derived directly from the observed data. An optimization routine was used to derive nine intracellular flows. The resulting dynamic model was found to be stable across a range of inputs suggesting that it can be perturbed and applied to other physiological states. Although nitrogen was generally in balance, leucine was in slight deficit compared to predicted needs for export protein synthesis, suggesting that an alternative source of leucine (e.g. peptides) was utilized. Simulations of varying glucagon concentrations indicated that an additional 5 mol/d of glucose could be synthesized at the reference substrate concentrations and blood flows. The increased glucose production was supported by increased removal from blood of lactate, glutamate, aspartate, alanine, asparagine, and glutamine. As glucose Output increased, ketone body and acetate release increased while CO2 release declined. The pattern of amino acids appearing in hepatic vein blood was affected by changes in amino acid concentration in portal vein blood, portal blood flow rate and glucagon concentration, with methionine and phenylalanine being the most affected of essential amino acids. Experimental evidence is insufficient to determine whether essential amino acids are affected by varying gluconeogenic demands. (C) 2004 Published by Elsevier Ltd.

A comparative evaluation of functions for partitioning nitrogen and amino acid intake between maintenance and growth in broilers

The results from three types of study with broilers, namely nitrogen (N) balance, bioassays and growth experiments, provided the data used herein. Sets of data on N balance and protein accretion (bioassay studies) were used to assess the ability of the monomolecular equation to describe the relationship between (i) N balance and amino acid (AA) intake and (ii) protein accretion and AA intake. The model estimated the levels of isoleucine, lysine, valine, threonine, methionine, total sulphur AAs and tryptophan resulting in zero balance to be 58, 59, 80, 96, 23, 85 and 32 mg/kg live weight (LW)/day, respectively. These estimates show good agreement with those obtained in previous studies. For the growth experiments, four models, specifically re-parameterized for analysing energy balance data, were evaluated for their ability to determine crude protein (CP) intake at maintenance and efficiency of utilization of CP intake for producing gain. They were: a straight line, two equations representing diminishing returns behaviour (monomolecular and rectangular hyperbola) and one equation describing smooth sigmoidal behaviour with a fixed point of inflexion (Gompertz). The estimates of CP requirement for maintenance and efficiency of utilization of CP intake for producing gain varied from 5.4 to 5.9 g/kg LW/day and 0.60 to 0.76, respectively, depending on the models.

Effect of prolonged intravenous glucose and essential amino acid infusion on nitrogen balance, muscle protein degradation and ubiquitin-conjugating enzyme gene expression in calves

Background: Intravenous infusions of glucose and amino acids increase both nitrogen balance and muscle accretion. We hypothesised that co-infusion of glucose ( to stimulate insulin) and essential amino acids (EAA) would act additively to improve nitrogen balance by decreasing muscle protein degradation in association with alterations in muscle expression of components of the ubiquitin-proteasome proteolytic pathway. Methods: We examined the effect of a 5 day intravenous infusions of saline, glucose, EAA and glucose + EAA, on urinary nitrogen excretion and muscle protein degradation. We carried out the study in 6 restrained calves since ruminants offer the advantage that muscle protein degradation can be assessed by excretion of 3 methyl-histidine and multiple muscle biopsies can be taken from the same animal. On the final day of infusion blood samples were taken for hormone and metabolite measurement and muscle biopsies for expression of ubiquitin, the 14-kDa E2 ubiquitin conjugating enzyme, and proteasome sub-units C2 and C8. Results: On day 5 of glucose infusion, plasma glucose, insulin and IGF-1 concentrations were increased while urea nitrogen excretion and myofibrillar protein degradation was decreased. Co-infusion of glucose + EAA prevented the loss of urinary nitrogen observed with EAA infusions alone and enhanced the increase in plasma IGF-1 concentration but there was no synergistic effect of glucose + EAA on the decrease in myofibrillar protein degradation. Muscle mRNA expression of the ubiquitin conjugating enzyme, 14-kDa E2 and proteasome sub-unit C2 were significantly decreased, after glucose but not amino acid infusions, and there was no further response to the combined infusions of glucose + EAA. Conclusion: Prolonged glucose infusion decreases myofibrillar protein degradation, prevents the excretion of infused EAA, and acts additively with EAA to increase plasma IGF-1 and improve net nitrogen balance. There was no evidence of synergistic effects between glucose + EAA infusion on muscle protein degradation or expression of components of the ubiquitin-proteasome proteolytic pathway.

Amino acid pairing preferences in parallel beta-sheets in proteins

Statistical approaches have been applied to examine amino acid pairing preferences within parallel beta-sheets. The main chain hydrogen bonding pattern in parallel beta-sheets means that, for each residue pair, only one of the residues is involved in main chain hydrogen bonding with the strand containing the partner residue. We call this the hydrogen bonded (HB) residue and the partner residue the non-hydrogen bonded (nHB) residue, and differentiate between the favorability of a pair and that of its reverse pair, e.g. Asn(HB)-Thr(nHB)versus Thr(HB)-Asn(nHB). Significantly (p < or = 0.000001) favoured pairings were rationalised using stereochemical arguments. For instance, Asn(HB)-Thr(nHB) and Arg(HB)-Thr(nHB) were favoured pairs, where the residues adopted favoured chi1 rotamer positions that allowed side-chain interactions to occur. In contrast, Thr(HB)-Asn(nHB) and Thr(HB)-Arg(nHB) were not significantly favoured, and could only form side-chain interactions if the residues involved adopted less favourable chi1 conformations. The favourability of hydrophobic pairs e.g. Ile(HB)-Ile(nHB), Val(HB)-Val(nHB) and Leu(HB)-Ile(nHB) was explained by the residues adopting their most preferred chi1 and chi2 conformations, which enabled them to form nested arrangements. Cysteine-cysteine pairs are significantly favoured, although these do not form intrasheet disulphide bridges. Interactions between positively and negatively charged residues were asymmetrically preferred: those with the negatively charged residue at the HB position were more favoured. This trend was accounted for by the presence of general electrostatic interactions, which, based on analysis of distances between charged atoms, were likely to be stronger when the negatively charged residue is the HB partner. The Arg(HB)-Asp(nHB) interaction was an exception to this trend and its favorability was rationalised by the formation of specific side-chain interactions. This research provides rules that could be applied to protein structure prediction, comparative modelling and protein engineering and design. The methods used to analyse the pairing preferences are automated and detailed results are available (http://www.rubic.rdg.ac.uk/betapairprefsparallel/).

Amino acid pairing preferences in parallel beta-sheets in proteins

Statistical approaches have been applied to examine amino acid pairing preferences within parallel beta-sheets. The main chain hydrogen bonding pattern in parallel beta-sheets means that, for each residue pair, only one of the residues is involved in main chain hydrogen bonding with the strand containing the partner residue. We call this the hydrogen bonded (HB) residue and the partner residue the non-hydrogen bonded (nHB) residue, and differentiate between the favourability of a pair and that of its reverse pair, e.g. Asn(HB)-Thr(nHB) versus Thr(HB)-Asn(nHB). Significantly (p <= 0.000001) favoured pairings were rationalised using stereochemical arguments. For instance, Asn(HB)-Thr(nHB) and Arg(HB)-Thr(nHB) were favoured pairs, where the residues adopted favoured chi(1) rotamer positions that allowed side-chain interactions to occur. In contrast, Thr(HB)-Asn(nHB) and Thr(HB)-Arg(nHB) were not significantly favoured, and could only form side-chain interactions if the residues involved adopted less favourable chi(1) conformations. The favourability of hydrophobic pairs e.g. Ile(HB)-Ile(nHB), Val(HB)-Val(nHB) and Leu(HB)-Ile(nHB) was explained by the residues adopting their most preferred chi(1) and chi(2) conformations, which enabled them to form nested arrangements. Cysteine-cysteine pairs are significantly favoured, although these do not form intrasheet disulphide bridges. Interactions between positively and negatively charged residues were asymmetrically preferred: those with the negatively charged residue at the HB position were more favoured. This trend was accounted for by the presence of general electrostatic interactions, which, based on analysis of distances between charged atoms, were likely to be stronger when the negatively charged residue is the HB partner. The Arg(HB)-Asp(nHB) interaction was an exception to this trend and its favourability was rationalised by the formation of specific side-chain interactions. This research provides rules that could be applied to protein structure prediction, comparative modelling and protein engineering and design. The methods used to analyse the pairing preferences are automated and detailed results are available (http:// www.rubic.rdg.ac.uk/betapairprefsparallel/). (c) 2005 Elsevier Ltd. All rights reserved.

Amino-acid cycling drives nitrogen fixation in the legume - Rhizobium symbiosis

The biological reduction of atmospheric N-2 to ammonium (nitrogen fixation) provides about 65% of the biosphere's available nitrogen. Most of this ammonium is contributed by legume rhizobia symbioses(1), which are initiated by the infection of legume hosts by bacteria (rhizobia), resulting in formation of root nodules. Within the nodules, rhizobia are found as bacteroids, which perform the nitrogen fixation: to do this, they obtain sources of carbon and energy from the plant, in the form of dicarboxylic acids(2,3). It has been thought that, in return, bacteroids simply provide the plant with ammonium. But here we show that a more complex amino-acid cycle is essential for symbiotic nitrogen fixation by Rhizobium in pea nodules. The plant provides amino acids to the bacteroids, enabling them to shut down their ammonium assimilation. In return, bacteroids act like plant organelles to cycle amino acids back to the plant for asparagine synthesis. The mutual dependence of this exchange prevents the symbiosis being dominated by the plant, and provides a selective pressure for the evolution of mutualism.