Comparison of electrostatic potential around proteins calculated from Amber and AM1 charges: application to mutants of prion protein

On the basis of arguments of complementary fit of shape and charge polarity or hydrophobicity, molecular electrostatic potentials (MEPs) around proteins are commonly used to deduce likely sites for interaction with ligands or other proteins, including for variations such as mutations. But protein MEPs calculated classically from fixed force field descriptions, including those with implicit solvent models such as in Delphi, do not allow for repolarization of protein residues within the protein system; hence, their representations are likely to be variably inaccurate. Linear-scaling methods now allow calculation of MEPs quantum mechanically for systems as large as proteins, and can account for polarization explicitly. Here we compare MEPs derived from AM1 charge distributions calculated by Mopac2000 with those from the classical Amber force field. Our models are mutants of prion protein (PrP), a protein with an unusually high number of charged residues. The results demonstrate that static point charges, as used in most current force fields, cannot reproduce the MEP of macromolecules. Also, it is not sufficient to account for the influence of nearby atoms connected by chemical bonds; the influence of nearby atoms in space is at least as important. Thus, further progress in the accuracy and wider applicability of force fields requires proper accounting for polarization. Mopac2000 calculations can provide the necessary data for checking new force fields and/or parameter fitting.

[1]  S. Prusiner,et al.  Glycosylation differences between the normal and pathogenic prion protein isoforms. , 1999, Proceedings of the National Academy of Sciences of the United States of America.

[2]  H. Berendsen,et al.  Molecular dynamics with coupling to an external bath , 1984 .

[3]  W. L. Jorgensen,et al.  Comparison of simple potential functions for simulating liquid water , 1983 .

[4]  W. Kabsch,et al.  Dictionary of protein secondary structure: Pattern recognition of hydrogen‐bonded and geometrical features , 1983, Biopolymers.

[5]  B. Roux,et al.  Implicit solvent models. , 1999, Biophysical chemistry.

[6]  Alistair P. Rendell,et al.  The Potassium Ion Channel: Comparison of Linear Scaling Semiempirical and Molecular Mechanics Representations of the Electrostatic Potential , 2001 .

[7]  Peter A. Kollman,et al.  AMBER, a package of computer programs for applying molecular mechanics, normal mode analysis, molecular dynamics and free energy calculations to simulate the structural and energetic properties of molecules , 1995 .

[8]  K Wüthrich,et al.  Prion protein NMR structure and familial human spongiform encephalopathies. , 1998, Proceedings of the National Academy of Sciences of the United States of America.

[9]  Jacopo Tomasi,et al.  Molecular Interactions in Solution: An Overview of Methods Based on Continuous Distributions of the Solvent , 1994 .

[10]  James J. P. Stewart,et al.  Application of localized molecular orbitals to the solution of semiempirical self‐consistent field equations , 1996 .

[11]  A. Burlingame,et al.  Site-specific characterization of the N-linked glycans of murine prion protein by high-performance liquid chromatography/electrospray mass spectrometry and exoglycosidase digestions. , 1999, Biochemistry.

[12]  K Wüthrich,et al.  NMR solution structure of the human prion protein. , 2000, Proceedings of the National Academy of Sciences of the United States of America.

[13]  S. Petersen,et al.  Surface and electrostatics of cutinases. , 1997, Methods in enzymology.

[14]  J. Gready,et al.  Molecular dynamics simulations of human prion protein: importance of correct treatment of electrostatic interactions. , 1999, Biochemistry.

[15]  P. Kollman,et al.  An all atom force field for simulations of proteins and nucleic acids , 1986, Journal of computational chemistry.

[16]  K Wüthrich,et al.  NMR structure of the bovine prion protein. , 2000, Proceedings of the National Academy of Sciences of the United States of America.

[17]  T. Simonson,et al.  Macromolecular electrostatics: continuum models and their growing pains. , 2001, Current opinion in structural biology.

[18]  S. Prusiner.  Novel proteinaceous infectious particles cause scrapie. , 1982, Science.

[19]  A L Burlingame,et al.  Glycosylinositol phospholipid anchors of the scrapie and cellular prion proteins contain sialic acid. , 1992, Biochemistry.

[20]  F. J. Luque,et al.  Theoretical Methods for the Description of the Solvent Effect in Biomolecular Systems. , 2000, Chemical reviews.

[21]  Alistair P. Rendell,et al.  Comparison of linear‐scaling semiempirical methods and combined quantum mechanical/molecular mechanical methods for enzymic reactions. II. An energy decomposition analysis , 2002, J. Comput. Chem..

[22]  Jonathan W. Essex,et al.  A review of protein-small molecule docking methods , 2002, J. Comput. Aided Mol. Des..

[23]  T. Halgren,et al.  Polarizable force fields. , 2001, Current opinion in structural biology.

[24]  Ruth Nussinov,et al.  Principles of docking: An overview of search algorithms and a guide to scoring functions , 2002, Proteins.

[25]  F E Cohen,et al.  Pathologic conformations of prion proteins. , 1998, Annual review of biochemistry.

[26]  T Darden,et al.  New tricks for modelers from the crystallography toolkit: the particle mesh Ewald algorithm and its use in nucleic acid simulations. , 1999, Structure.

[27]  Alistair P. Rendell,et al.  Comparison of semiempirical and ab initio QM decomposition analyses for the interaction energy between molecules , 2002 .

[28]  F. Cohen,et al.  Solution structure of Syrian hamster prion protein rPrP(90-231). , 1999, Biochemistry.

[29]  Alistair P. Rendell,et al.  Comparison of linear-scaling semiempirical methods and combined quantum mechanical/molecular mechanical methods applied to enzyme reactions , 2000 .

[30]  P. Kollman,et al.  A Second Generation Force Field for the Simulation of Proteins, Nucleic Acids, and Organic Molecules , 1995 .

[31]  K Wüthrich,et al.  NMR structures of three single-residue variants of the human prion protein. , 2000, Proceedings of the National Academy of Sciences of the United States of America.

[32]  Peng Wang,et al.  Ab initio calculations on phosphorus compounds. II. Effects of disubstitution on ligand apicophilicity in phosphoranes , 1993, J. Comput. Chem..

[33]  Robert J. Woods,et al.  Molecular Mechanical and Molecular Dynamic Simulations of Glycoproteins and Oligosaccharides. 1. GLYCAM_93 Parameter Development , 1995 .

[34]  K. Sharp,et al.  Calculating the electrostatic potential of molecules in solution: Method and error assessment , 1988 .

[35]  J. Gready,et al.  Molecular dynamics simulation of human prion protein including both N-linked oligosaccharides and the GPI anchor. , 2000, Glycobiology.

[36]  P. Kollman,et al.  Application of RESP charges to calculate conformational energies, hydrogen bond energies, and free energies of solvation , 1993 .

[37]  J. Gready,et al.  Investigating Enzyme Reaction Mechanisms with Quantum Mechanical-Molecular Mechanical Plus Molecular Dynamics Calculations , 1998 .

[38]  B. Caughey,et al.  Specific binding of normal prion protein to the scrapie form via a localized domain initiates its conversion to the protease‐resistant state , 1999, The EMBO journal.