PubMed Health⌕ Search

Biomedical subjects

Silvio C E Tosatto

Publications and source records attributed to Silvio C E Tosatto.

18 recordsLinked to original sources

Open and sustainable AI: challenges, opportunities and the road ahead in the life sciences.

Artificial intelligence (AI) has seen transformative breakthroughs in the life sciences, expanding possibilities to interpret biological information at an unprecedented capacity. To maximize return on growing investments and accelerate progress, it is urgent to address long-standing research challenges arising from the rapid adoption of AI methods. We review the erosion of trust in AI outputs driven by poor reusability and reproducibility, and highlight their impact on environmental sustainability. Furthermore, we discuss the fragmented components of the AI ecosystem and lack of guiding pathways to support open and sustainable AI model development. In response, this Perspective introduces practical open and sustainable AI recommendations mapped to over 300 ecosystem components and provides guiding implementation pathways. Our work connects researchers with relevant AI resources, facilitating the implementation of sustainable, reusable and reproducible AI. Built upon community consensus and aligned to existing efforts, these outputs will aid future policy development and structured pathways for guiding AI implementation.

Artificial Intelligence↗

GeomeTRe: accurate calculation of geometrical descriptors of tandem repeat proteins.

MOTIVATION: Structured tandem repeat proteins (STRPs) are characterized by preserved structural motifs arranged in a modular way. The structural and functional diversity of STRPs makes them particularly important for studying evolution and novel structure-function relationships, and ultimately for designing new synthetic proteins with specific functions. One crucial aspect of their classification is the estimation of geometrical parameters, which can provide better insight into their properties and the relationship between the spatial arrangement of repeated units and protein function. Calculating geometric descriptors for STRPs is challenging because naturally occurring repeats are not "perfect" and often contain insertions and deletions. Existing tools for predicting structural symmetry work well on simple cases but often fail for most natural proteins. RESULTS: Here, we present GeomeTRe, an algorithm that calculates geometrical descriptors such as curvature (yaw), twist (roll), and pitch for a protein structure with known repeat unit positions. The algorithm simulates the movement of consecutive units, identifies rotational axes, and calculates the corresponding Tait-Bryan angles. GeomeTRe's parameters can enhance STRP annotation and classification by identifying variations in geometric arrangements among different functional groups. The package is fast and suitable for processing large protein structure datasets when repeat region information (e.g. from RepeatsDB) is available. AVAILABILITY AND IMPLEMENTATION: GeomeTRe is available as a Python package; source code and documentation can be found at https://github.com/BioComputingUP/GeomeTRe.

Algorithms↗

MobiDB-lite 4.0: faster prediction of intrinsic protein disorder and structural compactness.

MOTIVATION: In recent years, many disorder predictors have been developed to identify intrinsically disordered regions (IDRs) in proteins, achieving high accuracy. However, it may be difficult to interpret differences in predictions across methods. Consensus methods offer a simple solution, highlighting reliable predictions while filtering out uncertain positions. Here, we present a new version of MobiDB-lite, a consensus method designed to predict long IDRs and classify them based on compositional biases and conformational properties. RESULTS: MobiDB-lite 4.0 pipeline was optimized to be ten times faster than the previous version. It now provides compactness annotations based on predicted apparent scaling exponent. The newly added features and disorder subclassifications allow the users to get a comprehensive insight into the protein's function and characteristics. MobiDB-lite 4.0 is integrated into the MobiDB and DisProt databases. A version without the compactness predictor is integrated into InterProScan, propagating MobiDB-lite annotations to UniProtKB. AVAILABILITY AND IMPLEMENTATION: The MobiDB-lite 4.0 source code and a Docker container are available from the GitHub repository: https://github.com/BioComputingUP/MobiDB-lite.

Intrinsically Disordered Proteins↗

The thioredoxin specificity of Drosophila GPx: a paradigm for a peroxiredoxin-like mechanism of many glutathione peroxidases.

Some members of the glutathione peroxidase (GPx) family have been reported to accept thioredoxin as reducing substrate. However, the selenocysteine-containing ones oxidise thioredoxin (Trx), if at all, at extremely slow rates. In contrast, the Cys homolog of Drosophila melanogaster exhibits a clear preference for Trx, the net forward rate constant, k'(+2), for reduction by Trx being 1.5x10(6) M(-1) s(-1), but only 5.4 M(-1) s(-1) for glutathione. Like other CysGPxs with thioredoxin peroxidase activity, Drosophila melanogaster (Dm)GPx oxidized by H(2)O(2) contained an intra-molecular disulfide bridge between the active-site cysteine (C45; C(P)) and C91. Site-directed mutagenesis of C91 in DmGPx abrogated Trx peroxidase activity, but increased the rate constant for glutathione by two orders of magnitude. In contrast, a replacement of C74 by Ser or Ala only marginally affected activity and specificity of DmGPx. Furthermore, LC-MS/MS analysis of oxidized DmGPx exposed to a reduced Trx C35S mutant yielded a dead-end intermediate containing a disulfide between Trx C32 and DmGPx C91. Thus, the catalytic mechanism of DmGPx, unlike that of selenocysteine (Sec)GPxs, involves formation of an internal disulfide that is pivotal to the interaction with Trx. Hereby C91, like the analogous second cysteine in 2-cysteine peroxiredoxins, adopts the role of a "resolving" cysteine (C(R)). Molecular modeling and homology considerations based on 450 GPxs suggest peculiar features to determine Trx specificity: (i) a non-aligned second Cys within the fourth helix that acts as C(R); (ii) deletions of the subunit interfaces typical of tetrameric GPxs leading to flexibility of the C(R)-containing loop. Based of these characteristics, most of the non-mammalian CysGPxs, in functional terms, are thioredoxin peroxidases.

Amino Acid Sequence↗

Improving the quality of protein structure models by selecting from alignment alternatives.

BACKGROUND: In the area of protein structure prediction, recently a lot of effort has gone into the development of Model Quality Assessment Programs (MQAPs). MQAPs distinguish high quality protein structure models from inferior models. Here, we propose a new method to use an MQAP to improve the quality of models. With a given target sequence and template structure, we construct a number of different alignments and corresponding models for the sequence. The quality of these models is scored with an MQAP and used to choose the most promising model. An SVM-based selection scheme is suggested for combining MQAP partial potentials, in order to optimize for improved model selection. RESULTS: The approach has been tested on a representative set of proteins. The ability of the method to improve models was validated by comparing the MQAP-selected structures to the native structures with the model quality evaluation program TM-score. Using the SVM-based model selection, a significant increase in model quality is obtained (as shown with a Wilcoxon signed rank test yielding p-values below 10(-15)). The average increase in TMscore is 0.016, the maximum observed increase in TM-score is 0.29. CONCLUSION: In template-based protein structure prediction alignment is known to be a bottleneck limiting the overall model quality. Here we show that a combination of systematic alignment variation and modern model scoring functions can significantly improve the quality of alignment-based models.

Computer Simulation↗

Spritz: a server for the prediction of intrinsically disordered regions in protein sequences using kernel machines.

Intrinsically disordered proteins have long stretches of their polypeptide chain, which do not adopt a single native structure composed of stable secondary and tertiary structure in the absence of binding partners. The prediction of intrinsically disordered regions in proteins from sequence is increasingly becoming of interest, as the presence of many such regions in the complete genome sequences are discovered and important functional roles are associated with them. We have developed a machine learning approach based on two support vector machines (SVM) to discriminate disordered regions from sequence. The SVM are trained and benchmarked on two sets, representing long and short disordered regions. A preliminary version of Spritz was shown to perform consistently well at the recent biannual CASP-6 experiment [Critical Assessment of Techniques for Protein Structure Prediction (CASP), 2004]. The fully developed Spritz method is freely available as a web server at http://distill.ucd.ie/spritz/ and http://protein.cribi.unipd.it/spritz/.

Artificial Intelligence↗

Fine-grained statistical torsion angle potentials are effective in discriminating native protein structures.

Modelling of drug targets requires the reliable selection of an accurate and representative structure from large ensembles of alternate models. Statistical potentials developed to discriminate native protein structures generally represent pairwise interactions between atoms, which are less sensitive to local conformational details. The discrimination of local distortions is therefore particularly difficult. Local interaction preferences, expressed through torsion angles, are rarely used, as some controversy exists in the literature regarding their discrimination power. The present study aims to benchmark the efficiency of different implementations of torsion angle propensities for selecting the native structure from ensembles of well-constructed decoys. Several statistical potentials derived from fine-grained discretisations of torsion angle space are constructed and evaluated. Results from a comparison with nine widely used statistical scoring functions show the torsion angle potentials to be more effective in recognising native structures and to improve with the number of torsion angles considered. These data suggest local structural propensities to be important for estimating the overall quality of native-like models.

Models, Chemical↗

A decoy set for the thermostable subdomain from chicken villin headpiece, comparison of different free energy estimators.

BACKGROUND: Estimators of free energies are routinely used to judge the quality of protein structural models. As these estimators still present inaccuracies, they are frequently evaluated by discriminating native or native-like conformations from large ensembles of so-called decoy structures. RESULTS: A decoy set is obtained from snapshots taken from 5 long (100 ns) molecular dynamics (MD) simulations of the thermostable subdomain from chicken villin headpiece. An evaluation of the energy of the decoys is given using: i) a residue based contact potential supplemented by a term for the quality of dihedral angles; ii) a recently introduced combination of four statistical scoring functions for model quality estimation (FRST); iii) molecular mechanics with solvation energy estimated either according to the generalized Born surface area (GBSA) or iv) the Poisson-Boltzmann surface area (PBSA) method. CONCLUSION: The decoy set presented here has the following features which make it attractive for testing energy scoring functions:1) it covers a broad range of RMSD values (from less than 2.0 A to more than 12 A);2) it has been obtained from molecular dynamics trajectories, starting from different non-native-like conformations which have diverse behaviour, with secondary structure elements correctly or incorrectly formed, and in one case folding to a native-like structure. This allows not only for scoring of static structures, but also for studying, using free energy estimators, the kinetics of folding;3) all structures have been obtained from accurate MD simulations in explicit solvent and after molecular mechanics (MM) energy minimization using an implicit solvent method. The quality of the covalent structure therefore does not suffer from steric or covalent problems. The statistical and physical effective energy functions tested on the set behave differently when native simulation snapshots are included or not in the set and when averaging over the trajectory is performed.

Animals↗

Functional insights from the structural modelling of a small Fe-hydrogenase.

Recently, a novel Fe-hydrogenase from a high rate of hydrogen producing Enterobacter cloacae strain IIT-BT08 was identified and partially characterized. This 147 residue protein was found to be much smaller than previously known Fe-hydrogenases, yet retaining a high catalytic activity. We predicted the structure of this protein and found it to be structurally similar to one of the two sub-domains containing the catalytic H-cluster so far jointly present in all other Fe-hydrogenases. This novel architecture allows a tentative explanation of protein function with the high rate of catalytic activity being due to a missing regulatory sub-domain, presumably allowing higher enzymatic activity at the cost of greater exposure to oxygen inactivation. This new insight may improve our understanding of the molecular and functional organization of other, more complex Fe-hydrogenases.

Amino Acid Sequence↗

Structural insights into the function of human caveolin 1.

Caveolin-1 (Cav-1) is emerging as the central protein controlling caveolae formation, caveolae trafficking, and cellular signalling. In particular, it is known that Cav-1 interacts and modulates the activity of several signalling proteins through the so-called caveolin scaffolding domain. In this paper, we used a bioinformatics approach to assess the validity of some long-standing structural features of Cav-1. We could confirm the existence of a membrane spanning region of Cav-1 and highlight an interesting pattern of palmitoylated cysteine residues explaining the structural features of the Cav-1 C-terminal region. Moreover, the scaffolding domain is predicted to have a different structure than previously reported.

Amino Acid Sequence↗

Functional interaction of phospholipid hydroperoxide glutathione peroxidase with sperm mitochondrion-associated cysteine-rich protein discloses the adjacent cysteine motif as a new substrate of the selenoperoxidase.

The mitochondrial capsule is a selenium- and disulfide-rich structure enchasing the outer mitochondrial membrane of mammalian spermatozoa. Among the proteins solubilized from the sperm mitochondrial capsule, we confirmed, by using a proteomic approach, the presence of phospholipid hydroperoxide glutathione peroxidase (PHGPx) as a major component, and we also identified the sperm mitochondrion-associated cysteine-rich protein (SMCP) and fragments/aggregates of specific keratins that previously escaped detection (Ursini, F., Heim, S., Kiess, M., Maiorino, M., Roveri, A., Wissing, J., and Flohé, L. (1999) Science 285, 1393-1396). The evidence for a functional association between PHGPx, SMCP, and keratins is further supported by the identification of a sequence motif of regularly spaced Cys-Cys doublets common to SMCP and high sulfur keratin-associated proteins, involved in bundling hair shaft keratin by disulfide cross-linking. Following the oxidative polymerization of mitochondrial capsule proteins, catalyzed by PHGPx, two-dimensional redox electrophoresis analysis showed homo- and heteropolymers of SMCP and PHGPx, together with other minor components. Adjacent cysteine residues in SMCP peptides are oxidized to cystine by PHGPx. This unusual disulfide is known to drive, by reshuffling oxidative protein folding. On this basis we propose that oxidative polymerization of the mitochondrial capsule is primed by the formation of cystine on SMCP, followed by reshuffling. Occurrence of reshuffling is further supported by the calculated thermodynamic gain of the process. This study suggests a new mechanism where selenium catalysis drives the cross-linking of structural elements of the cytoskeleton via the oxidation of a keratin-associated protein.

Amino Acid Motifs↗

Decomposing protein networks into domain-domain interactions.

UNLABELLED: The application of novel experimental techniques has generated large networks of protein-protein interactions. Frequently, important information on the structure and cellular function of protein-protein interactions can be gained from the domains of interacting proteins. We have designed a Cytoscape plugin that decomposes interacting proteins into their respective domains and computes a putative network of corresponding domain-domain interactions. To this end, the network graph of proteins has been extended by additional node and edge types for domain interactions, including different node and edge shapes and coloring schemes used for visualization. An additional plugin provides supplementary web links to Internet resources on domain function and structure. AVAILABILITY: Both Cytoscape plugins can be downloaded from http://www.cytoscape.org

Algorithms↗

The victor/FRST function for model quality estimation.

Scoring functions are widely used in the final step of model selection in protein structure prediction. This is of interest both for comparative modeling targets, where it is important to select the best model among a set of many good, "correct" ones, as well as for other (fold recognition or novel fold) targets, where the set may contain many incorrect models. A novel combination of four knowledge-based potentials recognizing different features of native protein structures is introduced and tested. The pairwise, solvation, hydrogen bond, and torsion angle potentials contain largely orthogonal information. Of these, the torsion angle potential is found to show the strongest correlation with model quality. Combining these features with a linear weighting function, it was possible to construct a robust energy function capable of discriminating native-like structures on several benchmarking sets. In a recent blind test (CAFASP-4 MQAP), the scoring function ranked consistently well and was able to reliably distinguish the correct template from an ensemble of high quality decoys in 52 of 70 cases (33 of 34 for comparative modeling). An executable version of the Victor/FRST function for Linux PCs is available for download from the URL http://protein.cribi.unipd.it/frst/.

Computational Biology↗

Application of MM/PBSA colony free energy to loop decoy discrimination: toward correlation between energy and root mean square deviation.

Accurate free energy estimation is needed in many predictive tasks. The molecular mechanics/Poisson-Boltzmann solvent accessible surface area (MM/PBSA) approach has proven to be accurate. However, the correlation between the estimated free energy and the distance (e.g., root mean square deviation [RMSD]) from the most stable conformation is hindered by the strong free energy dependence on minor conformational variations. In this paper, a protocol for MM/PBSA free energy estimation is designed and tested on several loop decoy sets. We show that further integration of MM/PBSA free energy estimator with the colony energy approach makes the correlation between the free energy and RMSD from the native structure apparent, for the test sets on which it could be applied. Our results suggest that (1) the MM/PBSA free energy estimator is able to detect native-like structures for most decoy sets, and (2) application of the colony energy approach greatly hampers the MM/energy strong dependence on minor conformational changes.

Amino Acids↗

The SSEA server for protein secondary structure alignment.

SUMMARY: We present a web server that computes alignments of protein secondary structures. The server supports both performing pairwise alignments and searching a secondary structure against a library of domain folds. It can calculate global and local secondary structure element alignments. A combination of local and global alignment steps can be used to search for domains inside the query sequence or help in the discrimination of novel folds. Both the SCOP and PDB fold libraries, clustered at 95 and 40% sequence identity, are available for alignment. AVAILABILITY: The web server interface is freely accessible to academic users at http://protein.cribi.unipd.it/ssea/. The executable version and benchmarking data are available from the same web page.

Algorithms↗

Simple consensus procedures are effective and sufficient in secondary structure prediction.

We have analyzed the performance of majority voting on minimal combination sets of three state-of-the-art secondary structure prediction methods in order to obtain a consensus prediction. Using three large benchmark sets from the EVA server, our results show a significant improvement in the average Q3 prediction accuracy of up to 1.5 percentage points by consensus formation. The application of an additional trivial filtering procedure for predicted secondary structure elements that are too short, does not significantly affect the prediction accuracy. Our analysis also provides valuable insight into the similarity of the results of the prediction methods that we combine as well as the higher confidence in consistently predicted secondary structure.

Computational Biology↗

MANIFOLD: protein fold recognition based on secondary structure, sequence similarity and enzyme classification.

We present a protein fold recognition method, MANIFOLD, which uses the similarity between target and template proteins in predicted secondary structure, sequence and enzyme code to predict the fold of the target protein. We developed a non-linear ranking scheme in order to combine the scores of the three different similarity measures used. For a difficult test set of proteins with very little sequence similarity, the program predicts the fold class correctly in 34% of cases. This is an over twofold increase in accuracy compared with sequence-based methods such as PSI-BLAST or GenTHREADER, which score 13-14% correct first hits for the same test set. The functional similarity term increases the prediction accuracy by up to 3% compared with using the combination of secondary structure similarity and PSI-BLAST alone. We argue that using functional and secondary structure information can increase the fold recognition beyond sequence similarity.

Algorithms↗

A divide and conquer approach to fast loop modeling.

We describe a fast ab initio method for modeling local segments in protein structures. The algorithm is based on a divide and conquer approach and uses a database of precalculated look-up tables, which represent a large set of possible conformations for loop segments of variable length. The target loop is recursively decomposed until the resulting conformations are small enough to be compiled analytically. The algorithm, which is not restricted to any specific loop length, generates a ranked set of loop conformations in 20-180 s on a desktop PC. The prediction quality is evaluated in terms of global RMSD. Depending on loop length the top prediction varies between 1.06 A RMSD for three-residue loops and 3.72 A RMSD for eight-residue loops. Due to its speed the method may also be useful to generate alternative starting conformations for complex simulations.

Algorithms↗