A highly parallelizable integral equation theory for three dimensional solvent distribution function: application to biomolecules.

Three dimensional (3D) hydration structure is informative to clarify the functions of hydrated waters around a protein. We develop a new approach to calculate 3D solvation structure with reasonable computational cost. In the present method, the total solvation structure is obtained using conventional one dimensional reference interaction site model (RISM) followed by integrating the 3D fragment data, which are evaluated around each atom (site) of solute. Thanks to this strategy, time-consuming 3D fast Fourier transformation, which is required in 3D-RISM theory, can be avoided and high-parallel performance is achieved. The method is applied to small molecular systems for comparison with 3D-RISM. The obtained results by the present method and by 3D-RISM show good agreement. The hydration structures for a large protein computed by the present method are also consistent with those obtained by x-ray crystallography.

[1]  Masahiro Kinoshita,et al.  Theoretical study for partial molar volume of amino acids in aqueous solution: Implication of ideal fluctuation volume , 2000 .

[2]  S. Sakaki,et al.  An integral equation theory for 3D solvation structure: A new procedure free from 3D Fourier transform , 2006 .

[3]  H Luecke,et al.  Structure of bacteriorhodopsin at 1.55 A resolution. , 1999, Journal of molecular biology.

[4]  M. P. Sears,et al.  A new efficient method for density functional theory calculations of inhomogeneous fluids , 2003 .

[5]  Shinya Honda,et al.  10 residue folded peptide designed by segment statistics. , 2004, Structure.

[6]  David Chandler,et al.  Optimized Cluster Expansions for Classical Fluids. II. Theory of Molecular Liquids , 1972 .

[7]  Fumio Hirata,et al.  Water molecules in a protein cavity detected by a statistical-mechanical theory. , 2005, Journal of the American Chemical Society.

[8]  Maxim V. Fedorov,et al.  Wavelet algorithm for solving integral equations of molecular liquids. A test for the reference interaction site model , 2004, J. Comput. Chem..

[9]  K. Nitta,et al.  Partial Molar Volumes and Adiabatic Compressibilities of Amino Acids in Dilute Aqueous Solutions at 5, 15, 25, 35, and 45 .degree.C , 1995 .

[10]  Benoît Roux,et al.  An Integral Equation To Describe the Solvation of Polar Molecules in Liquid Water , 1997 .

[11]  V. Lebedev,et al.  A QUADRATURE FORMULA FOR THE SPHERE OF THE 131ST ALGEBRAIC ORDER OF ACCURACY , 1999 .

[12]  D. Lévesque,et al.  Theoretical determination of the dielectric constant of a two dimensional dipolar fluid , 1981 .

[13]  Fumio Hirata,et al.  Self-consistent description of a metal–water interface by the Kohn–Sham density functional theory and the three-dimensional reference interaction site model , 1999 .

[14]  K Schulten,et al.  VMD: visual molecular dynamics. , 1996, Journal of molecular graphics.

[15]  H. Berendsen,et al.  Interaction Models for Water in Relation to Protein Hydration , 1981 .

[16]  N. Shimba,et al.  The pH-dependent structural variation of complementarity-determining region H3 in the crystal structures of the Fv fragment from an anti-dansyl monoclonal antibody. , 1999, Journal of molecular biology.

[17]  J. D. Talman,et al.  Numerical Fourier and Bessel transforms in logarithmic variables , 1978 .

[18]  平田 文男.  Molecular theory of solvation , 2003 .

[19]  Fumio Hirata,et al.  Molecular Theory of Solvation , 2004 .

[20]  Fumio Hirata,et al.  An extended rism equation for molecular polar fluids , 1981 .

[21]  H. Kandori.  Role of internal water molecules in bacteriorhodopsin. , 2000, Biochimica et biophysica acta.

[22]  W. L. Jorgensen,et al.  The OPLS [optimized potentials for liquid simulations] potential functions for proteins, energy minimizations for crystals of cyclic peptides and crambin. , 1988, Journal of the American Chemical Society.

[23]  I. Shimada,et al.  Conformational dynamics of complementarity-determining region H3 of an anti-dansyl Fv fragment in the presence of its hapten. , 2005, Journal of molecular biology.

[24]  M. Mizuguchi,et al.  Partial molar volumes and adiabatic compressibilities ofN-acetyl-DL-serinamide andN-acetyl-L-threonmamide in dilute aqueous solutions , 1997 .

[25]  A. Becke.  A multicenter numerical integration scheme for polyatomic molecules , 1988 .

[26]  K. Nitta,et al.  Partial Molar Volumes and Isentropic Compressibilities of Amino Acids in Dilute Aqueous Solutions , 1998 .

[27]  F. Hirata,et al.  Three-dimensional density profiles of water in contact with a solute of arbitrary shape: a RISM approach , 1998 .

[28]  P. Calmettes,et al.  Molecular solvent model for a cryptate solution in acetonitrile: A hypernetted chain study , 1994 .

[29]  Rob H. Bisseling,et al.  Parallel scientific computation - a structured approach using BSP and MPI , 2004 .

[30]  A. Zamyatnin,et al.  Amino acid, peptide, and protein volume in solution. , 1984, Annual review of biophysics and bioengineering.

[31]  Fumio Hirata,et al.  Potentials of mean force of simple ions in ambient aqueous solution. I. Three-dimensional reference interaction site model approach , 2000 .

[32]  Norio Yoshida,et al.  Selective ion binding by protein probed with the statistical mechanical integral equation theory. , 2007, The journal of physical chemistry. B.

[33]  Yuichi Harano,et al.  Theoretical study for partial molar volume of amino acids and polypeptides by the three-dimensional reference interaction site model , 2001 .

[34]  Rob H. Bisseling,et al.  Parallel Scientific Computation , 2004 .