Attractive hydration forces in DNA-dendrimer interactions on the nanometer scale

Single-molecule measurements of double-stranded DNA condensation

Electrostatically driven double-stranded DNA (dsDNA) condensation is critical in regulating many biological processes, including bacteriophage and virus replication and the packaging of chromosomal DNA in sperm heads. Here, we review single-molecule measurements of dsDNA condensed by cationic proteins, polypeptides, and small multivalent cations. Optical tweezers (OT) measurements of dsDNA collapsed by cationic condensing agents reveal a critical condensing force unique to each condensing agent that is tunable with condensing agent concentration and ionic strength. DNA globules visualized via atomic force microscopy, transmission electron microscopy, and cryoelectron microscopy reveal condensed dsDNA adopting several conformations including highly ordered toroids with a measurable central hole and, more recently, the maximally dense, yarn-ball-like structures observed with dsDNA condensed by the HIV-1 nucleocapsid protein. The combination of these approaches provides multifaceted insight into the shape and size of electrostatically condensed dsDNA globules and the kinetics of their formation and dissolution. We also review the physics of dsDNA condensation, including recent studies that show dsDNA globule size is tunable with ionic strength. Overall, this review provides important insights into understanding dsDNA condensate-regulated biological processes, as well as potential uses for gene delivery.

Grand canonical Monte Carlo and deep learning assisted enhanced sampling to characterize the distribution of Mg2+ and influence of the Drude polarizable force field on the stability of folded states of the twister ribozyme

Molecular dynamics simulations are crucial for understanding the structural and dynamical behavior of biomolecular systems, including the impact of their environment. However, there is a gap between the time scale of these simulations and that of real-world experiments. To address this problem, various enhanced simulation methods have been developed. In addition, there has been a significant advancement of the force fields used for simulations associated with the explicit treatment of electronic polarizability. In this study, we apply oscillating chemical potential grand canonical Monte Carlo and machine learning methods to determine reaction coordinates combined with metadynamics simulations to explore the role of Mg2+ distribution and electronic polarizability in the context of the classical Drude oscillator polarizable force field on the stability of the twister ribozyme. The introduction of electronic polarizability along with the details of the distribution of Mg2+ significantly stabilizes the simulations with respect to sampling the crystallographic conformation. The introduction of electronic polarizability leads to increased stability over that obtained with the additive CHARMM36 FF reported in a previous study, allowing for a distribution of a wider range of ions to stabilize twister. Specific interactions contributing to stabilization are identified, including both those observed in the crystal structures and additional experimentally unobserved interactions. Interactions of Mg2+ with the bases are indicated to make important contributions to stabilization. Notably, the presence of specific interactions between the Mg2+ ions and bases or the non-bridging phosphate oxygens (NBPOs) leads to enhanced dipole moments of all three moieties. Mg2+-NBPO interactions led to enhanced dipoles of the phosphates but, interestingly, not in all the participating ions. The present results further indicate the importance of electronic polarizability in stabilizing RNA in molecular simulations and the complicated nature of the relationship of Mg2+-RNA interactions with the polarization response of the bases and phosphates.

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.

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.

Optimization and computational studies evaluating molecular dynamics of EDA cored polymeric dendrimer

In this work we report the results acquired from molecular dynamics simulations as well as the optimization of different generations of polyamidoamine dendrimer. The analysis data revealed synthesized dendrimer as a suitable nanostructured candidate suitable for neutral as well as charged molecule delivery due to the presence of both electrostatic potential and van der Waals forces. The methyl ester terminating groups of half-generation dendrimers with characteristic IR peaks for carbonyl at 1670.41 cm⁻¹ tends to shift to 1514.17 cm⁻¹ on conversion to amide group of full-generation dendrimer. The study includes the usage of detailed analysis, demonstrating how molecular dynamics affect the dendrimer complexation. The present investigations provide an unprecedented insight into the computational and experimental system that may be of general significance for the clinical application of dendrimers.

New Insights into the Interactions of a DNA Oligonucleotide with mPEGylated-PAMAM by Circular Dichroism and Solution NMR

Dendrimers are well-defined, highly branched, synthetic three-dimensional molecules with a large number of reactive end groups. PAMAM dendrimers form stable complexes with DNA chemistries and constitute an important class of nonviral, cationic vectors in gene delivery. The aim of this study is to examine the interactions of a 12 bp DNA oligonucletide with PAMAM-G2 and mPEG- b-PAMAM-G3 having eight surface amine groups under physiological conditions, using constant DNA concentration but varying dendrimer concentration. 1D ³¹P NMR, 2D NOESY, and CD spectroscopic methods were employed to investigate the interactions between the dendrimer and the DNA. The CD experiments carried out with a constant DNA concentration of 25 μM and dendrimer concentrations from 0 to 100 μM indicated minimal change to the chirality of the DNA for both types of dendrimers. While the PAMAM-G2 dendrimer caused aggregation of the majority of the DNA, the 2D NMR data of the DNA with an mPEG- b-PAMAM-G3 dendrimer indicated general broadening of the 1D ³¹P peaks from the DNA phosphates, a small number of ¹H chemical shift perturbations (CSPs), and reduction of specific ¹H-¹H NOE intensities. These data suggest there is minimal structural alteration of the DNA in the complex and indicate preferential binding of the dendrimer to the central AATT region of the DNA sequence. The results herein are the first such results demonstrating a soluble DNA complex with the mPEG- b-PAMAM-G3 dendrimer analyzed by multidimensional NMR.

Probing the influence of cell surface polysaccharides on nanodendrimer binding to Gram-negative and Gram-positive bacteria using single-nanoparticle force spectroscopy

The safe use and design of nanoparticles (NPs) ask for a comprehensive interpretation of their potentially adverse effects on (micro)organisms. In this respect, the prior assessment of the interactions experienced by NPs in the vicinity of - and in contact with - complex biological surfaces is mandatory. It requires the development of suitable techniques for deciphering the processes that govern nano-bio interactions when a single organism is exposed to an extremely low dose of NPs. Here, we used atomic force spectroscopy (AFM)-based force measurements to investigate at the nanoscale the interactions between carboxylate-terminated polyamidoamine (PAMAM) nanodendrimers (radius ca. 4.5 nm) and two bacteria with very distinct surface properties, Escherichia coli and Lactococcus lactis. The zwitterionic nanodendrimers exhibit a negative peripheral surface charge and/or a positive intraparticulate core depending on the solution pH and salt concentration. Following an original strategy according to which a single dendrimer NP is grafted at the very apex of the AFM tip, the density and localization of NP binding sites are probed at the surface of E. coli and L. lactis mutants expressing different cell surface structures (presence/absence of the O-antigen of the lipopolysaccharides (LPS) or of a polysaccharide pellicle). In line with electrokinetic analysis, AFM force measurements evidence that adhesion of NPs onto pellicle-decorated L. lactis is governed by their underlying electrostatic interactions as controlled by the pH-dependent charge of the peripheral and internal NP components, and the negatively-charged cell surface. In contrast, the presence of the O-antigen on E. coli systematically suppresses the adhesion of nanodendrimers onto cells, may the apparent NP surface charge be determined by the peripheral carboxylate groups or by the internal amine functions. Altogether, this work highlights the differentiated roles played by surface polysaccharides in mediating NP attachment to Gram-positive and Gram-negative bacteria. It further demonstrates that the assessment of NP bioadhesion features requires a critical analysis of the electrostatic contributions stemming from the various structures composing the stratified cell envelope, and those originating from the bulk and surface NP components. The joint use of electrokinetics and AFM provides a valuable option for rapidly addressing the binding propensity of NPs to microorganisms, as urgently needed in NP risk assessments.

Collapse of DNA in packaging and cellular transport

The dawn of molecular biology and recombinant DNA technology arose from our ability to manipulate DNA, including the process of collapse of long extended DNA molecules into nanoparticles of approximately 100 nm diameter. This condensation process is important for the packaging of DNA in the cell and for transporting DNA through the cell membrane for gene therapy. Multivalent cations, such as natural polyamines (spermidine and spermine), were initially recognized for their ability to provoke DNA condensation. Current research is targeted on molecules such as linear and branched polymers, oligopeptides, polypeptides and dendrimers that promote collapse of DNA to nanometric particles for gene therapy and on the energetics of DNA packaging.

Three RNA Microenvironments Detected in Fluxional Gene Delivery Polyplex Nanoassemblies

Prototropic and solvatochromatic properties of fluorescein (FL) were employed to detect the presence of microenvironments in polyplexes consisting of polycationic polymer (POCP) and a fluorescein-conjugated RNA, the HIV-1 transactivation response element (TAR-FL). Results reveal new aspects of polyplex structure with respect to polyplex-bound RNA existing in the following local microenvironments: (a) RNA associated with the polyplex that experiences local pH changes in a manner dependent on POCP nitrogen to RNA phosphate ratio (N:P), (b) RNA experiencing relatively acidic local pH environment that remains constant in polyplexes formed after a charge-neutral ratio, and (c) RNA packed close enough to mediate fluorophore/fluorophore quenching. The magnitude of these changes observed as a function of POCP to nucleic acid N:P ratio is polymer dependent. Assessment of the different microenvironments can help elucidate the functional hierarchy of polyplex-bound oligonucleotides and additionally characterize POCPs based on the resulting local pH and solvent properties upon polyplex formation.

Rapid Exchange Between Free and Bound States in RNA-Dendrimer Polyplexes: Implications on the Mechanism of Delivery and Release

A combination of solution NMR, dynamic light scattering (DLS), and fluorescence quenching assays were employed to obtain insights into the dynamics and structural features of a polyplex system consisting of HIV-1 transactivation response element (TAR) and PEGylated generation 5 poly(amidoamine) dendrimer (G5-PEG). NMR chemical shift mapping and (13)C spin relaxation based dynamics measurements depict the polyplex system as a highly dynamic assembly where the RNA, with its local structure and dynamics preserved, rapidly exchanges ( Collapse

Key Words

• therapeutic rna delivery
• pamam
• pegylation
• tar
• dls
• nmr spin relaxation

Collapse

MESH Headings

• Dendrimers/chemistry
• HIV-1/chemistry
• Magnetic Resonance Spectroscopy/methods
• Polyamines/chemistry
• Polyethylene Glycols/chemistry
• RNA/chemistry
• Response Elements/genetics
• Transfection/methods

Collapse

Grants

• R01 AI066975 NIAID NIH HHS
• P50 GM103297 NIGMS NIH HHS
• R01 AI087463-01 NIAID NIH HHS
• R01 EB005028 NIBIB NIH HHS
• R01EB005028 NIBIB NIH HHS
• R01 GM089846 NIGMS NIH HHS

Collapse

Affiliation(s)

Anisha Shakya Department of Chemistry, University of Michigan, Ann Arbor, MI 48109-1055, USA
Casey A. Dougherty Department of Chemistry, University of Michigan, Ann Arbor, MI 48109-1055, USA
Yi Xue Department of Biochemistry and Chemistry, Duke University Medical Center, Durham, NC 27710, USA
Hashim M. Al-Hashimi Department of Biochemistry and Chemistry, Duke University Medical Center, Durham, NC 27710, USA
Mark M. Banaszak Holl Department of Chemistry, University of Michigan, Ann Arbor, MI 48109-1055, USA

Collapse

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.

Long-Range Lipid-Water Interaction as Observed by ATR-FTIR Spectroscopy

It is commonly assumed that the structure of water at a lipid-water interface is influenced mostly in the first hydration layer. However, recent results from different experimental methods show that perturbation extends through several hydration layers. Due to its low light penetration depth, attenuated total reflection Fourier transform infrared (ATR-FTIR) spectroscopy is specifically suited to study interlamellar water structure in multibilayers. Results obtained by this technique confirm the long-range water structure disturbance. Consequently, in confined membrane environments nearly all water molecules can be perturbed. It is important to note that the behavior of confined water molecules differs significantly in samples prepared in excess water and in partially hydrated samples. We show in what manner the interlamellar water perturbation is influenced by the hydration level and how it is sequentially modified with a step-by-step dehydration of samples either by water evaporation or by osmotic pressure. Our results also indicate that besides different levels of hydration the lipid-water interaction is modulated by different lipid headgroups and different lipid phases as well. Therefore, modification of interlamellar water properties may clarify the role of water-mediated effects in biological processes.

Molecular Modeling to Study Dendrimers for Biomedical Applications

Molecular modeling techniques provide a powerful tool to study the properties of molecules and their interactions at the molecular level. The use of computational techniques to predict interaction patterns and molecular properties can inform the design of drug delivery systems and therapeutic agents. Dendrimers are hyperbranched macromolecular structures that comprise repetitive building blocks and have defined architecture and functionality. Their unique structural features can be exploited to design novel carriers for both therapeutic and diagnostic agents. Many studies have been performed to iteratively optimise the properties of dendrimers in solution as well as their interaction with drugs, nucleic acids, proteins and lipid membranes. Key features including dendrimer size and surface have been revealed that can be modified to increase their performance as drug carriers. Computational studies have supported experimental work by providing valuable insights about dendrimer structure and possible molecular interactions at the molecular level. The progress in computational simulation techniques and models provides a basis to improve our ability to better predict and understand the biological activities and interactions of dendrimers. This review will focus on the use of molecular modeling tools for the study and design of dendrimers, with particular emphasis on the efforts that have been made to improve the efficacy of this class of molecules in biomedical applications.

Molecular modeling of polynucleotide complexes

Delivery of polynucleotides into patient cells is a promising strategy for treatment of genetic disorders. Gene therapy aims to either synthesize desired proteins (DNA delivery) or suppress expression of endogenous genes (siRNA delivery). Carriers constitute an important part of gene therapeutics due to limitations arising from the pharmacokinetics of polynucleotides. Non-viral carriers such as polymers and lipids protect polynucleotides from intra and extracellular threats and facilitate formation of cell-permeable nanoparticles through shielding and/or bridging multiple polynucleotide molecules. Formation of nanoparticulate systems with optimal features, their cellular uptake and intracellular trafficking are crucial steps for an effective gene therapy. Despite the great amount of experimental work pursued, critical features of the nanoparticles as well as their processing mechanisms are still under debate due to the lack of instrumentation at atomic resolution. Molecular modeling based computational approaches can shed light onto the atomic level details of gene delivery systems, thus provide valuable input that cannot be readily obtained with experimental techniques. Here, we review the molecular modeling research pursued on critical gene therapy steps, highlight the knowledge gaps in the field and providing future perspectives. Existing modeling studies revealed several important aspects of gene delivery, such as nanoparticle formation dynamics with various carriers, effect of carrier properties on complexation, carrier conformations in endosomal stages, and release of polynucleotides from carriers. Rate-limiting steps related to cellular events (i.e. internalization, endosomal escape, and nuclear uptake) are now beginning to be addressed by computational approaches. Limitations arising from current computational power and accuracy of modeling have been hindering the development of more realistic models. With the help of rapidly-growing computational power, the critical aspects of gene therapy are expected to be better investigated and direct comparison between more realistic molecular modeling and experiments may open the path for design of next generation gene therapeutics.

Ion competition in condensed DNA arrays in the attractive regime

Physical origin of DNA condensation by multivalent cations remains unsettled. Here, we report quantitative studies of how one DNA-condensing ion (Cobalt(3+) Hexammine, or Co(3+)Hex) and one nonDNA-condensing ion (Mg(2+)) compete within the interstitial space in spontaneously condensed DNA arrays. As the ion concentrations in the bath solution are systematically varied, the ion contents and DNA-DNA spacings of the DNA arrays are determined by atomic emission spectroscopy and x-ray diffraction, respectively. To gain quantitative insights, we first compare the experimentally determined ion contents with predictions from exact numerical calculations based on nonlinear Poisson-Boltzmann equations. Such calculations are shown to significantly underestimate the number of Co(3+)Hex ions, consistent with the deficiencies of nonlinear Poisson-Boltzmann approaches in describing multivalent cations. Upon increasing the concentration of Mg(2+), the Co(3+)Hex-condensed DNA array expands and eventually redissolves as a result of ion competition weakening DNA-DNA attraction. Although the DNA-DNA spacing depends on both Mg(2+) and Co(3+)Hex concentrations in the bath solution, it is observed that the spacing is largely determined by a single parameter of the DNA array, the fraction of DNA charges neutralized by Co(3+)Hex. It is also observed that only ∼20% DNA charge neutralization by Co(3+)Hex is necessary for spontaneous DNA condensation. We then show that the bath ion conditions can be reduced to one variable with a simplistic ion binding model, which is able to describe the variations of both ion contents and DNA-DNA spacings reasonably well. Finally, we discuss the implications on the nature of interstitial ions and cation-mediated DNA-DNA interactions.

Intermolecular forces between low generation PAMAM dendrimer condensed DNA helices: role of cation architecture

In recent years, dendriplexes, complexes of cationic dendrimers with DNA, have become attractive DNA delivery vehicles due to their well-defined chemistries. To better understand the nature of the forces condensing dendriplexes, we studied low generation poly(amidoamine) (PAMAM) dendrimer-DNA complexes and compared them to comparably charged linear arginine peptides. Using osmotic stress coupled with X-ray scattering, we have investigated the effect of molecular chain architecture on DNA-DNA intermolecular forces that determine the net attraction and equilibrium interhelical distance within these polycation condensed DNA arrays. In order to compact DNA, linear cations are believed to bind in DNA grooves and to interact with the phosphate backbone of apposing helices. We have previously shown a length dependent attraction resulting in higher packaging densities with increasing charge for linear cations. Hyperbranched polycations, such as polycationic dendrimers, presumably would not be able to bind to DNA and correlate their charges in the same manner as linear cations. We show that attractive and repulsive force amplitudes in PAMAM-DNA assemblies display significantly different trends than comparably charged linear arginines resulting in lower DNA packaging densities with increasing PAMAM generation. The salt and pH dependencies of packaging in PAMAM dendrimer-DNA and linear arginine-DNA complexes were also investigated. Significant differences in the force curve behaviour and salt and pH sensitivities suggest that different binding modes may be present in DNA condensed by dendrimers when compared to linear polycations.