We have developed a self-consistent description of an interface between a metal and a molecular liquid by combination of the density functional theory in the Kohn-Sham formulation ͑KS DFT͒ for the electronic structure, and the three-dimensional generalization of the reference interaction site model ͑3D RISM͒ for the classical site distribution profiles of liquid. The electron and classical subsystems are coupled in the mean field approximation. The procedure takes account of many-body effects of dense fluid on the metal-liquid interactions by averaging the pseudopotentials of liquid molecules over the classical distributions of the liquid. The proposed approach is substantially less time-consuming as compared to a Car-Parrinello-type simulation since it replaces molecular dynamics with the integral equation theory of molecular liquids. The calculation has been performed for pure water at normal conditions in contact with the ͑100͒ face cubic centered ͑fcc͒ surface of a metal roughly modeled after copper. The results are in good agreement with the Car-Parrinello simulation for the same metal model. The shift of the Fermi level due to the presence of water conforms with experiment. The electron distribution near an adsorbed water molecule is affected by dense water, and so the metal-water attraction follows the shapes of the metal effective electrostatic potential. For the metal model employed, it is strongest at the hollow site adsorption positions, and water molecules are adsorbed mainly at the hollow and bridge site positions rather than over metal atoms. Layering of water molecules near the metal surface is found. In the first hydration layer, adsorbed water molecules are oriented in parallel to the surface or tilted with hydrogens mainly outwards the metal. This orientation at the potential of zero charge agrees with experiment.

No abstract

We adapt the three-dimensional reference interaction site model (3D-RISM) to calculate the potentials of mean force for ion–molecular solution as a difference between the chemical potential of solvation of a cluster of solutes and of individual ones. The method yields the solvation structure around the cluster of solutes in detail. The solvation chemical potential is obtained for the three-dimensional hypernetted chain (3D-HNC) closure as well as for its partial linearization (3D-PLHNC approximation). The solvation chemical potential is obtained in a closed analytical form for both the 3D-HNC and 3D-PLHNC closures. The 3D-RISM integral equations are solved by using the supercell technique. A straightforward supercell treatment of ionic solute in polar molecular solvent leads to a big error in the potential of mean force as well as the solvation chemical potential, which for simple ions in water amounts to about 35 kcal/mol. We elaborated corrections to the 3D-RISM integral equations, alleviating the artifact of the supercell periodicity with an accuracy of 0.05 kcal/mol or better and restoring the long-range asymptotics of the solute–solvent correlation functions. The dielectrically consistent site–site RISM/HNC theory (DRISM/HNC) is employed for the solvent correlations to provide a proper description of the dielectric properties of solution. This allowed us to extend the description to solution at a finite salt concentration. We converge both the 3D-RISM and site–site DRISM integral equations by using the method of modified direct inversion in the iterative subspace. Owing to the proper initial guess of the correlation functions, iteration begins at once for a given temperature and full molecular charge, avoiding a gradual decrease of the temperature and increase of the site charges, which greatly reduces the computation time. We calculate and discuss the potentials of mean force for sodium chloride in ambient water at infinite dilution as well as at a finite concentration.

We present the three-dimensional molecular theory of solvation (also known as 3D-RISM) coupled with molecular dynamics (MD) simulation by contracting solvent degrees of freedom, accelerated by extrapolating solvent-induced forces and applying them in large multi-time steps (up to 20 fs) to enable simulation of large biomolecules. The method has been implemented in the Amber molecular modeling package, and is illustrated here on alanine dipeptide and protein G.

No abstract

We have modified the popular MM/PBSA or MM/GBSA approaches (molecular mechanics for a biomolecule, combined with a Poisson-Boltzmann or generalized Born electrostatic and surface area nonelectrostatic solvation energy) by employing instead the statistical-mechanical, three-dimensional molecular theory of solvation (also known as 3D reference interaction site model, or 3D-RISM-KH) coupled with molecular mechanics or molecular dynamics ( Blinov , N. ; et al. Biophys. J. 2010 ; Luchko , T. ; et al. J. Chem. Theory Comput. 2010 ). Unlike the PBSA or GBSA semiempirical approaches, the 3D-RISM-KH theory yields a full molecular picture of the solvation structure and thermodynamics from the first principles, with proper account of chemical specificities of both solvent and biomolecules, such as hydrogen bonding, hydrophobic interactions, salt bridges, etc. We test the method on the binding of seven biotin analogues to avidin in aqueous solution and show it to work well in predicting the ligand-binding affinities. We have compared the results of 3D-RISM-KH with four different generalized Born and two Poisson-Boltzmann methods. They give absolute binding energies that differ by up to 208 kJ/mol and mean absolute deviations in the relative affinities of 10-43 kJ/mol.

Tailoring the surface chemistry of metallic nanoparticles is generally a key step for their use in a wide range of applications. There are few examples of organic films covalently bound to metal nanoparticles. We demonstrate here that aryl films are formed on gold nanoparticles from the spontaneous reduction of diazonium salts. The structure and the bonding of the film is probed with surface-enhanced Raman scattering (SERS). Extinction spectroscopy and SERS show that a nitrobenzene film forms on gold nanoparticles from the corresponding diazonium salt. Comparison of the SERS spectrum with spectra computed from density functional theory models reveals a band characteristic of a Au-C stretch. The observation of this stretch is direct evidence of a covalent bond. A similar band is observed in high-resolution electron energy loss spectra of nitrobenzene layers on planar gold. The bonding of these types of films through a covalent interaction on gold is consistent with their enhanced stability observed in other studies. These findings provide motivation for the use of diazonium-derived films on gold and other metals in applications where high stability and/or strong adsorbate-substrate coupling are required.

We present a general-purpose model for biomolecular simulations at the molecular level that incorporates stochasticity, spatial dependence, and volume exclusion, using diffusing and reacting particles with physical dimensions. To validate the model, we first established the formal relationship between the microscopic model parameters (timestep, move length, and reaction probabilities) and the macroscopic coefficients for diffusion and reaction rate. We then compared simulation results with Smoluchowski theory for diffusion-limited irreversible reactions and the best available approximation for diffusion-influenced reversible reactions. To simulate the volumetric effects of a crowded intracellular environment, we created a virtual cytoplasm composed of a heterogeneous population of particles diffusing at rates appropriate to their size. The particle-size distribution was estimated from the relative abundance, mass, and stoichiometries of protein complexes using an experimentally derived proteome catalog from Escherichia coli K12. Simulated diffusion constants exhibited anomalous behavior as a function of time and crowding. Although significant, the volumetric impact of crowding on diffusion cannot fully account for retarded protein mobility in vivo, suggesting that other biophysical factors are at play. The simulated effect of crowding on barnase-barstar dimerization, an experimentally characterized example of a bimolecular association reaction, reveals a biphasic time course, indicating that crowding exerts different effects over different timescales. These observations illustrate that quantitative realism in biosimulation will depend to some extent on mesoscale phenomena that are not currently well understood.