David Abia

Two New Alleles of the abscisic aldehyde oxidase 3 Gene Reveal Its Role in Abscisic Acid Biosynthesis in Seeds

The abscisic aldehyde oxidase 3 (AAO3) gene product of Arabidopsis catalyzes the final step in abscisic acid (ABA) biosynthesis. An aao3-1 mutant in a Landsberg erecta genetic background exhibited a wilty phenotype in rosette leaves, whereas seed dormancy was not affected (Seo et al., 2000a). Therefore, it was speculated that a different aldehyde oxidase would be the major contributor to ABA biosynthesis in seeds (Seo et al., 2000a). Through a screening based on germination under high-salt concentration, we isolated two mutants in a Columbia genetic background, initially named sre2-1 and sre2-2 (for salt resistant). Complementation tests with different ABA-deficient mutants indicated that sre2-1 and sre2-2 mutants were allelic to aao3-1, and therefore they were renamed as aao3-2 and aao3-3, respectively. Indeed, molecular characterization of the aao3-2 mutant revealed a T-DNA insertional mutation that abolished the transcription of AAO3 gene, while sequence analysis of AAO3 in aao3-3 mutant revealed a deletion of three nucleotides and several missense mutations. Physiological characterization of aao3-2 and aao3-3 mutants revealed a wilty phenotype and osmotolerance in germination assays. In contrast to aao3-1, both aao3-2 and aao3-3 mutants showed a reduced dormancy. Accordingly, ABA levels were reduced in dry seeds and rosette leaves of both aao3-2 and aao3-3. Taken together, these results indicate that AAO3 gene product plays a major role in seed ABA biosynthesis.

Cross-Over between Discrete and Continuous Protein Structure Space: Insights into Automatic Classification and Networks of Protein Structures

Structural classifications of proteins assume the existence of the fold, which is an intrinsic equivalence class of protein domains. Here, we test in which conditions such an equivalence class is compatible with objective similarity measures. We base our analysis on the transitive property of the equivalence relationship, requiring that similarity of A with B and B with C implies that A and C are also similar. Divergent gene evolution leads us to expect that the transitive property should approximately hold. However, if protein domains are a combination of recurrent short polypeptide fragments, as proposed by several authors, then similarity of partial fragments may violate the transitive property, favouring the continuous view of the protein structure space. We propose a measure to quantify the violations of the transitive property when a clustering algorithm joins elements into clusters, and we find out that such violations present a well defined and detectable cross-over point, from an approximately transitive regime at high structure similarity to a regime with large transitivity violations and large differences in length at low similarity. We argue that protein structure space is discrete and hierarchic classification is justified up to this cross-over point, whereas at lower similarities the structure space is continuous and it should be represented as a network. We have tested the qualitative behaviour of this measure, varying all the choices involved in the automatic classification procedure, i.e., domain decomposition, alignment algorithm, similarity score, and clustering algorithm, and we have found out that this behaviour is quite robust. The final classification depends on the chosen algorithms. We used the values of the clustering coefficient and the transitivity violations to select the optimal choices among those that we tested. Interestingly, this criterion also favours the agreement between automatic and expert classifications. As a domain set, we have selected a consensus set of 2,890 domains decomposed very similarly in SCOP and CATH. As an alignment algorithm, we used a global version of MAMMOTH developed in our group, which is both rapid and accurate. As a similarity measure, we used the size-normalized contact overlap, and as a clustering algorithm, we used average linkage. The resulting automatic classification at the cross-over point was more consistent than expert ones with respect to the structure similarity measure, with 86% of the clusters corresponding to subsets of either SCOP or CATH superfamilies and fewer than 5% containing domains in distinct folds according to both SCOP and CATH. Almost 15% of SCOP superfamilies and 10% of CATH superfamilies were split, consistent with the notion of fold change in protein evolution. These results were qualitatively robust for all choices that we tested, although we did not try to use alignment algorithms developed by other groups. Folds defined in SCOP and CATH would be completely joined in the regime of large transitivity violations where clustering is more arbitrary. Consistently, the agreement between SCOP and CATH at fold level was lower than their agreement with the automatic classification obtained using as a clustering algorithm, respectively, average linkage (for SCOP) or single linkage (for CATH). The networks representing significant evolutionary and structural relationships between clusters beyond the cross-over point may allow us to perform evolutionary, structural, or functional analyses beyond the limits of classification schemes. These networks and the underlying clusters are available at http://ub.cbm.uam.es/research/ProtNet.php

Modifications in host cell cytoskeleton structure and function mediated by intracellular HIV-1 Tat protein are greatly dependent on the second coding exon

The human immunodeficiency virus type 1 (HIV-1) regulator Tat is essential for viral replication because it achieves complete elongation of viral transcripts. Tat can be released to the extracellular space and taken up by adjacent cells, exerting profound cytoskeleton rearrangements that lead to apoptosis. In contrast, intracellular Tat has been described as protector from apoptosis. Tat gene is composed by two coding exons that yield a protein of 101 amino acids (aa). First exon (1–72aa) is sufficient for viral transcript elongation and second exon (73–101 aa) appears to contribute to non-transcriptional functions. We observed that Jurkat cells stably expressing intracellular Tat101 showed gene expression deregulation 4-fold higher than cells expressing Tat72. Functional experiments were performed to evaluate the effect of this deregulation. First, NF-κB-, NF-AT- and Sp1-dependent transcriptional activities were greatly enhanced in Jurkat-Tat101, whereas Tat72 induced milder but efficient activation. Second, cytoskeleton-related functions as cell morphology, proliferation, chemotaxis, polarization and actin polymerization were deeply altered in Jurkat-Tat101, but not in Jurkat-Tat72. Finally, expression of several cell surface receptors was dramatically impaired by intracellular Tat101 but not by Tat72. Consequently, these modifications were greatly dependent on Tat second exon and they could be related to the anergy observed in HIV-1-infected T cells.

The amino terminal domain from Mrt4 protein can functionally replace the RNA binding domain of the ribosomal P0 protein

In Saccharomyces cerevisiae, the Mrt4 protein is a component of the ribosome assembly machinery that shares notable sequence homology to the P0 ribosomal stalk protein. Here, we show that these proteins can not bind simultaneously to ribosomes and moreover, a chimera containing the first 137 amino acids of Mrt4 and the last 190 amino acids from P0 can partially complement the absence of the ribosomal protein in a conditional P0 null mutant. This chimera is associated with ribosomes isolated from this strain when grown under restrictive conditions, although its binding is weaker than that of P0. These ribosomes contain less P1 and P2 proteins, the other ribosomal stalk components. Similarly, the interaction of the L12 protein, a stalk base component, is affected by the presence of the chimera. These results indicate that Mrt4 and P0 bind to the same site in the 25S rRNA. Indeed, molecular dynamics simulations using modelled Mrt4 and P0 complexes provide further evidence that both proteins bind similarly to rRNA, although their interaction with L12 displays notable differences. Together, these data support the participation of the Mrt4 protein in the assembly of the P0 protein into the ribosome and probably, that also of the L12 protein.

MM-ISMSA: An Ultrafast and Accurate Scoring Function for Protein-Protein Docking.

An ultrafast and accurate scoring function for protein-protein docking is presented. It includes (1) a molecular mechanics (MM) part based on a 12-6 Lennard-Jones potential; (2) an electrostatic component based on an implicit solvent model (ISM) with individual desolvation penalties for each partner in the protein-protein complex plus a hydrogen bonding term; and (3) a surface area (SA) contribution to account for the loss of water contacts upon protein-protein complex formation. The accuracy and performance of the scoring function, termed MM-ISMSA, have been assessed by (1) comparing the total binding energies, the electrostatic term, and its components (charge-charge and individual desolvation energies), as well as the per residue contributions, to results obtained with well-established methods such as APBSA or MM-PB(GB)SA for a set of 1242 decoy protein-protein complexes and (2) testing its ability to recognize the docking solution closest to the experimental structure as that providing the most favorable total binding energy. For this purpose, a test set consisting of 15 protein-protein complexes with known 3D structure mixed with 10 decoys for each complex was used. The correlation between the values afforded by MM-ISMSA and those from the other methods is quite remarkable (r(2) ∼ 0.9), and only 0.2-5.0 s (depending on the number of residues) are spent on a single calculation including an all vs all pairwise energy decomposition. On the other hand, MM-ISMSA correctly identifies the best docking solution as that closest to the experimental structure in 80% of the cases. Finally, MM-ISMSA can process molecular dynamics trajectories and reports the results as averaged values with their standard deviations. MM-ISMSA has been implemented as a plugin to the widely used molecular graphics program PyMOL, although it can also be executed in command-line mode. MM-ISMSA is distributed free of charge to nonprofit organizations.

Expression analysis revealing destabilizing mutations in phosphomannomutase 2 deficiency (PMM2-CDG)

Deficiency of phosphomannomutase (PMM2, MIM#601785) is the most common congenital disorder of glycosylation. Herein we report the genetic analysis of 22 Spanish PMM2 deficient patients and the functional analysis of 14 nucleotide changes in a prokaryotic expression system in order to elucidate their molecular pathogenesis. PMM2 activity assay revealed the presence of six protein changes with no enzymatic activities (p.R123Q, p.R141H, p.F157S, p.P184T, p.F207S and p.D209G) and seven mild protein changes with residual activities ranging from 16 to 54% (p.L32R, p.V44A p.D65Y, p.P113L p.T118S, p.T237M and p.C241S) and also one variant change with normal activity (p.E197A). The results obtained from Western blot analysis, degradation time courses of 11 protein changes and structural analysis of the PMM2 protein, suggest that the loss-of-function of most mutant proteins is based on their increased susceptibility to degradation or aggregation compared to the wild type protein, considering PMM2 deficiency as a conformational disease. We have identified exclusively catalytic protein change (p.D209G), catalytic protein changes affecting protein stability (p.R123Q and p.R141H), two protein changes disrupting the dimer interface (p.P113L and p.T118S) and several misfolding changes (p.L32R, p.V44A, p.D65Y, p.F157S, p.P184T, p.F207S, p.T237M and p.C241S). Our current work opens a promising therapeutic option using pharmacological chaperones to revert the effect of the characterized misfolding mutations identified in a wide range of PMM2 deficient patients.

Quantifying the evolutionary divergence of protein structures: The role of function change and function conservation

The molecular clock hypothesis, stating that protein sequences diverge in evolution by accumulating amino acid substitutions at an almost constant rate, played a major role in the development of molecular evolution and boosted quantitative theories of evolutionary change. These studies were extended to protein structures by the seminal paper by Chothia and Lesk, which established the approximate proportionality between structure and sequence divergence. Here we analyse how function influences the relationship between sequence and structure divergence, studying four large superfamilies of evolutionarily related proteins: globins, aldolases, P‐loop and NADP‐binding. We introduce the contact divergence, which is more consistent with sequence divergence than previously used structure divergence measures. Our main findings are: (1) Small structure and sequence divergences are proportional, consistent with the molecular clock. Approximate validity of the clock is also supported by the analysis of the clustering coefficient of structure similarity networks. (2) Functional constraints strongly limit the structure divergence of proteins performing the same function and may allow to identify incomplete or wrong functional annotations. (3) The rate of structure versus sequence divergence is larger for proteins performing different functions than for proteins performing the same function. We conjecture that this acceleration is due to positive selection for new functions. Accelerations in structure divergence are also suggested by the analysis of the clustering coefficient. (4) For low sequence identity, structural diversity explodes. We conjecture that this explosion is related to functional diversification. (5) Large indels are almost always associated with function changes. Proteins 2010.

Mobility of the native Bacillus subtilis conjugative plasmid pLS20 is regulated by intercellular signaling.

Horizontal gene transfer mediated by plasmid conjugation plays a significant role in the evolution of bacterial species, as well as in the dissemination of antibiotic resistance and pathogenicity determinants. Characterization of their regulation is important for gaining insights into these features. Relatively little is known about how conjugation of Gram-positive plasmids is regulated. We have characterized conjugation of the native Bacillus subtilis plasmid pLS20. Contrary to the enterococcal plasmids, conjugation of pLS20 is not activated by recipient-produced pheromones but by pLS20-encoded proteins that regulate expression of the conjugation genes. We show that conjugation is kept in the default “OFF” state and identified the master repressor responsible for this. Activation of the conjugation genes requires relief of repression, which is mediated by an anti-repressor that belongs to the Rap family of proteins. Using both RNA sequencing methodology and genetic approaches, we have determined the regulatory effects of the repressor and anti-repressor on expression of the pLS20 genes. We also show that the activity of the anti-repressor is in turn regulated by an intercellular signaling peptide. Ultimately, this peptide dictates the timing of conjugation. The implications of this regulatory mechanism and comparison with other mobile systems are discussed.

Functional characterization of novel genotypes and cellular oxidative stress studies in propionic acidemia.

Propionic acidemia (PA), caused by a deficiency of the mitochondrial biotin dependent enzyme propionyl-CoA carboxylase (PCC) is one of the most frequent organic acidurias in humans. PA is caused by mutations in either the PCCA or PCCB genes encoding the α- and β-subunits of the PCC enzyme which are assembled as an α6β6 dodecamer. In this study we have investigated the molecular basis of the defect in ten fibroblast samples from PA patients. Using homology modeling with the recently solved crystal structure of the PCC holoenzyme and a eukaryotic expression system we have analyzed the structural and functional effect of novel point mutations, also revealing a novel splice defect by minigene analysis. In addition, we have investigated the contribution of oxidative stress to cellular damage measuring reactive oxygen species (ROS) levels and apoptosis parameters in patient fibroblasts, as recent studies point to a secondary mitochondrial dysfunction as pathophysiological mechanism in this disorder. The results show an increase in intracellular ROS content compared to controls, correlating with the activation of the JNK and p38 signaling pathways. Highest ROS levels were present in cells harboring functionally null mutations, including one severe missense mutation. This work provides molecular insight into the pathogenicity of PA variants and indicates that oxidative stress may be a major contributing factor to the cellular damage, supporting the proposal of antioxidant strategies as novel supplementary therapy in this rare disease.

An interfering activity against lymphocytic choriomeningitis virus replication associated with enhanced mutagenesis

Previous studies have documented that, in the presence of the mutagenic base analogue 5-fluorouracil (FU), lymphocytic choriomeningitis virus (LCMV) that persisted in BHK-21 cells decreased its infectivity to a larger extent than intracellular viral RNA levels, prior to virus extinction. This observation, together with in silico simulations, led to the proposal of the lethal defection model of virus extinction. This model suggests the participation of defective-interfering genomes in the loss of infectivity by increased mutagenesis. Since LCMV naturally produces defective-interfering particles, it was important to show that a capacity to interfere is produced in association with FU treatment. Here, we document that BHK-21 cells persistently infected with LCMV grown in the presence of FU, but not in its absence, generated an interfering activity that suppressed LCMV infectivity. Interference was specific for LCMV and was sensitive to UV irradiation and its activity was dose- and time-dependent. The interfering preparations produced positive LCMV immunofluorescence and viral particles seen by electron microscopy when used to infect cells, despite some preparations being devoid of detectable infectivity. Interference did not involve significant increases of mutant spectrum complexity, as predicted by the lethal defection model. The results provide support for a specific interference associated with LCMV when the virus replicates in the presence of FU. The excess of interference relative to that observed in the absence of FU is necessary for LCMV extinction.