Atomically detailed folding simulation of the B domain of staphylococcal protein A from random structures

Abstract

The conformational space of the 10-55 fragment of the B-domain of staphylococcal protein A has been investigated by using the electrostatically driven Monte Carlo (EDMC) method. The ECEPP/3 (empirical conformational energy program for peptides) force-field plus two different continuum solvation models, namely SRFOPT (Solvent Radii Fixed with atomic solvation parameters OPTimized) and OONS (Ooi, Oobatake, Némethy, and Scheraga solvation model), were used to describe the conformational energy of the chain. After an exhaustive search, starting from two different random conformations, three of four runs led to native-like conformations. Boltzmann-averaged root-mean-square deviations (RMSD) for all of the backbone heavy atoms with respect to the native structure of 3.35 A and 4.54 A were obtained with SRFOPT and OONS, respectively. These results show that the protein-folding problem can be solved at the atomic detail level by an ab initio procedure, starting from random conformations, with no knowledge except the amino acid sequence. To our knowledge, the results reported here correspond to the largest protein ever folded from a random conformation by an initial-value formulation with a full atomic potential, without resort to knowledge-based information.