Search PubMed⌕ Search

Biomedical subjects

Joost Schymkowitz

Publications and source records attributed to Joost Schymkowitz.

11 recordsLinked to original sources

FoldX force field revisited, an improved version.

MOTIVATION: The FoldX force field was originally validated with a database of 1000 mutants at a time when there were few high-resolution structures. Here, we have manually curated a database of 5556 mutants affecting protein stability, resulting in 2484 highly confident mutations denominated FoldX stability dataset (FSD), represented in non-redundant X-ray structures with <2.5&#x2009;&#xc5; resolution, not involving duplicates, metals, or prosthetic groups. Using this database, we have created a new version of the FoldX force field by introducing pi stacking, pH dependency for all charged residues, improving aromatic-aromatic interactions, modifying the Ncap contribution and &#x3b1;-helix dipole, recalibrating the side-chain entropy of methionine, adjusting the H-bond parameters, and modifying the solvation contribution of tryptophan and others. RESULTS: These changes have led to significant improvements for the prediction of specific mutants involving the above residues/interactions and a statistically significant increase of FoldX predictions, as well as for the majority of the 20 aa. Removing all training sets data from FSD [Validation FoldX Stability Dataset (VFSD) dataset] resulted in improved predictions from R&#x2009;=&#x2009;0.693 (RMSE&#x2009;=&#x2009;1.277&#x2009;kcal/mol) to R&#x2009;=&#x2009;0.706 (RMSE&#x2009;=&#x2009;1.252&#x2009;kcal/mol) when compared with the previously released version. FoldX achieves 95% accuracy considering an error of &#xb1;0.85&#x2009;kcal/mol in prediction and an area under the curve&#x2009;=&#x2009;0.78 for the VFSD, predicting the sign of the energy change upon mutation. AVAILABILITY AND IMPLEMENTATION: FoldX versions 4.1 and 5.1 are freely available for academics at https://foldxsuite.crg.eu/.

Databases, Protein↗

PupaSuite: finding functional single nucleotide polymorphisms for large-scale genotyping purposes.

We have developed a web tool, PupaSuite, for the selection of single nucleotide polymorphisms (SNPs) with potential phenotypic effect, specifically oriented to help in the design of large-scale genotyping projects. PupaSuite uses a collection of data on SNPs from heterogeneous sources and a large number of pre-calculated predictions to offer a flexible and intuitive interface for selecting an optimal set of SNPs. It improves the functionality of PupaSNP and PupasView programs and implements new facilities such as the analysis of user's data to derive haplotypes with functional information. A new estimator of putative effect of polymorphisms has been included that uses evolutionary information. Also SNPeffect database predictions have been included. The PupaSuite web interface is accessible through http://pupasuite.bioinfo.cipf.es and through http://www.pupasnp.org.

Algorithms↗

SNPeffect v2.0: a new step in investigating the molecular phenotypic effects of human non-synonymous SNPs.

UNLABELLED: Single nucleotide polymorphisms (SNPs) constitute the most fundamental type of genetic variation in human populations. About 75 000 of these reported variations cause an amino acid change in the translated protein. An important goal in genomic research is to understand how this variability affects protein function, and whether or not particular SNPs are associated to disease susceptibility. Accordingly, the SNPeffect database uses sequence- and structure-based bioinformatics tools to predict the effect of non-synonymous SNPs on the molecular phenotype of proteins. SNPeffect analyses the effect of SNPs on three categories of functional properties: (1) structural and thermodynamic properties affecting protein dynamics and stability (2) the integrity of functional and binding sites and (3) changes in posttranslational processing and cellular localization of proteins. The search interface of the database can be used to search specifically for polymorphisms that are predicted to cause a change in one of these properties. Now based on the Ensembl human databases, the SNPeffect database has been remodeled to better fit an automatically updatable structure. The current edition holds the molecular phenotype of 74 567 nsSNPs in 23 426 proteins. AVAILABILITY: SNPeffect can be accessed through http://snpeffect.vib.be.

Algorithms↗

Protein aggregation and amyloidosis: confusion of the kinds?

Recent years have witnessed major advances in our understanding of the structural basis of protein aggregation on several fronts. Firstly, high-resolution structural information that remained elusive for many years was provided by a series of studies of amyloid fibers using NMR, X-ray crystallography and electron microscopy, thereby confirming earlier models based on lower resolution observations. Secondly, studies of the sequence determinants of protein aggregation culminated in the development of computer algorithms that predict aggregation-prone sequences with good accuracy, allowing the design of mutations that reduce aggregation. Thirdly, based on the first results from such predictions and on statistical analysis of naturally occurring aggregating sequences, a picture is emerging in which aggregation-prone sequences are capped by gatekeeper residues that oppose aggregation. In addition to their aggregation-opposing function, it seems that gatekeeper residues are also important in determining chaperone selectivity for strongly aggregating regions. Finally, recent computational and experimental work shows that preventing aggregation does not necessarily mean that amyloid formation is prevented and vice versa. Thus, although aggregation and amyloidosis correlate to a certain extent, they are different processes and should be treated as such.

Algorithms↗

Antigen binding and solubility effects upon the veneering of a camel VHH in framework-2 to mimic a VH.

Heavy chain only antibodies of camelids bind their antigens with a single domain, the VHH, which acquired adaptations relative to classical VHs to function in the absence of a VL partner. Additional CDR loop conformations, outside the canonical loop structures of VHs, broaden the repertoire of the antigen-binding site. The combined effects of part of the CDR3 that folds over the "former" VL binding site and framework-2 mutations to more hydrophilic amino acids, enhance the solubility of VHH domains and prevent VL pairing. cAbAn33, a VHH domain specific for the carbohydrate moiety of the variant surface glycoprotein of trypanosomes, has a short CDR3 loop that does not cover the former VL binding site as well as a VH-specific Trp47 instead of the VHH-specific Gly47. Resurfacing its framework-2 region (mutations Tyr37Val, Glu44Gly and Arg45Leu) to mimic that of a human VH restores the VL binding capacity. In solution, the humanised VHH behaves as a soluble, monomeric entity, albeit with reduced thermodynamic stability and affinity for its antigen. Comparison of the crystal structures of cAbAn33 and its humanised derivative reveals steric hindrance exerted by VHH-specific residues Tyr37 and Arg45 that prevent the VL domain pairing, whereas Glu44 and Arg45 are key elements to avoid insolubility of the domain.

Amino Acid Sequence↗

The FoldX web server: an online force field.

FoldX is an empirical force field that was developed for the rapid evaluation of the effect of mutations on the stability, folding and dynamics of proteins and nucleic acids. The core functionality of FoldX, namely the calculation of the free energy of a macromolecule based on its high-resolution 3D structure, is now publicly available through a web server at http://foldx.embl.de/. The current release allows the calculation of the stability of a protein, calculation of the positions of the protons and the prediction of water bridges, prediction of metal binding sites and the analysis of the free energy of complex formation. Alanine scanning, the systematic truncation of side chains to alanine, is also included. In addition, some reporting functions have been added, and it is now possible to print both the atomic interaction networks that constitute the protein, print the structural and energetic details of the interactions per atom or per residue, as well as generate a general quality report of the pdb structure. This core functionality will be further extended as more FoldX applications are developed.

Databases, Protein↗

Recognizing and defining true Ras binding domains II: in silico prediction based on homology modelling and energy calculations.

Considering the large number of putative Ras effector proteins, it is highly desirable to develop computational methods to be able to identify true Ras binding molecules. Based on a limited sequence homology among members of the Ras association (RA) and Ras binding (RB) sub-domain families of the ubiquitin super-family, we have built structural homology models of Ras proteins in complex with different RA and RB domains, using the FOLD-X software. A critical step in our approach is to use different templates of Ras complexes, in order to account for the structural variation among the RA and RB domains. The homology models are validated by predicting the effect of mutating hot spot residues in the interface, and residues important for the specificity of interaction with different Ras proteins. The FOLD-X calculated energies of the best-modelled complexes are in good agreement with previously published experimental data and with new data reported here. Based on these results, we can establish energy thresholds above, or below which, we can predict with 96% confidence that a RA/RB domain will or will not interact with Ras. This study shows the importance of in depth structural analysis, high quality force-fields and modelling for correct prediction. Our work opens the possibility of genome-wide prediction for this protein family and for others, where there is enough structural information.

Amino Acid Sequence↗

SNPeffect: a database mapping molecular phenotypic effects of human non-synonymous coding SNPs.

Single nucleotide polymorphisms (SNPs) are an increasingly important tool for genetic and biomedical research. However, the accumulated sequence information on allelic variation is not matched by an understanding of the effect of SNPs on the functional attributes or 'molecular phenotype' of a protein. Towards this aim we developed SNPeffect, an online resource of human non-synonymous coding SNPs (nsSNPs) mapping phenotypic effects of allelic variation in human genes. SNPeffect contains 31 659 nsSNPs from 12 480 human proteins. The current release of SNPeffect incorporates data on protein stability, integrity of functional sites, protein phosphorylation and glycosylation, subcellular localization, protein turnover rates, protein aggregation, amyloidosis and chaperone interaction. The SNP entries are accessible through both a search and browse interface and are linked to most major biological databases. The data can be displayed as detailed descriptions of individual SNPs or as an overview of all SNPs for a given protein. SNPeffect will be regularly updated and can be accessed at http://snpeffect.vib.be/.

Alleles↗

A systems biology perspective on protein structural dynamics and signal transduction.

The functional dynamics of signal transduction through protein interaction networks are determined both by network topology and by the signal processing properties of component proteins. In order to understand the emergent properties of signal transduction networks in terms of information processing, storage and decision making, we not only need to map the so-called 'interactome' but, perhaps more importantly, we also have to understand how the structural dynamics of constituent proteins shape non-linear responses through cooperativity and allostery. Several in silico methods have been developed to identify networks of cooperative residues in proteins and help infer their mode of action. Applying this type of analysis to important classes of modular signal transduction domains should, in principle, allow the function of these proteins to be abstracted in terms of their information processing characteristics, permitting better comprehension of the systemic properties of biological networks.

Binding Sites↗

Prediction of sequence-dependent and mutational effects on the aggregation of peptides and proteins.

We have developed a statistical mechanics algorithm, TANGO, to predict protein aggregation. TANGO is based on the physico-chemical principles of beta-sheet formation, extended by the assumption that the core regions of an aggregate are fully buried. Our algorithm accurately predicts the aggregation of a data set of 179 peptides compiled from the literature as well as of a new set of 71 peptides derived from human disease-related proteins, including prion protein, lysozyme and beta2-microglobulin. TANGO also correctly predicts pathogenic as well as protective mutations of the Alzheimer beta-peptide, human lysozyme and transthyretin, and discriminates between beta-sheet propensity and aggregation. Our results confirm the model of intermolecular beta-sheet formation as a widespread underlying mechanism of protein aggregation. Furthermore, the algorithm opens the door to a fully automated, sequence-based design strategy to improve the aggregation properties of proteins of scientific or industrial interest.

Algorithms↗

A comparative study of the relationship between protein structure and beta-aggregation in globular and intrinsically disordered proteins.

A growing number of proteins are being identified that are biologically active though intrinsically disordered, in sharp contrast with the classic notion that proteins require a well-defined globular structure in order to be functional. At the same time recent work showed that aggregation and amyloidosis are initiated in amino acid sequences that have specific physico-chemical properties in terms of secondary structure propensities, hydrophobicity and charge. In intrinsically disordered proteins (IDPs) such sequences would be almost exclusively solvent-exposed and therefore cause serious solubility problems. Further, some IDPs such as the human prion protein, synuclein and Tau protein are related to major protein conformational diseases. However, this scenario contrasts with the large number of unstructured proteins identified, especially in higher eukaryotes, and the fact that the solubility of these proteins is often particularly good. We have used the algorithm TANGO to compare the beta aggregation tendency of a set of globular proteins derived from SCOP and a set of 296 experimentally verified, non-redundant IDPs but also with a set of IDPs predicted by the algorithms DisEMBL and GlobPlot. Our analysis shows that the beta-aggregation propensity of all-alpha, all-beta and mixed alpha/beta globular proteins as well as membrane-associated proteins is fairly similar. This illustrates firstly that globular structures possess an appreciable amount of structural frustration and secondly that beta-aggregation is not determined by hydrophobicity and beta-sheet propensity alone. We also show that globular proteins contain almost three times as much aggregation nucleating regions as IDPs and that the formation of highly structured globular proteins comes at the cost of a higher beta-aggregation propensity because both structure and aggregation obey very similar physico-chemical constraints. Finally, we discuss the fact that although IDPs have a much lower aggregation propensity than globular proteins, this does not necessarily mean that they have a lower potential for amyloidosis.

Algorithms↗