Literature DB >> 24641469

Protein-triggered supramolecular disassembly: insights based on variations in ligand location in amphiphilic dendrons.

Diego Amado Torres1, Matteo Garzoni, Ayyagari V Subrahmanyam, Giovanni M Pavan, S Thayumanavan.   

Abstract

We use monodisperse dendrons that allow control over functional group presentation to investigate the influence of the location of a ligand on protein-induced disassembly and release of encapsulated small molecules. Based on both experiments and molecular dynamics simulations, we demonstrate that ligand location greatly influences release of guest molecules from the dendron-based supramolecular assembly. We show that a ligand moiety grafted to the dendron periphery is more accessible for the target protein in aqueous solution. On the other hand, the ligand moiety placed at the focal point or at the intermediate layer within the dendritic scaffold is less accessible, since it is surrounded by an environment rich in PEG chains, which hinders binding and even influences nonspecific interactions. We also demonstrate that the specific binding between one ligand and the target protein can destabilize the dendritic assembly. Furthermore, if more ligands are available, multivalent interactions are also possible with extravidin, which speed up disassembly and trigger the release of hydrophobic guests.

Entities:  

Mesh:

Substances:

Year:  2014        PMID: 24641469      PMCID: PMC4004214          DOI: 10.1021/ja500634u

Source DB:  PubMed          Journal:  J Am Chem Soc        ISSN: 0002-7863            Impact factor:   15.419


Introduction

Host systems that can spontaneously self-assemble and stably encapsulate guest molecules under one set of conditions, but disassemble and release the guest molecules when external conditions change, have been of great interest in supramolecular chemistry due to implications in a variety of biomedical applications.[1] While earlier studies have primarily focused on the former (i.e., self-assembly and binding), there has been a recent surge in interest in the latter features (disassembly and release in response to an environmental change). Also, the responsiveness of the assemblies has primarily focused on physical or chemical changes. For example, there has been a significant interest in systems that disassemble in response to chemical stimuli such as pH[2] or redox[3] variation as well as physical stimuli such as light,[4] temperature,[5] or a magnetic field.[6] While systems incorporating features sensitive to light and magnetic field have been designed to respond to external triggers for biological applications, pH- and redox-responsive systems have been designed to respond to the inherent imbalances observed in certain disease tissues. The anomalous pH or redox conditions in disease locations can be considered to be secondary imbalances in biology, as the primary imbalances are often the result of aberrant protein concentrations or enzymatic activity.[7] Therefore, there is a growing interest in developing supramolecular systems that respond to these primary factors in biology.[8] Protein-responsive systems can be broadly classified into two categories, viz. covalent modification of the hosts to disable their capacity to hold the guest molecules and noncovalent binding to the host assemblies to produce the same effect.[9,10] The former is often achieved by an enzyme-driven chemical reaction that modulates the host characteristics of the molecule. In the latter case, this is achieved due to a noncovalent binding interaction. While there have been several systems designed to be degraded or covalently modified by enzymes,[10] supramolecular assemblies that lose their host capacities due to noncovalent binding with proteins have been limited. Among the systems that undergo binding-induced disassembly, polyelectrolyte assemblies that can noncovalently bind a complementarily charged surface to cause a disassembly have attracted some attention.[11] Although nonspecific in its interaction, the simplicity of these systems has proved useful in applications such as separations.[12] It has also been clear that for a binding-induced approach to be useful in applications such as delivery and sensing, where specificity is critical, strategies that use specific ligand–protein interactions are needed.[13] Since dendritic macromolecules can be produced in high molecular weights, but with a great degree of control, these scaffolds have certain unique advantages for this strategy.[14] For example, the critical aggregate concentrations (CAC) of the dendrimer-based amphiphilic assemblies are low, an advantage that polymeric systems have over small molecule-based amphiphilic assemblies. Similarly, the control over functional group placements in dendritic architectures captures the advantage that small molecules have over their polymeric counterparts.[15] These features allow for structure–property relationship studies that unravel the factors that underlie the binding-induced supramolecular disassembly process. In this manuscript, we take advantage of this unique feature by incorporating ligands at specific locations within facially amphiphilic dendrons and interrogating the efficiency of the supramolecular disassembly and molecular release in response to a specific protein binding event. We use both experiments and molecular dynamics (MD) simulations synergistically to gain insights into the ligand–protein recognition-based disassembly event.

Results and Discussion

Molecular Design and Synthesis

Facially amphiphilic dendrons, containing a biaryl-based internal repeat unit and an aryl peripheral unit, bear both hydrophilic and hydrophobic functional groups as side chains in each of these repeat units.[16] The amphiphilic functional groups are placed at opposite faces of the biaryl backbone of the dendron such that these molecules are endowed with the capability to form micelle-like assemblies in aqueous solution and form inverse micelle-like assemblies in apolar solvents.[17] In our preliminary findings, we have shown that placing a ligand moiety at the hydrophilic face of the dendron can provide binding-induced disassembly of the micelle-like structure in the aqueous solvent.[13a] In that work, the ligand moiety was placed at the focal point of the dendron, as this structure is synthetically most easily accessible. However, the release efficiency of guests, following protein binding, was relatively moderate. In order to fully realize the potential of this process, it is critical that we understand the effect of the functional group placement, thus taking advantage of one of the most prominent features of dendritic architectures in this supramolecular process. Therefore, in this manuscript, we investigate the effect of the placement of a protein-specific ligand moiety at specific positions of a dendron upon the accessibility of the complementary protein to the supramolecular assembly and hence the effects on the concomitant guest release response (Figure 1). The relative encapsulation of functional groups, when placed in the periphery vs the focal point of the dendron is well-known.[18] Note however that the facially amphiphilic dendrons, utilized here, have two distinct characteristics in this context: (i) the ligand moieties placed in the hydrophilic face of the dendron is attached through a rather long ethyleneglycol linker. Therefore, the classical encapsulation by the dendritic backbone should not apply here; (ii) it has been previously shown using carboxylic acid based facially amphiphilic dendrimers that these functional groups are available for nonspecific binding to the positively charged surface of chymotrypsin.[12f] For these reasons, it is important that we investigate the effect of ligand placement on the binding-induced release.
Figure 1

Structure of the G1 and G2 containing ligands at specific locations. The CAC of each of these dendrons are listed in the table. The scheme is an illustration of the self-assembling dendrons with sequestered hydrophobic guest molecules, which are released in response to protein binding. The equilibrium illustrated in the illustration is to indicate the possible unimer-aggregate equilibrium.

Structure of the G1 and G2 containing ligands at specific locations. The CAC of each of these dendrons are listed in the table. The scheme is an illustration of the self-assembling dendrons with sequestered hydrophobic guest molecules, which are released in response to protein binding. The equilibrium illustrated in the illustration is to indicate the possible unimer-aggregate equilibrium. The CAC of even the first generation of these dendrons is substantially lower than the corresponding amphiphilic small molecule (μM compared to mM). Note that this 3 orders of magnitude difference in CAC is not accounted for by the simple difference in molecular weight of the amphiphilic molecules. Within each generation of dendrons however, the CAC gain has been relatively small, if any, especially when one accounts for the difference in molecular weights. Therefore, we focused on the G1 and G2 dendrons, which provide sufficient variations in the functional group placement. We targeted five different dendrons within these two generations, where the ligand placement is the key difference (see Figure 1 for the structures). Pentaethyleneglycol monomethyl ether (PEG) was chosen as the hydrophilic moiety not only for its ability to provide the requisite solubility in the aqueous phase but also for its propensity to present a charge neutral surface on the assembly that exhibits minimal nonspecific interactions.[19] Decyl moieties were incorporated as the hydrophobic components of the amphiphile. Biotin was chosen as the ligand, because of its well-established high affinity toward avidin.[20] A ligand can be incorporated on a G1 dendron either at the focal point or at the periphery. On a G2 dendron, a ligand can be incorporated on three different layers: (i) the focal point; (ii) the middle layer; and (iii) the periphery. To install a single ligand at a specific place, one of the PEG units was replaced by the ligand, attaching it to the dendron backbone via a hydrophilic linker. This placement allows the ligand to be exposed to the aqueous solution on the hydrophilic face of the assembly. Note that the focal point is synthetically the easiest place to attach a single functional group, since the focal point is a single and unique position in a dendron. In contrast, to install a single functional group at the periphery, one has to distinguish one of the peripheral monomers from otherwise identical ones within the dendron. The G1 dendron has two such identical positions in the periphery, while the G2 dendron has four such locations. Similarly, the middle layer of the G2 dendron contains two similar positions that need to be distinguished for specifically placing a ligand moiety. Accordingly, the degree of difficulty in placing a unique functional group in the dendron increases as we move from the focal point to the periphery and as we increase the generation. To achieve these selective functionalization, we and others have developed a variety of synthetic strategies that afford multifunctionalized dendrons and dendrimers.[15,21] In all our syntheses, we initially targeted dendrons containing an acetylene moiety, which will be used as the handle to ‘click’ the biotin moiety in the last steps of the syntheses. Overall, we made the syntheses of the dendrons modular in order to assemble the targeted dendrons in small number of steps. Accordingly, we first synthesized the biaryl, amphiphilic AB2 monomer 1, the propargyl-functionalized periphery unit 2, and the nonfunctionalized amphiphilic periphery unit 4 following the previously reported procedures.[13a] To achieve the targeted dendron 5, the biaryl monomer 1 was first monoalkylated with a periphery unit 2 to get the propargyl-functionalized scaffold 3 (Scheme 1) in 25% yield. The low isolated yield was because of the statistical distribution of mono- and disubstituted products obtained in this reaction. Molecule 3 was then treated with the amphiphilic peripheral monomer 4 under similar alkylation conditions to achieve the G1 dendron 5 containing the reactive propargyl moiety at its periphery. The G1 dendron 5 was then treated with the azide-modified biotin 6 under alkyne–azide click chemistry conditions[22] in the presence of cupric sulfate and sodium ascorbate to obtain G1-P in 70% yield. Similarly, the molecule G1-F was prepared by clicking the biotin azide to the propargyl moiety present in the biaryl repeat unit (see Supporting Information for synthetic details).
Scheme 1

Synthesis of G1-P, Installment of Ligand at the Periphery

The syntheses of the targeted G2 dendrons were achieved using a similar set of synthetic strategies, as shown in Scheme 2. To achieve the syntheses of these dendrons, the previously reported[13a] amphiphilic G1 dendron 7 was treated with the biaryl monomer 1, where the monoalkylated product 8 was separated from a statistical mixture in 24% yield. This molecule was then treated with the bromomethyl dendron 9 (obtained from the precursor to G1-F in one step) or 10 (obtained from 5 in one step) to obtain the G2 dendrons 11 or 12, containing the propargyl moiety at the middle layer or the periphery of the dendron, respectively. Copper-catalyzed Huisgen reaction of 11 and 12 with the azide-functionalized ligand 6 afforded the G2-M and G2-P dendrons, respectively. As with G1-F, the G2-F dendron was obtained by clicking the biotin azide to the propargyl moiety present in the biaryl repeat unit (see Supporting Information for synthetic details).
Scheme 2

Synthesis of G2-M and G2-P

Assembly and Disassembly

First, the self-assembly properties of these dendrons were studied by measuring their CACs using Nile red as the hydrophobic, spectroscopic probe (Figure S1). As anticipated, the CACs of final G1 and G2 dendrons were determined to be in the low micromolar range with values of 11.65 and 8.02 μM for G1-F and G1-P, respectively; 7.25, 7.08, and 7.01 μM for G2-F, G2-M, and G2-P, respectively. With the installment of a single ligand at different layers of a dendron and after knowing the lowest concentration at which they assemble into nanostructures, we investigated the response of the different dendritic aggregates in presence of the complementary protein extravidin and in presence of noncomplementary proteins with diverse pI values and molecular weights, viz. α-chymotrypsin (Chy, pI = 8.1–8.6), pepsin (Pep, pI = 2.9), and myoglobin (Myo, pI = 7.2), as shown in Figure 2.
Figure 2

DLS-based measure of size change of dendritic assemblies in aqueous phase upon interaction with different proteins. (a) G1-F, (b) G1-P, (c) G2-F, (d) G2-M, (e) G2-P, and (f) G1 and (g) G2 control dendrons.

Prior to analyzing the interaction between the dendritic assemblies and the proteins, we analyzed the size of the assemblies in aqueous phase using dynamic light scattering (DLS), with G1 dendron concentrations of 12.5 μM and G2 dendron concentrations of 10 μM (both above their respective CACs). The sizes of the assemblies were found to be in a few tens of nanometers, ranging from ∼30 to ∼200 nm, 44 nm for G1-F, 220 nm for G1-P, 29 nm for G2-F, 92 nm for G2-M, and 51 nm for G2-P. The reason for the variations in size with the subtle change in the position of the ligand is not clear. However, note that the replacement of a hydrophilic PEG chain in a G1 or G2 dendron by a less hydrophilic pendant biotin decreases the PEG density on the assembly surface, decreasing hydrophilicity and increasing the chance of nonspecific interactions with noncomplementary proteins. Such decrease in PEG density could be more relevant in G1, where substituting one of the three PEG chains in a dendron molecule could mean a reduction of up to 33%, while such a reduction in G2 is about 14%. It is possible that these differences confer changes in the way the assemblies pack together in solution, which results in size variations among the biotin-functionalized dendrimeric assemblies. It is interesting that the size of all these assemblies reduced to about ∼13 nm in presence of the protein, extravidin (2 μM) (Figure 2). A particular difference in disassembly among the biotin-functionalized dendrons was observed in the case of the G2-P that disassembled into smaller aggregates (∼7–8 nm) in presence of extravidin (Figure 2e). G1-P also presented some deviations, showing larger aggregates (∼28 nm) that could be formed by further aggregation of smaller protein–dendron complexes. DLS-based measure of size change of dendritic assemblies in aqueous phase upon interaction with different proteins. (a) G1-F, (b) G1-P, (c) G2-F, (d) G2-M, (e) G2-P, and (f) G1 and (g) G2 control dendrons. The decrease in the size of the assemblies was observed only in presence of the complementary protein, extravidin. In the presence of the noncomplementary proteins, no size decrease was observed, although a tendency of the biotin-functionalized dendrimeric assemblies to increase in size, forming larger aggregates, was noticed for G1-F and G1-P (Figure 2a,b). This aggregate enlargement did not occur when control dendrons with PEG replacing the ligand moiety were exposed to the same proteins, even after 15 h (Figure 2f,g). This suggested that the higher density of biotin functionalization in the G1 dendrons is the likely reason for such aggregate enlargement.

Release of Encapsulated Guest Molecules

Next, we investigated the host capabilities of the dendrons in the presence and absence of complementary and noncomplementary proteins. Specifically, we were interested in assessing the effect of incorporating the ligand moieties at different locations within the dendrons upon the disassembly-induced guest release from the dendrimer host. To investigate these differences, Nile red was encapsulated in the micelle-like nanoassemblies, and its release was triggered as a consequence of binding-induced disassembly upon exposure to extravidin. Nile red is a hydrophobic molecule that exhibits reduced fluorescence in water, unless it is sequestered in a hydrophobic pocket. Therefore, the reduction in fluorescence is a good indicator of the binding-induced disassembly event. As shown in Figures 3a,b, small differences in placement of the ligand in the dendron produced rather different responses. First, percentage of released dye was assessed after exposing 25 μM solutions of the dendrons to increasing concentrations of extravidin. We noticed that G1-P and G2-P responded to increasing concentrations of extravidin more than other dendrons studied. To further evaluate this behavior, we monitored the release profiles over time for all the dendritic assemblies upon exposing these assemblies to 14 μM extravidin (Figure 3b). In the first hour, the release in the G2-P assembly was as high as 65%, while that of the G1-P assembly was around 22% increasing to 40% after 3 h. Interestingly, the release from the G1-P assembly ultimately reached about 77%, which is comparable with the 81% observed for G2-P. In comparison, similar exposures to extravidin resulted in 35%, 25%, and 13% for G1-F, G2-F, and G2-M, respectively. The extent of release observed for G1-F and G2-F is consistent with our prior observation.[13a] Interestingly, the release from G2-M aggregates is comparable to the release percentages observed due to nonspecific interactions.
Figure 3

Fluorescence-based percentage of release of guest molecules from the amphiphilic dendritic assemblies: (a) with increasing concentration of extravidin, (b) time-dependent release in the presence of 14 μM of extravidin, (c) in buffer solutionin the absence of proteins, (d) in assemblies based on control dendrons when exposed to extravidin, (e) with increasing concentration of Chy, (f) with increasing concentration of Pep, (g) with increasing concentration of Myo. (h) Absorption-based percentage of release upon interaction of dendrons with 14 μM of Myo.

As a control experiment, the release of the Nile red from the dendritic assemblies was also monitored in the absence of any protein (Figure 3c). No discernible release (<10%) was seen in these dendritic assemblies. Similarly, control dendrons that lack the biotin ligand also did not exhibit appreciable dye release in the presence of extravidin (Figure 3d). These results show that the release profiles observed in Figure 3a,b are indeed due to the ligand–protein binding. Moreover, it is clear that among the second generation dendrons, G2-P assembly is the only one that releases the hydrophobic guests efficiently following the extravidin binding. To further test the selectivity in the systems toward the target protein, the biotin-functionalized dendrons were exposed to increasing concentrations of noncomplementary proteins, Chy, Pep, and Myo monitoring the change in Nile red fluorescence (Figure 3e–g). No significant change in the emission intensity was observed for any of the dendrons in the presence of Chy and Pep. However, while Myo did not exhibit any change in the fluorescence intensity in the G1-P, G2-M, and G2-P based assemblies, there was a significant change in fluorescence in the G1-F and G2-F based assemblies (Figure 3g and S3). Interestingly, these latter dendrons also exhibited much smaller release in response to extravidin (Figure 3b). It is noteworthy that Myo is a metalloprotein, and therefore the cofactors in metalloproteins could be simply quenching the fluorescence of the dye molecule without the need for releasing the contents from the amphiphilic assembly. In fact, such a phenomenon has been previously observed with polymer–surfactant coassemblies.[11c] Fluorescence-based percentage of release of guest molecules from the amphiphilic dendritic assemblies: (a) with increasing concentration of extravidin, (b) time-dependent release in the presence of 14 μM of extravidin, (c) in buffer solutionin the absence of proteins, (d) in assemblies based on control dendrons when exposed to extravidin, (e) with increasing concentration of Chy, (f) with increasing concentration of Pep, (g) with increasing concentration of Myo. (h) Absorption-based percentage of release upon interaction of dendrons with 14 μM of Myo. To test this possibility, we investigated the percentage of Nile red release in the presence of Myo by absorption spectroscopy. If it is simply a quenching phenomenon, there should be no change in the absorption spectrum since all dye molecules are still confined in the amphiphilic assembly without being released. Indeed, we noted that there was no change in the absorption spectrum over time, which suggested that the observed increase in release percentage, i.e., decrease in Nile red fluorescence, is likely due to quenching (Figure 3h). As a second step, we were interested in gaining insights into the mechanism for the observed fluorescence reduction in the presence of Myo, especially for G1-F and G2-F. Two limiting mechanisms are possible: (i) inherently different encapsulation stabilities among the dendritic assemblies, causing the dye to leak out of G1-F (or G2-F) and to move to the hydrophobic pockets in Myo, where the proximity between the metalloprotein cofactor and the dye molecule causes fluorescence quenching; or (ii) nonspecific interactions between assemblies formed by G1-F (or G2-F) and Myo to bring the encapsulated dye molecules in proximity to the metalloprotein cofactor favoring quenching. The mechanism (i) implies leakage of the cargo and “re-encapsulation” in the metalloprotein pockets. To test this possibility, we used a recently reported polymeric nanogel that has been well-established to have cross-link density-dependent encapsulation stabilities.[23] Nile red-encapsulated nanogels with 0%, 20%, and 50% cross-link densities were exposed to Myo. If mechanism (i) was possible, there should be a cross-link density-dependent leakage based quenching. In all cases, we found that the extent of quenching was quite independent of the cross-link density (Figure S4a,b). The quenching however was found to increase with increasing concentration of the Myo, as observed with G1-F and G2-F assemblies. These results are taken to suggest that, since the 0% cross-linked nanogel assembly is quite leaky, Myo itself is not capable of “encapsulating” Nile red molecules in its hydrophobic pockets. This rules out mechanism (i). The mechanism (ii) can be due to dynamic or static quenching. Since dynamic and static quenching differs in their temperature dependence, we exposed the assemblies G1-F and G2-F encapsulating Nile red to increasing concentrations of Myo at two different temperatures. The number of collisions between the dendritic assemblies and the protein was expected to increase at a higher temperature as in dynamic quenching, leading to an increase in quenching.[24a] On the other hand, in the case of static quenching, the nonspecific and weak dendritic assembly–protein complex would dissociate at a higher temperature, leading to a decrease in quenching.[24b]Figure S4c,d shows the Stern–Volmer plots for G1-F and G2-F, respectively, at 25 and 38 °C. The fact that quenching decreases with temperature supports the nonspecific and weakly bound complex hypothesis. Thus, it is interesting to note that G1-F and G2-F exhibit higher nonspecific interaction and, at the same time, do not exhibit significant release of the guest molecules in response to the target extravidin. On the other hand, the dendrons G1-P and G2-P that exhibited the highest release in response to extravidin did not exhibit any nonspecific quenching with Myo. These data and those from the previous sections indicate that if the release of hydrophobic guests is triggered by the specific extravidin–biotin interaction, the selective binding with the complementary protein can occur more easily when a biotin ligand is grafted at the periphery than when it is grafted at the middle layer or at the focal point of the dendron and that this makes the release of hydrophobic guests faster for G2-P than for the other constructs, as observed in Figure 3a,b. If disassembly is controlled by the same interaction, since the biotin–extravidin binding is considered irreversible, when extravidin binds to a dendron in the assembly a reduction in the size of assembly will occur after some time, while the system reorganizes and equilibrates in smaller size assemblies, as observed by DLS (Figure 2). Thus, although all systems reorganized into smaller size assemblies, the release was higher for the dendron systems with a ligand located at the periphery. The reason for this disparity was not entirely clear to us. Based on the release kinetics, which was also faster for the dendritic systems with a ligand at the periphery, we hypothesized that in these cases the assembly reorganization was drastic enough to produce a higher release. On the other hand, in the cases with a ligand at the middle layer and focal point the assemblies rearranged slower into smaller size structures, allowing for the encapsulated hydrophobic small molecules to still be accommodated in hydrophobic pockets. In fact, Figure 3b shows that even after 6 h of G2-M (25 μM) exposure to a constant concentration of extravidin (14 μM), the release was as low as 13%. At this point, some interesting questions, including why is the release from G2-P so high compared to the release from G2-M and G2-F remain still open. For this reason, we have employed MD simulations to gain additional insights into this process.

MD Simulations of the Single Dendrons in Solution

We have used MD simulations focusing on the second generation dendrons in order to understand why the different positioning of a biotin ligand within the dendron scaffold has such a strong effect on the final properties. In particular, we were interested in gaining insight on ligand accessibility from the external solution if it is tethered at the dendron periphery, at the focal point, or at the middle layer. First, it was important to understand what these molecules “look like” in solution. In fact, it is known that similar dendritic structures can undergo strong folding in solution,[25] so that if the biotin ligand is backfolded and surrounded by PEG in the experimental conditions the specific binding with extravidin will be unlikely. The entire simulation work was carried out with the AMBER 12 suite of programs.[26] Molecular models were created with three different functionalization points for biotin to understand how the individual dendrons arrange in solution. G2-P, G2-F, and G2-M dendron models were created and parametrized according to our similar studies on dendrons interacting with proteins.[27] Starting configurations of the dendrons were then immerged in a simulation box (Figure 4a) containing explicit water molecules (Figure S6 and details in the Supporting Information (SI)). All systems underwent 200 ns of MD in periodic boundary conditions at 25 °C (298 K) of temperature and 1 atm of pressure. During this time, all dendrons reached equilibrium with good stability. The root-mean-square displacement (RMSD) and the radius of gyration (Rg) data extracted from the MD simulations were used to assess the system’s equilibration. Computational details for the simulation procedure and data analysis are available in the SI.
Figure 4

MD simulation of G2 dendrons in water. (a) Starting configuration of G2-P immerged in a water box. G2-P scaffold is colored in black, hydrophobic decyl chains in red, hydrophilic PEG in blue, and biotin ligand in green. Oxygen atoms of water molecules are represented as transparent cyan spheres. (b) Final snapshot taken from the MD simulation of G2-P. During the MD run all dendrons undergo strong folding in solution assuming a globular shape (Rg = 9.4–10 Å). (c) RDF plots for the G2 dendrons. (d) Simplified model: the dendron–extravidin interaction is composed of two steps: first, unfolding of the biotin ligand (green) that becomes available for the protein at the dendron surface, and second, the specific binding between biotin and extravidin.

Figure 4a shows the starting configuration of G2-P immerged in water. Initially, all dendrons were constructed with all PEG chains on one side and all decyl chains on the other side. Molecular models for the G2 dendrons were also constructed with alternated chains to avoid configuration-dependent results. However, MD simulations of those systems suggested that the initial configuration does not have any impact on the shape and equilibrium configuration assumed by such small and flexible molecules in water, in terms of density distribution, radius of gyration, etc. (see Figure S7 and SI for details). These findings are consistent with the idea of treating these dendrons as facially amphiphilic structures.[28] MD simulation of G2 dendrons in water. (a) Starting configuration of G2-P immerged in a water box. G2-P scaffold is colored in black, hydrophobic decyl chains in red, hydrophilic PEG in blue, and biotin ligand in green. Oxygen atoms of water molecules are represented as transparent cyan spheres. (b) Final snapshot taken from the MD simulation of G2-P. During the MD run all dendrons undergo strong folding in solution assuming a globular shape (Rg = 9.4–10 Å). (c) RDF plots for the G2 dendrons. (d) Simplified model: the dendron–extravidin interaction is composed of two steps: first, unfolding of the biotin ligand (green) that becomes available for the protein at the dendron surface, and second, the specific binding between biotin and extravidin. The size of G2 in water, predicted by MD simulation, does not change substantially depending on the tethering position of biotin; the Rg in the three cases is 9.4–10 Å (Figure 4b). In general, the dendrons tend to compact the hydrophobic decyl chains at the core and to surround them with hydrophilic PEG. On the other hand, in terms of distribution of the biotin ligand the situation is different. The plots in Figure 4c report the radial distribution function g(r) of the biotin ligand calculated with respect to the dendrons center and expressed as a function of the dendron radius (Rg) for the cases where biotin is grafted at the periphery (G2-P: red), the middle layer (G2-M: blue), or at the focal point (G2-F: black). In general, the g(r) values give indication on the relative probability to find the biotin ligand at a certain distance from the dendron center, being the position of the g(r) maximum peak the most probable one. The biotin density going from the center to the surface is calculated at each simulation step, and the reported g(r) data are averaged in time over the equilibrated phase MD trajectories (the last 100 ns). Thus, high and sharp peaks in g(r) identify high biotin density regions, but they also indicate high localization, confinement and backfolding (namely, atoms that cannot move are counted at each step in the same region of space). On the contrary, flexible and fluctuating groups will have low and broad g(r) peaks. Figure 4c shows that at the equilibrium biotin distribution is very different for G2-P, G2-M, and G2-F. In particular, the biotin g(r) maximum peak for G2-P (red curve) corresponds well with its Rg indicating that, on average, the ligand availability on the surface is very good. On the contrary, the maximum peaks of black and blue g(r) curves at a distance r ∼ 0.5 Rg suggest that biotin is considerably more backfolded in the case of G2-F and G2-M. These data give indication on how much the biotin ligand is available at the surface of the dendron and thus also on the probability to have a specific binding with extravidin. In fact, we built a simplified model describing the dendron–extravidin specific interaction as composed of two phases, namely, the unfolding of the biotin ligand to make it available for the protein and then the specific biotin–avidin interaction (Figure 4d). According to this scheme, we can obtain information on the overall affinity of G2-P, G2-M, and G2-F for extravidin by evaluating the free energy of the dendron-protein binding process as: ΔAbind = ΔAspecific + ΔAunfold; where the specific biotin–avidin affinity is known experimentally (absolute free energy of binding ΔAspecific = −20.4 kcal mol–1)[29a,29b] and can be considered as a constant for all cases (all dendrons bear the same number of biotin ligands), ΔAunfold is the free energy necessary to drag out the biotin ligand from its backfolded state to make it available at the dendron surface. ΔAunfold values for the different cases can be extracted directly from the g(r)[29c] and depend on how much the ligand is backfolded within the dendron structure (details in the SI). In particular, if G2 dendrons are thought of as spheres with radius Rg, ΔAunfold = 0 for G2-P as in this case the biotin ligand is most probably available at the dendron surface (g(r) peak position coincides with Rg). On the other hand, for G2-F and G2-M ΔAunfold is an unfavorable term (>0) that depends on the level of backfolding, namely, the lower the biotin ligand availability at the surface (high level of backfolding), the higher the necessary free energy to make it available at the dendrons’ surface for avidin binding. Table 1 reports the free energy of binding values calculated for the different dendrons.
Table 1

Free Energy Values for the Dendron–Extravidin Interactiona

dendronΔAunfoldbΔAspecificcΔAbinddstatistical weighte
G2-P0–20.4–20.41
G2-M1.40–20.4–19.00.097
G2-F1.38–20.4–19.020.10

Data are expressed in kcal mol–1.

ΔAunfold depends on the level of ligand backfolding. It can be calculated for all cases directly from the g(r); see SI.

Experimental biotin–avidin affinity – ΔAspecific = −20.4 kcal mol–1.[28]a,b

Global dendron–avidin binding affinity was calculated as: ΔAbind = ΔAspecific + ΔAunfold.

Relative probability for the dendron–avidin specific interaction (statistical weight) can be calculated as exp(−ΔAbindkbT–1).

We also calculated the statistical weight for the different dendrons from the dendron–avidin affinity energies as exp(−ΔAbindkbT–1) (kbT = 0.593 kcal mol–1 at room temperature), which provides qualitative indication on the relative probability for a specific binding with extravidin in the case of G2-P, G2-M, and G2-F, depending on how much the biotin ligand is available at the surface. In particular, these data suggest that if we set the probability of having extravidin specifically bound to the biotin moiety in G2-P to 1, then the probability of having extravidin bound to the biotin moiety in G2-F and G2-M is reduced to ∼0.1 due to reduced ligand availability; a difference of 1 order of magnitude indicating that the extravidin specific binding to biotin ligands from the dendrons at the surface of G2-P aggregates is 10 times higher than in the case of G2-F and G2-M. We will come back to this important point in the next sections. Data are expressed in kcal mol–1. ΔAunfold depends on the level of ligand backfolding. It can be calculated for all cases directly from the g(r); see SI. Experimental biotin–avidin affinity – ΔAspecific = −20.4 kcal mol–1.[28]a,b Global dendron–avidin binding affinity was calculated as: ΔAbind = ΔAspecific + ΔAunfold. Relative probability for the dendron–avidin specific interaction (statistical weight) can be calculated as exp(−ΔAbindkbT–1).

Modeling the Specific Binding between G2-P and Avidin

Then we aimed at studying the specific binding between one dendron and extravidin. Since we did not have any information regarding the conformation assumed by the dendron during the binding with extravidin, we started from the unfolded configuration of one dendron (Figure 5d). Initially, the biotin ligand of one unfolded G2-P was superimposed to the native biotin present in chain A of the crystallographic structure for the extravidin tetramer (PDB: 1AVD) replacing it and thus obtaining the first G2-P+AVD molecular complex.
Figure 5

Dendron-extravidin specific binding. (a) When biotin is tethered at the periphery, it is more available for extravidin binding. On the contrary, when it is grafted at the middle layer (b) or at the focal point (c), it is surrounded by an environment rich in PEG, so that the deep penetration of the ligand in one of the binding pockets of extravidin is hindered. (d) Starting configuration of the G2-P+AVD molecular system. G2-P backbone is colored in black, hydrophobic decyl chains in red, hydrophilic PEG in blue, and biotin in green. Oxygen atoms of water molecules and Cl– counterions are represented as transparent cyan and purple spheres, respectively. Extravidin is represented as a black ribbon. (e) Last equilibrated snapshot taken from the MD simulation of the G2-P+AVD system. During the MD run the specific binding between biotin and extravidin remained very stable. Extravidin binding pocket is colored in yellow, and water molecules and counterions are not shown for clarity.

Dendron-extravidin specific binding. (a) When biotin is tethered at the periphery, it is more available for extravidin binding. On the contrary, when it is grafted at the middle layer (b) or at the focal point (c), it is surrounded by an environment rich in PEG, so that the deep penetration of the ligand in one of the binding pockets of extravidin is hindered. (d) Starting configuration of the G2-P+AVD molecular system. G2-P backbone is colored in black, hydrophobic decyl chains in red, hydrophilic PEG in blue, and biotin in green. Oxygen atoms of water molecules and Cl– counterions are represented as transparent cyan and purple spheres, respectively. Extravidin is represented as a black ribbon. (e) Last equilibrated snapshot taken from the MD simulation of the G2-P+AVD system. During the MD run the specific binding between biotin and extravidin remained very stable. Extravidin binding pocket is colored in yellow, and water molecules and counterions are not shown for clarity. When biotin is grafted at the middle layer (Figure 5b) or the focal point (Figure 5c) it is surrounded by an environment rich in PEG. This makes the specific binding with extravidin, which implies the deep penetration of the ligand inside one of the four binding pockets of tetrameric extravidin, very difficult compared to G2-P (Figure 5a). This difficulty was reflected also on modeling. In fact, it was not possible to create the initial complex with extravidin for G2-F and G2-M without incurring intractable distortions. For this reason and in light of the results on biotin availability discussed in the previous section, for what pertains to the modeling study of the interactions with the complementary protein, we focused only on G2-P. The starting configuration of the G2-P+AVD complex (Figure 5d) was immerged in a periodic box containing explicit water molecules and the minimum number of counterions necessary to guarantee the system neutrality (details in the SI). G2-P+AVD system was equilibrated for 200 ns of NPT MD simulation at 25 °C (298 K) and 1 atm of pressure. This time was sufficient to reach the equilibrium. During the MD run, the binding between biotin and extravidin remained very strong and stable (see SI for details). The interaction energy (ΔEbind) between G2-P and extravidin was extracted directly from the MD trajectories according to the MM-PBSA approach.[30] In particular, ΔEbind is the sum of the gas-phase in vacuo interaction energy (ΔEgas) and the solvation term (ΔEsol). In general, the more negative the ΔEbind value, the stronger the binding. In the case of a 1:1 binding between G2-P and extravidin, ΔEbind was calculated to be as strong as −93.9 ± 8.7 kcal mol–1 at the equilibrium.

Modeling G2-P Self-Assembly and Extravidin Binding-Induced Disassembly

In our preliminary communication, we had suggested that the change in hydrophilic–lipophilic balance upon protein binding likely drives the binding-induced disassembly.[13a] To further test this hypothesis, we first modeled the amphiphilic assembly of G2-P dendrons in water. Then, we aimed at understanding how the dendron assembly changes upon extravidin addition to the system. MD simulation allows for the study of hydrophobic aggregation in water, as it was recently reported in the case of the self-assembly of hydrophobic drug molecules and amphiphilic diblock copolymer micelles.[31] According to the same protocol, nine copies of G2-P dendron arranged on a plane were immerged in a simulation box filled with explicit water molecules. In particular, as a starting configuration for the G2-P dendron (Figure 6a), we chose the final equilibrated configuration produced by the simulation of G2-P+AVD binding (Figure 5e), with the biotin ligand initially extended in extravidin-bound conformation, so that it was also possible to add a bound extravidin to the central dendron without difficulties. We thus generated two molecular systems both containing nine G2-P dendrons solvated in water, where the central dendron is respectively unbound (Figure 6a: 9G2-P) or bound to one extravidin protein (Figure 6c; 9G2-P+AVD), in such a case the central dendron was simply replaced by the G2-P+AVD complex from Figure 5e.
Figure 6

G2-P self-assembly and extravidin-induced disassembly. (a) Initially, nine copies of unbound G2-P dendrons were immerged in a simulation box containing water (9G2-P system). (b) During 200 ns of MD simulation the nine dendrons aggregate and rearrange: hydrophilic PEG (blue) and biotins (green) surround the hydrophobic parts (red decyl chains). (c) The starting configuration for the G2-P+AVD system is the same of G2-P, but a tetrameric extravidin protein (AVD: black ribbons) is bound to the central G2-P dendron. Cl– ions (purple) were also added for neutralization. (d) The 9G2-P+AVD system was equilibrated for 200 ns of MD simulation, during which the specific binding between biotin (yellow) and AVD remains stable and the nine G2-P dendrons in the system self-assemble. (e) Rg plots for the 9G2-P assembly case (black curve). Rg plots are obtained also for the biotin ligands (green), PEG (blue), and hydrophobic decyl chains (red) to understand molecular reorganization in solution. (f) The central dendron (b,d: red) was used as a reference, and the self-assembly energy (ΔEass) was extracted for the MD simulations. ΔEass values were obtained for the 9G2-P (black) and 9G2-P+AVD (blue) systems and plotted as a function of simulation time. (g) The energetic difference ΔΔEass (red) identifies the G2-P self-assembly destabilization due to AVD binding.

Both systems were equilibrated for 200 ns of MD simulation at 25 °C and 1 atm pressure. During the simulation of the 9G2-P system, we found that the nine dendrons by themselves showed a strong propensity to self-assemble. The aggregate also reorganized in order to decrease the exposure of the hydrophobic regions (red) to the external solution. In general, Figure 6b shows that the red decyl chains converge at the core of the aggregate and are surrounded by blue PEG chains. Biotin ligands (green) are well exposed to the surface. Structural reorganization during self-assembly is also demonstrated by the Rg plots in Figure 6e obtained from the MD simulation of the 9G2-P system. Red decyl chains converge to a lower Rg value than PEG (blue) and biotin (green), indicating structural reorganization within the aggregate due to hydrophobic effects. G2-P self-assembly and extravidin-induced disassembly. (a) Initially, nine copies of unbound G2-P dendrons were immerged in a simulation box containing water (9G2-P system). (b) During 200 ns of MD simulation the nine dendrons aggregate and rearrange: hydrophilic PEG (blue) and biotins (green) surround the hydrophobic parts (red decyl chains). (c) The starting configuration for the G2-P+AVD system is the same of G2-P, but a tetrameric extravidin protein (AVD: black ribbons) is bound to the central G2-P dendron. Cl– ions (purple) were also added for neutralization. (d) The 9G2-P+AVD system was equilibrated for 200 ns of MD simulation, during which the specific binding between biotin (yellow) and AVD remains stable and the nine G2-P dendrons in the system self-assemble. (e) Rg plots for the 9G2-P assembly case (black curve). Rg plots are obtained also for the biotin ligands (green), PEG (blue), and hydrophobic decyl chains (red) to understand molecular reorganization in solution. (f) The central dendron (b,d: red) was used as a reference, and the self-assembly energy (ΔEass) was extracted for the MD simulations. ΔEass values were obtained for the 9G2-P (black) and 9G2-P+AVD (blue) systems and plotted as a function of simulation time. (g) The energetic difference ΔΔEass (red) identifies the G2-P self-assembly destabilization due to AVD binding. Similar aggregation occurred also in the 9G2-P+AVD system containing eight G2-P unbound dendrons surrounding a single G2-P+AVD complex (Figure 6c,d). However, interestingly in this case the aggregate was clearly less tight and less ordered at the equilibrium than in the 9G2-P case. We were interested in evaluating the native propensity of G2-P dendrons to self-assemble and in understanding how much the latter is affected when one specific binding with a complementary protein occurs at the surface of the aggregate. Thus, we obtained from the MD simulations the dendrons self-assembly energy (ΔEass) for the two systems. The latter was evaluated as the binding energy between the central dendron with the other eight surrounding ones in the systems (schemes in Figure 6b,d: red and blue dendrons, respectively).[15c] ΔEass was calculated as described in the previous section, taking into account for solute–solute and solute–solvent interactions (details on energetic analysis in the SI). ΔEass calculated values were further normalized per dendron (i.e., divided per 8), so that they become general indicators of how stable is the aggregation within the G2-P aggregates. Data reported in Figure 6f show that, at the equilibrium, the self-assembly energy triggering aggregation is as high as ΔEass = −32.5 ± 1.4 kcal mol–1 in absence of extravidin binding (9G2-P: black curve). On the other hand, in the 9G2-P+AVD system ΔEass is reduced to −16.8 ± 1.4 kcal mol–1 (blue curve). The destabilization induced by the specific binding with extravidin is represented by the energetic difference ΔΔEass = ΔEass (9G2-P+AVD) – ΔEass (9G2-P). Positive values for ΔΔEass indicate that, in general, when an AVD protein binds specifically to one biotin ligand at the surface of a G2-P aggregate, the self-assembly is locally destabilized, which could be the trigger for disassembly. In our case, the energetic destabilization due to extravidin binding converged to the value of ΔΔEass = +15.7 ± 1.7 kcal mol–1 (Figure 6g, red curve), indicating that when one AVD binding occurs, the stability of G2-P dendron self-assembly (ΔEass) is reduced to the half. Specific and nonspecific interaction at the interface between AVD and G2-P aggregates. (a,b) During the MD simulation the dendrons self-assemble and blue dendrons surround the red one, which is specifically bound to AVD (biotin in green). (c) The dendrons can interact with extravidin specifically (red) and also nonspecifically (blue). (d) Specific (red) and nonspecific (blue) interaction energies normalized per dendron expressed as a function of simulation time. (e) Nonspecific interactions with AVD (blue) are sensibly lower than the dendrons native self-assembly energy (green, the same 9G2-P values reported in black in Figure 6f), demonstrating why nonspecific interactions alone are energetically not strong enough to trigger disassembly. We also extracted useful information from the MD simulations regarding specific and nonspecific interactions (Figure 7). In fact, during the 200 ns of MD simulation of the 9G2-P+AVD system, we saw that all nine G2-P dendrons interact not only with each other but also with extravidin (Figure 7a–c). Importantly, while the central dendron (red) is specifically bound to the protein through the biotin ligand, the other eight interact with AVD nonspecifically. To obtain insights into the differences between specific and nonspecific interactions with AVD, we extracted the binding energy (ΔEbind) between the dendron aggregate and extravidin from the MD simulation of 9G2-P+AVD through the same approach adopted previously. ΔEbind measures the global interaction of the aggregate (composed of nine dendrons) for AVD, and it was calculated as ΔEbind = −162.5 ± 14.3 kcal mol–1, equivalent to ΔEbind = −18.1 ± 1.6 kcal mol–1 per dendron, at the equilibrium. Thus, this result shows that in general extravidin is attracted by G2-P aggregates, which suggests that as soon as extravidin is added to a solution containing dendron aggregates there will be a long-range trigger for molecular recognition.
Figure 7

Specific and nonspecific interaction at the interface between AVD and G2-P aggregates. (a,b) During the MD simulation the dendrons self-assemble and blue dendrons surround the red one, which is specifically bound to AVD (biotin in green). (c) The dendrons can interact with extravidin specifically (red) and also nonspecifically (blue). (d) Specific (red) and nonspecific (blue) interaction energies normalized per dendron expressed as a function of simulation time. (e) Nonspecific interactions with AVD (blue) are sensibly lower than the dendrons native self-assembly energy (green, the same 9G2-P values reported in black in Figure 6f), demonstrating why nonspecific interactions alone are energetically not strong enough to trigger disassembly.

How much of this attraction is due to the specific interaction (red dendron) and how much to nonspecific interactions? We obtained this information via a simple decomposition of the global ΔEbind energy on a per dendron basis. This analysis shows that most of the interaction is due to the specific binding between the red dendron in Figure 7 and AVD, that is as strong as ΔEbind = −91.6 ± 5.4 kcal mol–1. Such a high ΔEbind value is very close to that found for the 1:1 G2-P+AVD binding in Figure 5, which demonstrates that the presence of the other blue dendrons in the aggregate does not affect the strength of the specific binding. On the other hand, nonspecific interactions with extravidin are on average sensibly lower (ΔEbind = −8.8 ± 1.9 kcal mol–1). This result shows that in this case the strength of nonspecific interactions is 1 order of magnitude weaker than specific interactions and also less persistent and more discontinuous, as demonstrated by the standard deviation which is ∼20% of the ΔEbind average value (for specific interactions the latter is ∼6%). Interestingly, nonspecific interactions, even if present, are weaker than the intrinsic self-assembly energy of G2-P dendrons (Figure 7e: ∼9 vs. ∼33 kcal mol–1). For this reason, they are probably not strong enough to perturb the stability of the aggregates, which is consistent with our DLS evidence (Figure 2) showing no disassembly in the case of noncomplementary proteins.

Multivalent AVD Binding

All these results suggest that as soon as extravidin gets in contact with the surface of a G2-P aggregate, the formation of a specific binding will be an energetically favored event. In addition, it is worth noting that extravidin is a protein tetramer possessing four binding sites for biotin. Thus, after a first specific binding occurs between AVD with one biotin at the G2-P aggregate surface, the protein can find also other biotin ligands available in the neighborhood to establish many more specific bindings in a cooperative way. According to the so-called multivalent effect, this will be an energetically favored thus highly probable event.[32] To test this hypothesis, we performed the MD simulation, similar to that in Figure 5, but with four G2-P dendrimers bound to all of the biotin binding sites of tetrameric AVD (Figure 8a,b). The result of the energetic analysis for this case is in agreement with the multivalency principle. In fact, the interaction energy (ΔEbind) between AVD and all four G2-P dendrons extracted from this simulation and reported in Figure 8b is sensibly higher than that related to the binding of a single dendron reported in Figure 5 (multiplied per four, which equals to: ΔEbind = −375.6 ± 34.9 kcal mol–1). This large difference of ∼137 kcal mol–1 is a clear signal of multivalency, and it demonstrates that AVD is energetically favored to bind more biotins at the same time.
Figure 8

Multivalent AVD binding. (a) An additional system was simulated with four G2-P dendrons (yellow, green, magenta, and red) specifically bound to the four biotin binding sites of extravidin (PDB: 1AVD). (b) Binding energy extracted from this MD simulation indicates that the specific binding of AVD to multiple G2-P dendrons is energetically favored.

Multivalent AVD binding. (a) An additional system was simulated with four G2-P dendrons (yellow, green, magenta, and red) specifically bound to the four biotin binding sites of extravidin (PDB: 1AVD). (b) Binding energy extracted from this MD simulation indicates that the specific binding of AVD to multiple G2-P dendrons is energetically favored. The consequences of this behavior can be important. In fact, this tendency of AVD to bind more biotin ligands at the same time, when available, can in principle speed up the disassembly process. Conceptually, since AVD will tend to preserve its structure much more than the dendron aggregate, in the case of cooperative binding to multiple ligands, it is reasonable to think that the dendrons from the aggregate will adapt over extravidin, rather than the contrary (unlikely protein collapsing over the G2-P surface). This is consistent with a picture where G2-P aggregates are progressively degraded.

Two Different Proposed Mechanisms for Disassembly

One key factor allowing for multivalent binding is biotin availability, i.e., biotin ligands must be accessible at the aggregate surface and free to complete specific interactions with AVD. In fact, as we already discussed with Figure 4c,d, the chance of having multivalent binding between the dendrons with AVD will be extremely sensitive to biotin availability at the surface of the aggregates. Our MD simulations suggested that biotin availability is high for G2-P and low for G2-M and G2-F (Figure 4c,d). Thus, at the G2-P aggregates surface AVD will find many accessible binding spots for completing specific interactions. Moreover, after a first specific binding is established, the same AVD protein will be then energetically favored to bind more biotins from other G2-P dendrons. On the other hand, since the probability to have specific binding for G2-M and G2-F is reduced, even if a first specific interaction occurs between AVD and one biotin ligand, it is reasonable to hypothesize that the chance of having multivalent AVD binding at the surface of G2-M and G2-F aggregates will be even lower. Proposed mechanisms for disassembly. (a) Due to high biotin availability at the surface of G2-P (red) aggregates, the formation of multivalent AVD-biotin bindings is an energetically favored phenomenon. While completing multiple bindings, AVD would trigger the fast disassembly of the aggregates into small aggregates (AVD micelles) via exfoliation, leading to the release of hydrophobic fluorescent guest molecules in solution (green). (b) On the other hand, even if a less favored specific binding event occurs at the surface of G2-M and G2-F aggregates (blue), multivalent binding of the same AVD to multiple biotins is unlikely due to low ligand availability. The dendron aggregates would disassemble in larger and more ordered assemblies that would be able to retain the guests in their interior. (c) Starting configuration of an additional system containing one AVD (black ribbons) with four specific bindings with four G2-P dendrons (red) and surrounded by other 36 G2-P dendrons (transparent blue). The number of G2-P dendrons that can surround one AVD being in contact with the surface was calculated as ∼40 according to the Mansfield–Tomalia–Rakesh equation.[32] (d) Final snapshot of the MD simulation of the large aggregate in solution (oxygen atoms of water molecules in transparent cyan and counterions in purple). During the simulation all G2-P dendrons (transparent black) collapse and surround AVD (black). Specifically bound dendrons are colored in red. (e) Gyration radius of the large aggregate over simulation time (black). During the MD run the dendrimers first aggregate around AVD and then undergo rearrangement, surrounding the hydrophobic decyl chains (red) with the hydrophilic PEG chains (blue). Rg plots for biotin ligands are represented in green. In light of these results, we propose two possible mechanisms for the observed supramolecular disassembly (schematized in Figure 9). A first one for the G2-P case based on multivalent binding of AVD, leading to the rapid disassembly of the dendron aggregates in solution due to exfoliation (Figure 9a). A second one for G2-M and G2-F, where the aggregates would disassemble more slowly due to the destabilization (Figure 6) induced by the specific binding with AVD (Figure 9b). In principle, the first mechanism would result in the rapid production of smaller aggregates (since they are limited by AVD size) and in a higher level of hydrophobic guests release. In fact, the velocity of the process and the small size of the aggregates would not allow for the structural rearrangement necessary to retain the guest molecules. The second proposed mechanism would most likely produce larger size aggregates and lower levels of release.
Figure 9

Proposed mechanisms for disassembly. (a) Due to high biotin availability at the surface of G2-P (red) aggregates, the formation of multivalent AVD-biotin bindings is an energetically favored phenomenon. While completing multiple bindings, AVD would trigger the fast disassembly of the aggregates into small aggregates (AVD micelles) via exfoliation, leading to the release of hydrophobic fluorescent guest molecules in solution (green). (b) On the other hand, even if a less favored specific binding event occurs at the surface of G2-M and G2-F aggregates (blue), multivalent binding of the same AVD to multiple biotins is unlikely due to low ligand availability. The dendron aggregates would disassemble in larger and more ordered assemblies that would be able to retain the guests in their interior. (c) Starting configuration of an additional system containing one AVD (black ribbons) with four specific bindings with four G2-P dendrons (red) and surrounded by other 36 G2-P dendrons (transparent blue). The number of G2-P dendrons that can surround one AVD being in contact with the surface was calculated as ∼40 according to the Mansfield–Tomalia–Rakesh equation.[32] (d) Final snapshot of the MD simulation of the large aggregate in solution (oxygen atoms of water molecules in transparent cyan and counterions in purple). During the simulation all G2-P dendrons (transparent black) collapse and surround AVD (black). Specifically bound dendrons are colored in red. (e) Gyration radius of the large aggregate over simulation time (black). During the MD run the dendrimers first aggregate around AVD and then undergo rearrangement, surrounding the hydrophobic decyl chains (red) with the hydrophilic PEG chains (blue). Rg plots for biotin ligands are represented in green.

This hypothesis is consistent with all our experimental evidence. In fact, the high level of biotin availability of G2-P (Figure 4c,d) is compatible with the fast disassembly and hydrophobic guests release by G2-P aggregates in presence of extravidin, as shown by our DLS and fluorescence experiments. At the same time, our data demonstrate that G2-M and G2-F aggregates also disassemble in presence of extravidin (Figure 2), but more slowly. In addition, the final size of the disassembled aggregates is larger than that of G2-P (∼14 VS ∼7 nm), and no appreciable guest release is present during the disassembly of the G2-M and G2-F aggregates (Figure 3). Finally, we performed another MD simulation aiming at representing the final step of the disassembly mechanism based on multivalent binding and exfoliation represented in Figure 9a. First, we calculated how many G2-P dendrons, thought of as small spheres, would be necessary to surround completely an AVD molecule, thought of as a larger sphere, while being in direct contact with the protein surface. Using the gyration radii of a single G2-P dendron and AVD obtained from our MD simulations (respectively, Rg = 0.97 and 2.2 nm on average), we found that in our case this ratio equals to ∼40 G2-P dendrons for a single AVD according to the Mansfield–Tomalia–Rakesh equation.[33] Then we started from the final configuration obtained from the simulation of AVD specifically bound to four G2-P dendrons and have added other 36 equilibrated G2-P dendrons surrounding the complex (Figure 9c). The system was then immerged in a simulation box containing explicit solvent molecules and simulated for 150 ns in NPT conditions, at 25 °C of temperature and 1 atm of pressure. During this time, all 40 G2-P dendrons collapsed over the surface of AVD surrounding the protein (Figure 9d) and forming a real protein micelle.[34] Furthermore, during the MD simulation the entire aggregate underwent structural reorganization attempting to surround the surface with hydrophilic PEG, limiting as much as possible the exposure of the hydrophobic parts to water. This structural rearrangement is testified by the Rg plots of Figure 9e, demonstrating how during the MD simulation the PEG chains surround the surface of the aggregate. We also calculated the equilibrated size of the whole aggregate of Figure 9d, which resulted to be Rg = 3.3 nm. Considering the compact and globular nature of this large aggregate, the hydrodynamic size can be calculated as Rh ∼ 1.29 Rg.[25a] For this case the obtained size is thus Dh ∼ 8.5 nm, which is in good agreement with the size obtained for the disassembled aggregates in the case of G2-P (Figure 2).

Conclusions

Our research reveals that ligand placement on a supramolecular scaffold for binding-induced disassembly greatly impacts disassembly and release of encapsulated guest molecules, as we have shown from the high release difference observed, for example, between G2-P and G2-M. The best place to attach a ligand, looking for a protein triggered release from a dendritic micelle-like nanostructure, is the periphery. MD simulations show backfolding of the ligand, when attached to middle layer and focal point; and a better availability for protein binding, when the ligand is attached at the periphery. Incorporation of the biotin ligand in the dendrons gives the dendritic assemblies selectivity toward the target protein extravidin, regardless of the ligand position. Nonetheless, ligand positioning in the dendrons gives the assemblies sensitivity toward release upon binding of a target protein. MD simulations show that after AVD binding the stability of G2-P assembly is strongly affected, possible signal of AVD binding-induced disassembly. In addition, our computational efforts show that once a first specific binding between AVD and one G2-P dendron occurs, the multivalent binding of the same protein to other G2-P dendrons via specific biotin–extravidin interaction is an energetically favored event. This evidence allowed us to hypothesize two different mechanisms of disassembly induced by AVD binding. The divergence of the mechanism also allows us to explain the observed experimental differences in guest release that depend on the ligand location in the scaffold. A fast one for G2-P based on high biotin availability at the aggregate surface, multivalent interactions, and aggregates exfoliation. This mechanism leads to fast formation of small disassembled aggregates and to a high release of hydrophobic guests. A second mechanism for G2-M and G2-P, slower, based on AVD binding-induced aggregate destabilization and producing larger and more ordered aggregates in solution that are still capable of retaining the guest molecules in their interior. Substitution of a PEG unit in the dendrimer with a pendant biotin increased nonspecific interactions of the assemblies with proteins, which was seen as the formation of larger aggregates in solution. This became more evident when the density of PEG chains on the hydrophilic face was low, as it is in a G1 dendron compared to a G2. This, in turn, facilitates the formation of a weak complex with proteins, which was evidenced when a metalloprotein acted as a quencher, generating static quenching of the encapsulated fluorophore molecules. The combination of these results shows that G2 dendrons are better for selectivity and that the periphery is the best location for achieving binding-induced disassembly. Thus, G2-P exhibits most guest release and least nonspecific interactions with other proteins. Our MD simulations of G2-P show that, despite the general affinity between the aggregates and AVD, nonspecific interactions alone are too weak, and a specific binding (more than 1 order of magnitude stronger) is needed to trigger the release of the hydrophobic guest molecules. The research reported here gives a picture of how supramolecular disassembly and release might be largely affected by choosing a specific location for a trigger, rather than a random placement based on molecular architecture. Also, we have shown how controlled variations in PEG density could affect interactions of nanoparticles with proteins. Developing zwitterionic moieties,[35] which have the potential to circumvent the steric hindrance exhibited by lengthy PEG chains in amphiphilic systems, for use in protein-binding-induced disassembly is part of the ongoing focus in our laboratories.
  77 in total

1.  Validation of a novel molecular dynamics simulation approach for lipophilic drug incorporation into polymer micelles.

Authors:  Aliya O Kasimova; Giovanni M Pavan; Andrea Danani; Karine Mondon; Andrea Cristiani; Leonardo Scapozza; Robert Gurny; Michael Möller
Journal:  J Phys Chem B       Date:  2012-04-02       Impact factor: 2.991

2.  Self-cross-linked polymer nanogels: a versatile nanoscopic drug delivery platform.

Authors:  Ja-Hyoung Ryu; Reuben T Chacko; Siriporn Jiwpanich; Sean Bickerton; R Prakash Babu; S Thayumanavan
Journal:  J Am Chem Soc       Date:  2010-11-15       Impact factor: 15.419

3.  Functionally layered dendrimers: a new building block and its application to the synthesis of multichromophoric light-harvesting systems.

Authors:  William R Dichtel; Stefan Hecht; Jean M J Fréchet
Journal:  Org Lett       Date:  2005-09-29       Impact factor: 6.005

4.  Dendrimer analogues of linear molecules to evaluate energy and charge-transfer properties.

Authors:  Arpornrat Nantalaksakul; Raghunath Reddy Dasari; Tai-Sang Ahn; Rabih Al-Kaysi; Christopher J Bardeen; S Thayumanavan
Journal:  Org Lett       Date:  2006-07-06       Impact factor: 6.005

5.  Modeling the multivalent recognition between dendritic molecules and DNA: understanding how ligand "sacrifice" and screening can enhance binding.

Authors:  Giovanni M Pavan; Andrea Danani; Sabrina Pricl; David K Smith
Journal:  J Am Chem Soc       Date:  2009-07-22       Impact factor: 15.419

6.  Guest-release control in enzyme-sensitive, amphiphilic-dendrimer-based nanoparticles through photochemical crosslinking.

Authors:  Krishna R Raghupathi; Malar A Azagarsamy; S Thayumanavan
Journal:  Chemistry       Date:  2011-09-02       Impact factor: 5.236

7.  Self-assembling nanoprobes that display off/on 19F nuclear magnetic resonance signals for protein detection and imaging.

Authors:  Yousuke Takaoka; Takashi Sakamoto; Shinya Tsukiji; Michiko Narazaki; Tetsuya Matsuda; Hidehito Tochio; Masahiro Shirakawa; Itaru Hamachi
Journal:  Nat Chem       Date:  2009-09-23       Impact factor: 24.427

8.  Site-specific installation and study of electroactive units in every layer of dendrons.

Authors:  Malar A Azagarsamy; Kothandam Krishnamoorthy; Kulandaivelu Sivanandan; S Thayumanavan
Journal:  J Org Chem       Date:  2009-12-18       Impact factor: 4.354

9.  Double click: dual functionalized polymeric micelles with antibodies and peptides.

Authors:  Dianna P Y Chan; Shawn C Owen; Molly S Shoichet
Journal:  Bioconjug Chem       Date:  2013-01-04       Impact factor: 4.774

10.  HPV16 oncoproteins induce MMPs/RECK-TIMP-2 imbalance in primary keratinocytes: possible implications in cervical carcinogenesis.

Authors:  Laura Beatriz da Silva Cardeal; Enrique Boccardo; Lara Termini; Tatiana Rabachini; Maria Antonieta Andreoli; Celso di Loreto; Adhemar Longatto Filho; Luisa Lina Villa; Silvya Stuchi Maria-Engler
Journal:  PLoS One       Date:  2012-03-16       Impact factor: 3.240

View more
  13 in total

Review 1.  Experimental and theoretical investigations in stimuli responsive dendrimer-based assemblies.

Authors:  Mijanur Rahaman Molla; Poornima Rangadurai; Giovanni M Pavan; S Thayumanavan
Journal:  Nanoscale       Date:  2015-03-07       Impact factor: 7.790

2.  Smart Organic Two-Dimensional Materials Based on a Rational Combination of Non-covalent Interactions.

Authors:  Wei Bai; Ziwen Jiang; Alexander E Ribbe; S Thayumanavan
Journal:  Angew Chem Int Ed Engl       Date:  2016-08-04       Impact factor: 15.336

3.  A supramolecular dissociation strategy for protein sensing.

Authors:  Hui Wang; Jiaming Zhuang; Krishna R Raghupathi; S Thayumanavan
Journal:  Chem Commun (Camb)       Date:  2015-12-18       Impact factor: 6.222

4.  Protein-induced supramolecular disassembly of amphiphilic polypeptide nanoassemblies.

Authors:  Mijanur Rahaman Molla; Priyaa Prasad; S Thayumanavan
Journal:  J Am Chem Soc       Date:  2015-06-02       Impact factor: 15.419

5.  Role of Aromatic Interactions in Temperature-Sensitive Amphiphilic Supramolecular Assemblies.

Authors:  Oyuntuya Munkhbat; Matteo Garzoni; Krishna R Raghupathi; Giovanni M Pavan; S Thayumanavan
Journal:  Langmuir       Date:  2016-03-18       Impact factor: 3.882

Review 6.  Supramolecular Self-Associations of Amphiphilic Dendrons and Their Properties.

Authors:  Evgeny Apartsin; Anne-Marie Caminade
Journal:  Chemistry       Date:  2021-10-29       Impact factor: 5.020

7.  Stimuli-responsive nanomaterials for biomedical applications.

Authors:  Angela P Blum; Jacquelin K Kammeyer; Anthony M Rush; Cassandra E Callmann; Michael E Hahn; Nathan C Gianneschi
Journal:  J Am Chem Soc       Date:  2015-02-06       Impact factor: 15.419

8.  Consequences of chirality on the dynamics of a water-soluble supramolecular polymer.

Authors:  Matthew B Baker; Lorenzo Albertazzi; Ilja K Voets; Christianus M A Leenders; Anja R A Palmans; Giovanni M Pavan; E W Meijer
Journal:  Nat Commun       Date:  2015-02-20       Impact factor: 14.919

9.  Molecular glues for manipulating enzymes: trypsin inhibition by benzamidine-conjugated molecular glues.

Authors:  Rina Mogaki; Kou Okuro; Takuzo Aida
Journal:  Chem Sci       Date:  2015-03-18       Impact factor: 9.825

10.  Supramolecular disassembly of facially amphiphilic dendrimer assemblies in response to physical, chemical, and biological stimuli.

Authors:  Krishna R Raghupathi; Jing Guo; Oyuntuya Munkhbat; Poornima Rangadurai; S Thayumanavan
Journal:  Acc Chem Res       Date:  2014-06-17       Impact factor: 22.384

View more

北京卡尤迪生物科技股份有限公司 © 2022-2023.