Jinming
Sun‡
a,
Zi-You
Tian‡
a,
Jianbo
Liu‡
b,
Chuan
Wan
d,
Chuan
Dai
b,
Zhihong
Liu
b,
Yun
Xing
a,
Yujie
Wu
e,
Zhanfeng
Hou
a,
Wei
Han
c,
Feng
Yin
*ab,
Yuxin
Ye
*b and
Zigang
Li
*ab
aState Key Laboratory of Chemical Oncogenomics, School of Chemical Biology and Biotechnology, Peking University Shenzhen Graduate School, Shenzhen 518055, China. E-mail: lizg@pkusz.edu.cn
bPingshan Translational Medicine Center, Shenzhen Bay Laboratory, Shenzhen 518118, China
cDepartment of Chemistry, Faculty of Science, Hong Kong Baptist University, Kowloon Town, Hong Kong SAR, China
dCollege of Health Science and Environmental Engineering, Shenzhen Technology University, Shenzhen 518118, China
eOffice of Core Facilities, Shenzhen Bay Laboratory, Shenzhen 518118, China
First published on 9th August 2024
In nature, biochemical processes depend on polymorphism, a phenomenon by which discrete biomolecules can adopt specific conformations based on their environment. However, it is often difficult to explore the generation mechanism and achieve polymorphic control in artificial supramolecular assembly systems. In this work, we propose a feasible thought for exploring the transformation mechanism of polymorphism in peptide assembly from the perspective of thermodynamic regulation, which enables polymorphic composition to be limited by switchable intramolecular CH⋯π attraction within a certain temperature range. Combined with the density functional theory calculations, we obtained thermodynamic theoretical data supporting the conformation transition and the underlying polymorphism formation principle. Afterward, we properly designed the peptide to alter the probability of CH⋯π attraction occurring. Then, we selectively obtained a homogeneous morphological form with corresponding molecular conformation, which further demonstrated the important role of molecular conformational manipulation in polymorphism selection. This unique template-based strategy developed in this study may provide scientists with an additional line of thought to guide assembly paths in other polymorphic systems.
A great deal of progress has been made in supramolecular assembly in recent decades, particularly in the ability to arrange small molecules into sophisticated architectures.11–16 In nature, self-assembled systems are built using a variety of biomolecules, such as proteins, peptides, and nucleic acids, which usually show morphological and structural–functional diversity.17–22 This has aroused the aspiration of chemists to extend the supramolecular assembly area from specific small-molecular to complicated biomacromolecule systems in the laboratory. For example, mechanistic studies of glucose isomerase suggested that polymorph selection occurs at the earliest stage of the nucleation and exerts control over the system through site-directed mutagenesis.23 Various protein polymorphisms have been constructed by simply altering the chemical structure of the ligand.24 Moreover, the assembly of beta-amyloid peptides, closely related to neurodegenerative diseases, has been extensively explored.25–28
In essence, the polymorphism of biomacromolecules and peptides shares the same principle with self-assembly in nature which is brought about by the precisely located, elaborately balanced, and kT-related weak interaction among the subunits.29–32 Generally, the weak interactions contain hydrogen bonding, π–π stacking, and Coulomb attraction. Induced by these generic interactions, the organization of native proteins or peptides into fibers, micelles, or vesicles is obligatorily subject to analog influences.33 Exploring the driving force behind polymorphism, whether from detailed interactions or thermodynamic and kinetic factors, can provide valuable insights into the allosteric mechanisms of macromolecules.15 However, constrained by the relatively large molecular dimension and disordered chemical heterogeneity, especially the lack of stable secondary conformations, protein-directed generation of polymorphic or monomorphic assemblies remains a great challenge. In addition, the complex interactions are essential for addressing complicated cases, which have been difficult to predict with current theory. Furthermore, polymorphic biomacromolecules are influenced by thermodynamic forces.34 However, the distinct kinetic route and conditions often prefer a local rather than a global minimum, complicating the biomacromolecule as an ideal research subject.35 In contrast, the short peptide offers structural simplicity and robust self-assembly features which prove to be the ideal toolbox and building block to explore the insights of complex protein assemblies in nature.36 Specifically, short peptides with a stable secondary conformation could provide a unitary interaction model to dig out the picture of polymorphism. Benefiting from these advantages, short peptides can facilitate conformational/structural analysis, which has attracted widespread attention for their assembly and medical applications in recent years.37–39
Previously, we developed a chirality-induced helical (CIH) cyclic-pentapeptide system with superior assembly capacity. The types of amino acids in the side chain could manipulate the helicity of the short peptide, resulting in the enhanced self-assembly propensity for the CIH peptides.40,41 Benefiting from the discovery, it was realized that the CIH model could be capable of being finely tuned to explore the diversity of conformation with manipulation of amino acids. In this work, we reported an elaborate conformational polymorphism example based on the CIH peptide system. The peptide molecules could self-assemble into two coexisting and stable peptide crystals under appropriate conditions. The X-ray diffraction data and the DFT calculations revealed that a temperature-controlled intramolecular C–H⋯π interaction plays a pivotal role in switching the thermodynamic conformation of the peptide molecule, thereby adjusting its polymorphism. Moreover, through analysis of the structures and design of the sequence, we achieved the singular stabilization of a specific polymorphic form. In this study, we examined the CIH peptide polymorphic transformation process, affecting variables, and thermodynamically feasible theoretical basis, providing a conformation control model to understand the mechanisms of conformational polymorphism at the atom level.
Initially, a set of CIH peptides was employed to investigate self-assembling behavior by the ultrasonic method as previously reported (Fig. S2a†). However, peptides containing hydrophobic amino acids exhibited poor solubility and failed to undergo assembly at various concentrations. Subsequently, an alternative heating-annealing method was explored. Specifically, the peptide liquids were subjected to heating at 363 K for varying durations to facilitate the dissolution of the peptide suspension, and the annealing rate was kept constant at 0.1 K min−1. In this way, we obtained a two polymorph self-assembly product after heating for 240 minutes by using CIH peptide Ac-CAAIaS5(Ph)-NH2 (termed CIHAAIa), in which the second, third, and fourth residues are alanine (A), alanine (A), and, allo-isoleucine (Ia), respectively (Fig. 1a, b, and S2b†). Optical microscopy demonstrated that two polymorphs were generated concomitantly but had different morphologies (Fig. 1c). Scanning electron microscopy (SEM) further elucidated that the growth habits of these polymorphs were rods (termed polymorph-α) and blocks (termed polymorph-β), respectively, both exceeding 20 μm in size (Fig. 1d and e).
The crystals of two polymorphs were successfully separated manually for single-crystal X-ray diffraction (SCXRD) experiments according to their different morphologies. The crystal of polymorph-α was found to be orthorhombic with the space group of P212121, and crystals of polymorph-β were tetragonal with the space group of P41212 (Table S1†). Investigations of the crystal structure revealed that CIHAAIa monomers in two polymorphs were both in helical conformation maintained by four sets of hydrogen bonds (Fig. 2b), which were consistent with the CD analysis result (Fig. S1†). Further intermolecular interaction analysis showed that, in conformer-1 (the monomer conformation of CIHAAIa in polymorph-α), a C–H⋯π hydrogen bond was formed inside the benzene and side chain of the fourth residue (the distance from the γ hydrogen atom of the fourth residue to the benzene ring is 2.88 Å) (Fig. 2b, left). However, in conformer-2 (the monomer conformation of CIHAAIa in polymorph-β), the closest distance between the hydrogen atom on the fourth residue and the benzene ring was 3.46 Å (Fig. 2b, right), which exceeds the maximum distance of C–H⋯π contacts.43 Such an intramolecular interaction difference leads to distinct conformations of the side chain of the fourth residue in two conformers (Fig. 2c). Additionally, the C–H⋯π contact within the conformer-1 resulted in a slight rotation as well as a significant displacement of the benzene ring towards the fourth residue, consequently constraining the peptide in a more compact conformation than that in conformer-2 (Fig. 2c).
The differential monomer conformations, attributed to the presence or absence of C–H⋯π interaction, give rise to distinct packing arrangements of CIHAAIa molecules in the two polymorphs. In polymorph-α, the molecules assemble into one-dimensional cuboid columns along the a-axis, facilitated by the zipper-like hydrophobic interactions formed by the linker region, benzene group, acetyl group, and the side chains of the second and fourth residues (Fig. 3a). Subsequently, these cuboid columns utilize their amphiphilic surfaces to establish dry- and wet-interfaces with adjacent columns (Fig. 3b and S3†). The dry interface at the column's corners is relatively small and stems from the hydrophobic interactions along the b-axis, while the wet interface constitutes the major portion of the column interface, formed by the solvent-related hydrogen bond network (Fig. S3†). Finally, the dry- and wet-layers stacked along the c-axis to form rod-shaped crystals (Fig. 3c). In polymorph-α, the CIHAAIa molecule diverges from the conventional reliance on hydrogen bonds typically found in contiguous helical structures for its assembly. Instead, it assumes a role as a fundamental building block with a distinctive geometry. The growth of rod-shaped crystals is primarily facilitated by the hydrophobicity effect originating from the linker region and hydrophobic side chains along the a-axis (Fig. 3a and c).
In contrast to polymorph-α, the CIHAAIa molecule within polymorph-β adopts an assembly rule that optimally exploits its inherent helical conformation. In polymorph-β, the peptide molecules arrange themselves into a continuous straight α-helical column with head-to-tail interactions along the a- and b-axes, involving the acetyl groups and the carboxylic acid groups (Fig. 3d). The contiguous antiparallel helical columns further stack into a 2D lamellar structure through the zipper-type hydrophobic interaction, mediated by the linker region, benzene group, and hydrophobic side chains (Fig. 3e). This assembly is additionally reinforced by intramolecular C–H⋯π interactions between the benzene ring and the neighboring molecules (Fig. S4†). Subsequently, the adjacent 2D lamellas undergo 90-degree rotation stacking along the c-axis through hydrogen bonds facilitated by a hydrophilic interface and a solvent water molecule, ultimately resulting in the formation of block-shaped crystalline structures (Fig. 3f, g, and S4†).
The success in resolving distinct structures for CIHAAIa in two different morphologies encouraged us to explore the control of polymorphism by conformational transition. Temperature-dependent CD experiments of CIHAAIa revealed a conformational perturbation upon heating in an aqueous solution; the negative peak decreased significantly at 203 nm and 226 nm, and the positive peak increased significantly at 190 nm (Fig. 4a). This implied that the annealing temperature could potentially serve as a parameter for regulating conformational states, thereby influencing polymorphism. To confirm this hypothesis, we conducted a series of self-assembly experiments at varying heating temperatures, while maintaining a consistent annealing rate. The heating temperature was incrementally increased at intervals of 10 K, ranging from 333 K to 373 K. After each temperature interval, samples were collected for detailed characterization through SEM and Powder X-ray Diffraction (PXRD) analysis. The SEM results revealed the exclusive presence of rod-shaped crystals at 333 K and 343 K (Fig. 4b). Correspondingly, the PXRD results exhibited five distinct diffraction peaks at 10.8°, 14.2°, 18.9°, 21.7°, and 23.7°, aligning with the crystal faces (0, 1, 1), (1, 2, 1), (1, 3, 1), (2, 0, 3), and (2, 0, 4) in the simulated PXRD pattern of polymorph-α (Fig. 4c). Conversely, at 373 K, only tetragonal block-shaped crystals were observed (Fig. 4b), featuring three distinct diffraction peaks at 11.0°, 16.8°, and 22.4° in the PXRD results. These peaks corresponded to the crystal faces (0, 0, 8), (0, 0, 12), and (0, 0, 16) in the simulated PXRD pattern of polymorph-β (Fig. 4c). At 353 K and 363 K, both rod-shaped and block-shaped morphologies were observed (Fig. 4b), and all representative PXRD diffraction peaks characteristic of both polymorph-α and polymorph-β were detectable (Fig. 4c). These results indicated that CIHAAIa undergoes distinct stages of polymorph formation at specific temperatures. Taken together, this outcome underscores that the supramolecular assembly of CIHAAIa is a temperature-controlled process.
After achieving precise temperature control over polymorphism, we next contemplated whether such polymorphism regulation was due to the effect of temperature on the conformation state of the peptide monomer as we expected. Therefore, we executed Density Functional Theory (DFT) calculations of CIHAAIa in terms of conformational transformation. Two calculated conformations, iL1, and iL2, were located (Fig. 5a), and both closely resembling the experimental conformations, conformer-1 and conformer-2, respectively (the detailed calculations and benchmark discussion are provided in Section 4.1 of the ESI†). Initially, with reference to iL1, the Gibbs free energy difference at room temperature for iL2 in the aqueous phase was 0.4 kcal mol−1 (Fig. 5b), indicating that the conformational transformation from iL1 to iL2 is thermodynamically prohibited as evidenced by the positive ΔG value. This finding aligns with the experimental result that the preferred polymorph is determined by conformer-1 at lower temperatures, which is polymorph-α. Furthermore, AIM analysis revealed the presence of CH⋯π attraction44–47 in iL1 (Fig. S5†), while no such interaction was observed in iL2, suggesting that the thermodynamically favored conformation relies on this specific CH⋯π interaction.
Notably, the activation Gibbs free energy difference associated with the transition state (TS) (Fig. 5a) for the conformational shift from iL1 to iL2 was determined to be merely 6.4 kcal mol−1 (Fig. 5a). This low kinetic barrier implies that conformational changes between iL1 and iL2 can readily occur once the reaction is thermodynamically permissible.48,49 Subsequent thermodynamic calculations, incorporating temperature correction, demonstrated that the energy gap between these two conformations gradually diminishes as the temperature increases (Fig. S6†). Beyond the critical temperature of approximately 401 K, where ΔG equals 0, the conformational transformation becomes thermodynamically permissible, resulting in the equilibrium coexistence of both conformations (Fig. 5b). This is consistent with experimental observations of the simultaneous appearance of the two polymorphic forms. Of note, the error between the calculated critical temperature and the practical temperature was accepted, which was attributed to the precision related to the DFT calculations that was less than 1 kcal mol−1. Finally, the sign of ΔG switched as the temperature continued to increase,50 consistent with the dominance of the conformer 2 and related polymorph-β at elevated temperatures (Fig. 5b).
The DFT calculations of CIHAAIa revealed that the CH⋯π interaction, responsive to temperature variations, could thermodynamically influence the peptide conformation. This underscores the significance of CH⋯π interaction in controlling polycrystalline forms. Of note, such CH⋯π interaction emanated from the Cγ1 atom of the fourth allo-isoleucine side chain, characterized by two C–H bonds and one C–C bond (Cγ1–Cδ1) (Fig. 6a). We hypothesized that altering the number of available C–H bonds on this Cγ1 atom could manipulate the possibility of CH⋯π interaction, consequently altering ΔG to regulate the conformation distribution and attain control over polymorphic assembly. Two distinct peptide design pathways arise from this hypothesis: (i) mutating the fourth allo-isoleucine to valine to increase the number of C–H bonds on the Cγ1 atom to three, thereby increasing the possibility of CH⋯π interaction; (ii) reducing the C–H bond on the Cγ atom proximal to the benzene ring to just one by the substitution of the fourth allo-isoleucine with the more hindered leucine (Fig. 6a).
To validate our hypothesis, we synthesized CIHAAV and CIHAAL and conducted self-assembly experiments by the annealing method. SEM investigations revealed the exclusive formation of rod-shaped crystals for CIHAAV and block-shaped crystals for CIHAAL, respectively (Fig. 6b). SCXRD results indicated that both CIH peptides molecules exhibited helical conformations (Fig. 6c), consistent with CD analysis (Fig. S1†). Additionally, as anticipated, a 2.77 Å CH⋯π interaction was observed between the fourth valine side chain and benzene ring in CIHAAV, but no CH⋯π interaction was observed in CIHAAL (Fig. 6c). The single molecule X-ray structures of CIHAAV and CIHAAL were identical to the conformer-1 and conformer-2 of CIHAAIa, respectively (Fig. S7†). Furthermore, the space group, cell parameters, and molecular arrangements of the assemblies of CIHAAV and CIHAAL were consistent with polymorphs-α and -β of CIHAAIa, respectively (Table S2 and Fig. S8†). These results indicate that the modulation of CH⋯π interaction through amino acid substitution allows for the adjustment of CIH peptide conformation, consequently achieving uniformity in polymorphic assembly.
The temperature-dependent CD analysis revealed that the conformation of CIHAAV and CIHAAL remained unaltered with temperature variations (Fig. S9†). Moreover, the self-assembly form exhibited resilience to temperature changes, thereby indicating the inherent conformational stability of CIHAAV and CIHAAL. Therefore, the DFT calculations were employed to elucidate the thermodynamic stability of the conformations for CIHAAV and CIHAAL. For CIHAAV, three different conformations were located, which was definitely distinct from the isoleucine scenario. The preferred conformation V1 possessed an analogous CH⋯π interaction concerning the X-ray structure of CIHAAV (Fig. 6d). In addition, another two inferior conformations, which had no CH⋯π interaction, were less favored. In conformation V2, three severe H⋯H steric hindrances were induced within the valine side chain, the backbone of peptide, and the phenyl motif (Fig. 6d). In conformation V3, similar three serious H⋯H steric hindrances dominate the less favored consequence (Fig. 6d). These unfavorable elements in V1 and V2 contributed the higher Gibbs free energy with 1.3 and 2.0 kcal mol−1, underscoring that the reversal of the sign of ΔG with the less favored conformation demands higher temperature and less feasible under atmospheric conditions (Fig. 6d). Similarly, three distinct conformations were located for CIHAAL. As predicted, the preferred conformation L1 harbored less steric hindrance and no CH⋯π interaction was observed (Fig. 6e). The less favored conformation L2 was clearly understood that a dominant H⋯H repulsion was offered with 2.26 Å length, providing the higher 4.3 kcal mol−1 with Gibbs free energy (Fig. 6e). However, in another less favored conformation L3, besides an H⋯H steric hindrance with 2.3 Å was observed, a 2.8 Å CH⋯π interaction was also observed based on the DFT calculation (Fig. 6e). Nonetheless, the gap of the Gibbs free energy in terms of the preferred one was compensated to 2.6 kcal mol−1. Of note, the calculation related to the CIHAAL case was subtly different from the previous one (the detailed discussion is provided in Section 4.2 of the ESI†). Together, our DFT calculations are capable of providing molecular insight to reveal the essence of manipulating the polymorphism.
Footnotes |
† Electronic supplementary information (ESI) available. CCDC 2292337, 2292338, 2292339, and 2292343. For ESI and crystallographic data in CIF or other electronic format see DOI: https://doi.org/10.1039/d4sc02545h |
‡ J. S., Z. T., and J. L. contributed equally to this work. |
This journal is © The Royal Society of Chemistry 2024 |