Open Access Article
Vladimir
Kubyshkin
* and
Nediljko
Budisa
*
Biocatalysis Group, Institute of Chemistry, Technical University of Berlin, Müller-Breslau-Str. 10, Berlin 10623, Germany. E-mail: kubyshkin@win.tu-berlin.de; nediljko.budisa@tu-berlin.de
First published on 24th July 2017
Amide rotation of peptidyl–prolyl fragments is an important factor in backbone structure organization of proteins. Computational studies have indicated that this rotation preferentially proceeds through a defined transition-state structure (syn/exo). Here, we complement the computational findings by determining the amide bond rotation barriers for derivatives of the two symmetric proline analogues, meso and racemic pyrrolidine-2,5-dicarboxylic acids. The rotations around these residues represent syn/exo–syn/exo and anti/endo–syn/exo hybrid transition states for the meso and racemic diastereomer, respectively. The rotation barriers are lower for the former rotation by about 9 kJ mol−1 (aqueous medium), suggesting a strong preference for the syn/exo (clockwise) rotation over the anti/endo (anticlockwise) rotation. The results show that both hybrid rotation processes are enthalpically driven but respond differently to solvent polarity changes due to the different transition state dipole–dipole interactions.
![]() | ||
| Fig. 1 (A) The structure of proline (1), (B) two N-acyl-prolyl ground states, s-trans and s-cis, (C) four possible transition states for the amide rotation, the syn/exo transition state is most favored. The nomenclature is from ref. 21. | ||
The two ground states of the N-acyl prolyl amide bond are s-trans and s-cis (Fig. 1B), with s-trans conformation forming preferentially. The preference is about 3
:
1 in Ac-Pro-NHMe,15 4
:
1 in Ac-Pro-OH,16 5
:
1 in Ac-Pro-OMe,17 and 7
:
1 in GlyProGly-NH218 and AcGlyGlyProGlyGlyNH2 peptides (water, 298 K),19 which translates to approximately 4–5 kJ mol−1 in energetic preference. At the same time, the amide bond rotation should proceed through one of the four theoretically possible transition states, in which the amide conjugation should be lost (Fig. 1C). Several structures have successfully depicted the non-conjugated state trapped by a rigid molecular scaffold (so-called twisted amides).20 However, these elegant structures do not definitively model the amide transition state due to the following reason. Already in the mid-1990s the results from a theoretical study showed that the amide rotation preferably proceeds via the syn/exo transition states with the antiparallel orientation of the N → lone pair and carbonyl C → O dipoles,21 whereas in the covalent twisted amides this orientation rather resembles a parallel one. Only very recently metal coordination based twisted amides have been reported, and these demonstrated the antiparallel twist.22
Since the seminal theoretical reports, numerous computational studies have demonstrated that the amide rotation around proline has a preference for the syn/exo rotation scenario.23 These findings are further supported by recent theoretical24 and experimental25 studies of the catalytic mechanism of Pin1 peptidyl–prolyl cis–trans isomerase, and experimental demonstration of the self-catalytic mechanism of the amide rotation in nonpolar media.26
In the syn/exo transition state, the nitrogen lone pair is oriented on the same side and the carbonyl oxygen rotates into the opposite side with respect to the proline carbonyl group, whereas the upstream alkyl moiety is oriented outward. Curiously, despite the existing consensus among theoretical studies, a proper experimental model that addresses the energetic differences in the relevant rotation modes is lacking. While there is little doubt that the syn/exo transition state is the most relevant, it is not clear whether the anti/endo transition state can contribute to the rotation or whether the properties of the syn/exo conformation can be extrapolated to other transition scenarios. Therefore we designed an experimental model based on derivatives of pyrrolidine-2,5-dicarboxylic acid (Pdc). Pdc has two diastereomeric forms, meso (mPdc, with a symmetry plane, σ) and racemic (rPdc, with a symmetry axis, C2). We assumed that the amide ground states are equivalent in the N-acyl Pdc derivatives when the acyl-moiety is not chiral. At the same time, we expected the amide rotation to proceed via the syn/exo–syn/exo hybrid state for mPdc and the anti/endo–syn/exo hybrid state for rPdc derivatives (Fig. 2) (using the proline transition state nomenclature, Fig. 1C). We therefore synthesized a series of amides with non-chiral N-terminal acyl groups (Fig. 3). Here, we report the experimental amide bond rotation barriers for these two model situations. The results obtained for proline and pyrrolidine derivatives will also be shown for reference but will not be discussed in the context of this paper.
![]() | ||
| Fig. 2 Design of the study. Two diastereomeric forms of pyrrolidine-2,5-dicarboxylic acid yield different rotation modes around the N-acyl derivatives. | ||
![]() | ||
| Fig. 4 The pKa of the ammonium group in the pyrrolidine-2,5-dicarboxylic acid derivatives. The data on pyrrolidine and proline methyl esters are included for comparison. Standard error: ±0.10. | ||
Next, the X-ray crystal structures for both N-acetyl derivatives 2a and 3a demonstrate fully planar amide bonds (Fig. 5A).27 In the mPdc derivative 2a, the amide conformation breaks the symmetry of the structure and creates two enantiomers, which were present in a ratio of 4
:
2 in the cell unit. In the rPdc derivative 3a, the amide does not break the existing C2 symmetry of the structure, and the cell unit contains a racemate with a 2
:
2 ratio of the (R,R) and (S,S) enantiomers. The envelope conformation of the pyrrolidine ring in proline is usually described in terms of exo-/endo-pucker nomenclature (alternatively called up-/down-) depending on the direction of the C4-atom displacement from the plane of the ring relative to the carboxyl group orientation. Following this nomenclature, the mPdc derivative 2a adopts a hybrid of exo- and endo-pucker, whereas the rPdc derivative 3a is a hybrid of two endo-puckers. The crystal structure is consistent with the solution 1H NMR spectra, in which the two CH(CO2Me) resonances appear as a triplet and doublet of doublets in 2a and as two doublets in 3a (Fig. 5B).28 Importantly, these NMR spectral signatures persisted in all examined Pdc derivatives and in all employed NMR solvents (see ‘Polar effects’).
Based on this spatial organization argument we expect the steric requirements to be identical around the methyl group of the acetyl moiety, whereas the environment of the carbonyl group is different. These differences are very well described in the proline literature in terms of the pre-organization of the n → π* interaction between the carbonyl groups.29 As a result, some additional stability of about 1–3 kJ mol−1 can be expected for the mPdc derivatives relative to the rPdc-derived compounds irrespective from the polarity of the medium.17
Another aspect is the ψ-torsion due to the rotation of the carboxyl-group. For the proline derivative this rotation yields two states ψ1 ∼170° and ψ2 ∼−20° with different polarities. In 5a the first state is preferred by about 1.3–2.5 kJ mol−1.30 We modelled the carboxyl-group rotation in 2a–3a, and found that this preference remains nearly identical for the rPdc derivative 3a (ΔΔGψ1/ψ2∼ 1.6–2.7 kJ mol−1), whereas in 2a the two carboxymethyl groups are more proximal, and the rotation state difference becomes more prominent (ΔΔGψ1/ψ2∼ 3.7–4.7 kJ mol−1). This, however, produces only a little energy difference between the two ground states, and the resulting polarity differences are expected to become influential only in nonpolar media.30
Hereafter, we will not correct for these energies and will consider both ground states equivalent for simplicity of presentation.
![]() | ||
| Fig. 6 Amide rotation barriers for the N-acetyl derivatives 2a, 3a and 4a as a function of temperature, measured by 1H EXSY NMR in deuterium oxide. For details, see Tables S1 and S3.‡ | ||
For all three examined structures, the amide rotation is an enthalpically defined process with a small negative slope, −TΔS≠ (ΔS≠ > 0). The smallest entropic contribution was found for pyrrolidine derivative 4a (ΔS≠ = +3.2 ± 1.8 J mol−1 K−1), whereas for the mPdc derivative 2a, this value was somehow larger (ΔS≠ = +11.9 ± 0.9 J mol−1 K−1). We were not able to determine the entropic value for the rPdc derivative 3a with the same level of accuracy (ΔS≠ = +7.6 ± 9.0 J mol−1 K−1) due to the high relative experimental error for the very slow rates under the lower temperature conditions. Nonetheless, this lack of accuracy does not affect the overall conclusion, as the same tendencies were found for the N-acetylglycyl derivatives 2c–4c. In the latter case, the activation enthalpy was reduced by 5–7 kJ mol−1 and the entropy was reduced by 6–9 J mol−1 K−1 relative to the N-acetyl derivatives (see Tables S2 and S3‡).
A small positive activation entropy is also known for the amide rotation in the proline derivative 5a, where ΔS≠ has been reported to be +9.0 for cis → trans and +12.8 J mol−1 K−1 for trans → cis rotations, respectively.31,32 The positive entropy of the amide rotation process may originate from the higher molecular rigidity of the ground state (featured by the n → π* interactions, for instance),17,29 as well as the stronger polar interactions of the ground state with water (organization of the solvent). For example, negative activation entropy for the amide rotation in Ac-AA-OMe has been reported for some proline analogues with additional polar groups,33 nonetheless, for some other polar proline analogues this value was also positive.34
![]() | ||
| Fig. 7 Amide rotation barriers of N-acyl derivatives. The values were determined in deuterium oxide solution (for AcGly-derivatives, buffered at pH 7 with potassium phosphate buffer) at 298 or 310 K. See ESI Table S4‡ for details. | ||
Importantly, for all three types of tested compounds, the rotation barrier difference between the Pdc derivatives was about 8–10 kJ mol−1 (a factor of 25–55 at 298 K), with the barriers always lower for the mPdc derivative. This result confirms that the rotation barrier difference in this case originates from the transition state structures. It also highlights the fact that in both the syn/exo and anti/endo transition states shown in Fig. 1C, the alkyl group of the acetyl moiety moves outwards and cannot sterically interfere with the substituents in the pyrrolidine ring (see also Fig. 2). On the other hand, the carbonyl group of the amide fragment moves towards the pyrrolidine ring and may sterically interfere with a carboxymethyl group in rPdc derivatives and create polar interactions with the carboxymethyl groups in both rPdc and mPdc derivatives.
C–N ↔ −O–C
N+]. The Reichardt–Dimroth parameter (ET) is a solvatochromic value that quantifies the ability of a solvent to separate intramolecular charges. Therefore, plotting rotation barriers in the E≠–ET coordinates produces positive slopes due to the differences in the ground state stabilities.16
While the experimental slopes (Fig. 8) for the rPdc and pyrrolidine derivatives (3a and 4a) were nearly identical (0.076 ± 0.006 and 0.075 ± 0.010, respectively), the slope for the mPdc derivative (2a) was notably larger (0.088 ± 0.008). This result illustrates differences in the polar effects of the transition states. We performed semi-empirical modelling of the amide rotation, and found that the overall molecular dipole follows the orientation of the amide carbonyl group C
O in the pyrrolidine derivative 4a, and in the rPdc derivative 3a. In the latter compound the dipolar contribution of both carboxymethyl-moieties is mutually compensated for by the C2-symmetry of the residue (Fig. 9). As a result, the overall molecular dipole is defined by the orientation of the amide carbonyl group, and effectively, both 3a and 4a yield the same slope in the E≠–ET coordinates.
![]() | ||
| Fig. 8 Amide rotation barriers for the N-acetyl derivatives 2a, 3a and 4a measured in different solvents by 1H EXSY NMR (at 298 or 310 K). For details see ESI Table S4.‡ | ||
![]() | ||
Fig. 9 Schematic representation of the dipolar effects during the amide rotation in N-acetylated compounds. log P values represent partitioning of the compounds between octan-1-ol and water at 298 K. | ||
The orientation of the molecular dipole in 2a is more complex. The configuration of the two carboxymethyl substituents renders the compound more polar than 3a. For example, the Rf of mPdc derivatives are consistently lower compared with the equivalent rPdc derivatives (see the Experimental section). Higher polarity of the meso-diastereomer is also inferred from the octan-1-ol/water partitioning values as shown in Fig. 9. The amide rotation via the putative syn/exo–syn/exo hybrid transition state partially compensates for the intramolecular dipole–dipole interaction in the transition state, causing the transition state to be relatively favored in the nonpolar solvents. As a result, the slope in the E≠–ET coordinates increases for 2a. Similar dipolar compensation can be expected for the amide rotation via the syn/exo transition state in peptidyl–prolyl fragments, and this effect might be enhanced due to a higher polarity of the amide fragment versus the esters reported here. As a result, the syn/exo transition state dominance should increase in less polar environments.
This model system allowed us to demonstrate experimental features of the amide rotation relevant for peptidyl–prolyl groups: (1) the transitions are enthalpically driven and (2) differ in dipole–dipole compensation levels in the transition states and (3) the syn/exo transition scheme prevails in general over the anti/endo. These conclusions are consistent with the results of computational studies, which suggested the dominance of the syn/exo over the anti/endo transition state. At the same time, the other theoretically possible transition states, anti/exo and especially syn/endo, are expected to be disfavoured in peptide structures, as these would require the movement of a bulky upstream peptidyl-fragment under the pyrrolidine ring of proline.
:
1) as an eluent (Rf = 0.48 for rPdc and 0.27 for mPdc derivatives). The 1H NMR spectra exhibit a characteristic AB system for the diastereotopic benzyl CH2 group in the rPdc derivative, whereas this signal is a singlet in the mPdc derivative, where the CH2 moiety is not diastereotopic. 1H NMR (500 MHz, CDCl3), δ: Bn-(MeO)mPdc-OMe, 7.44–7.32 (m, 5H), 4.05 (s, 2H), 3.68 (s, 6H), 3.59 (br, 2H), 2.18 (m, 4H); Bn-(MeO)rPdc-OMe, 7.44–7.32 (m, 5H), 4.11 (d, J = 13 Hz, 1H), 3.92 (br d, J = 7 Hz, 2H), 3.91 (d, J = 13 Hz, 1H), 3.75 (s, 6H), 2.42 (m, 2H), 2.04 (m, 2H). The removal of the N-benzyl group was performed by stirring the substance with palladium on charcoal (10%) in methanol acidified with an equivalent amount of hydrochloric acid, under hydrogen pressure (40 atm), at 35 °C. The reaction was monitored by NMR. After the reaction was complete, the mixture was filtered, the solvent was removed under reduced pressure, and the residue was dissolved in water, filtered with activated charcoal and freeze-dried. The amino acid hydrochlorides were obtained as pinkish powders. 1H NMR (700 MHz, CD3OD), δ: HCl·H-(MeO)mPdc-OMe, 4.33 (t, J = 6 Hz, 2H), 3.85 (s, 6H), 2.40 (m, 2H), 2.15 (m, 2H); HCl·H-(MeO)rPdc-OMe, 4.22 (t, J = 6 Hz, 2H), 3.80 (s, 6H), 2.32 (m, 2H), 2.06 (m, 2H).
:
1) for 48 hours at room temperature. The solvent was gently removed under nitrogen gas flow, and the residue was dissolved in water (0.7 ml) and freeze-dried. The crude material was purified on a short silica gel column (about 10 g) using an ethyl acetate–methanol mixture (19
:
1) as an eluent (Rf = 0.49 for 2a and 0.63 for 3a).
O), 173.8 (C
O), 173.3 (C
O), 61.2 (CH), 60.1 (CH), 53.1 (CH3O), 52.8 (CH3O), 29.4 (CH2), 27.7 (CH2), 21.2 (CH3). HMRS (ESI-orbitrap): calcd for [M + H]+ C10H16NO5 230.1023 Da, found 230.1028 Th. M.p. 60–65 °C.
O), 174.0 (C
O), 173.7 (C
O), 61.1 (CH), 59.5 (CH), 53.4 (CH3O), 53.0 (CH3O), 29.2 (CH2), 27.0 (CH2), 21.1 (CH3). HMRS (ESI-orbitrap): calcd for [M + H]+ C10H16NO5 230.1023 Da, found 230.1030 Th. M.p. 95–100 °C.
The proline derivative Ac-Pro-OMe (5a) was prepared by mixing commercial racemic HCl·Pro-OMe with an equivalent amount of acetic anhydride in dichloromethane solution. The solvent was removed under reduced pressure, and the crude material was purified using a hexane–ethyl acetate gradient elution (from 1
:
1 to 0
:
1). The substance was obtained as clear oil, which crystallized upon long time storage at 4 °C. White solid, m.p. 35–40 °C. 1H NMR (700 MHz, D2O), δ (two rotamers, Ktrans/cis 4.95 ± 0.05 at 298 K): 4.62 (s-cis, dd, J = 8.8, 2.6 Hz) and 4.36 (s-trans, dd, J = 8.8, 4.6 Hz, 1H, α-CH), 3.72 (s-cis, s) and 3.67 (s, s-trans, 3H, CH3O), 3.59 and 3.55 (s-trans, two m) and 3.46 and 3.39 (s-cis, two m, 2H, δ-CH2), 2.28 and 2.15 (s-cis, two m) and 2.23 and 1.94 (s-trans, two m, 2H, β-CH2), 2.05 (s-trans, s) and 1.93 (s-cis, s, 3H, CH3C
O), 1.95 (s-trans, m) and 1.90 and 1.80 (s-cis, two m, 2H, γ-CH2).
N-Acetyl pyrrolidine 4a was prepared as follows. Pyrrolidine (0.25 ml, 3 mmol) was mixed with acetic anhydride (0.29 ml, 3 mmol, 1 equiv.) in dichloromethane (5 ml) (CAUTION: vigorous reaction!). The mixture was stirred at room temperature for 30 min. Sodium carbonate (0.3 g) was added, and this mixture was shaken for about 15 hours at room temperature. The mixture was filtered, and the solvent was removed under reduced pressure to give 4a as a clear liquid. Yield 0.26 g (2.3 mmol, 75%). 1H NMR (500 MHz, D2O), δ: 3.41 (t, J = 7.0 Hz, 2H), 3.28 (t, J = 7.0 Hz, 2H), 1.97 (s, 3H), 1.86 (m, 2H), 1.79 (m, 2H).
:
1) as an eluent (Rf = 0.41 for 2b and 0.51 for 3b).
O), 174.3 (C
O), 174.0 (C
O), 62.8 (CH), 60.8 (CH), 53.0 (CH3O), 52.6 (CH3O), 39.2 (C), 31.1 (CH2), 26.7 (CH3), 25.4 (CH2). HMRS (ESI-orbitrap): calcd for [M + H]+ C13H22NO5 272.1492 Da, found 272.1492 Th.
O), 175.1 (C
O), 174.8 (C
O), 61.9 (CH), 60.9 (CH), 53.4 (CH3O), 52.8 (CH3O), 39.1 (C), 30.3 (CH2), 26.5 (CH3), 24.8 (CH2). HMRS (ESI-orbitrap): calcd for [M + H]+ C13H22NO5 272.1492 Da, found 272.1494 Th. M.p. 80–85 °C.
The synthesis of 5b was fully reproduced37 starting from 0.5 g of the racemic HCl·H-Pro-OMe (3.0 mmol, 1/10 scale of the reported synthesis). The product was a colorless oil. Yield 0.60 g (2.8 mmol, 92%). 1H NMR (700 MHz, D2O), δ (s-trans rotamer only, Ktrans/cis = 46 ± 2 at 298 K): 4.34 (dd, J = 8.5, 6.0 Hz, 1H, α-CH), 3.79 (dt, J = 10.6, 6.5 Hz, 1H, δ-CH), 3.74 (dt, J = 10.6, 6.7 Hz, 1H, δ-CH), 3.67 (s, 3H, CH3O), 2.14 (m, 1H, β-CH), 1.98 (m, 1H, γ-CH), 1.92 (m, 1H, γ-CH), 1.79 (m, 1H, β-CH), 1.17 (s, 9H, (CH3)3C).
Compound 4b was prepared as follows. Pyrrolidine (1 ml, 12.2 mmol) was mixed with triethylamine (3.4 ml, 24.4 mmol, 2 equiv.) in dichloromethane (20 ml). Pivaloyl chloride (1.6 ml, 13.0 mmol, 1.07 equiv.) was added dropwise at room temperature (CAUTION: vigorous reaction!), and the mixture was stirred for an additional 30 min. The solution was washed with water (1 × 5 ml), 1
:
2 aqueous hydrochloric acid (1 × 5 ml) and 1 M sodium hydrogen carbonate (1 × 5 ml); subsequently, it was dried over sodium sulphate, filtered and evaporated. Pure 4b was obtained as a colorless solid after crystallization from dichloromethane. Yield 1.80 g (11.6 mmol, 95%). 1H NMR (700 MHz, D2O), δ: 3.62 (t, J = 6.1 Hz, 2H), 3.34 (t, J = 6.7 Hz, 2H), 1.86 (m, 2H), 1.73 (m, 2H), 1.16 (s, 9H). M.p. 57–63 °C.
:
2 aqueous hydrochloric acid (1 × 0.5 ml) and 1 M sodium hydrogen carbonate (1 × 0.5 ml); subsequently, it was dried over sodium sulphate, filtered and evaporated. The resulting oil was extracted with water (1.5 ml). The extract was filtered through a short cation-exchange column (1.5 ml), and the elution was performed with water (about 5 ml). The aqueous fractions were collected and freeze-dried. Final purification was accomplished by silica gel chromatography (about 10 g) using an ethyl acetate–methanol gradient elution (19
:
1–4
:
1–0
:
1). The resulting compounds were colorless glassy materials and were slightly acidic. Therefore, the NMR samples in deuterium oxide were buffered with potassium phosphate buffer (70 mM, pH 7) and contained 35 mM of the analytes.
O), 173.8 (C
O), 172.8 (C
O), 170.5 (C
O), 60.6 (CH), 59.9 (CH), 53.2 (CH3O), 52.9 (CH3O), 41.5 (NH-CH2), 29.4 (CH2), 27.2 (CH2), 21.6 (CH3). HRMS (ESI-orbitrap): calcd for [M + H]+ C12H19N2O6 287.1238 Da, found 230.1240 Th.
O), 174.0 (C
O), 173.5 (C
O), 170.4 (C
O), 60.0 (CH), 59.7 (CH), 53.5 (CH3O), 53.0 (CH3O), 41.5 (NH-CH2), 29.3 (CH2), 26.4 (CH2), 21.5 (CH3). HRMS (ESI-orbitrap): calcd for [M + H]+ C12H19N2O6 287.1238 Da, found 230.1239 Th.
The other two peptides 4c and 5c were prepared under analogous conditions, by employing N,N,N′,N‘-tetramethyl-O-(N-succinimidyl)uronium hexafluorophosphate (HSTU) as the coupling reagent. 4c, 1H NMR (700 MHz, D2O), δ: 3.93 (s, 2H), 3.40 (t, J = 6.9 Hz, 2H), 3.33 (t, J = 6.9 Hz, 2H), 1.98 (s, 3H), 1.90 (m, 2H), 1.80 (m, 2H). HRMS (ESI-orbitrap): calcd for [M + H]+ C8H15N2O2 171.1128 Da, found 171.1129 Th. 5c, 1H NMR (700 MHz, D2O), δ (s-trans rotamer only, Ktrans/cis = 6.15 ± 0.27 at 298 K): 4.41 (dd, J = 9.0, 4.4 Hz, 1H, α-CH), 4.03 (d, J = 17.3 Hz, 1H), 3.98 (d, J = 17.5 Hz, 1H), 3.68 (s, 3H, CH3O), 3.58 (m, 1H, δ-CH), 3.54 (m, 1H, δ-CH), 2.22 (m, 1H, β-CH), 1.98 (s, 3H), 1.97 (m, 2H, γ-CH2), 1.96 (m, 1H, β-CH). HRMS (ESI-orbitrap): calcd for [M + H]+ C10H17N2O4 229.1183 Da, found 229.1182 Th.
P values: 2a −0.74 ± 0.05, 3a −0.43 ± 0.03, 4a −0.20 ± 0.07, 5a −0.44 ± 0.05.
Footnotes |
| † The article is dedicated to the memory of Prof. Yuriy Kholin (1962–2017), who nurtured a generation of young talented Ukrainian chemists. |
| ‡ Electronic supplementary information (ESI) available: Tables S1–S5, copies of the NMR spectra for compounds. CCDC 1554475 and 1554476. For ESI and crystallographic data in CIF or other electronic format see DOI: 10.1039/c7ob01421j |
| This journal is © The Royal Society of Chemistry 2017 |