Bryden A. F.
Le Bailly
,
Liam
Byrne
,
Vincent
Diemer
,
Mohammadali
Foroozandeh
,
Gareth A.
Morris
and
Jonathan
Clayden
*
School of Chemistry, University of Manchester, Oxford Road, Manchester M13 9PL, UK. E-mail: clayden@man.ac.uk
First published on 21st January 2015
Although foldamers, by definition, are extended molecular structures with a well-defined conformation, minor conformers must be populated at least to some extent in solution. We present a quantitative analysis of these minor conformers for a series of helical oligomers built from achiral but helicogenic α-amino acids. By measuring the chain length dependence or chain position dependence of NMR or CD quantities that measure screw-sense preference in a helical oligomer, we quantify values for the decay constant of a conformational signal as it passes through the molecular structure. This conformational signal is a perturbation of the racemic mixture of M and P helices that such oligomers typically adopt by the inclusion of an N or C terminal chiral inducer. We show that decay constants may be very low (<1% signal loss per residue) in non-polar solvents, and we evaluate the increase in decay constant that results in polar solvents, at higher temperatures, and with more conformationally flexible residues such as Gly. Decay constants are independent of whether the signal originates from the N or the C terminus. By interpreting the decay constant in terms of the probability with which conformations containing a screw-sense reversal are populated, we quantify the populations of these alternative minor conformers within the overall ensemble of secondary structures adopted by the foldamer. We deduce helical persistence lengths for Aib polymers that allow us to show that in a non-polar solvent a peptide helix, even in the absence of chiral residues, may continue with the same screw sense for approximately 200 residues.
Oligomers of achiral quaternary amino acids, of which 2-aminoisobutyric acid (Aib) is the archetype, but which also include 1-aminocyclohexanecarboxylic acid (Ac6c, Fig. 1) along with 1-aminocyclopentanecarboxylic acid (Ac5c), diethylglycine (Deg) etc.22,23 form a class of foldamers that adopt specific compact conformations that take the form of hydrogen-bonded 310 helices.22,24–28 These helices are necessarily equal mixtures of left- and right-handed screw-sense conformers because (unlike natural peptides or nucleic acids) the monomers themselves are achiral.29 Nonetheless an equal population of the M and P conformers may be detectably perturbed by incorporation of one or more terminal chiral residues30–39 or interaction with chiral ligands.40–44 A helical chain of achiral monomers may be switched between a left- and a right-handed screw-sense preference either by configurational inversion of a single chiral residue at one terminus,45 or by reversible covalent40 or non-covalent46,47 binding of a chiral ligand to a terminal binding site. Such foldamers exhibit a form of switchability analogous to membrane-bound receptors, in which a conformational response to the binding of a ligand mediates signal transduction through the thickness of the cell membrane to mediate a remote biochemical outcome.13
Fig. 1 Remotely induced diastereoselective formation of diastereoisomeric foldamers 2a–d and the solvent-dependent chemical shift separation between their CH signals. |
We have exploited this ability of helical foldamers to communicate information by propagation of a conformational preference in some reactions exhibiting remote stereochemical control.48 Two chiral L-α-methylvaline (α-MeVal) residues at the terminus of a peptide-like oligomer of achiral quaternary amino acids are sufficient to induce essentially complete control over screw-sense preference which persists far enough to allow reactions to take place with 1,61 asymmetric induction and 88:12 diastereoselectivity. Similar reactions over shorter distances (1,46, 1,31 and 1,16 asymmetric induction) show successively increased levels of stereoselectivity. This chain-length dependence49 indicates that the ability of the helical foldamer to communicate stereochemical information is high, but not perfect (Fig. 1).
We may consider a conformational perturbation such as the screw-sense preference induced by the chiral N-terminal L-α-MeVal residues of 1 as a signal50–52 that propagates through the polymeric molecule without erosion only if that molecule truly adopts one single ‘specific compact conformation’.1 Any random intrusion of alternative conformations must attenuate the signal, allowing the local conformational preference within the helix to decay towards an equal population of left- and right-handed helices the further the polymer is extended away from the chiral terminus.
It is the occurrence of these alternative minor foldamer conformations in the solution state that we explore in this paper. We use NMR and other spectroscopic methods to quantify the attenuation of a conformational preference with increasing distance from a source of screw-sense induction located at either terminus of a series of Aib-containing foldamers. We show that a conformational signal is in many cases remarkably persistent, but its decrease in intensity on moving in either direction through the polymer may be characterized as a solvent-dependent exponential decay. We interpret the spatial decay constant of the signal in terms of the stochastic intrusion of screw-sense reversals into the three-dimensional structure of the oligomers, and use the results to quantify the population of minor contributors to the overall ensemble of conformations that describe the dynamic three-dimensional structure of the foldamer.
den = hei × (1 − f)n × ds | (1) |
We deduce that de decays by a factor of 0.8% per residue (i.e. f = 0.008) under the conditions of the reaction in THF at −50 °C.58–60 Given that when n = 0, dei = hei × ds, we can use de0 (extrapolated from Fig. 2a) = 90% and the reported value of he0 for (L-α-MeVal)2 of 99% (ref. 31) to deduce that the inherent diastereoselectivity of the reaction corresponds to a ds of ca. 95:5.
Constraints of temperature and solvent, which must be chosen to be compatible with the reaction in question, make diastereoselectivity an inconvenient observable quantity for the wider study helix of conformation. We therefore turned to the anisochronicity Δδ separating the signals arising from the major and minor diastereoisomers of the products of this addition reaction, which also allowed us to measure the conformational preference of the oligomers in alternative solvents.
We previously described a simple spectroscopic method for quantifying the relative proportions of the two screw-sense conformers, the helical excess (he), by measuring the anisochronicity of a pair of signals arising from diastereotopic reporter groups embedded in the helical chain.34 In the fast exchange regime, the chemical shift separation of two such anisochronous peaks is proportional to the local he experienced by the NMR reporter group, such that
he = Δδfast/Δδslow | (2) |
The same arguments must hold for pairs of signals arising from corresponding groups in diastereoisomeric structures. The detailed factors governing their anisochronicity will be more complicated to analyse, but these detailed conformational considerations are irrelevant to the extraction of an exponential decay factor. A decrease in the fast exchange regime anisochronicity of peaks arising from two diastereoisomers on increasing the length of the achiral segment separating the stereogenic centres that give rise to their diastereoisomeric relationship may be interpreted simply as a measure of that segment's decreasing ability to communicate stereochemical information between the relevant stereogenic centres.
Fig. 1 shows how the separation Δδ (= Δδfast of eqn (2)) between the CHN signals of the diastereoisomers (S,S,S) and (S,S,R)-2a–d decreases with increasing chain length at 23 °C in deuterated THF and at +23 °C or 40 °C in a more polar solvent, deuterated ethanol, and Fig. 2b plots anisochronicity against n (the length of the chain lying between the centres) for these solvents. Poor solubility prevented us studying these oligomers in chloroform. We fitted the points to a modification of the formula used earlier
Δδn = Δδi × (1 − f)n | (3) |
An alternative, simpler measure of signal decay may be obtained by inserting several NMR reporters, each containing a pair of diastereotopic groups giving rise to potentially anisochronous signals, at successive positions along a helical oligomer. Provided that the reporter groups are sufficiently far from the ends of the chain for their anisochronicity to be dictated solely by their local helical environment, the decrease in their anisochronicity along the chain gives a measure of signal decay. The challenge in such an approach is identifying which of several identical reporters is which. This is less difficult if the rate at which anisochronicity decreases is relatively high, because the signals from successive reporters are not superimposed upon one another, and a simplifying assumption can be made that anisochronicity can only decrease along the chain. We applied this principle to oligomer 3 (Fig. 3)65 which contains three glycine residues. These serve the dual purpose of firstly introducing greater flexibility than Aib, and therefore a steeper rate of decay of screw-sense preference, and secondly of providing pairs of diastereotopic protons as spectroscopic reporters of screw-sense preference. In ethanol and methanol, the three glycine ABX systems were resolved by Lorentz–Gauss transformation, while in CDCl3 and in CD2Cl2 poor dispersion prevented reliable assignment. In THF, however, the A and B coupled signals of the three ABX spin systems were identified using an absorption-mode 2D J-resolved experiment based on the PSYCHE method.66,67Fig. 4 plots the fast exchange anisochronicity Δδ between the A and B components of the AB systems68 against the position n of the relevant Gly residue relative to the N terminal Cbz-L-Phe at position i, and in each case a value for f was obtained using eqn (3).
Fig. 4 Anisochronicity Δδ in the three Gly residues of 3 in THF, C2D5OD (23 °C and 40 °C) and CD3OD. THF at 23 °C, f = 0.026; ethanol at 23 °C, f = 0.078; methanol at 23 °C, f = 0.14. |
This time, the additional conformational flexibility complicates the analysis of the decay constants. Gly introduces much more flexibility than Aib into a peptide chain,17 so we expect the values of f for Aib (fAib) and for Gly (fGly) to differ significantly. There are four Aib residues for every Gly in the chain, so the value f obtained from fitting eqn (3) to Fig. 4 is therefore:
f = 0.8fAib + 0.2fGly | (4) |
The values of f obtained exhibit a solvent dependence consistent with that noted for 2, with the more polar methanol allowing greater conformational flexibility than the less polar THF, and ethanol lying somewhere between. Using our previous estimates of values for fAib at +23 °C in THF, in ethanol (Fig. 2) and in methanol,31,34 we can further deduce values of fGly from eqn (4) indicating that glycine residues introduce a fall in Δδ at +23 °C of 13% in THF, 31% in C2D5OD or 49% in CD3OD, confirming that the Gly residues introduce significantly more conformational flexibility than the Aib residues.
We used this approach to explore the communication of conformational information through an Aib oligomer from the C terminus to the N terminus. The multiply labelled compounds 7 were synthesized using the strategy illustrated in Fig. 5. Labelled and unlabelled Aib trimers 4** and 4 were built by standard methods from a starting monomer that was either 100% doubly (4**) or 1.1% singly (4, i.e. natural abundance) 13C labelled in its Aib methyl groups. The two samples of trimer 4** and 4 were then mixed in ratios of 50:50 or 25:75 and these two samples sequentially ligated to the methyl ester of 4** to make a nonamer 5, three of whose residues are labelled with different abundances of 13C, as shown in Fig. 5. 5 was extended by ligation with an unlabelled Aib dimer, trimer or tetramer to make the three oligomers 6a–c, which were each coupled at the C terminus to the powerful chiral controlling residue L-AlaNHt-Bu32 to provide 12-mer 7a, 13-mer 7b and 14-mer 7c.
These three ‘frame-shifted’ oligomers between them contain a labelled probe at every one of the nine positions in the chain from i − 3 to i − 11, starting from the C terminal L-AlaNHt-Bu chiral controller at position i (Fig. 6). Only 7a was soluble in pure methanol, but all were soluble in a 9:1 mixture of methanol and chloroform. The relevant portions of the three 13C spectra are shown in Fig. 7. The abundance-labelled labels are evident in the successively smaller anisochronicities of pairs of peaks with successively smaller integrals (100:50:25). Similar spectra in THF and in CDCl3 showed significant peak overlap due to the much lower decay of anisochronicity in these solvents, but individual peaks were still evident in THF. In CDCl3 Lorentz–Gauss transformation assisted assignment of the 50% labelled signals, but the 25% labelled signals could not be located reliably.
Fig. 6 Multi-labelled oligomers 7a–c allow measurement of screw sense at each of the nine positions in the oligomer from i − 3 to i − 11. |
Fig. 7 Illustrative 13C NMR spectra of 7a (top), 7b (middle) and 7c (bottom) in 9:1 methanol:chloroform, with peaks arising from the 100%, 50% and 25% labels indicated by colours used in Fig. 6. |
Fig. 8 shows a plot of anisochronicity at each of the positions in the chain against the position (−n) of the label in the chain relative to the C-terminal alaninamide (residue i). Fitting an exponential decay (eqn (3), using positive values of n) to each plot gives signal decay constants f (the per-residue decrease in Δδ) of 6.1% in methanol, 5.4% in 9:1 methanol:chloroform, and 0.5% in THF. Earlier work had shown that aggregation effects may become important in solvents such as chloroform,69 but these were shown to be negligible in this case as the anisochronicity was independent of concentration over a 4-fold concentration range (4.6–18 μM). Changes in temperature induced a small variation in the value of f for the oligomers in methanol, which was 5.5% at 0 °C; 6.1% at 23 °C and 6.5% at 40 °C. These values correspond to an enthalpy difference of about 3 kJ mol−1 and an entropy difference of −13 J mol−1 K−1.
A further group of compounds 8a–d (Fig. 9), developed as part of a related project,70 had a better solubility profile than 7, and allowed comparisons of induced helical excess through four, eight, eleven or twelve achiral residues in methanol, chloroform and THF, as well as correlation with circular dichroism data. 1H NMR chemical shift separations Δδ between Ha and Hb at the two equivalent methylene groups of the azepine rings in 8a–d compounds at +23 °C are shown plotted against −n in Fig. 10, and an exponential decay (eqn (3)) fitted to these separations gave signal decay constants f of 1.3% in THF and 7.7% in MeOH. In CDCl3 the signals unfortunately showed no anisochronicity.
Fig. 9 Dibenzazepinocarbonyl-capped oligomers with screw-sense induction from a C terminal L-Ala residue. |
In the circular dichroism spectra of 8a–d the chromophore associated with the azepine biaryl (θmax at 249 nm) is well separated from those associated with the carbonyl groups in the chain (θmax < 220 nm), and thus we may expect molar ellipticity at 249 nm to be proportional71 to the induced helical excess at the N terminus of the chain. Fig. 10b shows a plot of this value against chain length. From a modification of eqn (3)
θn = θi × (1 − f)n | (5) |
The decay of screw-sense preference on lengthening a helical domain interposed between a screw-sense controller and reporter may be viewed in terms of the chance of helix reversals occurring in that domain. A domain containing a single reversal, or indeed any odd number of reversals, will give rise to opposite screw-sense preferences at the two termini of the domain, while zero or an even number of reversals will lead to identical screw senses at the termini. Thus a helix having P screw sense at one terminus will have P screw sense at the other if it contains zero or an even number of reversals and M screw sense if it contains an odd number of reversals. The likelihood of the two ends of the domain having the same screw sense may be calculated by considering the excess probability of an odd number over an even number of screw sense reversals, assuming an equal probability of screw-sense reversal at each location in the domain.
For a domain of N elements, the difference in probability between having an even number of elements in state α and an odd number in state α is65
(2p − 1)N | (6) |
hen = hei × (2p − 1)n | (7) |
Thus the empirical values for the signal decay constant f derived from eqn (3)–(5) may be interpreted structurally as a measure of the probability of screw-sense reversal, with f = 2 × (1 − p). In other words with p as the screw-sense fidelity, f/2 may be interpreted as the probability of a helix reversal occurring, and f, the signal decay constant, as the fractional decrease in helical excess as the screw-sense preference passes from one monomer to the next. Values for decay constant f and hence fidelity p obtained from Fig. 2, 4, 8 and 10 are assembled into Table 1.
Compound series | Monomer | Solvent | T/°C | Per-residue decrease in he fa/% | Screw-sense fidelity p | Average length of uniform screw sense L |
---|---|---|---|---|---|---|
a Determined by NMR unless otherwise stated. b Determined by circular dichroism. | ||||||
2 | Aib/Ac6c | THF | −50 | 0.8 | 0.996 | |
2 | Aib/Ac6c | THF | +23 | 0.1 | 0.9996 | |
7 | Aib | THF | +23 | 0.5 | 0.9975 | 200 |
8 | Aib | THF | +23 | 1.3 | 0.993 | |
8 | Aib | THF | +23 | 2.5b | 0.988 | |
2 | Aib/Ac6c | EtOH | +23 | 1.9 | 0.991 | 50 |
2 | Aib/Ac6c | EtOH | +40 | 2.5 | 0.987 | |
7 | Aib | 9:1 MeOH:CDCl3 | +23 | 5.4 | 0.973 | 18 |
7 | Aib | MeOH | 0 | 5.5 | 0.972 | 18 |
7 | Aib | MeOH | +23 | 6.1 | 0.970 | 16 |
7 | Aib | MeOH | +40 | 6.5 | 0.967 | 15 |
8 | Aib | MeOH | +23 | 7.7 | 0.962 | |
8 | Aib | MeOH | +23 | 7.6b | 0.962 | |
3 | Gly | THF | +23 | 13 | 0.935 | — |
3 | Gly | EtOH | +23 | 31 | 0.845 | — |
3 | Gly | MeOH | +23 | 49 | 0.76 | — |
This interpretation of the decay of he along a chain in terms of the likelihood of the two termini of a helical domain having the same or opposite screw sense makes it clear that the decay constant f for any monomer must be independent of whether the ‘signal’ is passing from the N to the C (as in 2 and 3) or the C to the N (as in 7 and 8) terminus. The values of p obtained for Aib in MeOH from 7 or 8 for a signal travelling from the C terminus to the N terminus (p = 0.970, 0.962, 0.962) are close to values previously calculated for a signal travelling from the N terminus to the C terminus (0.9735 (ref. 31,34,96), 0.965 (ref. 65)). It is also notable that the value of p obtained for Gly from 3 corresponds closely to a previous observation17,97 that Gly conducts a helical preference 59% as efficiently as Aib (i.e. (1 − fGly)/(1 − fAib) = 0.59; the corresponding value calculated from the data in Table 1 is 0.51/0.939 = 0.54).
The values in Table 1 clearly show that helical reversals are more likely at higher temperatures, and in more polar solvents. It may be deduced that whatever disruption of the hydrogen-bonding network is required to allow a helix reversal is more favourable in polar than non-polar solvents. The question of what structures these helix reversals adopt is still under investigation, but one important observation suggests that helix reversal does not require the formation of an unstructured random coil conformation. In such a conformation, oligomer NH groups would be exposed to the deuterated solvent and would therefore undergo hydrogen–deuterium exchange.98,99 Yet both the observed rate at which the NH signals of the oligomer exchange with D from the deuterated solvent, and the rate at which protonated signals returned on dissolution of deuterated oligomer in protonated solvent, are extremely slow. H/D exchange in CD3OD occurs on a time scale of hours to days at ambient temperature, while helical inversion occurs on the millisecond time scale,29,34,62–64 suggesting that random coil conformations make a negligible contribution to the ensemble of structures adopted by Aib oligomers. Thus we propose that the helical reversals that lead to signal decay are well-defined, local structural features which entail only a small perturbation of the 310 helical structure, and retain a well-defined intramolecularly hydrogen-bonded network. We are currently exploring experimentally the possible detailed structures for these helix reversals. As suggested by Kubasik,63 migration of a well-structured helical reversal along a chain, sequentially breaking and reforming only one or maybe two hydrogen bonds at a time, could provide a mechanism for the remarkably rapid inversion of a 310 helix.
The values also allow for the first time an estimate to be made, using the one-dimensional Ising model,95,100 of the average length over which hypothetical extended Aib homo-polymers might be expected to adopt a uniform screw sense. The Ising model, first used to describe ferromagnetic and antiferromagnetic arrays of atoms, assumes that the energetic cost of a state reversal (here a change in screw sense) in a linear array is determined solely by nearest neighbours. If the chance of a helix reversal at any one residue is (1 − p), the Boltzmann relation gives
(1 − p)/p = e−2J/kT | (8) |
According to the one-dimensional Ising model, the helical excess will decay exponentially over n segments as the nth power of tanh(J/kT). This factor is equivalent to our empirical decay factor (1 − f) (eqn (3)), giving as expected f = 2 × (1 − p). The persistence length L (the average distance between reversals) may be obtained from the Ising model by101
L = −1/ln[tanh(J/kT)] = −1/ln(2p − 1) = −1/ln(1 − f) | (9) |
Values for L are shown in Table 1 for selected data. A sufficiently soluble Aib polymer, in a non-polar solvent, can evidently be expected to maintain a preferred screw sense over distances measurable in hundreds of monomers.
One advantage of the Ising model is that it allows a simple estimate, from eqn (8), of the energetic cost 2J of a helix reversal. For example, in the case of 2a–d in THF at −50 °C (Fig. 2a), J = 10 kJ mol−1; for 7a–c in methanol at +40 °C it is 9 kJ mol−1, both considerably less than the typical strength of a hydrogen bond. This is consistent with the observation noted earlier that hydrogen–deuterium exchange is very slow in such systems, providing strong evidence that helix inversion simply involves a rearrangement of intramolecular hydrogen bonds rather than any reduction in their number. The Ising model assumes that the driving force for the retention of helicity is purely energetic; as the analysis of the temperature dependence of f for 7a–c in methanol in the preceding section shows, a substantial part of the driving force is in fact entropic.
Footnote |
† Electronic supplementary information (ESI) available: Synthesis and characterisation of all new compounds. See DOI: 10.1039/c4sc03944k |
This journal is © The Royal Society of Chemistry 2015 |