A HNP LATTICE MODEL OF PROTEINS

XH Wang,LX Zhang,DL Zhao

DOI: https://doi.org/10.3321/j.issn:1000-3304.2004.02.023

2004-01-01

Abstract:Effective interresidue contact energies for proteins in solutions are estimated from the number of residue-residue contacts observed in crystal structures of globular proteins ( http://www.rcsb.org/) by means of the Miyazawa-Jerniganm (MJ) interaction matrix. It is found that residues can be classified into three groups according to the contact energies: hydrophobic residues (H), neutral residues (N), and hydrophilic residues (P). Here we develop a modified HP lattice model, i.e. a HNP lattice model. In this model, the contact Hamiltonian is simply written asH = Sigma(i<j)E(Si-Sj)delta(r(ij), l) (S-i = H, or N, or P) (4)Where EH-H = -5.59RT, EH-N = -3.68RT, EH-P = -3.07RT, EN-N = 2.23RT, EN-P = -1.78RT, EP-P = -1.37RT, and r(ij) is the distance between residues i and j. We first calculate the conformational properties of the native states of short protein chains using 2-dimensional HNP lattice model and enumeration calculation method. The distributions of various types of residue-residue contacts and three types of residues of the native states are discussed. We also investigate the folding of the protein with a sequence of HHNHNPNHPPHPNPPHPHPPHHPHNH using the modified HP lattice model. Here the sequence is selected at random. The folded conformation is obtained in square lattice, and the folding energy is - 64.89 RT. Our calculation shows that the HNP lattice model may be more reasonable for simulation of the folding behavior of real proteins.