First-Principle Determination of Peptide Conformations in Solvents: Combination of Monte Carlo Simulated Annealing and RISM Theory

Abstract

This paper contributes to development of a microscopic approach to predicting stable conformations of proteins in solvent. We report results of the first attempt to combine Monte Carlo simulated annealing, a powerful conformational sampling technique, and the reference interaction site model (RISM) theory, a statistical-mechanical treatment for molecular fluids. In solvent the key function is the total energy defined as the sum of the conformational energy and the solvation free energy, and the RISM theory is employed to calculate the latter. Starting from an initial conformation given, our computer program samples many conformations and then finds the conformation with the minimum total energy. Met-enkephalin in the two different solvents, a model water and a simple, repulsive-potential system, are considered. In water the solvation free energy varies greatly from conformation to conformation, while in the simple solvent it remains almost unchanged against conformational changes. In water most of the con...