An ab initio distributed multipole study of the electrostatic potential around an undecapeptide cyclosporin derivative and a comparison with point charge electrostatic models

An SCF calculation has been performed on C63H113N11O12, a derivative of the immuno‐suppressive drug cyclosporin, using a 3‐21G basis set and a Direct SCF method. A distributed multipole analysis has been performed on the resulting charge density to give a set of multipoles at each atomic site, which are used to calculate the electrostatic potential around the molecule. The potential maxima and minima on the accessible surface of the molecule are compared with those predicted using the corresponding Mulliken charges, and also using a potential‐derived point‐charge model based on the force‐field of Kollman et al. The Mulliken charges give a misleading picture of the electrostatic potential around this peptide. The potential‐derived charges give results which are in far better agreement with the ab initio distributed multipole model, despite being derived from calculations on smaller molecules with different basis sets and geometries. The limitations of point‐charge models for describing the electrostatic interactions of polypeptides are discussed.

[1]  R. S. Mulliken.  Electronic Population Analysis on LCAO–MO Molecular Wave Functions. I , 1955 .

[2]  Armin Widmer,et al.  Peptide conformations. Part 31. The conformation of cyclosporin a in the crystal and in solution , 1985 .

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

[4]  J. Almlöf,et al.  Principles for a direct SCF approach to LICAO–MOab‐initio calculations , 1982 .

[5]  R. Bader,et al.  Quantum Theory of Atoms in Molecules–Dalton Revisited , 1981 .

[6]  H. Weber,et al.  Crystal and molecular structure of an iodo-derivative of the cyclic undecapeptide cyclosporin A. , 1976, Helvetica chimica acta.

[7]  P. Claverie,et al.  The exact multicenter multipolar part of a molecular charge distribution and its simplified representations , 1988 .

[8]  Warren J. Hehre,et al.  Computation of electron repulsion integrals involving contracted Gaussian basis functions , 1978 .

[9]  Mark A. Spackman,et al.  A SIMPLE QUANTITATIVE MODEL OF HYDROGEN-BONDING , 1986 .

[10]  Donald E. Williams,et al.  Lone-pair electronic effects on the calculated ab initio SCF-MO electric potential and the crystal structures of azabenzenes , 1983 .

[11]  Frank A. Momany,et al.  Determination of partial atomic charges from ab initio molecular electrostatic potentials. Application to formamide, methanol, and formic acid , 1978 .

[12]  Sarah L. Price,et al.  The electrostatic interactions in van der Waals complexes involving aromatic molecules , 1987 .

[13]  W. A. Sokalski,et al.  Correlated molecular and cumulative atomic multipole moments , 1987 .

[14]  R. Bonaccorsi,et al.  An approximate expression of the electrostatic molecular potential in terms of completely transferable group contributions , 1977 .

[15]  M. Paniagua,et al.  1.1.1 electrostatic description of molecular systems , 1985 .

[16]  Donald E. Williams,et al.  Representation of the molecular electrostatic potential by a net atomic charge model , 1981 .

[17]  Raymond J. Abraham,et al.  Charge calculations in molecular mechanics. III: Amino acids and peptides , 1985 .

[18]  P. Pulay.  Improved SCF convergence acceleration , 1982 .

[19]  Patrick W. Fowler,et al.  A model for the geometries of Van der Waals complexes , 1985 .

[20]  J. Pople,et al.  Self‐Consistent Molecular‐Orbital Methods. I. Use of Gaussian Expansions of Slater‐Type Atomic Orbitals , 1969 .

[21]  M. Alderton,et al.  Explicit formulae for the electrostatic energy, forces and torques between a pair of molecules of arbitrary symmetry , 1984 .

[22]  P. Kollman,et al.  An approach to computing electrostatic charges for molecules , 1984 .