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.

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.

Supramolecular peptide helix from a novel double turn forming peptide containing a beta-amino acid

A single-crystal X-ray diffraction study of the terminally protected tetrapeptide Boc-beta-Ala-Aib-Leu-Aib-OMe 1 (Aib: alpha-aminoisobutyric acid; beta-Ala: beta-Alanine) reveals that it adopts a new type of double turn structure which self-associates to form a unique supramolecular helix through intermolecular hydrogen bonds. Scanning electron microscopic studies show that peptide 1 exhibits amyloid-like fibrillar morphology in the solid state. (C) 2003 Elsevier Science Ltd. All rights reserved.

Chiral poly(aromatic amide ester) dendrimers bearing an amino acid derived C-3-symmetric core: synthesis and properties

The effect of hyperbranched macromolecular architectures (dendrimers) upon chirality has received significant attention in recent years in the light of the proposal of amplification of chirality. In particular, several studies have been carried out on the chiroptical properties of dendrimers that contain a chiral core and achiral branches in order to determine if the chirality of the central core can be transmitted to the distal. region of the macromolecule. In addition to interest of a pure academic nature, the presence of such chiral conformational order would be extremely useful in the development of asymmetric catalysts. In this paper, a novel class of chiral dendrimers is described - these perfect hyperbranched macromolecules have been prepared by a convergent route by the coupling of a chiral central core based upon tris(2-aminoethyl)amine and poly(aromatic amide ester) dendritic branches. The chiral properties of these dendrimers have been investigated by detailed optical rotation studies and circular dichroism analysis; the results of these studies are described herein. (C) Wiley-VCH Verlag GmbH Co.

Casein and whey exert different effects on plasma amino acid profiles, gastrointestinal hormone secretion and appetite

Protein, generally agreed to be the most satiating macronutrient, may differ in its effects on appetite depending on the protein source and variation in digestion and absorption. We investigated the effects of two milk protein types, casein and whey, on food intake and subjective ratings of hunger and fullness, and on postprandial metabolite and gastrointestinal hormone responses. Two studies were undertaken. The first study showed that energy intake from a buffet meal ad libitum was significantly less 90 min after a 1700 kJ liquid preload containing 48 g whey, compared with an equivalent casein preload (P<0.05). In the second study, the same whey preload led to a 28 % increase in postprandial plasma amino acid concentrations over 3 h compared with casein (incremental area under the curve (iAUC), P<0.05). Plasma cholecystokinin (CCK) was increased by 60 % (iAUC, P<0.005), glucagon-like peptide (GLP)-1 by 65 % (iAUC, P<0.05) and glucose-dependent insulinotropic polypeptide by 36 % (iAUC, P<0.01) following the whey preload compared with the casein. Gastric emptying was influenced by protein type as evidenced by differing plasma paracetamol profiles with the two preloads. Greater subjective satiety followed the whey test meal (P<0.05). These results implicate post-absorptive increases in plasma amino acids together with both CCK and GLP-1 as potential mediators of the increased satiety response to whey and emphasise the importance of considering the impact of protein type on the appetite response to a mixed meal.

Effects of genotype and environment on free amino acid levels in wheat grain: implications for acrylamide formation during processing

Acrylamide forms from free asparagine and reducing sugars during cooking, with asparagine concentration being the key parameter determining the formation in foods produced from wheat flour. In this study free amino acid concentrations were measured in the grain of varieties Spark and Rialto and four doubled haploid lines from a Spark x Rialto mapping population. The parental and doubled haploid lines had differing levels of total free amino acids and free asparagine in the grain, with one line consistently being lower than either parent for both of these factors. Sulfur deprivation led to huge increases in the concentrations of free asparagine and glutamine, and canonical variate analysis showed clear separation of the grain samples as a result of treatment (environment, E) and genotype (G) and provided evidence of G x E interactions. Low grain sulfur and high free asparagine concentration were closely associated with increased risk of acrylamide formation. G, E, and G x E effects were also evident in grain from six varieties of wheat grown at field locations around the United Kingdom in 2006 and 2007. The data indicate that progress in reducing the risk of acrylamide formation in processed wheat products could be made immediately through the selection and cultivation of low grain asparagme varieties and that further genetically driven improvements should be achievable. However, genotypes that are selected should also be tested under a range of environmental conditions.

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

Identification of genetic and phenotypic differences associated with prevalent and non-prevalent Salmonella Enteritidis phage types: analysis of variation in amino acid transport

In this study, differences at the genetic level of 37 Salmonella Enteritidis strains from five phage types (PTs) were compared using comparative genomic hybridization (CGH) to assess differences between PTs. There were approximately 400 genes that differentiated prevalent (4, 6, 8 and 13a) and sporadic (11) PTs, of which 35 were unique to prevalent PTs, including six plasmid-borne genes, pefA, B, C, D, srgC and rck, and four chromosomal genes encoding putative amino acid transporters. Phenotype array studies also demonstrated that strains from prevalent PTs were less susceptible to urea stress and utilized L-histidine, L-glutamine, L-proline, L-aspartic acid, gly-asn and gly-gln more efficiently than PT11 strains. Complementation of a PT11 strain with the transporter genes from PT4 resulted in a significant increase in utilization of the amino acids and reduced susceptibility to urea stress. In epithelial cell association assays, PT11 strains were less invasive than other prevalent PTs. Most strains from prevalent PTs were better biofilm formers at 37 degrees C than at 28 degrees C, whilst the converse was true for PT11 strains. Collectively, the results indicate that genetic and corresponding phenotypic differences exist between strains of the prevalent PTs 4, 6, 8 and 13a and non-prevalent PT11 strains that are likely to provide a selective advantage for strains from the former PTs and could help them to enter the food chain and cause salmonellosis.

Evidence for the complex relationship between free amino acid and sugar concentrations and acrylamide-forming potential in potato

Free amino acids and reducing sugars participate in the Maillard reaction during high-temperature cooking and processing. This results not only in the formation of colour, aroma and flavour compounds, but also undesirable contaminants, including acrylamide, which forms when the amino acid that participates in the reaction is asparagine. In this study, tubers of 13 varieties of potato (Solanum tuberosum), which had been produced in a field trial in 2010 and sampled immediately after harvest or after storage for 6 months, were analysed to show the relationship between the concentrations of free asparagine, other free amino acids, sugars and acrylamide-forming potential. The varieties comprised five that are normally used for crisping, seven that are used for French fry production and one that is used for boiling. Acrylamide formation was measured in heated flour, and correlated with glucose and fructose concentration. In French fry varieties, which contain higher concentrations of sugars, acrylamide formation also correlated with free asparagine concentration, demonstrating the complex relationship between precursor concentration and acrylamide-forming potential in potato. Storage of the potatoes for 6 months at 9°C had a significant, variety-dependent impact on sugar and amino acid concentrations and acrylamide-forming potential.

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/