Abstract
Despite the several available strategies to build complex supramolecular constructs, only a handful of different molecular knots have been synthesised so far. Here, in response to the quest for further designable topologies, we use Monte Carlo sampling and molecular dynamics simulations, informed by general principles of supramolecular assembly, as a discovery tool for thermodynamically and kinetically accessible knot types made of helical templates. By combining this approach with the exhaustive enumeration of molecular braiding patterns applicable to more general template geometries, we find that only few selected shapes have the closed, symmetric and quasiplanar character typical of synthetic knots. The corresponding collection of admissible topologies is extremely restricted. It covers all known molecular knots but it especially includes a limited set of novel complex ones that have not yet been obtained experimentally, such as 10_{124} and 15n_{41185}, making them privileged targets for future selfassembling experiments.
Introduction
In the past two decades there have been major breakthroughs in understanding the principles of directed selfassembly^{1,2,3,4,5,6,7,8,9,10}, as well as for the synthesis of molecules with complex topologies^{11,12,13,14}. Besides being appealing from a fundamental research perspective, molecular knots and links have versatile applications^{11,14,15,16}, from anion sensing and assisted catalysis^{17,18,19}, to seeding weaving structures^{20} to their potential as molecular nanocages^{21}, functionalizable nanoscale scaffolds^{22,23,24}, and the synthesis of polymers with special properties^{25,26,27,28,29,30,31}. The challenges encountered in designing and then obtaining knotted molecular architectures lie in the numerous concurrent physicochemical mechanisms that need to be balanced and steered. These include the length, thickness, intrinsic curvature and relative spatial orientation of the ligands and, above all, a judicious choice of the target topology^{32,33,34,35,36,37}. Not all knot types are, in fact, equally designable neither a priori nor in practice^{38}.
Until recently^{39}, the repertoire of addressable molecular topologies consisted of only the three simplest knot types: the trefoil^{40}, the figureofeight^{13} and the pentafoil knots^{41}, as well as the first fundamental link types^{15,35,42,43,44}. The 3_{1} and 5_{1} torus knots were assembled out of linear or circular double helicates whose crossings were templated by metal ions^{40,41}. The 4_{1} twist knot was instead assembled from a solution of organic flexiblerigid building blocks^{13}. Independently of the specific design strategy, the signature characteristics of these synthetic knots are two: cyclic symmetry and quasiplanar geometry.
In parallel with these experimental endeavours, theoretical and computational models have increasingly being used to optimise supramolecular design strategies and identify novel addressable topologies. To this end, our recent contribution was to devise general models based on selfassembling processes^{5,45,46} where a suitable shaping of the rigid building blocks boosted the spontaneous assembly of knotted or linked constructs^{38,47}. This minimalistic approach allowed for identifying a complex eightcrossing knot, the 8_{19} torus knot, as one of the most designable undiscovered knot types and hence the ideal next target of designed synthetic topologies. In the recent experiment of ref. ^{39}, Danon et al. have precisely succeeded in realising this very same type of 8_{19} molecular knot, thus underscoring the predictive power and reliability of general models of topological selfassembly.
At the same time, no simple design strategy exists for identifying additional privileged topologies obtainable with supramolecular techniques.
Motivated by these observations, here we use models and extensive stochastic simulations, informed with general principles of supramolecular assembly, for a systematic survey of thermodynamically and kinetically accessible knot types and use it to discover novel topologies that have not yet been obtained experimentally.
Our main finding is that the unsupervised selection procedure, besides recovering previously synthesised molecular knot types (3_{1}, 4_{1}, 5_{1} and 8_{19}), can indeed identify novel privileged topologies such as the 10_{124} and 15n_{41185} prime knots. These topologies have a higher complexity than those realised before and yet they share the same underlying architectural simplicity: they are quasiplanar, circularsymmetric and can be formed by only five templates of helical shape. To extend considerations beyond strictly helical templates, we survey the combinatorics of molecular braiding diagrams and show that the two mentioned topologies are part of a wider population of cyclicsymmetric knots that we propose as privileged addressable topologies for future molecular constructs. These include the 8_{18} topology that was reported experimentally^{48} after submission of this study.
Results
Model system
The templates, or building blocks, used for our selfassembling survey are modelled after helicates, which are commonly used in supramolecular synthetic chemistry^{17,40,43,49}.
Specifically, we consider helical fragments of unit radius, angular span α and pitch h. They are discretised as chains of beads with attractive patches at their termini, see Methods. The steric repulsion of the beads and the attractive patchy interactions allows the templates to attach one to another in a stringlike fashion^{50}.
The helical shape of the templates, with its nonplanar curved geometry, is ideally suited for our survey because it combines the simplicity of its mathematical parametrisation with the possibility to establish complex threedimensional structures either open or closed. For instance, prior to its elegant experimental realisation^{39}, the 8_{19} knot had been found to be designable precisely using helical templates in simulations of aspecific assembly processes^{38}.
Here, we instead consider assemblies formed under specific constraints, common to known molecular knot types, by helical templates with angular span and pitch varying respectively in the ranges 1.4π ≤ α ≤ 1.9π and 0.1 ≤ h ≤ 2.0, the unit of length being the helical radius. Combinations of templates with different chiralities are considered, too.
Monte Carlo discovery of novel addressable topologies
We proceeded in two steps, as sketched in Fig. 1. First, we explored the repertoire of thermodynamically accessible topologies of closed constructs formed by a fixed number of helical templates, n_{T}. Next, to identify the candidate shapes of feasible realisation, we select a posteriori those having the signature features of existing synthetic molecular knots: cyclic symmetry and quasiplanarity.
For the first step, we used a Monte Carlo procedure to explore the equilibrium conformational space accessible to constructs made of n_{T} = 3, 4 and 5, templates, as typical for known molecular knots. For a given value of n_{T}, we independently generated 1000 constructs for each template shape, i.e., (h, α), and possible combinations of left of righthanded template chiralities. The second step involved identifying the conformers with the sought cyclic symmetry. This was established by requiring that the constructs admitted one or more precise structural superpositions^{51} with their circular permutants, see Methods.
Figure 2 summarises the outcome of this phenomenologically informed discovery procedure. It presents the repertoire of the different types of closed, cyclicsymmetric knotted constructs found across the entire range of shapes and chiralities of the building blocks.
The graphical table of Fig. 2 displays the structural representatives of each knot type. Strikingly these are all inherently oblate, meaning that quasiplanarity is coopted by the cyclic symmetry constraint and does not need to be enforced additionally. Their Cartesian coordinates can be found in Supplementary Data 1.
Each representative in Fig. 2 is accompanied by a number of key properties, namely the braid representation of the knot, in both linear and circular form, and the topological designability score. We define the latter as the number of distinct template shapes, i.e., distinct points in the discretised (h,α) space that can assemble in that particular cyclic knot. This entropic measure is reflective of the designability of the knot that is the robustness to variations in template shapes. This ought to be a valuable quantitative element for guiding the design of novel molecular topologies. The regions in the (h, α) parameter space where the various knot types occurr, are highlighted in the topological state diagram of Supplementary Figs. 1 and 2.
The main result of Fig. 2 is that, across the wide range of template shapes and combinatorics of sampled constructs, the phenomenological selection is survived by only few privileged topologies.
These cover all knot types that have been synthesised so far and, above all, include different ones too. Entries in Fig. 2 that correspond to known molecular topologies are the 3_{1}, 4_{1}, 5_{1} and 8_{19} knots. The latter two are obtained with 5 and 4 templates, respectively, as in experimental realisations ^{39,41}. Interestingly, the 3_{1} knot is observed in different geometries, including the ideallike and the twisted one obtained with 3 and 4 templates, respectively, again as in the synthetic constructs^{40,52}. These instances are all torus knots^{53}.
The only nontorus knot in Fig. 2 is the 4_{1}, or figureofeight knot. Interestingly, one notes that this twist achiral knot can be established with four templates, either with the same or opposite chirality. The latter, racemic combination is more suitable because it covers a significantly wider region of parameter space, see designability score and Supplementary Figs. 1 and 2, and yields more planar constructs, see Supplementary Data 1 with coordinates. The 4_{1} topology has been experimentally obtained before, too, and with the same number of templates, though the building blocks were not helicates but flexible diblock modules^{13}. The 4_{1} instance in Fig. 2 therefore makes the important point that this topology would be realisable with nonnegligible probability also by using metal templating technique and helicates, but the latter should ideally be of opposite chirality. This is a feature that, to our knowledge, has not yet been systematically explored.
The fact that these previously obtained molecular knots are all included in the graphical table, gives confidence in the feasibility to synthesise in the future the two novel topologies appearing in Fig. 2.
These are the 10_{124} and 15n_{41185} knots. Though being made by only n_{T} = 5 templates (same chirality), like the 5_{1} synthetic knot, these topologies are more complex than all known synthetic knots obtained so far. In fact, it is arresting to observe that, out of the millions of prime and composite knots with of up to 15 crossings, those feasible with a handful of cyclically arranged helical templates are only 6.
The 10_{124} topology emerges as a particularly designable cyclic knot, because it can be obtained with tens of different template shapes (Supplementary Fig. 2). It is intriguing that this topology is among those that are only rarely populated in aspecific selfassemblies^{38}, i.e., without cyclic symmetry or quasiplanarity, and appears, with other instances, in energyminimizing arrangements of dipolar particles^{54}.
The 15crossing knot is an even more remarkable example of sophisticated interplay of geometry and topology, with as many as 15 crossings established by 5 templates only. It forms abundantly for a specific template shape, the one involved in the construct shown in Fig. 2, and more limitedly for a second template geometry, see Supplementary Fig. 2.
Of these two novel complex topologies, we therefore single out 10_{124} as the most promising one to be realised with current selfassembly strategies.
Kinetic accessibility of selfassembling knotted constructs
To test the kinetic accessibility of the topologies discovered by the Monte Carlo sampling, we performed hundreds of selfassembly simulations varying templates geometries. We used Langevin molecular dynamics (MD) to evolve dispersions of patchy helical templates in various conditions, see Fig. 3a, Supplementary Fig. 3 and Methods. To mimic helicate coordination by templating metal ions^{55}, we also considered adding several coordinating particles (same diameter of the beads) that do not bind to each other but can attract, and hence bring together, two or more templates.
By exploring combinations of varying number of templates and coordinating particles, from few to hundreds, see Supplementary Fig. 3, we established that the simplest selfassembling strategy for producing cyclicsymmetric quasiplanar constructs was to focus on systems with an equal number of templates and coordinating particles, n_{T} = 3, 4, 5 at a suitable density, see Methods.
With these systematic selfassembly simulations, we established that all privileged topologies listed in Fig. 2 except for the 5_{1} knot (because of its large contact angle, see Supplementary Fig. 4), are indeed kinetically accessible. Their representative structures formed by selfassembly are shown in Fig. 3b, while their cartesian coordinates are available in Supplementary Data 1. These conformers include the novel complex target topologies, 10_{124} and 15n_{41185}. This reinforces their general viability as new targets for molecular constructs with addressable topology.
Interestingly, besides the topologies of Fig. 2, we recurrently observed the formation of the 12n_{242} topology, as shown in Fig. 3b and Supplementary Table 1. This knotted structure is selfassembled from five templates and differs from the 15n_{41185} by a localised defect in the otherwise regular pattern of over and undercrossings. From its abundance, we conclude that the 12n_{242} topology, though not cyclicsymmetric, might be obtainable as a likely byproduct of the target 15n_{41185} knot.
Enumerative survey of cyclicsymmetric entangled constructs
The geometrical and topological repertoire spanned by the constructs of Figs. 2 and 3b can be recapitulated in terms of only two key parameters: the number of templates, n_{T}, which is equal to the number of crests on each edge of the corresponding braid, and the number of strands in the braid itself, n_{S}. The latter matches the number of complete turns around the cyclic symmetry axis, but is more general because it can be defined for noncyclic knots too. Both parameters have, in fact, proved useful in systematic explorations of possible synthetic topologies^{39}, arguably because they reflect different aspects of the practical difficulty of their realisation.
Classifying conformers in terms of n_{T} and n_{S} is an apt way to code their structure independently of the specific geometry of their templates. For this reason, a systematic search of the (n_{T}, n_{S}) parameter space can reveal additional addressable topologies besides those in Figs. 2 and 3b made of strictly helical building blocks or, else, rule out the existence of alternative ones.
Accordingly, we completed the survey with an exhaustive enumeration of braids for different combinations of (n_{T}, n_{S}). For practical reasons, we limited the combinatorial search up to n_{T} = 7 (see Supplementary Fig. 11 for n_{T} = 8, 9 and 10 too) and we assumed n_{T} ≥ n_{S} + 1, meaning that viable templates should cover less than a full turn when projected.
The results are shown in Fig. 4. Note that, since the search scheme based on braiding combinatorics is more general than the Monte Carlo exploration, the topologies listed for each (n_{T}, n_{S}) pair are not limited to those in Fig. 2. They include, for example, the 8_{18} and 10_{123} knots that, owing to their alternating character, are not realisable with a small number of templates if these have a strictly helical geometry. Moreover, when n_{T} and n_{S} are not mutually prime they include multicomponent constructs, such as catenanes and links. In such cases, the number of components corresponds to the greatest common divisor of n_{T} and n_{S}.
For both knots and links in Fig. 4, the number of projected crossings is equal to n_{T} × (n_{S} − 1). Note that this geometrical measure of complexity, which does not necessarily coincide with the socalled minimal crossing number, is identical for all topologies realisable with the same number of templates and strands, (n_{T}, n_{S}). For instance, various 8crossing knots are obtained for (n_{T} = 4, n_{S} = 3). The cyclic ones, shown in Fig. 4, are the 8_{19} and 8_{18} knots, with C_{4} symmetry, and the 3_{1} and 4_{1} ones, with C_{2} symmetry. Noncyclicsymmetric instances are listed but not shown in Fig. 4 and include nonminimal representations of 5_{1} and 5_{2} topologies. Similarly, for (n_{T} = 5, n_{S} = 4) one obtains various 15crossing knots: four cyclic ones, including the privileged 15n_{41185} topology, and many more acyclic ones, such as the 12n_{242} knot that recurred in selfassembling simulations, see Fig. 3b. Nonsymmetric knots are listed in Supplementary Figs. 8–11.
The topology indexing scheme of Fig. 4 has a twofold implication.
First, it systematically recapitulates the repertoire of designable topologies in terms of two parameters with direct bearings on the complexity of their practical realisation: the number of templates and the number of strands. The key emerging point is that ranking topologies in terms of realisation complexity, by the number of templates and braided strand, subverts the canonical order of nominal topological complexity^{56}. For instance, cyclicsymmetric realisations of knots with 7 minimal crossings, such as 7_{1}, requires more templates than prime knots with 8 or 10 minimal crossings (e.g., 8_{19} or 10_{124}). At the same time, complex knots with 8, 10, 14 and 15 crossings can be realised with fewer templates than topologies with 6 or 7 projected crossings. This stresses the fact that the expected difficulty of realisation does not necessarily parallel the nominal complexity, and hence reinforces a posteriori the necessity of systematic survey for judicious choices of target topologies^{57}.
Second, it provides a systematic route towards discovering new designable topologies by extending the range of n_{S} and n_{T} even beyond the cases considered here. This search ought not to be aimed only at the simplest types of novel knots, because complex threedimensional geometries would be even better suited for specific goals, such as realising molecular cages^{21,58}. Arguably, Sshaped^{34} or otherwise wavy building blocks ought to be more suited than constantcurvature ones (including straight and helical templates) to produce the more intricate entries in Fig. 4, and particularly those with the largest number of braid strands, n_{S}, for a given number of templates, n_{T}.
Discussion
In summary, we reported on a systematic discovery scheme for novel molecular topologies. The method combines Monte Carlo, molecular dynamics and braiding patterns enumeration and was used to singleout knot types that have the same signature features, notably cyclic symmetry and quasiplanarity, shared by all known synthetic molecular knots.
It was thus established that such repertoire of admissible topologies includes only a tiny fraction of all possible knot types. In particular, there are only 6 distinct cyclicsymmetric knot types that can be assembled with 5 or fewer helicatelike templates. Four of them, namely 3_{1}, 4_{1}, 5_{1} and 8_{19}, have been previously obtained in a remarkable progression of synthesis strategies spanning two decades; the other two, instead, are yet to be realised experimentally. These correspond to the 10_{124} and 15n_{41185} knot types. With their ten and 15 minimal crossings, respectively, these topologies largely surpass previous ones in terms of complexity. At the same time, their characteristics of cyclic symmetry, quasiplanarity and kinetic accessibility makes them ideal targets of future molecular designs efforts.
Finally, by exhaustive enumeration of braid patterns, we found that only few and specific types of additional knots might be addressable by using either a larger number of helical templates, or few templates but with wavy or Sshaped geometries^{34}, either rigid^{37} or with flexible termini^{13}.
The simplest of these privileged topologies appear in Fig. 4. The shown targets include the 8_{18} topology, which was reported experimentally^{48} after submission of this study. This molecular knot was obtained with a nonminimal geometry, corresponding to (n_{t} = 8, n_{s} = 3) in the scheme of Fig. 4, see Supplementary Fig. 11. This match confirms the predictive capabilities of the enumerative scheme of Fig. 4.
The in silico exploration and optimization of templates shapes ought to be valuable also in supramolecular DNA assembling strategies^{6,59,60,61} that are characterised by a good control of the local curvature of templates as well as of their spatial coordination and binding pattern.
Methods
Building blocks
As building blocks for selfassembling structures we consider rigid helical fragments. The helices project a circle of unit radius and their pitch h is varied in the [0.1;2.0] range at 0.1 increments. The angular span of the fragments, α is varied in the [1.4;1.9]π range in steps of 0.1π. Both possible chiralities χ are considered: right (+1) or left (−1). The helical trace of the fragments is described by the parametric equation {x = cos(αt); y = sin(αt); z = χht;} where t ∈ [0, 1]. On this helical trace we then placed the centres of spherical beads with diameter equal to σ = 1/3. The centres are equispaced at a distance in the [σ, 2^{1/6}σ] range, as needed to fit an integer number of beads along the contour. Following ref. ^{38} each template is functionalised with two small attractive patches (white spheres in Figures, enlarged for visual clarity) lying on the exposed surface of the terminal beads at the intersection point with the helical centerline.
Monte Carlo sampling
A Monte Carlo procedure based on crankshaft moves is used to explore the conformational space of closed constructs made of n_{T} = 3, 4, 5 templates. These constructs are formed by identical rigid templates, i.e., same (h, α) geometry, joined at their patchy ends in a circular fashion, see Fig. 1. For each template geometry and template number, all combinations of left and righthanded chiralities are considered. In the Monte Carlo scheme, bonded templates are joined by exactly superposing the patches at their termini. The connectivity of the constructs is preserved by using unrestricted crankshaft moves that are hinged at randomly chosen pairs of these superposed patches, see Fig. 1. Being their spatial superposition kept fixed at all times, the patches, and the templates too, cannot detach and therefore no patchpatch attractive potential is introduced. All configurations are accepted, except those with steric clashes (overlapping beads of different templates) or with large (>π/4) contacting angles between consecutive templates, which are rejected. For each combination of template shape, (h, α), and chirality, we sampled 1000 configurations spaced by 5 × 10^{4} Monte Carlo (crankshaft) moves, a timespan larger than the typical autocorrelation time for the considered systems.
From the set of sampled conformations we then selected a posteriori those, if any, with approximate cyclic symmetry. To this end, we structurally aligned the constructs with their circularly permuted variants, and kept track of the corresponding root mean square deviation (RMSD) of the alignments^{51}. The sought cyclicsymmetric constructs were identified from the atypically low RMSD, see Supplementary Fig. 5.
Molecular dynamics simulations
We used Langevin molecular dynamics simulations to study the spontaneous selfassembly of helical templates into complex constructs. The system used in the MC approach is here enriched by a number of particles in solution that do not bind to each other but can attract, and hence coordinate, the templates. The role of these coordinating particles is to mimic the effects of metal ion templating^{55}. Various simulations were carried out with a wide range of number (and relative size too) of templates and coordinating particles, from few to several hundreds, see Supplementary Fig. 3, as well as with different interaction parameters. These simulations were aimed at identifying the most suitable conditions yielding cyclicsymmetric quasiplanar constructs made of few templates. The choice of the specific parameters values was mostly guided by the need to make the assembly process computationally amenable.
The simplest conditions, those used for the production runs, are detailed below and correspond to systems with a small and equal number of templates and coordinating particles, n_{T} = 3, 4, 5, placed in periodic simulation boxes of side a few times larger than the templates’ projected diameters.
With exception of the small patches all other ones (templates beads and the coordinating particles) interact sterically through a truncated and shifted LennardJones potential:
where \(\epsilon\), the unit energy, is equal to the thermal energy of the system K_{B}T, σ is the particle diameter, d is the distance of the particles’ centres, C_{LJ} = 100, and θ is the Heaviside function.
The patches, instead, interact via an attractive Gaussian potential well,
where C_{p} = 25 and its range σ_{p} = 0.1σ, as in ref. ^{38}. The combination of these interaction terms results in an attractive directional potential well of depth ~18 K_{B}T between helical templates, see Supplementary Fig. 6.
The attraction of the coordinating particles and the templates’ beads is described with a Yukawalike potential
where C_{Y} = 8 and l_{Y} = σ. A repulsive version of this potential (same form but opposite sign) is introduced between the coordinating particles to keep them apart.
Simulations were carried out in a periodic cubic box of size L, with the LAMMPS software package^{62} with default values for the beads mass and viscous friction and integration time step equal to 0.006τ_{LJ}, where τ_{LJ} is the characteristic LennardJones time.
For each value of n_{T} and template geometry we carried out 20 independent simulations at fixed template number density ρ_{T} ≡ n_{T}/L^{3} = 0.0125, where L is in units of the helical radius. Note that the duration of each simulation was 264,000τ_{LJ}, i.e., sufficiently long to observe spontaneous binding and assembling of templates, but still shorter than their typical unbinding time (~4,000,000 τ_{LJ}, Supplementary Fig. 7).
Knot analysis
For most constructs, the patterns of crossings in a twodimensional projection was simple enough that the Alexander polynomial sufficed to identify the knot type^{63}. For more complicated patterns, the Dowker code was compared against a lookup table for prime knots of up to 16 crossings using the Knotscape software package of J. Hoste and M. Thistletwaite available from www.math.utk.edu/morwen/knotscape.html.
Data availability
The coordinates of the representative conformers shown in Figs. 2 and 3 are available as Supplementary Data 1.
The custom code for the Monte Carlo exploration is available from the authors. The code is made available freely for academic use and without user support.
References
 1.
Mallory, S. A., Valeriani, C. & Cacciuto, A. An active approach to colloidal selfassembly. Ann. Rev. Phys. Chem. 69, 59–79 (2017).
 2.
Furst, E. M. Directed selfassembly. Soft Matter 9, 9039–9045 (2013).
 3.
Sacanna, S. et al. Shaping colloids for selfassembly. Nat. Commun. 4, 1688 (2013).
 4.
Šarić, A. & Cacciuto, A. Selfassembly of nanoparticles adsorbed on fluid and elastic membranes. Soft Matter 9, 6677–6695 (2013).
 5.
Damasceno, P. F., Engel, M. & Glotzer, S. C. Predictive selfassembly of polyhedra into complex structures. Science 337, 453–457 (2012).
 6.
Aldaye, F. A., Palmer, A. L. & Sleiman, H. F. Assembling materials with dna as the giude. Science 321, 1795–1799 (2008).
 7.
Iacovella, C. R., Keys, A. S. & Glotzer, S. C. Selfassembly of softmatter quasicrystals and their approximants. Proc. Natl Acad. Sci. USA 108, 0935–20940 (2011).
 8.
Romano, F. & Sciortino, F. Colloidal selfassembly: patchy from the bottom up. Nat. Mater. 10, 171 (2011).
 9.
Chen, Q., Bae, S. C. & Granick, S. Directed selfassembly of a colloidal kagome lattice. Nature 469, 381 (2011).
 10.
Sciortino, F., Giacometti, A. & Pastore, G. Phase diagram of janus particles. Phys. Rev. Lett. 103, 237801 (2009).
 11.
Sauvage, J.P. & DietrichBuchecker, C. Catenanes, Rotaxanes and Knots. (VCHWiley, Weinheim, Germany, 1999).
 12.
Ayme, J.F., Beves, J. E., Campbell, C. J. & Leigh, D. A. Template synthesis of molecular knots. Chem. Soc. Rev. 42, 1700–1712 (2013).
 13.
Ponnuswamy, N., Cougnon, F. B. L., Dan Pantos, G. & Sanders, J. K. M. Homochiral and meso figure eight knots and a solomon link. J. Am. Chem. Soc. 136, 8243–8251 (2014).
 14.
GilRamirez, G., Leigh, D. A. & Stephens, A. J. Catenanes: fifty years of molecular links. Angew. Chem. Int. Ed. 54, 6110–6150 (2015).
 15.
Chichak, K. S. et al. Molecular borromean rings. Science 304, 1308–1312 (2004).
 16.
Horner, K. E., Miller, M. A., Steed, J. W. & Sutcliffe, P. M. Knot theory in modern chemistry. Chem. Soc. Rev. 45, 6432–6448 (2016).
 17.
Marcos, V. et al. Allosteric initiation and regulation of catalysis with a molecular knot. Science 352, 1555–1559 (2016).
 18.
Ayme, J.F. et al. Strong and selective anion binding within the central cavity of molecular knots and links. J. Am. Chem. Soc. 137, 9812–9815 (2015).
 19.
Bilbeisi, R. A. et al. [CH…anion] interactions mediate the templation and anion binding properties of topologically nontrivial metalorganic structures in aqueous solutions. Chem. Sci. 7, 2524–2531 (2016).
 20.
Liu, Y. et al. Weaving of organic threads into a crystalline covalent organic framework. Science 351, 365–369 (2016).
 21.
Riddell, I. A. et al. Cation and anionexchanges induce multiple distinct rearrangements within metallosupramolecular architectures. J. Am. Chem. Soc. 136, 9491–9498 (2014).
 22.
Chichak, K. S., Peters, A. J., Cantrill, S. J. & Stoddart, J. F. Nanoscale borromeates. J. Org. Chem. 70, 7956–7962 (2005).
 23.
Yates, C. R., Benítez, D., Khan, S. I. & Stoddart, J. F. Hexafunctionalized borromeates using olefin cross metathesis. Org. Lett. 9, 2433–2436 (2007).
 24.
Lukin, O. & Vögtle, F. Knotting and threading of molecules: chemistry and chirality of molecular knots and their assemblies. Angew. Chem. Int. Ed. 44, 1456–1477 (2005).
 25.
Yamamoto, T. & Tezuka, Y. Cyclic polymers revealing topology effects upon selfassemblies, dynamics and responses. Soft Matter 11, 7458–7468 (2015).
 26.
Cao, P.F., Mangadlao, J. & Advincula, R. A trefoil knotted polymer produced through ring expansion. Angew. Chem. Int. Ed. 54, 5127–5131 (2015).
 27.
Cao, P.F., Mangadlao, J. D., de Leon, A., Su, Z. & Advincula, R. C. Catenated poly (εcaprolactone) and poly (llactide) via ringexpansion strategy. Macromolecules 48, 3825–3833 (2015).
 28.
Yamamoto, T., Yagyu, S. & Tezuka, Y. Lightand heattriggered reversible linear–cyclic topological conversion of telechelic polymers with anthryl end groups. J. Am. Chem. Soc. 138, 3904–3911 (2016).
 29.
Sato, H., Aoki, D. & Takata, T. Synthesis and star/linear topology transformation of a mechanically linked abc terpolymer. ACS Macro Lett. 5, 699–703 (2016).
 30.
Cao, P.F., Rong, L.H., Mangadlao, J. D. & Advincula, R. C. Synthesizing a trefoil knotted block copolymer via ringexpansion strategy. Macromolecules 50, 1473–1481 (2017).
 31.
Wu, Q. et al. Poly[n]catenanes: synthesis of molecular interlocked chains. Science 358, 1434–1439 (2017).
 32.
Prakasam, T. et al. Simultaneous selfassembly of a [2]catenane, a trefoil knot, and a solomon link from a simple pair of ligands. Angew. Chem. Int. Ed. 52, 9862–9862 (2013).
 33.
GilRamirez, G. et al. Tying a molecular overhand knot of single handedness and asymmetric catalysis with the corresponding pseudod3symmetric trefoil knot. J. Am. Chem. Soc. 138, 13159–13162 (2016).
 34.
Beves, J. E., Danon, J. J., Leigh, D. A., Lemonnier, J.F. & VitoricaYrezabal, I. J. A solomon link through an interwoven molecular grid. Angew. Chem. Int. Ed. 54, 7555–7559 (2015).
 35.
Cougnon, F. B. L., Ponnuswamy, N., Pantos, G. D. & Sanders, J. K. M. Molecular motion of donoracceptor catenanes in water. Org. Biomol. Chem. 13, 2927–2930 (2015).
 36.
Meyer, C. D. et al. The dynamic chemistry of molecular borromean rings and solomon knots. Chem. Eur. J. 16, 12570–12581 (2010).
 37.
Schouwey, C. et al. Selfassembly of a giant molecular solomon link from 30 subcomponents. Angew. Chem. Int. Ed. 53, 11261–11265 (2014).
 38.
Polles, G., Marenduzzo, D., Orlandini, E. & Micheletti, C. Selfassembling knots of controlled topology by designing the geometry of patchy templates. Nat. Commun. 6, 6423 (2015).
 39.
Danon, J. J. et al. Braiding a molecular knot with eight crossings. Science 355, 159–162 (2017).
 40.
DietrichBuchecker, C. O. & Sauvage, J. P. A synthetic molecular trefoil knot. Angew. Chem. Int. Ed. 28, 189–192 (1989).
 41.
Ayme, J.F. et al. A synthetic molecular pentafoil knot. Nat. Chem. 4, 15–20 (2011).
 42.
Nierengarten, J. F., DietrichBuchecker, C. O. & Sauvage, J. P. Synthesis of a doubly interlocked [2]catenane. J. Am. Chem. Soc. 116, 375–376 (1994).
 43.
Leigh, D. A., Pritchard, R. G. & Stephens, A. J. A star of david catenane. Nat. Chem. 6, 978–982 (2014).
 44.
Pentecost, C. D. et al. A molecular solomon link. Angew. Chem. Int. Ed. 46, 218–222 (2007).
 45.
Romano, F. & Sciortino, F. Patterning symmetry in the rational design of colloidal crystals. Nat. Commun. 3, 975 (2012).
 46.
Ye, X. et al. Competition of shape and interaction patchiness for selfassembling nanoplates. Nat. Chem. 5, 466–473 (2013).
 47.
Polles, G., Orlandini, E. & Micheletti, C. Optimal selfassembly of linked constructs and catenanes via spatial confinement. ACS Macro Lett. 5, 931–935 (2016).
 48.
Kim, D. H. et al. Coordinationdriven selfassembly of a molecular knot comprising sixteen crossings. Angew. Chem. Int. Ed. 57, 5669–5673 (2018).
 49.
DietrichBuchecker, C. O., Sauvage, J. P. & Kern, J. M. Templated synthesis of interlocked macrocyclic ligands: the catenands. J. Am. Chem. Soc. 106, 3043–3045 (1984).
 50.
Tavares, J. M., Rovigatti, L. & Sciortino, F. Quantitative description of the selfassembly of patchy particles into chains and rings. J. Chem. Phys. 137, 044901 (2012).
 51.
Kabsch, W. A discussion of the solution for the best rotation to relate two sets of vectors. Acta Cryst. A 34, 827–828 (1978).
 52.
Ayme, J.F. et al. Lanthanide template synthesis of a molecular trefoil knots. J. Am. Chem. Soc. 136, 13142–13145 (2014).
 53.
Adams, C. C. The Knot Book (Freeman, 1994).
 54.
Miller, M. A. & Wales, D. J. Novel structural motifs in clusters of dipolar spheres: Knots, links, and coils. J. Phys. Chem. B 109, 23109–23111 (2005).
 55.
Forgan, R. S., Sauvage, J.P. & Stoddart, J. F. Chemical topology: complex molecular knots, links, and entanglements. Chem. Rev. 111, 5354–5364 (2011).
 56.
Rolfsen, D. Knots and Links. (Mathematics Lecture Series 7, Publish or Perish, Inc., Houston, Texas, 1976).
 57.
Chambron, J.C. & Sauvage, J.P. Topologically complex molecules obtained by transition metal templation: it is the presentation that determines the synthesis strategy. New J. Chem. 37, 49–57 (2013).
 58.
Leenders, S. H. A. M., GramageDoroia, R., de Bruin, B. & Reek, J. N. H. Transition metal catalysis in confined spaces. Chem. Soc. Rev. 44, 433–448 (2015).
 59.
Rothemund, P. W. K. Folding DNA to create nanoscale shapes and patterns. Nature 440, 297–302 (2006).
 60.
Han, D., Pal, S., Liu, Y. & Yan, H. Folding and cutting dna into reconfigurable topological nanostructures. Nat. Nanotech. 5, 712–717 (2010).
 61.
Zhang, F., Nangreave, J., Liu, Y. & Yan, H. Structural dna nanotechnology: state of the art and future perspective. J. Am. Chem. Soc. 136, 11198–11211 (2014).
 62.
Plimpton, S. J. Fast parallel algorithms for shortrange molecular dynamics. J. Comput. Phys. 117, 1–19 (1995).
 63.
Orlandini, E. & Whittington, S. G. Statistical topology of closed curves: some applications in polymer physics. Rev. Mod. Phys. 79, 611–642 (2007).
Acknowledgements
We acknowledge support from the Italian Ministry of Education. We thank Guido Polles for valuable discussions, technical help and scientific advice.
Author information
Affiliations
Contributions
M.M., E.O. and C.M. designed and performed the research and wrote the manuscript.
Corresponding author
Ethics declarations
Competing interests
The authors declare no competing interests.
Additional information
Publisher's note: Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Electronic supplementary material
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Marenda, M., Orlandini, E. & Micheletti, C. Discovering privileged topologies of molecular knots with selfassembling models. Nat Commun 9, 3051 (2018). https://doi.org/10.1038/s4146701805413z
Received:
Accepted:
Published:
Further reading

A molecular endless (74) knot
Nature Chemistry (2021)

Tying different knots in a molecular strand
Nature (2020)
Comments
By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.