Force Biased Molecular Dynamics Simulation Study of Effect of Dendrimer Generation on Interaction with DNA

Potential of Mean Force for DNA Wrapping Around a Cationic Nanoparticle

Sharp bending and wrapping of DNA around proteins and nanoparticles (NPs) has been of extensive research interest. Here, we present the potential of mean force (PMF) for wrapping a DNA double helix around a cationic NP using coarse-grained models of a double-stranded DNA and a cationic NP. Starting from a NP wrapped around by DNA, the PMF was calculated along the distance between the center of the NP and one end of the DNA molecule. A relationship between the distance and the extent of DNA wrapping is used to calculate the PMF as a function of DNA wrapping around a NP. In particular, the PMF was compared for two DNA sequences of (AT)₂₅/(AT)₂₅ and (AC)₂₅/(GT)₂₅, for which the persistence lengths are different by ∼10 nm. The simulation results provide solid evidence of the thermodynamic preference for complex formation of a cationic NP with more flexible DNA over the less flexible DNA. Furthermore, we estimated the elastic energy of DNA bending, which was in good order-of-magnitude agreement with the theoretical prediction of elastic rods. This work suggests that the variation of sequence-dependent DNA flexibility can be utilized in DNA nanotechnologies, in which the position and dynamics of NPs are regulated on large-scale DNA structures, or the structural transformation of DNA is triggered by the sequence-dependent binding of NPs.

Mechanistic Understanding from Molecular Dynamics in Pharmaceutical Research 2: Lipid Membrane in Drug Design

We review the use of molecular dynamics (MD) simulation as a drug design tool in the context of the role that the lipid membrane can play in drug action, i.e., the interaction between candidate drug molecules and lipid membranes. In the standard "lock and key" paradigm, only the interaction between the drug and a specific active site of a specific protein is considered; the environment in which the drug acts is, from a biophysical perspective, far more complex than this. The possible mechanisms though which a drug can be designed to tinker with physiological processes are significantly broader than merely fitting to a single active site of a single protein. In this paper, we focus on the role of the lipid membrane, arguably the most important element outside the proteins themselves, as a case study. We discuss work that has been carried out, using MD simulation, concerning the transfection of drugs through membranes that act as biological barriers in the path of the drugs, the behavior of drug molecules within membranes, how their collective behavior can affect the structure and properties of the membrane and, finally, the role lipid membranes, to which the vast majority of drug target proteins are associated, can play in mediating the interaction between drug and target protein. This review paper is the second in a two-part series covering MD simulation as a tool in pharmaceutical research; both are designed as pedagogical review papers aimed at both pharmaceutical scientists interested in exploring how the tool of MD simulation can be applied to their research and computational scientists interested in exploring the possibility of a pharmaceutical context for their research.

Phase equilibria and conformational behavior of dendrimers in porous media: Towards chromatographic analysis of dendrimers

HYPOTHESIS

The intricate entropy-enthalpy interplay of dendrimers confined in pores affects their conformation and retention in the porous stationary phase. This work aims at providing important insights into its impacts on partitioning and chromatographic separation in both size-exclusion chromatography (SEC) and interaction chromatography (IC) regimes.

SIMULATIONS

Using Monte Carlo (MC) simulations, we investigated the bulk-pore phase equilibria and the conformational behavior of flexible dendrimers differing in generation, in spacer length and in fraction of modified terminal groups interacting differently with pore walls than the majority building units.

FINDINGS

With increasing interaction strength, a distinct transition from a roughly spherical shape caused by simultaneous interactions with two walls to an ellipsoidal (or even disklike) conformation tenaciously adhering to only one wall was observed for moderately confined dendrimers. The strongly deformed dendrimers subjected to severe confinement gain high energy and the samples differing in the degree of modification become chromatographically discernable thanks to large energy differences. Consequently, our results suggest that the column fillings with fairly narrow pores which are ineffective in SEC, are highly efficient separation media for dendrimer studies by IC above the critical adsorption point (CAP). Overall, our simulations reveal useful information for advancing and optimizing experimental liquid chromatography studies of dendrimers.

Effect of DNA Flexibility on Complex Formation of a Cationic Nanoparticle with Double-Stranded DNA

We present extensive molecular dynamics simulations of a cationic nanoparticle and a double-stranded DNA molecule to discuss the effect of DNA flexibility on the complex formation of a cationic nanoparticle with double-stranded DNA. Martini coarse-grained models were employed to describe double-stranded DNA molecules with two different flexibilities and cationic nanoparticles with three different electric charges. As the electric charge of a cationic nanoparticle increases, the degree of DNA bending increases, eventually leading to the wrapping of DNA around the nanoparticle at high electric charges. However, a small increase in the persistence length of DNA by 10 nm requires a cationic nanoparticle with a markedly increased electric charge to bend and wrap DNA around. Thus, a more flexible DNA molecule bends and wraps around a cationic nanoparticle with an intermediate electric charge, whereas a less flexible DNA molecule binds to a nanoparticle with the same electric charge without notable bending. This work provides solid evidence that a small difference in DNA flexibility (as small as 10 nm in persistence length) has a substantial influence on the complex formation of DNA with proteins from a biological perspective and suggests that the variation of sequence-dependent DNA flexibility can be utilized in DNA nanotechnology as a new tool to manipulate the structure of DNA molecules mediated by nanoparticle binding.

Mechanistic Understanding From Molecular Dynamics Simulation in Pharmaceutical Research 1: Drug Delivery

In this review, we outline the growing role that molecular dynamics simulation is able to play as a design tool in drug delivery. We cover both the pharmaceutical and computational backgrounds, in a pedagogical fashion, as this review is designed to be equally accessible to pharmaceutical researchers interested in what this new computational tool is capable of and experts in molecular modeling who wish to pursue pharmaceutical applications as a context for their research. The field has become too broad for us to concisely describe all work that has been carried out; many comprehensive reviews on subtopics of this area are cited. We discuss the insight molecular dynamics modeling has provided in dissolution and solubility, however, the majority of the discussion is focused on nanomedicine: the development of nanoscale drug delivery vehicles. Here we focus on three areas where molecular dynamics modeling has had a particularly strong impact: (1) behavior in the bloodstream and protective polymer corona, (2) Drug loading and controlled release, and (3) Nanoparticle interaction with both model and biological membranes. We conclude with some thoughts on the role that molecular dynamics simulation can grow to play in the development of new drug delivery systems.

PAMAM dendrimer: a pH-controlled nanosponge

Using fully atomistic molecular dynamics simulation that are several hundred nanoseconds long, we demonstrate the pH-controlled sponge action of PAMAM dendrimer. We show how at varying pH levels, the PAMAM dendrimer acts as a wet sponge; at neutral or low pH levels, the dendrimer expands noticeably and the interior of the dendrimer opens up to host several hundreds to thousands of water molecules depending on the generation number. Increasing the pH (i.e., going from low pH to high pH) leads to the collapse of the dendrimer size, thereby expelling the inner water, which mimics the ‘sponge’ action. As the dendrimer size swells up at a neutral pH or low pH due to the electrostatic repulsion between the primary and tertiary amines that are protonated at this pH, there is dramatic increase in the available solvent accessible surface area (SASA), as well as solvent accessible volume (SAV).

Polymer-Nucleic Acid Interactions

Gene therapy is an important therapeutic strategy in the treatment of a wide range of genetic disorders. Polymers forming stable complexes with nucleic acids (NAs) are non-viral gene carriers. The self-assembly of polymers and nucleic acids is typically a complex process that involves many types of interaction at different scales. Electrostatic interaction, hydrophobic interaction, and hydrogen bonds are three important and prevalent interactions in the polymer/nucleic acid system. Electrostatic interactions and hydrogen bonds are the main driving forces for the condensation of nucleic acids, while hydrophobic interactions play a significant role in the cellular uptake and endosomal escape of polymer-nucleic acid complexes. To design high-efficiency polymer candidates for the DNA and siRNA delivery, it is necessary to have a detailed understanding of the interactions between them in solution. In this chapter, we survey the roles of the three important interactions between polymers and nucleic acids during the formation of polyplexes and summarize recent understandings of the linear polyelectrolyte-NA interactions and dendrimer-NA interactions. We also review recent progress optimizing the gene delivery system by tuning these interactions.

Binding of single stranded nucleic acids to cationic ligand functionalized gold nanoparticles

The interactions of nanoparticles (NPs) with single stranded nucleic acids (NAs) have important implications in gene delivery, and nanotechnological and biomedical applications. Here, the complexation of cationic ligand functionalized gold nanoparticles with single stranded deoxyribose nucleic acid (DNA) and ribonucleic acid (RNA) are examined using all atom molecular dynamics simulations. The results indicated that complexation depends mostly on charge of nanoparticle, and, to lesser extent, sequence and type of nucleic acid. For cationic nanoparticles, electrostatic interactions between charged ligands and the nucleic acid backbone dominate binding regardless of nanoparticle charge. Highly charged nanoparticles bind more tightly and cause compaction of the single-stranded NAs through disruption of intrastrand π-π stacking and hydrogen bonding. However, poly-purine strands (polyA-DNA, polyA-RNA) show less change in structure than poly-pyrimidine strands (polyT-DNA, polyU-RNA). Overall, the results show that control over ssNA structure may be achieved with cationic NPs with a charge of more than 30, but the extent of the structural changes depends on sequence.

Computational and experimental approaches for investigating nanoparticle-based drug delivery systems

Most therapeutic agents suffer from poor solubility, rapid clearance from the blood stream, a lack of targeting, and often poor translocation ability across cell membranes. Drug/gene delivery systems (DDSs) are capable of overcoming some of these barriers to enhance delivery of drugs to their right place of action, e.g. inside cancer cells. In this review, we focus on nanoparticles as DDSs. Complementary experimental and computational studies have enhanced our understanding of the mechanism of action of nanocarriers and their underlying interactions with drugs, biomembranes and other biological molecules. We review key biophysical aspects of DDSs and discuss how computer modeling can assist in rational design of DDSs with improved and optimized properties. We summarize commonly used experimental techniques for the study of DDSs. Then we review computational studies for several major categories of nanocarriers, including dendrimers and dendrons, polymer-, peptide-, nucleic acid-, lipid-, and carbon-based DDSs, and gold nanoparticles. This article is part of a Special Issue entitled: Membrane Proteins edited by J.C. Gumbart and Sergei Noskov.

On the Possibility of Facilitated Diffusion of Dendrimers Along DNA

We investigate the electrostatics, energetics, and dynamics of dendrimer-DNA interactions that mimic protein-DNA complexes as a means to design facilitated mechanisms by which dendrimers can slide and search DNA for targets. By using all-atom molecular dynamics simulations, we calculated the free energy profiles of dendrimer-binding around the DNA via umbrella sampling. We also calculated electrostatic interaction maps in comparison to proteins, as well as the dynamical changes induced by DNA-dendrimer interactions via NMR-measurable order parameters. Our results show that for dendrimers to go around DNA, there is a free-energy barrier of 8.5 kcal/mol from the DNA major groove to DNA minor groove, with a minimum in the major groove. This barrier height makes it unlikely for an all-amine dendrimer to slide along DNA longitudinally, but following a helical path may be possible along the major groove. Comparison of the nonbonded interaction energy and the interaction free-energy profiles reveal a considerable entropic cost as the dendrimer binds to DNA. This is also supported by the mobility patterns obtained from NMR-measurable order parameter values, which show a decreased mobility of the dendrimer N-H bond vectors in the DNA-binding mode.

Association of the anti-tuberculosis drug rifampicin with a PAMAM dendrimer

The association of the anti-tuberculosis drug rifampicin (RIF) with a 4th-generation poly(amidoamine) (G4-PAMAM) dendrimer was investigated by means of molecular dynamics simulations. The RIF load capacity was estimated to be around 20 RIF per G4-PAMAM at neutral pH. The complex formed by 20 RIF molecules and the dendrimer (RIF20-PAMAM) was subjected to 100 ns molecular dynamics (MD) simulations at two different pH conditions (neutral and acidic). The complex was found to be significantly more stable in the simulation at neutral pH compared to the simulation at low pH in which the RIF molecules were rapidly and almost simultaneously expelled to the solvent bulk. The high stability of the RIF-PAMAM complex under physiological pH and the rapid release of RIF molecules under acidic medium provide an interesting switch for drug targeting since the Mycobacterium resides within acidic domains of the macrophage. Altogether, these results suggest that, at least in terms of stability and pH-dependent release, PAMAM-like dendrimers may be considered suitable drug delivery systems for RIF and derivatives.

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.

Transport of nucleosides in the vcCNT facilitated by sodium gradients from molecular dynamics simulations

Nucleosides are required for DNA and RNA synthesis, and the nucleoside adenosine has a function in a variety of signaling processes. Nucleosides require a specialized class of integral membrane proteins, known as nucleoside transporters (NTs), for specific transport across cell membranes. NTs are also important determinants for the transport of nucleoside-derived drugs across cell membranes. Recently, the crystal structure of the vcCNT (Vibrio cholerae Concentrative Nucleoside Transporter) was reported. Here we perform molecular dynamics (MD) simulations for the vcCNT structure in the presence of various sodium gradients, since CNTs are sodium-coupled transporters. The results highlight the important role of sodium bound to the vcCNT in the transport of uridine. Our MD simulations show that, without NaCl, uridine remains stable in the binding pocket of the vcCNT. In the presence of 20 mM NaCl, uridine moves from the binding pocket and approaches the entrance of the intracellular side. In the presence of 100 mM NaCl, uridine passes through most part of the entrance and approaches the intracellular side. The polar/charged amino acids in the binding pocket are important in the transport process. They first "fix" the ribose and allow the uracil base of uridine to approach the entrance of the intracellular side, and then "release" the ribose to allow uridine to move freely into the intracellular side coupled with the movement of sodium ions and HP1b. Finally, we propose a detailed mechanism of the nucleoside transport from the binding pocket to the intracellular side of the vcCNT.