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.     

A molecular mechanics and semiempirical conformational analysis of the herbicide diuron inhibitor of photosystem II

The potential energy surface (PES) for the herbicide diuron (DCMU), a photosystem II inhibitor, has been extensively investigated using the quantum-mechanical semiempirical molecular orbital methods AM1 and PM3 and molecular mechanics method. A detailed conformational search has been carried out which revealed the occurrence of four genuine minimum energy structures. The relative stability of the conformers and rotational barriers to conformational interconversion were evaluated using distinct theoretical approaches. The results showed that thetrans form of the diuron molecule is more stable than thecis form in all methods, and so it may possibly be the biologically active isomer.     

Exploring the similarities between potential smoothing and simulated annealing

Simulated annealing and potential function smoothing are two widely used approaches for global energy optimization of molecular systems. Potential smoothing as implemented in the diffusion equation method has been applied to study partitioning of the potential energy surface (PES) for N‐Acetyl‐Ala‐Ala‐N‐Methylamide (CDAP) and the clustering of conformations on deformed surfaces. A deformable version of the united‐atom OPLS force field is described, and used to locate all local minima and conformational transition states on the CDAP surface. It is shown that the smoothing process clusters conformations in a manner consistent with the inherent structure of the undeformed PES. Smoothing deforms the original surface in three ways: structural shifting of individual minima, merging of adjacent minima, and energy crossings between unrelated minima. A master equation approach and explicit molecular dynamics trajectories are used to uncover similar features in the equilibrium probability distribution of CDAP minima as a function of temperature. Qualitative and quantitative correlations between the simulated annealing and potential smoothing approaches to enhanced conformational sampling are established. © 2000 John Wiley & Sons, Inc. J Comput Chem 21: 531–552, 2000     

Proton mobility in protonated peptides: a joint molecular orbital and RRKM study

The mobile proton model was critically evaluated by using purely theoretical models which include quantum mechanical calculations to determine stationary points on the potential energy surface (PES) of a model compound, and Rice-Ramsperger-Kassel-Marcus (RRKM) calculations to determine the rate constants of various processes (conformational changes, proton transfer reactions) which occur during mass analysis of protonated peptides. Extensive mapping of the PES of protonated N-formylglycinamide resulted in various minima which were stabilized by one or more of the following types of interaction: internal hydrogen bond, charge transfer interaction, charge delocalization, and ring formation. The relative energies of most of the investigated minima are less then 20 kcal mol(-1) compared with the most stable species. More importantly, the relative energies of the transition structures connecting these minima are fairly low, allowing facile transitions among the energetically low-lying species. It is demonstrated that a path can be found leading from the energetically most stable species, protonated on an amide oxygen, to the structure from which the energetically most favorable fragmentation occurs. It is also shown that the added proton can sample all protonation sites prior to fragmentation. The RRKM calculations applied the results of ab initio computations (structures, energetics, vibrational frequencies) to the reactions (internal rotations, proton transfers) occurring in protonated N-formylglycinamide, and clearly lend additional evidence to the mobile proton model. Based on the results of the PES search on protonated N-formylglycinamide, we also comment on the mechanism proposed by Arnot et al. (Arnot D, Kottmeier D, Yates N, Shabanowitz J, Hunt D F. 42(nd) ASMS Conference on Mass Spectrometry, 1994; 470) and Reid et al. (Reid G E, Simpson R J, O'Hair R A J. J. Am. Soc. Mass Spectrom. 1998; 9:945) for the formation of b(2)(+) ions. According to the high level ab initio results, the mechanism relying on amide oxygen protonated species seems to be less feasible than the one which involves N-protonated species.     

1,5-Diamino-1H-1,2,3,4-tetrazolium picrate: X-ray molecular and crystal structures and ab initio MO calculations

The crystal and molecular structures of 1,5-diamino-1H-1,2,3,4-tetrazolium picrate (DATP) were determined by X-ray diffraction analysis. The tetrazolium cation in DATP has a structure with protonated N4 atom of the ring. Two amino groups in the cation are found to be rather different. The 5-amino group lies in the plane of the tetrazole ring and valence angles around the N atom are close to 120°, which indicates sp² hybridization of atomic orbital of the nitrogen atom. In contrast, valence angles around the N atom of the 1-amino group are close to tetrahedral angle, which suggests sp³ hybridization. The exocyclic C-N bond in the cation is substantially shorter than that in 1,5-diaminotetrazole. The obtained results indicate a conjugation between the π-system of the tetrazole ring and the 5-amino group. The results of ab initio calculations of electronic structure and relative stability for various tautomeric forms of protonated 1,5-diaminotetrazole using MP2/6-31G* and B3LYP/6-31G* levels of theory are in a good agreement with X-ray data and show that there are differences in σ-electron overlap populations for the C-N bonds in the cation in DATP, while π-electrons are delocalized.     

Conformational properties of penicillins: quantum chemical calculations and molecular dynamics simulations of benzylpenicillin

Herein, we present theoretical results on the conformational properties of benzylpenicillin, which are characterized by means of quantum chemical calculations (MP2/6-31G* and B3LYP/6-31G*) and classical molecular dynamics simulations (5 ns) both in the gas phase and in aqueous solution. In the gas phase, the benzylpenicillin conformer in which the thiazolidine ring has the carboxylate group oriented axially is the most favored one. Both intramolecular CH. O and dispersion interactions contribute to stabilize the axial conformer with respect to the equatorial one. In aqueous solution, a molecular dynamics simulation predicts a relative population of the axial:equatorial conformers of 0.70:0.30 in consonance with NMR experimental data. Overall, the quantum chemical calculations as well as the simulations give insight into substituent effects, the conformational dynamics of benzylpenicillin, the frequency of ring-puckering motions, and the correlation of side chain and ring-puckering motions.     

The rotational barrier in ethane: a molecular orbital study

The energy change on each Occupied Molecular Orbital as a function of rotation about the C-C bond in ethane was studied using the B3LYP, mPWB95 functional and MP2 methods with different basis sets. Also, the effect of the ZPE on rotational barrier was analyzed. We have found that σ and π energies contribution stabilize a staggered conformation. The σ(s) molecular orbital stabilizes the staggered conformation while the stabilizes the eclipsed conformation and destabilize the staggered conformation. The π(z) and molecular orbitals stabilize both the eclipsed and staggered conformations, which are destabilized by the π(v) and molecular orbitals. The results show that the method of calculation has the effect of changing the behavior of the energy change in each Occupied Molecular Orbital energy as a function of the angle of rotation about the C-C bond in ethane. Finally, we found that if the molecular orbital energy contribution is deleted from the rotational energy, an inversion in conformational preference occurs.     

Semi-empirical (PM3 and AM1) and ab initio molecular orbital calculations of 1,2,4-oxadiazoles, 4,5-dihydro-1,2,4-oxadiazoles and 4,4-di-n-butyl-2-phenylbenzo-1,3-oxazine

Molecular orbital calculations using semi-empirical (PM3 and AM1) and ab initio (HF/6-31G) types have been carried out on several 3,5-disubstituted 1,2,4-oxadiazoles 1a–d, 5-n-butyl-3,5-diaryl-4,5-dihydro-1,2,4-oxadiazoles 2a–d, and 4,4-di-n-butyl-2-phenylbenzo-1,3-oxazine 3. A comparison of the results by the two computational procedures has been made. Transformation of the oxadiazole ring to 4,5-dihydro-1,2,4-oxadiazole having both aryl and n-butyl groups at C-5 exhibited interesting conformational features. Also, examination of 1,3-oxazine 3 gave an idea about the structure of this compound. The rotational barrier of each phenyl group in 1a and 1d has been calculated using the ab initio method HF/6-31G(d).     

Molecular dynamics, density functional and second-order Møller–Plesset theory study of the structure and conformation of acetylcholine in vacuo and in solution

 Molecular mechanics minimizations based on the CVFF force field and molecular dynamics simulation for a time of 2.5 ns were performed to examine the conformational behaviour and the molecular motion of acetylcholine in vacuo and in aqueous solution. Five low-lying conformations, namely the TT, TG, GG, G*G and GT, were obtained from molecular mechanics computations with the GT structure as the absolute minimum. Molecular dynamics trajectories in vacuo and in water show that only four (GT, GG, G*G and TG) and three (TG, TT and GT) conformations are present in the simulation time, respectively. Density functional B3LYP and second-order M?ller–Plesset (MP2) methods were then used to study all the five lowest-lying conformers of acetylcholine neurotransmitter in vacuo and in water by the polarizable continuum model approach. The B3LYP and MP2 computations show that in the gas phase all minima lie in a narrow range of energy with the G*G conformer as the most stable one. The relative minima GG, GT, TG and TT are located at 1.1 (3.3), 1.8 (4.2), 2.1 (4.5) and 4.3 (7.3) kcal/mol above the absolute one at the B3LYP (MP2) level. The preferred conformation in water is the TG. Solvation reduces the relative energy differences between the five minima in both computations. Received: 4 April 2001 / Accepted: 5 July 2001 / Published online: 30 October 2001     

四唑及其衍生物的理论研究（6）

对7个肼基四唑衍生物进行了HF／6－31G水平的几何构型全优化及MP2／6－31G／／HF／6－31G水平的总能量计算．结果表明，四唑环基本上取平面构型，肼基不与环共面．键长与重叠布居之间不存在平行关系．环上N1和N4原子荷负电多、利于质子化．1H．肼基四唑的HOMO－LUMO能级差（△E）较2H式的大，中性分子的△E教相应的负离子的大．还计算研究了7个标题物的红外光谱并由此求得了它们的热力学性质．     

Ab initio calculations and molecular mechanics (MM3) force field development for ammonium and protonated aliphatic amines

The geometries and vibrational frequencies of 11 training molecules containing the ammonium ion moiety were calculated at the MP2/6-31+G* level of theory. Various torsional energy profiles were also calculated using this basis set. From those ab initio calculations, a molecular mechanics (MM3) force field was developed using our Parameter Analysis and Refinement Toolkit System (PARTS). Using this set of parameters, the MM3 force field was found to well reproduce the molecular geometries and vibrational spectra for the all training molecules. CPU time was reduced from days to seconds. The availability of this new force field dramatically increases the feasibility of the computer-assisted drug design involving ammonium and protonated amino groups. © 1997 by John Wiley & Sons, Inc. J Comput Chem 18 : 1371–1391, 1997     

Thermodynamic conformational analysis and structural stability of the nicotinic analgesic ABT-594

This work presents a theoretical study of the nicotinic analgesic ABT-594. We describe its neutral (precursor) and protonated (active) forms in vacuum and aqueous solution at the MP2/cc-pVDZ level. A conformational analysis is performed on the two torsional angles describing the orientation of the azetidinyl group and the azetidinylmethoxy moiety. To account for entropic effects, a thermostatistical study of conformational populations at physiological temperature is carried out. In the neutral form, conformer I is found as the most populated in vacuum and solution. Here, the nitrogen of the azetidinyl group is far from the electron pairs of the oxygen and the pyridinic nitrogen. In the protonated form, conformer VIII is the most stable in vacuum and solution. Now, the additional proton on the azetidinyl group is oriented toward the electron lone pairs of oxygen. The structural stability of conformers I and VIII is considered through the atoms in molecules theory. The conformer I, in the neutral forms, is stabilized by an intramolecular hydrogen bond. The preference of conformer VIII in the protonated forms is explained by the higher strength of its intramolecular hydrogen bond over the cation-pi interaction found in conformer I. The effect of the interaction energy with the receptor on the conformational preferences of protonated ABT-594 is simulated. The result is that the population of conformers associated to the rotation of the azetidinyl group increases. So, the molecule can easily adopt the optimal internitrogen separation for interaction with the receptor.     

Molecular orbital studies of polyethylene deformation

Young's modulus E for polyethylene in the chain direction is calculated with molecular orbital theory applied to n-alkanes C₃H₈ through n-C₁₃H₂₈ and analyzed with the cluster-difference method. Semiempirical CNDO, MNDO, and AM1 models and ab initio HF/STO-3G, HF/6-31G, HF/6-31G*, and MP2/6-31G* models are used. Cluster-difference results, when extrapolated to infinite chain length, give E in good agreement with moduli evaluated with molecular cluster or crystal orbital methods, provided minimal basis sets are employed. E decreases from 495 GPa (CNDO) to 336 GPa (MP2/6-31G*) as the level of theory is improved, consistent with established behaviors of the various models. Our calculations do not reproduce earlier molecular cluster or crystal orbital results, which gave E < 330 GPa. The most rigorous MP2/6-31G* model is known to overestimate force constants by ∼ 11%; the scaled modulus E = 299 GPa is in good accord with E = 306 GPa from recent calculations based on experimental vibration frequencies. © 1996 John Wiley & Sons, Inc.     

可乐定分子构象异构和互变异构的理论研究

采用HF/6-31G*方法, 对氨基型和亚胺型可乐定中性分子以及质子化的可乐定分子进行了构象分析. 之后采用B3LYP/6-31G**方法对处于势能面上的能量极小构象异构体进行全自由度几何优化和频率分析. Onsager反应场溶剂模型用于水相的计算. 结果表明, 在气相和水相中可乐定中性分子主要以亚胺型存在, 这同已有实验结果一致. 进一步, 寻找构象异构化过渡态和氨基型亚胺型互变异构化过渡态, 探讨质子化和水溶剂化效应对异构化过程的几何结构和能量的影响. 为了考察氯苯的共轭效应对可乐定互变异构体稳定性的影响, 在B3LYP/6-31G**水平上, 研究了2-氨基-2-咪唑啉的互变异构化反应机理.     

Ab initio and DFT studies of the molecular structures and vibrational spectra of succinonitrile

The Molecular structure, conformational stability and vibrational frequencies of succinonitrile NCCH2CH2CN have been investigated with ab initio and density functional theory (DFT) methods implementing the standard 6-311++G* basis set. The potential energy surfaces (PES) have been explored at DFT-B3LYP, HF and MP2 levels of theory. In agreements with previous experimental results, the molecule was predicted to exist in equilibrium mixture of trans and gauche conforms with the trans form being slightly lower in energy. The vibrational frequencies and the corresponding vibrational assignments of succinonitrile in both C2h and C2 symmetry were examined theoretically and the calculated Infrared and Raman spectra of the molecule were plotted. Observed frequencies for normal modes were compared with those calculated from normal mode coordinate analysis carried out on the basis of ab initio and DFT force fields using the standard 6-311++G* basis set of the theoretical optimized geometry. Theoretical IR intensities and Raman activities are reported.     

Development of polyphosphate parameters for use with the AMBER force field

Accurate force fields are essential for reproducing the conformational and dynamic behavior of condensed-phase systems. The popular AMBER force field has parameters for monophosphates, but they do not extend well to polyphorylated molecules such as ADP and ATP. This work presents parameters for the partial charges, atom types, bond angles, and torsions in simple polyphosphorylated compounds. The parameters are based on molecular orbital calculations of methyldiphosphate and methyltriphosphate at the RHF/6-31+G* level. The new parameters were fit to the entire potential energy surface (not just minima) with an RMSD of 0.62 kcal/mol. This is exceptional agreement and a significant improvement over the current parameters that produce a potential surface with an RMSD of 7.8 kcal/mol to that of the ab initio calculations. Testing has shown that the parameters are transferable and capable of reproducing the gas-phase conformations of inorganic diphosphate and triphosphate. Also, the parameters are an improvement over existing parameters in the condensed phase as shown by minimizations of ATP bound in several proteins. These parameters are intended for use with the existing AMBER 94/99 force field, and they will permit users to apply AMBER to a wider variety of important enzymatic systems.     

Transiting the molecular potential energy surface along low energy pathways: The TRREAT algorithm

The Transition Rapidly exploring Random Eigenvector Assisted Tree (TRREAT) algorithm is introduced to perform searches along low curvature pathways on a potential energy surface (PES). The method combines local curvature information about the PES with an iterative Rapidly exploring Random Tree algorithm (LaValle, Computer Science Department, Iowa State University, 1998, TR98–11) that quickly searches high‐dimensional spaces for feasible pathways between local minima. Herein, the method is applied to identifying conformational changes of molecular systems using Cartesian coordinates while avoiding a priori definition of collective variables. We analyze the pathway identification problem for alanine dipeptide, cyclohexane and glycine using nonreactive and reactive forcefields. We show how TRREAT‐identified pathways can be used as valuable input guesses for double‐ended methods such as the Nudged Elastic Band when ascertaining transition state energies. This method can be utilized to improve/extend the reaction databases that lie at the core of automatic chemical reaction mechanism generator software currently developed to build kinetic models of chemical reactions. © 2013 Wiley Periodicals, Inc.     

Ab initio and DFT calculations on the initial step in thiamin catalysis

The diphosphate ester (ThDP) of thiamin (vitamin B₁) is an important cofactor of enzymes within the carbohydrate metabolism. The initial reaction step shared by all ThDP-dependent enzymes is the deprotonation of the C2–H of the thiazolium ring. The replacement of the 4′-amino group by a hydroxyl one in the pyrimidine ring leads to the oxy-ThDP analogue which is known as an antagonist in thiamin catalysis.

Ab initio and DFT calculations on the MP2/6-31G* and B3LYP/6-31G* level were performed to study the proton relay function in thiamin and oxythiamin systems. Both MP2 and B3LYP calculations show significant differences of the reaction coordinate of the ylide formation in the systems. Tautomers, protonated and deprotonated species of both systems show different trends regarding their stability. The influence of correlation effects on the results is discussed by comparison with the HF-SCF/6-31G* data. Frequency calculations on the B3LYP/6-31G* level were performed to characterize the minima and transition state structures, respectively.     

A mass spectrometric and molecular orbital study of H2O loss from protonated tryptophan and oxidized tryptophan derivatives

Protonated N-acetyltryptophan, oxindolylalanine (a mono-oxidized derivative of tryptophan), and N-acetyloxindolylalanine, as well as several di- and tripeptide derivatives containing oxindolylalanine, undergo a range of fragmentation reactions in the gas phase, including the loss of water. In order to elucidate the sites of water loss within these ions, and to determine the mechanisms associated with these processes, we have conducted a series of experiments employing multistage tandem mass spectrometry (MS/MS and MS(3)) in a quadrupole ion trap mass spectrometer, regiospecific structural labeling, and independent solution-phase syntheses of proposed product ion structures, coupled with the use of molecular orbital calculations at the B3LYP/6-31G* level of theory. We demonstrate that the loss of H(2)O from the amide carbonyl group of protonated N-acetyltryptophan O-methyl ester occurs via a "side-chain-backbone" neighboring group reaction to yield a protonated carboline derivative. In contrast, the loss of water from the O-methyl ester of protonated oxindolylalanine results in the formation of a tricyclic structure by "backbone-side-chain" nucleophilic attack from the amino nitrogen to the C2 position of the indole ring. The O-methyl ester of protonated N-acetyloxindolylalanine was found to dissociate via the loss of water from both possible sites, i.e. from the side-chain indolyl oxygen and the backbone amide carbonyl group. An estimate of the relative preference for water loss from each site was obtained from the abundances of product ions formed from MS(3) analysis of regiospecifically labeled derivatives of N-acetyloxindolylalanine, and from the results of molecular orbital calculations. These studies indicate the absence of a characteristic 'signature' ion or neutral loss for peptides containing oxindolylalanine residues under low-energy ion trap CID conditions.