Abstract
We present a new method for the calculation of solvent accessible surface areas at the atomic and residue levels, which we call parameter optimized surfaces (POPS-A and POPS-R ). Atomic and residue areas (the latter simulated with a single sphere centered at the C(alpha)s atom for amino acids and at the P atom for nucleotides) have been optimized versus accurate all-atoms methods. We concentrated on an analytical formula for the approximation of solvent accessibilities. The formula is simple, easily derivable and fast to compute, therefore it is practical for use in molecular dynamics simulations as an approximation to the first solvation shell. The residue based approach POPS-R has been derived as a useful tool for the analysis of large macromolecular assemblies like the ribosome, and is especially suited for use in refinement of low resolution structures. The structures of the 70S, 50S and 30S ribosomes have been analyzed in detail and most of the interactions within the subunits and at their interfaces were clearly identified. Some interesting differences between 30S alone and within the 70S have been highlighted. Owing to the presence of the P-tRNA in the 70S ribosome, localized conformational rearrangements occur within the subunits, exposing Arg and Lys residues to negatively charged binding sites of P-tRNA. POPS-R also allows for estimates of the loss of free energy of solvation upon complex formation, particularly useful in designing new protein-RNA complexes and in suggesting more focused experimental work.