Searching for signatures of cold adaptations in modern and archaic humans: hints from the brown adipose tissue genes.

Adaptation to low temperatures has been reasonably developed in the human species during the colonization of the Eurasian landmass subsequent to Out of Africa migrations of anatomically modern humans. In addition to morphological and cultural changes, also metabolic ones are supposed to have favored human isolation from cold and body heat production and this can be hypothesized also for most Neandertal and at least for some Denisovan populations, which lived in geographical areas that strongly experienced the last glacial period. Modulation of non-shivering thermogenesis, for which adipocytes belonging to the brown adipose tissue are the most specialized cells, might have driven these metabolic adaptations. To perform an exploratory analysis aimed at looking into this hypothesis, variation at 28 genes involved in such functional pathway was investigated in modern populations from different climate zones, as well as in Neandertal and Denisovan genomes. Patterns of variation at the LEPR gene, strongly related to increased heat dissipation by mitochondria, appeared to have been shaped by positive selection in modern East Asians, but not in Europeans. Moreover, a single potentially cold-adapted LEPR allele, different from the supposed adaptive one identified in Homo sapiens, was found also in Neandertal and Denisovan genomes. These findings suggest that independent mechanisms for cold adaptations might have been developed in different non-African human groups, as well as that the evolution of possible enhanced thermal efficiency in Neandertals and in some Denisovan populations has plausibly entailed significant changes also in other functional pathways than in the examined one.

High-throughput SNP discovery in the rabbit (Oryctolagus cuniculus) genome by next-generation semiconductor-based sequencing.

The European rabbit (Oryctolagus cuniculus) is a domesticated species with one of the broadest ranges of economic and scientific applications and fields of investigation. Rabbit genome information and assembly are available (oryCun2.0), but so far few studies have investigated its variability, and massive discovery of polymorphisms has not been published yet for this species. Here, we sequenced two reduced representation libraries (RRLs) to identify single nucleotide polymorphisms (SNPs) in the rabbit genome. Genomic DNA of 10 rabbits belonging to different breeds was pooled and digested with two restriction enzymes (HaeIII and RsaI) to create two RRLs which were sequenced using the Ion Torrent Personal Genome Machine. The two RRLs produced 2 917 879 and 4 046 871 reads, for a total of 280.51 Mb (248.49 Mb with quality >20) and 417.28 Mb (360.89 Mb with quality >20) respectively of sequenced DNA. About 90% and 91% respectively of the obtained reads were mapped on the rabbit genome, covering a total of 15.82% of the oryCun2.0 genome version. The mapping and ad hoc filtering procedures allowed to reliably call 62 491 SNPs. SNPs in a few genomic regions were validated by Sanger sequencing. The Variant Effect Predictor Web tool was used to map SNPs on the current version of the rabbit genome. The obtained results will be useful for many applied and basic research programs for this species and will contribute to the development of cost-effective solutions for high-throughput SNP genotyping in the rabbit.

Exploring copy number variation in the rabbit (Oryctolagus cuniculus) genome by array comparative genome hybridization.

The European rabbit (Oryctolagus cuniculus) is relevant in a large spectrum of fields: it is a livestock, a pet, a biomedical model and a biotechnology tool, a wild resource and a pest. The sequencing of the rabbit genome has opened new perspectives to study this lagomorph at the genome level. We herein investigated for the first time the O. cuniculus genome by array comparative genome hybridization (aCGH) and established a first copy number variation (CNV) genome map in this species comprising 155 copy number variation regions (CNVRs; 95 gains, 59 losses, 1 with both gain and loss) covering ~0.3% of the OryCun2.0 version. About 50% of the 155 CNVRs identified spanned 139 different protein coding genes, 110 genes of which were annotated or partially annotated (including Major Histocompatibility Complex genes) with 277 different gene ontology terms. Many rabbit CNVRs might have a functional relevance that should be further investigated.

A first comparative map of copy number variations in the sheep genome.

We carried out a cross species cattle-sheep array comparative genome hybridization experiment to identify copy number variations (CNVs) in the sheep genome analysing ewes of Italian dairy or dual-purpose breeds (Bagnolese, Comisana, Laticauda, Massese, Sarda, and Valle del Belice) using a tiling oligonucleotide array with ~385,000 probes designed on the bovine genome. We identified 135 CNV regions (CNVRs; 24 reported in more than one animal) covering ~10.5 Mb of the virtual sheep genome referred to the bovine genome (0.398%) with a mean and a median equal to 77.6 and 55.9 kb, respectively. A comparative analysis between the identified sheep CNVRs and those reported in cattle and goat genomes indicated that overlaps between sheep and both other species CNVRs are highly significant (P<0.0001), suggesting that several chromosome regions might contain recurrent interspecies CNVRs. Many sheep CNVRs include genes with important biological functions. Further studies are needed to evaluate their functional relevance.

Molecular mechanism of thioredoxin regulation in photosynthetic A2B2-glyceraldehyde-3-phosphate dehydrogenase.

Chloroplast glyceraldehyde-3-phosphate dehydrogenase (GAPDH) is a light-regulated, NAD(P)H-dependent enzyme involved in plant photosynthetic carbon reduction. Unlike lower photosynthetic organisms, which only contain A(4)-GAPDH, the major GAPDH isoform of land plants is made up of A and B subunits, the latter containing a C-terminal extension (CTE) with fundamental regulatory functions. Light-activation of AB-GAPDH depends on the redox state of a pair of cysteines of the CTE, which can form a disulfide bond under control of thioredoxin f, leading to specific inhibition of the NADPH-dependent activity. The tridimensional structure of A(2)B(2)-GAPDH from spinach chloroplasts, crystallized in the oxidized state, shows that each disulfide-containing CTE is docked into a deep cleft between a pair of A and B subunits. The structure of the CTE was derived from crystallographic data and computational modeling and confirmed by site-specific mutagenesis. Structural analysis of oxidized A(2)B(2)-GAPDH and chimeric mutant [A+CTE](4)-GAPDH revealed that Arg-77, which is essential for coenzyme specificity and high NADPH-dependent activity, fails to interact with NADP in these kinetically inhibited GAPDH tetramers and is attracted instead by negative residues of oxidized CTE. Other subtle changes in catalytic domains and overall conformation of the tetramers were noticed in oxidized A(2)B(2)-GAPDH and [A+CTE](4)-GAPDH, compared with fully active A(4)-GAPDH. The CTE is envisioned as a redox-sensitive regulatory domain that can force AB-GAPDH into a kinetically inhibited conformation under oxidizing conditions, which also occur during dark inactivation of the enzyme in vivo.

K-Fold: a tool for the prediction of the protein folding kinetic order and rate.

UNLABELLED: K-Fold is a tool for the automatic prediction of the protein folding kinetic order and rate. The tool is based on a support vector machine (SVM) that was trained on a data set of 63 proteins, whose 3D structure and folding mechanism are known from experiments already described in the literature. The method predicts whether a protein of known atomic structure folds according to a two-state or a multi-state kinetics and correctly classifies 81% of the folding mechanisms when tested over the training set of the 63 proteins. It also predicts as a further option the logarithm of the folding rate. To the best of our knowledge, the tool discriminates for the first time whether a protein is characterized by a two state or a multiple state kinetics, during the folding process, and concomitantly estimates also the value of the constant rate of the process. When used to predict the logarithm of the folding rate, K-Fold scores with a correlation value to the experimental data of 0.74 (with a SE of 1.2). AVAILABILITY: http://gpcr.biocomp.unibo.it/cgi/predictors/K-Fold/K-Fold.cgi. SUPPLEMENTARY INFORMATION: http://gpcr.biocomp.unibo.it/~emidio/K-Fold/K-Fold_help.html.

The effect of backbone on the small-world properties of protein contact maps.

In the last years, small-world behavior has been extensively described for proteins, when they are represented by the undirected graph defined by the inter-residue protein contacts. By adopting this representation it was possible to compute the average clustering coefficient (C) and characteristic path length (L) of protein structures, and their values were found to be similar to those of graphs characterized by small-world topology. In this comment, we analyze a large set of non-redundant protein structures (1753) and show that by randomly mimicking the protein collapse, the covalent structure of the protein chain significantly contributes to the small-world behavior of the inter-residue contact graphs. When protein graphs are generated, imposing constraints similar to those induced by the backbone connectivity, their characteristic path lengths and clustering coefficients are indistinguishable from those computed using the real contact maps showing that L and C values cannot be used for 'protein fingerprinting'. Moreover we verified that these results are independent of the selected protein representations, residue composition and protein secondary structures.

Terlipressin infusion induces ischemia of breast skin in a cirrothic patient with hepatorenal syndrome.

We report a case of a 65-year-old woman with hepatitis C virus-related decompensated cirrhosis with hepatorenal syndrome, treated by high dose of terlipressin. Few hours after the highest dose was started, the patient complained burning pain in breasts, followed by the development of extensive bilateral cyanosis of breast's skin. When terlipressin was immediately stopped, pain and skin cyanosis rapidly disappeared. The peculiarity of our case is that cyanosis did not develop in common peripheral sites (e.g. fingers, toes, etc.) but in an atypical area, as skin of the breasts. Probably, this particular behaviour could be explained by the anatomical position of her large size breasts, that resulting as an extremely sloping and stretching region thus filling the maximum effect of gravity.

Predicting the insurgence of human genetic diseases associated to single point protein mutations with support vector machines and evolutionary information.

MOTIVATION: Human single nucleotide polymorphisms (SNPs) are the most frequent type of genetic variation in human population. One of the most important goals of SNP projects is to understand which human genotype variations are related to Mendelian and complex diseases. Great interest is focused on non-synonymous coding SNPs (nsSNPs) that are responsible of protein single point mutation. nsSNPs can be neutral or disease associated. It is known that the mutation of only one residue in a protein sequence can be related to a number of pathological conditions of dramatic social impact such as Alzheimer's, Parkinson's and Creutzfeldt-Jakob's diseases. The quality and completeness of presently available SNPs databases allows the application of machine learning techniques to predict the insurgence of human diseases due to single point protein mutation starting from the protein sequence. RESULTS: In this paper, we develop a method based on support vector machines (SVMs) that starting from the protein sequence information can predict whether a new phenotype derived from a nsSNP can be related to a genetic disease in humans. Using a dataset of 21 185 single point mutations, 61% of which are disease-related, out of 3587 proteins, we show that our predictor can reach more than 74% accuracy in the specific task of predicting whether a single point mutation can be disease related or not. Our method, although based on less information, outperforms other web-available predictors implementing different approaches. AVAILABILITY: A beta version of the web tool is available at http://gpcr.biocomp.unibo.it/cgi/predictors/PhD-SNP/PhD-SNP.cgi

DCCD-sensitive proton permeability of bacterial photosynthetic membranes. Cross-reconstitution studies with purified bovine heart Fo subunits.

The DCCD-sensitive proton permeability of chromatophores, from a green strain of Rhodobacter Capsulatus is potentiometrically detected following the proton release induced by a transmembrane diffusion potential imposed by a valinomycin-mediated potassium influx with a procedure already used for bovine heart submitochondrial particles (ESMP) and vesicles from Escherichia coli (Zanotti et al. (1994) Eur. J. Biochem. 222, 733-741). In the photosynthetic system, addition of increasing amounts of DCCD inhibits, with a similar titre, both proton permeability and MgATP-dependent ATPase activity as detected in the dark. The titre for 50% inhibition coincides with that obtained measuring proton permeability and ATP hydrolysis in ESMP. Upon removal of F1, the passive proton permeability is much less sensitive to DCCD in chromatophores than in USMP, suggesting that in chromatophores the F1-Fo interaction shapes the DCCD-sensitive proton conducting pathway. Addition of the purified mitochondrial FoI-PVP and oligomycin sensitivity-conferring (OSCP) proteins to the F1 stripped chromatophores restored the sensitivity of proton permeability to DCCD detected in untreated chromatophores. Analysis of the binding of 14C[DCCD] on F1 stripped chromatophores shows that the increase of DCCD sensitivity of proton permeability, caused by addition of mitochondrial Fo proteins, is related to an increase of the binding of the inhibitor to subunit c of Fo sector of ATP synthase complex.

Characterization of 9-aminoacridine interaction with chromatophore membranes and modelling of the probe response to artificially induced transmembrane delta pH values.

We analyze the adsorption of the fluorescent monoamine 9-aminoacridine to the membrane phase of photosynthetic chromatophores, in the physiological interval of pH values ranging from 5.5 to 8.5 and at ionic strengths of 0.005 and 0.150 M. The interaction of the probe with the membrane phase is described with S-shaped isotherms of the Hill type and is modulated by electrostatic effects as modelled with the Gouy-Chapman-Boltzman theory. This description is consistent with different values of the surface change density of the chromatophore membranes decreasing from about 1.3 x 10(-3) to about 0.5 x 10(-3) e-/A2, on changing the pH from 8.5/7.5 to 6.5/5.5, respectively. Furthermore we show that, when the free concentrations of the probe in the inner and outer vesicle compartments are computed from the adsorbing isotherms at the proper pH values, the model considering the equilibrium distribution of the neutral monoamine following the onset of a delta pH is sufficient to describe the dependence of the artificially induced transmembrane delta pH values on the observed quenching of the probe fluorescence.

Light-dark adaptation of bacteriorhodopsin in triton-treated purple membrane.

Solubilization of purple membrane with Triton X-100 yields Triton micelles containing bacteriorhodopsin monomers. The absorption maximum of dark-adapted solubilized bacteriorhodopsin is blue-shifted to 549 nm. Light adaption increases the absorbance by 4% and shifts the absorption maximum to 553 nm, i.e., the extent of light adaptation is considerably less than in intact purple membrane. Extraction of dark-adapted bacteriorhodopsin in Triton yields a 13-cis- to all-trans-retinal ratio of 58 : 42 which changes after light adaptation to 38 : 62. It has been shown by Sperling et al. (Sperling, W., Carl, P., Rafferty, Ch.N. and Dencher, N.A. (1977) Biophys. Struct. Mech. 3, 79-94) that light adaptation in intact purple membrane occurs through a branching of the 13-cis photoreaction cycle, so that part of the pigment during each cycle crosses over into the all-trans photoreaction cycle. We explain the decreased extent of light adaptation in solubilized bacteriorhodopsin by assuming a significant back reaction from the all-trans to the 13-cis cycle. This assumption predicts a wavelength dependence of the extent of light adaptation, which is born out by experiment.

Topology prediction for helical transmembrane proteins at 86% accuracy.

Previously, we introduced a neural network system predicting locations of transmembrane helices (HTMs) based on evolutionary profiles (PHDhtm, Rost B, Casadio R, Fariselli P, Sander C, 1995, Protein Sci 4:521-533). Here, we describe an improvement and an extension of that system. The improvement is achieved by a dynamic programming-like algorithm that optimizes helices compatible with the neural network output. The extension is the prediction of topology (orientation of first loop region with respect to membrane) by applying to the refined prediction the observation that positively charged residues are more abundant in extra-cytoplasmic regions. Furthermore, we introduce a method to reduce the number of false positives, i.e., proteins falsely predicted with membrane helices. The evaluation of prediction accuracy is based on a cross-validation and a double-blind test set (in total 131 proteins). The final method appears to be more accurate than other methods published: (1) For almost 89% (+/-3%) of the test proteins, all HTMs are predicted correctly. (2) For more than 86% (+/-3%) of the proteins, topology is predicted correctly. (3) We define reliability indices that correlate with prediction accuracy: for one half of the proteins, segment accuracy raises to 98%; and for two-thirds, accuracy of topology prediction is 95%. (4) The rate of proteins for which HTMs are predicted falsely is below 2% (+/-1%). Finally, the method is applied to 1,616 sequences of Haemophilus influenzae. We predict 19% of the genome sequences to contain one or more HTMs. This appears to be lower than what we predicted previously for the yeast VIII chromosome (about 25%).

Transmembrane helices predicted at 95% accuracy.

We describe a neural network system that predicts the locations of transmembrane helices in integral membrane proteins. By using evolutionary information as input to the network system, the method significantly improved on a previously published neural network prediction method that had been based on single sequence information. The input data were derived from multiple alignments for each position in a window of 13 adjacent residues: amino acid frequency, conservation weights, number of insertions and deletions, and position of the window with respect to the ends of the protein chain. Additional input was the amino acid composition and length of the whole protein. A rigorous cross-validation test on 69 proteins with experimentally determined locations of transmembrane segments yielded an overall two-state per-residue accuracy of 95%. About 94% of all segments were predicted correctly. When applied to known globular proteins as a negative control, the network system incorrectly predicted fewer than 5% of globular proteins as having transmembrane helices. The method was applied to all 269 open reading frames from the complete yeast VIII chromosome. For 59 of these, at least two transmembrane helices were predicted. Thus, the prediction is that about one-fourth of all proteins from yeast VIII contain one transmembrane helix, and some 20%, more than one.

Quantum mechanical analysis of oxygenated and deoxygenated states of hemocyanin: theoretical clues for a plausible allosteric model of oxygen binding.

In this work with ab initio computations, we describe relevant interactions between protein active sites and ligands, using as a test case arthropod hemocyanins. A computational analysis of models corresponding to the oxygenated and deoxygenated forms of the hemocyanin active site is performed using the Density Functional Theory approach. We characterize the electron density distribution of the binding site with and without bound oxygen in relation to the geometry, which stems out of the crystals of three hemocyanin proteins, namely the oxygenated form from the horseshoe crab Limulus polyphemus, and the deoxygenated forms, respectively, from the same source and from another arthropod, the spiny lobster Panulirus interruplus. Comparison of the three available crystals indicate structural differences at the oxygen binding site, which cannot be explained only by the presence and absence of the oxygen ligand, since the geometry of the ligand site of the deoxygenated Panulirus hemocyanin is rather similar to that of the oxygenated Limulus protein. This finding was interpreted in the frame of a mechanism of allosteric regulation for oxygen binding. However, the cooperative mechanism, which is experimentally well documented, is only partially supported by crystallographic data, since no oxygenated crystal of Panulirus hemocyanin is presently available. We address the following question: is the local ligand geometry responsible for the difference of the dicopper distance observed in the two deoxygenated forms of hemocyanin or is it necessary to advocate the allosteric regulation of the active site conformations in order to reconcile the different crystal forms? We find that the difference of the dicopper distance between the two deoxygenated hemocyanins is not due to the small differences of ligand geometry found in the crystals and conclude that it must be therefore stabilized by the whole protein tertiary structure.

Tryptophanyl fluorescence lifetime distribution of hyperthermophilic beta-glycosidase from molecular dynamics simulation: a comparison with the experimental data.

A molecular dynamics simulation approach has been utilized to understand the unusual fluorescence emission decay observed for beta-glycosidase from the hyperthermophilic bacterium Solfolobus sulfotaricus (Sbeta gly), a tetrameric enzyme containing 17 tryptophanyl residues for each subunit. The tryptophanyl emission decay of Sbeta gly results from a bimodal distribution of fluorescence lifetimes with a short-lived component centered at 2.5 ns and a long-lived one at 7.4 ns (Bismuto E, Nucci R, Rossi M, Irace G, 1999, Proteins 27:71-79). From the examination of the trajectories of the side chains capable of causing intramolecular quenching for each tryptophan microenvironment and using a modified Stern-Volmer model for the emission quenching processes, we calculated the fluorescence lifetime for each tryptophanyl residue of Sbeta gly at two different temperatures, i.e., 300 and 365 K. The highest temperature was chosen because in this condition Sbeta gly evidences a maximum in its catalytic activity and is stable for a very long time. The calculated lifetime distributions overlap those experimentally determined. Moreover, the majority of trytptophanyl residues having longer lifetimes correspond to those originally identified by inspection of the crystallographic structure. The tryptophanyl lifetimes appear to be a complex function of several variables, such as microenvironment viscosity, solvent accessibility, the chemical structure of quencher side chains, and side-chain dynamics. The lifetime calculation by MD simulation can be used to validate a predicted structure by comparing the theoretical data with the experimental fluorescence decay results.

Prediction of the transmembrane regions of beta-barrel membrane proteins with a neural network-based predictor.

A method based on neural networks is trained and tested on a nonredundant set of beta-barrel membrane proteins known at atomic resolution with a jackknife procedure. The method predicts the topography of transmembrane beta strands with residue accuracy as high as 78% when evolutionary information is used as input to the network. Of the transmembrane beta-strands included in the training set, 93% are correctly assigned. The predictor includes an algorithm of model optimization, based on dynamic programming, that correctly models eight out of the 11 proteins present in the training/testing set. In addition, protein topology is assigned on the basis of the location of the longest loops in the models. We propose this as a general method to fill the gap of the prediction of beta-barrel membrane proteins.

A neural network approach to evaluate fold recognition results.

Fold recognition techniques assist the exploration of protein structures, and web-based servers are part of the standard set of tools used in the analysis of biochemical problems. Despite their success, current methods are only able to predict the correct fold in a relatively small number of cases. We propose an approach that improves the selection of correct folds from among the results of two methods implemented as web servers (SAMT99 and 3DPSSM). Our approach is based on the training of a system of neural networks with models generated by the servers and a set of associated characteristics such as the quality of the sequence-structure alignment, distribution of sequence features (sequence-conserved positions and apolar residues), and compactness of the resulting models. Our results show that it is possible to detect adequate folds to model 80% of the sequences with a high level of confidence. The improvements achieved by taking into account sequence characteristics open the door to future improvements by directly including such factors in the step of model generation. This approach has been implemented as an automatic system LIBELLULA, available as a public web server at http://www.pdg.cnb.uam.es/servers/libellula.html.

Can functional regions of proteins be predicted from their coding sequences? The case study of G-protein coupled receptors.

A filter based on a set of unsupervised neural networks trained with a winner-take-all strategy discloses signals along the coding sequences of G-protein coupled receptors. By comparing with the existing experimental data it appears that these signals correlate with putative functional domains of the proteins. After protein alignment within subfamilies, signals cluster in protein regions which, according to the presently available experimental results, are described as possible functional domains of the folded proteins. The mapping procedure reveals characteristic regions in the coding sequences common and/or characteristic of the receptor subtype. This is particularly noticeable for the third cytoplasmic loop, which is likely to be involved in the molecular coupling of all the subfamilies with G-proteins. The results indicate that our mapping can highlight intrinsic representative features of the coding sequences which, in the case of G-protein coupled receptors, are characteristic of protein functional regions and suggest a possible application of the filter for predicting functional determinants in proteins starting from the coding sequence.

A high diffusion coefficient for coenzyme Q10 might be related to a folded structure.

We measured the lateral diffusion of different coenzyme Q homologues and analogues in model lipid vesicles using the fluorescence collisional quenching technique with pyrene derivatives and found diffusion coefficients in the range of 10(-6) cm2/s. Theoretical diffusion coefficients for these highly hydrophobic components were calculated according to the free volume theory. An important parameter in the free volume theory is the relative dimension between diffusant and solvent: a molecular dynamics computer simulation of the coenzymes yielded their most probable geometries and volumes and revealed surprisingly similar sizes of the short and long homologues, due to a folded structure of the isoprenoid chain in the latter, with a length for coenzyme Q10 of 21 A. Using this information we were able to calculate diffusion coefficients in the range of 10(-6) cm2/s, in good agreement with those found experimentally.