A computational model for the quantitative prediction of protein thermostability has been developed by means of the Volsurf method. A data set of 22 enzymes of reported thermostability in water systems, for the most part coming from thermophilic and hyperthermophilic organisms, has been built up. Molecular descriptors of the protein surface have been calculated and their role in the stabilization of the macromolecule has been analyzed by a multivariate statistical approach. The resulting regression model has shown a good predictivity and it has been able to quantitatively identify some structural requirements correlated with protein stability. The method can be the basis for a new computational support tool in rational protein design, which is complementary to the existing methods based on the sequence analysis.