Xavier de la Cruz

PMUT: a web-based tool for the annotation of pathological mutations on proteins

PMUT allows the fast and accurate prediction (approximately 80% success rate in humans) of the pathological character of single point amino acidic mutations based on the use of neural networks. The program also allows the fast scanning of mutational hot spots, which are obtained by three procedures: (1) alanine scanning, (2) massive mutation and (3) genetically accessible mutations. A graphical interface for Protein Data Bank (PDB) structures, when available, and a database containing hot spot profiles for all non-redundant PDB structures are also accessible from the PMUT server.

PupasView: a visual tool for selecting suitable SNPs, with putative pathological effect in genes, for genotyping purposes

We have developed a web tool, PupasView, for the selection of single nucleotide polymorphisms (SNPs) with potential phenotypic effect. PupasView constitutes an interactive environment in which functional information and population frequency data can be used as sequential filters over linkage disequilibrium parameters to obtain a final list of SNPs optimal for genotyping purposes. PupasView is the first resource that integrates phenotypic effects caused by SNPs at both the translational and the transcriptional level. PupasView retrieves SNPs that could affect conserved regions that the cellular machinery uses for the correct processing of genes (intron/exon boundaries or exonic splicing enhancers), predicted transcription factor binding sites and changes in amino acids in the proteins for which a putative pathological effect is calculated. The program uses the mapping of SNPs in the genome provided by Ensembl. PupasView will be of much help in studies of multifactorial disorders, where the use of functional SNPs will increase the sensitivity of the identification of the genes responsible for the disease. The PupasView web interface is accessible through and through .

Cerebrospinal fluid-derived circulating tumour DNA better represents the genomic alterations of brain tumours than plasma

Cell-free circulating tumour DNA (ctDNA) in plasma has been shown to be informative of the genomic alterations present in tumours and has been used to monitor tumour progression and response to treatments. However, patients with brain tumours do not present with or present with low amounts of ctDNA in plasma precluding the genomic characterization of brain cancer through plasma ctDNA. Here we show that ctDNA derived from central nervous system tumours is more abundantly present in the cerebrospinal fluid (CSF) than in plasma. Massively parallel sequencing of CSF ctDNA more comprehensively characterizes the genomic alterations of brain tumours than plasma, allowing the identification of actionable brain tumour somatic mutations. We show that CSF ctDNA levels longitudinally fluctuate in time and follow the changes in brain tumour burden providing biomarkers to monitor brain malignancies. Moreover, CSF ctDNA is shown to facilitate and complement the diagnosis of leptomeningeal carcinomatosis.

The (In)dependence of alternative splicing and gene duplication

Alternative splicing (AS) and gene duplication (GD) both are processes that diversify the protein repertoire. Recent examples have shown that sequence changes introduced by AS may be comparable to those introduced by GD. In addition, the two processes are inversely correlated at the genomic scale: large gene families are depleted in splice variants and vice versa. All together, these data strongly suggest that both phenomena result in interchangeability between their effects. Here, we tested the extent to which this applies with respect to various protein characteristics. The amounts of AS and GD per gene are anticorrelated even when accounting for different gene functions or degrees of sequence divergence. In contrast, the two processes appear to be independent in their influence on variation in mRNA expression. Further, we conducted a detailed comparison of the effect of sequence changes in both alternative splice variants and gene duplicates on protein structure, in particular the size, location, and types of sequence substitutions and insertions/deletions. We find that, in general, alternative splicing affects protein sequence and structure in a more drastic way than gene duplication and subsequent divergence. Our results reveal an interesting paradox between the anticorrelation of AS and GD at the genomic level, and their impact at the protein level, which shows little or no equivalence in terms of effects on protein sequence, structure, and function. We discuss possible explanations that relate to the order of appearance of AS and GD in a gene family, and to the selection pressure imposed by the environment.

Autoacetylation regulates P/CAF nuclear localization.

Acetylation is a posttranslational modification that alters the biological activities of proteins by affecting their association with other proteins or DNA, their catalytic activities, or their subcellular distribution. The acetyltransferase P/CAF is autoacetylated and acetylated by p300 in vivo. P/CAF autoacetylation is an intramolecular or intermolecular event. Intramolecular acetylation targets five lysines within the nuclear localization signal at the P/CAF C terminus. We analyzed how the subcellular distribution of P/CAF is regulated by intramolecular autoacetylation and found that a P/CAF mutant lacking histone acetyltransferase activity accumulated primarily in the cytoplasm. This cytoplasmic fraction of P/CAF is enriched for nonautoacetylated P/CAF. In addition, P/CAF deacetylation by HDAC3 and in a minor degree by HDAC1, HDAC2, or HDAC4 leads to cytoplasmic accumulation of P/CAF. Importantly, our data show that P/CAF accumulates in the cytoplasm during apoptosis. These results reveal the molecular mechanism of autoacetylation control of P/CAF nuclear translocation and suggest a novel pathway by which P/CAF activity is controlled in vivo.

Classical molecular interaction potentials: Improved setup procedure in molecular dynamics simulations of proteins

The latest version of the classical molecular interaction potential (CMIP) has the ability to predict the position of crystallographic waters in several proteins with great accuracy. This article analyzes the ability of the CMIP functional to improve the setup procedure of the molecular system in molecular dynamics (MD) simulations of proteins. To this end, the CMIP strategy is used to include both water molecules and counterions in different protein systems. The structural details of the configurations sampled from trajectories obtained using the CMIP setup procedure are compared with those obtained from trajectories derived from a standard equilibration process. The results show that standard MD simulations can lead to artifactual results, which are avoided using the CMIP setup procedure. Because the CMIP is easy to implement at a low computational cost, it can be very useful in obtaining reliable MD trajectories. Proteins 2001;45:428–437.

Genome-wide analysis reveals that Smad3 and JMJD3 HDM co-activate the neural developmental program

Neural development requires crosstalk between signaling pathways and chromatin. In this study, we demonstrate that neurogenesis is promoted by an interplay between the TGFβ pathway and the H3K27me3 histone demethylase (HDM) JMJD3. Genome-wide analysis showed that JMJD3 is targeted to gene promoters by Smad3 in neural stem cells (NSCs) and is essential to activate TGFβ-responsive genes. In vivo experiments in chick spinal cord revealed that the generation of neurons promoted by Smad3 is dependent on JMJD3 HDM activity. Overall, these findings indicate that JMJD3 function is required for the TGFβ developmental program to proceed.

Scoring by Intermolecular Pairwise Propensities of Exposed Residues (SIPPER): A New Efficient Potential for Protein―Protein Docking

A detailed and complete structural knowledge of the interactome is one of the grand challenges in Biology, and a variety of computational docking approaches have been developed to complement experimental efforts and help in the characterization of protein-protein interactions. Among the different docking scoring methods, those based on physicochemical considerations can give the maximum accuracy at the atomic level, but they are usually computationally demanding and necessarily noisy when implemented in rigid-body approaches. Coarser-grained knowledge-based potentials are less sensitive to details of atomic arrangements, thus providing an efficient alternative for scoring of rigid-body docking poses. In this study, we have extracted new statistical potentials from intermolecular pairs of exposed residues in known complex structures, which were then used to score protein-protein docking poses. The new method, called SIPPER (scoring by intermolecular pairwise propensities of exposed residues), combines the value of residue desolvation based on solvent-exposed area with the propensity-based contribution of intermolecular residue pairs. This new scoring function found a near-native orientation within the top 10 predictions in nearly one-third of the cases of a standard docking benchmark and proved to be also useful as a filtering step, drastically reducing the number of docking candidates needed by energy-based methods like pyDock.

Ligand-induced changes in the binding sites of proteins

Classical molecular interaction potentials, in conjunction with other theoretical techniques, are used to analyze the dependence of the binding sites of representative proteins on the bound ligand. It is found that the ligand bound introduces in general small structural perturbations at the binding site of the protein. However, such small structural changes can lead to important alterations in the recognition pattern of the protein. The impact of these findings in docking procedures is discussed.

An Atomistic View to the Gas Phase Proteome

Extended all-atom molecular dynamics simulations on all protein metafolds have been performed to obtain a complete picture of the gas phase proteome. The structural atlas of the gas phase proteome obtained here shows an unexpected maintenance of the global and local structure and of the general deformability pattern upon transfer to the gas phase under electrospray conditions. Despite a general compression, the solution structure can be easily very well recognized from the gas phase one, and most structural details, such as secondary structure, are well preserved upon vaporization. Rehydration of the gas phase protein leads in most cases to a very fast transition from gas phase to solution structure. Overall, our massive analysis (over 4 micros in solution and over 12 micros in the gas phase) demonstrates that solution-like structures can be determined by using mass spectroscopy and related techniques to obtain fast approximations to the solution structure.