An ant algorithm for the conformational analysis of flexible molecules

Originally, the ant system was developed for optimization in discrete search spaces such as the traveling salesman problem. We detail our adaptation of the algorithm to optimization in the continuous search space of conformational analysis. The parameters of the algorithm were tuned using a simple test molecule, undecane, and a drug molecule, imatinib. The algorithm is further tested on four more drug or drug-like molecules, on vitamin A and on alanine tetrapeptide.     

Ionization-loss stimulated Raman spectroscopy for conformational probing of flexible molecules

The approach of studying structural and dynamical properties of flexible molecules is of substantial interest, as it allows decoding the shapes and intrinsic properties of isolated molecular constituents, which have an influence on the selectivity and functionality in biological processes. Combining quantum computation methods with double resonance or infrared hole burning techniques, mainly covering hydride stretch vibrations, recently led to great progress in understanding the structure of a variety of biological building blocks. Measurements of spectra in the lower frequency range, with relatively compact and convenient laser sources, still pose major challenges. For this reason, the method of ionization-loss stimulated Raman spectroscopy (ILSRS) has been developed and applied for monitoring the spectral features of the 2-phenylethanol prototype. The bands observed in the Raman spectra of its two conformers uniquely identify their structures and are in accord with anharmonic results obtained by density functional theory calculations. These findings point to future opportunities for ILSRS as a powerful conformational probe and set new standards for detailed interrogation of structure and intra- and inter-molecular interactions.     

Applications of simulated annealing to the conformational analysis of flexible molecules

We describe in this article our solution to the global minimum problem which uses the simulated annealing algorithm of Kirkpatrick. This method is a Metropolis (e-ΔE/kT) Monte Carlo sampling of conformation space with simultaneous constraint of the search by lowering the temperature T so that the search converges on the global minimum. The Anneal-Conformer program has been extensively tested with peptides and organic molecules using either the Amber or MM2 force fields. A history file of the simulated annealing process allows reconstruction of the random walk in conformation space for subsequent examination. Thus plots of distance and dihedral angle changes during the search for the global minimum can be examined to deduce molecular shape and flexibility. A separate program Conf-Gen reads the history file and extracts all low energy conformations visited during the run.     

Influence of vibrational energy flow on isomerization of flexible molecules: incorporating non-Rice-Ramsperger-Kassel-Marcus kinetics in the simulation of dipeptide isomerization

The conformational isomerization of a dipeptide, N-acetyl-tryptophan methyl amide (NATMA), is studied computationally by including important dynamical corrections to Rice-Ramsperger-Kassel-Marcus (RRKM) theory for the transition rate between pairs of isomers. The dynamical corrections arise from incomplete or sluggish vibrational energy flow in the dipeptide, a property suggested by the mode-selective chemistry that has been observed by Dian et al. [J. Chem. Phys. 120, 133 (2004)]. We compute the extent and rate of vibrational energy flow in NATMA quantum mechanically using local random matrix theory, which we then use to correct the RRKM theory rates. The latter rates are then introduced into a master equation to study the population dynamics of the dipeptide. Incomplete or slow vibrational energy flow is found to enhance the conformational selectivity of NATMA over RRKM estimates.     

Cascade of coil-globule conformational transitions of single flexible polyelectrolyte molecules in poor solvent

We show that hydrophobic flexible polyelectrolyte molecules of poly(2-vinylpyridine) and poly(methacryloyloxyethyl dimethylbenzylammonium chloride) are trapped and frozen due to adsorption on the mica surface, and the observed AFM single molecule structures reflect the molecular conformation in solution. An increase of the ionic strength of the solution induces the cascade of abrupt conformational transitions due to the intrachain segregation from elongated coil to compact globule conformations through intermediate pearl necklace-globule conformations with different amounts of beads per chain. The length of the necklaces and the number of beads decrease, while the diameter of beads increases with the increase of ionic strength. Coexistence at the same time of extended coils, necklaces with different amounts of beads, and compact globules indicates the cascade of the first-order-type phase transitions.     

Roles of polymer-solvent complexes in conformational ordering of polymer molecules dispersed in gels

On gelation of crystallizable stereoregular polystyrenes (isotactic IPS and syndiotactic SPS), highly ordered skeletal conformations of particular types characteristic of the respective polymer are formed. The conformational orderliness is evaluated in the absolute scale through quantitative analysis of infrared absorption bands due to the respective ordered conformation. The kinetic and/or thermodynamic behaviors of gelation process and the accompanied conformational ordering are strongly influenced by solvent. Solvent effect on the aggregation state formed in SPS gels has been investigated by means of small angle neutron scattering and discussed in relation to the gelation rate. In some solvents, two different ordering processes, gelation and crystallization, proceed in a competitive way. On the basis of these experimental facts, the roles of specific polymer-solvent interactions in gels are considered from the viewpoint of molecular-level structure.     

Unusual conformational isomerization of oxygen-containing silacyclohexanes

Russian Journal of Organic Chemistry -     

Symmetry of nonrigid molecules and isomerization processes

Starting from the concept of a totally flexible molecule, rigidity is introduced by restricting progressively the allowed permutations. Each nonrigid isomer is then characterized by its Longuet-Higgins group ? of allowed permutation inversions. For a given ? all nonrigidity types correspond to all possible choices of the symmetry subgroup ? of the skeleton. The structure of ? also allows a characterization of the isomerization processes studied in chemical kinetics. In both situations, isomerization mechanisms may be associated with generators and are, in the simplest situations, represented by Cayley graphs.     

ELECT++: Faster conformational search method for docking flexible molecules using molecular similarity

We have developed a program, ELECT++ (Effective LEssening of Conformations by Template molecules in C++), to speed up the conformational search for small flexible molecules using the similar property principle. We apply this principle to molecular shape and, importantly, to molecular flexibility. After molecules in a database are clustered according to flexibility and shape (FCLUST++), additional reagents are generated to screen the conformational space of molecules in each cluster (TEMPLATE++). We call these representative reagents of each cluster template reagents. Template reagents and clustered reagents produce, after reaction, template molecules and clustered molecules, respectively (tREACT++). The conformations of a template molecule are searched in the context of a macromolecular target. Acceptable conformational choices are then applied to all molecules in its cluster, thus effectively biasing conformational space to speed up conformational searches (tSEARCH++). In our incremental search method, it is necessary to calculate the root-mean-square deviations (RMSD) matrix of distances between different conformations of the same molecule to reduce the number of conformations. Instead of calculating the RMSD matrix for all molecules in a cluster, the RMSD matrix of a template molecule is chosen as a reference and applied to all the molecules in its cluster. We demonstrate that FCLUST++ clusters the primary amine reagents from the Available Chemicals Directory (ACD) successfully. The program tSEARCH++ was applied to dihydrofolate reductase with virtual molecules generated by tREACT++ using clustered primary amine reagents. The conformational search by the program tSEARCH++ was about 4.8 times faster than by SEARCH++, with an acceptable range of errors. © 1998 John Wiley & Sons, Inc. J Comput Chem 19: 1834–1852, 1998     

Combination of molecular dynamics method and 3D-RISM theory for conformational sampling of large flexible molecules in solution

We have developed an algorithm for sampling the conformational space of large flexible molecules in solution, which combines the molecular dynamics (MD) method and the three-dimensional reference interaction site model (3D-RISM) theory. The solvent-induced force acting on solute atoms was evaluated as the gradient of the solvation free energy with respect to the solute-atom coordinates. To enhance the computation speed, we have applied a multiple timestep algorithm based on the RESPA (Reversible System Propagator Algorithm) to the combined MD/3D-RISM method. By virtue of the algorithm, one can choose a longer timestep for renewing the solvent-induced force compared with that of the conformational update. To illustrate the present MD/3D-RISM simulation, we applied the method to a model of acetylacetone in aqueous solution. The multiple timestep algorithm succeeded in enhancing the computation speed by 3.4 times for this model case. Acetylacetone possesses an intramolecular hydrogen-bonding capability between the hydroxyl group and the carbonyl oxygen atom, and the molecule is significantly stabilized due to this hydrogen bond, especially in gas phase. The intramolecular hydrogen bond was kept intact during almost entire course of the MD simulation in gas phase, while in the aqueous solutions the bond is disrupted in a significant number of conformations. This result qualitatively agrees with the behavior on a free energy barrier lying upon the process for rotating a torsional degree of freedom of the hydroxyl group, where it is significantly reduced in aqueous solution by a cancellation between the electrostatic interaction and the solvation free energy.     

Experimental determination of conformational isomerization energy thresholds in serotonin

Stimulated emission pumping-population transfer (SEP-PT) and hole-filling (SEP-HF) spectroscopies were used to determine the energy thresholds to isomerization between thirteen reactant-product conformer pairs in the biomolecule serotonin (SERO). Serotonin is a close structural analog of tryptamine (TRA), differing in having a hydroxyl group in the 5 position of the indole ring. A previous spectroscopic study (LeGreve; et al. J. Am. Chem. Soc. 2007, 129 (13), 4028) identified eight conformational isomers of SERO, whose interconversion involves motion of the 3-ethylamine side chain, the 5-OH group, or both. In the cases in which only an ethylamine side chain reorientation occurred, the barriers were found to be similar to, but systematically somewhat smaller than, those in TRA, which has been studied by similar methods (Dian; et al. Science 2004, 303 (5661), 1169; Clarkson; et al. J. Chem. Phys. 2005, 122 (21), Art. No. 214311). In most cases, the experimental thresholds are well reproduced by calculated transition states separating the conformational wells; however, tunneling effects may artificially reduce the thresholds observed for isomerization of SERO(A,Gpy(out)) and SERO(B,Gpy(up)) into SERO(C,Gph(out)). The A --> A' isomerization involving only the OH rotation from anti to syn was found to be 721-761 cm-1, in accordance with the calculated classical barrier. For isomerizations in which the ethylamine side chain reorients as does the OH group, the barriers to isomerization were consistent with sequential rather than concerted motion of both groups. Finally, some evidence for mode-specific effects in the product quantum yields near threshold is presented.     

The dynamics of conformational isomerization in flexible biomolecules. I. Hole-filling spectroscopy of N-acetyl tryptophan methyl amide and N-acetyl tryptophan amide

The conformational isomerization dynamics of N-acetyl tryptophan methyl amide (NATMA) and N-acetyl tryptophan amide (NATA) have been studied using the methods of IR-UV hole-filling spectroscopy (HFS) and IR-induced population transfer spectroscopy (IR-PTS), which were developed for this purpose. Single conformations of these molecules were selectively excited in well-defined NH stretch fundamentals. This excess energy was used to drive conformational isomerization. By carrying out the infrared excitation early in a supersonic expansion, the excited molecules were recooled into their zero-point levels, partially refilling the hole created in the ground state population of one of the conformers, and creating gains in population in other conformers. These changes in population were detected using laser-induced fluorescence downstream in the expansion. In HFS, the IR wavelength is fixed and the UV laser tuned in order to determine where the population went following selective infrared excitation. In IR-PTS, the UV is fixed to monitor the population of a given conformation, and the IR is tuned to record the IR-induced changes in the population of the monitored conformer. Besides demonstrating the capability of the experiment to change the downstream conformational population distribution, the IR-PTS scans were used to extract two quantitative results: (i) The fractional populations of the conformers in the absence of the infrared, and (ii) the isomerization quantum yields for each of the six unique amide NH stretch fundamentals (three conformers each with two amide groups). The method for obtaining quantum yields is described in detail. In both NATMA and NATA, the quantum yields show modest conformational specificity, but only a hint of vibrational mode specificity. The prospects for the hole-filling technique for providing insight into energy flow in large molecules are discussed, leaving a more detailed theoretical modeling to the adjoining paper [Evans et al. J. Chem. Phys. 120, 148 (2004)].     

Near-infrared-induced conformational isomerization of 2-chloroethanol in an argon matrix

The infrared spectrum of 2-chloroethanol in a low-temperature Ar matrix shows marked intensity changes in the O---H and C---Cl stretching reg     

Torsional isomerization of biologically active bicyclic molecules

Rotational isomerization of bipyridines C₅H₄N? C₅H₄N was studied by CNDO /2, PPP /CI , and CNDO /CI methods. It is shown that CNDO /2 overestimates the angle of rotation ? between the pyridine rings ca. two times. The angle ? was determined for 2,2′-bipyridine by means of correlation of the theoretical (CNDO /2) and experimental dipole moment. It was also found from the correlation between the theoretical and experimental UV spectra. It is shown that there is an explicit dependence of the results upon the distance between heteroatoms (PPP /CI ). It has been found that the CNDO /CI method correctly predicts the value of the rotational angles and their sequence in bipyridines.     

Os络合物催化烯烃加氢和异构化反应

综述了近年来锇络合物用于催化烯烃加氢和异构化反应的研究进展。 Os催化剂在H2分子和转移加氢二个方面用于烯烃加氢反应均表现出较高的活性和选择性。因此它有望成为有机合成中的一个强有力的工具。     

Picosecond isomerization in photochromic ruthenium-dimethyl sulfoxide complexes

The complexes [Ru(tpy)(bpy)(dmso)](OSO(2)CF(3))(2) and trans-[Ru(tpy)(pic)(dmso)](PF(6)) (tpy is 2,2':6',2' '-terpyridine, bpy is 2,2'-bipyridine, pic is 2-pyridinecarboxylate, and dmso is dimethyl sulfoxide) were investigated by picosecond transient absorption spectroscopy in order to monitor excited-state intramolecular S-->O isomerization of the bound dmso ligand. For [Ru(tpy)(bpy)(dmso)](2+), global analysis of the spectra reveals changes that are fit by a biexponential decay with time constants of 2.4 +/- 0.2 and 36 +/- 0.2 ps. The first time constant is assigned to relaxation of the S-bonded (3)MLCT excited state. The second time constant represents both excited-state relaxation to ground state and excited-state isomerization to form O-[Ru(tpy)(bpy)(dmso)](2+). In conjunction with the S-->O isomerization quantum yield (Phi(S)(-->)(O) = 0.024), isomerization of [Ru(tpy)(bpy)(dmso)](2+) occurs with a time constant of 1.5 ns. For trans-[Ru(tpy)(pic)(dmso)](+), global analysis of the transient spectra reveals time constants of 3.6 +/- 0.2 and 118 +/- 2 ps associated with these two processes. In conjunction with the S-->O isomerization quantum yield (Phi(S)(-->)(O) = 0.25), isomerization of trans-[Ru(tpy)(pic)(dmso)](+) occurs with a time constant of 480 ps. In both cases, the thermally relaxed excited states are assigned as terpyridine-localized (3)MLCT states. Electronic state diagrams are compiled employing these data as well as electrochemical, absorption, and emission data to describe the reactivity of these complexes. The data illustrate that rapid bond-breaking and bond-making reactions can occur from (3)MLCT excited states formed from visible light irradiation.     

Pathways of conformational isomerization of 4-phenyl-1,3-dioxane

A study of the potential energy surface of 4-phenyl-1,3-dioxane by non-empirical quantum chemical RHF/STO-3G and 6-31G(d) approximations reveals six energetically inequivalent pathways of conformational isomerization of equatorial and axial chair forms.