Biomedical subjects

Julián Echave

Publications and source records attributed to Julián Echave.

8 recordsLinked to original sources

PubMed2006

Quaternary structure constraints on evolutionary sequence divergence.

The structurally constrained protein evolution (SCPE) model simulates protein divergence considering protein structure explicitly. The model is based on the observation that protein structure is more conserved during evolution than the sequences encoding for that structure. In the previous work, the SCPE model considered only the tertiary structure. Here we show that the performance of the model is enhanced when the oligomeric structure is taken into account. Our results agree with recent evolutionary studies of oligomeric proteins, which show that conservation of the quaternary structure imposes additional constraints on sequence divergence. The incorporation of protein-protein interactions into protein evolution models may be important in the study of quaternary protein structures and complex protein assemblies.

Algorithms↗

PubMed2006

Evolutionary conservation of protein backbone flexibility.

Internal protein dynamics is essential for biological function. During evolution, protein divergence is functionally constrained: properties more relevant for function vary more slowly than less important properties. Thus, if protein dynamics is relevant for function, it should be evolutionary conserved. In contrast with the well-studied evolution of protein structure, the evolutionary divergence of protein dynamics has not been addressed systematically before, apart from a few case studies. X-Ray diffraction analysis gives information not only on protein structure but also on B-factors, which characterize the flexibility that results from protein dynamics. Here we study the evolutionary divergence of protein backbone dynamics by comparing the C(alpha) flexibility (B-factor) profiles for a large dataset of homologous proteins classified into families and superfamilies. We show that C(alpha) flexibility profiles diverge slowly, so that they are conserved at family and superfamily levels, even for pairs of proteins with nonsignificant sequence similarity. We also analyze and discuss the correlations among the divergences of flexibility, sequence, and structure.

Amino Acid Sequence↗

PubMed2005

Gamma carbonic anhydrase like complex interact with plant mitochondrial complex I.

We report the identification by two hybrid screens of two novel similar proteins, called Arabidopsis thaliana gamma carbonic anhydrase like1 and 2 (AtgammaCAL1 and AtgammaCAL2), that interact specifically with putative Arabidopsis thaliana gamma Carbonic Anhydrase (AtgammaCA) proteins in plant mitochondria. The interaction region that was located in the N-terminal 150 amino acids of mature AtgammaCA and AtgammaCA like proteins represents a new interaction domain. In vitro experiments indicate that these proteins are imported into mitochondria and are associated with mitochondrial complex I as AtgammaCAs. All plant species analyzed contain both AtgammaCA and AtgammaCAL sequences indicating that these genes were conserved throughout plant evolution. Structural modeling of AtgammaCAL sequences show a deviation of functionally important active site residues with respect to gammaCAs but could form active interfaces in the interaction with AtgammaCAs. We postulate a CA complex tightly associated to plant mitochondrial complex.

Arabidopsis↗

PubMed2004

Generality of the structurally constrained protein evolution model: assessment on representatives of the four main fold classes.

The Structurally Constrained Protein Evolution (SCPE) model simulates protein evolution by introducing random mutations into the evolving sequences and selecting them against too much structural perturbation. Given a single protein structure, the SCPE model can be used to obtain a whole set of site-dependent amino acid substitution matrices. The set of SCPE substitution matrices for a given protein family can be seen as an independent-sites model of evolution for that family. Thus, these matrices can be compared with other substitution-matrix-based models of evolution. So far, SCPE has been tested only on left-handed parallel beta helix (LbetaH) proteins. Here, we address the question of generality by assessing the SCPE model on representatives of the four main classes of folds: alpha, beta, alpha+beta, and alpha/beta. We compare with other models using the likelihood ratio test with parametric bootstrapping. We show that SCPE performs better than the popular JTT model for all cases considered. Furthermore, by considering the relative contributions of mutation and selection, we found that the key to the success of the SCPE model is the selection step.

Algorithms↗

PubMed2004

The structurally constrained protein evolution model accounts for sequence patterns of the LbetaH superfamily.

BACKGROUND: Structure conservation constrains evolutionary sequence divergence, resulting in observable sequence patterns. Most current models of protein evolution do not take structure into account explicitly, being unsuitable for investigating the effects of structure conservation on sequence divergence. To this end, we recently developed the Structurally Constrained Protein Evolution (SCPE) model. The model starts with the coding sequence of a protein with known three-dimensional structure. At each evolutionary time-step of an SCPE simulation, a trial sequence is generated by introducing a random point mutation in the current coding DNA sequence. Then, a "score" for the trial sequence is calculated and the mutation is accepted only if its score is under a given cutoff, lambda. The SCPE score measures the distance between the trial sequence and a given reference sequence, given the structure. In our first brief report we used a "global score", in which the same reference sequence, the ancestral one, was used at each evolutionary step. Here, we introduce a new scoring function, the "local score", in which the sequence accepted at the previous evolutionary time-step is used as the reference. We assess the model on the UDP-N-acetylglucosamine acyltransferase (LPXA) family, as in our previous report, and we extend this study to all other members of the left-handed parallel beta helix fold (LbetaH) superfamily whose structure has been determined. RESULTS: We studied site-dependent entropies, amino acid probability distributions, and substitution matrices predicted by SCPE and compared with experimental data for several members of the LbetaH superfamily. We also evaluated structure conservation during simulations. Overall, SCPE outperforms JTT in the description of sequence patterns observed in structurally constrained sites. Maximum Likelihood calculations show that the local-score and global-score SCPE substitution matrices obtained for LPXA outperform the JTT model for the LPXA family and for the structurally constrained sites of class i of other members within the LbetaH superfamily. CONCLUSION: We extended the SCPE model by introducing a new scoring function, the local score. We performed a thorough assessment of the SCPE model on the LPXA family and extended it to all other members of known structure of the LbetaH superfamily.

Acyltransferases↗

PubMed2004

Dynactins p25 and p27 are predicted to adopt the LbetaH fold.

Dynactin is a multimeric protein essential for the minus-end-directed transport driven by microtubule-based motor dynein. The pointed-end subcomplex in dynactin contains p62, p27, p25, and Arp11 subunits, and is thought to participate in interactions with membranous cargoes. We used sequence and structure prediction analysis to study dynactins p25 and p27. Here we present evidence that strongly supports that dynactins p27 and p25 contain the isoleucine-patch motif and adopt the left-handed parallel beta-helix fold. The structural models we obtained could contribute to the understanding of the complex interactions that dynactins are able to establish with cargo particles, microtubules or other dynactin subunits.

Amino Acid Sequence↗

PubMed2004

Gamma carbonic anhydrases in plant mitochondria.

Three genes from Arabidopsis thaliana with high sequence similarity to gamma carbonic anhydrase (gammaCA), a Zn containing enzyme from Methanosarcina thermophila (CAM), were identified and characterized. Evolutionary and structural analyses predict that these genes code for active forms of gammaCA. Phylogenetic analyses reveal that these Arabidopsis gene products cluster together with CAM and related sequences from alpha and gamma proteobacteria, organisms proposed as the mitochondrial endosymbiont ancestor. Indeed, in vitro and in vivo experiments indicate that these gene products are transported into the mitochondria as occurs with several mitochondrial protein genes transferred, during evolution, from the endosymbiotic bacteria to the host genome. Moreover, putative CAM orthologous genes are detected in other plants and green algae and were predicted to be imported to mitochondria. Structural modeling and sequence analysis performed in more than a hundred homologous sequences show a high conservation of functionally important active site residues. Thus, the three histidine residues involved in Zn coordination (His 81, 117 and 122), Arg 59, Asp 61, Gin 75, and Asp 76 of CAM are conserved and properly arranged in the active site cavity of the models. Two other functionally important residues (Glu 62 and Glu 84 of CAM) are lacking, but alternative amino acids that might serve to their roles are postulated. Accordingly, we propose that photosynthetic eukaryotic organisms (green algae and plants) contain gammaCAs and that these enzymes codified by nuclear genes are imported into mitochondria to accomplish their biological function.

Amino Acid Sequence↗

PubMed2003

Sequence determinants of quaternary structure in lumazine synthase.

Riboflavin, an essential cofactor for all organisms, is biosynthesized in plants, fungi and microorganisms. The penultimate step in the pathway is catalyzed by the enzyme lumazine synthase. One of the most distinctive characteristics of this enzyme is that it is found in different species in two different quaternary structures, pentameric and icosahedral, built from practically the same structural monomeric unit. In fact, the icosahedral structure is best described as a capsid of twelve pentamers. Despite this noticeable difference, the active sites are virtually identical in all structurally studied members. Furthermore, the main regions involved in the catalysis are located at the interface between adjacent subunits in the pentamer. Thus, the two quaternary forms of the enzyme must meet similar structural requirements to achieve their function, but, at the same time, they should differ in the sequence traits responsible for the different quaternary structures observed. Here, we present a combined analysis that includes sequence-structure and evolutionary studies to find the sequence determinants of the different quaternary assemblies of this enzyme. A data set containing 86 sequences of the lumazine synthase family was recovered by sequence similarity searches. Seven of them had resolved three-dimensional structures. A subsequent phylogenetic reconstruction by maximum parsimony (MP) allowed division of the total set into two clusters in accord with their quaternary structure. The comparison between the patterns of three-dimensional contacts derived from the known three-dimensional structures and variation in sequence conservation revealed a significant shift in structural constraints of certain positions. Also, to explore the changes in functional constraints between the two groups, site-specific evolutionary rate shifts were analyzed. We found that the positions involved in icosahedral contacts suffer a larger increase in constraints than the rest. We found eight sequence sites that would be the most important icosahedral sequence determinants. We discuss our results and compare them with previous work. These findings should contribute to refinement of the current structural data, to the design of assays that explore the role of these positions, to the structural characterization of new sequences, and to initiation of a study of the underlying evolutionary mechanisms.

Base Sequence↗