New protein parameters are reported for the all-atom empirical energy function in the CHARMM program. The parameter evaluation was based on a self-consistent approach designed to achieve a balance between the internal (bonding) and interaction (nonbonding) terms of the force field and among the solvent−solvent, solvent−solute, and solute−solute interactions. Optimization of the internal parameters used experimental gas-phase geometries, vibrational spectra, and torsional energy surfaces supplemented with ab initio results. The peptide backbone bonding parameters were optimized with respect to data for N-methylacetamide and the alanine dipeptide. The interaction parameters, particularly the atomic charges, were determined by fitting ab initio interaction energies and geometries of complexes between water and model compounds that represented the backbone and the various side chains. In addition, dipole moments, experimental heats and free energies of vaporization, solvation and sublimation, molecular volume...

/pdf/all-atom-empirical-potential-for-molecular-modeling-and-2nds0dxw12.pdf

All-atom empirical potential for molecular modeling and dynamics studies of proteins.

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-

The OPLS [optimized potentials for liquid simulations] potential functions for proteins, energy minimizations for crystals of cyclic peptides and crambin.

IA Basic Concepts.- 1 The Importance of Hydrogen Bonds.- 1.1 Historical Perspective.- 1.2 The Importance of Hydrogen Bonds in Biological Structure and Function.- 1.3 The Role of the Water Molecules.- 1.4 Significance of Small Molecule Crystal Structural Studies.- 1.5 The Structural Approach.- 2 Definitions and Concepts.- 2.1 Definition of the Hydrogen Bond - Strong and Weak Bonds.- 2.2 Hydrogen-Bond Configurations: Two- and Three-Center Hydrogen Bonds Bifurcated and Tandem Bonds.- 2.3 Hydrogen Bonds Are Very Different from Covalent Bonds.- 2.4 The van der Waals Radii Cut-Off Criterion Is Not Useful.- 2.5 The Concept of the Hydrogen-Bond Structure.- 2.6 The Importance of ? and ? Cooperativity.- 2.7 Homo-, Anti- and Heterodromic Patterns.- 2.8 Hydrogen Bond Flip-Flop Disorder: Conformational and Configurational.- 2.9 Proton-Deficient Hydrogen Bonds.- 2.10 The Excluded Region.- 2.11 The Hydrophobic Effect.- 3 Experimental Studies of Hydrogen Bonding.- 3.1 Infrared Spectroscopy and Gas Electron Diffraction.- 3.2 X-Ray and Neutron Crystal Structure Analysis.- 3.3 Treatment of Hydrogen Atoms in Neutron Diffraction Studies.- 3.4 Charge Density and Hydrogen-Bond Energies.- 3.5 Neutron Powder Diffraction.- 3.6 Solid State NMR Spectroscopy.- 4 Theoretical Calculations of Hydrogen-Bond Geometries.- 4.1 Calculating Hydrogen-Bond Geometries.- 4.2 Ab-Initio Molecular Orbital Methods.- 4.3 Application to Hydrogen-Bonded Complexes.- 4.4 Semi-Empirical Molecular Orbital Methods.- 4.5 Empirical Force Field or Molecular Mechanics Methods.- 5 Effect of Hydrogen Bonding on Molecular Structure.- IB Hydrogen-Bond Geometry.- 6 The Importance of Small Molecule Structural Studies.- 6.1 Problems Associated with the Hydrogen-Bond Geometry.- 6.2 The Hydrogen Bond Can Be Described Statistically.- 6.3 The Problems of Measuring Hydrogen-Bond Lengths and Angles in Small Molecule Crystal Structures.- 7 Metrical Aspects of Two-Center Hydrogen Bonds.- 7.1 The Metrical Properties of O-H *** O Hydrogen Bonds.- 7.1.1 Very Strong and Strong OH *** O Hydrogen Bonds Occur with Oxyanions, Acid Salts, Acid Hydrates, and Carboxylic Acids.- 7.1.2 OH *** O Hydrogen Bonds in the Ices and High Hydrates.- 7.1.3 Carbohydrates Provide the Best Data for OH ... O Hydrogen Bonds: Evidence for the Cooperative Effect.- 7.2 N-H *** O Hydrogen Bonds.- 7.3 N-H *** N Hydrogen Bonds.- 7.4 O-H *** N Hydrogen Bonds.- 7.5 Sequences in Lengths of Two-Center Hydrogen Bonds.- 7.6 H/D Isotope Effect.- 8 Metrical Aspects of Three- and Four-Center Hydrogen Bonds.- 8.1 Three-Center Hydrogen Bonds.- 8.2 Four-Center Hydrogen Bonds.- 9 Intramolecular Hydrogen Bonds.- 10 Weak Hydrogen-Bonding Interactions Formed by C-H Groups as Donors and Aromatic Rings as Acceptors.- 11 Halides and Halogen Atoms as Hydrogen-Bond Acceptors.- 12 Hydrogen-Bond Acceptor Geometries.- II Hydrogen Bonding in Small Biological Molecules.- 13 Hydrogen Bonding in Carbohydrates.- 13.1 Sugar Alcohols (Alditols) as Model Cooperative Hydrogen-Bonded Structures.- 13.2 Influence of Hydrogen Bonding on Configuration and Conformation in Cyclic Monosaccharides.- 13.3 Rules to Describe Hydrogen-Bonding Patterns in Monosaccharides.- 13.4 The Water Molecules Link Hydrogen-Bond Chains into Nets in the Hydrated Monosaccharide Crystal Structures.- 13.5 The Disaccharide Crystal Structures Provide an Important Source of Data About Hydrogen-Bonding Patterns in Polysaccharides.- 13.6 Hydrogen Bonding in the Tri- and Tetrasaccharides Is More Complex and Less Well Defined.- 13.7 The Hydrogen Bonding in Polysaccharide Fiber Structures Is Poorly Defined.- 14 Hydrogen Bonding in Amino Acids and Peptides: Predominance of Zwitterions.- 15 Purines and Pyrimidines.- 15.1 Bases Are Planar and Each Contains Several Different Hydrogen-Bonding Donor and Acceptor Groups.- 15.2 Many Tautomeric Forms Are Feasible But Not Observed.- 15.3 ?-Bond Cooperativity Enhances Hydrogen-Bonding Forces.- 15.4 General, Non-Base-Pairing Hydrogen Bonds.- 16 Base Pairing in the Purine and Pyrimidine Crystal Structures.- 16.1 Base-Pair Configurations with Purine and Pyrimidine Homo-Association.- 16.2 Base-Pair Configurations with Purine-Pyrimidine Hetero-Association: the Watson-Crick Base-Pairs.- 16.3 Base Pairs Can Combine to Form Triplets and Quadruplets.- 17 Hydrogen Bonding in the Crystal Structures of the Nucleosides and Nucleotides.- 17.1 Conformational and Hydrogen-Bonding Characteristics of the Nucleosides and Nucleotides.- 17.2 A Selection of Cyclic Hydrogen-Bonding Patterns Formed in Nucleoside and Nucleotide Crystal Structures.- 17.3 General Hydrogen-Bonding Patterns in Nucleoside and Nucleotide Crystal Structures.- III Hydrogen Bonding in Biological Macromolecules.- 18 O-H *** O Hydrogen Bonding in Crystal Structures of Cyclic and Linear Oligoamyloses: Cyclodextrins, Maltotriose, and Maltohexaose.- 18.1 The Cyclodextrins and Their Inclusion Complexes.- 18.2 Crystal Packing Patterns of Cyclodextrins Are Determined by Hydrogen Bonding.- 18.3 Cyclodextrins as Model Compounds to Study Hydrogen-Bonding Networks.- 18.4 Cooperative, Homodromic, and Antidromic Hydrogen-Bonding Patterns in the ?-Cyclodextrin Hydrates.- 18.5 Homodromic and Antidromic O-H *** O Hydrogen-Bonding Systems Analyzed Theoretically.- 18.6 Intramolecular Hydrogen Bonds in the ?-Cyclodextrin Molecule are Variable - the Induced-Fit Hypothesis.- 18.7 Flip-Flop Hydrogen Bonds in ?-Cyclodextrin * 11 H2O.- 18.8 From Flip-Flop Disorder to Ordered Homodromic Arrangements at Low lbmperature: The Importance of the Cooperative Effect.- 18.9 Maltohexaose Polyiodide and Maltotriose - Double and Single Left-Handed Helices With and Without Intramolecular O(2) *** O(3?) Hydrogen Bonds.- 19 Hydrogen Bonding in Proteins.- 19.1 Geometry of Secondary-Structure Elements: Helix, Pleated Sheet, and Turn.- 19.2 Hydrogen-Bond Analysis in Protein Crystal Structures.- 19.3 Hydrogen-Bonding Patterns in the Secondary Structure Elements.- 19.4 Hydrogen-Bonding Patterns Involving Side-Chains.- 19.5 Internal Water Molecules as Integral Part of Protein Structures.- 19.6 Metrical Analysis of Hydrogen Bonds in Proteins.- 19.7 Nonsecondary-Structure Hydrogen-Bond Geometry Between Main-Chains, Side-Chains and Water Molecules.- 19.8 Three-Center (Bifurcated) Bonds in Proteins.- 19.9 Neutron Diffraction Studies on Proteins Give Insight into Local Hydrogen-Bonding Flexibility.- 19.10 Site-Directed Mutagenesis Gives New Insight into Protein Thermal Stability and Strength of Hydrogen Bonds.- 20 The Role of Hydrogen Bonding in the Structure and Function of the Nucleic Acids.- 20.1 Hydrogen Bonding in Nucleic Acids is Essential for Life.- 20.2 The Structure of DNA and RNA Double Helices is Determined by Watson-Crick Base-Pair Geometry.- 20.3 Systematic and Accidental Base-Pair Mismatches: "Wobbling" and Mutations.- 20.4 Noncomplementary Base Pairs Have a Structural Role in tRNA.- 20.5 Homopolynucleotide Complexes Are Stabilized by a Variety of Base-Base Hydrogen Bonds - Three-Center (Bifurcated) Hydrogen Bonds in A-Tracts.- 20.6 Specific Protein-Nucleic Acid Recognition Involves Hydrogen Bonding.- IV Hydrogen Bonding by the Water Molecule.- 21 Hydrogen-Bonding Patterns in Water, Ices, the Hydrate Inclusion Compounds, and the Hydrate Layer Structures.- 21.1 Liquid Water and the Ices.- 21.2 The Hydrate Inclusion Compounds.- 21.3 Hydrate Layer Structures.- 22 Hydrates of Small Biological Molecules: Carbohydrates, Amino Acids, Peptides, Purines, Pyrimidines, Nucleosides and Nucleotides.- 23 Hydration of Proteins.- 23.1 Characterization of "Bound Water" at Protein Surfaces - the First Hydration Shell.- 23.2 Sites of Hydration in Proteins.- 23.3 Metrics of Water Hydrogen Bonding to Proteins.- 23.4 Ordered Water Molecules at Protein Surfaces - Clusters and Pentagons.- 24 Hydration of Nucleic Acids.- 24.1 Two Water Layers Around the DNA Double Helix.- 24.2 Crystallographically Determined Hydration Sites in A-, B-, Z-DNA. A Statistical Analysis.- 24.3 Hydration Motifs in Double Helical Nucleic Acids.- 24.3.1 Sequence-Independent Motifs.- 24.3.2 Sequence-Dependent Motifs.- 24.4 DNA Hydration and Structural Transitions Are Correlated: Some Hypotheses.- 25 The Role of Three-Center Hydrogen Bonds in the Dynamics of Hydration and of Structure Transition.- References.- Refcodes.

Hydrogen Bonding in Biological Structures

A molecular mechanics force field implemented in the Sybyl program is described along with a statistical evaluation of its efficiency on a variety of compounds by analysis of internal coordinates and thermodynamic barriers. The goal of the force field is to provide good quality geometries and relative energies for a large variety of organic molecules by energy minimization. Performance in protein modeling was tested by minimizations starting from crystallographic coordinates for three cyclic hexapeptides in the crystal lattice with rms movements of 0.019 angstroms, 2.06 degrees, and 6.82 degrees for bond lengths, angles, and torsions, respectively, and an rms movement of 0.16 angstroms for heavy atoms. Isolated crambin was also analyzed with rms movements of 0.025 angstroms, 2.97 degrees, and 13.0 degrees for bond lengths, angles, and torsions respectively, and an rms movement of 0.42 angstroms for heavy atoms. Accuracy in calculating thermodynamic barriers was tested for 17 energy differences between conformers, 12 stereoisomers, and 15 torsional barriers. The rms errors were 0.8, 1.7, and 1.13 kcal/mol, respectively, for the three tests. Performance in general purpose applications was assessed by minimizing 76 diverse complex organic crystal structures, with and without randomization by coordinate truncation, with rms movements of 0.025 angstroms, 2.50 degrees, and 9.54 degrees for bond lengths, angles and torsions respectively, and an average rms movement of 0.192 angstroms for heavy atoms.

Validation of the general purpose Tripos 5.2 force field

Carp Muscle Calcium-binding Protein II. STRUCTURE DETERMINATION AND GENERAL DESCRIPTION

Structure of acetylacetone by electron diffraction

A procedure has been developed for the determination of molecular structure by electron diffraction which yields accurate intensity data and obviates the necessity for visual examination of the diffraction photographs. The theory for computing radial distribution curves has been extended to permit accurate curves to be obtained from scattering data covering only a restricted range of angle. From this method, it is possible to obtain not only equilibrium distances but also the probability distributions for the vibrational motion between pairs of atoms in a molecule. The procedure has been applied to CCl4 and CO2 and when comparisons may be made with spectroscopic results, satisfactory agreement is obtained.

Internal Motion and Molecular Structure Studies by Electron Diffraction

We have investigated the heretofore unknown unimolecular decomposition pathway of the explosive molecule diaminodinitroethylene (DADNE). With the use of DFT methods, whose accuracy has been calibrated by means of ab initio calculations (MP2, MP4, G2) on a simpler but related molecule, nitroethylene, we have been able to characterize the entire decomposition reaction pathway. Importantly, we find that the reaction is initiated by a nitro-to-nitrite rearrangement with a calculated energy barrier of magnitude 59.1 kcal/mol obtained by use of B3LYP (59.7 kcal/mol B3P86) which is very close to the experimental activation energy of 58 kcal/mol. We have been able to characterize every step in the decomposition reaction path leading to fragments NO, HONO, CO, NH2, and HNC. These may interact to yield final stable products, CO, N2, H2O with an energy release that on average is adequate to initiate two additional DADNE molecular decompositions, and thus, sustain a chain reaction. The structural parameters we have c...

Proposed Mechanism of 1,1-Diamino-Dinitroethylene Decomposition: A Density Functional Theory Study

We describe a kernel energy method (KEM) for applying quantum crystallography to large molecules, with an emphasis on the calculation of the molecular energy of peptides. The computational difficulty of representing the system increases only modestly with the number of atoms. The calculations are carried out on modern parallel supercomputers. By adopting the approximation that a full biological molecule can be represented by smaller “kernels” of atoms, the calculations are greatly simplified. Moreover, collections of kernels are, from a computational point of view, well suited for parallel computation. The result is a modest increase in computational time as the number of atoms increases, while retaining the ab initio character of the calculations. We describe a test of our method, and establish its accuracy using 15 different peptides of biological interest. © 2005 Wiley Periodicals, Inc. Int J Quantum Chem, 2005

Kernel energy method illustrated with peptides

The diffraction patterns of both silica and germania glasses are consistent with a structure in which nearly all of the atoms belong to tridymite-like regions of up to about 20 angstroms or more that are bonded efficiently together in a manner analogous to that found in twinned crystals.

https://science.sciencemag.org/content/sci/179/4069/177.full.pdf

Jerome Karle

Papers

Structure of acetylacetone by electron diffraction

Internal Motion and Molecular Structure Studies by Electron Diffraction

Proposed Mechanism of 1,1-Diamino-Dinitroethylene Decomposition: A Density Functional Theory Study

Kernel energy method illustrated with peptides

Crystalline ordering in silica and germania glasses.