Chiral expression from molecular to macroscopic level via pH modulation in terbium coordination polymers

Chiral expression from the molecular to macroscopic level is common in biological systems, but is difficult to realise for coordination polymers (CPs). The assembly of homochiral CPs in both crystalline and helical forms can provide a bridge for understanding the relationship between the molecular and macroscopic scales of chirality. Herein, we report homochiral helices of [Tb(R- or S-pempH)3]∙2H2O (R - or S -1) (pempH2 = (1-phenylethylamino)methylphosphonic acid) and their crystalline counterparts (R - or S -3), which are formed at different pH of the reaction mixtures under hydrothermal conditions. By combining the experiments and molecular simulations, we propose that the formation of helices of R -1 or S -1 occurs via a hierarchical self-assembly route, which involves twisted packing due to the geometric incompatibility of the different types of chains. The observed chiral transcription from molecules to morphologies is significant for understanding bio-related self-assembly processes on the nano- to macro-scale.

1 I n biological systems, chirality is widely expressed from small molecules such as L-amino acids and D-sugars to macroscopic helices such as double helical DNA and α-helical proteins. The helical assembly of large biomolecules not only enables their biological functions required to sustain life, but also allows for the development of novel applications such as bio-computing, biocatalysis, biomedicine and material science [1][2][3][4][5][6][7][8] . On the other hand, DNA and protein crystallisation is increasingly demanding in order to elucidate the structure-function relationships at the atomic and molecular scale in the biological systems. Despite the tremendous difficulties in assembling the helical bio-motifs into three-dimensional crystalline lattices, significant progress has been witnessed in obtaining diffraction-quality protein, peptide and DNA crystals and in establishing their roles as catalytic materials, drug delivery vehicles and chemical reaction vessels [9][10][11][12][13] .
DNA and protein crystallisation requires the enhancement of interactions between helical motifs through non-covalent or covalent bondings. One of the most promising routes is metalmediated head-to-tail assembly by taking the advantage of the metal-ligand interactions at the ends of the growing crystal 14 . In this case, the resulted crystals can be viewed as 'biocoordination polymers' in which the helical biomacromolecules are connected by metal ions into polymers with ordered structures.
In non-biological systems, coordination polymers (CPs) refer to a broad range of metal-organic compounds with one-, twoand three-dimensional structures, composed of metal ions as nodes and small organic ligands as linkages. These materials usually possess precise structures that can be designed and constructed through careful selection of metal ions and ligands 15,16 . As a subclass of CPs, homochiral CPs are very attractive in the fields of enantioselective catalysis and separation, chiral sensing, nonlinear optics and multifunctional materials [17][18][19][20] . In contrast to the biomacromolecules that are difficult to crystallise, most homochiral CPs appear as crystalline materials with well-defined structures, and the chirality is typically expressed only at the molecular level. Examples of homochiral CPs with helical morphologies are extremely rare [21][22][23][24][25][26] , despite the great efforts devoted to the fabrication of artificial chiral nanoarchitectures in order to mimic the natural biological systems [27][28][29][30][31][32][33] . By utilising chiral amino acids, Tang et al. obtained homochiral Ag(I)/ cysteine helical nanobelts in which the chirality transcription occurs from cysteine molecule to the assembly entities 21 . The proposed mechanism involved the merging of nonhelical nanobelts of Ag(I)/cysteine layers through lateral attachment, which developed into hierarchical helices with a specific twist direction. Huang et al. isolated right-handed Ca-cholate helical nanoribbons, which were further used as templates for the fabrication of helical inorganic nanomaterials 22 . The formation mechanism was supposed to involve the twist of supramolecular layers composed of cholate bilayer strips connected via calcium-carboxyl coordination. In all these cases, crystalline forms of the same materials were not obtained. Therefore, structural illustration of these helices was not very clear due to the absence of single crystal data, although theoretical calculations were conducted to propose the structures. The scarce examples demonstrate the challenges remaining in the construction of homochiral helical CPs. More challenging is the assembly of homochiral CPs in both crystalline and helical forms, which can provide a bridge for understanding the relationship between the molecular and macroscopic scales of chirality.
In nature, the macroscopic scale biomolecules contain supramolecular nanostructures, hierarchically organised via weak interactions such as π-π stacking, hydrogen bonding and hydrophobic interactions. Hence, the selection of suitable organic ligands and the control of weak interactions between the chain or layer motifs could be a key for the construction of homochiral CPs in crystalline and/or helical morphologies.
The optically active amino phosphonic acids are known to play important roles in biological activities 34,35 . They can dissociate in aqueous solutions with the release of one or two protons, depending on the pH of the media. Each phosphonate group can bind up to nine metal atoms with versatile coordination modes 36,37 . Modification of the organic groups can provide additional sites for metal coordination, hydrogen bonding and/or π-π stacking. Previously, we found that enantiopure R-or S-(1phenylethylamino)methylphosphonic acid (R-or S-pempH 2 ) can react with metal ions to form layered 38,39 or nanotubular 40,41 structures. Considering that trivalent lanthanide ions show high kinetic lability and a lack of stereochemical preference, herein we report on the reactions of Tb 3+ nitrate and R-or S-pempH 2 under hydrothermal conditions. Interestingly, homochiral coordination polymers of [Tb(R-or S-pempH) 3 ]•2H 2 O with both helical (R-, S-1) and crystalline (R-, S-3) morphologies are obtained by pH modulation of the reaction mixture. Results based on both the experiments and molecular simulations reveal that the formation of helices of R-1 or S-1 occurs via a hierarchical self-assembly route, which involves twisted packing due to the geometric incompatibility of the different types of chains.

Results
Synthesis and characterisation of homochiral helices. An aqueous solution of Tb(NO 3 ) 3 and R-pempH 2 ( Fig. 1) in molar ratio 1:5 was adjusted to pH 3.1 using 0.5 M NaOH and heated hydrothermally at 120°C for 2 days. After cooling to room temperature, both the flocculent precipitates on the top and the powder precipitates on the bottom of the vessel were obtained. The flocculent precipitates were collected manually, and the other precipitates were collected by suction filtration. Powder X-ray diffraction (PXRD) measurements confirmed that the two precipitates were the same material, which is referred to hereafter as R-1. When S-pempH 2 was used as the starting material, the final product, referred to as S-1, showed an identical PXRD pattern to that of R-1 ( Supplementary Fig. 1). Figure 2 shows the scanning electron microscope (SEM) images of the as-synthesised flocculent precipitates of R-1 and S-1. R-1 exhibited a pure right-handed helical morphology, whereas S−1 showed a pure left-handed helical morphology. The diameters of the helices were ca. 6-9 μm and the lengths ranged from 100 to 300 μm. The pitches of the helices were~20-30 μm, and the pitch angles ranged from 70 to 80°for most helices ( Supplementary Fig. 2).
The mirror image relationship of R-1 and S-1 was reflected not only in their morphologies, but also in their optical properties. The infrared spectra (IR) of both R-1 and S-1 revealed sharp peaks in the 900-1200 cm −1 region, which were attributed to the stretching vibrations of the -PO 3 group ( Supplementary Fig. 3). Solid-state circular dichroism (CD) spectra revealed Cotton effects centred at 262 nm, originating from the ligand; opposing symmetries were observed for R-1 and S-1, which indicated that the chirality was transferred from the ligand to the inner structure of the helices ( Fig. 2e Table 2) and EDX analyses, the molecular formulae of the homochiral helices were proposed to be h-Tb(R-pempH) 3 •2H 2 O for R-1 and h-Tb(S-pempH) 3 •2H 2 O for S-1 (h means helical). The number of water molecules was confirmed by thermal analyses ( Supplementary  Fig. 9).
Influence of reaction conditions on the helix formation. To determine the effects of the reaction conditions on the formation of the helices, a systematic study was performed on the Tb/R-pempH 2 system by varying the pH, reaction time, temperature, as well as the cations and anions in the reaction mixture. First, Tb(NO 3 ) 3 was allowed to react with R-pempH 2 under hydrothermal conditions (120°C, 2 d) at different pH values, which were appropriately adjusted using 0.5 M NaOH (Supplementary Fig. 10). A transparent solution without any precipitates was obtained when the pH was sufficiently low (1.5-2.4). At pH 2.5-2.7, block-like colourless crystals of (H 3 O) [Tb 3 (R-pempH 2 ) 2 (R-pempH) 7 ][Tb 3 (R-pempH 2 )(R-pempH) 8 ] (NO 3 ) 4 •11H 2 O (referred to hereafter as R-2) formed (Supplementary Fig. 11), together with a small amount of white precipitate of R-1. Higher pH value (2.7-2.9) resulted in a mixture of block-like crystals of R-2 and helices of R-1. When the pH was 3.0-3.2, only helices of R-1 were obtained. When the pH was further increased to 3.3-3.6, both helices of R-1 and rod-like crystals of c-Tb(R-pempH) 3 •2H 2 O (referred to hereafter as R-3, c means crystalline) were isolated. Finally, at pH 3.7-4.5, only R-3 was isolated. These results clearly indicated that the pH of the reaction mixture was essential for the formation of the helices. To investigate whether helices of R-1 could be crystallised, we added additional 0.5 M NaOH solution to the helices in mother liquid (obtained at pH = 3.1) until the final pH was about 4.5, and put the autoclave back into the oven at 120°C. After 2 days, the white powdered helices were transformed into the rod-like crystals of R-3. Therefore, we propose that the structure of the helices of R-1 is closely related to the structures of R-2 and R-3.
Next, hydrothermal reactions of Tb(NO 3 ) 3 and R-pempH 2 (pH 3.1) were conducted at 120°C for different periods of time (1-24 h). As shown in Supplementary Figs. 12 and 13, helices of R-1 were present even at 1 h. When the reaction time was increased, both the length and the diameter of the helices increased. When the reaction time was sufficiently long, e.g., after 8 h, the morphology of the final products did not change significantly.
We next asked whether the presence of metal ions other than Na + could affect the final products. To test this, we conducted the hydrothermal reactions of Tb(NO 3 ) 3 and R-pempH 2 (pH = 3.1) at 120°C for 2 d, and used other alkali or alkaline earth metal hydroxides, such as KOH and Ca(OH) 2 , to adjust the pH (instead of NaOH) or added Ba(NO 3 ) 2 /Sr(NO 3 ) 2 to the reaction mixture. The SEM images and the PXRD patterns of the final products showed that generally the cations had little influence on the formation of the helices (Supplementary Figs. 14 and 15). However, the addition of a particular cation did affect the morphology of the helices. Supplementary Fig. 14 clearly shows that helices formed in the presence of Ba 2+ had smaller helical pitches.
The effect of anions was also studied. Different terbium salts, such as TbCl 3 , Tb(OAc) 3 and Tb(ClO 4 ) 3 , were used as representative to react with R-pempH 2 (pH = 3.1, adjusted using 0.5 M NaOH) at 120°C for 2 d. Surprisingly, no helices were found in the final products, as shown in Supplementary Fig. 16. In other words, among the tested anions, only NO 3 − could induce the formation of helices.
Finally, hydrothermal reactions of Tb(NO 3 ) 3 and R-pempH 2 (pH 3.1) were performed for 20 h but at different temperatures (80-160°C). Supplementary Fig. 17 showed that helices of R-1 were obtained when the reaction temperature reached 100-140°C, but the pitches of the helices were significantly different.
Apparently, helices of R-1 can form in a suitable temperature range. The reaction time and addition of alkali or alkaline earth metal ions could affect the size and morphology of the helices. The helix formation is pH-sensitive and anion-dependent. Notably, helices formed at ca. pH 3.1, whereas block-or rodlike crystals of R-2 and R-3 were isolated below or above this pH. An accurate determination of the crystal structures can help reveal information about the formation mechanism of the helices. We next examined the crystal structures of the compounds we obtained.
The N and C atoms are omitted for clarity. b Side and top views of the chains in structures R-2, R-2′ and R-3, which contain left-handed triple strands. The pitches and diameters of each helical chain are also given. c Packing diagrams of structures R-2, R-2′ and R-3 viewed along the crystallographic bor c-axis. The interval of adjacent chains are defined by the distances between the relative 2 1 axis along the b-direction in R-2 and R-2′ or 6 5 axis along the c-direction in R-3 crystallographically distinguished NO 3 − anions (N19, N20, N21 and N22) and one H 3 O + (O11W) exist as charge-balancing counterions between the chains. The H 3 O + cations are hydrogen bonded to the phosphonate oxygen atoms from chain-II and lattice water molecules (Fig. 4). The NO 3 − anions are involved in the hydrogen bond network with the protonated phosphonate oxygen (O18, O27 and O30), amino nitrogen, -CH 2and phenyl groups as well as lattice water molecules (Supplementary Table 5).
Structural description of homochiral block-like crystals of R-2′. Interestingly, compound R-2 underwent a solid state structural transformation soon after exposure to air at room temperature, forming Tb 3 (R-pempH 2 ) 2 (R-pempH) 7 Table 8). The equivalent Tb atoms are triply bridged by two μ 3 -O(P) and one O-P-O units into a helical chain running along the c-axis (named as chain-III) (Fig. 3a). The Tb···Tb separation is 4.043(1) Å. The Tb···Tb···Tb angle is 174.8 (1)°. Chain-III also contains left-handed triple helical strands (Fig. 3b). Remarkably, although the pitch of the helical chains in R-3 is only slightly smaller than those in R-2 (24.16 Å vs. 24.28 Å in R-2), there is a significant reduction in the chain diameter (16.71 Å vs. 18.11-18.98 Å in R-2). The adjacent helical chains are packed into a 3D supramolecular framework via hydrogen bond interactions (Supplementary Table 9) with an interchain distance of 15.824 Å (Fig. 3c).
Comparison of the rod-like crystals of R-3 and helices of R-1. The chemical composition of R-3 is identical to that of R-1, suggesting that the two materials could be closely related to each other. Indeed, the IR spectra of R-1 and R-3 were identical, and the PXRD patterns of R-1 and R-3 were similar, as shown in Fig. 5. The reflections of (h00) for R-1 appeared as a series of strong and equidistant peaks with h up to 5 due to preferential orientation, which indicated that the helices possessed long-range intermolecular order and had a similar structure as that of the rod-like crystals of R-3. The diffraction peaks can be indexed by using TOPAS 4.2 programme 42 , giving a set of unit cell parameters with space group P6 5 , a = 16.28 Å, c = 24.24 Å and V = 5564.9 Å 3 for R-1 ( Supplementary Fig. 30). Compared with those in R-3, all peaks for R-1 helices were shifted to lower angles and the cell volume was expanded (Supplementary Table 10). As such, a larger interchain distance can exist in the helices of R-1 than in the crystals of R-3. The average distance between two layers assembled by the 1D helical chains in the bc plane was estimated to be 14.1 Å in the helices of R-1, much larger than that in the crystals of R-3 (13.7 Å). This difference may originate from the distortion of the 1D chains during the self-assembly process. According to the experimental data, we conclude that the helices of R-1 have a high degree of crystallinity, with almost the same structure at the molecular level as that of rod-like crystals of R-3, but at the macroscopic level different from R-3. To the best of our knowledge, this is the first example of coordination polymers that show both crystalline and helical morphologies.
Characterisation of intermediate self-assembly states. In order to get more insights into the formation mechanism of the chiral morphology, we monitored the solid products after hydrothermal reactions of Tb(NO 3 ) 2 and R-pempH 2 (pH = 3.1) at 120°C for different period of time. The PXRD measurements revealed that   the 0-min product (before hydrothermal reaction) contained only the ligand of R-pempH 2 . After 10 min, helices of R-1 appeared together with the ligand, evidenced by the emergence of a peak at a lower angle of 6.183°. The two phases, e.g., R-1 and the solid R-pempH 2 , coexisted in the reaction mixture even after 1 h. When the reaction time reached 2 h and above, pure phases of R-1 can be observed (Supplementary Fig. 31). The IR spectra can recognise the presence of un-coordinated nitrate anions by the peak at ca. 1385 cm −1 , which was significant in the case of pure R-2. As shown in Supplementary Fig. 32, the intensity of this peak was markedly increased when the reaction time reached 4-6 h, indicating that the NO 3 − anion and hence the positively charged chains (chain-I, -II in R-2) could be involved in the helices of R-1. The increase of the peak intensity cannot be identified when the reaction time was less than 2 h or more than 8 h, possibly due to the interference of the ligand and/or the presence of neutral chains in the helices of R-1. The results suggested that the amount of positively charged chains in R-1 could be very small compared with that of the neutral ones. This may explain the fact that both the PXRD pattern and the chemical composition of the helices of R-1 are close to those of compound R-3.
The SEM images give a clear visualisation about the morphology of the solid products. The 0-min product contained sheets of R-pempH 2 covered by some amorphous nano-particles of size 100-300 nm (Supplementary Fig. 33). The nano-particles without helical morphologies were also observed in the 10-min product (Fig. 6), although PXRD result suggested the helix formation after 10 min of reaction. In the 15-min product, aggregates of both nano-particles and nanorods were found. Still, the nanorods did not show helical morphologies, but their Tb/P ratio was close to 1/3. For comparison, the Tb/P ratio in the nano-particle area was ca. 1/7 (Supplementary Fig. 34).
Interestingly, small helices appeared in the 20-min product like an actinian, the widths and lengths of which were ca. 150-500 nm and 2-5 μm, respectively. The Tb/P atomic ratio was ca. 1/3 for the helices, but ca. 0/3 for the bottom area without helices ( Supplementary Fig. 35). After reacting for 30 min, aggregates of helices were observed like the bird's nest. Compared with the 20min product, however, the lengths of the helices (ca. 2.05-11.78 μm) become longer, and the widths (ca. 0.19-0.97 μm) become wider ( Supplementary Fig. 36). For the 1-h product, starfish-like aggregates of helices can be recognised together with the separate ones. Meanwhile, the lengths (ca. 6.37-119.45 μm) and the widths (ca. 0.19-2.26 μm) of the helices increased further, from nanometre to micrometre scale. After 2-8 h of hydrothermal reactions, the helices of R-1 appeared separately, and the average length and width of the helices increased with increasing reaction time. Notably, the widths of the helices do not change significantly after 8 h of reactions, although the lengths of the helices increased continuously (Fig. 6

and Supplementary Figs. 37 and 38).
Noting that the growth of helices occurred on the surface of the ligand, the un-dissolved R-pempH 2 could play two roles during the helix formation: (1) it served as a buffer against the pH change in solution during the self-assembly process. The pH of a saturated solution of R-pempH 2 was ca. 3.5 ( Supplementary  Fig. 39). The coordination of R-pempH 2 with Tb 3+ would release the protons, and then decreased the pH in solution (Supplementary Fig. 40). Thus, the dissolve of solid R-pempH 2 could help maintaining the pH of solution in a suitable range (ca. 3.1) during helix formation. (2) It served as nucleation centres at which the helices of R-1 grew.
From the above experimental results, we can conclude that the helices of R-1 formed on the surface of the un-dissolved ligand of R-pempH 2 , forming nano-particles first, then the nanorods, and finally the helices of R-1. The growth of the helices followed a hierarchical process with the length direction growing much faster than the width direction. Both the neutral chains of R-3 and the positively charged chains of R-2 were involved in the helix formation process.

Molecular modelling of the formation mechanism of helices.
To better understand the underlying mechanism of the helix formation, we applied the all-atom molecular dynamics simulation to investigate the interaction energy among Tb 3+ , R-pempH − , R-pempH 2 , and also used the coarse-grained Brownian dynamics (BD) simulation to study the self-assembly of chain-II and chain-III under different cases. As shown in Fig. 7, the interaction energy between Tb 3+ and R-pempH − was much larger (about 15 times) than that between R-pempH − and R-pempH − (or R-pempH 2 ), indicating that the growth rate along the chain direction should be faster than that along the side direction. However, in the presence of NO 3 − , since there existed hydrogen bonds among NO 3 − , R-pempH − and R-pempH 2 , the interaction energy increased a lot (about three times). As a result, the difference of the growth rate between the chain direction and side direction would not become very obvious and the length and diameter of assembly would be comparable, which may help explain the block crystal in R-2 system.
To better clarify the experimental observations, we further used Brownian dynamics simulation to investigate the system from the mesoscopic view (Fig. 8). In the case of pure chain-III, rod-like chains with one or several monomers in the diameter were observed (Fig. 9a). While in the case of pure chain-II, due to the strong side interaction, block-like aggregates with similar length and height were observed (Fig. 9b). These results were in accordance with the inference by all-atom simulation, and consisted with the experimental findings. More importantly, in the case of the mixture of chain-III and chain-II (the ratio is about 4:1), the curved or twisted chains were found (Fig. 9c). Since the molecular symmetries of chain-III (P6 5 ) and chain-II (P2 1 ) are different, when chain-II bound to chain-III along the axis, the assembled chain would become a bit curved instead of linear growth in R-3 system. Moreover, since the side interaction between chain-II and chain-II was much larger than that between chain-III and chain-II/chain-III, chain-II preferred to aggregate with each other, which further made the growth along the chain direction more curved. Notably, the curved chain does not mean that it is chiral. But the curvature of the chain must be the prerequisite for the helix formation; in other words, pure linear chains can never form helices. As shown in Supplementary  Fig. 41, the correlation function of the linear chain was totally different from that of the curved chain and the helical chain. Recently Grason et al. 43 also demonstrated theoretically that chiral filaments may occur due to the frustration of inter-filament spacing when there exist many self-twisting curved chains or bundles. We should mention that due to the limitation of present computing technology, the length scale of the obtained chains in the simulations was about tens of (at most one hundred) nanometres, which was far from the length scale of the helices (tens of micrometres) in the experiments. Thus the chains here were still not the direct evidence of macroscopic helices. Nevertheless, the present simulation provides some hints or possible explanations for the reason why the macroscopic scales of chirality could be observed in this coordination polymer system.
Additionally, if there did not exist the NO 3 − ions in the system, namely, the side interaction between chain-II and chain-II was the same as that between chain-III and chain-II/chain-III, the aggregation of chain-II in the assembly would not become obvious, and finally the assembly was nearly linear chain instead of twisted one (Fig. 9d), which indicated the importance of NO 3 − in the chiral assembly and was in good agreement with the experimental findings.
Besides, we also found that if the mixture ratio of chain-III and chain-II was lower (e.g., 1:1), since the side interaction between chain-II and chain-II was large, block-like aggregates were again observed (Fig. 9e). On the contrary, if the mixture ratio of chain-III and chain-II was much greater (16:1), since there existed no enough chains-II to induce the curved growth of chains-III, the twist or curve of the assembly was not obvious (Fig. 9f), indicating that the ratio of chain-III and chain-II was also of great importance in helices formation. Considering that pH is related to the ratio of chain-III and chain-II in the system, our simulation results here could also clarify the role of pH in the experiments.

Discussion
We described the generation of homochiral helices of R-1 or S-1 under hydrothermal conditions when the pH of the reaction mixture was about 3.1. Below or above this pH, crystalline compounds R-2 (or S-2) and R-3 (or S-3) were isolated. Such a pH-responsive morphology change is reminiscent of biological systems, where a pH is often crucial for conformational changes  Chain-III Chain-II Fig. 8 The setup of Brownian dynamics simulations. a Schematic illustration of the CG models for R-pempH − (R-pempH 2 ) molecule and Tb 3+ . Snapshots of CG model for chain-II b and chain-III c with six Tb 3+ ions and eighteen R-pempH − (R-pempH 2 ) molecules. To better differentiate the chain-II and chain-III, the bead that can interact with Tb 3+ is used as pink and lime one, respectively. d The initial conformation of the self-assembly system, where the chains-II (and/or chains-III) are distributed uniformly in the simulation box positively charged chains, which coexist in the same reaction solution under suitable pH condition (ca. 3.1). The NO 3 − ions, which serve as spacers and hydrogen bond acceptors, are responsible for the twist of the chain aggregates, and hence play an important role in the formation of the helices of R-1.
Based on both the experimental and simulation results, the formation mechanism of the helices of R-1 can be proposed as illustrated in Fig. 10. Assembly starts with the formation of Tb/Rpemp complexes, and fragments of 1D helical chains are constructed via coordination interactions between the Tb 3+ ions and phosphonate ligands. The degree of protonation of the phosphonate groups and thus the composition of the helical chains in solution are dependent on the pH of the reaction mixture. At pH ca. 3.1, both neutral (chain-III in R-3) and positively charged chains (chain-I and chain-II in R-2) coexist in the same reaction solution with the former being the dominant species. The binding of the positively charged chains to the neutral ones triggers the twist of chain growth in pure R-3 system. Since chains of R-2 prefer to aggregate with each other, the growth along the chain direction becomes more curved and twisted. The twisted chains further assemble into hierarchical bundles, like the artificial supramolecular helices 45,46 . The assembly rate along the length of the helix is faster than that of the width, thus leading to R-1 with a 1D helical morphology.
The helices of R-1 are right-handed, whereas the Tb(R-pempH) 3 chain (chain-III) has a left-handed helical conformation. Same phenomenon was observed in biomolecules such as collagen, a right-handed helical superhelix composed of three polypeptides where each has a left-handed helical conformation 47 . The chirality of the chains is determined by the stereo configuration of the phosphonate ligand, whereas the transfer of chiral information from molecular to the macroscopic level can be attributed to the weak interchain interactions. Clearly, the pendant amino and phenyl groups of the phosphonate ligands and nitrate ions play non-trivial roles in the successful construction of helices of R-1 and S-1.
In summary, we report homochiral CPs with formulae [Tb(Ror S-pempH) 3 ]•2H 2 O with both helical and crystalline forms. The geometric incompatibility of the coexisting chain types in the same reaction mixture leads to the formation of helices. This work provides new insight into the design and construction of homochiral CPs with helical morphologies. Furthermore, the observed chiral transcription from molecules to morphologies is also significant for understanding bio-related self-assembly processes on the nano-to macro-scale.

Methods
Materials and physical measurements. R-or S-(1-phenylethylamino)methylphosphonic acid (pempH 2 ) was prepared according to methods reported in the literature 38 . All other starting materials were of reagent grade and were used as received from commercial sources without further purification. Elemental analyses for C, H and N were carried out on a PE 240C analyser. IR spectra were recorded with a Bruker Tensor 27 spectrometer using KBr discs. The pH value was measured by a Sartorius PB-10 pH metre. Thermal analyses were performed under nitrogen in the temperature range of 25-800°C at a heating rate of 5°C min −1 on a METTLER TOLEDO TGA/DSC 1 instrument. PXRD data were collected using a Bruker D8 advance diffractometer. SEM measurements were performed on SHI-MADZU SSX-550. The UV-Vis spectra were recorded on Perkin Elmer Lambda Preparation of rod-like crystals of S-3. Crystals of S-3 were obtained using the same procedure as for R-3 except that S-pempH 2 was used as the starting material. Elemental analysis calcd (%) for C 27 48 . Adsorption corrections were applied. The structures were solved using direct methods and refined on F 2 using full-matrix least-squares using SHELXTL 49 . Anisotropic temperature factors were used to refine all atoms excluding hydrogen. All hydrogen atoms bound to carbon were refined isotropically in riding mode; hydrogen atoms of water molecules were detected via experimental electron density and then refined isotropically with reasonable restriction of O-H bond distances and H-O-H angles. For R-3 and S-3, only one lattice water molecule with 0.5 occupancy could be determined from difference Fourier maps due to poor diffraction data or lattice solvent disorder even though the samples were sealed in mother solution for data collection. And the number of water molecules in lattice was determined by thermal analyses and elemental analyses. The crystallographic data are given in Supplementary Table 3, and the selected bond lengths and angles are shown in Supplementary Tables 4, 6, 8 and 11.
All-atom molecular dynamics simulations. The all-atom MD simulations were performed by using Gromacs 5.0.4 package 50 with the Amber force field 51 and the TIP3P water model 52 in the NVT ensemble. The force field parameter for R-pempH − and R-pempH 2 was built by using Antechamber tool 53 . The parameter for Tb 3+ was adopt by Li et al. 54 and that for NO 3 − was adopted by Liu et al. 55 During the simulation, the temperature was coupled at 300 K using Nosé-Hoover method 56 . The Particle Mesh Ewald (PME) method was used to calculate the electrostatic interactions and the cutoff of Lennard-Jones (LJ) interaction was 1.2 nm. The periodic boundary conditions were applied in all three dimensions. The time step was chosen as 2 fs and each simulation was at least conducted for 50 ns.
Brownian dynamics simulations. In order to better understand the distinct assembly under different cases, the BD simulations 57 were also used. Actually, since there existed a large number of chain-II and/or chain-III in this system, it was far beyond the computing ability of present all-atom molecular simulation. Here, the Tb 3+ was treated as one (green) bead (i.e., Tb bead) and the R-pempH − /R-pempH 2 molecule was simplified as two beads, where the ochre bead (i.e., P1 bead) can interact with each other to promote the growth along the side direction and the pink (for chain-II) or lime (for chain-III) bead (i.e., P2 bead) can interact with Tb bead to promote the growth along the chain direction (Fig. 9a). The harmonic spring interaction U s ¼ k s ðl i;iþ1 À l 0 Þ 2 was applied between connected beads in the polymers and the binding sites of particle beads, where k s = 3000 k B T 0 r 0 −2 , l 0 = 0.4 r 0 . In order to depict the attractive interaction of the Tb bead-P2 bead, the P1 bead-P1 bead, the Lennard-Jones potential was used (σ = 0.4 r 0 ). According to the allatom simulation, the interaction energy between Tb 3+ and R-pempH − was fifteen times as large as that between R-pempH − and R-pempH − (R-pempH 2 ) molecules. Here for the sake of simplicity, we set ϵ p1−p1 /k B T 0 = 1.0 and ϵ Tb−p2 /k B T 0 = 15.0. In particular, when in the presence of NO 3 − , the all-atom simulation results showed that the interaction energy between R-pempH − and R-pempH 2 molecules was three times as large as that in the absence of NO 3 − . Thus, we set the ϵ p1−p1 /k B T 0 = 3.0 in chain-II when in the presence of NO 3 − . Besides, the shifted Lennard−Jones potential, cutoff at 2 1/6 σ, was used to model the repulsive interaction between the other beads (ϵ/k B T 0 = 1.0, σ = 0.4 r 0 ). Since there are six Tb 3+ ions and eighteen R-pempH − /R-pempH 2 molecules in each pitch of chain-II/chain-III, the monomer for the assembly in the simulations was chosen as six Tb beads with eighteen ochre beads and 18 pink or lime beads (see Fig. 9b, c). In addition, the initial coordinate of Tb beads and P1 beads (in the monomer) was obtained by using the coordinate of Tb 3+ ion and P atom (of R-pempH − molecule) in experimental crystal data, respectively. Moreover, Tb beads and P1 beads in each monomer were treated as rigid ones in BD simulation, namely, the relative position of Tb beads and P1 beads in the monomer was kept fixed, which was due to the following two reasons. First, as shown in Supplementary Fig. 42, the all-atom MD simulation results showed that the relative position of Tb 3+ ion and P atom changed very little since the peak of the Tb-P distance and the P-Tb-P angle was very sharp in the probability distribution functions. Actually, the values of these parameters in all-atom MD simulations were close to those in the experiments (see the caption in Supplementary Fig. 42). Second, if we did not constrain the relative position of Tb beads and P1 beads in the monomer, the molecular structure obtained in BD simulation could be totally different from that in the all-atom MD simulation, and more importantly, differed from the unit in the experimental structure (e.g., chain-II and chain-III cannot be distinguished). As a result, here we treat the Tb beads and P1 beads in each monomer as rigid ones in the BD simulations.
All BD simulations were performed in the NVT ensemble by using the LAMMPS package (15 May 2015) 58 . During the simulation, as did in the experiments, the temperature was firstly coupled at 1.3 T 0 (about 390 K or 120°C) in the initial 100,000 τ and then decreased to 1.0 T 0 (about 300 K, i.e., room temperature) using the Langevin thermostat. The time step was 0.008 τ, and the data were collected every 80 τ, with the total simulation time lager than 400,000 τ.