Excluded volume and ion-ion correlation effects on the ionic atmosphere around B-DNA: theory, simulations, and experiments

Theory and simulation of ligand functionalized nanoparticles - a pedagogical overview

Synthesizing reconfigurable nanoscale synthons with predictive control over shape, size, and interparticle interactions is a holy grail of bottom-up self-assembly. Grand challenges in their rational design, however, lie in both the large space of experimental synthetic parameters and proper understanding of the molecular mechanisms governing their formation. As such, computational and theoretical tools for predicting and modeling building block interactions have grown to become integral in modern day self-assembly research. In this review, we provide an in-depth discussion of the current state-of-the-art strategies available for modeling ligand functionalized nanoparticles. We focus on the critical role of how ligand interactions and surface distributions impact the emergent, pre-programmed behaviors between neighboring particles. To help build insights into the underlying physics, we first define an "ideal" limit - the short ligand, "hard" sphere approximation - and discuss all experimental handles through the lens of perturbations about this reference point. Finally, we identify theories that are capable of bridging interparticle interactions to nanoscale self-assembly and conclude by discussing exciting new directions for this field.

Java application for cytoskeleton filament characterization (JACFC)

JACFC is a Java web application (http://neuronanobiophysics.utsa.edu/) that provides both experts and non-experts in the field suitable tools for elucidating the molecular mechanisms modulating the electrical signal propagation, stability, and bundle formation of microtubules and actin filaments under different molecular (wild type, isoforms, mutants) and environmental (physiological and pathological) conditions. This acknowledgment might reveal the potential role of cytoskeleton filaments in neuronal activities, including molecular-level processing of information and neural regeneration. Molecular understanding of the polyelectrolyte properties of bionanowires, is also crucial for development of reliability, highly functioning small devices with biotechnological applications such as bionanosensors and computing bionanoprocessors.

Size and charge correlations in spherical electric double layers: a case study with fully asymmetric mixed electrolytes within the solvent primitive model

Size and charge correlations in spherical electric double layers are investigated through Monte Carlo simulations and density functional theory, through a solvent primitive model representation. A fully asymmetric mixed electrolyte is used for the small ions, whereas the solvent, apart from being a continuum dielectric, is also treated as an individual component. A partially perturbative density functional theory is adopted here, and for comparison, a standard canonical ensemble Monte Carlo simulation is used. The hard-sphere free energy is treated within a weighted density approach and the residual ionic contribution is estimated through perturbation around the uniform density. The results from both methods corroborate each other quantitatively over a wide range of physical parameters. The importance of structural correlations is envisaged through the size and charge asymmetry of the supporting electrolytes that includes the solvent as a component.

Size and charge correlations in spherical electric double layers are investigated through Monte Carlo simulations and density functional theory, through a solvent primitive model representation.

The Self-Adaptation Ability of Zinc Oxide Nanoparticles Enables Reliable Cancer Treatments

Optimal procedures for reliable anti-cancer treatments involve the systematic delivery of zinc oxide nanoparticles, which spread through the circulatory system. The success of these procedures may largely depend on the NPs' ability of self-adapting their physicochemical properties to overcome the different challenges facing at each stage on its way to the interior of a cancerous cell. In this article, we combine a multiscale approach, a unique nanoparticle model, and available experimental data to characterize the behavior of zinc oxide nanoparticles under different vessels rheology, pH levels, and biological environments. We investigate their ability to prevent aggregation, allow prolonged circulation time in the bloodstream, avoid clearance, conduct themselves through the capillarity system to reach damaged tissues, and selectively approach to target cancerous cells. Our results show that non-functionalized spherical zinc oxide nanoparticles with surface density N = 5.89 × 10-6 mol/m2, protonation and deprotonation rates pKa = 10.9 and pKb = -5.5, and NP size in the range of 20-50 nm are the most effective, smart anti-cancer agents for biomedical treatments.

A Java Application to Characterize Biomolecules and Nanomaterials in Electrolyte Aqueous Solutions

The electrostatic, entropic and surface interactions between a macroion (nanoparticle or biomolecule), surrounding ions and water molecules play a fundamental role in the behavior and function of colloidal systems. However, the molecular mechanisms governing these phenomena are still poorly understood. One of the major limitations in procuring this understanding is the lack of appropriate computational tools. Additionally, only experts in the field with an extensive background in programming, who are trained in statistical mechanics, and have access to supercomputers are able to study these systems. To overcome these limitations, in this article, we present a free, multiplatform, portable Java software, which provides experts and non-experts in the field an easy and efficient way to obtain an accurate molecular characterization of electrical and structural properties of aqueous electrolyte mixture solutions around both cylindrical- and spherical-like rigid macroions under multiple conditions. These properties include the normalized ions and water density profile distributions, the mean electrostatic potential, the integrated charge, the zeta potential, the electrostatic potential energy, the particle crowding entropy energy, the ion-ion electrostatic direct correlation energy, and the ionic potential of mean force. The Java software does not require outstanding skills and comes with detailed user-guide documentation. The application is based on the so-called Classical Density Functional Theory Solver (CSDFTS), which was successfully applied to a variety of rod-like biopolymers, rigid-like globular proteins, nanoparticles, and nano-rods. CSDFTS implements several electrolyte and macroion models, uses different levels of approximation and takes advantage of high performance Fortran90 routines and optimized libraries. These features enable the software to run on single processor computers at low-to-moderate computational cost depending on the computer performance, the grid resolution, and the characterization of the macroion and the electrolyte solution, among other factors. As a unique feature, the software comes with a graphical user interface (GUI) that allows users to take advantage of the visually guided setup of the required input data to properly characterize the system and configure the solver. Several examples on nanomaterials and biomolecules are provided to illustrate the use of the GUI and the solver performance.

Predicting Monovalent Ion Correlation Effects in Nucleic Acids

Ion correlation and fluctuation can play a potentially significant role in metal ion-nucleic acid interactions. Previous studies have focused on the effects for multivalent cations. However, the correlation and fluctuation effects can be important also for monovalent cations around the nucleic acid surface. Here, we report a model, gMCTBI, that can explicitly treat discrete distributions of both monovalent and multivalent cations and can account for the correlation and fluctuation effects for the cations in the solution. The gMCTBI model enables investigation of the global ion binding properties as well as the detailed discrete distributions of the bound ions. Accounting for the ion correlation effect for monovalent ions can lead to more accurate predictions, especially in a mixed monovalent and multivalent salt solution, for the number and location of the bound ions. Furthermore, although the monovalent ion-mediated correlation does not show a significant effect on the number of bound ions, the correlation may enhance the accumulation of monovalent ions near the nucleic acid surface and hence affect the ion distribution. The study further reveals novel ion correlation-induced effects in the competition between the different cations around nucleic acids.

A multi-scale approach to describe electrical impulses propagating along actin filaments in both intracellular and in vitro conditions

An accurate and efficient characterization of the polyelectrolyte properties for cytoskeleton filaments are key to the molecular understanding of electrical signal propagation, bundle and network formation, as well as their potential nanotechnological applications. In this article, we introduce an innovative multi-scale approach able to account for the atomistic details of a protein molecular structure, its biological environment, and their impact on electrical impulses propagating along wild type actin filaments. The formulation includes non-trivial contributions to the ionic electrical conductivity and capacitance coming from the diffuse part of the electrical double layer of G-actins. We utilize this monomer characterization in a non-linear inhomogeneous transmission line prototype model to account for the monomer–monomer interactions, dissipation and damping perturbations along the filament length. A novel, simple, accurate, approximate analytic expression has been obtained for the transmission line model. Our results reveal the propagation of electrical signal impulses in the form of solitons for the range of voltage stimulus and electrolyte solutions typically present for intracellular and in vitro conditions. The approach predicts a lower electrical conductivity with higher linear capacitance and non-linear accumulation of charge for intracellular conditions. Our results show a significant influence of the voltage input on the electrical impulse shape, attenuation and kern propagation velocity. The filament is able to sustain the soliton propagation at almost constant kern velocity for the in vitro condition, whereas the intracellular condition displays a remarkable deceleration. Additionally, the solitons are narrower and travel faster at higher voltage input. As a unique feature, this multi-scale theory is able to account for molecular structure conformation (mutation) and biological environment (protonations/deprotonations) changes often present in pathological conditions. It is also applicable to other highly charged rod-like polyelectrolytes with relevance in biomedicine and biophysics.

An innovative analytic solution accounting for the molecular structure, its biological environment, and their impact on electrical impulses along microfilaments.

Electrical double layer properties of spherical oxide nanoparticles

The accurate characterization of the electrical double layer properties of nanoparticles is of fundamental importance for optimizing their physicochemical properties for specific biotechnological and biomedical applications. In this article, we use classical solvation density functional theory and a surface complexation model to investigate the effects of the pH and the nanoparticle size on the structural and electrostatic properties of an electrolyte solution surrounding a spherical silica oxide nanoparticle. The formulation has been particularly useful for identifying dominant interactions governing the ionic driving force at a variety of pH levels and nanoparticle sizes. As a result of the energetic interplay displayed between electrostatic potential, ion-ion correlation and particle crowding effects on the nanoparticle surface titration, rich, non-trivial ion density profiles and mean electrostatic potential behavior have been found.

The Role of Correlation and Solvation in Ion Interactions with B-DNA

The ionic atmospheres around nucleic acids play important roles in biological function. Large-scale explicit solvent simulations coupled to experimental assays such as anomalous small-angle x-ray scattering can provide important insights into the structure and energetics of such atmospheres but are time- and resource intensive. In this article, we use classical density functional theory to explore the balance among ion-DNA, ion-water, and ion-ion interactions in ionic atmospheres of RbCl, SrCl2, and CoHexCl3 (cobalt hexamine chloride) around a B-form DNA molecule. The accuracy of the classical density functional theory calculations was assessed by comparison between simulated and experimental anomalous small-angle x-ray scattering curves, demonstrating that an accurate model should take into account ion-ion correlation and ion hydration forces, DNA topology, and the discrete distribution of charges on the DNA backbone. As expected, these calculations revealed significant differences among monovalent, divalent, and trivalent cation distributions around DNA. Approximately half of the DNA-bound Rb(+) ions penetrate into the minor groove of the DNA and half adsorb on the DNA backbone. The fraction of cations in the minor groove decreases for the larger Sr(2+) ions and becomes zero for CoHex(3+) ions, which all adsorb on the DNA backbone. The distribution of CoHex(3+) ions is mainly determined by Coulomb and steric interactions, while ion-correlation forces play a central role in the monovalent Rb(+) distribution and a combination of ion-correlation and hydration forces affect the Sr(2+) distribution around DNA. This does not imply that correlations in CoHex solutions are weaker or stronger than for other ions. Steric inaccessibility of the grooves to large CoHex ions leads to their binding at the DNA surface. In this binding mode, first-order electrostatic interactions (Coulomb) dominate the overall binding energy as evidenced by low sensitivity of ionic distribution to the presence or absence of second-order electrostatic correlation interactions.

Polar-solvation classical density-functional theory for electrolyte aqueous solutions near a wall

A precise description of the structural and dielectric properties of liquid water is critical to understanding the physicochemical properties of solutes in electrolyte solutions. In this article, a mixture of ionic and dipolar hard spheres is considered to account for water crowding and polarization effects on ionic electrical double layers near a uniformly charged hard wall. As a unique feature, solvent hard spheres carrying a dipole at their centers were used to model water molecules at experimentally known concentration, molecule size, and dipolar moment. The equilibrium ionic and dipole density profiles of this electrolyte aqueous model were calculated using a polar-solvation classical density-functional theory (PSCDFT). These profiles were used to calculate the charge density distribution, water polarization, dielectric permittivity function, and mean electric potential profiles as well as differential capacitance, excess adsorptions, and wall-fluid surface tension. These results were compared with those corresponding to the pure dipolar model and unpolar primitive solvent model of electrolyte aqueous solutions to understand the role that water crowding and polarization effects play on the structural and thermodynamic properties of these systems. Overall, PSCDFT predictions are in agreement with available experimental data.

Ion-ion correlation, solvent excluded volume and pH effects on physicochemical properties of spherical oxide nanoparticles

One major source of complexity in the implementation of nanoparticles in aqueous electrolytes arises from the strong influence that biological environments has on their physicochemical properties. A key parameter for understanding the molecular mechanisms governing the physicochemical properties of nanoparticles is the formation of the surface charge density. In this article, we present an efficient and accurate approach that combines a recently introduced classical solvation density functional theory for spherical electrical double layers with a surface complexation model to account for ion-ion correlation and excluded volume effects on the surface titration of spherical nanoparticles. We apply the proposed computational approach to account for the charge-regulated mechanisms on the surface chemistry of spherical silica (SiO2) nanoparticles. We analyze the effects of the nanoparticle size, as well as pH level and electrolyte concentration of the aqueous solution on the nanoparticle's surface charge density and Zeta potential. We validate our predictions for 580Å and 200Å nanoparticles immersed in acid, neutral and alkaline mono-valent aqueous electrolyte solutions against experimental data. Our results on mono-valent electrolyte show that the excluded volume and ion-ion correlations contribute significantly to the surface charge density and Zeta potential of the nanoparticle at high electrolyte concentration and pH levels, where the solvent crowding effects and electrostatic screening have shown a profound influence on the protonation/deprotonation reactions at the liquid/solute interface. The success of this approach in describing physicochemical properties of silica nanoparticles supports its broader application to study other spherical metal oxide nanoparticles.

Multivalent ion-mediated nucleic acid helix-helix interactions: RNA versus DNA

Ion-mediated interaction is critical to the structure and stability of nucleic acids. Recent experiments suggest that the multivalent ion-induced aggregation of double-stranded (ds) RNAs and DNAs may strongly depend on the topological nature of helices, while there is still lack of an understanding on the relevant ion-mediated interactions at atomistic level. In this work, we have directly calculated the potentials of mean force (PMF) between two dsRNAs and between two dsDNAs in Co(NH3)6 (3+) (Co-Hex) solutions by the atomistic molecular dynamics simulations. Our calculations show that at low [Co-Hex], the PMFs between B-DNAs and between A-RNAs are both (strongly) repulsive. However, at high [Co-Hex], the PMF between B-DNAs is strongly attractive, while those between A-RNAs and between A-DNAs are still (weakly) repulsive. The microscopic analyses show that for A-form helices, Co-Hex would become 'internal binding' into the deep major groove and consequently cannot form the evident ion-bridge between adjacent helices, while for B-form helices without deep grooves, Co-Hex would exhibit 'external binding' to strongly bridge adjacent helices. In addition, our further calculations show that, the PMF between A-RNAs could become strongly attractive either at very high [Co-Hex] or when the bottom of deep major groove is fixed with a layer of water.