Discovering privileged topologies of molecular knots with self-assembling models

Unraveling the Influence of Topology and Spatial Confinement on Equilibrium and Relaxation Properties of Interlocked Ring Polymers

We use Langevin dynamics simulations to study linked ring polymers in channel confinement. We address the in- and out-of-equilibrium behavior of the systems for varying degrees of confinement and increasing topological and geometrical complexity of the interlocking. The main findings are three. First, metric observables of different link topologies collapse onto the same master curve when plotted against the crossing number, revealing a universal response to confinement. Second, the relaxation process from initially stretched states is faster for more complex links. We ascribe these properties to the interplay of several effects, including the dependence of topological friction on the link complexity. Finally, we show that transient forms of geometrical entanglement purposely added to the initial stressed state can leave distinctive signatures in force-spectroscopy curves. The insight provided by the findings could be leveraged in single-molecule nanochannel experiments to identify geometric entanglement within topologically linked rings.

Quantum-inspired encoding enhances stochastic sampling of soft matter systems

Quantum advantage in solving physical problems is still hard to assess due to hardware limitations. However, algorithms designed for quantum computers may engender transformative frameworks for modeling and simulating paradigmatically hard systems. Here, we show that the quadratic unconstrained binary optimization encoding enables tackling classical many-body systems that are challenging for conventional Monte Carlo. Specifically, in self-assembled melts of rigid lattice ring polymers, the combination of high density, chain stiffness, and topological constraints results in divergent autocorrelation times for real-space Monte Carlo. Our quantum-inspired encoding overcomes this problem and enables sampling melts of lattice rings with fixed curvature and compactness, unveiling counterintuitive topological effects. Tackling the same problems with the D-Wave quantum annealer leads to substantial performance improvements and advantageous scaling of sampling computational cost with the size of the self-assembled ring melts.

Can Polymer Helicity Affect Topological Chirality of Polymer Knots?

We investigate the effect of helicity in isolated polymers on the topological chirality of their knots with computer simulations. Polymers are described by generic worm-like chains (WLC), where helical conformations are promoted by chiral coupling between segments that are neighbors along the chain contour. The sign and magnitude of the coupling coefficient u determine the sense and strength of helicity. The corrugation of the helix is adjusted via the radius R of a spherical, hard excluded volume around each WLC segment. Open and compact helices are, respectively, obtained for R that is either zero or smaller than the length of the WLC bond, and R that is a few times larger than the bond length. We use a Monte Carlo algorithm to sample polymer conformations for different values of u, spanning the range from achiral polymers to chains with well-developed helices. Monitoring the average helix torsion and fluctuations of chiral order as a function of u, for two very different chain lengths, demonstrates that the coil-helix transition in this model is not a phase transition but a crossover. Statistical analysis of conformations forming the simplest chiral knots, 3₁, 5₁, and 5₂, demonstrates that topological mirror symmetry is broken─knots formed by helices with a given sense prefer one handedness over the other. For the 3₁ and 5₁ knots, positive helical sense favors positive handedness. Intriguingly, an opposite trend is observed for 5₂ knots, where positive helical sense promotes negative handedness. We argue that this special coupling between helicity and topological chirality stems from a generic mechanism: conformations where some of the knot crossings are found in "braids" formed by two tightly interwoven sections of the polymer.

Circular Polycatenanes: Supramolecular Structures with Topologically Tunable Properties

Polycatenanes, macrochains of topologically interlocked rings with unique physical properties have recently gained considerable interest in supramolecular chemistry, biology, and soft matter. Most of the work has been, so far, focused on linear chains and on their variety of conformational properties compared to standard polymers. Here we go beyond the linear case and show that, by circularizing such macrochains, one can exploit the topology of the local interlockings to store twist in the system, significantly altering its metric and local properties. Moreover, by properly defining the twist (Tw) and writhe (Wr) of these macrorings we show the validity of a relation equivalent to the Călugăreanu-White-Fuller theorem Tw+Wr=const, originally proved for ribbonlike structures such as double stranded DNA. Our results suggest that circular polycatenanes with storable and tunable twist can form a new category of highly designable multiscale structures with potential applications in supramolecular chemistry and material science.

Density-tunable pathway complexity in a minimalistic self-assembly model

An open challenge in self-assembly is learning how to design systems that can be conditionally guided towards different target structures depending on externally-controlled conditions. Using a theoretical and numerical approach, here we discuss a minimalistic self-assembly model that can be steered towards different types of ordered constructs at the equilibrium by solely tuning a facile selection parameter, namely the density of building blocks. Metadynamics and Langevin dynamics simulations allow us to explore the behavior of the system in and out of equilibrium conditions. We show that the density-driven tunability is encoded in the pathway complexity of the system, and specifically in the competition between two different main self-assembly routes. A comprehensive set of simulations provides insight into key factors allowing to make one self-assembling pathway prevailing on the other (or vice versa), determining the selection of the final self-assembled products. We formulate and validate a practical criterion for checking whether a specific molecular design is predisposed for such density-driven tunability of the products, thus offering a new, broader perspective to realize and harness this facile extrinsic control of conditional self-assembly.

Knotting matters: orderly molecular entanglements

Entangling strands in a well-ordered manner can produce useful effects, from shoelaces and fishing nets to brown paper packages tied up with strings. At the nanoscale, non-crystalline polymer chains of sufficient length and flexibility randomly form tangled mixtures containing open knots of different sizes, shapes and complexity. However, discrete molecular knots of precise topology can also be obtained by controlling the number, sequence and stereochemistry of strand crossings: orderly molecular entanglements. During the last decade, substantial progress in the nascent field of molecular nanotopology has been made, with general synthetic strategies and new knotting motifs introduced, along with insights into the properties and functions of ordered tangle sequences. Conformational restrictions imparted by knotting can induce allostery, strong and selective anion binding, catalytic activity, lead to effective chiral expression across length scales, binding modes in conformations efficacious for drug delivery, and facilitate mechanical function at the molecular level. As complex molecular topologies become increasingly synthetically accessible they have the potential to play a significant role in molecular and materials design strategies. We highlight particular examples of molecular knots to illustrate why these are a few of our favourite things.

Molecular communications in complex systems of dynamic supramolecular polymers

Supramolecular polymers are composed of monomers that self-assemble non-covalently, generating distributions of monodimensional fibres in continuous communication with each other and with the surrounding solution. Fibres, exchanging molecular species, and external environment constitute a sole complex system, which intrinsic dynamics is hard to elucidate. Here we report coarse-grained molecular simulations that allow studying supramolecular polymers at the thermodynamic equilibrium, explicitly showing the complex nature of these systems, which are composed of exquisitely dynamic molecular entities. Detailed studies of molecular exchange provide insights into key factors controlling how assemblies communicate with each other, defining the equilibrium dynamics of the system. Using minimalistic and finer chemically relevant molecular models, we observe that a rich concerted complexity is intrinsic in such self-assembling systems. This offers a new dynamic and probabilistic (rather than structural) picture of supramolecular polymer systems, where the travelling molecular species continuously shape the assemblies that statistically emerge at the equilibrium.

The dynamic structure of supramolecular polymers is challenging to determine both in experiments and in simulations. Here the authors use coarse-grained molecular models to provide a comprehensive analysis of the molecular communication in these complex molecular systems.

Topological Friction and Relaxation Dynamics of Spatially Confined Catenated Polymers

We study catenated ring polymers confined inside channels and slits with Langevin dynamics simulations and address how the contour position and size of the interlocked or physically linked region evolve with time. We show that the catenation constraints generate a drag, or topological friction, that couples the contour motion of the interlocked regions. Notably, the coupling strength decreases as the interlocking is made tighter, but also shorter, by confinement. Though the coupling strength differs for channel and slit confinement, the data outline a single universal curve when plotted against the size of the linked region. Finally, we study how the relaxation kinetics changes after one of the rings is cut open and conclude that considering interlocked circular polymers is key for isolating the manifestations of topological friction. The results ought to be relevant for linked biomolecules in experimental or biological confining conditions.

Topological and physical links in soft matter systems

Linking, or multicomponent topological entanglement, is ubiquitous in soft matter systems, from mixtures of polymers and DNA filaments packedin vivoto interlocked line defects in liquid crystals and intertwined synthetic molecules. Yet, it is only relatively recently that theoretical and experimental advancements have made it possible to probe such entanglements and elucidate their impact on the physical properties of the systems. Here, we review the state-of-the-art of this rapidly expanding subject and organize it as follows. First, we present the main concepts and notions, from topological linking to physical linking and then consider the salient manifestations of molecular linking, from synthetic to biological ones. We next cover the main physical models addressing mutual entanglements in mixtures of polymers, both linear and circular. Finally, we consider liquid crystals, fluids and other non-filamentous systems where topological or physical entanglements are observed in defect or flux lines. We conclude with a perspective on open challenges.

Piecewise-linear embeddings of knots and links with rotoinversion symmetry

This article describes the simplest members of an infinite family of knots and links that have achiral piecewise-linear embeddings in which linear segments (sticks) meet at corners. The structures described are all corner- and stick-2-transitive - the smallest possible for achiral knots.

A molecular endless (74) knot

Current strategies for the synthesis of molecular knots focus on twisting, folding and/or threading molecular building blocks. Here we report that Zn(II) or Fe(II) ions can be used to weave ligand strands to form a woven 3 × 3 molecular grid. We found that the process requires tetrafluoroborate anions to template the assembly of the interwoven grid by binding within the square cavities formed between the metal-coordinated criss-crossed ligands. The strand ends of the grid can subsequently be joined through within-grid alkene metathesis reactions to form a topologically trivial macrocycle (unknot), a doubly interlocked [2]catenane (Solomon link) and a knot with seven crossings in a 258-atom-long closed loop. This 7₄ knot topology corresponds to that of an endless knot, which is a basic motif of Celtic interlace, the smallest Chinese knot and one of the eight auspicious symbols of Buddhism and Hinduism. The weaving of molecular strands within a discrete layer by anion-template metal-ion coordination opens the way for the synthesis of other molecular knot topologies and to woven polymer materials.

Effects of turn-structure on folding and entanglement in artificial molecular overhand knots

The length and constitution of spacers linking three 2,6-pyridinedicarboxamide units in a molecular strand influence the tightness of the resulting overhand (open-trefoil) knot that the strand folds into in the presence of lanthanide(iii) ions. The use of β-hairpin forming motifs as linkers enables a metal-coordinated pseudopeptide with a knotted tertiary structure to be generated. The resulting pseudopeptide knot has one of the highest backbone-to-crossing ratios (BCR)—a measure of knot tightness (a high value corresponding to looseness)—for a synthetic molecular knot to date. Preorganization in the crossing-free turn section of the knot affects aromatic stacking interactions close to the crossing region. The metal-coordinated pseudopeptide knot is compared to overhand knots with other linkers of varying tightness and turn preorganization, and the entangled architectures characterized by NMR spectroscopy, ESI-MS, CD spectroscopy and, in one case, X-ray crystallography. The results show how it is possible to program specific conformational properties into different key regions of synthetic molecular knots, opening the way to systems where knotting can be systematically incorporated into peptide-like chains through design.

Spacers linking 2,6-pyridinedicarboxamide units influence the tightness of the corresponding lanthanide-coordinated overhand knot. β-Hairpin forming motifs generate a metal-coordinated pseudopeptide with a knotted tertiary structure.

Isogonal weavings on the sphere: knots, links, polycatenanes

Mathematical knots and links are described as piecewise linear - straight, non-intersecting sticks meeting at corners. Isogonal structures have all corners related by symmetry (`vertex'-transitive). Corner- and stick-transitive structures are termed regular. No regular knots are found. Regular links are cubic or icosahedral and a complete account of these (36 in number) is given, including optimal (thickest-stick) embeddings. Stick 2-transitive isogonal structures are again cubic and icosahedral and also encompass the infinite family of torus knots and links. The major types of these structures are identified and reported with optimal embeddings. The relevance of this work to materials chemistry and biochemistry is noted.

Tying different knots in a molecular strand

The properties of knots are exploited in a range of applications, from shoelaces to the knots used for climbing, fishing and sailing¹. Although knots are found in DNA and proteins², and form randomly in other long polymer chains^3,4, methods for tying⁵ different sorts of knots in a synthetic nanoscale strand are lacking. Molecular knots of high symmetry have previously been synthesized by using non-covalent interactions to assemble and entangle molecular chains^6-15, but in such instances the template and/or strand structure intrinsically determines topology, which means that only one type of knot is usually possible. Here we show that interspersing coordination sites for different metal ions within an artificial molecular strand enables it to be tied into multiple knots. Three topoisomers-an unknot (0₁) macrocycle, a trefoil (3₁) knot^6-15, and a three-twist (5₂) knot-were each selectively prepared from the same molecular strand by using transition-metal and lanthanide ions to guide chain folding in a manner reminiscent of the action of protein chaperones¹⁶. We find that the metal-ion-induced folding can proceed with stereoinduction: in the case of one knot, a lanthanide(III)-coordinated crossing pattern formed only with a copper(I)-coordinated crossing of particular handedness. In an unanticipated finding, metal-ion coordination was also found to translocate an entanglement from one region of a knotted molecular structure to another, resulting in an increase in writhe (topological strain) in the new knotted conformation. The knot topology affects the chemical properties of the strand: whereas the tighter 5₂ knot can bind two different metal ions simultaneously, the looser 3₁ isomer can bind only either one copper(I) ion or one lutetium(III) ion. The ability to tie nanoscale chains into different knots offers opportunities to explore the modification of the structure and properties of synthetic oligomers, polymers and supramolecules.

Circuit Topology Analysis of Polymer Folding Reactions

Circuit topology is emerging as a versatile measure to classify the internal structures of folded linear polymers such as proteins and nucleic acids. The topology framework can be applied to a wide range of problems, most notably molecular folding reactions that are central to biology and molecular engineering. In this Outlook, we discuss the state-of-the art of the technology and elaborate on the opportunities and challenges that lie ahead.

Identifying knot types of polymer conformations by machine learning

We investigate the use of artificial neural networks (NNs) as an alternative tool to current analytical methods for recognizing knots in a given polymer conformation. The motivation is twofold. First, it is of interest to examine whether NNs are effective at learning the global and sequential properties that uniquely define a knot. Second, knot classification is an important and unsolved problem in mathematical and physical sciences, and NNs may provide insights into this problem. Motivated by these points, we generate millions of polymer conformations for five knot types: 0, 3_{1}, 4_{1}, 5_{1}, and 5_{2}, and we design various NN models for classification. Our best model achieves a five-class classification accuracy of above 99% on a polymer of 100 monomers. We find that the sequential modeling ability of recurrent NNs is crucial for this result, as it outperforms feed-forward NNs and successfully generalizes to differently sized conformations as well. We present our methods and suggest that deep learning may be used in specific applications of knot detection where some error is permissible. Hopefully, with further development, NNs can offer an alternative computational method for knot identification and facilitate knot research in mathematical and physical sciences.

Macrocyclic multinuclear metal complexes acting as catalysts for organic synthesis

Recent progress in homogeneous catalysis with macrocyclic multinuclear metal complexes (categories A–C) is overviewed.

Molecular knot with nine crossings: Structure and electronic properties from density functional theory computation

The electronic structure of a molecule with nine-crossing composite knots 9₇³ link denoted by the Alexander-Briggs notation (complex-1) are studied by means of theoretical methods (DFT). The most interesting feature of this kind of molecules is their capability to capture anion spices inside the cage. Stability and chemical reactivity were evaluated taking advantage of the criteria chemical hardness and chemical potential. The simulation of the infrared spectra is also included and shows the characteristic signal of the molecule in a range 1000-1600 cm^-1. The frontier molecular orbitals were also analyzed. Whereas the capability to capture chlorine ion into the cavity of the complex-1 is explored by means the analysis of bond energy. Also, the electron density distribution of the chlorine complex was studied by means the quantum theory of atoms in molecules (QTAIM) formalism in order to stablish its bonding properties as well as the electron transfer between chlorine ion and complex-1 which was approached by the natural bonding orbital (NBO) and Hirshfeld charge. Ours results revels semiconductor behaviors for both compounds.

Stereoselective Synthesis of Molecular Square and Granny Knots

We report on the stereoselective synthesis of both molecular granny and square knots through the use of lanthanide-complexed overhand knots of specific handedness as three-crossing “entanglement synthons”. The composite knots are assembled by combining two entanglement synthons (of the same chirality for a granny knot; of opposite handedness for a square knot) in three synthetic steps: first, a CuAAC reaction joins together one end of each overhand knot. Ring-closing olefin metathesis (RCM) then affords the closed-loop knot, locking the topology. This allows the lanthanide ions necessary for stabilizing the entangled conformation of the synthons to subsequently be removed. The composite knots were characterized by ¹H and ¹³C NMR spectroscopy and mass spectrometry and the chirality of the knot stereoisomers compared by circular dichroism. The synthetic strategy of combining building blocks of defined stereochemistry (here overhand knots of Λ- or Δ-handed entanglement) is reminiscent of the chiron approach of using minimalist chiral synthons in the stereoselective synthesis of molecules with multiple asymmetric centers.

Topological entanglement of interlocked knotted-unknotted polymer rings

Topological entanglements in biopolymers could drive them to certain internal statics and dynamics with important implications for biological functions. In this study, by means of molecular dynamics simulations, we demonstrate that the minimal crossing pattern of a braid plays a major role in its structural and dynamical properties; the braid consists of a knotted ring and an interlocked entwined unknotted polymer ring. In particular, we show that depending on the bending rigidity of the chains, the conformational energy of the braid can be either lower or higher than the unlocked polymer rings. Additionally, we find that a non-identical crossing pattern in the braid could distinctly enforce concerted internal conformational fluctuations between the interlocked rings.