Solvation effect on conformations of 1,2:dimethoxyethane: charge-dependent nonlinear response in implicit solvent models

Ion Solvation Free Energy Calculation Based on Ab Initio Molecular Dynamics Using a Hybrid Solvent Model

Free energy calculation of small molecules or ion species in aqueous solvent is one of the most important problems in electrochemistry study. Although there are many previous approaches to calculate such free energies, they are based on ab initio methods and suffer from various limitations and approximations. In the current work, we developed a hybrid approach based on ab initio molecular dynamics (AIMD) simulations to calculate the ion solvation energy. In this approach, a small water cluster surrounding the central ion is used, and implicit solvent model is applied outside the water cluster. A dynamic potential well is used during AIMD to keep the water cluster together. Quasi-harmonic approximation is used to calculate the entropy contribution, while the total energy average is used to calculate the enthalpy term. The obtained solvation voltages of the bulk metal agree with experiments within 0.3 eV, and the simulation results for the solvation energies of gaseous ions are close to the experimental observations. Besides the free energies, radial pair distribution functions and coordination numbers of hydrated cations are also obtained. The remaining challenges of this method are also discussed.

Explicit solvent theory of salt-induced dielectric decrement

Within the framework of an electrolyte model composed of structured solvent molecules and salt ions coupled by electrostatic and hard-core interactions, we characterize the physical mechanism behind salt-induced dielectric decrement.

Near-infrared dyes for two-photon absorption in the short-wavelength infrared: strategies towards optical power limiting

The recent advances in the field of two-photon absorbing chromophores in the short-wavelength infrared spectral range (SWIR 1100–2500 nm) are summarized, highlighting the development of optical power limiting devices in this spectral range.

Analytical Theory for Sequence-Specific Binary Fuzzy Complexes of Charged Intrinsically Disordered Proteins

Intrinsically disordered proteins (IDPs) are important for biological functions. In contrast to folded proteins, molecular recognition among certain IDPs is "fuzzy" in that their binding and/or phase separation are stochastically governed by the interacting IDPs' amino acid sequences, while their assembled conformations remain largely disordered. To help elucidate a basic aspect of this fascinating yet poorly understood phenomenon, the binding of a homo or heterodimeric pair of polyampholytic IDPs is modeled statistical mechanically using cluster expansion. We find that the binding affinities of binary fuzzy complexes in the model correlate strongly with a newly derived simple "joint sequence charge decoration" parameter readily calculable from the pair of IDPs' sequence charge patterns. Predictions by our analytical theory are in essential agreement with coarse-grained explicit-chain simulations. This computationally efficient theoretical framework is expected to be broadly applicable to rationalizing and predicting sequence-specific IDP-IDP polyelectrostatic interactions.

Hybrid Solvation Model with First Solvation Shell for Calculation of Solvation Free Energy

We present a hybrid solvation model with first solvation shell to calculate solvation free energies. This hybrid model combines the quantum mechanics and molecular mechanics methods with the analytical expression based on the Born solvation model to calculate solvation free energies. Based on calculated free energies of solvation and reaction profiles in gas phase, we set up a unified scheme to predict reaction profiles in solution. The predicted solvation free energies and reaction barriers are compared with experimental results for twenty bimolecular nucleophilic substitution reactions. These comparisons show that our hybrid solvation model can predict reliable solvation free energies and reaction barriers for chemical reactions of small molecules in aqueous solution.

A new lattice Monte Carlo simulation for dielectric saturation in ion-containing liquids

We develop a new, rapid method for the lattice Monte Carlo simulation of ion-containing liquids that accounts for the effects of the reorganization of solvent dipoles under external electrostatic fields. Our results are in reasonable agreement with the analytical solutions to the dielectric continuum theory of Booth for single ions, ion pairs, and ionic cross-links. We also illustrate the substantial disparity between the dielectric functions for like and unlike charges on the nanometer scale. Our simulation rationalizes the experimental data for the dependence of the bulk dielectric value of water on ion concentrations in terms of saturated dipoles near ions.

Communication: modeling charge-sign asymmetric solvation free energies with nonlinear boundary conditions

We show that charge-sign-dependent asymmetric hydration can be modeled accurately using linear Poisson theory after replacing the standard electric-displacement boundary condition with a simple nonlinear boundary condition. Using a single multiplicative scaling factor to determine atomic radii from molecular dynamics Lennard-Jones parameters, the new model accurately reproduces MD free-energy calculations of hydration asymmetries for: (i) monatomic ions, (ii) titratable amino acids in both their protonated and unprotonated states, and (iii) the Mobley "bracelet" and "rod" test problems [D. L. Mobley, A. E. Barber II, C. J. Fennell, and K. A. Dill, "Charge asymmetries in hydration of polar solutes," J. Phys. Chem. B 112, 2405-2414 (2008)]. Remarkably, the model also justifies the use of linear response expressions for charging free energies. Our boundary-element method implementation demonstrates the ease with which other continuum-electrostatic solvers can be extended to include asymmetry.

Modified Poisson–Boltzmann equations for characterizing biomolecular solvation

Methods for computing electrostatic interactions often account implicitly for the solvent, due to the much smaller number of degrees of freedom involved. In the Poisson–Boltzmann (PB) approach the electrostatic potential is obtained by solving the Poisson–Boltzmann equation (PBE), where the solvent region is modeled as a homogeneous medium with a high dielectric constant. PB however is not exempt of problems. It does not take into account for example the sizes of the ions in the atmosphere surrounding the solute, nor does it take into account the inhomogeneous dielectric response of water due to the presence of a highly charged surface. In this paper we review two major modifications of PB that circumvent these problems, namely the size-modified PB (SMPB) equation and the Dipolar Poisson–Boltzmann Langevin (DPBL) model. In SMPB, steric effects between ions are accounted for with a lattice gas model. In DPBL, the solvent region is no longer modeled as a homogeneous dielectric media but rather as an assembly of self-orienting interacting dipoles of variable density. This model results in a dielectric profile that transits smoothly from the solute to the solvent region as well as in a variable solvent density that depends on the charges of the solute. We show successful applications of the DPBL formalism to computing the solvation free energies of isolated ions in water. Further developments of more accurately modified PB models are discussed.

A novel implicit solvent model for simulating the molecular dynamics of RNA

Although molecular dynamics simulations can be accelerated by more than an order of magnitude by implicitly describing the influence of the solvent with a continuum model, most currently available implicit solvent simulations cannot robustly simulate the structure and dynamics of nucleic acids. The difficulties become exacerbated especially for RNAs, suggesting the presence of serious physical flaws in the prior continuum models for the influence of the solvent and counter ions on the nucleic acids. We present a novel, to our knowledge, implicit solvent model for simulating nucleic acids by combining the Langevin-Debye model and the Poisson-Boltzmann equation to provide a better estimate of the electrostatic screening of both the water and counter ions. Tests of the model involve comparisons of implicit and explicit solvent simulations for three RNA targets with 20, 29, and 75 nucleotides. The model provides reasonable agreement with explicit solvent simulations, and directions for future improvement are noted.

Gradient Models in Molecular Biophysics: Progress, Challenges, Opportunities

In the interest of developing a bridge between researchers modeling materials and those modeling biological molecules, we survey recent progress in developing nonlocal-dielectric continuum models for studying the behavior of proteins and nucleic acids. As in other areas of science, continuum models are essential tools when atomistic simulations (e.g. molecular dynamics) are too expensive. Because biological molecules are essentially all nanoscale systems, the standard continuum model, involving local dielectric response, has basically always been dubious at best. The advanced continuum theories discussed here aim to remedy these shortcomings by adding features such as nonlocal dielectric response, and nonlinearities resulting from dielectric saturation. We begin by describing the central role of electrostatic interactions in biology at the molecular scale, and motivate the development of computationally tractable continuum models using applications in science and engineering. For context, we highlight some of the most important challenges that remain and survey the diverse theoretical formalisms for their treatment, highlighting the rigorous statistical mechanics that support the use and improvement of continuum models. We then address the development and implementation of nonlocal dielectric models, an approach pioneered by Dogonadze, Kornyshev, and their collaborators almost forty years ago. The simplest of these models is just a scalar form of gradient elasticity, and here we use ideas from gradient-based modeling to extend the electrostatic model to include additional length scales. The paper concludes with a discussion of open questions for model development, highlighting the many opportunities for the materials community to leverage its physical, mathematical, and computational expertise to help solve one of the most challenging questions in molecular biology and biophysics.

Context and force field dependence of the loss of protein backbone entropy upon folding using realistic denatured and native state ensembles

The loss of conformational entropy is the largest unfavorable quantity affecting a protein's stability. We calculate the reduction in the number of backbone conformations upon folding using the distribution of backbone dihedral angles (ϕ,ψ) obtained from an experimentally validated denatured state model, along with all-atom simulations for both the denatured and native states. The average loss of entropy per residue is TΔS(BB)(U-N) = 0.7, 0.9, or 1.1 kcal·mol(-1) at T = 298 K, depending on the force field used, with a 0.6 kcal·mol(-1) dispersion across the sequence. The average equates to a decrease of a factor of 3-7 in the number of conformations available per residue (f = Ω(Denatured)/Ω(Native)) or to a total of f(tot) = 3(n)-7(n) for an n residue protein. Our value is smaller than most previous estimates where f = 7-20, that is, our computed TΔS(BB)(U-N) is smaller by 10-100 kcal mol(-1) for n = 100. The differences emerge from our use of realistic native and denatured state ensembles as well as from the inclusion of accurate local sequence preferences, neighbor effects, and correlated motions (vibrations), in contrast to some previous studies that invoke gross assumptions about the entropy in either or both states. We find that the loss of entropy primarily depends on the local environment and less on properties of the native state, with the exception of α-helical residues in some force fields.

Fast Analytical Methods for Macroscopic Electrostatic Models in Biomolecular Simulations

We review recent developments of fast analytical methods for macroscopic electrostatic calculations in biological applications, including the Poisson-Boltzmann (PB) and the generalized Born models for electrostatic solvation energy. The focus is on analytical approaches for hybrid solvation models, especially the image charge method for a spherical cavity, and also the generalized Born theory as an approximation to the PB model. This review places much emphasis on the mathematical details behind these methods.

Extended structures in RNA folding intermediates are due to nonnative interactions rather than electrostatic repulsion

RNA folding occurs via a series of transitions between metastable intermediate states for Mg(2+) concentrations below those needed to fold the native structure. In general, these folding intermediates are considerably less compact than their respective native states. Our previous work demonstrates that the major equilibrium intermediate of the 154-residue specificity domain (S-domain) of the Bacillus subtilis RNase P RNA is more extended than its native structure. We now investigate two models with falsifiable predictions regarding the origins of the extended intermediate structures in the S-domains of the B. subtilis and the Escherichia coli RNase P RNA that belong to different classes of P RNA and have distinct native structures. The first model explores the contribution of electrostatic repulsion, while the second model probes specific interactions in the core of the folding intermediate. Using small-angle X-ray scattering and Langevin dynamics simulations, we show that electrostatics plays only a minor role, whereas specific interactions largely account for the extended nature of the intermediate. Structural contacts in the core, including a nonnative base pair, help to stabilize the intermediate conformation. We conclude that RNA folding intermediates adopt extended conformations due to short-range, nonnative interactions rather than generic electrostatic repulsion of helical domains. These principles apply to other ribozymes and riboswitches that undergo functionally relevant conformational changes.

Structural parameters and vibrational spectra of a series of zinc meso-phenylporphyrins: a DFT and experimental study

The influences of meso-phenyl substitution on the geometric structure and vibrational spectra have been studied by DFT calculation (B3LYP/6-31G(d)) and experiment on a series of zinc porphyrins (ZnTPP: zinc 5,10,15,20-tetraphenylporphyrin; ZnTrPP: zinc 5,10,15-triphenylporphryin; ZnDPP: zinc 5,15-dipenylporphyirn; ZnMPP: zinc 5-monophenylporphyrin; ZnP: zinc porphine). Calculation indicates that meso-phenyl substitution gives rise to slight out-of-plane distortion but significant in-plane distortion, especially for the configuration around C(m) atom, to zinc porphyrin. The assignment of experimental vibrational spectra was proposed mainly on the basis of calculation. Different shifting tendency upon meso-phenyl substitution is observed for different structure-sensitive bands, such as the shifting of nu(2), nu(3), nu(6), and nu(8) modes toward higher frequencies and nu(4) and nu(28) modes toward lower frequencies, upon meso-phenyl substitution. This is attributed primarily to in-plane nuclear reorganization effect (IPNR), though the contribution from out-of-plane distortion cannot be excluded completely. Analysis on vibrational structure reveals that asymmetric meso-phenyl substitution, especially the 5,15-diphenyl substitution of ZnDPP, brings about asymmetric vibrational displacement, or even splitting of vibrational structure to normal modes involving mainly the motion of meso-C(m). This is ascribed to the reduction of symmetry of porphyrin skeleton caused by asymmetric meso-phenyl substitution.

Fluctuation in electrolyte solutions: the self energy

We address the issue of the self energy of the mobile ions in electrolyte solutions within a general Gaussian renormalized fluctuation theory using a field-theoretic approach. We introduce the Born radii of the ions in the form of a charge distribution allowing for different Born radii between the cations and anions. The model thus automatically yields a theory free of divergences and accounts for the solvation of the ions at the level of continuous dielectric media. In an inhomogeneous dielectric medium, the self energy is in general position dependent and differences in the self energy between cations and anions can give rise to local charge separation in a macroscopically neutral system. Treating the Born radius a as a smallness parameter, we show that the self energy can be split into an O(a(-1)) nonuniversal contribution and an O(a0) universal contribution that depends only on the ion concentration, valency, and the spatially varying dielectric constant. For a weakly inhomogeneous dielectric medium, the nonuniversal part of the self energy is shown to have the form of the Born energy with the local dielectric constant. This self energy is incorporated into the Poisson-Boltzmann equation as a simple means of including this local fluctuation effect in a mean-field theory. We illustrate the phenomenon of charge separation by considering cations and anions of difference sizes and valencies in a periodic dielectric medium.

A single molecule as a dielectric medium

For three molecules with weak or negligible charge overlap, we prove that the three-body interaction energy obtained from quantum perturbation theory (to leading order) fits a dielectric model with a nonlocal electronic screening function. The electronic charge cloud of each molecule acts as a dielectric medium for the interaction of the remaining two with the nonlocal dielectric function epsilon(r,r') obtained by O. S. Jenkins and K. L. C. Hunt [J. Chem. Phys. 119, 8250 (2003)], by considering the charge redistribution induced in a single molecule by an external perturbation. The dielectric function depends parametrically on the coordinates of the nuclei, within the Born-Oppenheimer approximation. We also prove that the force on each nucleus in molecule A depends on intramolecular dielectric screening within A. The potential from the charge distribution of B, screened by C acting as a dielectric medium, is further screened linearly within A; and similarly, with the roles of B and C reversed. In addition, the potential due to the unperturbed charge distribution of B and the potential due to the unperturbed charge distribution of C, acting simultaneously, are screened nonlinearly within A. The results show that nonlocal dielectric theory holds on the molecular level, provided that the overlap of the electronic charge distributions is weak.

Image Charge Methods for a Three-Dielectric-Layer Hybrid Solvation Model of Biomolecules

This paper introduces a three dielectric layer hybrid solvation model for treating electrostatic interactions of biomolecules in solvents using the Poisson-Boltzmann equation. In this model, an interior spherical cavity will contain the solute and some explicit solvent molecules, and an intermediate buffer layer and an exterior layer contain the bulk solvent. A special dielectric permittivity profile is used to achieve a continuous dielectric transition from the interior cavity to the exterior layer. The selection of this special profile using a harmonic interpolation allows an analytical solution of the model by generalizing the classical Kirkwood series expansion. Discrete image charges are used to speed up calculations for the electrostatic potential within the interior and buffer layer regions. Semi-analytical and least squares methods are used to construct an accurate discrete image approximation for the reaction field due to solvent with or without salt effects. In particular, the image charges obtained by the least squares method provide accurate approximations to the reaction field independent of the ionic concentration of the solvent. Numerical results are presented to validate the accuracy and effectiveness of the image charge methods.

Computing ion solvation free energies using the dipolar Poisson model

A new continuum model is presented for computing the solvation free energies of cations in water. It combines in a single formalism based on statistical thermodynamics the Poisson model for electrostatics with the Langevin dipole model to account for nonuniform water dipole distribution around the ions. An excellent match between experimental and computed solvation free energies is obtained for 10 monovalent and divalent ions.

Langevin-Debye model for nonlinear electrostatic screening of solvated ions

Ion-ion interactions are central to numerous phenomena in geology, biology, and material science, but current understanding is based on linear theories of limited physical applicability. Using the Langevin-Debye model, we investigate how nonlinear dielectric saturation alters screening of electrostatic interactions between solvated charges and show that charge screening in liquids strongly depends on the charge magnitudes but negligibly on the ion radii. Qualitatively different universal screening functions are predicted for ions with charges of the same or opposite signs.

Quantifying the structural requirements of the folding transition state of protein A and other systems

The B-domain of protein A is a small three-helix bundle that has been the subject of considerable experimental and theoretical investigation. Nevertheless, a unified view of the structure of the transition-state ensemble (TSE) is still lacking. To characterize the TSE of this surprisingly challenging protein, we apply a combination of psi analysis (which probes the role of specific side-chain to side-chain contacts) and kinetic H/D amide isotope effects (which measures hydrogen-bond content), building upon previous studies using mutational phi analysis (which probes the energetic influence of side-chain substitutions). The second helix is folded in the TSE, while helix formation appears just at the carboxy and amino termini of the first and third helices, respectively. The experimental data suggest a homogenous yet plastic TS with a native-like topology. This study generalizes our earlier conclusion, based on two larger alpha/beta proteins, that the TSEs of most small proteins achieve approximately 70% of their native state's relative contact order. This high percentage limits the degree of possible TS heterogeneity and requires a reevaluation of the structural content of the TSE of other proteins, especially when they are characterized as small or polarized.

Influence of nonlinear electrostatics on transfer energies between liquid phases: charge burial is far less expensive than Born model

The widely used Born model describes the electrostatic response of continuous media using static dielectric constants. However, when applied to a liquid environment, a comparison of Born model predictions with experimental values (e.g., transfer free energies and pK(a) shifts) found that agreement is only achieved by using physically unrealistic dielectric constants for proteins, lipids, etc., and/or equally unrealistic atomic radii. This leads to questions concerning the physical origins for this failure of the Born model. We partially resolve this question by applying the Langevin-Debye (LD) model of a continuous distribution of point, polarizable dipoles, a model that contains an added dependence of the electrostatic response on the solvent's optical dielectric constant and both gas- and liquid-phase dipole moments, features absent in the Born model to which the LD model reduces for weak fields. The LD model is applied to simple representations of three biologically relevant systems: (i) globular proteins, (ii) lipid bilayers, and (iii) membrane proteins. The linear Born treatment greatly overestimates both the self-energy and the transfer free energy from water to hydrophobic environments (e.g., a protein interior). By using the experimental dielectric constant, the energy cost of charge burial in either globular or membrane proteins of the Born model is reduced by almost 50% with the nonlinear theory as is the pK(a) shift, and the shifts agree well with experimental trends.

Benchmarking implicit solvent folding simulations of the amyloid beta(10-35) fragment

A pathogenetic feature of Alzhemier disease is the aggregation of monomeric beta-amyloid proteins (Abeta) to form oligomers. Usually these oligomers of long peptides aggregate on time scales of microseconds or longer, making computational studies using atomistic molecular dynamics models prohibitively expensive and making it essential to develop computational models that are cheaper and at the same time faithful to physical features of the process. We benchmark the ability of our implicit solvent model to describe equilibrium and dynamic properties of monomeric Abeta(10-35) using all-atom Langevin dynamics (LD) simulations, since Alphabeta(10-35) is the only fragment whose monomeric properties have been measured. The accuracy of the implicit solvent model is tested by comparing its predictions with experiment and with those from a new explicit water MD simulation, (performed using CHARMM and the TIP3P water model) which is approximately 200 times slower than the implicit water simulations. The dependence on force field is investigated by running multiple trajectories for Alphabeta(10-35) using the CHARMM, OPLS-aal, and GS-AMBER94 force fields, whereas the convergence to equilibrium is tested for each force field by beginning separate trajectories from the native NMR structure, a completely stretched structure, and from unfolded initial structures. The NMR order parameter, S2, is computed for each trajectory and is compared with experimental data to assess the best choice for treating aggregates of Alphabeta. The computed order parameters vary significantly with force field. Explicit and implicit solvent simulations using the CHARMM force fields display excellent agreement with each other and once again support the accuracy of the implicit solvent model. Alphabeta(10-35) exhibits great flexibility, consistent with experiment data for the monomer in solution, while maintaining a general strand-loop-strand motif with a solvent-exposed hydrophobic patch that is believed to be important for aggregation. Finally, equilibration of the peptide structure requires an implicit solvent LD simulation as long as 30 ns.