Cristina Marino-Buslje

Protein Conformational Diversity Correlates with Evolutionary Rate

Native state of proteins is better represented by an ensemble of conformers in equilibrium than by only one structure. The extension of structural differences between conformers characterizes the conformational diversity of the protein. In this study, we found a negative correlation between conformational diversity and protein evolutionary rate. Conformational diversity was expressed as the maximum root mean square deviation (RMSD) between the available conformers in Conformational Diversity of Native State database. Evolutionary rate estimations were calculated using 16 different species compared with human sharing at least 700 orthologous proteins with known conformational diversity extension. The negative correlation found is independent of the protein expression level and comparable in magnitude and sign with the correlation between gene expression level and evolutionary rate. Our findings suggest that the structural constraints underlying protein dynamism, essential for protein function, could modulate protein divergence.

Sequence evolution of the intrinsically disordered and globular domains of a model viral oncoprotein.

In the present work, we have used the papillomavirus E7 oncoprotein to pursue structure-function and evolutionary studies that take into account intrinsic disorder and the conformational diversity of globular domains. The intrinsically disordered (E7N) and globular (E7C) domains of E7 show similar degrees of conservation and co-evolution. We found that E7N can be described in terms of conserved and coevolving linear motifs separated by variable linkers, while sequence evolution of E7C is compatible with the known homodimeric structure yet suggests other activities for the domain. Within E7N, inter-residue relationships such as residue co-evolution and restricted intermotif distances map functional coupling and co-occurrence of linear motifs that evolve in a coordinate manner. Within E7C, additional cysteine residues proximal to the zinc-binding site may allow redox regulation of E7 function. Moreover, we describe a conserved binding site for disordered domains on the surface of E7C and suggest a putative target linear motif. Both homodimerization and peptide binding activities of E7C are also present in the distantly related host PHD domains, showing that these two proteins share not only structural homology but also functional similarities, and strengthening the view that they evolved from a common ancestor. Finally, we integrate the multiple activities and conformations of E7 into a hierarchy of structure-function relationships.

The residue 179 is involved in product specificity of the Bacillus circulans DF 9R cyclodextrin glycosyltransferase

Cyclodextrin glycosyltransferases (CGTases) are important enzymes in biotechnology because of their ability to produce cyclodextrin (CD) mixtures from starch whose relative composition depends on enzyme source. A multiple alignment of 46 CGTases and Shannon entropy analysis allowed us to find differences and similarities that could be related to product specificity. Interestingly, position 179 has Gly in all the CGTases except in that from Bacillus circulans DF 9R which possesses Gln. The absence of a side chain at that position has been considered as a strong requirement for substrate binding and cyclization process. Therefore, we constructed two mutants of this enzyme, Q179L and Q179G. The activity and kinetic parameters of Q179G remained unchanged while the Q179L mutant showed a different CDs ratio, a lower catalytic efficiency, and a decreased ability to convert starch into CDs. We show that position 179 is involved in CGTase product specificity and must be occupied by Gly—without a side chain—or by amino acid residues able to interact with the substrate through hydrogen bonds in a way that the cyclization process occurs efficiently. These findings are also explained on the basis of a structural model.

I-COMS: Interprotein-COrrelated Mutations Server

Interprotein contact prediction using multiple sequence alignments (MSAs) is a useful approach to help detect protein–protein interfaces. Different computational methods have been developed in recent years as an approximation to solve this problem. However, as there are discrepancies in the results provided by them, there is still no consensus on which is the best performing methodology. To address this problem, I-COMS (interprotein COrrelated Mutations Server) is presented. I-COMS allows to estimate covariation between residues of different proteins by four different covariation methods. It provides a graphical and interactive output that helps compare results obtained using different methods. I-COMS automatically builds the required MSA for the calculation and produces a rich visualization of either intraprotein and/or interprotein covariating positions in a circos representation. Furthermore, comparison between any two methods is available as well as the overlap between any or all four methodologies. In addition, as a complementary source of information, a matrix visualization of the corresponding scores is made available and the density plot distribution of the inter, intra and inter+intra scores are calculated. Finally, all the results can be downloaded (including MSAs, scores and graphics) for comparison and visualization and/or for further analysis.

Kin-Driver: a database of driver mutations in protein kinases

Somatic mutations in protein kinases (PKs) are frequent driver events in many human tumors, while germ-line mutations are associated with hereditary diseases. Here we present Kin-driver, the first database that compiles driver mutations in PKs with experimental evidence demonstrating their functional role. Kin-driver is a manual expert-curated database that pays special attention to activating mutations (AMs) and can serve as a validation set to develop new generation tools focused on the prediction of gain-of-function driver mutations. It also offers an easy and intuitive environment to facilitate the visualization and analysis of mutations in PKs. Because all mutations are mapped onto a multiple sequence alignment, analogue positions between kinases can be identified and tentative new mutations can be proposed for studying by transferring annotation. Finally, our database can also be of use to clinical and translational laboratories, helping them to identify uncommon AMs that can correlate with response to new antitumor drugs. The website was developed using PHP and JavaScript, which are supported by all major browsers; the database was built using MySQL server. Kin-driver is available at: http://kin-driver.leloir.org.ar/

The Arabidopsis DNA Polymerase δ Has a Role in the Deposition of Transcriptionally Active Epigenetic Marks, Development and Flowering

DNA replication is a key process in living organisms. DNA polymerase α (Polα) initiates strand synthesis, which is performed by Polε and Polδ in leading and lagging strands, respectively. Whereas loss of DNA polymerase activity is incompatible with life, viable mutants of Polα and Polε were isolated, allowing the identification of their functions beyond DNA replication. In contrast, no viable mutants in the Polδ polymerase-domain were reported in multicellular organisms. Here we identify such a mutant which is also thermosensitive. Mutant plants were unable to complete development at 28°C, looked normal at 18°C, but displayed increased expression of DNA replication-stress marker genes, homologous recombination and lysine 4 histone 3 trimethylation at the SEPALLATA3 (SEP3) locus at 24°C, which correlated with ectopic expression of SEP3. Surprisingly, high expression of SEP3 in vascular tissue promoted FLOWERING LOCUS T (FT) expression, forming a positive feedback loop with SEP3 and leading to early flowering and curly leaves phenotypes. These results strongly suggest that the DNA polymerase δ is required for the proper establishment of transcriptionally active epigenetic marks and that its failure might affect development by affecting the epigenetic control of master genes.

Activating Mutations Cluster in the “Molecular Brake” Regions of Protein Kinases and Do Not Associate with Conserved or Catalytic Residues

Mutations leading to activation of proto‐oncogenic protein kinases (PKs) are a type of drivers crucial for understanding tumorogenesis and as targets for antitumor drugs. However, bioinformatics tools so far developed to differentiate driver mutations, typically based on conservation considerations, systematically fail to recognize activating mutations in PKs. Here, we present the first comprehensive analysis of the 407 activating mutations described in the literature, which affect 41 PKs. Unexpectedly, we found that these mutations do not associate with conserved positions and do not directly affect ATP binding or catalytic residues. Instead, they cluster around three segments that have been demonstrated to act, in some PKs, as “molecular brakes” of the kinase activity. This finding led us to hypothesize that an auto inhibitory mechanism mediated by such “brakes” is present in all PKs and that the majority of activating mutations act by releasing it. Our results also demonstrate that activating mutations of PKs constitute a distinct group of drivers and that specific bioinformatics tools are needed to identify them in the numerous cancer sequencing projects currently underway. The clustering in three segments should represent the starting point of such tools, a hypothesis that we tested by identifying two somatic mutations in EPHA7 that might be functionally relevant.

The role of the N-terminal tail for the oligomerization, folding and stability of human frataxin.

The N‐terminal stretch of human frataxin (hFXN) intermediate (residues 42–80) is not conserved throughout evolution and, under defined experimental conditions, behaves as a random‐coil. Overexpression of hFXN56–210 in Escherichia coli yields a multimer, whereas the mature form of hFXN (hFXN81–210) is monomeric. Thus, cumulative experimental evidence points to the N‐terminal moiety as an essential element for the assembly of a high molecular weight oligomer. The secondary structure propensity of peptide 56–81, the moiety putatively responsible for promoting protein–protein interactions, was also studied. Depending on the environment (TFE or SDS), this peptide adopts α‐helical or β‐strand structure. In this context, we explored the conformation and stability of hFXN56–210. The biophysical characterization by fluorescence, CD and SEC‐FPLC shows that subunits are well folded, sharing similar stability to hFXN90–210. However, controlled proteolysis indicates that the N‐terminal stretch is labile in the context of the multimer, whereas the FXN domain (residues 81–210) remains strongly resistant. In addition, guanidine hydrochloride at low concentration disrupts intermolecular interactions, shifting the ensemble toward the monomeric form. The conformational plasticity of the N‐terminal tail might impart on hFXN the ability to act as a recognition signal as well as an oligomerization trigger. Understanding the fine‐tuning of these activities and their resulting balance will bear direct relevance for ultimately comprehending hFXN function.

MIToS.jl: mutual information tools for protein sequence analysis in the Julia language

Motivation MIToS is an environment for mutual information analysis and a framework for protein multiple sequence alignments (MSAs) and protein structures (PDB) management in Julia language. It integrates sequence and structural information through SIFTS, making Pfam MSAs analysis straightforward. MIToS streamlines the implementation of any measure calculated from residue contingency tables and its optimization and testing in terms of protein contact prediction. As an example, we implemented and tested a BLOSUM62-based pseudo-count strategy in mutual information analysis. Availability and Implementation The software is totally implemented in Julia and supported for Linux, OS X and Windows. Its freely available on GitHub under MIT license: http://mitos.leloir.org.ar . Contacts diegozea@gmail.com or cmb@leloir.org.ar. Supplementary information Supplementary data are available at Bioinformatics online.

How Does pH Fit in with Oscillating Polar Growth

Polar growth in root hairs and pollen tubes is an excellent model for investigating plant cell size regulation. While linear plant growth is historically explained by the acid growth theory, which considers that auxin triggers apoplastic acidification by activating plasma membrane P-type H+-ATPases (AHAs) along with cell wall relaxation over long periods, the apoplastic pH (apopH) regulatory mechanisms are unknown for polar growth. Polar growth is a fast process mediated by rapid oscillations that repeat every ∼20-40s. In this review, we explore a reactive oxygen species (ROS)-dependent mechanism that could generate oscillating apopH gradients in a coordinated manner with growth and Ca2+ oscillations. We propose possible mechanisms by which apopH oscillations are coordinated with polar growth together with ROS and Ca2+ waves.