Prediction of protein conformation in water and on surfaces by Monte Carlo simulations using united-atom method

Abstract

The united-atom method has been used to model an avian pancreatic polypeptide (APP) in water and the adsorption process of an albumin subdomain (AS) onto graphite surface to observe the capability of this lumped modelling approach to generate structures observed in protein data bank (PDB) and from atomistic modelling. The subdomain structure of a protein is simplified by the united-atom approximation where the side chains and peptide groups are represented by lumped spheres. The total potential energy of the adsorption process involves the interaction between these lumped spheres by means of virtual bond chain interaction and the interaction of the spheres with the graphite surface by means of Lennard-Jones potential. The protein/polypeptide structure has been perturbed by Monte Carlo with energy minimisation to obtain the global minimum. Results on the APP in water showed a near-to-experimental PDB conformation revealing the two α-helix structures of this small protein molecule with the root mean square deviation among carbon backbone atoms of 5.9 Å. Protein adsorption on biosurfaces has been made by modelling AS, which has 60 amino acids. The surface is graphite, which is characterised by its hydrophobicity. Graphite was chosen because of its widely used applications in certain implants that interact with blood. Our simulation results showed final conformation close to that obtained by atomistic modelling. It also proved that the whole pattern of intramolecular hydrogen bonds was distorted. The model also demonstrated the random conformation of the original α-helix secondary structures of AS consistent with experimental and atomistic results. While atomistic simulation works well for simulating individual small proteins, the united-atom model is more efficient when simulating macromolecular and multiple protein adsorption where time and limiting computer capacity are key factors.