Kinetics and equilibria of cis/trans isomerization of backbone amide bonds in peptoids

The biological activities of N-substituted glycine oligomers (peptoids) have been the subject of extensive research. As compared to peptides, both the cis and trans conformations of the backbone amide bonds of peptoids can be significantly populated. Thus, peptoids are mixtures of configurational isomers, with the number of isomers increasing by a factor of 2 with each additional monomer residue. Here we report the results of a study of the kinetics and equilibria of cis/trans isomerization of the amide bonds of N-acetylated peptoid monomers, dipeptoids, and tripeptoids by NMR spectroscopy. Resonance intensities indicate the cis conformation of the backbone amide bonds of the peptoids studied is more populated than is generally the case for the analogous secondary amide bond to proline residues in acyclic peptides. Rate constants were measured by inversion-magnetization transfer techniques over a range of temperatures, and activation parameters were derived from the temperature dependence of the rate constants. The rate of cis/trans isomerization by rotation around the amide bonds in the peptoids studied is generally slower than that around amide bonds to proline residues and takes place on the NMR inversion-magnetization transfer time scale only by rotation around the amide bond to the C-terminal peptoid residue.     

Puckering transition of 4-substituted proline residues

The puckering transition of 4-substituted proline residues by electron-withdrawing groups, i.e., 4(R)-hydroxy-L-proline (Hyp) and 4(R)-fluoro-L-proline (Flp) residues, with trans and cis prolyl peptide bonds was studied by adiabatic optimizations along the torsion angle chi1 of the prolyl ring at the HF/6-31+G(d) level. By analyzing the potential energy surface and local minima, it is observed that the puckering transition of the prolyl ring for Hyp and Flp residues proceeds from a down-puckered conformation to an up-puckered one through the transition state with an envelope form having the N atom at the top of envelope and not a planar one for both trans and cis conformers, which is the same as found for the unsubstituted proline residue. At HF/6-31+G(d) and B3LYP/6-311++G(d,p) levels, the structures of the backbone and prolyl ring for local minima of Ac-Hyp-NHMe and Ac-Flp-NHMe are quite similar to those of Ac-Pro-NHMe. However, the relative stability of the up-puckered conformation to the down-puckered one is increased for Ac-Hyp-NHMe with the cis imide bond and for Ac-Flp-NHMe with the trans and cis imide bonds. In particular, the 4(R)-substitution by hydroxy and fluorine groups has brought some structural changes in the prolyl ring of the transition states and the changes in barriers for the puckering transition. The puckering transitions for Ac-Hyp-NHMe and Ac-Flp-NHMe are proven to be predominantly electronically driven by analyzing the electronic and enthalpic contributions to the barriers, as seen for Ac-Pro-NHMe.     

Conformational preferences of non-prolyl and prolyl residues

The conformational study on Ac-Ala-NHMe (the alanine dipeptide) and Ac-Pro-NHMe (the proline dipeptide) is carried out using ab initio HF and density functional methods with the self-consistent reaction field method to explore the differences in the backbone conformational preference and the cis-trans isomerization for the non-prolyl and prolyl residues in the gas phase and in the solutions (chloroform and water). For the alanine and proline dipeptides, with the increase of solvent polarity, the populations of the conformation tC with an intramolecular C(7) hydrogen bond significantly decrease, and those of the polyproline II-like conformation tF and the alpha-helical conformation tA increase, which is in good agreement with the results from circular dichroism and NMR experiments. For both the dipeptides, as the solvent polarity increases, the relative free energy of the cis conformer to the trans conformer decreases and the rotational barrier to the cis-trans isomerization increases. It is found that the cis-trans isomerization proceeds in common through only the clockwise rotation with omega' approximately +120 degrees about the non-prolyl and prolyl peptide bonds in both the gas phase and the solutions. The pertinent distance d(N...H-N(NHMe)) can successfully describe the increase in the rotational barriers for the non-prolyl and prolyl trans-cis isomerization as the solvent polarity increases and the higher barriers for the non-prolyl residue than for the prolyl residue, as seen in experimental and calculated results. By analysis of the contributions to rotational barriers, the cis-trans isomerization for the non-prolyl and prolyl peptide bonds is proven to be entirely enthalpy driven in the gas phase and in the solutions. The calculated cis populations and rotational barriers to the cis-trans isomerization for both the dipeptides in chloroform and/or water accord with the experimental values.     

End-to-end vs interior loop formation kinetics in unfolded polypeptide chains

The conformational search for favorable intramolecular interactions during protein folding is limited by intrachain diffusion processes. Recent studies on the dynamics of loop formation in unfolded polypeptide chains have focused on loops involving residues near the chain ends. During protein folding, however, most contacts are formed between residues in the interior of the chain. We compared the kinetics of end-to-end loop formation (type I loops) to the formation of end-to-interior (type II loops) and interior-to-interior loops (type III loops) using triplet-triplet energy transfer from xanthone to naphthylalanine. The results show that formation of type II and type III loops is slower compared to type I loops of the same size and amino acid sequence. The rate constant for type II loop formation decreases with increasing overall chain dimensions up to a limiting value, at which loop formation is about 2.5-fold slower for type II loops compared to type I loops. Comparing type II loops of different loop size and amino acid sequence shows that the ratio of loop dimension over total chain dimension determines the rate constant for loop formation. Formation of type III loops is 1.7-fold slower than formation of type II loops, indicating that local chain motions are strongly coupled to motions of other chain segments which leads to faster dynamics toward the chain ends. Our results show that differences in the kinetics of formation of type I, type II, and type III loops are mainly caused by differences in internal flexibility at the different positions in the polypeptide chain. Interactions of the polypeptide chain with the solvent contribute to the kinetics of loop formation, which are strongly viscosity-dependent. However, the observed differences in the kinetics of formation of type I, type II, and type III loops are not due to the increased number of peptide-solvent interactions in type II and type III loops compared to type I loops as indicated by identical viscosity dependencies for the kinetics of formation of the different types of loops.     

cis,cis- and trans,trans-ceratospongamide, new bioactive cyclic heptapeptides from the Indonesian red alga Ceratodictyon spongiosum and symbiotic sponge Sigmadocia symbiotica

Chemical investigation of the marine red alga (Rhodophyta) Ceratodictyon spongiosum containing the symbiotic sponge Sigmadocia symbiotica collected from Biaro Island, Indonesia, yielded two isomers of a new and bioactive thiazole-containing cyclic heptapeptide, cis,cis-ceratospongamide (1) and trans, trans-ceratospongamide (2). Isolation of these peptides was assisted by bioassay-guided fractionation using a brine shrimp toxicity assay (Artemia salina). The structures of the ceratospongamides, which each consist of two L-phenylalanine residues, one (L-isoleucine)-L-methyloxazoline residue, one L-proline residue, and one (L-proline)thiazole residue, were established through extensive NMR spectroscopy, including (1)H-(13)C HMQC-TOCSY, and (1)H-(15)N HMBC experiments, as well as chemical degradation and chiral analysis. cis,cis- and trans,trans-ceratospongamide are stable conformational isomers of the two proline amide bonds. Molecular modeling of these two ceratospongamide isomers showed the trans, trans isomer to be quite planar, whereas the cis,cis isomer has a more puckered overall conformation. trans,trans-Ceratospongamide exhibits potent inhibition of sPLA(2) expression in a cell-based model for antiinflammation (ED(50) 32 nM), whereas the cis,cis isomer is inactive. trans,trans-Ceratospongamide was also shown to inhibit the expression of a human-sPLA(2) promoter-based reporter by 90%.     

Puckering transition of proline residue in water

The puckering transition of the proline residue with trans and cis prolyl peptide bonds was explored by optimizations along the torsion angle chi1 of the prolyl ring using quantum-chemical methods in water. By analyzing the potential energy surfaces and local minima in water, it is observed that the puckering transition of the proline residue proceeds from a down-puckered conformation to an up-puckered one and vice versa through the transition state with an envelope form having the N atom at the top of the envelope and not a planar one, as seen in the gas phase, although the backbone conformations are different in the gas phase and in water. The barriers to the puckering transition DeltaGup-->down are estimated to be 3.12 and 3.00 kcal/mol for trans and cis conformers at the B3LYP/6-311++G(d,p) level of theory in water, respectively, which are about 1.7 kcal/mol higher than those in the gas phase. Out of 2197 prolines from the 241 high-resolution PDB chains, four transition-state-like structures with the envelope ring puckering are identified. Three of them have the trans prolyl peptide bonds and one has the cis one. The favorable or steric interactions by neighboring residues may be responsible for the stabilization of these transition-state-like ring structures in the proteins.     

Conformational preferences of proline oligopeptides

A conformational study on the terminally blocked proline oligopeptides, Ac-(Pro)(n)()-NMe(2) (n = 2-5), is carried out using the ab initio Hartree-Fock level of theory with the self-consistent reaction field method in the gas phase and in solutions (chloroform, 1-propanol, and water) to explore the preference and transition between polyproline II (PPII) and polyproline I (PPI) conformations depending on the chain length, the puckering, and the solvent. The mean differences in the free energy per proline of the up-puckered conformations relative to the down-puckered conformations for both diproline and triproline increases for the PPII-like conformations and decreases for the PPI-like conformations as the solvent polarity increases. These calculated results indicate that the PPII-like structures have preferentially all-down puckerings in solutions, whereas the PPI-like structures have partially mixed puckerings. The free energy difference per proline residue between the PPII- and PPI-like structures decreases as the proline chain becomes longer in the gas phase but increases as the proline chain becomes longer in solutions and the solvent polarity increases. In particular, our calculated results indicate that each of the proline oligopeptides can exist as an ensemble of conformations with the trans and cis peptide bonds in solutions, although the PPII-like structure with all-trans peptide bonds is dominantly preferred, which is reasonably consistent with the previously observed results. In diproline Ac-(Pro)(2)-NMe(2), the rotational barrier to the cis-to-trans isomerization for the first prolyl peptide bond increases as the solvent polarity increases, whereas the rotational barrier for the second prolyl peptide bond does not show the monotonic increase as the solvent polarity increases. When the rotational barriers for these two prolyl peptide bonds were compared, it could be deduced that the conformational transition from PPI with the cis peptide bond to PPII with the trans peptide bond is initiated at the C-terminus and proceeds to the N-terminus in water. This is consistent with the results from NMR experiments on polyproline in D(2)O but opposite to the results from enzymatic hydrolysis kinetics experiments on polyproline.     

Conformational profile of a proline-arginine hybrid

The intrinsic conformational preferences of a new nonproteinogenic amino acid have been explored by computational methods. This tailored molecule, named ((β)Pro)Arg, is conceived as a replacement for arginine in bioactive peptides when the stabilization of folded turn-like conformations is required. The new residue features a proline skeleton that bears the guanidilated side chain of arginine at the C(β) position of the five-membered pyrrolidine ring, in either a cis or a trans orientation with respect to the carboxylic acid. The conformational profiles of the N-acetyl-N'-methylamide derivatives of the cis and trans isomers of ((β)Pro)Arg have been examined in the gas phase and in solution by B3LYP/6-31+G(d,p) calculations and molecular dynamics simulations. The main conformational features of both isomers represent a balance between geometric restrictions imposed by the five-membered pyrrolidine ring and the ability of the guanidilated side chain to interact with the backbone through hydrogen bonds. Thus, both cis- and trans-((β)Pro)Arg exhibit a preference for the α(L) conformation as a consequence of the interactions established between the guanidinium moiety and the main-chain amide groups.     

Electronic control of amide cis-trans isomerism via the aromatic-prolyl interaction

The cis-trans isomerization of prolyl amide bonds results in large structural and functional changes in proteins and is a rate-determining step in protein folding. We describe a novel electronic strategy to control cis-trans isomerization, based on the demonstration that interactions between aromatic residues and proline are tunable by aromatic electronics. A series of peptides of sequence TXPN, X = Trp, pyridylalanine, pentafluorophenylalanine, or 4-Z-phenylalanine derivatives (Z = electron-donating, electron-withdrawing, or electron-neutral substituents), was synthesized and Ktrans/cis analyzed by NMR. Electron-rich aromatic residues stabilized cis amide bond formation, while electron-poor aromatics relatively favored trans amide bond formation. A Hammett correlation between aromatic electronics and cis-trans isomerization was observed. These results indicate that the interaction between aromatic residues and proline, which is observed to stabilize cis amide bonds and is also a general stabilizing interaction ubiquitous in proteins and protein-protein complexes, is not stabilized exclusively by a classical hydrophobic effect. To a large extent, the aromatic-prolyl interaction is driven and controllable by an electronic effect between the aromatic ring pi-electrons and the proline ring, consistent with a C-H-pi interaction as the key stabilizing force. The aromatic-prolyl interaction is electronically tunable by 0.9 kcal/mol and is enthalpic in nature. In addition, by combining aromatic ring electronics and stereoelectronic effects using 4-fluoroprolines, we demonstrate broad tuning (2.0 kcal/mol) of cis-trans isomerism in tetrapeptides. We demonstrate a simple tetrapeptide, TWflpN, that exhibits 60% cis amide bond and adopts a type VIa1 beta-turn conformation.     

Conformational preferences of beta- and gamma-aminated proline analogues

Quantum mechanical calculations have been used to investigate how the incorporation of an amino group to the Cbeta- or Cgamma-positions of the pyrrolidine ring affects the intrinsic conformational properties of the proline. Specifically, a conformational study of the N-acetyl-N'-methylamide derivatives of four isomers of aminoproline, which differ not only in the beta- or gamma-position of the substituent but also in its cis or trans relative disposition, has been performed. To further understand the role of the intramolecular hydrogen bonds between the backbone carbonyl groups and the amino side group, a conformational study was also performed on the corresponding four analogues of (dimethylamino)proline. In addition, the effects of solvation on aminoproline and (dimethylamino)proline dipeptides have been evaluated using a self-consistent reaction field model, and considering four different solvents (carbon tetrachloride, chloroform, methanol and water). Results indicate that the incorporation of the amino substituent into the pyrrolidine ring affects the conformational properties, with backbone...side chain intramolecular hydrogen bonds detected when it is incorporated in a cis relative disposition. In general, the incorporation of the amino side group tends to stabilize those structures where the peptide bond involving the pyrrolidine nitrogen is arranged in cis. The aminoproline isomer with the substituent attached to the Cgamma-position with a cis relative disposition is the most stable in the gas phase and in chloroform, methanol and water solutions. Replacement of the amino side group by the dimethylamino substituent produces significant changes in the potential energy surfaces of the four investigated (dimethylamino)proline-containing dipeptides. Thus, these changes affect not only the number of minima, which increases considerably, but also the backbone and pseudorotational preferences. In spite of these effects, comparison of the conformational preferences, i.e., the more favored conformers, calculated for different isomers of aminoproline and (dimethylamino)proline dipeptides showed a high degree of consistency for the two families of compounds.     

Substituted 2-azabicyclo[2.1.1]hexanes as constrained proline analogues: implications for collagen stability

Among the proteinogenic amino acids, only proline is a secondary amine and only proline has a saturated ring. Electronegative substituents on C-4 (that is, C(gamma)) have a substantial effect on the trans/cis ratio of the prolyl peptide bond and the pucker of the pyrrolidine ring. 2-Azabicyclo[2.1.1]hexane is, in essence, a proline analogue with two C(gamma) atoms, one in each of the two prevalent ring puckers of proline. Here, 2-azabicyclo[2.1.1]hexane analogues of 2S-proline, (2S,4S)-4-hydroxyproline, and (2S,4S)-4-fluoroproline residues were synthesized, and their trans/cis ratios were shown to be invariant in a particular solvent. Thus, the substitution of a proline residue on C-4 affects the trans/cis ratio by altering the pucker of its pyrrolidine ring. This finding has implications for the conformation of collagen, which has an abundance of 2S-proline and (2S,4R)-4-hydroxyproline residues, and can be stabilized by (2S,4R)-4-fluoroproline and (2S,4S)-4-fluoroproline residues.     

Conformational preferences of pseudoproline residues

The conformational study on N-acetyl-N'-methylamides of oxazolidine and thiazolidine residues (Ac-Oxa-NHMe and Ac-Thz-NHMe) is carried out using ab initio HF and density functional B3LYP methods with the self-consistent reaction field method to explore the effects of the replacement of the C(gamma)H(2) group in the prolyl ring by oxygen or sulfur atoms on the conformational preferences and prolyl cis-trans isomerization in the gas phase and in solution (chloroform and water). As the solvent polarity increases, the conformations C with the C7 intramolecular hydrogen bonds become depopulated, the PPII- or PPI-like conformations F become more populated, and the cis populations increase for both Oxa and Thz dipeptides, as found for the Pro dipeptide, although the populations of backbone conformations and puckerings are different in pseudoproline and proline dipeptides. As the increase of solvent polarity, the populations of the trans/up conformations decrease for Oxa and Thz dipeptides, but they increase for the Pro dipeptide. It is found that the cis-trans isomerization proceeds through the anticlockwise rotation with omega' approximately -60 degrees about the oxazolidyl peptide bond and the clockwise rotation with omega' approximately +120 degrees about the thiazolidyl peptide bond in the gas phase and in solution, whereas the clockwise rotation is preferred for the prolyl peptide bond. The pertinent distance d(N...H-N(NHMe)) and the pyramidality of the prolyl nitrogen can describe the role of this hydrogen bond in stabilizing the transition state structure but the lower rotational barriers for Oxa and Thz dipeptides than those for the Pro dipeptide, which is observed from experiments, cannot be rationalized. The calculated cis populations and rotational barriers to the cis-trans isomerization for both Oxa and Thz dipeptides in chloroform and/or water are consistent with the experimental values.     

CIS-Trans isomerism of a,2-bistrichlorosilylethylene

Summary The existence of cis and trans isomers of 1,2-bis-rrichlorosilylethylene was detected spectrographically. It was noted that the frequencies of the double bonds of the cis and trans isomers of 1,2-bistrichlorosilyethylene, are unusually far apart as compared with the frequencies found for the cis and trans isomers of other compounds.     

First-principle computational study on the full conformational space of L-threonine diamide, the energetic stability of cis and trans isomers

First-principle computations were carried out on the conformational space of trans and cis peptide bond isomers of HCO-Thr-NH2. Using the concept of multidimensional conformational analysis (MDCA), geometry optimizations were performed at the B3LYP/6-31G(d) level of theory, and single-point energies as well as thermodynamic functions were calculated at the G3MP2B3 level of theory for the corresponding optimized structures. Two backbone Ramachandran-type potential energy surfaces (PESs) were computed, one each for the cis and trans isomers, keeping the side chain at the fully extended orientation (chi1=chi2=anti). Similarly, two side chain PESs for the cis and trans isomers were generated for the (phi=psi=anti) orientation corresponding to approximately the betaL backbone conformation. Besides correlating the relative Gibbs free energy of the various stable conformations with the number of stabilizing hydrogen bonds, the process of trans-->cis isomerization is discussed in terms of intrinsic stabilities as measured by the computed thermodynamic functions.     

Conformational analysis of the complexes of copper(ii) with two identical α-amino acid ligands,relative stability of the cis and trans isomers

Conclusions The energy closeness of the cis and trans isomers of the (Ala)₂Cu complexes, and consequently the possibility of the presence of the cis isomers in the equilibrium mixture, was shown. The differences in the energies of the cis and trans isomers increase in the complexes with the N-mono- and N,N-dimethyl-substituted alanine. The coexistence of both isomers is thereby possible for the monosubstituted compounds, but the disubstituted cis isomers are virtually unrealized. According to the conformational calculations, the formation of the optically active compounds containing ligands of one configurational series is more preferred among the complexes considered.Translated from Izvestiya Akademii Nauk SSSR, Seriya Khimicheskaya, No. 11, pp. 2484–2487, November, 1987.     

Mechanism of photoisomerization of 1-naphthyl-2-phenylethylenes in organic glasses

Photoisomerization reactions of cis isomers of 1-beta-naphthyl-2-phenylethylene, an o-methylated homolog and 1-alpha-naphthyl-2-phenylethylene in organic glasses at liquid nitrogen temperature were studied. Reactions were followed by changes in UV-absorption spectra of the irradiated samples. Formation of an unstable trans-photoproduct was detected only with the o-methylated homolog. The results are consistent with high regioselective Hula-twist photoisomerization at the benzylic sites of the three compounds examined. Calculated data on relative energies of the conformers of both the cis and the trans isomers are in agreement with the suspected conformational population of the starting materials and the photoproducts.     

The mechanism of cis-trans isomerization of prolyl peptides by cyclophilin

The mechanism of cis-trans isomerization of prolyl peptides catalyzed by cyclophilin (CyP) was studied computationally via molecular dynamics (MD) simulations of the transition state (TS) and the cis and trans forms of the ground state (GS), when bound to CyP and when free in aqueous solution. The MD simulations include four enzyme-bound species of tetrapeptide (Suc-Ala-XC([double bond]O)-NPro-Phe-pNA; X = Gly, Trp, Ala, and Leu). In water, the prolyl amide bond is favorably planar with the presence of conformers exhibiting +/-20 degrees twist of the C-N dihedral. In the active site a hydrogen bond between the cis-prolyl amide carbonyl O and the backbone amide N-H of Asn102 retains the 20 degrees twist of the C-N dihedral. The TS structure is characterized by a 90 degrees twist of the amide C-N bond and a more favorable interaction with Asn102 due to the shorter distance between Asn102(HN) and the amide carbonyl O. The conformational change of cis --> TS also involves pyramidalization of the amide N, which results in the formation of a hydrogen bond between the amide N and the guanidino group of Arg55. Both Asn102 and Arg55 are held in the same position in CyP.cis-isomer as in CyP.TS. In the ligand-free CyP the Arg55 guanidino group is highly disorganized and Asn102 is displaced 1 A from the position in the ligand-bound CyP. Thus, the organization of Arg55 and Asn102 occurs upon substrate binding. The geometrical complimentarity of the organized enzyme structure to the TS structure is a result of preferential binding of the proline N and the amide carbonyl of the TS compared to that of GS. However, the N-terminal part (Suc-Ala) becomes repositioned in the TS such that two hydrogen bonds disappear, one hydrogen bond appears and two other hydrogen bonds becomes weaker on the conversion of CyP.cis to CyP.TS. During this conversion, total hydrophobic contact between enzyme and the peptide is preserved. Thus, the interaction energies of GS and TS with enzyme are, as a whole, much alike. This does not support the contention that TS is bound more tightly than GS by K(m)/K(TS) = 10(6) in the cis --> trans reaction. Repositioning of the N-terminal part of the peptide on CyP.TS formation becomes more pronounced when the substrate X residue is changed from Gly < Trp < Ala < Leu. We propose that the larger turning of the N-terminus is responsible for the larger value of the experimentally observed Delta S(++) and Delta H(++), which sum up to little change in Delta G(++). The positioning of the Arg55 and the degree of 20 degrees twist of the amide C-N bond are considered as criteria for Near Attack Conformers (NACs) in cis-trans isomerization. NACs account for approximately 30% of the total GS populations of the cis-isomer. Similar NAC populations were observed with four different substrates. This is consistent with the insensitivity of enzymatic activity to the nature of the X residue. Also, the NAC population in CyP.trans-AAPF was comparable to that in CyP.cis-AAPF, in accord with similar experimentally measured rates of the cis --> trans and trans --> cis reaction in CyP. These NACs, found in CyP.cis and CyP.trans, resemble only one of the four possible TS configurations in the water reaction. The identity of this TS structure (syn/exo) is in accord with experimentally determined KIE values in the enzymatic reaction. However, the geometry of the active site was also complementary to another TS structure (anti/exo) that was not detected in the active site by the same KIE measurements, implying that the geometrical fitness of the TS cannot be a single determining factor for enzymatic reactions.     

Conformational analysis of dehydrodidemnin B (aplidine) by NMR spectroscopy and molecular mechanics/dynamics calculations.

Dehydrodidemnin B (DDB or aplidine), a potent antitumoral natural product currently in phase II clinical trials, exists as an approximately 1:1 mixture of two slowly interconverting conformations. These are sufficiently long-lived so as to allow their resolution by HPLC. NMR spectroscopy shows that this phenomenon is a consequence of restricted rotation about the Pyr-Pro(8) terminal amide bond of the molecule's side chain. The same technique also indicates that the overall three-dimensional structures of both the cis and trans isomers of DDB are similar despite the conformational change. Molecular dynamics simulations with different implicit and explicit solvent models show that the ensembles of three-dimensional structures produced are indeed similar for both the cis and trans isomers. These studies also show that hydrogen bonding patterns in both isomers are alike and that each one is stabilized by a hydrogen bond between the pyruvyl unit at the terminus of the molecule's side chain and the Thr(6) residue situated at the junction betwen the macrocycle and the molecule's side chain. Nevertheless, each conformational isomer forms this hydrogen bond using a different pyruvyl carbonyl group: CO(2) in the case of the cis isomer and CO(1) in the case of the trans isomer.     

Conformational preferences of alpha-substituted proline analogues

DFT calculations at the B3LYP/6-31+G(d,p) level have been used to investigate how the replacement of the alpha hydrogen by a more sterically demanding group affects the conformational preferences of proline. Specifically, the N-acetyl-N'-methylamide derivatives of L-proline, L-alpha-methylproline, and L-alpha-phenylproline have been calculated, with both the cis/trans isomerism of the peptide bonds and the puckering of the pyrrolidine ring being considered. The effects of solvation have been evaluated by using a Self-Consistent Reaction Field model. As expected, tetrasubstitution at the alpha carbon destabilizes the conformers with one or more peptide bonds arranged in cis. The lowest energy minimum has been found to be identical for the three compounds investigated, but important differences are observed regarding other energetically accessible backbone conformations. The results obtained provide evidence that the distinct steric requirements of the substituent at C (alpha) may play a significant role in modulating the conformational preferences of proline.