Weak interaction-steered evolution of polyoxovanadate-based metal-organic polyhedra from transformation via interlock to packing

Yn-Hu Wng ,Ke-Wei Tong ,Cho-Qin Chen ,Jing Du ,Peng Yng,*

a College of Chemistry and Chemical Engineering,Advanced Catalytic Engineering Research Center of the Ministry of Education,Hunan University,Changsha 410082,China

b College of Chemistry and Materials Science,Testing and Analysis Center,Hebei Normal University,Shijiazhuang 050024,China

Keywords: Polyoxometalates Polyoxovanadates Metal-organic polyhedron Interlock Weak intermolecular forces

ABSTRACT A bottleneck in biomimetic synthesis consists in the full copy of,for example,the hierarchical structure of proteins directed by weak interactions.By contrast with covalent bonds bearing definite orientation and high stability,weak intermolecular forces within a continuous dynamic equilibrium can be hardly tamed for molecular design.In this endeavor,a ligand-dominated strategy that embodies tunable electrostatic repulsion and π···π stacking was first employed to shape polyoxovanadate-based metal-organic polyhedra (VMOPs).Structural evolution involving transformation,interlock,and discovery of an unprecedented prototype of the Star of David was hence achievable.Not only as a handy tool for the primary structural control over VMOPs,these weak forces allow for an advanced management on the spatial distribution of such manmade macromolecules as well as the associated physicochemical behaviors,representing an ideal model for simulating and interpreting the conformation-function relationship of proteins.

Weak intermolecular forces (e.g.,hydrogen bonds,π···πstacking,and van der Waals interactions) are vital to master the spatial structure of biomacromolecules as well as their associated functions [1].Typical examples include the complementary base pairing in DNA and theα-helix/β-sheet in proteins.With the upsurge of biomimetic synthesis,chemists are wild about a higher-order molecular designviaan ensemble of both covalent bonds and weak forces [2-4].It is in this context that the research of metalorganic frameworks (MOFs) and metal-organic polyhedra (MOPs)has achieved encouraging progress [5-7].Thanks to the guidance by reticular chemistry,their fundamental coordination scaffolds are now programmable,on the bench orinsilicocomputations,viaan elaborate combination of secondary building units and organic linkers of various topologies [7,8].Yet,predictions on the exact structural details are not always reliable as the existence of sophisticated intermolecular weak interactions.It has been established that,along with the ligand extension and resultant enlarged porosity,MOFs prefer to adopt framework interpenetration spontaneously (usually driven by hydrogen bonds andπ···πstacking),in order to get rid of the extra energy that causes metastability [9,10].In contrast to MOFs of infinite networks,the discrete MOPs apparently suffer a lot more from the weakly interacting constraints inside and outside their cage-like structures [11,12].Construction of aesthetically pleasing architectures,including distorted and interlocked MOPs are,of course,dominated by coordination interactions,yet the contribution of weak intermolecular forces cannot be ignored [9-12].As a consequence,it would be challenging and rewarding to understand and take advantage of such weak forces as an effective tool to steer the self-assembly of MOPs.

Amongst MOPs built upon polynuclear metallic vertices,polyoxovanadate-based metal-organic polyhedra (VMOPs),in particular,represent a subset with great developing prospects [13-16].Facilitated by an atom-by-atom control over both the addenda and heterogroup,the alkoxo-vanadates are capable of being customized into three-(e.g.,[VIV6O6(OCH3)9(SO4)(COO)3]2-,{V6S}) [17],four-(e.g.,[VIV4VVO9Cl(COO)4]2-,{V5Cl}) [18],and five-connected (e.g.,[WVIVIV5O11(SO4)(COO)5]3-,{V5W}) [19] nodes (Fig.S1 in Supporting information).Of these,in combination with a rich reservoir of edge and panel linkers,a large family of Platonic and Archimedean solids have been established with tetrahedral,octahedral,cubic,and spherical configurations [17-29].In this research frontier,Wang,Su,and Fangetal.have greatly advanced the heuristic findings of Hartl [13] and Zaworotkoetal.[14,30],especially in the structural engineering of modulesviacoordination interactions [17,20-29].From the research perspective of weak forces,on the other side,VMOPs could be considered as viable candidates on account of their tunable charge distribution,adaptive pore size,and distinct template effect.Indeed,the emergence of supramolecular recognition and isomerization in VMOPs delivers a positive signal to exploit the underlying weak interactions in such a system [26,31].Referenced by the seminal achievements in the related fields of polyoxometalates and molecular cages,weak intermolecular forces are nowadays valued for shaping high-nuclear clusters [32],as well as for directing the steric configuration of manmade macromolecules [33-44],which are of central importance in biological modeling synthesis.

Adoption of a ligand-dominated strategy that embodies tunable electrostatic repulsion andπ···πstacking,we attempted to tap the potential of such weak forces on the structural regulation of VMOPs.In our work,a series of dipolar ligands with varied lengths and angles were pre-synthesized to scaffold VMOPs towards a metastable state,which would set the stage for watching the possible weak interactions formed spontaneously.In addition to evaluating their type and strength,such forces were hired for a multi-scale structural design from basic geometryviabimolecular interlock to spatial packing,in a manner analogous to the hierarchical proteins evolved in nature.

By means of one-pot reactions under solvothermal conditions,a group of five linear dicarboxylic ligands that vary in length,namely H2AzDC (azobenzene-4,4′-dicarboxylic acid),H2EDC(4,4′-ethynylenedibenzoic acid),H2TPDC ([p-terphenyl]-4,4′′-dicarboxylic acid),H2QPDC ([1,1′:4′,1′′:4′′,1′′′-quaterphenyl]-4,4′′′-dicarboxylic acid),H2PEBDC {(4,4′-(1,4-phenylenebis(ethyne-2,1-diyl))dibenzoic acid)},were selected to assemble with thein-situformed {V6S} for building VMOPs of the supposed tetrahedral archetype (see Supporting information).Beyond our expectation,the tetramethylammonium (TMA) salts of three anticipated “tetrahedra” (VMOTs) and two unprecedented “capsules” (VMOCs) were obtained respectively:

(TMA)8{[V6O6(OCH3)9(SO4)]4(AzDC)6}·16CH3OH·15DMF

(TMA-VMOT-1);

(TMA)8{[V6O6(OCH3)9(SO4)]4(EDC)6}·xCH3OH·yDMF

(TMA-VMOT-2) [31];

(TMA)8{[V6O6(OCH3)9(SO4)]4(TPDC)6}·17CH3OH·11DMF

(TMA-VMOT-3);

(TMA)4{[V6O6(OCH3)9(SO4)]2(QPDC)3}·3CH3OH·3DMF

(TMA-VMOC-4);

(TMA)4{[V6O6(OCH3)9(SO4)]2(PEBDC)3}·13CH3OH·13DMF

(TMA-VMOC-5).

In general,the applied ligand-elongation strategy could largely vary the pore size but not the topology of backbones,as previously indicated in the study of IRMOFs [45].As expected,two regular tetrahedra ofVMOT-1andVMOT-2were received using AzDC and EDC as respective linkers.In a typical tetrahedral geometry,four {V6S} motifs serving as three-connected vertices are bridged by six organic ligands as edges,of which the length fall into the range of 13.0-13.6 ˚A (Figs.1a and b).If displacing the linker with a prolonged TPDC (15.7 ˚A),a concave tetrahedron ofVMOT-3was obtained (Fig.1c).It has been noted that the six edges bend inwards the body center,resulting in curves with an average radian of 0.14π(Fig.S2a in Supporting information).Further elongation of the ligands to QPDC (20.0 ˚A) and PEBDC (20.7 ˚A) led to a drastic conversion from the classical tetrahedron to an unprecedented capsule-style cage.In the structure ofVMOC-4andVMOC-5,two {V6S} units occupy the two poles and are vertically bridged by three ligands,with an idealized point group symmetry ofD3h(Figs.1d and e).On the contrary to the bending mode inVMOT-3,the edges herein bulge outwards in the middle with an average radian of 0.24π(VMOC-4) and 0.18π(VMOC-5),respectively(Figs.S2b and c in Supporting information).It follows from the above that a rare structural evolution of VMOPs,from tetrahedronviaa concave intermediate to capsule,was achieved as a function of ligand length (Fig.1f).To explain this from the standpoint of weak interactions,a systematic analysis on their electrostatic potential (ESP) diagrams was carried out [46].

Fig.1.Transformation steered by ligand extension of varied electrostatic potential.(a-e) Combined polyhedral/stick representation of VMOT-1～3,VMOC-4,and VMOC-5.Color code: {VO6},turquoise octahedra;S,yellow;O,red;N,blue;C,gray.Hydrogen atoms are omitted for clarity.(f) Structural evolution from tetrahedron (VMOT-1 and VMOT-2) via a concave intermediate (VMOT-3) to capsule (VMOC-4 and VMOC-5).(g) Electrostatic potential map (isovalue 0.01) of the representative ligands,and (h) of the corresponding VMOT-BDC,VMOT-1,VMOT-3,as well as (i) VMOC-5.Blue: electron accumulation;red: electron depletion.

As depicted in Fig.1g,the ESP of several representative dicarboxylic acid-based ligands were calculated.From the perspective of electric potential distribution,the more negative ESP values are predominantly located on the carboxylic functions,while the ESP values in the middle of the ligands are relatively positive.As the length of ligands being extended from BDC (terephthalate) to PEBDC,the nucleophilicity of the central aromatic parts decreases correspondingly.So,a progressive separation of ESP was clearly evidenced in the linear ligands being used.As shown in Fig.1h,the ESP surfaces of VMOT regions are primarily concentrated on the semi-exposed sulfate heterogroups.Therefore,it is speculated that the electrostatic repulsion between the sulfate functions and peripheral aromatic domains is essential for supporting the tetrahedral configuration.From the referencedVMOT-BDC({[V6O6(OCH3)9(SO4)]4(BDC)6},Fig.S3a in Supporting information)[17] toVMOT-3,such repulsive forces are gradually weakened till the formation of a pressure-bearing concave state with inadequate inner support.Upon further stretch of the edges,the assumed tetrahedra might become unstable and be forced to shrink into a compact capsule-like archetype (Fig.1i).Accordingly,the opening angles in the three-connected {V6S} are shifted from a diverging pattern (e.g.,θ=29.5° inVMOT-1,Table S1 in Supporting information) to a convergent fashion (e.g.,θ=12.0° inVMOC-5),in an effort to gather the ligands and restore the internal electrostatic forces.As demonstrated by the thermogravimetric analysis,the decomposition ofTMA-VMOT-1andTMA-VMOC-4was ended at ≥500°C (Figs.S4a and b in Supporting information),whileTMA-VMOT-3was completely destructed at a lower temperature (468°C,Fig.S4c in Supporting information),testifying its metastability from another view.

Inspired by these findings,it interested us more to verify whether such electrostatic forces are able to steer the structural transformation on purpose.Previous research has experienced the deformation of VMOPsviathe substitution of sulfate ({V6S}) by a sterically demanding phenylphosphate ({V6P}) within a confined space [27-29].Different from the above,the void inVMOT-1～3are large enough to ignore the steric hindrance of heterogroups,while the distribution of ESP on {V6P} should be varied.In this way,the replacement of {V6S}→{V6P} was proceeded in the scenario of as-made compounds.Of prime interest,the original tetrahedra ofVMOT-1～3were converted into their respective capsulestyle counterpartsVMOC-P-1～3(Figs.2a-c).Besides,the backbone of the derivedVMOC-P-4andVMOC-P-5remain nearly unchanged(Figs.2d and e).The formula of their TMA salts are listed below:

(TMA)4{[V6O6(OCH3)9(PhPO3)]2(AzDC)3}·6CH3OH·3DMF

(TMA-VMOC-P-1);

(TMA)4{[V6O6(OCH3)9(PhPO3)]2(EDC)3}·7CH3OH·5DMF

(TMA-VMOC-P-2);

(TMA)4{[V6O6(OCH3)9(PhPO3)]2(TPDC)3}·11CH3OH·7DMF

(TMA-VMOC-P-3);

(TMA)4{[V6O6(OCH3)9(PhPO3)]2(QPDC)3}·3CH3OH·4DMF

(TMA-VMOC-P-4);

(TMA)4{[V6O6(OCH3)9(PhPO3)]2(PEBDC)3}·3CH3OH·5DMF

(TMA-VMOC-P-5).

Apparently,the ESP analysis showcased a reduced nucleophilicity on the heterogroup after the exchange of {V6S}→{V6P}(Fig.2f).As a result,the electrostatic repulsion in the cage could be further turned down,thereby aiding the conversion from tetrahedron (VMOT-1～3) to capsule (VMOC-P-1～3),which is in full agreement with the experimental data.It should also be stressed that,in the previous study on vertice replacement,the tetrahedral backbone with relatively short edges are maintained inVMOT-P-BDC({[V6O6(OCH3)9(PhPO3)]4(BDC)6},Figs.S3b and c in Supporting information) andVMOT-PNDC({[V6O6(OCH3)9(PhPO3)]4(NDC)6},NDC=naphthalene-2,6-dicarboxylate,Fig.S3d in Supporting information) [27,47].Thus it can be seen that,a manual switch of configuration is accessible in VMOPs on condition of prolonged linkers being used.By virtue of the regulation on the difference of ESP either on the edge or vertice,the processable weak electrostatic repulsion could assume more responsibility for the structural engineering of such macromolecules.To visualize the region of non-covalent interactions withinVMOC-P-3,we applied the independent gradient model based on the Hirshfeld partition (IGMH) method [48] to process the calculated electron densities in the Multiwfn code [49].The van der Waals interactions were colored by green and imaged with a Visual Molecular Dynamics (VMD) 1.9.3 software [50].From the IGMH map (Fig.2g),it can be seen that a thin and broad isosurface appears between the phenyl group in {V6P} and nearby aromatic linkers,which ideally exhibitsπ···πand C-H···πinteractions.Through the substitution of {V6S}→{V6P},additional van der Waals interactions could be generated between the vertices and linkers,which would further weaken the repulsion between them and hence in favor of a capsule-like architecture.A lower and more-stable energy state after the structural transformation was again proved by thermogravimetric analysis.As exemplified byTMA-VMOT-3→TMA-VMOC-P-3,the final decomposition temperature was markedly increased from 468°C to 501°C (Fig.S4d in Supporting information).Meanwhile,the emergence ofπ···πstacking was also defined in the crystal structure of all {V6P}-based capsules (e.g.,centroid-to-centroid distance of 3.92 ˚A inVMOC-P-3,Fig.2h),reminding us to exam theπ···πinteractions appeared in the whole system.

Fig.2.Transformation steered by vertice replacement of tunable electrostatic potential and IGMH analyses.(a-c) Structural transformation from the tetrahedral (VMOT-1～3)to capsule-like archetype (VMOC-P-1～3) via vertice replacement.(d,e) Vertice exchange from VMOC-4 and VMOC-5 to VMOC-P-4 and VMOC-P-5.(f) Electrostatic potential map (isovalue 0.01) of {V6S} and {V6P}.(g) The sign(λ2)ρ colored isosurfaces of δginter=0.001 a.u.of VMOC-P-3 corresponding to IGMH analyses (The color bar represents the common interpretation of coloring method of mapped function sign(λ2)ρ in IGMH maps).(h) The π···π interactions (dotted green line) in VMOC-P-3.

A detailed review on the structure ofVMOT-3disclosed a facedirected interlock of two concave tetrahedra,resulting in a shuttleshaped supramolecular assembly (Figs.3a and b).Of the two pieces being interpenetrated,one is flipped over 180° and then hooked with the other on the central belt.The view from the top pictures a hexagram made up of two stagger basal planes,which are glued by six groups ofπ···πstacking (Fig.3c).Specifically,the phenyl rings in the middle of the three linkers (purple) adopt faceto-face interactions with the corresponding aromatic groups next to the vertices (brown),with the centroid-to-centroid distances ranging from 3.76 ˚A to 3.91 ˚A.Based on the semi-stable tetrahedron ofVMOT-3,the presence of bimolecular interpenetration is understandable as to spontaneously fill the large pores and improve their mutual stability.Being calculated by the Multiwfn code [48] using the Bondi van der Waals radius [51],the pore volume inVMOT-3is compressed from 1337.98 ˚A3(individual) to 635.63 ˚A3(interlocked),almost halving the central space (Fig.S5 in Supporting information).Fig.3d presents thesign(λ2)ρcolored IGMHδginterisosurfaces of 0.001 a.u.,the two individual tetrahedra are respectively defined as separate fragments.Notably,evidentπ···πinteractions are positioned in the middle aromatic area of ligands,as there exists a flat and broadδginterisosurface.It is apparent that,with the increasing number of aromatic groups,the ESP distribution in the middle of the linkers would be weakened,whereas the intermolecularπ···πinteractions get strengthened instead,which jointly power the interpenetration of such tetrahedra.As another kind of weak forces found in the VMOP system,the role ofπ···πstacking in the management of bimolecular aggregation states deserves further investigation.

Fig.3.Structural and IGMH analyses of the face-and edge-interlock assembly.Structural representation of (a) the individual VMOT-3 and (b) its bimolecular face-interlock assembly.(c) The π···π interactions in the face-interlock assembly of VMOT-3.(d) The sign(λ2)ρ colored isosurfaces of δginter=0.001 a.u.corresponding to IGMH analyses for the face-interlock assembly of VMOT-3.Structural representation of e) the individual VMOC-6 and (f) its bimolecular edge-interlock assembly.(g) The π···π interactions in the edge-interlock assembly of VMOC-6.(h) The sign(λ2)ρ colored isosurfaces of δginter=0.0005 a.u.corresponding to IGMH analyses for the edge-interlock assembly of VMOC-6.

In an effort to introduceπ···πstacking into the family of VMOCs,H2SPDC (4,4′-sulfonyldibenzoic acid) was employed as linker with a flexible dihedral angle.Considering the electronwithdrawing sulfonyl group in the middle of the edge (i.e.,reduced nucleophilicity on the aromatic area,Fig.S6 in Supporting information),a capsule-like cage was envisioned initially because of the relatively weak electrostatic repulsion within it.Moreover,the central void in the prospective capsule might be enlarged due to the folding nature of edges,thereby creating favorable conditions forπ···πinteractions.As expected,a capsule with the formula(TMA)4{[V6O6(OMe)9(SO4)]2(SPDC)3}·6CH3OH·5DMF (TMA-VMOC-6) was harvested successfully.In the structure ofVMOC-6,the three edges display a bulgy configuration with the dihedral angles fall in the range of 106.00° to 113.23°,making room for a larger interior space (Fig.3e).In order to liberate the excess energy,two of such capsules are cross-linked perpendicularly to afford an edgeinterlock mode,which is stabilized by four sets ofπ···πinteractions (Fig.3f).For eachVMOC-6,the edge being entangled penetrates into the domain that is enclosed by two non-locked linkers of the other cage.Every aromatic ring of the former interacts with the two phenyls of the latter,with an average centroid-tocentroid distance of 4.20 ˚A (Fig.3g).As illustrated in Fig.3h,the IGMH analyses showcased the interactions ofπ···πstacking are mainly concentrated in the center of the dimer,and are much enhanced compared to that in the face-interlockVMOT-3.This tells us that,by adjusting the folding angle of edges,additional sites forπ···πinteractions inside the capsule would be provided,giving rise to a tighter interlocked pattern of narrow space (65.66 ˚A3,Fig.S7 in Supporting information).From the above,it is confirmed thatπ···πstacking could be engaged to customize the mechanically interlocked VMOPs of varying archetypes,thereby broadening the scope of design from individual molecule to supramolecular aggregates.

What if the two weak forces are simultaneously amplified in one system? In other words,a synergistic strategy of ligand stretch and vertice replacement was taken to severely weaken the electrostatic repulsion.At the same time,expanding further the central cavity to maximize the satisfaction ofπ···πstacking and framework interpenetration.To meet these requirements,a customized ligand of H2DBIDC (1,4-bis(5-carboxy-1Hbenzimidazole-2-yl) benzene) was involved in the synthesis,and a“hexagram”-like prototype (Fig.4a) was obtained for the first time:(TMA)12{[V6O6(OCH3)9(PhPO3)]6(DBIDC)9}·40CH3OH·15DMF (TMAVMOH-P-1).

Fig.4.Structural and IGMH analyses of the Star of David.(a) The assembly process of VMOH-P-1.(b) The face-edge-interlock mode in VMOH-P-1.(c) The π···π interactions in VMOH-P-1.(d) The sign(λ2)ρ colored isosurfaces of δginter=0.001 a.u.corresponding to IGMH analyses for VMOH-P-1.

In view of the prolonged linker (19.0 ˚A) and {V6P} vertice of reduced nucleophilicity,a capsule-like cage could be foreseen (Fig.S8 in Supporting information).However,the V-shaped coordination of DBIDC significantly distorts the pristine capsule,and finally breaks it from one side.Then the two arms being thrown out are united with another {V6P} to re-form a basket-like assembly (Fig.4b).The remained connection point on the handle of each basket are occupied by a shared Z-shaped DBIDC,which pierces through the subface of both and bounds them together.It is obvious that the open basket encompasses a larger nanoscale central void,wherein the distance from the bottom center to apex (the outmost carbon atom in {V6P}) is 15.58 ˚A (Fig.S9 in Supporting information).In this context,a mixed face-edge-interlock pattern directed byπ···πstacking was generated.In the center ofVMOH-P-1,two sets of face-to-face aromatic interactions are located,with the centroid-tocentroid distance of 3.80 ˚A (Fig.4c).According to the IGMH analyses shown in Fig.4d,installation of aromatic groups on both the vertices and edges would make for more intense intramolecular weak interactions presented inVMOH-P-1.From theδginterisosurface,in addition to the identifiedπ···πstacking amongst the linkers in the center,the C-H···πinteractions appear as well between the ligands and heterogroups around vertices.It is certain that a modification on both the ligand and vertice would promote the communication of weak forces withinVMOH-P-1,and further orchestrate the face-and edge-interlock modes in one system.As all said above,the integration of both electrostatic repulsion andπ···πstacking is full of unlimited potential,especially in the development of spectacular architectures beyond the classical polyhedra of high symmetries.

To deepen our understanding in the weak interaction-steered VMOPs,the research focus was moved from the primary structures towards their distribution in the three-dimensional (3D) space.On top of this is disclosure of the relationship between the molecular spatial arrangement and associated physicochemical properties.Directed by electrostatic repulsion andπ···πstacking,the evolution of VMOPs is reflected not only in their structural types but also the adjustable electronegativity.From capsuleviatetrahedron to hexagram,the stoichiometric ratio of verticesVSlinkers are 2:3,4:6,and 6:9,respectively,corresponding to a gradually enhanced negative charge of -4,-8,and -12.Herein,the {V6P}-based capsule (TMA-VMOC-P-2),tetrahedron (TMA-VMOT-P-NDC),and hexagram (TMA-VMOH-P-1) were selected for dye adsorption.Similar to the behavior of reported VMOPs [22,26,28,47],TMA-VMOCP-2(Fig.5a),TMA-VMOT-P-NDC(Fig.5b),andTMA-VMOH-P-1(Fig.5c) exclusively adsorbed the cationic MB+(methylene blue)rather than the neutral SD0(Sudan I,Fig.S10 in Supporting information) or anionic AO-(acid orange,Fig.S11 in Supporting information).It was noted that VMOPs bearing higher negative charges(VMOH-P-1＞VMOT-P-NDC＞VMOC-P-2) behaved a faster adsorption rate and larger loading capacity (Fig.5d),which could be attributed to the cation-exchange mechanism between MB+and TMA+molecules.

Fig.5.Relationship of spatial packing and dye adsorption performance.(a-c) The UV-vis spectra of MB+ in ethanol in the presence of TMA-VMOC-P-2,TMA-VMOT-P-NDC,and TMA-VMOH-P-1 at various time intervals (inset,photographs of the color of dye solutions before and after adsorption).(d) The comparison of adsorption performance of selected VMOPs.(e) The crystal packing diagram of TMA-VMOC-P-2 ordered by π···π stacking.(f) The crystal packing diagram of TMA-VMOC-P-5 without π···π interactions.(g) The UV-vis spectra of MB+ in ethanol in the presence of TMA-VMOC-P-5 at various time intervals,as well as the adsorption performance compared to that of TMAVMOC-P-2.

It is noteworthy that the equally charged capsules inTMAVMOC-P-1～5,on the other hand,adopt two types of packing modes.Along with the ligand extension from AzDC to PEBDC the corresponding radian of edges gradually increases from 0.12πto 0.20π(Figs.S12-S14 in Supporting information).ForTMA-VMOCP-1～3,the neighboring capsules are ordered by two or three sets ofπ···πstacking (e.g.,centroid-to-centroid distance of 4.15 ˚A inTMA-VMOC-P-2),and are further cascaded into a wave-like 1D chain.Such chains are packed in parallel with each other,resulting in a supramolecular network with zigzag channels (Fig.5e).By sharp contrast,in the spatial arrangement ofTMA-VMOC-P-4andTMA-VMOC-P-5,their significantly bulgy edges might exert additional steric constraints to diminish the probability ofπ···πinteractions amongst themselves.As a result,the adjacent capsules are arranged above and below the same layer alternatively,and a honeycomb-like framework with hexagonal channels was achieved instead (Fig.5f).Interestingly,a better adsorption performance ofTMA-VMOC-P-5was reached as compared to that ofTMA-VMOCP-2(Fig.5g),albeit the channels in both are large enough to store MB+molecules.One possible explanation for this may be the distinct microenvironment forπ···πinteractions between the capsule and MB+molecules.Since most of the aromatic groups are taken up for theπ···πstacking of capsules in between,MB+are unlikely to enter the channels inTMA-VMOC-P-2without appropriate binding sites (Fig.S15 in Supporting information).InTMA-VMOC-P-5,however,the three aromatic linkers are all free and face to the center of the hexagons,which is propitious to accommodate MB+viaπ···πinteractions (Fig.S16 in Supporting information).From the above,weak intermolecular forces can be taken as chemist’s right hand to challenge the spatial control over macromolecules of higher difficulty.Therefore,the weak force-directed VMOPs are on the right track for a complete simulation of protein’s multi-level structures,as well as the associated conformation-function relationship.

In response to ligand design,the spontaneous electrostatic repulsion andπ···πstacking were magnified and captured together with the structural evolution of VMOPs.These weak interactions were,in turn,hired to direct the structural transformation,interlock,and even to develop a Star-of-David prototype of VMOPs,either alone or in collaboration.But more than that was the successful customization on the packing mode has been proved to be competent for a protein-mimicking,hierarchical management on the self-assembly of VMOPs,from coordination framework geometriesviabimolecular aggregation states to spatial distribution characteristics.Our work promises a feasible route for the design of VMOPs by electrostatic repulsion andπ···πstacking,in addition to the conventional modular matching.Also,being dramatically responsive to the weak intermolecular forces,such a system represents an ideal model for simulating and elucidating the conformation-function mechanism in proteins at a molecular scale.

Declaration of competing interest

The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.

Acknowledgments

This research was financially supported by the National Natural Science Foundation of China (No.22001066),the Natural Science Foundation of Hunan Province (Nos.2021JJ40049 and 2022JJ20007),and the Science and Technology Innovation Program of Hunan Province (No.2022RC1115).J.Du acknowledges the Science and Technology Project of Hebei Education Department (No.QN2023049),Science Foundation of Hebei Normal University (No.L2023B51),and TianHe-2 (LvLiang,China) Cloud Computing Center for support.

Supplementary materials

Supplementary material associated with this article can be found,in the online version,at doi:10.1016/j.cclet.2023.109066.