Searching and optimizing structure ensembles for complex flexible sugars

NMR restrictions are suitable to specify the geometry of a molecule when a single well-defined global free energy minimum exists that is significantly lower than other local minima. Carbohydrates are quite flexible, and therefore, NMR observables do not always correlate with a single conformer but instead with an ensemble of low free energy conformers that can be accessed by thermal fluctuations. In this communication, we describe a novel procedure to identify and weight the contribution to the ensemble of local minima conformers based on comparison to residual dipolar couplings (RDCs) or other NMR observables, such as scalar couplings. A genetic algorithm is implemented to globally minimize the R factor comparing calculated RDCs to experiment. This is done by optimizing the weights of different conformers derived from the exhaustive local minima conformational search program, fast sugar structure prediction software (FSPS). We apply this framework to six human milk sugars, LND-1, LNF-1, LNF-2, LNF-3, LNnT, and LNT, and are able to determine corresponding population weights for the ensemble of conformers. Interestingly, our results indicate that in all cases the RDCs can be well represented by only a few most important conformers. This confirms that several, but not all of the glycosidic linkages in histo-blood group "epitopes" are quite rigid.     

Spanning set of silica cluster isomer topologies from QTAIM

Structural and chemical properties of the building block of silica nanowires, (SiO(2))(6), are investigated with the theory of atoms and molecules (QTAIM). Twenty-five conformers have been analyzed, ten of which have not been reported before. We extend the silica (SiO(2))(6) topology phase space using QTAIM; the Poincaré-Hopf topological sum rules are applied and used to identify the spanning set of topologies, and this includes finding eight new distinct topologies that satisfy the Poincaré-Hopf relation. A simple phase diagram of the solutions of the Poincaré-Hopf relation is created with the aid of a new classification scheme to determine the boundary between topological stability and instability. Sum rules are then found to be applicable to any set of isomers. We determine that O-O bonding interactions exist for the silica (SiO(2))(6) conformers in regions where the energy surface is flattest. In addition, we identify unstable local minima in the topology of the charge density in order to further compare conformer instabilities. We quantify the dimensionality of a molecule using the Poincaré-Hopf relation instead of Euclidean geometry. This quantum topological definition of geometry shows that the four most energetically stable (SiO(2))(6) conformers are quantified as two-dimensional within the new quantum topology.     

Preferred conformers and photochemical (lambda > 200 nm) reactivity of serine and 3,3-dideutero-serine in the neutral form

A systematic investigation of the conformational potential energy surface of neutral serine [HOCH2CHNH2COOH] and 3,3-dideutero-serine [HOCD2CHNH2COOH] was undertaken, revealing the existence of 61 different minima. The structures and vibrational spectra of the most stable conformers, which were estimated to have relative energies within 7 kJ mol(-1) and account for ca. 93% of the total conformational population at room temperature, were calculated at both the MP2 and DFT/BLYP levels of theory with the 6-311++G(d,p) basis-set and used to interpret the spectroscopic data obtained for the compounds isolated in low-temperature inert matrixes. The assignment of the main spectral infrared features observed in the range 4000-400 cm(-1) to the most stable conformers of serine was undertaken. In addition, UV irradiation (lambda > 200 nm) of the matrix-isolated compounds was also performed, leading to decarboxylation, which was found to be strongly dependent on the conformation assumed by the reactant molecule.     

Ab initio Study on Ionization Energies of 3-Amino-1-propanol

Fourteen conformers of 3-amino-1-propanol as the minima on the potential energy surface are examined at the MP2/6-311++G** level. Their relative energies calculated at B3LYP, MP3 and MP4 levels of theory indicated that two most stable conformers display the in-tramolecular OH…N hydrogen bonds. The vertical ionization energies of these conformers calculated with ab initio electron propagator theory in the P3/aug-cc-pVTZ approximation are in agreement with experimental data from photoelectron spectroscopy. Natural bond orbital analyses were used to explain the differences of IEs of the highest occupied molec-ular ortibal of conformers. Combined with statistical mechanics principles, conformational distributions at various temperatures are obtained and the temperature dependence of pho-toelectron spectra is interpreted.     

Camdas: An automated conformational analysis system using molecular dynamics

We present an automated conformational analysis program, CAMDAS (Conformational Analyzer with Molecular Dynamics And Sampling). CAMDAS performs molecular dynamics (MD) calculations for a target molecule and samples conformers from the trajectory of the MD. The program then evaluates the similarities between each of the sampled conformers in terms of the root- mean-square deviations of the atomic positions, clusters similar conformers, and finally prints out the clustered conformers. This MD-based conformational analysis is a broadly used method, and CAMDAS is intended to provide a convenient framework for the method. CAMDAS has the ability to find the representative conformers automatically from an arbitrarily given structure of the molecule. The accuracy of the program was examined using N- acetylalanine-N-methylamide, and the obtained result was consistent with that of the systematic search method. In the test calculation of cyclodecane, CAMDAS could identify most of the known conformers and their conformational enantiomers by examining only 5000 conformers. In addition, the potential-scaled method, which we have developed previously as an accelerating technique for MD, could find two additional conformers of cyclodecane that have not been reported. CAMDAS presents a convenient way to find the energetically possible conformers of a molecule, which is needed especially in the early stage of drug design.     

The molecular conformation of butyl acrylate – A vibrational spectroscopy and computational study

Potential energy curves allow us to probe the flexibility of butyl acrylate molecule. We found two local minima belonging to C_s symmetry for s-cis and s-trans conformations. The optimized geometries at RHF and DFT using extended basis set are in good agreement with electron diffraction data of methyl acrylate, and the s-cis is the most stable form. Ab initio calculations of harmonic frequencies at the fully optimized geometries of the planar s-cis and s-trans conformers have been performed at the DFT//B3LYP/6-311+G** level of theory. Scaled harmonic frequencies and potential energy distribution are used for the assignment of the experimental IR and Raman bands. We noticed an acceptable qualitative agreement between the experimental and the computed spectra.     

Solvent-assisted intramolecular proton transfer in thiopurinol: application of M06-2X functional

All available conformers of tisopurine as an important pharmaceutical molecule are optimized and frequency calculations calculated at M06-2X/6-311++G(2d,2p) level of theory. These conformers are classified in 22 different tautomers, tautomer Z showing the most stable tautomer in the gas phase. Effects of four different solvents on the most stable conformer of each tautomer is calculated. Solvents cause stabilization of all conformers and relative solvent stabilization is as follows: water > DMSO > acetone > toluene. Energy profile for such stabilization is illustrated and mechanism of proton transfer studied at the same level of theory. Solvent-assisted proton transfer performed when water and methanol used as solvents. Results indicate that explicit solvent effect has much more stabilization on tautomerization processes compared to implicit solvent effect.     

Conformational transitions of glycine induced by vibrational excitation of the O-H stretch

Vibrational energy flow and conformational transitions following excitation of the OH stretching mode of the most stable conformer of glycine are studied by classical trajectories. "On the fly" simulations with the PM3 semiempirical electronic structure method for the potential surface are used. Initial conditions are selected to correspond to the ν=1 excitation of the OH stretch. The main findings are: (1) An an equilibrium-like ratio is established between the populations of the 3 lowest-lying conformers after about 10 picoseconds. (2) There is a high probability throughout the 150 ps of the simulations for finding the molecule in geometries far from the equilibrium structures of the lowest-energy conformers. (3) Energy from the initial excited OH (ν=1) stretch flows preferentially to 5 other vibrational modes, including the bending motion of the H atom. (4) RRK theory yields conformational transition rates that deviate substantially from the classical trajectory results. Possible implication of these results for vibrational energy flow and conformational transitions in small biological molecules are discussed.     

Structure of Isolated 1,4-Butanediol: Combination of MP2 Calculations, NBO Analysis, and Matrix-Isolation Infrared Spectroscopy

Theoretical calculations at the MP2 level, NBO and AIM analysis, and matrix-isolation infrared spectroscopy have been used to investigate the structure of the isolated molecule of 1,4-butanediol (1,4-BDO). Sixty-five structures were found to be minima on the potential energy surface, and the three most stable forms are characterized by a folded backbone conformation leading to the formation of an intramolecular H-bond. To better characterize the intramolecular interactions and particularly the hydrogen bonds, natural bond orbital analysis (NBO) was performed for the four most stable conformers, and was further complemented with an atoms-in-molecules (AIM) topological analysis. Infrared spectra of 1,4-BDO isolated in low-temperature argon and xenon matrixes show a good agreement with a population-weighted mean theoretical spectrum, and the spectral features of the conformers expected to be trapped in the matrixes were observed experimentally. Annealing the xenon matrix from 20 to 60 K resulted in significant spectral changes, which were interpreted based on the barriers to intramolecular rotation. An estimation of the intramolecular hydrogen bond energy was carried out following three different methodologies.     

Probing furanose ring conformation by gas-phase computational methods: energy profile and structural parameters in methyl beta-D-arabinofuranoside as a function of ring conformation

The potential energy surface of methyl beta-D-arabinofuranoside (3) has been studied by ab initio molecular orbital (HF/6-31G) and density functional theory (B3LYP/6-31G) calculations via minimization of the 10 possible envelope conformers. The partial potential energy surface identified that the global minimum and lowest energy northern conformer was E(2). In the HF calculations, (2)E was the most stable southern conformer, while the density functional theory methods identified (4)E as the local minimum in this hemisphere. Additional calculations at higher levels of theory showed that the B3LYP-derived energies of many of the envelope conformers of 3 are dependent upon the basis set used. It has also been demonstrated that B3LYP/6-31+G//B3LYP/6-31G single point energies are essentially the same as those obtained from full geometry optimizations at the B3LYP/6-31+G level. The northern and southern minima of the B3LYP/6-31+G surface are, respectively, the E(2) and (2)E conformers. The B3LYP/6-31G geometries were used to study the relationship between ring conformation and various structural parameters including bond angles, dihedral angles, bond lengths, and interatomic distances.     

Molecular Self‐Recognition: Rotational Spectra of the Dimeric 2‐Fluoroethanol Conformers

Fluoroalcohols show competitive formation of intra‐ and intermolecular hydrogen bonds, a property that may be crucial for the protein‐altering process in a fluoroalcohol/water solution. In this study, we examine the intra‐ and intermolecular interactions of 2‐fluoroethanol (FE) in its dimeric conformers by using rotational spectroscopy and ab initio calculations. Three pairs of homo‐ and heterochiral dimeric FE conformers are predicted to be local minima at the MP2/6‐311++G(d,p) level of theory. They are solely made of the slightly distorted most stable G+g?/G?g+ FE monomer units. Jet‐cooled rotational spectra of four out of the six predicted dimeric conformers were observed and unambiguously assigned for the first time. All four observed dimeric conformers have compact geometries in which the fluoromethyl group of the acceptor tilts towards the donor and ensures a large contact area. Experimentally, the insertion of the O? H group of one FE subunit into the intramolecular O? H???F bond of the other was found to lead to a higher stabilisation than the pure association through an intermolecular O? H???O? H link. The hetero‐ and homochiral combinations were observed to be preferred in the inserted and the associated dimeric conformers, respectively. The experimental rotational constants and the stability ordering are compared with the ab initio calculations at the MP2 level with the 6‐311++G(d,p) and aug‐cc‐pVTZ basis sets. The effects of fluorination and the competing inter‐ and intramolecular hydrogen bonds on the stability of the dimeric FE conformers are discussed.     

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.     

Do computed crystal structures of nonpolar molecules depend on the electrostatic interactions? The case of tetracene

We present a strategy for comparing the global properties of competing potential models. By systematically sampling the potential energy surface of crystalline tetracene, we assess how the number, energy and structure of its minima are modified by switching on (or off) the Coulombic interactions. The increased complexity of the Coulombic potential leads to a more "rugged" potential energy surface with a larger number of minima, but the effect is not large. In fact, we find a subset of minima stable only in presence of the Coulombic interactions, a smaller subset stable only in their absence, and a large majority stable in both cases. Among these, there is a very good, but not perfect, correlation between the energies and the structures computed with and without the electrostatic interactions. Although electrostatic interactions play a role even in a rigid nonpolar molecule such as tetracene, they are not as crucial as often believed, because altering the electrostatic model (or switching it off completely) leads, in most cases, to equivalent results.     

Structural and vibrational characterisation of 3-amino-1-propanol a concerted SCF-MO ab initio, Raman and infrared (matrix isolation and liquid phase) spectroscopy study

Results obtained for the isolated and liquid 3-amino-1-propanol by a concerted molecular orbital and vibrational spectroscopic approach are reported. The relative energies and both structural and vibrational data of the different conformers of the studied compound were calculated using the extended 6-31G* basis set both at the HF-SCF and MP2 ab initio levels of theory and the theoretical results used to interpret Raman and infrared experimental data. In the gaseous phase and for the molecule isolated in an Argon matrix, monomeric 3-amino-1-propanol exists as a mixture of conformers, the first and second lowest energy forms corresponding to conformers which exhibit an intramolecular OH-N hydrogen bond (forms I and II). On the other hand, in the pure liquid, where intermolecular H-bonding occurs, the monomeric unit within the aggregates assumes a conformation similar to that of the third most stable form found for the isolated molecule situation (form III), which is characterised by having a weak intramolecular NH-O bond. The experimental data obtained for the pure liquid also reveals the presence of monomeric form I in this phase, a result that is in consonance with the strongly stabilizing OH-N intramolecular hydrogen bond that is present in this conformer.     

Spectroscopic evidence for gas-phase formation of successive beta-turns in a three-residue peptide chain

We report the first gas-phase spectroscopic study of a three-residue model of a peptide chain, Ac-Phe-Gly-Gly-NH2 (Ac = acetyl), using the IR/UV double resonance technique. The existence of at least five different conformers under supersonic expansion conditions is established, most of them exhibiting rather strong intramolecular H-bonds. One of the most populated conformers, however, exhibits a different H-bonding network characterized by two weak H-bonds. Comparison of the amide A and I/II experimental data with density functional theory calculations carried out on a series of selected conformations enables us to assign this conformer to two successive beta-turns along the peptide chain, the two H-bonds being of C10 type, i.e., each of them closing a 10-atom ring in the molecule. The corresponding form is found to be more stable than the 310 helix secondary structure (not observed), presumably because of specific effects due to the glycine residues.     

Potential energy and free energy landscapes

Familiar concepts for small molecules may require reinterpretation for larger systems. For example, rearrangements between geometrical isomers are usually considered in terms of transitions between the corresponding local minima on the underlying potential energy surface, V. However, transitions between bulk phases such as solid and liquid, or between the denatured and native states of a protein, are normally addressed in terms of free energy minima. To reestablish a connection with the potential energy surface we must think in terms of representative samples of local minima of V, from which a free energy surface is projected by averaging over most of the coordinates. The present contribution outlines how this connection can be developed into a tool for quantitative calculations. In particular, stepping between the local minima of V provides powerful methods for locating the global potential energy minimum, and for calculating global thermodynamic properties. When the transition states that link local minima are also sampled we can exploit statistical rate theory to obtain insight into global dynamics and rare events. Visualizing the potential energy landscape helps to explain how the network of local minima and transition states determines properties such as heat capacity features, which signify transitions between free energy minima. The organization of the landscape also reveals how certain systems can reliably locate particular structures on the experimental time scale from among an exponentially large number of local minima. Such directed searches not only enable proteins to overcome Levinthal's paradox but may also underlie the formation of "magic numbers" in molecular beams, the self-assembly of macromolecular structures, and crystallization.     

Structure and dynamics of the aniline-argon complex as derived from its potential energy surface

The structure and dynamics of the van der Waals (vdW) complex of aniline (An) with argon (Ar) are studied using ab initio methods. The inversion potential of the aniline-argon (AnAr) complex perturbed by the weak vdW interaction is calculated taking into account subtle corrections from the zero-point energy of the vdW modes and from the frequency shifts of the An normal modes modified by the complexation. The intermolecular potential energy surface (PES) of the AnAr complex is determined by performing a large-scale computation of the interaction energy and the fitting of the analytical many-body expansion to the set of single-point interaction energies. The PES determined shows two deep local minima corresponding to the anti and syn AnAr conformers. The difference in the energies of these two minima is only 15 cm-1, but it is sufficient to localize the inversion wave functions and to form the two conformers. In the conformer anti (syn) of lower (higher) energy, Ar is bound to the An ring opposite (adjacent) the amino-hydrogens. In the additional local minima higher in energy, Ar approaches the aniline ring between the C-H bonds near its plane. An additional local minimum is located opposite of nitrogen between the two N-H bonds. The high-energy minima are, however, too flat to form stable conformers. The perturbation of the interaction of Ar with the phenyl ring by the NH2 group is described by the vdW hole, which is responsible for unusually strong intermode mixing in the excited intermolecular vibrational states. The analysis of these states calculated for the ground (S0) as well as the first excited electronic state (S1) resolves difficulties faced earlier with the assignment of the observed vibronic bands of AnAr.     

A full conformational space analysis of bilirubin

Ab initio methods were utilized in a gas-phase systematic conformational search of bilirubin conformers. The whole molecule was divided into four fragments. Most stable conformers of them were employed to build 196 conformers of the complete bilirubin molecule. Initial geometries were optimized using HF/3-21G level of theory and the minimum energy conformers were then reoptimized at B3LYP/6-31G(d) level. Ridge-tile conformer was the most stable one, in perfect agreement with X-ray data. We found that while tetrapyrrole backbone shows some flexibility, propionic acid side chains have a greater influence in bilirubin conformation because they can interact through different hydrogen bond patterns with the backbone and between them.     

Inherent structures of crystalline tetracene

We have systematically sampled the potential energy surface of crystalline tetracene to identify its local minima. These minima represent all possible stable configurations and constitute the "inherent structures" of the system. The crystal is described in terms of rigid molecules with Coulombic and atom-atom interactions. Hundreds of distinct minima are identified, mostly belonging to the space groups P (triclinic) and P2(1)/c (monoclinic), with a variety of structural arrangements. The deepest minimum corresponds to the high temperature-low pressure polymorph. This is the only polymorph with a completely described X-ray structure, which is satisfactorily described by the calculations. The next deep minimum is likely to correspond to the low temperature-high pressure polymorph, which has been experimentally identified but not yet fully described.