A Multiscale Scheme for the Simulation of Conformational and Solution Properties of Different Dendrimer Molecules

Temperature-Transferable Coarse-Grained Models for Volumetric Properties of Poly(lactic Acid)

A new coarse-grained (CG) model, for which each monomer is mapped as one bead at its center of mass, was developed for simulating the volumetric properties of the polylactide (PLA) bulk. The three bonded CG potentials are first parametrized against the strain energies of the dimer, trimer, and tetramer, and the nonbonded CG potentials are then optimized to match the melt densities of the decamer. With the derived CG potentials, molecular dynamics (MD) simulations are found to reproduce thermal expansion and glass transition. By rescaling the dihedral and nonbonded potentials with temperature-independent factors, the glass transition temperature (Tg) is also satisfactorily restored with little modifications on the volumetric expansive coefficients at both rubbery and glassy states. Therefore, the finally optimized CG potentials exhibit excellent temperature transferability, as rationalized by the Simha-Boyer relation. Furthermore, it is confirmed that the dihedral torsions and nonbonded interactions play key roles in glass transition. Also, the simulated bulk moduli and conformational properties in a wide temperature range compare well with the referenced data. The proposed multiscale scheme has great potential in simulating thermo-mechanical properties of PLA and other polymers.

An Enhanced Scheme for Multiscale Modeling of Thermomechanical Properties of Polymer Bulks

While multiscale modeling significantly enhances the capability of molecular simulations of polymer systems, it is well realized that the systematically derived coarse-grained (CG) models generally underestimate the thermomechanical properties. In this work, a charge-based mapping scheme has been adopted to include explicit electrostatic interactions and benchmarked against two typical polymers, atactic poly(methyl methacrylate) (PMMA) and polystyrene (PS). The CG potentials are parameterized against the oligomer bulks of nine monomers per chain to match the essential structural features and the two basic pressure-volume-temperature (PVT) properties, which are obtained from the all-atomistic (AA) molecular dynamics (MD) simulations at a single elevated temperature. The so-parameterized CG potentials are extended with the MD method to simulate the two polymer bulks of one hundred monomers per chain over a wide temperature range. Without any scaling, all the simulated results, including mass densities and bulk moduli at room temperature, thermal expansion coefficients at rubbery and glassy states, and glass transition temperatures (T_g), compare well with the corresponding experimental data. The proposed scheme not only contributes to realistically simulating various thermomechanical properties of both apolar and polar polymers but also allows for directly simulating their electrical properties.

Exploring the Structural Properties of Positively Charged Peptide Dendrimers

We report a combined experimental and computational approach to study the structural behavior of positively charged peptide dendrimers. Third-generation dendrimers containing combinations of positive/neutral amino acid residues in the different dendrimer generations were synthesized and their overall size evaluated using diffusion NMR. Molecular dynamics simulations were performed to obtain a comprehensive description of the molecular-level phenomena substantiating the structural differences observed. Comparison of the results presented with previous findings reveals a striking charge-dependent tendency in these systems, where the simple number and placement of charged amino acids in the sequence allows an extensive control over the exhibited structural features. Indeed, we observe that peptide dendrimers bearing progressively higher amounts of charged residues are characterized by an increasing structural plasticity, with a myriad of conformational states equally accessible to them. On the other hand, dendrimers containing only small amounts of charged residues evidence, to some extent, a characteristic structural rigidity.

HBP Builder: A Tool to Generate Hyperbranched Polymers and Hyperbranched Multi-Arm Copolymers for Coarse-grained and Fully Atomistic Molecular Simulations

Computer simulation has been becoming a versatile tool that can investigate detailed information from the microscopic scale to the mesoscopic scale. However, the crucial first step of molecular simulation is model building, particularly for hyperbranched polymers (HBPs) and hyperbranched multi-arm copolymers (HBMCs) with complex and various topological structures. Unlike well-defined polymers, not only the molar weight of HBPs/HBMCs with polydispersity, but the HBPs/HBMCs with the same degree of polymerization (DP) and degree of branching (DB) also have many possible topological structures, thus making difficulties for user to build model in molecular simulation. In order to build a bridge between model building and molecular simulation of HBPs and HBMCs, we developed HBP Builder, a C language open source HBPs/HBMCs building toolkit. HBP Builder implements an automated protocol to build various coarse-grained and fully atomistic structures of HBPs/HBMCs according to user's specific requirements. Meanwhile, coarse-grained and fully atomistic output structures can be directly employed in popular simulation packages, including HOOMD, Tinker and Gromacs. Moreover, HBP Builder has an easy-to-use graphical user interface and the modular architecture, making it easy to extend and reuse it as a part of other program.

Ligand-Receptor Interaction-Mediated Transmembrane Transport of Dendrimer-like Soft Nanoparticles: Mechanisms and Complicated Diffusive Dynamics

Nearly all nanomedical applications of dendrimer-like soft nanoparticles rely on the functionality of attached ligands. Understanding how the ligands interact with the receptors in cell membrane and its further effect on the cellular uptake of dendrimer-like soft nanoparticles is thereby a key issue for their better application in nanomedicine. However, the essential mechanism and detailed kinetics for the ligand-receptor interaction-mediated transmembrane transport of such unconventional nanoparticles remain poorly elucidated. Here, using coarse-grained simulations, we present the very first study of molecular mechanism and kinetics behaviors for the transmembrane transport of dendrimer-like soft nanoparticles conjugated with ligands. A phase diagram of interaction states is constructed through examining ligand densities and membrane tensions that allows us to identify novel endocytosis mechanisms featured by the direct wrapping and the penetration-extraction vesiculation. The results provide an in-depth insight into the diffusivity of receptors and dendrimer in the membrane plane and demonstrate how the ligand density influences receptor diffusion and uptake kinetics. It is interesting to find that the ligand-conjugated dendrimers present superdiffusive behaviors on a membrane, which is revealed to be driven by the random fluctuation dynamics of the membrane. The findings facilitate our understanding of some recent experimental observations and could establish fundamental principles for the future development of such important nanomaterials for widespread nanomedical applications.

Molecular dynamics simulation studies of hyperbranched polyglycerols and their encapsulation behaviors of small drug molecules

Computer simulation could disclose more details about the conformations of HPGs and their encapsulation behaviors of guest molecules.

Unique dynamical approach of fully wrapping dendrimer-like soft nanoparticles by lipid bilayer membrane

Wrapping dendrimer-like soft nanoparticles by cell membrane is an essential event in their endocytosis in drug and gene delivery, but this process remains poorly elucidated. Using computer simulations and theoretical analysis, we report the detailed dynamics of the process in which a lipid bilayer membrane fully wraps a dendrimer-like soft nanoparticle. By constructing a phase diagram, we firstly demonstrate that there exist three states in the interaction between a dendrimer and a lipid bilayer membrane, i.e., penetration, penetration and partial wrapping, and full wrapping states. The wrapping process of dendrimer-like nanoparticles is found to take a unique approach where the penetration of the dendrimer into the membrane plays a significant role. The analysis of various energies within the system provides a theoretical justification to the state transition observed from simulations. The findings also support recent experimental results and provide a theoretical explanation for them. We expect that these findings are of immediate interest to the study of the cellular uptake of dendrimer-like soft nanoparticles and can prompt the further application of this class of nanoparticles in nanomedicine.

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).

Brownian dynamics simulations on CPU and GPU with BD_BOX

There has been growing interest in simulating biological processes under in vivo conditions due to recent advances in experimental techniques dedicated to study single particle behavior in crowded environments. We have developed a software package, BD_BOX, for multiscale Brownian dynamics simulations. BD_BOX can simulate either single molecules or multicomponent systems of diverse, interacting molecular species using flexible, coarse-grained bead models. BD_BOX is written in C and employs modern computer architectures and technologies; these include MPI for distributed-memory architectures, OpenMP for shared-memory platforms, NVIDIA CUDA framework for GPGPU, and SSE vectorization for CPU.