Molecular properties from conformational ensembles. 1. Dipole moments of molecules with multiple internal rotations

We present a novel method for constructing the stable conformational space of small molecules with many rotatable bonds that uses our iterative stochastic elimination (ISE) algorithm, a robust stochastic search method capable of finding ensembles of best solutions for large combinatorial problems. To validate the method, we show that ISE reproduces the best conformers found in a fully exhaustive search, as well as compare computed dipole moments to experimental values, based on molecular ensembles and their Boltzmann distributions. Results were also compared to the alternative molecular dynamics and simulated annealing methods. Our results clarify that many low energy conformations may be required to reproduce molecular properties, while single low energy conformers or ensembles of low energy conformers cannot account for the experimental properties of flexible molecules. Whereas ISE well reproduces conformations that are not separated by very large energy barriers, it has not been successful in reproducing conformations of strained molecules.     

Ab initio Hartree-Fock and density functional theory study on molecular structures, energies, and vibrational frequencies of conformations of 2-hydroxy-3-nitropyridine and 3-hydroxy-2-nitropyridine

The optimized molecular structures, vibrational frequencies and corresponding vibrational assignments of conformations of 2-hydroxy-3-nitropyridine and 3-hydroxy-2-nitropyridine molecules have been investigated using ab initio Hartree-Fock (HF) and density functional theory (B3LYP) methods with 6-311++G(d,p) basis set. The comparison of the experimental and calculated spectra of the molecules have shown that they exist in two conformations with the two OH bond angles (110 degrees and 250 degrees ) respective to the CO bond in the ground state and their energy curves having two minimums have been drawn.     

A new method for fast and accurate derivation of molecular conformations

During molecular simulations, three-dimensional conformations of biomolecules are calculated from the values of their bond angles, bond lengths, and torsional angles. In this paper we study how to efficiently derive three-dimensional molecular conformations from the values of torsional angles. This case is of broad interest as torsional angles greatly affect molecular shape and are always taken into account during simulations. We first review two widely used methods for deriving molecular conformations, the simple rotations scheme and the Denavit-Hartenberg local frames method. We discuss their disadvantages which include extensive bookkeeping, accumulation of numerical errors, and redundancies in the local frames used. Then we introduce a new, fast, and accurate method called the atomgroup local frames method. This new method not only eliminates the disadvantages of earlier approaches but also provides lazy evaluation of atom positions and reduces the computational cost. Our method is especially useful in applications where many conformations are generated or updated such as in energy minimization and conformational search.     

酰基吡唑啉酮化合物的计算、合成和抑菌活性

利用Sybyl 6.3软件包中的分子力学程序对2个4-酰基-2-吡唑啉-5-酮类化合物分子进行了系统研究,得到了它们的最低能量构象,发现它们具有适宜的配位环境,可能是一类新型的三齿配体。 合成和表征了1-苯基-3-甲基-4-(邻-甲酰基苯甲酸)-2-吡唑啉-5-酮 和1-苯基-3-甲基-4-(顺-乙酰基丁烯酸)-2-吡唑啉-5-酮化合物。 抑菌活性测定结果表明它们对革兰氏阳性和阴性菌均有较好的抑制作用,为进一步研究其配合物的生理活性提供了参考依据。     

Comparison of double-ended transition state search methods

While a variety of double-ended transition state search methods have been developed, their relative performance in characterizing complex multistep pathways between structurally disparate molecular conformations remains unclear. Three such methods (doubly-nudged elastic band, a string method, and a growing string method) are compared for a series of benchmarks ranging from permutational isomerizations of the seven-atom Lennard-Jones cluster (LJ(7)) to highly cooperative LJ(38) and LJ(75) rearrangements, and the folding pathways of two peptides. A database of short paths between LJ(13) local minima is used to explore the effects of parameters and suggest reasonable default values. Each double-ended method was employed within the framework of a missing connection network flow algorithm to construct more complicated multistep pathways. We find that in our implementation none of the three methods definitively outperforms the others, and that their relative effectiveness is strongly system and parameter dependent.     

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.     

Protein structure prediction aided by geometrical and probabilistic constraints

Database-assisted ab initio protein structure prediction methods have exhibited considerable promise in the recent past, with several implementations being successful in community-wide experiments (CASP). We have employed combinatorial optimization techniques toward solving the protein structure prediction problem. A Monte Carlo minimization algorithm has been employed on a constrained search space to identify minimum energy configurations. The search space is constrained by using radius of gyration cutoffs, the loop backbone dihedral probability distributions, and various secondary structure packing conformations. Simulations have been carried out on several sequences and 1000 conformations have been initially generated. Of these, 50 best candidates have then been selected as probable conformations. The search for the optimum has been simplified by incorporating various geometrical constraints on secondary structural elements using distance restraint potential functions. The advantages of the reported methodology are its simplicity, and modifiability to include other geometric and probabilistic restraints.     

On searching in, sampling of, and dynamically moving through conformational space of biomolecular systems: A review

Methods to search for low-energy conformations, to generate a Boltzmann-weighted ensemble of configurations, or to generate classical-dynamical trajectories for molecular systems in the condensed liquid phase are briefly reviewed with an eye to application to biomolecular systems. After having chosen the degrees of freedom and method to generate molecular configurations, the efficiency of the search or sampling can be enhanced in various ways: (i) efficient calculation of the energy function and forces, (ii) application of a plethora of search enhancement techniques, (iii) use of a biasing potential energy term, and (iv) guiding the sampling using a reaction or transition pathway. The overview of the available methods should help the reader to choose the combination that is most suitable for the biomolecular system, degrees of freedom, interaction function, and molecular or thermodynamic properties of interest.     

Computational protocols for prediction of solute NMR relative chemical shifts. a case study of L-tryptophan in aqueous solution

In this study, we have applied two different spanning protocols for obtaining the molecular conformations of L-tryptophan in aqueous solution, namely a molecular dynamics simulation and a molecular mechanics conformational search with subsequent geometry re-optimization of the stable conformers using a quantum mechanically based method. These spanning protocols represent standard ways of obtaining a set of conformations on which NMR calculations may be performed. The results stemming from the solute-solvent configurations extracted from the MD simulation at 300 K are found to be inferior to the results stemming from the conformations extracted from the MM conformational search in terms of replicating an experimental reference as well as in achieving the correct sequence of the NMR relative chemical shifts of L-tryptophan in aqueous solution. We find this to be due to missing conformations visited during the molecular dynamics run as well as inaccuracies in geometrical parameters generated from the classical molecular dynamics simulations.     

Conformational analysis of 14- and 16-membered unsaturated oxa, thia, and selena macrocyclic ligands

Crystallographic results retrieved from the Cambridge Structural Database (CSD) have been used to perform systematic conformational analyses of the free and metal-coordinated ligands 1,4,8,11-tetraoxa/thia/selenacyclotetradecane, 1,5,9,13-tetraoxa/thia/selenacyclohexadecane, and their derivatives. Conformational classifications established using symmetry-modified Jarvis-Patrick cluster analysis have been displayed in torsional space by principal component analysis (PCA) plots. Relative molecular mechanics energies of free macrocycles in the observed conformations are compared with the cluster populations and the effect of metal coordination investigated. While the conformations of the free macrocycles in the solid state are those with low molecular mechanics energies, substantial reorganization is required to produce tetradentate chelating conformers, particularly with S donors. With the latter, several such conformers occur, while one predominates with O donors for both macrocycles, in agreement with molecular mechanics predictions. Factors influencing the conformations adopted are discussed.     

An exploration of a novel strategy for superposing several flexible molecules

Summary This paper describes a computational strategy for the superposition of a set of flexible molecules. The combinatorial problems of searching conformational space and molecular matching are reduced drastically by the combined use of simulated annealing methods and cluster analysis. For each molecule, the global minimum of the conformational energy is determined by annealing and the search trajectory is retained in a history file. All the significantly different low-energy conformations are extracted by cluster analysis of data in the history file. Each pair of molecules, in each of their significantly different conformations, is then matched by simulated annealing, using the difference-distance matrix as the objective function. A set of match statistics is then obtained, from which the best match taken from all different conformations can be found. The molecules are then superposed either by reference to a base molecule or by a consensus method. This strategy ensures that as wide a range of conformations as possible is considered, but at the same time that the smallest number of significantly different conformations is used. The method has been tested on a set of six angiotensin II antagonists with between 7–11 rotatable bonds.     

Possible conformations of the seven-membered rings and relative stabilities of the C10 epimers oftrans,trans-,trans,cis-, andcis,trans-linked pseudoguai-11(13)-EN-8,12-olides

A calculation has been made by the molecular mechanics method of the CIO-epimers of trans,trans-, trans,cis-, and cis,trans-linked pseudoguai-11(13)-en-8,12-olides. Possible conformations in these compounds and probabilities of their realization have been determined.     

Theoretical investigation of the coexistence of α and β-nitric acid trihydrates (NAT) molecular conformations

In this letter we document the possibility of the existence of a second molecular configuration for nitric acid trihydrates. Density functional theory (DFT) methods have been used for studying the nitric acid trihydrates α and β-NAT conformations, their spectroscopic and thermodynamics properties and dipole moments have been calculated. This study describes the gas–solid phase transition of the NAT and provides two possible pathways for the molecular structure transformation between α and β-NAT.     

The relationship between peptide plane rotation (PPR) and similar conformations

Currently, several energy functions and conformational search methods have been developed that are based on the observed distribution of phi and psi angles in protein structures. The definition of phi and psi angles is directly related to the orientation of the peptide plane (CA CO NH CA). Starting from one conformation and rotating a single peptide plane, the angles psi for one residue and phi for the consecutive residue that are linked by the peptide plane, display a continuous range of values within one global conformation. When peptide plane rotation is analyzed in several different conformations generated from a restricted conformation database, a large number of these conformations are related. Based on these observations, a new simplified all-atom representation for protein folding simulations is presented where only one torsion angle variable is required for each residue. The underlying theme of this article is that conformational search methods using phi and psi torsion space, search through many redundant conformations. These conformations are related by anticorrelated torsion changes of peptide plane rotations. ©1999 John Wiley & Sons, Inc. J Comput Chem 20: 947–955, 1999     

Study of potassium O,O'-Dibutyldithiophosphate combining DFT, 31P CP/MAS NMR and infrared spectroscopy

Dithiophosphates are used in many different industrial applications. To explain their functions and properties in these applications, a fundamental understanding on a molecular level is needed. Potassium O, O'-Dibutyldithiophosphate and its anion have been investigated by means of a combination of DFT and (31)P CP/MAS NMR and infrared spectroscopy. Several low-energy conformations were studied by DFT. Three different conformations with significantly different torsion angles of the O-C bond relative to the O-P-O plane were selected for further studies of infrared frequencies and (31)P NMR chemical-shift tensors. A good agreement between theoretical and experimental results was obtained, especially when the IR spectra or (31)P chemical shift tensor parameters of all three conformations were added, indicating that, because of the low energy difference between the conformations, the molecules are rapidly fluctuating between them.     

On the conformational structure of cyclohexadecane,the corresponding ketone,the 1,9-dione and the corresponding ketals

The title compounds have been prepared, and their structures have been studied using molecular mechanics calculations, together with dipole moment measurements, and thermodynamic and spectroscopic data reported in the literature. It is concluded that the parent hydrocarbon exists to a very large extent in a “square” diamondoid conformation (D_2d symmetry). Other conformations (rectangular, nearly D₂, and puckered, (D_2d) are considerably higher in energy. The ketone consists of a mixture of both square and rectangular conformations with the CO group in several different positions. The monoethylene ketal exists primarily in a conformation which is square, with the ketal group on a corner. The diketone, diketal, and monoketal of the dione are considered in the light of the above, and conclusions are drawn regarding their conformations. The dipole moments of the difunctional compounds are in agreement with the conformational conclusions reached by the calculations.     

A comparison of a direct search method and a genetic algorithm for conformational searching

We present results from the application of two conformational searching methods: genetic algorithms (GA) and direct search methods for finding low energy conformations of organic molecules. GAs are in a class of biologically motivated optimization methods that evolve a population of individuals in which individuals who are more “fit” have a higher probability of surviving into subsequent generations. The parallel direct search method (PDS) is a type of pattern search method that uses an adaptive grid to search for minima. Both methods found energies equal to or lower than the energy of the relaxed crystal structure in all cases, at a relatively small cost in CPU time. We suggest that either method would be a good candidate to find 3-D conformations in a large scale screening application. © 1996 by John Wiley & Sons, Inc.     

Investigation of the dynamics of the conformational transitions of dimeric isoquinoline alkaloids by nmr spectroscopy and molecular mechanics

The conformational dynamics of the dimeric isoquinoline alkaloid turconidine has been investigated by the methods of NMR spectroscopy and molecular mechanics. The main conformational states and the trajectories of their mutual transitions have been established. A high lability of the molecule and the predominant state of the dimer in three of the eight possible conformations with respect to the ether bond have been shown.Translated from Khimiya Prirodnykh Soedinenii, No. 2, pp. 285–289, March-April, 1995. Original article submitted November 28, 1994.     

Protein structure predictions by parallel simulated annealing molecular dynamics using genetic crossover

We propose a conformational search method to find a global minimum energy structure for protein systems. The simulated annealing is a powerful method for local conformational search. On the other hand, the genetic crossover can search the global conformational space. Our method incorporates these attractive features of the simulated annealing and genetic crossover. In the previous works, we have been using the Monte Carlo algorithm for simulated annealing. In the present work, we use the molecular dynamics algorithm instead. To examine the effectiveness of our method, we compared our results with those of the normal simulated annealing molecular dynamics simulations by using an α-helical miniprotein. We used genetic two-point crossover here. The conformations, which have lower energy than those obtained from the conventional simulated annealing, were obtained.     

Hydrogen bonded oligohydrazide foldamers and their recognition for saccharides

This paper describes the synthesis and characterization of the first series of hydrogen bonding-driven hydrazide foldamers and their recognition for alkyl saccharides in chloroform. Oligomers 1, 2-4, 5, 6, and 7, which contain one, two, four, six, or twelve repeated dibenzoyl hydrazide residues, respectively, have been prepared. The rigid and planar conformations of 1 and 2 or 4 have been established with X-ray analysis and (1)H NMR spectroscopy, whereas the folding and helical conformations of 5-7 have been evidenced by the 1D and 2D (1)H NMR and IR spectroscopy and molecular mechanics calculations. Molecular mechanics calculations also revealed that 5, 6, and 7 possess a rigid cavity with size of ca. 10.6 to 11.1 A, and half of the carbonyl groups in the folding conformations are orientated inwardly inside the cavity. (1)H NMR and CD experiments revealed that 5-7 efficiently complex alkylated mono- and disaccharides 32-35 in chloroform. The association constants (K(assoc)) of the complexes have been determined with the (1)H NMR and fluorescent titration methods. The energy-minimized conformation of 6.34 has been obtained with molecular mechanics calculation. The hydrazide-based folding structures described here represent novel examples of hydrogen bonding-driven foldamers that act as artificial receptors for selective molecular recognition.