Abstract: A complete set of intermolecular potential functions has been developed for use in computer simulations of proteins in their native environment. Parameters are reported for 25 peptide residues as well as the common neutral and charged terminal groups. The potential functions have the simple Coulomb plus Lennard-Jones form and are compatible with the widely used models for water, TIP4P, TIP3P, and SPC. The parameters were obtained and tested primarily in conjunction with Monte Carlo statistical mechanics simulations of 36 pure organic liquids and numerous aqueous solutions of organic ions representative of subunits in the side chains and backbones of proteins. Bond stretch, angle bend, and torsional terms have been adopted from the AMBER united-atom force field. As reported here, further testing has involved studies of conformational energy surfaces and optimizations of the crystal structures for four cyclic hexapeptides and a cyclic pentapeptide. The average root-mean-square deviation from the X-ray structures of the crystals is only 0.17 A for the atomic positions and 3% for the unit cell volumes. A more critical test was then provided by performing energy minimizations for the complete crystal of the protein crambin, including 182 water molecules that were initially placed via a Monte Carlo simulation. The resultant root-mean-square deviation for the non-hydrogen atoms is still ca. 0.2 A and the variation in the errors for charged, polar, and nonpolar residues is small. Improvement is apparent over the AMBER united-atom force field which has previously been demonstrated to be superior to many alternatives. Computer simulations are undoubtedly destined to became an increasingly important means for investigating the structures and dynamics of biomolecular systems.' At the heart of such theoretical calculations are the force fields that describe the interatomic interactions and the mechanics of deformations of the molecules.* There is also little doubt that there will be a continual evolution in force fields with added complexity and improved performance paralleling the availability of computer resources. Our own efforts in this area over the last few years have resulted in the OPLS potential functions for proteins whose development and performance are summarized here. These potential functions have a simple form and they have been parametrized directly to reproduce experimental thermodynamic and structural data on fluids. Consequently, they are computationally efficient and their description of proteins in solution or crystalline environments should be superior to many alterantives that have been developed with limited condensed-phase data. The latter point is pursued here primarily through calculations on the crystal structures for four cyclic hexapeptides, a cyclic pentapeptide, and the protein crambin. Improvements are apparent in comparison to the AMBER united-atom force field3 which has previously been shown to be superior to many alternative^.^ (1) Beveridge, D. L., Jorgensen, W. L., Eds. Ann. N.Y. Acad. Sci. 1986, 482. ( 2 ) For reviews, see: (a) Levitt, M. Annu. Reu. Biophys. Eioeng. 1982, 11, 251. (b) McCammon, J. A. Rep. Prog. Phys. 1984, 47, 1. (3) Weiner, S. J.; Kollman, P. A.; Case, D. A,; Singh, U. C.; Ghio, C.; Alagona, G.; Profeta, S.; Weiner, P. J. Am. Chem. SOC. 1984, 106, 765. Parametrization The peptide residues of proteins contain readily identifiable organic subunits such as amides, hydrocarbons, alcohols, thioethers, etc. In view of this and since data are available on the corresponding pure organic liquids, our approach to developing a force field for proteins was to build it up from parameters demonstrated to yield good descriptions of organic liquids. U1timately, the force field would need to treat both intramolecular terms for bond stretches, angle bends, and torsions, as well as the intermolecular and intramolecular nonbonded interactions. The latter are generally accepted to be the most difficult part of the problem and have been our focus.3 A simple, computationally efficient form was chosen to represent the nonbonded interactions through Coulomb and Lennard-Jones terms interacting between sites centered on nuclei (eq 1). Thus, the intermolecular inter-