How neocarcerand octacid4 self-assembles with guests into irreversible noncovalent complexes and what accelerates the assembly

Nature

How neocarcerand octacid4 self-assembles with guests into irreversible noncovalent complexes and what accelerates the assembly"


Play all audios:

Loading...

ABSTRACT Cram’s supramolecular capsule Octacid4 can irreversibly and noncovalently self-assemble with small-molecule guests at room temperature, but how they self-assemble and what


accelerates their assembly remain poorly understood. This article reports 81 distinct Octacid4•guest self-assembly pathways captured in unrestricted, unbiased molecular dynamics simulations.


These pathways reveal that the self-assembly was initiated by the guest interaction with the cavity portal exterior of Octacid4 to increase the portal collisions that led to the portal


expansion for guest ingress, and completed by the portal contraction caused by the guest docking inside the cavity to impede guest egress. The pathways also reveal that the self-assembly was


accelerated by engaging populated host and guest conformations for the exterior interaction to increase the portal collision frequency. These revelations may help explain why the presence


of an exterior binding site at the rim of the enzyme active site is a fundamental feature of fast enzymes such as acetylcholinesterase and why small molecules adopt local minimum


conformations when binding to proteins. Further, these revelations suggest that irreversible noncovalent complexes with fast assembly rates could be developed—by engaging populated host and


guest conformations for the exterior interactions—for materials technology, data storage and processing, molecular sensing and tagging, and drug therapy. SIMILAR CONTENT BEING VIEWED BY


OTHERS A SELF-COMPLEMENTARY MACROCYCLE BY A DUAL INTERACTION SYSTEM Article Open access 26 September 2022 MACROMOLECULAR ARCHITECTURES CONSTRUCTED BY BISCALIX[5]ARENE–[60]FULLERENE


HOST–GUEST INTERACTIONS Article 18 November 2022 COMPLEX STRUCTURES ARISING FROM THE SELF-ASSEMBLY OF A SIMPLE ORGANIC SALT Article 10 February 2021 INTRODUCTION Irreversible bimolecular


complexes are highly desirable because (1) the residence of one molecule inside the cavity of another molecule is permanent until the breakdown of the cavity-containing molecule, (2) the


complex formation is one-to-one stoichiometric, and (3) for drug therapy, the permanent residence and 1:1 stoichiometry confer the cavity-residing molecule with desired high metabolic


stability, high potency-to-mass ratio, low-acting dose, and low off-target activity1,2,3,4. These complexes generally refer to irreversible covalent complexes that are self-assembled by two


molecules between which there is a covalent bond resulting from the self-assembly process. The desirability of irreversible covalent complexes is apparent from the drug action mechanisms of


aspirin5, penicillin6, and sotorasib7, all of which involve an irreversible covalent complex. Notably, sotorasib was approved in May 2021


(https://www.fda.gov/news-events/press-announcements/fda-approves-first-targeted-therapy-lung-cancer-mutation-previously-considered-resistant-drug) as the first-in-class personalized


treatment for a lung-cancer mutation previously considered resistant to drug therapy due to sotorasib’s unique capability to clinically block the function of KRASG12C (an enzyme mutant


responsible for ~13% non-small-cell lung cancers) via irreversible complexation. However, the irreversible covalent complexes from in situ cysteine conjugation are limited by the infrequent


presence of the noncatalytic cysteine in a protein cavity, and inhibiting KRASG12C offers treatment for only a subset of cancer patients. A paradigm shift is needed for the design of


irreversible bimolecular complexes. In terms of both intrinsic binding from the thermodynamic perspective and constrictive binding from the kinetic perspective8,9,10, the irreversible


complexes include irreversible noncovalent complexes that are self-assembled by two molecules between which there is no covalent bond. Here, the intrinsic binding is the complexation


governed by intermolecular interactions between the two molecules and between the solvent and each of the two molecules, while the constrictive binding is the complexation controlled by the


thermal energy required for the guest to overcome the steric hindrance from the host during the assembly or disassembly process, as exemplified below by Cram’s supramolecular capsules known


as carcerand11, hemicarcerand12, and neocarcerand13, 14. The hallmark of the constrictive binding is the heightened energy barrier for the assembly or disassembly of a host•guest complex.


This barrier can make the dissembled or assembled molecules kinetically stable, namely, it takes a long time to assemble the two disassembled molecules or disassemble the two assembled


molecules if the barrier for the conversion is heightened. When the barrier is extremely high, assembly or disassembly requires covalent-bond making or breaking, respectively. For example, a


guest can be noncovalently trapped in the cavity of the two bowl-shaped fragments that are rim-to-rim tethered by linkers (of a host known as carcerand) when the guest is present in the


reaction medium for the tethering reaction; the disassembly of the resulting carcerand complex is not allowed unless the host is broken by excessive heat11. When the barrier is moderately


high, assembly or disassembly requires annealing. For example, a guest can enter or exit the cavity of a host known as hemicarcerand once the bimolecular system is heated to expand the


cavity portal, and the guest remains inside the cavity once the system is cooled to contract the portal; the disassembly of the hemicarcerand complex is disallowed without heating12, 15.


When the barrier is slightly high, assembly or disassembly can occur spontaneously and slowly at an ambient temperature. For some host•guest complexes, the disassembly can be disabled at the


ambient temperature by the complexation that subsequently heightens the barrier for disassembly. For example, a guest can enter the cavity of a host known as neocarcerand at room


temperature and remain in the cavity at room temperature unless the portal is opened by heating because the docking of the guest at the cavity induces a host conformational change that


consequently closes all cavity portals14. The applications of carcerand and hemicarcerand are however limited because their complexes cannot be formed in situ and adiabatically. Although


neocarcerand can form its complexes in situ and adiabatically13, 14, 16, the applications of neocarcerand complexes are also limited due to their slow complexation rates. This underscores


the need to understand how two molecules self-assemble into an irreversible noncovalent complex and what accelerates their assembly as these high-level questions hold the key to designing


irreversible noncovalent bimolecular complexes with fast complexation rates for broad applications. To promote irreversible noncovalent bimolecular complex design, this article reports 81


distinct pathways of the irreversible noncovalent self-assembly of neocarcerand Octacid4 with three known guests13—1,4-dioxane (dioxane), _p_-xylene (xylene), and naphthalene (Fig. 1 and


Table S1). These pathways were captured in multiple distinct, independent, unrestricted, unbiased, and classical isobaric–isothermal molecular dynamics (MD) simulations at a high time


resolution with an aggregated simulation time exceeding 3.761664 milliseconds at 298–363 K, rendering the structural and kinetic information needed to answer the two high-level questions and


guide the design of irreversible noncovalent complexes that can be formed in situ and adiabatically with desired kinetics for materials technology, data storage and processing, molecular


sensing and tagging, and drug therapy. RESULTS THE CHALLENGE OF CAPTURING SELF-ASSEMBLY PATHWAYS In view of the current state of computational work on guest/ligand-binding


pathways17,18,19,20, unrestricted and unbiased MD simulations of the self-assembly of Octacid4 with its known guests are challenging to perform because the complexation times of


Octacid4•guest were estimated to be a few minutes (or hours for the bulkiest guest naphthalene) according to NMR experiments13. These complexation times are many orders of magnitude longer


than current MD simulation times that are on the order of milliseconds. To demonstrate the challenge of capturing the Octacid4•guest self-assembly pathways, multiple distinct, independent,


unrestricted, unbiased, and classical isobaric–isothermal MD simulations (hereafter shortened to simulations) were performed using the fully deprotonated, _apo_ Octacid4 that was surrounded


by eight neutralizing sodium cations, 10 xylenes to mimic the use of the guest in 10 equivalent excess in the NMR experiments13, 60 NaCl molecules to approximate the ionic strength of the


experimental conditions13, and 2359 TIP3P water molecules21 to mimic the experimental aqueous solution13. Indeed, no autonomous complexation was observed in any set of 20 14,251.6-ns


simulations at 298, 340, 363, or 370 K (Table S1) because the simulation time was many orders of magnitude shorter than the experimentally estimated complexation time for xylene. USING THE


PHASE-TRANSFER CATALYST TO CAPTURE THE SELF-ASSEMBLY According to the reported NMR experiments13, the experimentally observed self-assembly of Octacid4 with dioxane, xylene, or naphthalene


in a sodium borate buffer at pH 9 is mechanistically driven by the sodium cation as a phase-transfer catalyst22 that chelates with the carboxylates on the Octacid4 surface and consequently


accumulates the guest on the host surface via the cation–π interaction23 or the sodium chelation. The guest accumulation on the host surface is similar to immersing the sodium-chelated


Octacid4 in a neat guest solution. Simulating the latter can substantially accelerate the self-assembly according to the law of mass action, enabling determination of relative complexation


times of different guests with Octacid4 for mechanistic insights into the self-assembly. Simulating the latter also allows the use of linear-regression analysis to examine the convergency


and internal consistency of the MD simulations. A goodness of fit (_r_2) of <0.70 for the natural logarithm plot of the host population versus the simulation time indicates problematic


simulations. This is because the first-order rate (viz., the exponential decay of the host population over the simulation time) is expected for simulating the self-assembly of Octacid4 with


its guest in large excess. Accordingly, 40 simulations were performed using the fully deprotonated, _apo_ Octacid4 that was surrounded by eight neutralizing sodium cations and 150 xylenes.


Here, the number of xylenes was arbitrarily chosen. Gratifyingly, all 40 14,251.6-ns simulations with the fully deprotonated, _apo_ Octacid4 with 150 xylenes captured the self-assembly of


Octacid4 with xylene at 298 K. Additional simulations were performed using each of the five variations: (1) replacing the fully deprotonated, _apo_ Octacid4 with the octa-anionic Octacid4


possessing five water molecules inside the cavity, (2) increasing the number of xylenes to 250, (3) replacing 150 xylenes with 150 dioxanes, (4) replacing 150 xylenes with 150 naphthalenes,


or (5) changing the Berendsen thermostat to the Langevin thermostat (Table S1). All these simulations captured the self-assembly event. For the simulations with the water-bound Octacid4, all


water molecules inside the cavity had a high propensity to interact with the sodium cations and the carboxylates outside the cavity, and egression of all five water molecules occurred prior


to the ingression of xylene. By contrast, no complexation was observed under the same conditions if the fully deprotonated, _apo_ Octacid4 surrounded by 8 sodium cations was replaced with


the fully protonated, _apo_ Octacid4 without any cations (Table S1). All simulations described hereafter used the fully deprotonated, _apo_ Octacid4 surrounded with 8 sodium cations, 150


guests, and the Berendsen thermostat. CONFORMATIONAL CHARACTERIZATION OF THE OCTACID4•XYLENE SELF-ASSEMBLY PATHWAYS Visual inspections of the 40 self-assembly pathways of Octacid4•xylene at


298 K revealed that all pathways comprised three common steps (Fig. 2 and Videos S1–S6 and S11–S16). In Step 1, xylene entered the cavity portal exterior space that was confined by two


aromatic linkers of the host; as apparent from the representative videos, there were substantial collisions between xylene and the equatorial portal in the linker region of Octacid4 that led


to the portal expansion for guest ingress. In Step 2, xylene passed one methyl group through the equatorial portal, then the phenyl group through the portal with the phenyl plane


perpendicular to the axial axis (viz., the axis passing two axial portals in the bowl-shaped region), and last the other methyl group through the equatorial portal. In Step 3, xylene rotated


~90° to keep its phenyl plane parallel to the axial axis. This rotation caused the equatorial portal contraction and impeded guest egress, according to the reported structural analyses of


_apo_ Octacid4 and the Octacid4•xylene complex14. KINETIC CHARACTERIZATION OF THE OCTACID4•XYLENE SELF-ASSEMBLY PATHWAYS A survival analysis of the 40 14,251.6-ns simulations that all


captured the self-assembly of Octacid4•xylene at 298 K showed xylene’s mean complexation time to be 1022 ns (95% confidence interval: 750–1392 ns). Here the complexation time was defined as


the first time instant at which xylene was inside the host cavity and had its phenyl plane parallel to the axial axis of the host. This xylene orientation was found in the most populated


Octacid4•xylene conformation14. To dissect the self-assembly kinetics, the duration of the first step is herein termed priming time, and the duration of the last two steps is termed


ingression time. These names are used because during the first step, the host and guest conformations are primed, through self-selection and conformational rearrangements, for guest ingress,


and because during the last two steps, the guest enters the host cavity and rotates ~90° to form a complex. The priming time was defined as a time period from the beginning of the MD


simulation to the last time instant at which the distance between any guest hydrogen atoms and any host methylene hydrogen atoms was greater than 2.6 Å. The hydrogen–hydrogen distance cutoff


(abbreviated as HH cutoff) was set at 2.6 Å for the following reasons. Visual inspection of all 40 Octacid4•xylene pathways revealed that xylene was outside of the Octacid4 cavity as long


as the center-of-mass distance cutoff (abbreviated as COM cutoff) for Octacid4 and its guest was ≥7 Å. However, xylene could have its terminal methyl group contact the host cavity portal to


slightly enter the host cavity portal in some of the 40 Octacid4•xylene pathways at the COM cutoff of 7 or 8 Å. Xylene could also be relatively away from the portal in some pathways at the


COM cutoff of ≥8 Å. Therefore, rather than using the COM cutoff, the HH cutoff of 2.6 Å was used to avoid both the methyl group contacting the portal (which consequently shortens the


ingression time) and the guest being away from the portal (which consequently lengthens the ingression time). The ingression time was defined as a time period from the last time instant at


which the distance between any guest hydrogen atoms and any host methylene hydrogen atoms was greater than 2.6 Å to the first time instant at which the guest rotated ~90° inside the cavity.


The complexation time is now a sum of the priming and ingression times. This dissection reveals that the ingression time (8–934 ps; Table S2) is a fraction of the priming time (12–14,199 ns;


Table S2), and hence the contribution of the ingression time to the complexation time is inconsequential. In other words, the guest ingress is so fast that the complexation time is


determined primarily by the priming time. This conclusion is independent of the use of the HH cutoff of 2.6 Å because changes of this cutoff only slightly affect the ingression time that is


several orders of magnitude shorter than the priming time. For example, as apparent from Table S2, the priming and complexation times determined from the HH cutoff of 2.6 Å are identical to


those determined from the COM cutoff of 10 Å, and the average ingression times from the HH cutoff versus the COM cutoff are 95 versus 105 ps. RELATIVE COMPLEXATION TIMES OF DIFFERENT GUESTS


WITH OCTACID4 Multiple simulations were performed for dioxane or naphthalene under the same simulation conditions as those for xylene. These simulations showed that the self-assembly of


Octacid4 with dioxane or naphthalene was much slower than that of Octacid4•xylene according to the priming time defined using either the HH cutoff of 2.6 Å or the COM cutoff of 8 Å for


dioxane and 10 Å for naphthalene (Table S2). The use of the HH cutoff was according to the Octacid4•guest-pathway analysis, which revealed that dioxane and naphthalene were close enough


(without slightly entering the host cavity portal) to the portal at the HH cutoff of 2.6 Å. The use of two different COM cutoffs was because, unlike dioxane that was close enough to the


portal at the COM cutoff of 8 Å, naphthalene could have its β-carbon atoms contact the portal at the COM cutoff of 8 Å or be relatively away from the portal in some pathways at the COM


cutoff of ≥8 Å. Consistent with the complexation kinetics of xylene described above, the priming and complexation times of dioxane and naphthalene determined from the HH cutoff are also


identical to those determined from the COM cutoff (Table S2), and the average ingression times determined from the HH cutoff versus COM cutoff are 29/18 versus 29/17 ps for dioxane at


298/340 K and 88/670/7790 versus 88/668/7790 ps for naphthalene at 298/340/363 K (Table S2). For dioxane, 100 6320-ns simulations captured two self-assembly pathways at 298 K (priming and


ingression times using the HH cutoff of 2.6 Å: 3668 and 5681 ns and 6 and 52 ps, respectively; Table S2), and 40 12,640-ns simulations at 340 K captured 24 self-assembly pathways (priming


and ingression times using the HH cutoff of 2.6 Å: 34–2816 ns, 6–50 ps; Table S2). For naphthalene, 100 6320-ns simulations at 298 K captured one self-assembly pathway (priming and


ingression times using the HH cutoff of 2.6 Å: 5721 ns and 88 ps, respectively; Table S2), 100 6320-ns simulations at 340 K captured four self-assembly pathways (priming and ingression times


using the HH cutoff of 2.6 Å: 1865–4389 ns, and 16–2189 ps, respectively; Table S2), and 100 7900-ns simulations at 363 K captured ten self-assembly pathways (priming and ingression times


using the HH cutoff of 2.6 Å: 651–7234 ns and 51–31,035 ps, respectively; Table S2). The ingression times of naphthalene were substantially longer than those of xylene and dioxane, but these


longer ingression times of naphthalene were still a small portion of the naphthalene priming times, confirming that the ingression time is so short that the complexation time is governed


largely by the priming time. All pathways of dioxane and naphthalene shared the three common steps of the Octacid4•xylene pathways (Figs. 3 and 4), except for more profound collisions of


dioxane or naphthalene with the equatorial portal than those of xylene in the first step (Videos S7–S10 and S17–S18) and subtle differences in the second step noted as follows. For dioxane,


the second step involved the passing of the dioxane oxygen through the equatorial portal during which the oxygen was in the energetically less stable half-chair conformation, then the


four-methylene portion of dioxane through the portal during which the methylene portion was in the energetically stable chair conformation, and last the other dioxane oxygen through the


portal during which the oxygen was again in the half-chair conformation (Figs. 3 and 5 and Videos S7 and S17). Interestingly, the Octacid4•dioxane self-assembly process followed the mutually


induced fitting mechanism that is akin to the mechanism for a reported synthetic complex24. For naphthalene, the passing of the two head β-carbon atoms, then the α-carbon portion, and


finally the tail β-carbon atoms of the guest through the host portal was completed at the second step (Fig. 4 and Videos S9–S10). EFFECT OF CONFORMATIONAL STABILITY ON COMPLEXATION TIME To


understand why the complexation time of xylene is much shorter than those of dioxane and naphthalene, conformational analyses of the 81 pathways were performed and revealed the involvement


of three clusters of the Octacid4 conformations during the Octacid4•xylene self-assembly process at 298 K (Fig. 6). The most-populated cluster (population: 21/40) had two nearly orthogonal


aromatic linkers that strongly interacted with xylene according to the intermolecular interactions depicted by the noncovalent interaction gradient isosurfaces25 (Fig. 6a and Videos S1–S2


and S11–S12), and the assembly involving this cluster was captured mainly at the early stage of the simulations. The less-populated cluster (population: 12/40) had two nearly parallel,


face-to-face aromatic linkers that moderately interacted with xylene (Fig. 6b and Videos S3–S4 and S13–S14), and the complexation involving this cluster was captured at the intermediate


stage. The least-populated cluster (population: 7/40) had two nearly-coplanar aromatic linkers that weakly interacted with xylene (Fig. 6c and Videos S5–S6 and S15–S16), and the complexation


involving this cluster was captured at the late stage. In contrast to xylene that had the attraction from the two nearly orthogonal aromatic linkers in its top-5 fastest pathways at 298 K


(complexation times: 12–41 ns; Table S2; Figs. 2 and 6a), dioxane and naphthalene had the attraction from the two nearly parallel aromatic linkers in their fastest pathways at 298 K


(complexation times: 3668 and 5729 ns; Table S2; Figs. 3 and 4). Consistent with the nature of π–π interactions26, the Octacid4 conformation with two nearly orthogonal aromatic linkers


corresponded to the 5th most populated conformation of _apo_ Octacid4 in water, but the one with two nearly parallel, face-to-face aromatic linkers corresponded to none of the top-10 most


populated conformations of the aqueous _apo_ Octacid4. These results demonstrate the effect of the host conformational stability on complexation time. More importantly, the results reveal


that the self-assembly of Octacid4 with its guest is governed by the conformational complementarity between the two molecules not only during the ingression time but also during the priming


time, and that the host or guest molecule can adopt, at a cost of lengthening the complexation time, an unpopulated conformation to proceed the self-assembly along a repulsion-free path when


its populated conformation is not complementary to the conformation of its partner. OPPOSING EFFECTS OF THE SODIUM CATION ON COMPLEXATION TIME According to the conformational analysis of


the 43 self-assembly pathways of Octacid4 with dioxane (2 pathways), xylene (40 pathways), and naphthalene (1 pathway) at 298 K, the sodium cation that chelated with the Octacid4 carboxylate


either coordinated with the guest oxygen atom or formed the cation–π interaction with the guest aromatic ring. However, the chelation with the oxygen atom of the guest that subsequently


entered the host cavity was observed in one of the two pathways of dioxane (Table S2); the transient cation–π interaction with the aromatic ring of the cavity-entering guest was observed in


only 20 of the 41 self-assembly pathways of xylene and naphthalene (Table S2). These observations indicate that the cavity-entering guest either does not interact at all or does not strongly


interact with the carboxylate-chelated sodium cation, so that the cavity-entering guest is not trapped at the host linker region. Instead, the cavity-entering guest interacts with an


immobilized guest, which is trapped in the linker region due to its strong interaction with the carboxylate-chelated sodium cation, via the π–π interaction for xylene or naphthalene or the


van der Waals interaction for dioxane. Because the π–π interaction26 and the van der Waals interaction are generally weaker than the cation–π interaction23, these weak interactions enable


the immobilized guest to usher the cavity-entering guest into the cavity, revealing the role of the sodium cation in shortening the complexation time of Octacid4 with dioxane, xylene, and


naphthalene. The conformational analysis of the 43 pathways also identified two small but interesting clusters of Octacid4 conformations that were derived a priori from the MD simulations


(Fig. 7 and Videos S11–S18). In one cluster with a population of 11/43 (Table S2) that was associated mainly with the host conformations with two nearly parallel linkers, the carboxylates


from two nearby linkers of Octacid4 formed a bidentate coordination with a sodium cation (Fig. 7a and Videos S11–S14 and S17–S18). This bidentate coordination rigidified a pair of the


bidendate linkers and blocked one of the four equatorial portals of the host. In the other cluster with a population of 3/43 (Table S2) that was associated mainly with the host conformations


with two nearly coplanar linkers, all four linkers of Octacid4 were involved in the bidentate coordination, resulting in rigidification of two pairs of the bidendate linkers and blockage of


two of the four equatorial portals (Fig. 7b and Videos S15–S16). According to the energy minimization and frequency calculations of the two simulation-derived conformations with one or two


pairs of the bidendate linkers using the Gaussian 16 program and HF/6-31 G* or B3LYP/6-31 G*, the minimization-derived conformations closely resembled the simulation-derived conformations


(heavy-atom root-mean-square deviations for one or two pairs of the bidendate linkers: 0.75 or 0.32 Å for HF/6-31 G* and 0.75 or 0.34 Å for B3LYP/6-31 G*), respectively, and no imaginary


frequencies were found for the minimization-derived conformations, suggesting that the conformations with one or two pairs of the bidendate linkers derived a priori from the simulations were


local minimal conformations. Because the host conformations with the nearly parallel or nearly coplanar linkers were associated with long complexation times as described above, these


results indicate that the sodium cation also played a role in lengthening the complexation time of Octacid4 with dioxane, xylene, and naphthalene and that the lengthening role is outweighed


by the shortening role of the sodium cation as the populations of the host conformations with the bidendate linkers are lower than those of the host conformations without the bidendate


linkers (Table S2). DISCUSSION INTERNAL AND EXTERNAL CONSISTENCIES OF THE SELF-ASSEMBLY PATHWAYS The 81 self-assembly pathways described above exhibited internal or external consistencies as


follows. (1) The natural logarithm of the Octacid4 population versus the simulation time exhibited a linear relationship with _r_2 of 0.98 for all 40 simulations that captured the


Octacid4•xylene self-assembly, 0.85 for the first 20 of the 40 simulations, and 0.96 for the last 20 of the 40 simulations (Fig. 8). These _r_2 values indicate the convergency and internal


consistency of the 40 simulations for which the first-order self-assembly rate is expected as explained above. (2) As described above, the simulations using the fully deprotonated Octacid4


under various conditions all captured the self-assembly of Octacid4 with dioxane, xylene, and naphthalene, but no autonomous complexation was observed under the same conditions if the fully


deprotonated host was replaced by the fully protonated host. These results are consistent with the use of the Octacid4 solution containing the sodium borate buffer at pH 9 to detect the


Octacid4•guest complexation in the NMR experiments13 and with the phase-transfer catalysis undergirding those NMR experiments. (3) One key finding of the present work is that the guest


ingress is so fast that complexation time is determined primarily by the priming time. This finding is consistent with the report that many dense-phase reactions can be considered as gated


reactions in that the rate of a local reaction is governed largely by the initial formation of a permissive atomic arrangement (viz., determined mainly by the system priming) within which


the local transformation can proceed relatively rapidly27. (4) Conformational analyses of all 43 pathways of xylene, dioxane, and naphthalene at 298 K showed that Octacid4 adopted


exclusively a cluster of V-shaped conformations to gulp its guests. These V-shaped conformations have the mean C17a–C6b distance (the distance between two diphenoxymethane carbon atoms that


control the width of the cavity portal as shown in Fig. 1) of 7.0 Å (95% confidence interval: 7.0–7.1 Å) for the entrance portal and the corresponding mean C17a–C6b distance of 5.6 Å (95%


confidence interval: 5.5–5.7 Å) for the opposing portal (Table S3). These mean distances are consistent with the reported V-shaped conformation proposed for the sliding-door mechanism for


the gating of hemicarcerands15 that are closely related to Octacid4. (5) For the self-assembly at 298 K, 12 of the 40 captured pathways for xylene and all captured pathways for dioxane and


naphthalene had two nearly parallel linkers that channeled the guest into the cavity (Table S2). These linker channels are consistent with the report that an antechamber formed by two


parallel linkers of a hemicarcerand played a role in the gating of hemicarcerands15. (6) Relative to the noncovalent interaction gradient isosurfaces25 of dioxane that showed the stronger


intramolecular repulsion in the half-chair conformation than that in the chair conformation (Fig. 5), the isosurfaces of the Octacid4 complexes in the 43 self-assembly pathways of xylene


(40), dioxane (2), and naphthalene (1) at 298 K revealed no repulsion (or only a trace of repulsion for naphthalene) but increasingly strong attraction between Octacid4 and its guest in the


self-assembly process (Figs. 2–4). These observations are consistent with the report that the Corey–Pauling–Koltun space-filling models of xylene, dioxane, and naphthalene could be pushed


(or pushed with effort for naphthalene) through portals of the space-filling model of Octacid413, indicating that a repulsion-free path exists for Octacid4 to self-assemble with xylene,


dioxane, and naphthalene. (7) Most importantly, the capturing of the self-assembly of Octacid4 with dioxane, xylene, and naphthalene in the simulations at 298 K is consistent with the


Octacid4 complexation with these guests detected in the NMR experiments at 298 K13. HOW OCTACID4 SELF-ASSEMBLES WITH ITS GUESTS AND WHAT ACCELERATES THEIR ASSEMBLY Collectively, the 81


pathways indicate that the irreversible noncovalent self-assembly process of Octacid4 with dioxane, xylene, and naphthalene was (1) initiated by the guest interaction with the cavity portal


exterior of the host to increase the collisions of the portal that led to the portal expansion for guest ingress, (2) accelerated by engaging populated conformations of host and guest for


the exterior intermolecular interaction to increase the portal collision frequency, and (3) completed by the portal contraction caused by the guest docking inside the cavity to impede guest


egress. This type of self-assembly is a process of stepwise conformational rearrangements of the two molecules—kinetically starting from the pre-assembly state (at which the guest interacts


with the exterior of the host) to the assembly and post-assembly states (at which the guest interacts partially and fully with the interior of the host, respectively)—to gradually increase


the attraction between the two molecules along a repulsion–free path until the maximal attraction is reached. The kinetics of this process was governed primarily by the complementarity


between the two molecules during the priming time (on the order of microseconds or longer) rather than the ingression time (on the order of picoseconds to nanoseconds). HETEROGENEITY OF THE


OCTACID4•XYLENE SELF-ASSEMBLY PATHWAYS The 40 Octacid4•xylene self-assembly pathways derived from the converged MD simulations (as evident from Fig. 8) demonstrate the heterogeneity of the


self-assembly pathways. As apparent from Table S2, at 298 K, xylene can enter the cavity of the three distinct clusters of the host conformations shown in Fig. 6 and Videos S1–S6. In each of


these clusters, xylene can also enter the cavity of the host conformations with or without bidendate linkers (Videos S1–S6 and S11–S16), indicating that the Octacid4•xylene pathways can be


shortened or lengthened by the sodium cation (Table S2). This heterogeneity is akin to the heterogeneity of protein-folding pathways28. It points to a limitation of the simulation protocol


used in this work because the protocol captured only the fast self-assembly pathways of Octacid4•dioxane and Octacid4•naphthalene at 298 K. This calls for the development of a new simulation


protocol to avoid oversimplification of the self-assembly process such as equating the process to a few fast self-assembly pathways for Octacid4•dioxane and Octacid4•naphthalene. As


apparent from the Octacid4•xylene self-assembly pathways, it is not a few complexation times (from the fast pathways) but the mean complexation time (from the fast, intermediate, and slow


pathways) that offers insight into the delicate balance between the sodium cation’s roles in lengthening and shortening the complexation time of Octacid4•xylene. Given the known theoretical


work on gated reactions27 and the finding of the present work that the ingression time is so short that the complexation time is governed largely by the priming time, it is worth noting the


need to avoid overcomplication of the self-assembly process, such as simulating the Octacid4•guest self-assembly using the octa-anionic Octacid4 possessing five water molecules inside the


cavity or a more complicated host with a variable number of water molecules inside the cavity, for at least two reasons. First, it has not been determined experimentally whether the


egression of water molecules that was observed in the present work will actually occur under the NMR experiment conditions once the water-bound Octacid4 is surrounded by a layer of


water-insoluble xylene or naphthalene. Second, the water molecule is much smaller than dioxane, xylene, and naphthalene. The ingression and egression times of water are hence much shorter


than those of the three guests, and binding of water molecules to the host cavity is likely opportunistic rather than intrinsic or constrictive. The inclusion of water molecules inside the


cavity will affect the Octacid4•guest self-assembly but not likely to the extent that the mean complexation time will be significantly altered. DESIGN OF IRREVERSIBLE NONCOVALENT COMPLEXES


WITH FAST ASSEMBLY RATES Despite the heterogeneity issue as described above, the realization that the self-assembly initiation by the exterior intermolecular interaction and the control of


the complexation time largely by the priming time is an advance and may help explain why the presence of a cryptic exterior site for substrate or inhibitor binding at the rim of the


enzyme-active site is a fundamental feature of fast enzymes such as acetylcholinesterase29,30,31,32,33,34,35,36 and why it is advantageous for the proapoptotic multidomain BAK protein to


have a noncanonical BH3-binding groove37 that abuts the canonical one for constrictive complexation of BAK—due to the steric hindrance from R88 and Y89 on the edge of the canonical


groove38—with proapoptotic BH3-only proteins. The realization that the self-assembly acceleration by the adoption of populated conformations for the exterior intermolecular interaction is an


advance and may help explain why small molecules reportedly prefer to adopt local minimum conformations when binding to proteins39. Further, these realizations suggest that irreversible


noncovalent complexes with fast assembly rates could be designed by accounting for the complementarity between guest and host both of which adopt populated conformations for the exterior


intermolecular interactions during the priming time. This design strategy may facilitate a paradigm shift from irreversible covalent complex design to irreversible noncovalent complex design


for materials technology, data storage and processing, molecular sensing and tagging, and drug therapy—especially the personalized drug therapies for cancer patients with somatic mutations


other than the KRASG12C mutation7. METHODS MOLECULAR DYNAMICS SIMULATION The fully deprotonated, _apo_ Octacid4 neutralized with 8 sodium ions (or an Octacid4 in a different configuration,


such as the fully protonated Octacid4, as listed in Table S1) was manually solvated with 150/250 copies of a guest (xylene, dioxane, or naphthalene) using PyMOL V1.7.0.3 (https://pymol.org)


and tLEaP of the AmberTools 16 package (University of California, San Francisco) and then energy-minimized for 100 cycles of steepest-descent minimization followed by 900 cycles of


conjugate-gradient minimization to remove close van der Waals contacts using SANDER of the AMBER 11 package (University of California, San Francisco), FF12MClm40, and a cutoff of 8.0 Å for


noncovalent interactions. The tLEaP input file for building the fully protonated Octacid4 and the Cartesian coordinates of the energy-minimized Octacid4•guest (in all configurations as


listed in Table S1) are provided in Data S1 and S2, respectively. The energy-minimized system was slowly heated to 298/340/363/370 K in 30 steps under constant temperature and constant


volume, and then equilibrated for 106 timesteps under constant temperature of 298/340/363/370 K and constant pressure of 1 atm employing isotropic molecule-based scaling. Finally, a set of


20/40/100 distinct, independent, unrestricted, unbiased, and classical isobaric–isothermal MD simulations was performed for the resulting system using PMEMD of the AMBER 14/16/18/20 package


(University of California, San Francisco), FF12MClm40, and a periodic boundary condition at 1 atm and 298/340/363/370 K. All simulations used (i) a dielectric constant of 1.0, (ii) the


Berendsen coupling algorithm41 for thermostat and barostat, (iii) the particle mesh Ewald method42 to calculate electrostatic interactions of two atoms at a separation of >8 Å, (iv) Δ_t_ 


= 1.00 fs of the standard-mass time40, 43, (v) the SHAKE-bond-length constraint applied to all bonds involving hydrogen, (vi) a protocol to save the image closest to the middle of the


“primary box” to the restart and trajectory files, (vii) a formatted restart file, (viii) the revised alkali-ion parameters44, (ix) a cutoff of 8.0 Å for noncovalent interactions, (x) a


uniform 10-fold reduction in the atomic masses of the entire simulation system (both solute and solvent)40, 43, (xi) NTWX = 100 steps for coordinates’ output, and (xii) default values of all


other inputs of PMEMD. Available in the Supporting Information of Ref. 43, FF12MClm is a revised AMBER protein forcefield with no parameterization for any Octacid4•guest complexes40. This


forcefield is able to (1) capture the experimentally observed exponential decay of the non-native state population of fast-folding proteins over simulation time with _r_2 > 0.90 and (2)


fold these proteins with agreements between simulated and experimental folding times within factors of 0.6–1.445. FF12MClm was used in this study to investigate the noncovalent self-assembly


of small-molecule guests with Octacid4 whose aromatic linkers can flip between left‐ and right‐handed configurations and usher the guest into the host cavity. This was because of the


effectiveness of FF12MClm in simulating the experimentally observed flipping between left‐ and right‐handed configurations for C14–C38 of bovine pancreatic trypsin inhibitor in solution40


and because of the need to compress the simulation time (viz., speed up simulations) by a factor of 101/2 through 10-fold uniform reduction of the system mass. While the hydrogen mass


repartitioning scheme can also speed up simulations, it was not used in this study because it would affect dynamic properties of the system46. The forcefield parameters for the fully


deprotonated Octacid4, dioxane, xylene, and naphthalene are available in the Supplementary Information of Ref. 14. The forcefield parameters for the neutral Octacid4 were developed using a


published procedure14 and provided in Data S1. The ab initio calculations for developing the fully protonated Octacid4 forcefield parameters were performed using Gaussian 98 (Revision A.7;


Gaussian, Inc. Wallingford, CT). The chemical structures of HC2 and HCD2 and the assembly of HC2 into the fully protonated Octacid4 are shown in Fig. S1. All MD simulations were performed


using a dedicated cluster of 100 12-core Apple Mac Pros with Intel Westmere (2.40/2.93 GHz) and computers at the University of Minnesota Supercomputing Institute and the Mayo Clinic


high-performance computing facility at the University of Illinois Urbana-Champaign National Center for Supercomputing Applications. SURVIVAL ANALYSIS The mean complexation time and its 95%


confidence interval for the Octacid4•xylene self-assembly at 298 K was obtained from the 40 individual complexation times of xylene listed in Table S2 using the parametric survival function


[the Surreg() function] implemented in the R survival package Version 3.2.047. NONCOVALENT INTERACTION GRADIENT ISOSURFACE All noncovalent interaction gradient isosurfaces were generated


using the NCIPLOT program (Version 4)48 with keyword RANGE (3, –0.1 to –0.02, –0.02 to 0.02, and 0.02–0.1 au) for Fig. 5 or keywords LIGAND (4.0 Å) and RANGE (3, –0.1 to –0.02, –0.02 to


0.02, and 0.02–0.1 au) for Figs. 2–4 and 6. GAUSSIAN ENERGY MINIMIZATION AND FREQUENCY CALCULATIONS All energy minimization and frequency calculations were performed using Gaussian 16


(Revision C.01; Gaussian, Inc. Wallingford, CT) and HF/6-31 G* or B3LYP/6-31 G* on computers from the Mayo Clinic high-performance computing facility at the University of Illinois


Urbana-Champaign National Center for Supercomputing Applications. Each frequency calculation was done using %mem = 40 Gb and %nproc = 16. The Cartesian coordinates of the Octacid4


conformations (with one or two pairs of bidendate linkers) before and after the minimization are provided in Data S3. VIDEO The time series of the Cartesian coordinates of all heavy atoms of


Octacid4 and the guest (that was initially outside of the Octacid4 cavity and then entered and remained inside the cavity) was extracted from the output file of an MD simulation and saved


to a concatenated PDB file using the CPPTRAJ of the AMBER 14/16/18 package. The guest involved in the complexation was identified through the analysis of eight intermolecular distances


between a guest atom (C3 for xylene, O1 for dioxane, and C1 for naphthalene) and the eight diphenoxymethane carbon atoms of Octacid4. The eight distances for each of 150 guests were


calculated from the output file using the CPPTRAJ, and the guest was considered (confirmed by visual inspection) to be inside the Octacid4 cavity if all eight distances were <9.0 Å. The


concatenated PDB file was then converted to a set of PNG files using the PyMOL V1.7.0.3 showing the top or side view of the bimolecular complex with the guest in the stick-and-ball model and


the host in the stick model. Before the PDB-to-PNG conversion, each complex conformation was superimposed over the proceeding one using the PyMOL. The resulting PNG files were lastly


converted to a video file using Adobe Photoshop 2021 (https://www.adobe.com). DATA AVAILABILITY Fig. S1, Tables S1–S3, Videos S1–S18, and Data S1–S3 are provided in the Supplementary


Information. Other relevant data are available from the corresponding author upon reasonable request. REFERENCES * Smith, A. J., Zhang, X., Leach, A. G. & Houk, K. N. Beyond picomolar


affinities: quantitative aspects of noncovalent and covalent binding of drugs to proteins. _J. Med. Chem._ 52, 225–233 (2009). Article  CAS  PubMed  PubMed Central  Google Scholar  * Singh,


J., Petter, R. C., Baillie, T. A. & Whitty, A. The resurgence of covalent drugs. _Nat. Rev. Drug Discov._ 10, 307–317 (2011). Article  CAS  PubMed  Google Scholar  * Pang, Y.-P. Insect


acetylcholinesterase as a target for effective and environmentally safe insecticides. _Adv. Insect Physiol._ 46, 435–494 (2014). Article  Google Scholar  * Vita, E. 10 years into the


resurgence of covalent drugs. _Future Med. Chem._ 13, 193–210 (2021). Article  CAS  PubMed  Google Scholar  * Vane, J. R. & Botting, R. M. The mechanism of action of aspirin. _Thromb.


Res._ 110, 255–258 (2003). Article  CAS  PubMed  Google Scholar  * Yocum, R. R., Rasmussen, J. R. & Strominger, J. L. The mechanism of action of penicillin. Penicillin acylates the


active site of Bacillus stearothermophilus D-alanine carboxypeptidase. _J. Biol. Chem._ 255, 3977–3986 (1980). Article  CAS  PubMed  Google Scholar  * Canon, J. et al. The clinical


KRAS(G12C) inhibitor AMG 510 drives anti-tumour immunity. _Nature_ 575, 217–223 (2019). Article  CAS  PubMed  Google Scholar  * Quan, M. L. C. & Cram, D. J. Constrictive binding of large


guests by a hemicarcerand containing 4 portals. _J. Am. Chem. Soc._ 113, 2754–2755 (1991). Article  CAS  Google Scholar  * Warmuth, R. & Yoon, J. Recent highlights in hemicarcerand


chemistry. _Acc. Chem. Res._ 34, 95–105 (2001). Article  CAS  PubMed  Google Scholar  * Liu, F., Helgeson, R. C. & Houk, K. N. Building on Cram’s legacy: stimulated gating in


hemicarcerands. _Acc. Chem. Res._ 47, 2168–2176 (2014). Article  CAS  PubMed  PubMed Central  Google Scholar  * Cram, D. J., Karbach, S., Kim, Y. H., Baczynskyj, L. & Kalleymeyn, G. W.


Shell closure of 2 cavitands forms carcerand complexes with components of the medium as permanent guests. _J. Am. Chem. Soc._ 107, 2575–2576 (1985). Article  CAS  Google Scholar  * Tanner,


M. E., Knobler, C. B. & Cram, D. J. Hemicarcerands permit entrance to and egress from their inside phases with high structural recognition and activation free-energies. _J. Am. Chem.


Soc._ 112, 1659–1660 (1990). Article  CAS  Google Scholar  * Yoon, J. Y. & Cram, D. J. The first water-soluble hermicarceplexes. _Chem. Commun_., 10.1039/A607353K (1997). * McFerrin, K.


G. & Pang, Y.-P. How the water-soluble hemicarcerand incarcerates guests at room temperature decoded with modular simulations. _Commun. Chem_. 4, 26 (2021). * Houk, K. N., Nakamura, K.,


Sheu, C. M. & Keating, A. E. Gating as a control element in constrictive binding and guest release by hemicarcerands. _Science_ 273, 627–629 (1996). Article  CAS  PubMed  Google Scholar


  * Jankowska, K. I., Pagba, C. V., Chekler, E. L. P., Deshayes, K. & Piotrowiak, P. Electrostatic docking of a supramolecular host-guest assembly to cytochrome c probed by bidirectional


photoinduced electron transfer. _J. Am. Chem. Soc._ 132, 16423–16431 (2010). Article  CAS  PubMed  Google Scholar  * Dixon, T., Lotz, S. D. & Dickson, A. Predicting ligand binding


affinity using on- and off-rates for the SAMPL6 SAMPLing challenge. _J. Comput. Aided Mol. Design_ 32, 1001–1012 (2018). Article  CAS  Google Scholar  * You, W. L., Tang, Z. Y. & Chang,


C. E. A. Potential mean force from umbrella sampling simulations: what can we learn and what is missed? _J. Chem. Theory Comput._ 15, 2433–2443 (2019). Article  CAS  PubMed  PubMed Central 


Google Scholar  * Miao, Y. L., Bhattarai, A. & Wang, J. N. Ligand Gaussian accelerated molecular dynamics (LiGaMD): characterization of ligand binding thermodynamics and kinetics. _J.


Chem. Theory Comput._ 16, 5526–5547 (2020). Article  CAS  PubMed  PubMed Central  Google Scholar  * Limongelli, V. Ligand binding free energy and kinetics calculation in 2020. _WIREs Comput.


Mol. Sci_. 10, e1455 (2020). * Jorgensen, W. L., Chandreskhar, J., Madura, J. D., Impey, R. W. & Klein, M. L. Comparison of simple potential functions for simulating liquid water. _J.


Chem. Phys._ 79, 926–935 (1983). Article  CAS  Google Scholar  * Starks, C. M. Phase-transfer catalysis. I. Heterogeneous reactions involving anion transfer by quaternary ammonium and


phosphonium salts. _J. Am. Chem. Soc._ 93, 195–199 (1971). Article  CAS  Google Scholar  * Ma, J. C. & Dougherty, D. A. The cation-pi interaction. _Chem. Rev._ 97, 1303–1324 (1997).


Article  CAS  PubMed  Google Scholar  * Sawada, T., Hisada, H. & Fujita, M. Mutual induced fit in a synthetic host-guest system. _J. Am. Chem. Soc._ 136, 4449–4451 (2014). Article  CAS 


PubMed  Google Scholar  * Johnson, E. R. et al. Revealing noncovalent interactions. _J. Am. Chem. Soc._ 132, 6498–6506 (2010). Article  CAS  PubMed  PubMed Central  Google Scholar  * Hunter,


C. A. & Sanders, J. K. M. The nature of π-π interactions. _J. Am. Chem. Soc._ 112, 5525–5534 (1990). Article  CAS  Google Scholar  * Northrup, S. H. & Mccammon, J. A. Gated


reactions. _J. Am. Chem. Soc._ 106, 930–934 (1984). Article  CAS  Google Scholar  * Fersht, A. R. On the simulation of protein folding by short time scale molecular dynamics and distributed


computing. _Proc. Natl Acad. Sci. USA_ 99, 14122–14125 (2002). Article  CAS  PubMed  PubMed Central  Google Scholar  * Changeux, J. P. Responses of acetylcholinesterase from Torpedo


marmorata to salts and curarizing drugs. _Mol. Pharmacol._ 2, 369–392 (1966). CAS  PubMed  Google Scholar  * Taylor, P. & Lappi, S. Interaction of fluorescence probes with


acetylcholinesterase. The site and specificity of propidium binding. _Biochemistry_ 14, 1989–1997 (1975). Article  CAS  PubMed  Google Scholar  * Berman, H. A., Yguerabide, J. & Taylor,


P. Fluorescence energy transfer on acetylcholinesterase: spatial relationship between peripheral site and active center. _Biochemistry_ 19, 2226–2235 (1980). Article  CAS  PubMed  Google


Scholar  * Berman, H. A. & Leonard, K. Interaction of tetrahydroaminoacridine with acetylcholinesterase and butyrylcholinesterase. _Mol. Pharmacol._ 41, 412–418 (1992). CAS  PubMed 


Google Scholar  * Radic, Z., Reiner, E. & Taylor, P. Role of the peripheral anionic site on acetylcholinesterase: inhibition by substrates and coumarin derivatives. _Mol. Pharmacol._ 39,


98–104 (1991). CAS  PubMed  Google Scholar  * Pang, Y.-P., Quiram, P., Jelacic, T., Hong, F. & Brimijoin, S. Highly potent, selective, and low cost bis-tetrahydroaminacrine inhibitors


of acetylcholinesterase: Steps toward novel drugs for treating Alzheimer’s disease. _J. Biol. Chem._ 271, 23646–23649 (1996). Article  CAS  PubMed  Google Scholar  * Kousba, A. A., Sultatos,


L. G., Poet, T. S. & Timchalk, C. Comparison of chlorpyrifos-oxon and paraoxon acetylcholinesterase inhibition dynamics: potential role of a peripheral binding site. _Toxicol. Sci._ 80,


239–248 (2004). Article  CAS  PubMed  Google Scholar  * Auletta, J. T., Johnson, J. L. & Rosenberry, T. L. Molecular basis of inhibition of substrate hydrolysis by a ligand bound to the


peripheral site of acetylcholinesterase. _Chem. Biol. Interact._ 187, 135–141 (2010). Article  CAS  PubMed  PubMed Central  Google Scholar  * Ye, K. et al. Characterization of an


alternative BAK binding site for BH3 peptides. _Nat. Commun._ 11, 3301 (2020). Article  CAS  PubMed  PubMed Central  Google Scholar  * Pang, Y.-P. et al. Bak conformational changes induced


by ligand binding: Insight into BH3 domain binding and Bak homo-oligomerization. _Sci. Rep._ 2, 257 (2012). Article  PubMed  PubMed Central  Google Scholar  * Wang, Q. & Pang, Y.-P.


Preference of small molecules for local minimum conformations when binding to proteins. _PLoS ONE_ 2, e820 (2007). Article  PubMed  PubMed Central  Google Scholar  * Pang, Y.-P. FF12MC: a


revised AMBER forcefield and new protein simulation protocol. _Proteins_ 84, 1490–1516 (2016). Article  CAS  PubMed  PubMed Central  Google Scholar  * Berendsen, H. J. C., Postma, J. P. M.,


van Gunsteren, W. F., Di Nola, A. & Haak, J. R. Molecular dynamics with coupling to an external bath. _J. Chem. Phys._ 81, 3684–3690 (1984). Article  CAS  Google Scholar  * Darden, T.


A., York, D. M. & Pedersen, L. G. Particle mesh Ewald: An N log(N) method for Ewald sums in large systems. _J. Chem. Phys._ 98, 10089–10092 (1993). Article  CAS  Google Scholar  * Pang,


Y.-P. Low-mass molecular dynamics simulation for configurational sampling enhancement: more evidence and theoretical explanation. _Biochem. Biophys. Rep._ 4, 126–133 (2015). PubMed  PubMed


Central  Google Scholar  * Joung, I. S. & Cheatham, T. E. Determination of alkali and halide monovalent ion parameters for use in explicitly solvated biomolecular simulations. _J. Phys.


Chem. B_ 112, 9020–9041 (2008). Article  CAS  PubMed  PubMed Central  Google Scholar  * Pang, Y.-P. How fast fast-folding proteins fold in silico. _Biochem. Biophys. Res. Commun._ 492,


135–139 (2017). Article  CAS  PubMed  Google Scholar  * Feenstra, K. A., Hess, B. & Berendsen, H. J. C. Improving efficiency of large time-scale molecular dynamics simulations of


hydrogen-rich systems. _J. Comput. Chem._ 20, 786–798 (1999). Article  CAS  Google Scholar  * Therneau, T. M. & Grambsch, P. M. _Modeling Survival_ _Data: Extending the Cox Model_


(Springer-Verlag, 2000). * Contreras-Garcia, J. et al. NCIPLOT: a program for plotting noncovalent interaction regions. _J. Chem. Theory Comput._ 7, 625–632 (2011). Article  CAS  PubMed 


PubMed Central  Google Scholar  Download references ACKNOWLEDGEMENTS This work was supported by the US Army Research Office (W911NF-16-1-0264) and the Mayo Foundation for Medical Education


and Research. Responsibility for the information and views in this study lies entirely with the author. The author acknowledges the computing resources provided by the University of


Minnesota Supercomputing Institute and the Mayo Clinic high-performance computing facility at the University of Illinois Urbana-Champaign National Center for Supercomputing Applications.


AUTHOR INFORMATION AUTHORS AND AFFILIATIONS * Computer-Aided Molecular Design Laboratory, Mayo Clinic, Rochester, MN, USA Yuan-Ping Pang Authors * Yuan-Ping Pang View author publications You


can also search for this author inPubMed Google Scholar CONTRIBUTIONS Conceptualization, methodology, investigation, visualization, funding acquisition, project administration, supervision,


and writing: Y.P.P. CORRESPONDING AUTHOR Correspondence to Yuan-Ping Pang. ETHICS DECLARATIONS COMPETING INTERESTS The author declares no competing interests. PEER REVIEW PEER REVIEW


INFORMATION _Communications Chemistry_ thanks Zhiyong Zhang, Baoqing Zhou, and Ferran Feixas for their contribution to the peer review of this work. Peer reviewer reports are available.


ADDITIONAL INFORMATION PUBLISHER’S NOTE Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations. SUPPLEMENTARY INFORMATION PEER


REVIEW FILE DESCRIPTION OF ADDITIONAL SUPPLEMENTARY FILES SUPPLEMENTAL MATERIAL VIDEO S1 VIDEO S2 VIDEO S3 VIDEO S4 VIDEO S5 VIDEO S6 VIDEO S7 VIDEO S8 VIDEO S9 VIDEO S10 VIDEO S11 VIDEO S12


VIDEO S13 VIDEO S14 VIDEO S15 VIDEO S16 VIDEO S17 VIDEO S18 DATA S1 DATA S2 DATA S3 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/. Reprints and permissions ABOUT THIS ARTICLE CITE THIS ARTICLE Pang, YP. How neocarcerand Octacid4 self-assembles with guests into irreversible


noncovalent complexes and what accelerates the assembly. _Commun Chem_ 5, 9 (2022). https://doi.org/10.1038/s42004-022-00624-4 Download citation * Received: 20 July 2021 * Accepted: 30


December 2021 * Published: 20 January 2022 * DOI: https://doi.org/10.1038/s42004-022-00624-4 SHARE THIS ARTICLE Anyone you share the following link with will be able to read this content:


Get shareable link Sorry, a shareable link is not currently available for this article. Copy to clipboard Provided by the Springer Nature SharedIt content-sharing initiative


Trending News

Jalisco officials say this indigenous leader died in a car crash. His community isn't buying it

The Indigenous Wixárika community in northern Jalisco is raising doubts about the official explanation that a vehicle ac...

Microfluidic electronic paper | Nature Photonics

Pigmented inks in a microfluidic structure provide a new approach for fabricating bright, colourful electronic paper wit...

What we’ve learnt from building africa’s biggest genome library

The human genome was first sequenced in 2003 by multiple research centres across the world. The breakthrough was hailed ...

At an average of 1,200 passengers per day, maya train numbers far from eventual target

Will the Maya Train railroad eventually become a popular mode of travel between beach resorts, colonial cities and archa...

Over 500 chiapas residents flee cartel violence into guatemala

Hundreds of Mexican families fleeing cartel violence in Mexico’s southern state of Chiapas have sought refuge across the...

Latests News

How neocarcerand octacid4 self-assembles with guests into irreversible noncovalent complexes and what accelerates the assembly

ABSTRACT Cram’s supramolecular capsule Octacid4 can irreversibly and noncovalently self-assemble with small-molecule gue...

Response of pacific-sector antarctic ice shelves to the el niño/southern oscillation

ABSTRACT Satellite observations over the past two decades have revealed increasing loss of grounded ice in West Antarcti...

Risks and rewards in 2016: a year in review (opinion)

[embedded content] What’s life without a little risk? In order to achieve great things, we must continue to push ourselv...

How rangers can bridge 14-year gap and go where no scottish club has gone before

The hand of history will be on Rangers’ shoulders as they embark on their Europa League last 16 tie against Red Star Bel...

National apprenticeship week 2015: deputy prime minister's message

News story NATIONAL APPRENTICESHIP WEEK 2015: DEPUTY PRIME MINISTER'S MESSAGE Nick Clegg marks National Apprentices...

Top