MONTE CARLO SIMULATION ON AQUEOUS SOLUTION OF BIOLOGICAL MOLECULES——ALANINE SOLUTION

The Monte Carlo simulation is performed for a cluster consisting of a neutral alanine molecule surrounded by 56 water molecules. The average water-water and alanine-water interaction energies are found. The space surrounding the alanine molecule is divided into three regions, where the central atoms are C’ atom for the earboxyl group region, N atom for the amino group region and Catom for the methyl group region. The water-water and water-alanine interaction energies as functions of the distance between the oxygen atom in a water molecule and the central atom in each region are calculated. In each region the orientational correlation function for the dipole moments of water molecules, the radial distribution function for the oxygen and hydrogen atoms of water molecules are evaluated. In addition, the numbers of water molecules in the first solvation shells of each region and of entire alanine molecule are also counted.