Effects of Electronegativity and Hydration Energy on the Selective Adsorption of Heavy Metal Ions by Synthetic NaX Zeolite

The adsorption capacity of synthetic NaX zeolite for Pb²⁺, Cd²⁺, Cu²⁺ and Zn²⁺ in single and multi-component systems were investigated. The effects of electronegativity and hydration energy on the selective adsorption, as well as potential selective adsorption mechanism of the NaX zeolite for Pb²⁺, Cd²⁺, Cu²⁺ and Zn²⁺ were also discussed. The maximum adsorption capacity order of the heavy metals in the single system was Pb²⁺ > Cd²⁺ > Cu²⁺ > Zn²⁺, and this could be related to their hydration energy and electronegativity. The values of the separation factors (α) and affinity constant (K_EL) in different binary systems indicated that Pb²⁺ was preferentially adsorbed, and Zn²⁺ presented the lowest affinity for NaX zeolite. The selective adsorption capacities of the metals were in the order, Pb²⁺ > Cd²⁺ ≈ Cu²⁺ > Zn²⁺. The trend for the selective adsorption of NaX zeolite in ternary and quaternary systems was consistent with that in the binary systems. Pb²⁺ and Cu²⁺ reduced the stability of the Si-O-Al bonds and the double six-membered rings in the NaX framework, due to the high electronegativity of Pb²⁺ and Cu²⁺ than that of Al³⁺. The selective adsorption mechanism of NaX zeolite for the high electronegative metal ions could mainly result from the negatively charged O in the Si-O-Al structure of the NaX zeolite, hence heavy metal ions with high electronegativity display a strong affinity for the electron cloud of the oxygen atoms in the Si-O-Al. This study could evaluate the application and efficiency of zeolite in separating and recovering certain metal ions from industrial wastewater.

Thermodynamic Integration in 3n Dimensions without Biases or Alchemy for Protein Interactions

Thermodynamic integration (TI), a powerful formalism for computing Gibbs free energy, has been implemented for many biophysical processes with alchemical schemes that require delicate human efforts to choose/design biasing potentials for sampling the desired biophysical events and to remove their artifactitious consequences afterwards. Theoretically, an alchemical scheme is exact but practically, an unsophisticated implementation of this exact formula can cause error amplifications. Small relative errors in the input parameters can be amplified many times in their propagation into the computed free energy [due to subtraction of similar numbers such as (105 ± 5)‒(100 ± 5) = 5 ± 7]. In this paper, we present an unsophisticated implementation of TI in 3n dimensions (3nD) (n=1,2,3…) for the potential of mean force along a 3nD path connecting one state in the bound state ensemble to one state in the unbound state ensemble. Fluctuations in these 3nD are integrated in the bound and unbound state ensembles but not along the 3nD path. Using TI3nD, we computed the standard binding free energies of three protein complexes: trometamol in Salmonella effector SpvD (n=1), biotin-avidin (n=2), and Colicin E9 endonuclease with cognate immunity protein Im9 (n=3). We employed three different protocols in three independent computations of E9-Im9 to show TI3nD's robustness. We also computed the hydration energies of ten biologically relevant compounds (n=1 for water, acetamide, urea, glycerol, trometamol, ammonium and n=2 for erythritol, 1,3-propanediol, xylitol, biotin). Each of the 15 computations is accomplishable within one (for hydration) to ten (for E9-Im9) days on an inexpensive GPU workstation. The computed results all agree with the available experimental data.

Small molecule hydration energy and entropy from 3D-RISM

Implicit solvent models offer an attractive way to estimate the effects of a solvent environment on the properties of small or large solutes without the complications of explicit simulations. One common test of accuracy is to compute the free energy of transfer from gas to liquid for a variety of small molecules, since many of these values have been measured. Studies of the temperature dependence of these values (i.e. solvation enthalpies and entropies) can provide additional insights into the performance of implicit solvent models. Here, we show how to compute temperature derivatives of hydration free energies for the 3D-RISM integral equation approach. We have computed hydration free energies of 1123 small drug-like molecules (both neutral and charged). Temperature derivatives were also used to calculate hydration energies and entropies of 74 of these molecules (both neutral and charged) for which experimental data is available. While direct results have rather poor agreement with experiment, we have found that several previously proposed linear hydration free energy correction schemes give good agreement with experiment. These corrections also provide good agreement for hydration energies and entropies though simple extensions are required in some cases.

Predicting hydration energies for multivalent ions

We have predicted the free energy of hydration for 40 monovalent and multivalent cations and anions using density functional theory and the implicit solvent model COnductor like Screening MOdel for Real Solvents (COSMO-RS) at the Becke-Perdew (BP)/Triple zeta valence with polarization functions (TZVP) level. Agreement with experimental data for monovalent and divalent ions is good and shows no significant systematic errors. Predictions are noticeably better than with standard COSMO. The agreement with experimental data for trivalent and tetravalent ions is slightly worse and shows systematic errors. Our results indicate that quantum chemical calculations combined with COSMO-RS solvent treatment is a reliable method for treating multivalent ions in solution, provided one hydration shell of explicit water molecules is included for metal cations. The accuracy is not high enough to allow absolute predictions of hydration energies but could be used to investigate trends for several ions, thanks to the low computational cost, in particular for ligand exchange reactions.

Development of polarizable models for molecular mechanical calculations. 4. van der Waals parametrization

In the previous publications of this series, we presented a set of Thole induced dipole interaction models using four types of screening functions. In this work, we document our effort to refine the van der Waals parameters for the Thole polarizable models. Following the philosophy of AMBER force field development, the van der Waals (vdW) parameters were tuned for the Thole model with linear screening function to reproduce both the ab initio interaction energies and the experimental densities of pure liquids. An in-house genetic algorithm was applied to maximize the fitness of "chromosomes" which is a function of the root-mean-square errors (RMSE) of interaction energy and liquid density. To efficiently explore the vdW parameter space, a novel approach was developed to estimate the liquid densities for a given vdW parameter set using the mean residue-residue interaction energies through interpolation/extrapolation. This approach allowed the costly molecular dynamics simulations be performed at the end of each optimization cycle only and eliminated the simulations during the cycle. Test results show notable improvements over the original AMBER FF99 vdW parameter set, as indicated by the reduction in errors of the calculated pure liquid densities (d), heats of vaporization (H(vap)), and hydration energies. The average percent error (APE) of the densities of 59 pure liquids was reduced from 5.33 to 2.97%; the RMSE of H(vap) was reduced from 1.98 to 1.38 kcal/mol; the RMSE of solvation free energies of 15 compounds was reduced from 1.56 to 1.38 kcal/mol. For the interaction energies of 1639 dimers, the overall performance of the optimized vdW set is slightly better than the original FF99 vdW set (RMSE of 1.56 versus 1.63 kcal/mol). The optimized vdW parameter set was also evaluated for the exponential screening function used in the Amoeba force field to assess its applicability for different types of screening functions. Encouragingly, comparable performance was observed when the optimized vdW set was combined with the Thole Amoeba-like polarizable model, particularly for the interaction energy and liquid density calculations. Thus, the optimized vdW set is applicable to both types of Thole models with either linear or Amoeba-like screening functions.

Cystic fibrosis transmembrane conductance regulator. Physical basis for lyotropic anion selectivity patterns

The cystic fibrosis transmembrane conductance regulator (CFTR) Cl channel exhibits lyotropic anion selectivity. Anions that are more readily dehydrated than Cl exhibit permeability ratios (P(S)/P(Cl)) greater than unity and also bind more tightly in the channel. We compared the selectivity of CFTR to that of a synthetic anion-selective membrane [poly(vinyl chloride)-tridodecylmethylammonium chloride; PVC-TDMAC] for which the nature of the physical process that governs the anion-selective response is more readily apparent. The permeability and binding selectivity patterns of CFTR differed only by a multiplicative constant from that of the PVC-TDMAC membrane; and a continuum electrostatic model suggested that both patterns could be understood in terms of the differences in the relative stabilization of anions by water and the polarizable interior of the channel or synthetic membrane. The calculated energies of anion-channel interaction, derived from measurements of either permeability or binding, varied as a linear function of inverse ionic radius (1/r), as expected from a Born-type model of ion charging in a medium characterized by an effective dielectric constant of 19. The model predicts that large anions, like SCN, although they experience weaker interactions (relative to Cl) with water and also with the channel, are more permeant than Cl because anion-water energy is a steeper function of 1/r than is the anion-channel energy. These large anions also bind more tightly for the same reason: the reduced energy of hydration allows the net transfer energy (the well depth) to be more negative. This simple selectivity mechanism that governs permeability and binding acts to optimize the function of CFTR as a Cl filter. Anions that are smaller (more difficult to dehydrate) than Cl are energetically retarded from entering the channel, while the larger (more readily dehydrated) anions are retarded in their passage by "sticking" within the channel.

Relaxation study of complex formation between monovalent cations and cyclic polyethers

Equilibrium quotients for the 1:1 complexes of dibenzo-30-crown-10 with Na(+), K(+), Rb(+), Cs(+), NH(4) (+), and T1(+) in methanol solution have been determined spectrophotometrically. The observed selectivity pattern is discussed in terms of the hydration energy of the cations and the cation-ligand binding energy. The rate of complex formation is diffusion-controlled. Relaxation amplitude data are consistent with a mechanism that involves conformational change of the crown compound.