Dissecting the role of low-complexity regions in the evolution of vertebrate proteins

Low-complexity regions (LCRs) in proteins are tracts that are highly enriched in one or a few aminoacids. Given their high abundance, and their capacity to expand in relatively short periods of time through replication slippage, they can greatly contribute to increase protein sequence space and generate novel protein functions. However, little is known about the global impact of LCRs on protein evolution. We have traced back the evolutionary history of 2,802 LCRs from a large set of homologous protein families from H.sapiens, M.musculus, G.gallus, D.rerio and C.intestinalis. Transcriptional factors and other regulatory functions are overrepresented in proteins containing LCRs. We have found that the gain of novel LCRs is frequently associated with repeat expansion whereas the loss of LCRs is more often due to accumulation of amino acid substitutions as opposed to deletions. This dichotomy results in net protein sequence gain over time. We have detected a significant increase in the rate of accumulation of novel LCRs in the ancestral Amniota and mammalian branches, and a reduction in the chicken branch. Alanine and/or glycine-rich LCRs are overrepresented in recently emerged LCR sets from all branches, suggesting that their expansion is better tolerated than for other LCR types. LCRs enriched in positively charged amino acids show the contrary pattern, indicating an important effect of purifying selection in their maintenance. We have performed the first large-scale study on the evolutionary dynamics of LCRs in protein families. The study has shown that the composition of an LCR is an important determinant of its evolutionary pattern.

Contribution of S6K1/MAPK signaling pathways in the response to oxidative stress: activation of RSK and MSK by hydrogen peroxide

Cells respond to different kind of stress through the coordinated activation of signaling pathways such as MAPK or p53. To find which molecular mechanisms are involved, we need to understand their cell adaptation. The ribosomal protein, S6 kinase 1 (S6K1), is a common downstream target of signaling by hormonal or nutritional stress. Here, we investigated the initial contribution of S6K1/MAPK signaling pathways in the cell response to oxidative stress produced by hydrogen peroxide (H2O2). To analyze S6K1 activation, we used the commercial anti-phospho-Thr389-S6K1 antibody most frequently mentioned in the bibliography. We found that this antibody detected an 80-90 kDa protein that was rapidly phosphorylated in response to H2O2 in several human cells. Unexpectedly, this phosphorylation was insensitive to both mTOR and PI3K inhibitors, and knock-down experiments showed that this protein was not S6K1. RSK and MSK proteins were candidate targets of this phosphorylation. We demonstrated that H2O2 stimulated phosphorylation of RSK and MSK kinases at residues that are homologous to Thr389 in S6K1. This phosphorylation required the activity of either p38 or ERK MAP kinases. Kinase assays showed activation of RSK and MSK by H2O2. Experiments with mouse embryonic fibroblasts from p38 animals" knockout confirmed these observations. Altogether, these findings show that the S6K1 signaling pathway is not activated under these conditions, clarify previous observations probably misinterpreted by non-specific detection of proteins RSK and MSK by the anti-phospho-Thr389-S6K1 antibody, and demonstrate the specific activation of MAPK signaling pathways through ERK/p38/RSK/MSK by H2O2.

Elastic properties and secondary structure formation of single-stranded DNA at monovalent and divalent salt conditions

Single-stranded DNA (ssDNA) plays a major role in several biological processes. It is therefore of fundamental interest to understand how the elastic response and the formation of secondary structures are modulated by the interplay between base pairing and electrostatic interactions. Here we measure force-extension curves (FECs) of ssDNA molecules in optical tweezers set up over two orders of magnitude of monovalent and divalent salt conditions, and obtain its elastic parameters by fitting the FECs to semiflexible models of polymers. For both monovalent and divalent salts, we find that the electrostatic contribution to the persistence length is proportional to the Debye screening length, varying as the inverse of the square root of cation concentration. The intrinsic persistence length is equal to 0.7 nm for both types of salts, and the effectivity of divalent cations in screening electrostatic interactions appears to be 100-fold as compared with monovalent salt, in line with what has been recently reported for single-stranded RNA. Finally, we propose an analysis of the FECs using a model that accounts for the effective thickness of the filament at low salt condition and a simple phenomenological description that quantifies the formation of non-specific secondary structure at low forces.

PhylomeDB v4: zooming into the plurality of evolutionary histories of a genome

Phylogenetic trees representing the evolutionary relationships of homologous genes are the entry point for many evolutionary analyses. For instance, the use of a phylogenetic tree can aid in the inference of orthology and paralogy relationships, and in the detection of relevant evolutionary events such as gene family expansions and contractions, horizontal gene transfer, recombination or incomplete lineage sorting. Similarly, given the plurality of evolutionary histories among genes encoded in a given genome, there is a need for the combined analysis of genome-wide collections of phylogenetic trees (phylomes). Here, we introduce a new release of PhylomeDB (http://phylomedb.org), a public repository of phylomes. Currently, PhylomeDB hosts 120 public phylomes, comprising >1.5 million maximum likelihood trees and multiple sequence alignments. In the current release, phylogenetic trees are annotated with taxonomic, protein-domain arrangement, functional and evolutionary information. PhylomeDB is also a major source for phylogeny-based predictions of orthology and paralogy, covering >10 million proteins across 1059 sequenced species. Here we describe newly implemented PhylomeDB features, and discuss a benchmark of the orthology predictions provided by the database, the impact of proteome updates and the use of the phylome approach in the analysis of newly sequenced genomes and transcriptomes.

Analysis of meiotic recombination in 22q11.2, a region that frequently undergoes deletions and duplications

Background: The 22q11.2 deletion syndrome is the most frequent genomic disorder with an estimated frequency of 1/4000 live births. The majority of patients (90%) have the same deletion of 3 Mb (Typically Deleted Region, TDR) that results from aberrant recombination at meiosis between region specific low-copy repeats (LCRs). Methods: As a first step towards the characterization of recombination rates and breakpoints within the 22q11.2 region we have constructed a high resolution recombination breakpoint map based on pedigree analysis and a population-based historical recombination map based on LD analysis. Results: Our pedigree map allows the location of recombination breakpoints with a high resolution (potential recombination hotspots), and this approach has led to the identification of 5 breakpoint segments of 50 kb or less (8.6 kb the smallest), that coincide with historical hotspots. It has been suggested that aberrant recombination leading to deletion (and duplication) is caused by low rates of Allelic Homologous Recombination (AHR) within the affected region. However, recombination rate estimates for 22q11.2 region show that neither average recombination rates in the 22q11.2 region or within LCR22-2 (the LCR implicated in most deletions and duplications), are significantly below chromosome 22 averages. Furthermore, LCR22-2, the repeat most frequently implicated in rearrangements, is also the LCR22 with the highest levels of AHR. In addition, we find recombination events in the 22q11.2 region to cluster within families. Within this context, the same chromosome recombines twice in one family; first by AHR and in the next generation by NAHR resulting in an individual affected with the del22q11.2 syndrome. Conclusion: We show in the context of a first high resolution pedigree map of the 22q11.2 region that NAHR within LCR22 leading to duplications and deletions cannot be explained exclusively under a hypothesis of low AHR rates. In addition, we find that AHR recombination events cluster within families. If normal and aberrant recombination are mechanistically related, the fact that LCR22s undergo frequent AHR and that we find familial differences in recombination rates within the 22q11.2 region would have obvious health-related implications.

ORMDL proteins are a conserved new family of endoplasmic reticulum membrane proteins

Background: Annotations of completely sequenced genomes reveal that nearly half of the genes identified are of unknown function, and that some belong to uncharacterized gene families. To help resolve such issues, information can be obtained from the comparative analysis of homologous genes in model organisms. Results: While characterizing genes from the retinitis pigmentosa locus RP26 at 2q31-q33, we have identified a new gene, ORMDL1, that belongs to a novel gene family comprising three genes in humans (ORMDL1, ORMDL2 and ORMDL3), and homologs in yeast, microsporidia, plants, Drosophila, urochordates and vertebrates. The human genes are expressed ubiquitously in adult and fetal tissues. The Drosophila ORMDL homolog is also expressed throughout embryonic and larval stages, particularly in ectodermally derived tissues. The ORMDL genes encode transmembrane proteins anchored in the endoplasmic reticulum (ER). Double knockout of the two Saccharomyces cerevisiae homologs leads to decreased growth rate and greater sensitivity to tunicamycin and dithiothreitol. Yeast mutants can be rescued by human ORMDL homologs. Conclusions: From protein sequence comparisons we have defined a novel gene family, not previously recognized because of the absence of a characterized functional signature. The sequence conservation of this family from yeast to vertebrates, the maintenance of duplicate copies in different lineages, the ubiquitous pattern of expression in human and Drosophila, the partial functional redundancy of the yeast homologs and phenotypic rescue by the human homologs, strongly support functional conservation. Subcellular localization and the response of yeast mutants to specific agents point to the involvement of ORMDL in protein folding in the ER.

Nota. Detección molecular de Heterobasidion annosum s.s. en claras de Pinus nigra en el Pirineo Aragonés

En diciembre de 2007, en la Sierra de Loarre (Huesca) observamos signos de Newrobasidion annosum sensu lato asociados a daños de defoliación y mortalidad en una repoblación de Pinus nigra. Se ha identificado la especie biológica de aislamientos recogidos de este patógeno mediante técnicas moleculares, y se ha corroborado mediante cruces con aislamientos homocariontes, que se trata de H. annosum sensu stricto. Esta especie biológica ataca al género Pinus y puede causar daños sobre Pinus nigra.

Evolution and cellular function of monothiol glutaredoxins: involvement in iron-sulphur cluster assembly

A number of bacterial species, mostly proteobacteria, possess monothiol glutaredoxins homologous to the Saccharomyces cerevisiae mitochondrial protein Grx5, which is involved in iron–sulphur cluster synthesis. Phylogenetic profiling is used to predict that bacterial monothiol glutaredoxins also participate in the iron–sulphur cluster (ISC) assembly machinery, because their phylogenetic profiles are similar to the profiles of the bacterial homologues of yeast ISC proteins. High evolutionary cooccurrence is observed between the Grx5 homologues and the homologues of the Yah1 ferredoxin, the scaffold proteins Isa1 and Isa2, the frataxin protein Yfh1 and the Nfu1 protein. This suggests that a specific functional interaction exists between these ISC machinery proteins. Physical interaction analyses using low-definition protein docking predict the formation of strong and specific complexes between Grx5 and several components of the yeast ISC machinery. Two-hybrid analysis has confirmed the in vivo interaction between Grx5 and Isa1. Sequence comparison techniques and cladistics indicate that the other two monothiol glutaredoxins of S. cerevisiae, Grx3 and Grx4, have evolved from the fusion of a thioredoxin gene with a monothiol glutaredoxin gene early in the eukaryotic lineage, leading to differential functional specialization. While bacteria do not contain these chimaeric glutaredoxins, in many eukaryotic species Grx5 and Grx3/4-type monothiol glutaredoxins coexist in the cell.

Newly described human polyomaviruses Merkel Cell, KI and WU are present in urban sewage and may represent potential environmental contaminants

Recently, three new polyomaviruses (KI, WU and Merkel cell polyomavirus) have been reported to infect humans. It has also been suggested that lymphotropic polyomavirus, a virus of simian origin, infects humans. KI and WU polyomaviruses have been detected mainly in specimens from the respiratory tract while Merkel cell polyomavirus has been described in a very high percentage of Merkel cell carcinomas. The distribution, excretion level and transmission routes of these viruses remain unknown. Here we analyzed the presence and characteristics of newly described human polyomaviruses in urban sewage and river water in order to assess the excretion level and the potential role of water as a route of transmission of these viruses. Nested-PCR assays were designed for the sensitive detection of the viruses studied and the amplicons obtained were confirmed by sequencing analysis. The viruses were concentrated following a methodology previously developed for the detection of JC and BK human polyomaviruses in environmental samples. JC polyomavirus and human adenoviruses were used as markers of human contamination in the samples. Merkel cell polyomavirus was detected in 7/8 urban sewage samples collected and in 2/7 river water samples. Also one urine sample from a pregnant woman, out of 4 samples analyzed, was positive for this virus. KI and WU polyomaviruses were identified in 1/8 and 2/8 sewage samples respectively. The viral strains detected were highly homologous with other strains reported from several other geographical areas. Lymphotropic polyomavirus was not detected in any of the 13 sewage neither in 9 biosolid/sludge samples analyzed. This is the first description of a virus isolated from sewage and river water with a strong association with cancer. Our data indicate that the Merkel cell polyomavirus is prevalent in the population and that it may be disseminated through the fecal/urine contamination of water. The procedure developed may constitute a useful tool for studying the excreted strains, prevalence and transmission of these recently described polyomaviruses.

Periodontal health and esthetic results in impacted teeth exposed by apically positioned flap technique

Objectives: This study evaluates the periodontal health status and the esthetic results of teeth subjected to orthodontic traction, after their exposure by an apically positioned flap. Study design: Fifteen patients were included in the study, ages between 11 and 28 years old. The fenestrated teeth and their homologous contralateral normally erupted teeth, used as control, were evaluated. Results: Statistically significant differences were found in the position of the gingival margin (p = 0.005), with an average distance between cemento-enamel junction (CEJ) and gingival margin of 2.47 mm (SD 1.19) in control teeth and of 1 mm (SD 1.31) in the operated teeth, and in the depth of palatal probing (p = 0.031), with 2.1 mm (SD 0.9) for the experimental teeth and 1.7 mm (SD 0.8) for the control teeth. The gingival index, the bleeding during probing and the probing depth did not show statistically significant differences. The patient"s subjective esthetic evaluation was more favorable for the control teeth in most of the cases. Conclusions: The surgical approach for the impacted teeth by means of the apically positioned flap resulted to be a predictable technique allowing the maintenance of the periodontal health on a long-term basis.

Microscopic-macroscopic approach for binding energies with the Wigner-Kirkwood method

The semiclassical Wigner-Kirkwood ̄h expansion method is used to calculate shell corrections for spherical and deformed nuclei. The expansion is carried out up to fourth order in ̄h. A systematic study of Wigner-Kirkwood averaged energies is presented as a function of the deformation degrees of freedom. The shell corrections, along with the pairing energies obtained by using the Lipkin-Nogami scheme, are used in the microscopic-macroscopic approach to calculate binding energies. The macroscopic part is obtained from a liquid drop formula with six adjustable parameters. Considering a set of 367 spherical nuclei, the liquid drop parameters are adjusted to reproduce the experimental binding energies, which yields a root mean square (rms) deviation of 630 keV. It is shown that the proposed approach is indeed promising for the prediction of nuclear masses.

The structure of human tRNALys3 anticodon bound to HIV genome is stabilized by modified nucleosides and adjacent mismatch base pairs

Replication of human immunodeficiency virus (HIV) requires base pairing of the reverse transcriptase primer, human tRNA(Lys3), to the viral RNA. Although the major complementary base pairing occurs between the HIV primer binding sequence (PBS) and the tRNA's 3'-terminus, an important discriminatory, secondary contact occurs between the viral A-rich Loop I, 5'-adjacent to the PBS, and the modified, U-rich anticodon domain of tRNA(Lys3). The importance of individual and combined anticodon modifications to the tRNA/HIV-1 Loop I RNA's interaction was determined. The thermal stabilities of variously modified tRNA anticodon region sequences bound to the Loop I of viral sub(sero)types G and B were analyzed and the structure of one duplex containing two modified nucleosides was determined using NMR spectroscopy and restrained molecular dynamics. The modifications 2-thiouridine, s(2)U(34), and pseudouridine, Psi(39), appreciably stabilized the interaction of the anticodon region with the viral subtype G and B RNAs. The structure of the duplex results in two coaxially stacked A-form RNA stems separated by two mismatched base pairs, U(162)*Psi(39) and G(163)*A(38), that maintained a reasonable A-form helix diameter. The tRNA's s(2)U(34) stabilized the interaction between the A-rich HIV Loop I sequence and the U-rich anticodon, whereas the tRNA's Psi(39) stabilized the adjacent mismatched pairs.

Elastic properties and secondary structure formation of single-stranded DNA at monovalent and divalent salt conditions

Single-stranded DNA (ssDNA) plays a major role in several biological processes. It is therefore of fundamental interest to understand how the elastic response and the formation of secondary structures are modulated by the interplay between base pairing and electrostatic interactions. Here we measure force-extension curves (FECs) of ssDNA molecules in optical tweezers set up over two orders of magnitude of monovalent and divalent salt conditions, and obtain its elastic parameters by fitting the FECs to semiflexible models of polymers. For both monovalent and divalent salts, we find that the electrostatic contribution to the persistence length is proportional to the Debye screening length, varying as the inverse of the square root of cation concentration. The intrinsic persistence length is equal to 0.7 nm for both types of salts, and the effectivity of divalent cations in screening electrostatic interactions appears to be 100-fold as compared with monovalent salt, in line with what has been recently reported for single-stranded RNA. Finally, we propose an analysis of the FECs using a model that accounts for the effective thickness of the filament at low salt condition and a simple phenomenological description that quantifies the formation of non-specific secondary structure at low forces.

Synthesis of Janus compounds for the recognition of G-U mismatched nucleobase pairs

The design and synthesis of two Janus-type heterocycles with the capacity to simultaneously recognize guanine and uracyl in G-U mismatched pairs through complementary hydrogen bond pairing is described. Both compounds were conveniently functionalized with a carboxylic function and efficiently attached to a tripeptide sequence by using solid-phase methodologies. Ligands based on the derivatization of such Janus compounds with a small aminoglycoside, neamine, and its guanidinylated analogue have been synthesized, and their interaction with Tau RNA has been investigated by using several biophysical techniques, including UV-monitored melting curves, fluorescence titration experiments, and 1H NMR. The overall results indicated that Janus-neamine/guanidinoneamine showed some preference for the +3 mutated RNA sequence associated with the development of some tauopathies, although preliminary NMR studies have not confirmed binding to G-U pairs. Moreover, a good correlation has been found between the RNA binding affinity of such Janus-containing ligands and their ability to stabilize this secondary structure upon complexation.

Fast and economic immobilization methods described for non-commercial Pseudomonas lipases

There is an increasing interest to seek new enzyme preparations for the development of new products derived from bioprocesses to obtain alternative bio-based materials. In this context, four non-commercial lipases from Pseudomonas species were prepared, immobilized on different low-cost supports, and examined for potential biotechnological applications. Results: To reduce costs of eventual scaling-up, the new lipases were obtained directly from crude cell extracts or from growth culture supernatants, and immobilized by simple adsorption on Accurel EP100, Accurel MP1000 and Celite (R) 545. The enzymes evaluated were LipA and LipC from Pseudomonas sp. 42A2, a thermostable mutant of LipC, and LipI. 3 from Pseudomonas CR611, which were produced in either homologous or heterologous hosts. Best immobilization results were obtained on Accurel EP100 for LipA and on Accurel MP1000 for LipC and its thermostable variant. Lip I. 3, requiring a refolding step, was poorly immobilized on all supports tested ( best results for Accurel MP1000). To test the behavior of immobilized lipases, they were assayed in triolein transesterification, where the best results were observed for lipases immobilized on Accurel MP1000. Conclusions: The suggested protocol does not require protein purification and uses crude enzymes immobilized by a fast adsorption technique on low-cost supports, which makes the method suitable for an eventual scaling up aimed at biotechnological applications. Therefore, a fast, simple and economic method for lipase preparation and immobilization has been set up. The low price of the supports tested and the simplicity of the procedure, skipping the tedious and expensive purification steps, will contribute to cost reduction in biotechnological lipase-catalyzed processes.