Marcelo Marucho

Biomolecular electrostatics and solvation: a computational perspective

An understanding of molecular interactions is essential for insight into biological systems at the molecular scale. Among the various components of molecular interactions, electrostatics are of special importance because of their long-range nature and their influence on polar or charged molecules, including water, aqueous ions, proteins, nucleic acids, carbohydrates, and membrane lipids. In particular, robust models of electrostatic interactions are essential for understanding the solvation properties of biomolecules and the effects of solvation upon biomolecular folding, binding, enzyme catalysis, and dynamics. Electrostatics, therefore, are of central importance to understanding biomolecular structure and modeling interactions within and among biological molecules. This review discusses the solvation of biomolecules with a computational biophysics view toward describing the phenomenon. While our main focus lies on the computational aspect of the models, we provide an overview of the basic elements of biomolecular solvation (e.g. solvent structure, polarization, ion binding, and non-polar behavior) in order to provide a background to understand the different types of solvation models.

Simulations of RNA interactions with monovalent ions.

RNA folding and binding reactions are mediated by interactions with ions that make up the surrounding aqueous electrolytic milieu. Although Mg(2+) ions are often implicated as being crucial for RNA folding, it is known that folding is feasible in high concentrations of monovalent alkali-halide salts. Experiments have yielded important information regarding the salt dependence of RNA stability. Recent work has shown that molecular simulations based on explicit representations of solvent molecules and monovalent ions can also provide useful insights regarding the ionic atmospheres around model RNA systems. These insights can help rationalize intriguing observations regarding the dependence of RNA stability on cation type providing one pays attention to important considerations that go into the proper design of molecular simulations. These issues are discussed in detail and the methods are applied to an A-form RNA and B-form DNA sequence. The results of these simulations are compared to previous work on a kissing-loop system with analogous sequence. In particular, ionic atmospheres obtained from molecular simulations are compared to those obtained using the nonlinear Poisson Boltzmann model for continuum electrostatics for these three different nucleic acid systems. The comparisons indicate reasonable agreement in terms of coarse-grained observables such as the numbers of counterions accumulated around the solutes. However, details of the ionic atmospheres, captured in terms of spatial free energy density profiles, are quite different between the two approaches. These comparisons suggest the need for improvements in continuum models to capture sequence-specific effects, ion-ion correlation, and the effects of partial dehydration of ions.

Ionic asymmetry and solvent excluded volume effects on spherical electric double layers: a density functional approach.

In this article, we present a classical density functional theory for electrical double layers of spherical macroions that extends the capabilities of conventional approaches by accounting for electrostatic ion correlations, size asymmetry, and excluded volume effects. The approach is based on a recent approximation introduced by Hansen-Goos and Roth for the hard sphere excess free energy of inhomogeneous fluids [J. Chem. Phys. 124, 154506 (2006); Hansen-Goos and Roth, J. Phys.: Condens. Matter 18, 8413 (2006)]. It accounts for the proper and efficient description of the effects of ionic asymmetry and solvent excluded volume, especially at high ion concentrations and size asymmetry ratios including those observed in experimental studies. Additionally, we utilize a leading functional Taylor expansion approximation of the ion density profiles. In addition, we use the mean spherical approximation for multi-component charged hard sphere fluids to account for the electrostatic ion correlation effects. These approximations are implemented in our theoretical formulation into a suitable decomposition of the excess free energy which plays a key role in capturing the complex interplay between charge correlations and excluded volume effects. We perform Monte Carlo simulations in various scenarios to validate the proposed approach, obtaining a good compromise between accuracy and computational cost. We use the proposed computational approach to study the effects of ion size, ion size asymmetry, and solvent excluded volume on the ion profiles, integrated charge, mean electrostatic potential, and ionic coordination number around spherical macroions in various electrolyte mixtures. Our results show that both solvent hard sphere diameter and density play a dominant role in the distribution of ions around spherical macroions, mainly for experimental water molarity and size values where the counterion distribution is characterized by a tight binding to the macroion, similar to that predicted by the Stern model.

Adsorption of Soft and Hard Proteins onto OTCEs under the Influence of an External Electric Field

The adsorption behavior of hard and soft proteins under the effect of an external electric field was investigated by a combination of spectroscopic ellipsometry and molecular dynamics (MD) simulations. Optically transparent carbon electrodes (OTCE) were used as conductive, sorbent substrates. Lysozyme (LSZ) and ribonuclease A (RNase A) were selected as representative hard proteins, whereas myoglobin (Mb), α-lactalbumin (α-LAC), bovine serum albumin (BSA), glucose oxidase (GOx), and immunoglobulin G (IgG) were selected to represent soft proteins. In line with recent publications from our group, the experimental results revealed that while the adsorption of all investigated proteins can be enhanced by the potential applied to the electrode, the effect is more pronounced for hard proteins. In contrast with the incomplete monolayers formed at open-circuit potential, the application of +800 mV to the sorbent surface induced the formation of multiple layers of protein. These results suggest that this effect can be related to the intrinsic polarizability of the protein (induction of dipoles), the resulting surface accessible solvent area (SASA), and structural rearrangements induced upon the incorporation on the protein layer. The described experiments are critical to understand the relationship between the structure of proteins and their tendency to form (under electric stimulation) layers with thicknesses that greatly surpass those obtained at open-circuit conditions.

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

The ionic atmosphere around a nucleic acid regulates its stability in aqueous salt solutions. One major source of complexity in biological activities involving nucleic acids arises from the strong influence of the surrounding ions and water molecules on their structural and thermodynamic properties. Here, we implement a classical density functional theory for cylindrical polyelectrolytes embedded in aqueous electrolytes containing explicit (neutral hard sphere) water molecules at experimental solvent concentrations. Our approach allows us to include ion correlations as well as solvent and ion excluded volume effects for studying the structural and thermodynamic properties of highly charged cylindrical polyelectrolytes. Several models of size and charge asymmetric mixtures of aqueous electrolytes at physiological concentrations are studied. Our results are in good agreement with Monte Carlo simulations. Our numerical calculations display significant differences in the ion density profiles for the different aqueous electrolyte models studied. However, similar results regarding the excess number of ions adsorbed to the B-DNA molecule are predicted by our theoretical approach for different aqueous electrolyte models. These findings suggest that ion counting experimental data should not be used alone to validate the performance of aqueous DNA-electrolyte models.

Biodirected Synthesis and Nanostructural Characterization of Anisotropic Gold Nanoparticles

Gold nanoparticles with anisotropic structures have tunable absorption properties and diverse bioapplications as image contrast agents, plasmonics, and therapeutic-diagnostic materials. Amino acids with electrostatically charged side chains possess inner affinity for metal ions. Lysine (Lys) efficiently controlled the growing into star-shape nanoparticles with controlled narrow sizes (30-100 nm) and produced in high yields (85-95%). Anisotropic nanostructures showed tunable absorbance from UV to NIR range, with extraordinary colloidal stability (-26 to -42 mV) and surface-enhanced Raman scattering properties. Advanced electron microscopy characterization through ultra-high-resolution SEM, STEM, and HR-TEM confirmed the size, nanostructure, crystalline structure, and chemical composition. Molecular dynamics simulations revealed that Lys interacted preferentially with Au(I) through the -COOH group instead of their positive side chains with a binding free energy (BFE) of 3.4 kcal mol(-1). These highly monodisperse and colloidal stable anisotropic particles prepared with biocompatible compounds may be employed in biomedical applications.

Adsorption and catalytic activity of glucose oxidase accumulated on OTCE upon the application of external potential

This article describes the adsorption of glucose oxidase (GOx) onto optically transparent carbon electrodes (OTCE) under the effect of applied potential and the analysis of the enzymatic activity of the resulting GOx/OTCE substrates. In order to avoid electrochemical interferences with the enzyme redox center, control electrochemical experiments were performed using flavin adenine dinucleotide (FAD) and GOx/OTCE substrates. Then, the enzyme adsorption experiments were carried out as a function of the potential applied (ranged from the open circuit potential to +950mV), the pH solution, the concentration of enzyme, and the ionic strength on the environment. The experimental results demonstrated that an increase in the adsorbed amount of GOx on the OTCE can be achieved when the potential was applied. Although the increase in the adsorbed amount was examined as a function of the potential, a maximum enzymatic activity was observed in the GOx/OTCE substrate achieved at +800mV. These experiments suggest that although an increase in the amount of enzyme adsorbed can be obtained by the application of an external potential to the electrode, the magnitude of such potential can produce detrimental effects in the conformation of the adsorbed protein and should be carefully considered. As such, the article describes a simple and rational approach to increase the amount of enzyme adsorbed on a surface and can be applied to improve the sensitivity of a variety of biosensors.

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

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.

MPBEC, a Matlab Program for Biomolecular Electrostatic Calculations

One of the most used and efficient approaches to compute electrostatic properties of biological systems is to numerically solve the Poisson-Boltzmann (PB) equation. There are several software packages available that solve the PB equation for molecules in aqueous electrolyte solutions. Most of these software packages are useful for scientists with specialized training and expertise in computational biophysics. However, the user is usually required to manually take several important choices, depending on the complexity of the biological system, to successfully obtain the numerical solution of the PB equation. This may become an obstacle for researchers, experimentalists, even students with no special training in computational methodologies. Aiming to overcome this limitation, in this article we present MPBEC, a free, cross-platform, open-source software that provides non-experts in the field an easy and efficient way to perform biomolecular electrostatic calculations on single processor computers. MPBEC is a Matlab script based on the Adaptative Poisson Boltzmann Solver, one of the most popular approaches used to solve the PB equation. MPBEC does not require any user programming, text editing or extensive statistical skills, and comes with detailed user-guide documentation. As a unique feature, MPBEC includes a useful graphical user interface (GUI) application which helps and guides users to configure and setup the optimal parameters and approximations to successfully perform the required biomolecular electrostatic calculations. The GUI also incorporates visualization tools to facilitate users pre- and post- analysis of structural and electrical properties of biomolecules.