Stereoelectronic tuning of the structure and stability of the trp cage miniprotein

Proline residues are critical structural elements in proteins, defining turns, loops, secondary structure boundaries, and polyproline helices. Control of proline conformation therefore may be used to define protein structure and stability. 4-Substituted proline derivatives may be used to control proline ring pucker, which correlates with protein main chain conformation. To examine the use of proline conformational restriction to tune globular protein stability, a series of peptides derived from the trp cage miniprotein was synthesized. Proline at residue 12 of the trp cage miniprotein, which adopts a Cgamma-exo ring pucker in the NMR structure, was replaced with 4-substituted proline derivatives, including 4R derivatives favoring a Cgamma-exo ring pucker and 4S derivatives favoring a Cgamma-endo ring pucker. Eight trp cage peptides were synthesized, five of which included residues that are not commercially available, without requiring any solution phase chemistry. Analysis of the trp cage peptides by circular dichroism and NMR indicated that the structure and stability of the trp cage miniprotein was controllable based on the conformational bias of the proline derivative. Replacement of Pro12 with 4S-substituted proline derivatives that favor the Cgamma-endo ring pucker destabilized the trp cage, while replacement of Pro12 with 4R-substituted proline derivatives that favor a Cgamma-exo ring pucker resulted in increased alpha-helicity and thermal stability of the trp cage. The most stable trp cage derivatives contained benzoates of 4R-hydroxyproline, which also exhibited the most pronounced stereoelectronic effects in TYProxN model peptides. Overall, the stability of the trp cage was tunable by over 50 degrees C depending on the identity of the proline side chain at residue 12.     

Conformers of gaseous proline

Accurate geometries, relative energies, rotational and quartic centrifugal distortion constants, dipole moments, harmonic vibrational frequencies, and infrared intensities were determined from ab initio electronic structure calculations for eighteen conformers of the neutral form of the amino acid L-proline. Only four conformers have notable population at low and moderate temperature. The second most stable conformer is only 2+/-2 kJ mol(-1) above the global minimum, while the third and fourth conformers are nearly degenerate and have an excess energy of 7+/-2 kJ mol(-1) relative to the global minimum. All four conformers have one hydrogen bond: N.HO in the lower energy pair of conformers, and NH.O in the higher energy pair of conformers. The conformer pairs differ only in their ring puckering. The relative energies of the conformers include corrections for valence electron correlation, extrapolated to the complete basis set limit, as well as core correlation and relativistic effects. Structural features of the pyrrolidine ring of proline are discussed by using the concept of pseudorotation. The accurate rotational and quartic centrifugal distortion constants as well as the vibrational frequencies and infrared intensities should aid identification and characterization of the conformers of L-proline by rotational and vibrational spectroscopy, respectively. Bonding features of L-proline, especially intramolecular hydrogen bonds, were investigated by the atoms-in-molecules (AIM) technique.     

DFT‐Aided Vibrational Circular Dichroism Spectroscopy Study of (−)‐S‐cotinine

The implementation of a strategy comprising the use of vibrational circular dichroism spectroscopy and DFT calculations allows determination of the solution‐state conformational distribution in (?)‐S‐cotinine, giving further proof of the extra conformer‐discriminating potential of this experimental technique, which may offer unique molecular fingerprints of subtly dissimilar molecular conformers of chiral samples. Natural bond orbital electronic structure calculations of the rotational barrier height between the two main conformers of the species indicate that hyperconjugative effects are the key force governing the conformational equilibrium. The negligible effect of the solvent’s polarity over both structure and conformational energy profile supports this result.     

Conformational equilibrium of bis(trifluoromethanesulfonyl) imide anion of a room-temperature ionic liquid: Raman spectroscopic study and DFT calculations

The structure of bis(trifluoromethanesulfonyl) imide (TFSI-) in the liquid state has been studied by means of Raman spectroscopy and DFT calculations. Raman spectra of 1-ethyl-3-methylimidazolium (EMI+) TFSI- show relatively strong bands arising from TFSI- at about 398 and 407 cm(-1). Interestingly, the 407 cm(-1) band, relative to the 398 cm(-1) one, is appreciably intensified with raising temperature, suggesting that an equilibrium is established between TFSI- conformers in the liquid state. According to DFT calculations followed by normal frequency analyses, two conformers of C2 and C1 symmetry, respectively, constitute global and local minima, with an energy difference 2.2-3.3 kJ mol(-1). The wagging omega-SO2 vibration appears at 396 and 430 cm(-1) for the C1 conformer and at 387 and 402 cm(-1) for the C2 one. Observed Raman spectra over the range 380-440 cm(-1) were deconvoluted to extract intrinsic bands of TFSI- conformers, and the enthalpy of conformational change from C2 to C1 was evaluated. The enthalpy value is in good agreement with that obtained by theoretical calculations. We thus conclude that a conformational equilibrium is established between the C1 and C2 conformers of TFSI- in the liquid EMI+TFSI-, and the C2 conformer is more favorable than the C1 one.     

Molecular dynamics simulations of a nucleoside analogue of 1,4-dihydro-4-oxoquinoline-3-carboxylic acid synthesized as a potential antiviral agent: Conformational studies in vacuum and in water

Conformational analysis of nucleosides may have direct applications to the structure–activity relationship (SAR) studies and in the design of new drug candidates. Although conformational analysis may be accessed in many different ways, in this work it was performed using molecular dynamics (MD) simulation in order to study the dynamic behavior of a nucleoside derivative of 1,4-dihydro-4-oxoquinoline-3-carboxylic acid, synthesized by our group as a potential antiviral agent. The MD simulation was carried out during 10 ns in vacuum and in a box of water at two different temperatures (i.e., 300 and 600 K) using the AMBER force field. The in vacuum MD simulation results are in agreement with the crystallographic structure and with the DFT calculations of the nucleoside, revealing the anti conformer as the more stable one. The simulation in water, however, shows that both conformers may exist at 300 K, the temperature of the in vivo and in vitro assays, revealing that both the syn and anti conformers should be considered in a MD simulation study of the inhibitor–enzyme complex. Simulations are also in agreement with the NOE experiment, which shows that the anti conformer is the preferential one in DMSO-d6 solution at 298 K.     

Conformational space of the Pseudosaccharin allyl ether 3-(allyloxy)-1,2-benzisothiazole 1,1-dioxide in gas phase and in rare gas matrices

In this work, the conformational space of the pseudosaccharyl ether 3-(allyloxy)-1,2-benzisothiazole 1,1-dioxide (ABID) has been studied by means of infrared spectroscopy and density functional theory (DFT) calculations. Five different low energy conformers (TSk, TC, GSk, GSk' and GC, with relative energies of 0.00, 1.97, 2.00, 3.82 and 6.02 kJ mol(-1), respectively) were found on the DFT(B3LYP)/6-311++G(3df,3pd) potential energy surface of the molecule, all of them differing in the conformation of the allyl substituent. According to the calculations, in the gaseous phase all conformers are significantly populated (TSk:TC:GSk:GSk':GC = 47%:16%:18%:12%:7%, at 350 K). In the cryogenic matrices, however, only the TSk and TC conformers exist due to isomerization from the higher energy gauche forms to the most stable trans isomers during deposition of the matrix (conformational cooling). The observed conformational cooling is in consonance with the low calculated energy barriers for the GSk --> TSk, GSk' --> TSk and GC --> TC isomerization processes. Results from annealing experiments in krypton matrix doubtlessly show that in this matrix the order of stability of the TSk and TC conformers is reversed, with the more planar TC form becoming the most stable conformer.     

Conformational stability of cyclobutanol from temperature dependent infrared spectra of xenon solutions, r0 structural parameters, ab initio calculations and vibrational assignment

Variable temperature (-55 to -100 degrees C) studies of the infrared spectra (4000-400 cm(-1)) of cyclobutanol, c-C4H7OH dissolved in liquid xenon have been carried out. The infrared spectrum (4000-100 cm(-1)) of the gas has also been recorded. From these data two of the four possible stable conformers have been confidently identified and their order of stabilities has been experimentally determined where the first indicator is for the position of attachment of the hydroxyl group on the bent cyclobutyl ring (Eq=equatorial or Ax=axial) and the second one (t=trans, g=gauche) is the relative position of the hydroxyl rotor, i.e. rotation around the ring C-O bond. The enthalpy difference between the most stable Eq-t conformer and the second most stable rotamer, Eq-g, has been determined to be 200+/-50 cm(-1) (2.39+/-0.60 kJ/mol). This experimentally determined order is consistent with the order of stability predicted by ab initio calculations Eq-t>Eq-g>Ax-g>Ax-t. Evidence was obtained for the third conformer Ax-g which is predicted by ab initio calculations to be less stable by more than 650cm(-1) than the Eq-t form. The percentage of each conformer at ambient temperature is estimated to be Eq-t (50%), Eq-g (47%) and Ax-g (3%). The conformational stabilities, harmonic force fields, infrared intensities, Raman activities, depolarization ratios and vibrational frequencies have been obtained for all of the conformers from MP2(full)/6-31G(d) ab initio calculations. The optimized geometries and conformational stabilities have been obtained from ab initio calculations utilizing several different basis sets up to MP2(full)/aug-cc-pVTZ and from density functional theory calculations by the B3LYP method. By utilizing previously reported microwave rotational constants for the Eq-t conformer combined with ab initio MP2(full)/6-311+G(d,p) predicted structural values, adjusted r0 parameters have been obtained. The determined heavy atom structural parameters for the Eq-t conformer are: the distances C1-C4=1.547(5) angstroms, C4-C6=1.552(5)angstroms, C-O=1.416(5) angstroms and angles angleC6C4C1=86.6(5) degrees , angleC4C1C5=88.9(5) degrees and angleC6C5C1C4=22.8(5) degrees . The results are discussed and compared to the corresponding properties of some similar molecules.     

A sobering assessment of small‐molecule force field methods for low energy conformer predictions

We have carried out a large scale computational investigation to assess the utility of common small‐molecule force fields for computational screening of low energy conformers of typical organic molecules. Using statistical analyses on the energies and relative rankings of up to 250 diverse conformers of 700 different molecular structures, we find that energies from widely used classical force fields (MMFF94, UFF, and GAFF) show unconditionally poor energy and rank correlation with semiempirical (PM7) and Kohn–Sham density functional theory (DFT) energies calculated at PM7 and DFT optimized geometries. In contrast, semiempirical PM7 calculations show significantly better correlation with DFT calculations and generally better geometries. With these results, we make recommendations to more reliably carry out conformer screening.     

Conformational behavior of cis-2-methoxy, cis-2-methylthio, and cis-2-methylselenocyclohexanol: a theoretical and experimental investigation

Studies on the conformational equilibria of 2-methoxy, 2-methylthio, and 2-methylselenocyclohexanol are reported. Dynamic NMR spectroscopy experiments at 203-210 K were performed, which provided the percentages of each conformer in equilibrium. Theoretical calculations using the B3LYP method and aug-cc-pvdz basis set were applied to determine the differences in energy between the conformers. The analysis of the potential energy surface of each conformer showed the presence of two rotamers. Natural bond orbital analysis provided an explanation of which factors are driving the rotamer and conformer preferences.     

Investigation of the two-step spin crossover complex Fe[5-NO2-sal-(1,4,7,10)] using density functional theory

A quantum chemical study of the Fe[5-NO2-sal-(1,4,7,10)] ((1,10-bis(5-nitrosalicylaldehyde)-1,4,7,10-tetra-ezdecane-O,O',N,N',N' ',N' ')iron(II)) molecule was performed using density functional theory (DFT). Starting from the different X-ray crystallographic structures, geometry optimizations have been performed. These calculations confirmed the conformational isomerism of this complex in each spin states of the molecule ((1)A(1g) and (5)T(2g)). Each employed DFT method (B3LYP, B3LYP*, BP86, HCTH407) reproduced correctly the structural differences between the two calculated conformers when compared to the experimental structures. Furthermore, electronic polarizabilities have been calculated in each spin state and for each conformer. These calculations revealed a higher polarizability in the singlet state in agreement with the measured higher dielectric constant in this state.     

Conformation and stereodynamics of 1,2-diaryltetrahydropyrimidine and of its five- and seven-membered ring analogs

The 1-(2-nitrophenyl)-2-(2-methylphenyl)-1,4,5,6-tetrahydropyrimidine and its five- and seven-membered ring analogs were synthesized and their conformational properties investigated by low temperature NMR spectroscopy and DFT theoretical calculations. Restricted rotation of the aryl substituents were observed in all cases and the corresponding barriers determined. In the case of the six-membered ring derivative the additional conformers resulting from a ring inversion process were also detected.     

Evidence of conformational equilibrium of 1-ethyl-3-methylimidazolium in its ionic liquid salts: Raman spectroscopic study and quantum chemical calculations

Raman spectra of liquid 1-ethyl-3-methylimidazolium (EMI+) salts, EMI(+)BF4-, EMI(+)PF6-, EMI(+)CF3SO3-, and EMI(+)N(CF3SO2)2-, were measured over the frequency range 200-1600 cm(-1). In the range 200-500 cm(-1), we found five bands originating from the EMI+ ion at 241, 297, 387, 430, and 448 cm(-1). However, the 448 cm(-1) band could hardly be reproduced by theoretical calculations in terms of a given EMI+ conformer, implying that the band originates from another conformer. This is expected because the EMI+ involves an ethyl group bound to the N atom of the imidazolium ring, and the ethyl group can rotate along the C-N bond to yield conformers. The torsion energy for the rotation was then theoretically calculated. Two local minima with an energy difference of ca. 2 kJ mol(-1) were found, suggesting that two conformers are present in equilibrium. Full geometry optimizations followed by normal frequency analyses indicate that the two conformers are those with planar and nonplanar ethyl groups against the imidazolium ring plane, and the nonplanar conformer is favorable. It elucidates that bands at 241, 297, 387, and 430 cm(-1) mainly originate from the nonplanar conformer, whereas the 448 cm(-1) band does originate from the planar conformer. Indeed, the enthalpy for conformational change from nonplanar to planar EMI+ experimentally obtained by analyzing band intensities of the conformers at varying temperatures is practically the same as that evaluated by theoretical calculations. We thus conclude that the EMI+ ion exists as either a nonplanar or planar conformer in equilibrium in its liquid salts.     

Diastereomer configurations from joint experimental-computational analysis

The potential of the approach combining nuclear magnetic resonance (NMR) spectroscopy, relaxed grid search (RGS), molecular dynamics (MD) simulations, and quantum mechanical (QM) calculations for the determination of diastereomer configurations is demonstrated using four diastereomers of a trisubstituted epoxide. Since the change in configuration of the chiral center is expected to change the distribution of conformer populations (including those of side-chain rotamers), changes in NMR parameters [chemical shifts, J couplings, and nuclear Overhauser effects (NOEs)] are expected. The method therefore relies on (1) identification of possible conformations in each diastereomer using relaxed grid search analysis and MD simulations; (2) geometry optimizations of conformers selected from step (1), followed by calculations of their relative energies (populations) using QM methods; (3) calculations of averaged NMR parameters using QM methods; (4) matching calculated and experimental values of NMR parameters of diastereomers. The diastereomer configurations are considered resolved, if three NMR parameters different in nature, chemical shifts, J couplings, and NOEs, are in agreement. A further advantage of this method is that full structural and dynamics characterization of each of the diastereomers is achieved based on the joint analysis of experimental and computational data.     

A combined experimental and theoretical study on the conformational behavior of a calix[6]arene

An experimental and theoretical study on the conformational behavior of the 1,3,5-OMe-2,4,6-OCH(2)CONHOH-p-tert-butylcalix[6]arene has been carried out. In particular, semiempirical (AM1) and density functional theory (DFT) calculations have been performed in order to identify the possible conformers. The obtained results show that the cone structure is the most stable conformer at any level of theory, even if significant differences have been obtained for the other species. The inclusion of solvent effect, through a continuum model, also points out the relevant role played by the solvent in the stabilization of the cone structure in solution. These latter results have been confirmed by NMR experiments, which clearly show the presence of only the cone conformer in a polar solvent, such as DMSO. Finally, (1)H and (13)C NMR spectra on model systems, i.e., two successive phenol rings (Ar(1)-CH(2)-Ar(2)), have been computed at the DFT level and compared with the experimental spectra of the complete molecule. The results show an overall good agreement with the experimental data, thus leading to an unambiguous assignment of the experimental spectra.     

Impact of azaproline on amide cis-trans isomerism: conformational analyses and NMR studies of model peptides including TRH analogues

The beta-turn is a well-studied motif in both proteins and peptides. Four residues, making almost a complete 180 degree-turn in the direction of the peptide chain, define the beta-turn. Several types of the beta-turn are defined according to Phi and Psi torsional angles of the backbone for residues i + 1 and i + 2. One special type of beta-turn, the type VI-turn, usually contains a proline with a cis-amide bond at residue i + 2. In an aza-amino acid, the alpha-carbon of the amino acid is changed to nitrogen. Peptides containing azaproline (azPro) have been shown to prefer the type VI beta-turn both in crystals and in organic solvents by NMR studies. MC/MD simulations using the GB/SA solvation model for water explored the conformational preferences of azPro-containing peptides in aqueous systems. An increase in the conformational preference for the cis-amide conformer of azPro was clearly seen, but the increased stability was relatively minor with respect to the trans-conformer as compared to previous suggestions. To test the validity of the calculations in view of the experimental data from crystal structures and NMR in organic solvents, [azPro(3)]-TRH and [Phe(2), azPro(3)]-TRH were synthesized, and their conformational preferences were determined by NMR in polar solvents as well as the impact of the azPro substitution on their biological activities.     

Anomalous conformational change in 1-butyl-3-methylimidazolium tetrafluoroborate-D2O mixtures

We have investigated the effect of deuterated water on the conformational equilibrium between the gauche and trans conformers of the [bmim] cation in mixtures of water and 1-butyl-3-methylimidazolium tetrafluoroborate ([bmim][BF(4)]), an ionic liquid, at room temperature. A comparison of the results obtained from solutions made with H(2)O and with D(2)O highlights an anomalous conformational change in the D(2)O solution showing an extended N-shaped behavior. The gauche conformer of the [bmim] cation in D(2)O increased up to x = ~50 (D(2)O mol %); however, it decreased up to higher water concentrations of x = ~85 before again increasing drastically toward x = ~100. We provide spectroscopic evidence that the anomalous conformational dynamics of the [bmim] cation in D(2)O is directly related to the H/D exchange reaction of the C-H group at position 2 of the imidazolium ring.     

Circular dichroism,proton magnetic resonance and conformation of steroid 4-en-3-ones, 4, 9-dien-3-ones and 4, 9, 11-trien-3-ones

Room and low-temperature CD of the title dienones and trienones (3 to 9) clearly show that a conformational equilibrium occurs between two ring A half-chair conformers. The relative stability of the two conformers depends on substitution of ring A. These results may be extended to other series of steroids and in particular to 4-en-3-ones. CD of these enones (14 to 28) is consistent with a conformational equilibrium between the known quasi-cis-quasi-trans conformers. The controversial conformational behaviour of 2 β-substituted 19-nor-4-en-3-ones is explained by a dynamic equilibrium in solution between the two afore-mentioned conformers rather than by single twist or deformed boat conformer.     

Matrix isolation and low temperature solid state FTIR spectroscopic study of alpha-furil

Alpha-furil [C(4)H(3)O-C(=O)-C(=O)-C(4)H(3)O] has been isolated in argon and xenon matrices and studied by FTIR spectroscopy, supported by DFT(B3LYP)/6-311++G(d,p) calculations. The obtained spectra were fully assigned and revealed the presence in the matrices of three different conformers, all of them exhibiting skewed conformations around the intercarbonyl bond with the two C(4)H(3)O-C(=O) fragments nearly planar. The three conformers differ in the orientation of the furan rings relative to the carbonyl groups: the most stable conformer, I (C(2) symmetry; O=C-C=O intercarbonyl dihedral equal to 153.1 degrees), has both furan rings orientated in such a way that one of their beta-hydrogen atoms approaches the oxygen atom of the most distant carbonyl group, forming two H-C=C-C-C=O six-membered rings; the second most stable conformer, II (C(1) symmetry; O=C-C=O intercarbonyl dihedral equal to 126.9 degrees ), has one furan ring orientated as in I, while the second furan group is rotated by ca. 180 degrees (resulting in an energetically less favourable H-C=C-C=O five-membered ring); in the third conformer, III (C(2) symmetry; O=C-C=O dihedral equal to 106.2 degrees ), both furan rings assume the latter orientation relative to the dicarbonyl group. The theoretical calculations predicted the two higher energy forms being 5.85 and 6.22 kJ mol(-1) higher in energy than the most stable form, respectively, and energy barriers for conformational interconversion higher than 40 kJ mol(-1). These barriers are high enough to prevent observation of conformational isomerization for the matrix isolated compound. The three possible conformers of alpha-furil were also found to be present in CCl(4) solution, as well as in a low temperature neat amorphous phase of the compound prepared from fast condensation of its vapour onto a suitable 10 K cooled substrate. On the other hand, in agreement with the available X-ray data [S. C. Biswas, S. Ray and A. Podder, Chem. Phys. Lett., 1987, 134, 541], the IR spectra obtained for the neat low temperature crystalline state reveals that, in this phase, alpha-furil exists uniquely in its most stable conformational state, I.     

Structure and conformational equilibrium of thiacalix[4]arene by density functional theory

Density functional theory calculations for the structure and conformational equilibrium of thiacalix[4]arene are reported. The conformational equilibrium of thiacalix[4]arene, a heterocalixarene in which the phenol groups are bridged by sulphur atoms is compared to the conformational equilibrium of calix[4]arene. Thiacalix[4]arene conformational energies relative to the cone conformer (ΔE's) are reduced in comparison with calix[4]arene. This conformational change is in qualitative agreement with recent NMR spectroscopy measurements of the conformational equilibrium for a tetraethylether of thiacalix[4]arene in a CDCl₃ solution which indicates an enhanced chemical exchange of thiacalixarene conformers in comparison with similar methylene bridged structures. Density functional theory results for the structure of thiacalix[4]arene are in good agreement with recent X-ray diffraction measurements. The electrostatic potentials in the cone conformers of thiacalix[4]arene and calix[4]arene suggest that their complexation or recognition abilities can be significantly different. Dipole moments of the four thiacalix[4]arene conformers are in the order: cone>1,2-alternate>partial-cone>1,3-alternate.