Brownian Dynamics Simulations of Dendrimers under Shear Flow

Drug-dendrimer complexes and conjugates: Detailed furtherance through theory and experiments

Dendritic nanovectors-based drug delivery has gained significant attention in the past couple of decades. Dendrimers play a crucial role in deciding the solubility of sparingly soluble drug molecules and help in improving pharmacokinetics. A few important steps in drug delivery through dendrimers, such as drug encapsulation, formulation, and target-specific delivery, play an important role in deciding the fate of a drug molecule. It is also of prime importance to understand the interactions between a drug molecule and dendrimers at atomistic levels to decode the mechanism of action of drug-dendrimer complexes and their reliability in terms of drug delivery. Colossal progress in current experimental and computational approaches in the field has resulted in a vast amount of data that needs to be curated to be further implemented efficiently. Improved computational power has led to greater accuracy and prompt predictions of properties of drug-dendrimer complexes and their mechanism of action. The current review encapsulates the pioneering work in the field, experimental achievements in terms of drug delivery, and newer computational techniques employed in the advancement of the field.

QSPR modelling for intrinsic viscosity in polymer-solvent combinations based on density functional theory

Linear and nonlinear quantitative structure-property relationship (QSPR) models were developed based on a dataset with 65 polymer-solvent combinations. Seven quantum chemical descriptors, dipole moment, hardness, chemical potential, electrophilicity index, total energy, HOMO and LUMO orbital energies, were calculated with density functional theory at the B3LYP/6-31 G(d) level for polymers and solvents. Considering the strong correlation between intrinsic viscosity and weight, size, shape as well as topological structure of polymers and solvents, topological descriptors were also applied in this work. Meanwhile, the most appropriate polymer structure representation was investigated by considering 1-5 monomeric repeating units. The molecular descriptors were first screened by using the genetic algorithms-multiple linear regression (GA-MLR), with coefficient of determinations (r2) of 0.78 and 0.83 for the training set and the prediction set, respectively. The support vector machine model (SVM) model based on the selected descriptors subset showed a r2 value of 0.95 for the training set and 0.93 for the prediction set. All statistical results suggest that the established QSPR models have good predictability. Furthermore, a new test set obtained from the literature was used for further validation. The r2 values were 0.81 for the MLR model and 0.90 for the SVM model.

Static Structure Factor and Viscoelastic Properties of Dendrimer Grafted Nanoparticles in Solution

The theory for the dynamics of multiscale branched polymeric structures is applied to understand the dendrimer-grafted nanoparticles in a dilute solution. The multiscale nature of dendrimer-grafted nanoparticles arises due to larger beads for the nanoparticles and the smaller beads for the polymeric structure connected through the harmonic springs. The multiscale generalized Gaussian structure approach allows us to study several viscoelastic properties: (i) storage and loss moduli and (ii) intrinsic viscosity. The influence of nanoparticles in the dendrimer structure is reflected in low and intermediate frequency regimes of the viscoelastic relaxation moduli. The increase in the size and the number fraction of nanoparticle shows an anomalous enhancement in the relaxation moduli. The increase in number fraction of nanoparticle in dendrimer-grafted nanoparticles decreases the transition frequency between solid- and liquid-like viscoelastic region. The intrinsic viscosity of dendrimer-grafted nanoparticles increases with increasing the size of nanoparticle. The inclusion of hydrodynamic interactions facilitates the dynamics of dendrimer-grafted nanoparticles. The Kratky plot of the static structure factor of all conformation of dendrimer-grafted nanoparticles is also analyzed as a function of number fraction and the size of the nanoparticles. At low wavenumbers, all conformations of dendrimer-grafted nanoparticles show a universal behavior. The compactness of dendrimer-grafted nanoparticles increases with the increase in number fraction and the size of the nanoparticles.

Hybrid Polycarbosilane-Siloxane Dendrimers: Synthesis and Properties

A series of carbosilane dendrimers of the 4th, 6th, and 7th generations with a terminal trimethylsilylsiloxane layer was synthesized. Theoretical models of these dendrimers were developed, and equilibrium dendrimer conformations obtained via molecular dynamics simulations were in a good agreement with experimental small-angle X-ray scattering (SAXS) data demonstrating molecule monodispersity and an almost spherical shape. It was confirmed that the glass transition temperature is independent of the dendrimer generation, but is greatly affected by the chemical nature of the dendrimer terminal groups. A sharp increase in the zero-shear viscosity of dendrimer melts was found between the 5th and the 7th dendrimer generations, which was qualitatively identical to that previously reported for polycarbosilane dendrimers with butyl terminal groups. The viscoelastic properties of high-generation dendrimers seem to follow some general trends with an increase in the generation number, which are determined by the regular branching structure of dendrimers.

Dendrimer Dynamics: A Review of Analytical Theories and Molecular Simulation Methods

The theoretical study of dendrimers is reviewed, considering both analytical approaches and molecular simulation methods. We discuss the effect of molecular symmetry on the degeneracy of the relaxation times, and then the calculation of observable quantities, in particular the intrinsic viscosity, and then the viscoelastic complex modulus and the dynamic structure factor, in comparison with the available experimental data. In particular, the maximum intrinsic viscosity with increasing molar mass is analyzed in some detail. The approximations and/or assumptions of the adopted methods are also described in connection with analogous results for polymer of a different topology, in particular linear and star polymers.

Orientational Relaxation of Poly(propylene imine) Dendrimers at Different pH

The dilute solution dynamics of poly(propylene imine) (PPI) dendrimers is investigated at three different solution pH through molecular dynamics (MD) simulations. The dynamics of PPI dendrimers is characterized by both global and local relaxations that occur at different time and length scales. While the global dynamics may be described in terms of rotational diffusion, the local motion may be characterized through orientational relaxation dynamics measured in terms of the time autocorrelation function (ACF), second-order orientational ACF, and the spin-lattice relaxation rate. The global motion of dendrimers decreases with an increase in the size from high pH to low pH with increasing generations of growth. The results reveal that the segments at low pH relax faster than those at high pH, and the local mobility of the segments near the periphery is higher than the core segments. This observation is also evident from the spectral density and spin-lattice relaxation rate. High values of the spectral density at higher frequencies imply higher segmental mobility of the dendrimer at low pH relative to that at high pH. A shift in the maximum of the spin-lattice relaxation rate toward lower frequencies with decreasing generations indicates the dependence of local mobility on the topological distance of the segment from the periphery at all pH conditions.

Effects of generation number, spacer length and temperature on the structure and intramolecular dynamics of siloxane dendrimer melts: molecular dynamics simulations

The structure and properties of polysiloxane dendrimer melts are studied by extensive atomistic molecular dynamics simulations. Two homologous series differing in the spacer length are considered. In the first series the dendrimer spacers are the shortest ones, comprising only one oxygen atom, while in the second series the spacers consist of two oxygen atoms with the silicon atom in between. Melts of the dendrimers from the 3rd up to the 6th generation number are modelled in a wide temperature range from 273 to 600 K. A comparative study of the macroscopic melt characteristics such as the melt density and thermal expansion coefficients is performed for the two series. Analysis of the dendrimer structure in melts and in the isolated state shows that intermolecular interactions and interpenetration of dendrimer molecules in melts hardly affect the dendrimer interior organization. However, the presence of neighboring molecules significantly slows down their intramolecular dynamics in melts in comparison with that of isolated dendrimers. An increasing generation number causes an increase of the radius of the dendrimer interior region unavailable for neighboring molecules, which starts to exceed the length of the peripheral interpenetration layer for high-generation dendrimers; this fact could lead to different mechanisms of melt dynamics for lower and higher generation dendrimers.

Effect of pH on Size and Internal Structure of Poly(propylene imine) Dendrimers: A Molecular Dynamics Simulation Study

The behavior of poly(propylene imine) dendrimers at three different solution pH is investigated through molecular dynamics (MD) simulations in explicit solvent. MD simulations provide an insight into the conformational properties of dendrimers via the evaluation of their size, shape, radial density distribution, static structure factor, and scattering intensity. The size of the dendrimer increases from high solution pH to low pH. The internal structure of the dendrimer is quantified in terms of the radial atomic density profile and the terminal amine group density distribution. While the radial atomic density distribution shifts away from the core of the dendrimer with decreasing pH, a significant degree of back-folding of the outer generations is observed at high pH for higher generations of growth. Results from the structure factor and scattering intensity indicate two types of conformational transitions: (i) as a function of the solution pH, where the dendrimer evolves from an expanded structure at low pH to a highly compact one at high pH (except for higher generations), and (ii) with increasing generations, where the open structure of the dendrimer at lower generations transforms to a compact structure at higher generations at both high and low pH, characterized by a change in the fractal dimension.

A Comparative Study of Intramolecular Mobility of Single Siloxane and Carbosilane Dendrimers via Molecular Dynamics Simulations

A comparative analysis of intramolecular dynamics of four types of isolated dendrimers from the fourth to the seventh generations belonging to the siloxane and carbosilane families, differing in spacer length, core functionality, and the type of chemical bonds, has been performed via atomic molecular dynamics simulations. The average radial and angular positions of all Si branching atoms of various topological layers within the dendrimer interior, as well as their variations, have been calculated, and the distributions of the relaxation times of their radial and angular motions have been found. It has been shown that the dendrons of all the dendrimers elongate from the center and decrease in a solid angle with an increasing generation number. The characteristic relaxation times of both angular and radial motions of Si atoms are of the order of a few nanoseconds, and they increase with an increasing generation number and decrease with temperature, with the angular relaxation times being larger than the radial ones. The relaxation times in the carbosilanes are larger than those in the siloxanes. The rotational angle dynamics of the carbosilane dendrimers show that the chain bending is mainly realized via trans-gauche transitions in the Si branching bonds.

Molecular dynamics simulations of single siloxane dendrimers: Molecular structure and intramolecular mobility of terminal groups

Molecular dynamics simulations of two types of isolated siloxane dendrimers of various generations (from the 2nd to the 8th) have been performed for temperatures ranging from 150 K to 600 K. The first type of dendrimer molecules has short spacers consisting of a single oxygen atom. In the dendrimers of the second type, spacers are longer and comprised of two oxygen atoms separated by a single silicon atom. A comparative analysis of molecular macroscopic parameters such as the gyration radius and the shape factor as well as atom distributions within dendrimer interior has been performed for varying generation number, temperature, and spacer length. It has been found that the short-spacer dendrimers of the 7th and 8th generations have a stressed central part with elongated bonds and deformed valence angles. Investigation of the time evolution of radial displacements of the terminal Si atoms has shown that a fraction of the Si groups have a reduced mobility. Therefore, rather long time trajectories (of the order of tens of nanoseconds) are required to study dendrimer intramolecular dynamics.

Hyperbranched polymer stars with Gaussian chain statistics revisited

Conformational properties of regular dendrimers and more general hyperbranched polymer stars with Gaussian statistics for the spacer chains between branching points are revisited numerically. We investigate the scaling for asymptotically long chains especially for fractal dimensions df = 3 (marginally compact) and df = 2.5 (diffusion limited aggregation). Power-law stars obtained by imposing the number of additional arms per generation are compared to truly self-similar stars. We discuss effects of weak excluded-volume interactions and sketch the regime where the Gaussian approximation should hold in dense solutions and melts for sufficiently large spacer chains.

Mechanochemistry of Topological Complex Polymer Systems

Although existing since the concept of macromolecules, polymer mechanochemistry is a burgeoning field which attracts great scientific interest in its ability to bias conventional reaction pathways and its potential to fabricate mechanoresponsive materials. We review here the effect of topology on the mechanical degradation of polymer chains and the activation of mechanophores in polymer backbones. The chapter focuses on both experimental and theoretical work carried out in the past 70 years. After a general introduction (Sect. 1), where the concept, the history, and the application of polymer mechanochemistry are briefly described, flow fields to study polymer mechanochemistry are discussed (Sect. 2), results of mechanochemistry study are presented for linear polymers (Sect. 3), cyclic polymers (Sect. 4), graft polymers (Sect. 5), star-shaped polymers (Sect. 6), hyperbranched polymers and dendrimers (Sect. 7), and systems with dynamic topology (Sect. 8). Here we focus on (1) experimental results involving the topological effect on the coil-to-stretch transition and the fracture of the polymer chains, (2) the underlying mechanisms and the key factor that determines the mechanical stability of the macromolecules, (3) theoretical models that relate to the experimental observations, and (4) rational design of mechanophores in complex topology to achieve multiple activations according to the existing results observed in chain degradation.

Computational design principles for the discovery of bioactive dendrimers: [s]-triazines and other examples

INTRODUCTION

Chemistry yields dendrimers of many classes and compositions. Translating this synthetic success to bioactivity is significantly aided by the use of computational modeling and our knowledge of the three-dimensional shapes of these macromolecules.

AREAS COVERED

This review discusses the lessons learned during the investigations of [s]-triazine dendrimers. Specifically, the article focuses on the evolving role that computational models have taken in the exploration of these macromolecules. These lessons, furthermore, can be generalized across many dendrimer classes.

EXPERT OPINION

Computational models and the resulting structural data from molecular dynamics simulations provide insights into: shape, solvent penetration, shielding of biolabile linkers, and the density of hydrophobic patches. These models have evolved from artistic representations, through bases for rationalization, to hypothesis-generating tools that drive synthesis. With further advances expected in both software and hardware the answer to the question, 'What does a specific dendrimer look like in solution?' is becoming increasingly clear. Moreover, the authors believe that answer to this question lies at the heart of the design of bioactive dendrimers.

Structure of polyamidoamide dendrimers up to limiting generations: a mesoscale description

The polyamidoamide (PAMAM) class of dendrimers was one of the first dendrimers synthesized by Tomalia and co-workers at Dow. Since its discovery the PAMAMs have stimulated many discussions on the structure and dynamics of such hyperbranched polymers. Many questions remain open because the huge conformation disorder combined with very similar local symmetries have made it difficult to characterize experimentally at the atomistic level the structure and dynamics of PAMAM dendrimers. The higher generation dendrimers have also been difficult to characterize computationally because of the large size (294,852 atoms for generation 11) and the huge number of conformations. To help provide a practical means of atomistic computational studies, we have developed an atomistically informed coarse-grained description for the PAMAM dendrimer. We find that a two-bead per monomer representation retains the accuracy of atomistic simulations for predicting size and conformational complexity, while reducing the degrees of freedom by tenfold. This mesoscale description has allowed us to study the structural properties of PAMAM dendrimer up to generation 11 for time scale of up to several nanoseconds. The gross properties such as the radius of gyration compare very well with those from full atomistic simulation and with available small angle x-ray experiment and small angle neutron scattering data. The radial monomer density shows very similar behavior with those obtained from the fully atomistic simulation. Our approach to deriving the coarse-grain model is general and straightforward to apply to other classes of dendrimers.

Energetic and conformational aspects of dendrimer overcharging by linear polyelectrolytes

Extensive Brownian dynamics simulations of conformational changes accompanying the overcharging of a dendrimer by an oppositely charged long linear polyelectrolyte (LPE) have been carried out. The simulated results have been compared with the predictions of the Nguen and Shklovskii correlation theory [Physica A 293, 324 (2001)] for impenetrable charged spherical macroion. Dendrimer overcharging is caused by the spatial correlations between the "excess" of the LPE charges adsorbed onto its surface. The simulated LPE-length dependence of the corresponding "correlation" energy is in agreement with the theoretical predictions. Maximum of the LPE adsorption occurs at some critical LPE length N{ch};{c} , and the first order phase transition from completely coiled conformation to the conformation with released tails takes place. The phase transition is accompanied by the drastic increase in the relative fluctuations of the polyelectrolyte size. Upon increasing the linear-chain length above N{ch};{c} , the one-long-tail conformation becomes energetically preferable; the exchange time between the long-tail conformation and the short-tail conformation is very large.

Structure of poly(propyl ether imine) dendrimer from fully atomistic molecular dynamics simulation and by small angle x-ray scattering

We study the structure of carboxylic acid terminated neutral poly(propyl ether imine) (PETIM) dendrimer from generations 1-6 (G1-G6) in a good solvent (water) by fully atomistic molecular dynamics (MD) simulations. We determine as a function of generation the structural properties such as radius of gyration, shape tensor, asphericity, fractal dimension, monomer density distribution, and end-group distribution functions. The sizes obtained from the MD simulations have been validated by small angle x-ray scattering experiment on dendrimer of generations 2-4 (G2-G4). A good agreement between the experimental and theoretical value of radius of gyration has been observed. We find a linear increase in radius of gyration with the generation. In contrast, Rg scales as approximately Nx with the number of monomers. We find two distinct exponents depending on the generations, x=0.47 for G1-G3 and x=0.28 for G3-G6, which reveal their nonspace filling nature. In comparison with the amine terminated poly(amidoamine) (PAMAM) dendrimer, we find that Rg of Gth generation PETIM dendrimer is nearly equal to that of (G+1)th generation of PAMAM dendrimer as observed by Maiti et al. [Macromolecules 38, 979 (2005)]. We find substantial back folding of the outer subgenerations into the interior of the dendrimer. Due to their highly flexible nature of the repeating branch units, the shape of the PETIM dendrimer deviates significantly from the spherical shape and the molecules become more and more spherical as the generation increases. The interior of the dendrimer is quite open with internal cavities available for accommodating guest molecules, suggesting the use of PETIM dendrimer for guest-host applications. We also give a quantitative measure of the number of water molecules present inside the dendrimer.

Soft matter with soft particles

In this review we present a summary of recent progress achieved in examining the equilibrium and dynamical properties of concentrated solutions of two novel kinds of soft matter systems: and starburst molecules known as . The two systems share a host of interesting properties. The both consist of highly branched polymers, they allow for tuning of their properties through modification of the macromolecular architecture and they are both representatives of a quite novel class of colloidal particles, termed . On the other hand there are also important differences, reflecting the fundamental difference in their architecture. It will be shown that a combination of scattering techniques and rheology with computer simulations and analytical methods from the realm of theoretical physics can shed light on the unusual properties of such systems. In this fashion, new ways appear for the manipulation of soft matter systems under external influences and promising perspectives for the fabrication of new materials are opened up and the versatility in manipulating soft matter is underlined.

Monte Carlo calculations for the intrinsic viscosity of several dendrimer molecules

We have performed Monte Carlo simulations to reproduce the intrinsic viscosity corresponding to different generation of several types of dendrite molecules: polyamidoamine dendrimers with an ethylendiamine core, polypropylene-imine with a diaminobutane core, and monodendrons and tridendrons of polybenzylether. With this end, we have employed coarse-grained idealizations of the molecules constituted by only two beads in each repeat unit (one in a branching or end unit and one intermediate along the repeat unit) and a simple hard-sphere potential between non-neighboring beads. Our goal is to investigate if this simple model is able to provide a reasonable description of some differences between these systems that have been observed experimentally, in particular, the location of the maximum in the intrinsic viscosity as a function of the generation number. Experimental radii of gyration in a given solvent are reproduced by a fit of the hard-sphere potential diameter. Subsequently, intrinsic viscosities are calculated by the variational approach of Fixman, which yields an accurate lower-bound value with an additional hydrodynamic interaction parameter (the friction radius of the beads). The results show a pronounced variation of the maximum location with the value of the friction radius and the structural details that cannot be mimicked with simpler models. The initial conformations for the Monte Carlo procedure are taken from atomistic configurations thermalized by means of a molecular dynamics.