Tandem mass spectrometry-based proteomics is currently in great demand of computational methods that facilitate the elimination of likely false positives in peptide and protein identification. In the last few years, a number of new peptide identification programs have been described, but scores or other significance measures reported by these programs cannot always be directly translated into an easy to interpret error rate measurement such as the False Discovery Rate (FDR). In this work we used Generalized Lambda Distributions (GLDs) to model frequency distributions of database search scores computed by MASCOT, X!TANDEM with k-score plugin, OMSSA and InsPecT. From these distributions, we could successfully estimate p-values and false discovery rates with high accuracy. From the set of peptide assignments reported by any of these engines, we also defined a generic protein scoring scheme that enabled accurate estimation of protein-level p-values by simulation of random score distributions, which was also found to yield good estimates of protein-level FDR. The performance of these methods was evaluated by searching four freely available datasets ranging from 40,000 to 285,000 MS/MS spectra.