A new method for obtaining optimized parameters for semiempirical methods has been developed and applied to the modified neglect of diatomic overlap (MNDO) method. The method uses derivatives of calculated values for properties with respect to adjustable parameters to obtain the optimized values of parameters. The large increase in speed is a result of using a simple series expression for calculated values of properties rather than employing full semiempirical calculations. With this optimization procedure, the rate-determining step for parameterizing elements changes from the mechanics of parameterization to the assembling of experimental reference data.

Optimization of parameters for semiempirical methods I. Method

Scaling factors for obtaining fundamental vibrational frequencies, low-frequency vibrations, zero-point vibrational energies (ZPVE), and thermal contributions to enthalpy and entropy from harmonic frequencies determined at 19 levels of theory have been derived through a least-squares approach. Semiempirical methods (AM1 and PM3), conventional uncorrelated and correlated ab initio molecular orbital procedures [Hartree−Fock (HF), Moller−Plesset (MP2), and quadratic configuration interaction including single and double substitutions (QCISD)], and several variants of density functional theory (DFT:  B-LYP, B-P86, B3-LYP, B3-P86, and B3-PW91) have been examined in conjunction with the 3-21G, 6-31G(d), 6-31+G(d), 6-31G(d,p), 6-311G(d,p), and 6-311G(df,p) basis sets. The scaling factors for the theoretical harmonic vibrational frequencies were determined by a comparison with the corresponding experimental fundamentals utilizing a total of 1066 individual vibrations. Scaling factors suitable for low-frequency vib...

Harmonic Vibrational Frequencies: An Evaluation of Hartree−Fock, Møller−Plesset, Quadratic Configuration Interaction, Density Functional Theory, and Semiempirical Scale Factors

A significant modification to the additive all-atom CHARMM lipid force field (FF) is developed and applied to phospholipid bilayers with both choline and ethanolamine containing head groups and with both saturated and unsaturated aliphatic chains. Motivated by the current CHARMM lipid FF (C27 and C27r) systematically yielding values of the surface area per lipid that are smaller than experimental estimates and gel-like structures of bilayers well above the gel transition temperature, selected torsional, Lennard-Jones and partial atomic charge parameters were modified by targeting both quantum mechanical (QM) and experimental data. QM calculations ranging from high-level ab initio calculations on small molecules to semiempirical QM studies on a 1,2-dipalmitoyl-sn-phosphatidylcholine (DPPC) bilayer in combination with experimental thermodynamic data were used as target data for parameter optimization. These changes were tested with simulations of pure bilayers at high hydration of the following six lipids: ...

Update of the CHARMM All-Atom Additive Force Field for Lipids: Validation on Six Lipid Types

MNDO/AM1‐type parameters for twelve elements have been optimized using a newly developed method for optimizing parameters for semiempirical methods. With the new method, MNDO‐PM3, the average difference between the predicted heats of formation and experimental values for 657 compounds is 7.8 kcal/mol, and for 106 hypervalent compounds, 13.6 kcal/mol. For MNDO the equivalent differences are 13.9 and 75.8 kcal/mol, while those for AM1, in which MNDO parameters are used for aluminum, phosphorus, and sulfur, are 12.7 and 83.1 kcal/mol, respectively. Average errors for ionization potentials, bond angles, and dipole moments are intermediate between those for MNDO and AM1, while errors in bond lengths are slightly reduced.

Optimization of parameters for semiempirical methods II. Applications

Several modifications that have been made to the NDDO core-core interaction term and to the method of parameter optimization are described. These changes have resulted in a more complete parameter optimization, called PM6, which has, in turn, allowed 70 elements to be parameterized. The average unsigned error (AUE) between calculated and reference heats of formation for 4,492 species was 8.0 kcal mol−1. For the subset of 1,373 compounds involving only the elements H, C, N, O, F, P, S, Cl, and Br, the PM6 AUE was 4.4 kcal mol−1. The equivalent AUE for other methods were: RM1: 5.0, B3LYP 6–31G*: 5.2, PM5: 5.7, PM3: 6.3, HF 6–31G*: 7.4, and AM1: 10.0 kcal mol−1. Several long-standing faults in AM1 and PM3 have been corrected and significant improvements have been made in the prediction of geometries.

/pdf/optimization-of-parameters-for-semiempirical-methods-v-597lxcumix.pdf

Optimization of parameters for semiempirical methods V: modification of NDDO approximations and application to 70 elements.

Studies have found that mutant, misfolded superoxide dismutase [Cu–Zn] (SOD1) can convert wild type SOD1 (wtSOD1) in a prion-like fashion, and that misfolded wtSOD1 can be propagated by release and uptake of protein aggregates. In developing a prion-like mechanism for this propagation of SOD1 misfolding we have previously shown how enervation of the SOD1 electrostatic loop (ESL), caused by the formation of transient non-obligate SOD1 oligomers, can lead to an experimentally observed gain of interaction (GOI) that results in the formation of SOD1 amyloid-like filaments. It has also been shown that freedom of ESL motion is essential to catalytic function. This work investigates the possibility that restricting ESL mobility might not only compromise superoxide catalytic activity but also serve to promote the peroxidase activity of SOD1, thus implicating the formation of SOD1 oligomers in both protein misfolding and in protein oxidation.

A model for gain of function in superoxide dismutase.

A study of the aluminum hydride reduction of unsaturated cyclic epoxides

The matrix metalloproteinase family has been a pharmaceutical target for most of the last three decades, but success has been hampered by unwanted side effects caused by lack of selectivity, poor oral bioavailability and decreased potency in vivo. The surface-expressed metalloproteinases ADAM10 and ADAM17, the latter also referred to as TACE, play important roles in various physiological processes, especially involving tissue repair and development. Because of its role in the release of the cytokine TNF-α TACE has been a key target for pharmaceutical intervention in the treatment of rheumatoid arthritis. An extensive body of structural activity data has been developed for a series of small molecule inhibitors of TACE based on a sulfonamide scaffold containing key acetylenic substituents. We have undertaken an extensive molecular modeling study of select members of this ligand group to better understand the structural nuances involved in the development of ever more potent TACE inhibitors, and identify those elements of structure-based design that would enhance the selectivity of such inhibitors for TACE over ADAM10. Results include the identification of a flexible loop, comparable to that found in other MMPs that plays a subtle, yet significant, role in determining inhibitor potency.

Acetylenic inhibitors of ADAM10 and ADAM17: in silico analysis of potency and selectivity.

While synthesis projects are not uncommon in the undergraduate organic curriculum, the problems posed are usually not representative of what chemists actually do. The sudden dislocation to online l...

Eamonn F. Healy

Papers

A model for gain of function in superoxide dismutase.

A study of the aluminum hydride reduction of unsaturated cyclic epoxides

Acetylenic inhibitors of ADAM10 and ADAM17: in silico analysis of potency and selectivity.

Tips and Tools for Teaching Organic Synthesis Online

A mechanism of action for small heat shock proteins