Melissa A.
Reith
a,
Sinan
Kardas
bc,
Chiel
Mertens
a,
Mathieu
Fossépré
b,
Mathieu
Surin
*b,
Jan
Steinkoenig
*a and
Filip E.
Du Prez
*a
aPolymer Chemistry Research Group, Centre of Macromolecular Chemistry (CMaC), Department of Organic and Macromolecular Chemistry, Faculty of Sciences, Ghent University, Krijgslaan 281 S4-bis, Ghent B-9000, Belgium. E-mail: filip.duprez@ugent.be; jan.steinkonig@ugent.be
bLaboratory of Chemistry of Novel Materials, Center of Innovation in Materials and Polymers (CIRMAP), University of Mons - UMONS, Place du Parc 20, Mons B-7000, Belgium. E-mail: mathieu.surin@umons.ac.be
cInstitute for Complex Molecular Systems, Eindhoven University of Technology, P.O. Box 513, Eindhoven 5600 MB, The Netherlands
First published on 21st June 2021
Macromolecules found in Nature display a precise control over the primary as well as higher ordered architectures. To mimic the folding found in Nature, we herein demonstrate the design and characterization of single-chain nanoparticles that are formed by the folding of sequence-defined macromolecules with metal ions. The study showcases the influence of the loop size of such precision macromolecules on their relative hydrodynamic radius. The sequence-defined structures are fabricated using thiolactone chemistry, where two picolyl moieties are installed forming a valuable ligand system for subsequent metal complexation. Next, metal ions such as Ni(II) and Cu(II) ions are introduced to fold the unimers into sequence-defined single-chain nanoparticles (SD-SCNPs). After proving the successful complexation using a trimer, a systematic study is conducted altering the distance between the respective ligands by incorporating variable numbers of non-functionalized spacer units. Finally, the loop size formation of the SD-SCNPs is evidenced by DOSY measurements. The result indicates that the positioning of the ligands plays a crucial role on the compaction process and, more specifically, on the final size of the SD-SCNP. In addition, molecular dynamics (MD) simulations show the effects of the sequence and Ni(II) complexation on the structure and compaction of the SD-SCNPs, and highlight the differences of the nanoparticles’ shape when varying the number of spacer units. Finally, the system is further expanded to a dodecamer and even a heptadecamer with drastically decreased hydrodynamic radii after compaction.
Thus, sequence-defined macromolecules, a vastly advancing field in polymer science, have been developed to attain Nature's precision synthetically with the possibility to extend the applicable monomeric library to non-natural compounds. For instance, the preparation of highly defined artificial biomacromolecules was assisted by earlier Nobel-Prize-awarded methodology such as the Merrifield synthesis. By adapting the protocols, sequence-defined synthetic macromolecules can be prepared using solid-phase approaches.3,7 Furthermore, other routes were developed, such as liquid-phase synthesis2,8 and single unit monomer insertion (SUMI)9–12 where diverse chemical platforms allow for efficient monomer installation. These monodisperse and precise macromolecules are used in bioinspired applications such as information storage13,14 and encryption,15,16 as well as enzymatic-like catalysis.17 Recent reviews are highlighting the advances on the field of discrete synthetic macromolecules focusing on their properties and applications.18–20
Earlier attempts to mimic protein folding were performed using linear polymeric chains with dispersity ranging from approx. 1.121 up to 2.23.22 The introduction of functional moieties along the backbone enables intramolecular compaction upon a stimulus. This subclass of materials referred to as single-chain nanoparticles (SCNPs), emerged in the early 2000s.23,24 A diverse library of folding chemistries including irreversible covalent linkage formation25–27 reversible covalent bonds,28–30 supramolecular interactions,31,32 and degradable (self-immolative) bonds33 was described to form SCNPs. Due to their specific characteristics, SCNPs find application in many fields, such as drug delivery,34–37 catalysis,25,38,39 imaging,40–43 sensors44 and responsive materials.45,46 The morphology and dynamics of these nanostructures have been revealed by combining experimental and molecular simulation studies, mainly in solution.47–49 In particular, it has been reported that the sparse morphology of current SCNPs in solution resembles that of intrinsically disordered proteins.49 Consequently, several strategies have been proposed for controlling the SCNP folding49–51 but the tuning of their internal structure remains an unsolved issue. Nevertheless, the morphology of SCNPs in solution is expected to influence their activity and selectivity, in particular for catalysis applications.25,51,52
In addition to SCNPs, other monodisperse structures have been studied, more specifically foldamers and bioactive macrocycles. Macrocycles such as cyclic peptides are the first class combining precise primary structure with a higher ordered architecture. Here, loops are formed mainly via head-to-tail cyclization,53 but they can also emerge from side chain-to-side chain,54 head-to-side chain or side chain-to-tail.55 Cyclic peptides have found broad applications from drug delivery56 to nanomaterials.57 However, macrocycles are not limited to natural or biological macromolecules but cyclization of linear synthetic sequence-defined macromolecules has also been studied, as highlighted in recent papers by the groups of Alabi58 and Meier.59
The second important class of monodisperse, synthetic macromolecules to be highlighted in this section is the one of foldamers, which collapse spontaneously into a secondary structure (e.g. helix). Foldamers can be derived from biotic (such as peptoids60 or peptide nucleic acids61) and abiotic (e.g. oligo-phenyl-ethylenes62 and aromatic oligoamides63) building blocks. In the work of Zuckermann, the binding affinity of a peptoidic two-helix bundle towards zinc was studied.64
In the current work, we will focus on the metallo-induced folding of synthetic sequence-defined structures for which we investigated the impact of the distance between two functional handles on the folding behaviour. The relevance of such a study was outlined in different perspectives. In 2011, Lutz and Sawamoto emphasized the importance of controlling primary and secondary structures.65 More recently, Meier and Barner-Kowollik reiterated the influence of spatial precision in macromolecular design.66
A careful scan of existing synthetic protocols for the design of sequence-defined macromolecules with the option for a precise installation of suitable moieties for subsequent folding resulted in the choice for the well-known thiolactone-based protocol.7 Here, picolyl moieties could be easily installed, forming a valuable ligand system for subsequent metal complexation (vide infra). Thus, the obtained precisely designed single-chain nanoparticles allowed us to investigate the influence of different loop sizes – being the distance between two functional cross-linking units in a chain – on the size of the formed nanoparticle. The synthesis is outlined as follows: (i) an iterative two-step solid-phase approach using thiolactone chemistry to synthesize sequence-defined macromolecules; (ii) a selected metal-induced intra-molecular complexation to compact the structures and investigate their behaviour (Scheme 1).
Taking inspiration from the work of Chattopadhyay and co-workers,67 a picolyl moiety was chosen to design a N,S-ligand, which can construct stable N2S2 donor sets with different metal species (e.g. Cu(II) and Ni(II)68 ions). Subsequently, we studied the effect of an increasing loop size on the folding behaviour by incorporating from one up to fifteen spacer units between the functional handles.
Fig. 2 (a) DOSY (500 MHz) spectrum of S1,7. (b) DOSY (500 MHz) spectrum of its folded structure. Dashed lines display diffusion coefficients of the corresponding macromolecules. |
Fig. 2a showcases the DOSY spectrum of sequence S1,7. All other DOSY measurements (i.e. S1,2–S1,6) can be found in the ESI.† The dashed line highlights the diffusion coefficient of the macromolecule. Here, the Stokes–Einstein relation71 provides an hydrodynamic radius (Rh) of 1.36 nm. Comparing the diffusion traces’ maxima of Fig. 2a and b, a shift of the diffusion coefficient upon the complex formation was observed. The shift resulted in a decrease of the hydrodynamic radius to 1.30 nm, indicating the intramolecular folding and thereby the formation of a sequence-defined oligomeric single-chain nanoparticle.
To emphasize folding effects for all sequences, we calculated the compaction ratio, i.e. the ratio between hydrodynamic radii of both folded and non-folded sequences (see Fig. 3, eqn (I)). The hydrodynamic radius’ numerical value of the folded oligomer was divided through the numerical value of its non-complexed analogue to obtain the size decrease as a percentage. It was expected that with an increasing number of spacer units between the picolyl handles, a decrease of the hydrodynamic radius would be observed, as the end-to-end vector would be decreased with an increasing loop size. Experimentally, a minimum value of the compaction ratio was observed for an intermediate number of spacer units (see Fig. 3). Most significantly, the compaction ratio of S1,2 (i.e. both functional moieties in close proximity) was found the highest of all heptamers (8%). Thus, the formed SCNPs were the smallest. Then, the compaction ratio decreased from 8% (for S1,2) to 3% (for S1,4) and increased again to 5% (for S1,7). The experiments were repeated, and the results showed a reproducibility with an error below 1%. An explanation for the minimum was not found.
In addition to the set of heptamers with two picolyl functionalities, an oligomer with only one picolyl moiety (S7) was synthesized to underpin that both picolyl moieties are necessary for the intramolecular folding. As previously established, Ni(ClO4)2 was added to a 0.01 molar methanol solution of the sequence. The Rh obtained through DOSY measurements increased after the addition of the Ni(II) ions from 1.28 nm to 1.32 nm, indicating the reaction of the ligand functionality with the metal ion but not the formation of SCNPs (see Fig. S36†). Since dimerization could in principle also occur with this oligomer, a reference sequence without any picolyl functionalities was also synthesized with the aim to demonstrate that only those units, and thus no other functionalities along the oligomeric structure, are responsible for the folding. However, this hypothesis could not be confirmed via a similar DOSY experiment as a result of the insolubility of this oligomer in methanol.
Moreover, the SCNP was studied during temperature dependent DOSY measurements (Fig. S37†) revealing that the complexed form unfolds upon heating to 100 °C. For the study, the complex was dissolved in DMF and heated in steps of 20 °C from 25 °C to 100 °C. A DOSY spectrum was recorded at each equilibrated temperature. Afterwards, the system was cooled down in the same intervals, scanning a second spectrum for each of those steps. After the heating cycle was completed, an increase in the size of the SCNP was observed, indicating the expansion of the SCNP and thus the reversion of the complexation. In addition to the DOSY measurements, the successful complexation was confirmed using UV–Vis measurements (Fig. S38†) and 1H NMR spectra (Fig. S12†) as described for the model study using the trimer. Furthermore, the complexed structure was observed during high resolution ESI-MS measurements (Fig. S40†). Observing the heptamer in HR-ESI MS was more challenging. A small excess of the Ni(II) ions was therefore added to shift the equilibrium towards the complexed form. Consequently, additional salt adduct formation could be observed as indicated in Table S13.† Further, two species were observed in a LC-ESI-MS spectrum, similar to the trimer study (Fig. S6†). Finally, the oligomer length was further increased up to a dodecamer (sequence S1,12) and heptadecamer (sequence S1,17) to showcase the folding for larger structures. The previously described protocol was adapted to ensure full conversions by prolonging the reaction time of the elongation steps. The picolyl functionalities were placed at both ends of the oligomer. The folding, again induced with Ni(ClO4)2 resulted in a compaction ratio of 18% for the dodecamer and 24% for the heptadecamer (Fig. 3). The observed ratios indicated that with an increasing length of the oligomers, the value for the compaction ratio increased as well and was comparable to the ones obtained for disperse polymeric single-chain nanoparticles, e.g. 13% observed from Murray and Fulton30 or 22% recorded by Steinkoenig et al.72
Fig. 4 (a) Chemical structure of the model sequence S1,7. Molecular model of the heptamer is constructed by joining nine fragments. The corresponding fragment numbers are also shown. The functional handles are depicted in magenta. (b) Left: last snapshots of MD simulations of S1,7 for non-complexed (top) and complexed (bottom) structures, in sticks representation. The two functional handles are shown in magenta, with nitrogen and sulphur atoms in blue and yellow, respectively; the nickel ion is depicted as a cyan sphere. (b) Middle: corresponding 2D distance maps (for the fragment number, see Fig. 4a). (b) Right: the 2D difference distance map is obtained by subtracting the 2D distance map of the non-complexed sequence from that of the complexed sequence. The regions with large conformational changes, corresponding to complexation sites, are represented by large distances (yellowish regions) and marked with magenta squares. (c) 2D difference distance maps of the different sequences. From left to right: S1,2/S1,4/S1,7/S1,12/S1,17. The regions with large conformational changes, corresponding to complexation sites, are represented by large distances (yellowish regions) and marked with magenta squares. |
Further, the change in conformation upon compaction was supported by the corresponding Root-Mean-Square Fluctuation plots (ESI Fig. S47†), which showed a bell-shaped profile, characteristic of flexible polymeric systems. Therefore, the non-complexed heptamers were highly flexible at the termini and unstructured in the middle of the chain, while the complexation with the nickel ion reduced their conformational freedom, leading to a more structured overall arrangement at the functional picolyl handles.
To evaluate the effect of the complexation on the compactness of heptamers, the radii of gyration (Rg) were estimated from MD simulations (Fig. 5a). It is important to note that the hydrodynamic radii (Rh), determined from experimental diffusion coefficients (from DOSY experiments), were based on a model that assumes that SCNPs are hard spheres in solution, giving access to the ‘apparent’ sizes of the SCNPs without information about the absolute size of SCNPs.69 In this particular context, radii of gyration and shape descriptors determined from sufficiently long MD simulations gave access to a counterpoint to hydrodynamic radii in terms of sizes and shape of heptamers. For the non-complexed heptamers (black dots), the average value of Rg of the S1,4 sequence was low compared to those found for the S1,2 and S1,7 sequences (see Fig. 5a). However, while the average values of Rg of S1,2 and S1,4 were not influenced much by the complexation, the S1,7 sequence had a lower average value of Rg in its complexed form. Therefore, the simulations suggested that the S1,7 sequence was the most compacted upon complexation. These observations were consistent with the folding effect represented by the 2D distance maps (Fig. 4 for S1,7 and ESI Fig. S43, S44† for S1,2 and S1,4, respectively). However, estimates of Rg were averaged over the entire set of MD conformations (see values and standard deviations in the ESI†). The deviation from the ideal spherical shape was not included in such a measurement. Thus, the main challenge was to characterize the shape of heptamers in non-complexed and complexed forms more precisely.
Consequently, the shape characteristics of the heptamers were analysed.73–75 They were described by the ratios λz/λx and λy/λx of principal moments of the gyration tensor, referred to as shape descriptors (see details in the ESI†). The evolution of shape descriptors with the number of spacer units is shown in Fig. 5b. For non-complexed heptamers, all the ratios of eigenvalues were comparable, regardless of the spacer length, which indicated that the heptamers were elongated in one direction and the number of spacer units did not influence the degree of elongation. On the other hand, for the compacted heptamers, the ratios of the eigenvalues converged to a similar value when increasing the spacer length. Thus, the results suggested that, upon complexation, the heptamers with an increasing number of spacer units tended to a spherical shape, which was consistent with the last snapshots of MD simulations (see for example Fig. 4 left for S1,7 and Fig. S43, Fig. S44† for S1,2 and S1,4, respectively).
The asphericity parameter b has previously been used to describe the shape of both folded and unfolded polymeric systems,76 and thus the average values for each sequence were calculated. Since b is a linear combination of the eigenvalues of the gyration tensor (see details in the ESI†), a similar behaviour as in λz/λx and λy/λx was observed in the average value of the asphericity parameter b as a function of the spacer length (Fig. 5c). Indeed, for the non-complexed heptamers, the values of the asphericity parameter b were large, indicating that the shapes of heptamers were far from a spherical symmetry. In contrast, for the complexed heptamers, the asphericity b decreased as the spacer length increases, again pointing out that the structures evolved from elongated (for close functional handles) towards more compact conformations when the number of spacer units was increased.
The overall change in the solvent-accessible surface area of heptamers upon complexation was an additional effect of local conformational rearrangements. On an average, the heptamers gained around 1230 Å2 of solvent accessibility upon complexation (Table S15†). This could be related to a more structured overall arrangement of the heptamers upon compaction, leading the structures to get more exposed to the surrounding methanol molecules. However, the accessibility of complexed heptamers was not influenced by the spacer length. Although the complexed structures adopted more compact conformations as the spacer length increased, they had similar percentages of their surface areas that were exposed to the surrounding solvent indicating that the most compact structure (i.e. S1,7) had likely a more complex surface roughness compared to those of the more extended structures (i.e. S1,2 and S1,4).
To investigate the effect of the spacer length on their folding behaviour, we further performed MD simulations for the dodecamer (S1,12) and the heptadecamer (S1,17), containing two functional handles separated by 10 and 15 spacer units, respectively. MD simulations show the same qualitative trends as those observed for complexation experiments. Average values of radii of gyration (Fig. 5a) suggested that upon complexation the compaction increases with an increasing number of spacer units in the sequence. Indeed, the 2D distance maps indicated that the conformational changes of oligomers upon complexation are more pronounced when increasing the length of the spacer (see Fig. 4c). As for the heptamers, complexation leads to a more structured overall arrangement at the picolyl functional handles for the longer oligomers (see ESI Fig. S43 and 44†). Furthermore, analysis of the shape descriptors and the asphericity parameter b revealed that the morphologies of SCNPs of S1,12 and S1,17 (Fig. 5b and c) are closer to spherical objects (lower asphericity) than the non-complexed oligomers. This effect is more pronounced as the spacer length increases. However, though some local compact globular shapes are formed in the 3D structure, SCNPs still exhibit extended portions (loop-like structures), which lead to non-spherical morphologies.
Footnote |
† Electronic supplementary information (ESI) available. See DOI: 10.1039/d1py00229e |
This journal is © The Royal Society of Chemistry 2021 |