Simulations of Dendrimers with Flexible Spacer Chains and Explicit Counterions under Low and Neutral pH Conditions

Binding mechanisms in dendrimer-surfactant complexes

Molecular dynamics simulations were employed to investigate the impact of interactions between dendritic polyeclectrolytes and amphiphilic surfactants on the supramolecular complex formation. We recognize two crucial parameters that govern association of surfactants within dendrimers: surfactant hydrophobicity, ε^{*}, and dendrimer generation, G. We find that depending on the values of ε^{*} and G encapsulation of surfactants by dendrimers is either noncooperative or cooperative. The noncooperative binding is characterized by absorption of surfactants as unimers, whereas in cooperative binding absorption of unimers is followed by aggregate formation through hydrophobic attractions between the surfactant tails. Our results provide guidelines for controlled encapsulation of guest molecules in dendrimer-based guest-host complexes.

Spatial segregation of mixed-sized counterions in dendritic polyelectrolytes

Langevin dynamics simulations are utilized to study the structure of a dendritic polyelectrolyte embedded in two component mixtures comprised of conventional (small) and bulky counterions. We vary two parameters that trigger conformational properties of the dendrimer: the reduced Bjerrum length, [Formula: see text], which controls the strength of electrostatic interactions and the number fraction of the bulky counterions, [Formula: see text], which impacts on their steric repulsion. We find that the interplay between the electrostatic and the counterion excluded volume interactions affects the swelling behavior of the molecule. As compared to its neutral counterpart, for weak electrostatic couplings the charged dendrimer exists in swollen conformations whose size remains unaffected by [Formula: see text]. For intermediate couplings, the absorption of counterions into the pervaded volume of the dendrimer starts to influence its conformation. Here, the swelling factor exhibits a maximum which can be shifted by increasing [Formula: see text]. For strong electrostatic couplings the dendrimer deswells correspondingly to [Formula: see text]. In this regime a spatial separation of the counterions into core-shell microstructures is observed. The core of the dendrimer cage is preferentially occupied by the conventional ions, whereas its periphery contains the bulky counterions.

Dendritic polyelectrolytes revisited through the Poisson-Boltzmann-Flory theory and the Debye-Hückel approximation

The properties of a dendritic polyelectrolyte in equilibrium with a reservoir of monovalent salts are investigated using the cell model and the Poisson-Boltzmann-Flory theory. Within this approach we use the Debye-Hückel approximation to solve the Poisson-Boltzmann equation and minimize the semi-grand potential of the system with respect to the size of the molecule which enables us to inspect its conformations as well as the electric field, the ionic density profile, the overall charge density, the effective charge of the dendrimer and the osmotic pressure based on their response to the salt concentration and the dendrimer charge. The model predicts pronounced trapping of salt ions, a local charge neutrality and a zero electric field in the volume of the molecule as well as oscillations of the density profiles and the electric field in the vicinity of the dendrimer-bulk interface. As a result of ion trapping and screening of Coulomb interactions monovalent salts are found to have a minor effect on the size of the dendrimer. Specifically, the dendrimer exists in slightly swollen states as compared to the neutral molecule which indicates that the conformational properties of the polyelectrolyte depend weakly on monovalent salts. These observations harmonise with the equilibrium behavior of the dendrimer pressure, the internal pressure and the bulk pressure, respectively.

Dendritic polyelectrolytes as seen by the Poisson–Boltzmann–Flory theory

The conformational and electrostatic properties of dendritic polyelectrolytes accompanied by counterions are investigated using the Poisson–Boltzmann–Flory theory.

Computer Simulation Studies on the pH-Responsive Self-Assembly of Amphiphilic Carboxy-Terminated Polyester Dendrimers in Aqueous Solution

This paper investigates the pH-responsive self-assembly of an amphiphilic carboxyl-terminated polyester dendrimer, H20-COOH, in aqueous solution using the dissipative particle dynamics method. The electrostatic interactions were described by introducing the explicit interaction between the smeared charges on ionized polymer beads and the counterions. The results show that the self-assemblies could change from unimolecular micelles, microphase-separated small micelles, wormlike micelles, sheetlike micelles, and small vesicles to large vesicles with the decrease in the degree of ionization (α) of carboxylic acid groups. In addition, the detailed self-assembly mechanisms and the molecular packing models have also been disclosed for each self-assembly stages. Interestingly, the wormlike micelles are found to change from linear to branched when α decreases from 0.182 to 0.109. The current work might serve as a comprehensive understanding on the effect of carboxylic acid groups on the self-assembly behaviors of dendritic polymers.

Polyelectrolyte brushes in external fields: molecular dynamics simulations and mean-field theory

A mean-field model is developed to predict the layer-thickness of sparse and salt-free polyelectrolyte brushes, exposed to an external electric field, which attracts the polyelectrolytes to the substrate. In molecular dynamics simulations, it is shown that a fraction of polymers collapses entirely to screen the charge of the substrate. The remaining brushes are then treated as field-free brushes at reduced grafting density. The mean-field model may thus be applied to field-free brushes, both in their osmotic and their weak charge regimes. It yields simple, closed equations for the brush height and for the partition of counterions in- and outside the brushes, and accurately reproduces simulation data of the collapse of the brushes during the crossover between both charge regimes.

Adsorption of branched and dendritic polymers onto flat surfaces: a Monte Carlo study

Using Monte Carlo simulations based on the bond fluctuation model we study the adsorption of starburst dendrimers with flexible spacers onto a flat surface. The calculations are performed for various generation number G and spacer length S in a wide range of the reduced temperature τ as the measure of the interaction strength between the monomers and the surface. Our simulations indicate a two-step adsorption scenario. Below the critical point of adsorption, τc, a weakly adsorbed state of the dendrimer is found. Here, the dendrimer retains its shape but sticks to the surface by adsorbed spacers. By lowering the temperature below a spacer-length dependent value, τ*(S) < τc, a step-like transition into a strongly adsorbed state takes place. In the flatly adsorbed state the shape of the dendrimer is well described by a mean field model of a dendrimer in two dimensions. We also performed simulations of star-polymers which display a simple crossover-behavior in full analogy to linear chains. By analyzing the order parameter of the adsorption transition, we determine the critical point of adsorption of the dendrimers which is located close to the critical point of adsorption for star-polymers. While the order parameter for the adsorbed spacers displays a critical crossover scaling, the overall order parameter, which combines both critical and discontinuous transition effects, does not display simple scaling. The step-like transition from the weak into the strong adsorbed regime is confirmed by analyzing the shape-anisotropy of the dendrimers. We present a mean-field model based on the concept of spacer adsorption which predicts a discontinuous transition of dendrimers due to an excluded volume barrier. The latter results from an increased density of the dendrimer in the flatly adsorbed state which has to be overcome before this state is thermodynamically stable.

Two universality classes for random hyperbranched polymers

We grow AB2 random hyperbranched polymer structures in different ways and using different simulation methods. In particular we use a method of ad hoc construction of the connectivity matrix and the bond fluctuation model on a 3D lattice. We show that hyperbranched polymers split into two universality classes depending on the growth process. For a "slow growth" (SG) process where monomers are added sequentially to an existing molecule which strictly avoids cluster-cluster aggregation the resulting structures share all characteristic features with regular dendrimers. For a "quick growth" (QG) process which allows for cluster-cluster aggregation we obtain structures which can be identified as random fractals. Without excluded volume interactions the SG model displays a logarithmic growth of the radius of gyration with respect to the degree of polymerization while the QG model displays a power law behavior with an exponent of 1/4. By analyzing the spectral properties of the connectivity matrix we confirm the behavior of dendritic structures for the SG model and the corresponding fractal properties in the QG case. A mean field model is developed which explains the extension of the hyperbranched polymers in an athermal solvent for both cases. While the radius of gyration of the QG model shows a power-law behavior with the exponent value close to 4/5, the corresponding result for the SG model is a mixed logarithmic-power-law behavior. These different behaviors are confirmed by simulations using the bond fluctuation model. Our studies indicate that random sequential growth according to our SG model can be an alternative to the synthesis of perfect dendrimers.

Charged dendrimers in trivalent salt solutions under the action of DC electric fields

The structural properties and electrophoretic mobility of charged dendrimers in 3:1 electrolyte solutions subjected to direct current electric fields are studied using molecular dynamics simulations. The simulated dendrimer size is studied in zero fields and found to scale as R(g) ∼ N(0.29). The dendrimers exhibit shape distortions when the applied electric field is larger than some critical value, which scales with the number of dendrimer monomers as E(z,crit) ∼ N(0.39(6)). Families of curves, such as the curves of the square of radius of gyration, the asphericity, the degree of prolateness, and the electrophoretic mobility of dendrimers, are shown to collapse to single, master curves in electric fields through appropriate scaling. This reflects the fractal characteristics of these systems. The density profile of the surface monomers and salt cations reveals two pronounced combination effects between the polarization of dendrimer complexes and stripping-off of the condensed salt cations from the dendrimer surface.

The escape of a charged dendrimer from an oppositely charged planar surface

Many of the envisioned applications of dendrimers revolve around placing these molecules at and removing them from charged interfaces. Herein, we provide a prescription for the conditions needed to release a charged dendrimer from an oppositely charged flat substrate. Identifying an effective segment step length that reflects the intramolecular repulsions due to excluded volume and electrostatics, as well as the dendrimer's branching, provides the essential concept leading to an analytical prediction for the boundary between captured and free molecules. We find that this effective step length obeys trends similar to those predicted for linear chains, but is modified by the dendrimer's connectivity. Moreover, the boundary predicted for the capture of linear chains holds for dendrimers once this effective step length is employed. Monte Carlo computer simulations of coarse-grained model dendrimers escaping from charged surfaces validate these findings. The simulations consider generations 2 through 6 with a range of lengths between the branch points, as well as a range of solution ionic strengths and surface charge densities.

Theoretical and computational studies of dendrimers as delivery vectors

It is a great challenge for nanomedicine to develop novel dendrimers with maximum therapeutic potential and minimum side-effects for drug and gene delivery. As delivery vectors, dendrimers must overcome lots of barriers before delivering the bio-agents to the target in the cell. Extensive experimental investigations have been carried out to elucidate the physical and chemical properties of dendrimers and explore their behaviors when interacting with biomolecules, such as gene materials, proteins, and lipid membranes. As a supplement of the experimental techniques, it has been proved that computer simulations could facilitate the progress in understanding the delivery process of bioactive molecules. The structures of dendrimers in dilute solutions have been intensively investigated by monomer-resolved simulations, coarse-grained simulations, and atom-resolved simulations. Atomistic simulations have manifested that the hydrophobic interactions, hydrogen-bond interactions, and electrostatic attraction play critical roles in the formation of dendrimer-drug complexes. Multiscale simulations and statistical field theories have uncovered some physical mechanisms involved in the dendrimer-based gene delivery systems. This review will focus on the current status and perspective of theoretical and computational contributions in this field in recent years. (275 references).

Monte Carlo simulations of charged dendrimer-linear polyelectrolyte complexes and explicit counterions

We study complexes composed of one dendrimer of generation G = 4 (G4 dendrimer) with N(t) = 32 charged terminal groups and an oppositely charged linear polyelectrolyte accompanied by neutralizing counterions in an athermal solvent using Monte Carlo simulations based on the bond fluctuation model. In our study both the full Coulomb potential and the excluded volume interactions are taken into account explicitly with the reduced temperature τ and the chain length N(ch) as the main simulation parameters. Our calculations indicate that there exist three temperature ranges that determine the behavior of such complexes. At τ(complex) stable charged dendrimer-linear polyelectrolyte complexes are formed first, which are subsequently accompanied by selective counterion localization within the complex interior at τ(loc) ≤ τ(complex), and counterion condensation as temperature is further decreased below τ(cond) < τ(loc). In particular, we observe that condensation takes place exclusively on the excess charges in the complex and thus no condensation is observed at the compensation point (N(ch) = N(t)), irrespective of τ. For N(ch) ≠ N(t) the complex is overally charged. Furthermore, we discuss the size and structure of the dendrimer and the linear polyelectrolyte within the complex, as well as spatial distributions of monomers and counterions. Conformations of the chain in the bound state are analysed in terms of loops, trains, and tails.