Stochastic searches for lactone and cycloalkene conformers

The stochastic search method coupled to MM3(92) has been used to locate as many conformers as possible for 6- to 11-membered ring lactones, trans-cycloalkenes, and cis-cycloalkenes. A comparison was carried out between the conformers of lactones and cycloalkenes of the same ring size for each force field. These comparisons were carried out by means of (1) conformational distances, defined as the rms deviation between the dihedral angles of the conformers being compared, and (2), substitution, in which a lactone was transformed into an olefin and vice versa, trying to keep as much as possible the initial geometry, followed by reoptimization. It is found that cycloalkenes and lactones share many common characteristics. The thermodynamic information provided by MM3 was used to study the dependence of conformer population upon (1) temperature and (2) total number of conformers. © 1993 John Wiley & Sons, Inc.     

Conformational energetics of 1,3-dichloropropane as predicted by several calculations methods

Two semiempirical methods (MNDO and AM1), a molecular mechanics technique (MM2) and two ab initio approaches (6–31G* full optimization and 3–21G/6–31G*) were used to calculate the ordering of and energy difference between conformers in 1,3-dichloropropane. The semiempirical methods did not order the conformers properly or predict correct energy differences. Both ab initio methods ordered the conformers and predicted energy differences correctly, with the 6–31G* full optimization performing slightly better. The MM2 results were presented for calculations involving a force field with no hydrogens and a full force field of all atoms. The full force field properly ordered the conformers but did not correctly predict the energy differences. The nonhydrogen field ordered the conformers based on the Cl…Cl nonbonded distance. The data show that conformer stability is not a simple matter of maximizing the Cl…Cl nonbonded distance, but is also related to some other stabilizing interaction(s).     

Molecular mechanics (MM4) calculations on carbonyl compounds part I: aldehydes

Aliphatic aldehydes have been studied with the aid of the MM4 force field. The structures, moments of inertia, vibrational spectra, conformational energies, barriers to internal rotation, and dipole moments have been examined for six compounds (nine conformations). MM4 parameters have been developed to fit the indicated quantities to the wide variety of experimental data. Ab initio (MP2) and density functional theory (B3LYP) calculations have been used to augment and/or replace experimental data, as appropriate. Because more, and to some extent, better, data have become available since MM3 was developed, it was anticipated that the overall accuracy of the information calculated with MM4 would be better than with MM3. The best single measure of the overall accuracy of a force field is the accuracy to which the moments of inertia of a set of compounds (from microwave spectroscopy) can be reproduced. For all of the 20 moments (seven conformations) experimentally known for the aldehyde compounds, the MM4 rms error is 0.30%, while with MM3, the most accurate force field presently available, the rms error over the same set is 1.01%. The calculation of the vibrational spectra was also improved overall. For the four aldehydes that were fully analyzed (over a total of 78 frequencies), the rms errors with MM4 and MM3 are 18 and 38 cm^?1, respectively. These improvements came from several sources, but the major ones were separate parameters involving the carbonyl carbon for formaldehyde, the alkyl aldehydes and the ketones, and new crossterms featured in the MM4 force field that are not present in the MM3 version. © 2001 John Wiley & Sons, Inc. J Comput Chem 22: 1396–1425, 2001     

Accurate assessment of the strain energy in a protein-bound drug using QM/MM X-ray refinement and converged quantum chemistry

An ongoing question regarding the energetics of protein‐ligand binding has been; what is the strain energy that a ligand pays (if any) when binding to its protein target? The traditional method to estimate strain energy uses force fields to calculate the energy difference between the ligand bound conformation and its nearest local minimum/global minimum on the gas‐phase or aqueous phase potential energy surface. This makes the implicit assumption that the underlying force field as well as the reference crystal structure is accurate. Herein, we use ibuprofen as a test case and compare MMFF and ab initio QM methods to identify the local and global minimum conformations. Nine low energy conformations were identified with HF/6‐31G* geometry optimization in vacuo. We also obtained highly accurate relative energies for ibuprofen's conformational energy surface based on M06/aug‐cc‐pVXZ (X = D and T), MP2/aug‐cc‐pVXZ (X = D and T) and the MP2/CBS method (with and without solvent corrections). Moreover, we curate and re‐refine the ibuprofen‐protein complex (PDB 2BXG) using QM/MM X‐ray refinement approaches (HF/6‐31G* was the QM method and the MM model was the AMBER force field ff99sb), which were compared with the low energy conformers to calculate the strain energy. The result indicates that there was an 88% reduction in ibuprofen conformation strain using the QM/MM refined structure versus the original PDB ibuprofen conformations. Furthermore, our results indicate that, due to its inherent limitations in estimating electrostatic interactions, force fields are not suitable to gauge strain energy for charged drug molecules like ibuprofen. The present work offers a carefully validated conformational potential energy surface for a drug molecule as well as a reliable QM/MM re‐refined X‐ray structure that can be used to test current structure‐based drug design approaches. © 2011 Wiley Periodicals, Inc. J Comput Chem, 2011     

Some Comments on the Conformations of Methyl Phenyl Sulfides,Sulfoxides and Sulfones

¹³C-NMR. and He (Iα) photoelectron spectra of alkyl phenyl sulfides, sulfoxides and sulfones have been used to probe how their conformations depend on the size of the alkyl groups R. The results are interpreted to indicate that in the sulfides the S, R-bond is twisted out of the planar conformation with increasing size of the alkyl group, whereas in the sulfoxides and sulfones the preferred conformation with the S, R-bond perpendicular to the phenyl group plane seems to be independent of the size of R. These conclusions are in agreement with previous work on the alkyl phenyl sulfides and provide strong support for previous conjectures concerning the preferred conformations of alkyl phenyl sulfoxides and sulfones.     

Calculation of number and free energy of the conformers of linear alkanes with medium and long chains. Implications for catalysis

A method for calculating the number of rotamers of a linear alkane and of the number of rotamers with a given number of gauche conformations along the chain as a function of the total number of atoms in the chain, using general equations, is presented. A graphical method for generating individual rotamers was applied to the homologs up to decane, which has 1134 rotamers. The steric energies calculated by molecular mechanics (MM2 force field) were used as measures of the heat of formation for the coiled conformations relative to the anti conformer for each molecule, whereas the statistical entropy differences were calculated for classes of coiled rotamers grouped by the number of gauche bonds and steric energy. The free energy values calculated from these components show that already at 400 K hexane exists preferentially in conformations containing gauche bonds. For larger chains the free energy advantage for the coiled chains increases very steeply. The implications for the question of reactions of linear alkanes occurring on the surface or inside the channels of small- and medium-pore zeolites are briefly examined.     

Molecular mechanics study of conformational trends in simple alcohols and ethers. II. Intramolecular hydrogen bonding

In this article we present an extension of our modified MM2(80) force field MM2MOD in which a potential function for hydrogen bonding in alcohols and ethers is included. The results of applying MM2(85), MM2(87), and MM2MOD on ethylene glycol, 2-methoxyethanol and 1,3-propanediol are reported and compared with available experimental data and ab initio results. It is concluded that hydrogen bonding plays an important role in determining the molecular conformations of these systems. © 1992 by John Wiley & Sons, Inc.     

Searching for conformers of nine- to twelve-ring hydrocarbons on the MM2 and MM3 energy surfaces: Stochastic search for interconversion pathways

The stochastic search method was employed to find as many conformers on the MM2 and MM3 energy surfaces as possible for cyclic saturated hydrocarbons with ring sizes from 9 through 12. The number found was 8 MM2 (8 MM3) for 9 rings, 18 MM2 (16 MM3) for 10 rings, 40 MM2 (29 MM3) for 11 rings, and 111 MM2 (90 MM3) for 12 rings. A measure of similarity between pairs of conformers of a compound, called conformational distance, is described. It was used to correlate similar MM2 and MM3 conformers. It was discovered that some conformers on each energy surface are not close to minima on the other surface in rings larger than 9. On refinement with the other optimizer, they changed considerably—going downhill to other previously found minima on the other energy surface or (in a few cases) going to minima which had not been found by direct searches. Conformational distance was also employed as an indication of which pairs of MM2 (or MM3) conformers are likely to interconvert rapidly. A new stochastic procedure of using small kicks was used to search for the most likely interconversion processes among the conformers. There is fairly good agreement between the most facile pathways located by it and unusually short conformational distances. Several additional 12-ring conformers (not found with previous methods) were located through application of this small kick procedure.     

Alcohols,ethers, carbohydrates,and related compounds. IV. Carbohydrates

Ab initio calculations [B3LYP/6-311++G(2d,2p)] have been carried out on 84 conformations of 12 different sugars (hexoses), in both pyranose and furanose forms, with the idea of generating a data base for carbohydrate structural energies that may be used for developing the predictive value of molecular mechanics calculations for carbohydrates. The average value for the apparent gas phase anomeric effect for a series of 31 pairs of pyranose conformations was found to be 1.83 kcal/mol (vs. 2.67 kcal/mol with a smaller basis set used in earlier calculations). In developing MM4 to reproduce these data, it was necessary first to have good energies for simple alcohols and ethers, together with an adequate treatment of hydrogen bonding, and then to include the anomeric effect, and the ethylene glycol type system, as was previously recognized. It was also found that the so-called delta-2 effect, long recognized in carbohydrates, must be explicitly included, in order to obtain acceptable results. When a force field that included all of these items as developed from the small molecules based on the MM4 hydrocarbon force field was applied without any parameter adjustment to the set of hexopyranose and furanose conformations mentioned earlier, the E(beta) - E(alpha) was found to have an average value of 1.88 kcal/mol, versus 1.74 for the quantum calculations. The signed average and RMS deviations of the MM4 from the QM results were +0.15 and 0.87 kcal/mol.     

1, 3-二氟丙烷构象和电子结构的从头算和 密度泛函理论计算

用从头算和密度泛函数理论研究了1,3-二氟丙烷的构象和电子结构。在计算的各种理论水平下,GG构象是最稳定构象,AG构象次之。利用分子力场的非键作用定量化旁式效应,由MM2力场的非键参数计算的结果较为合理。对于GG和AG构象,在HF/6-31G^*和HF/6-31+G^*^*水平预测的构象分布与实验值接近。     

Molecular mechanics conformational analysis of cyclononane using the RIPS method and comparison with quantum-mechanical calculations

The recently reported Random Incremental Pulse Search (RIPS) technique has been used to probe the conformational energy surface of cyclononane. The stochastic method permits searching of the potential energy surface for all minimum-energy conformations. The search located all previously reported structures together with three additional conformations that were not found by earlier, primitive searching techniques. Two of these structures are high-nergy skew forms, and the third is a low-energy conformer that should contribute significantly to the overall equilibrium set of cyclononane conformations. The global minimum has been found to be the D₃ symmetrical twist chair-boat (TBC) form in accordance with previous studies. The newly discovered low-energy structure, which lies only 2.2 kcal/mol above the global minimum, has been designated twist chair-twist chair (TCTC). The two higher energy conformers are skewed chair-chair (SCC) and skewed boat-boat (SBB) forms that are 5.7 kcal/mol and 10.4 kcal/mol above the global minimum, respectively. The seven reported conformations were reanalyzed quantum mechanically (AM 1), and a comparison between MM 2 and AM 1 results is presented.     

Molecular mechanics calculations on the differentiation of diastereomeric complexes ofcis-decalin withβ-cyclodextrin

Molecular mechanics calculations with the latest available version of Allinger's MM2 force field (MM2(91)) on the diastereomeric complexes of both enantiomeric conformations ofcis-decalin with -cyclodextrin show a small preference (1.67 kJ mol^–1) for one of them, in agreement with the available¹³C-NMR results. Calculations were found to be sensitive to the procedure used.     

Conformational Analysis of C-Disaccharides using Molecular Mechanics Calculations

ABSTRACT

Relaxed-residue energy maps based on the MM3 force field were computed for the methyl glycosides of eight C-linked D-glucosyl disaccharides: the two-bond axial-equatorial linked disaccharides β-kojibioside [(1→2)α–], β-nigeroside [(1→3)α–] and β-maltose [(1→4)α–], the two-bond equatorial-equatorial linked disaccharides β-sophoroside [(1→2)β–], β–laminarabioside [(1→3)β-], β–cellobioside [(1→4)β–] and the three-bond-linked (1→6) disacharides C-isomaltoside and C-gentiobioside. Optimized structures were calculated on a 20° grid spacing of the torsional angles about the C-glycosidic bonds and the final isoenergy surfaces were based on 11664 conformations, for the two-bond-linked disaccharides and 69984 conformations for the three-bond-linked disaccharides. Boltzmann-weighted ³J coupling constants were calculated and compared to the experimental values. They are satisfactory except for maltose where hydrogen bonds cause an over-estimation of the energy differences between the conformers. The energy maps are similar to maps of the corresponding O-disaccharides, but there are differences in the locations and the relative energies of the minima. The preferred conformations of the C-glycosidic bonds are as if they were conforming to the exo-anomeric effect but are closer to staggered conformations than shown by the MM3 results for the O-linkages.     

Molecular mechanics calculations of conjugated amides and an ab initio investigation of acrylamide and its β-amino derivative: Conformational analysis and rotational barriers

Conformations and rotational barriers in a series of conjugated primary and tertiary amides have been analyzed by a modified MM2(91) force field, which treats the amide nitrogen as part of the conjugated system by redefining the atom type for the nitrogen. Ab initio molecular orbital calculations at the MP2/6-31G* level have been performed on the stable conformers and transition structures of acrylamide and β-trans-aminoacrylamide. The results have been used, with published experimental and computational data, to generate parameters for the MM2 force field. The force field has been applied to various conjugated amides, such as reduced nicotinamide adenine dinucleotide (NADH) and NAD⁺ analogues, nicotinamide, urea, vinylogous urea derivatives, and nucleic acid bases. The fundamental difference between primary and tertiary conjugated amides with respect to both conformation and barrier is highlighted. © 1996 by John Wiley & Sons, Inc.     

Conformational preferences in alkylbenzenes and aryl-alkylamines: A comparative study using CAMSEQ,MM2 and molecular dynamics methods

Conformational preferences in a series of alkylbenzenes and protonated and neutral aryl-alkylamines of biological interest have been examined. First, a general picture was obtained for the 25 compounds in the series by means of the CAMSEQ empirical potential software system. This provided solution as well as vacuum data. A low-energy “folded” conformation (alkyl chain coiled toward aromatic ring) was observed in every case. The presence of an amino nitrogen atom in the alkyl chain did not significantly influence conformational preference. Secondly, a group of 14 compounds, representative of the subgroupings within the main series, was selected and the previously established (CAMSEQ) folded and extended minimum positions were further examined by means of a modified MM2 program. Energy differences between conformers were also calculated, and the presence of a stable folded form was confirmed. A feature common to many of the folded conformations is the positioning over the ring of a hydrogen atom from the terminal group of the chain. The MM2 folded/extended energy differences were in some cases smaller than those determined by CAMSEQ, with the “extended” form in many cases being about 1 or 2 kcal/mol more stable. Thirdly, three representative compounds from the series were examined by means of a molecular dynamics program which permitted the sampling of conformational space throughout the transition from extended to folded forms. This method gave energy differences between folded and extended conformations which agreed with the corresponding MM2 differences.     

Conformational properties of the macrocyclic trichothecene mycotoxin verrucarin A in solution

Phase-sensitive nuclear Overhauser enhancement spectroscopy (NOESY) experiments, (3)J couplings and computational molecular modeling (MM2* and MMFF force fields) were employed to examine the conformational properties of verrucarin A in chloroform solutions. The MMFF force field calculations resulted in a family of 12 low-energy structures along with their populations, the latter being determined by the NMR analysis of molecular flexibility in solution(NAMFIS) deconvolution analysis. The concluded model was capable of reproducing successfully the experimental NOESY cross-peak volumes and the proton-coupling constants. Among the 12 conformers, the one which was similar to the structure of verrucarin A in the solid state was the predominant accounting for 75% of the total relative population, although other low-energy conformations contributed to a lesser degree in order to explain the experimental data.     

Molecular mechanics conformational analysis of tylosin

The conformations of the 16-membered macrolide antibiotic tylosin were studied with molecular mechanics (AMBER* force field) including modelling of the effect of the solvent on the conformational preferences (GB/SA). A Monte Carlo conformational search procedure was used for finding the most probable low-energy conformations. The present study provides complementary data to recently reported analysis of the conformations of tylosin based on NMR techniques. A search for the low-energy conformations of protynolide, a 16-membered lactone containing the same aglycone as tylosin, was also carried out, and the results were compared with the observed conformation in the crystal as well as with the most probable conformations of the macrocyclic ring of tylosin. The dependence of the results on force field was also studied by utilizing the MM3 force field. Some particular conformations were computed with the semiempirical molecular orbital methods AM1 and PM3.     

Peptide free‐energy profile is strongly dependent on the force field: Comparison of C96 and AMBER95

The C96 and AMBER95 force fields were compared with small model peptides Ac‐(Ala)_n‐NMe (Ac = CH₃CO, NMe = NHCH₃, n=2 and 3) in vacuo and in TIP3P water by computing the free‐energy profiles using multicanonical molecular dynamics method. The C96 force field is a modified version of the AMBER95 force field, which was adjusted to reproduce the energy difference between extended β‐ and constrained α‐helical energies for the alanine tetrapeptide, obtained by the high level ab initio MO method. The slight modification resulted in a large difference in the free energy profiles. The C96 force field prefers relatively extended conformers, whereas the AMBER95 force field favors turn conformations. © 2000 John Wiley & Sons, Inc. J Comput Chem 21: 748–762, 2000     

COMPUTER SIMULATIONS ON UNPERTURBED CONFIGURATIONAL DIMENSIONS OF POLY (METHYL ACRYLATE)

A full-relaxation optimization of molecule and the popular MM2 force field are em-ployed to obtain the geometry parameters and the conformational energy surface of a mesoor a racemic dyad of poly(methyl acrylate) (PMA) with a specified carbonyl-bond orien-tation in side-groups. It is found that the conformational energy maps calculated hereconsiderably differ from those calculated with the rigid molecular model as reported in theearlier studies. The g~- state cannot be omitted in the obtained contour maps. Two impor-tant conformers tg~- and g~(-t) with energy minima were newly detected for a racemic dyad.The analysis on the conformations with energy minima confirmed that the ester groupsare not always perpendicular to the plane defined by the two adjacent skeletal bonds andmay change their relative orientations to meet the requirement of lower energies duringthe conformational state transition. Instead of the early way of adjusting the interactionenergy parameters to fit the experimental data, we attempt to predict unperturbed chaindimensions via the reliable force field and the configurational statistical mechanics. Theproposed scheme with three rotational states identified from the contour maps allowed usto satisfactorily reproduce the experimental dimensions of random PMA chains.     

MM3 Modeling of Ribose and 2-Deoxyribose Ring Puckering

ABSTRACT

Conformational energy maps for furanosyl and pyranosyl rings of α- and β-ribose and 2-deoxyribose were generated with the molecular mechanics program MM3. For the furanosyl tautomers, low-energy Northern and Southern minima were found. For the pyranosyl rings, lowest-energy minima corresponded to chair forms. Hydrogen-hydrogen coupling constants for the minimal-energy conformers were calculated based on Karplus equations. Computational and experimental results indicate that several tautomeric forms of ribose and 2-deoxyribose exist in multiple conformations in solution. Ring conformations of related crystal structures have energies within ～2 kcal/mol of the calculated global minima.