Search

The three-deimensional distribution functions (3D-DFs) and potentials of mean force (PMFs) of small neutral molecules inside the two aquaporin channels, AQP1 and GlpF, are calculated based on the 3D-RISM theory, the statistical mechanics theory of molecular liquids, in order to investigate the permeability of those ligands through the channels. The ligands investigated are neon (Ne), carbon dioxide (CO2), nitric oxide (NO), ammonia (NH3), urea, and glycerol. Neon shows continuous distribution throughout the channel pore in AQP1 as is the case of water, although the PMF of Ne at the selective filter (SF) region is higher than that of water, indicating that the stability of molecules in the channel is determined not only by their size, but also by the charge distribution. The ligand molecules, CO2, NO, urea, and glycerol, have a large barrier in PMF at the SF region in AQP1, indicating that the channel is not permeable by those ligands. On the other hand, NH3has only a small activation barrier, ∼2.5 kJ/mol, to be overcome. Therefore, our theory predicts that a NH3molecule can be permeated through the AQP1 channel. In GlpF, all the ligands have negative PMF throughout the channel pore except for glycerol, which has a small barrier at the SF area, ∼2.1 kJ/mol. The barrier can be readily overcome by the thermal motion. So, our results are quite consistent with the experiments for urea and glycerol, for which the corresponding data are available. The results obtained by the 3D-RISM theory show striking differences from those obtained by the MD simulations, especially in the case of GlpF. Possible causes of the difference in the results between the two methods are discussed. [ABSTRACT FROM AUTHOR]

Selective ion binding by human lysozyme and its mutants is probed with the three-dimensional interaction site model theory which is the statistical mechanical integral equation theory. Preliminary and partial results of the study have been already published (Yoshida, N. et al. J. Am. Chem. Soc. 2006, 128, 12042−12043). The calculation was carried out for aqueous solutions of three different electrolytes, CaCl2, NaCl, and KCl, and for four different mutants of the human lysozyme:  wild type, Q86D, A92D, and Q86D/A92D, which have been studied experimentally. The discussion of this article focuses on the cleft that consists of amino acid residues from Q86 to A92. For the wild type of protein in the aqueous solutions of all the electrolytes studied, there are no distributions observed for the ions inside the cleft. The Q86D mutant shows essentially the same behavior with that of the wild type. The A92D mutant shows strong binding ability to Nain the recognition site, which is in accord with the experimental results. There are two isomers of the Q86D/A92D mutant, e.g., apo-Q86D/A92D and holo-Q86D/A92D. Although both isomers exhibit the binding ability to the Naand Ca2ions, the holo isomer shows much greater affinity compared with the apo isomer. Regarding the selective ion binding of the holo-Q86D/A92D mutant, it shows greater affinity to Ca2than to Na, which is also consistent with the experimental observation. [ABSTRACT FROM AUTHOR]