Atomic radii play important roles in scientific research. The covalent radii of atoms, ionic radii of ions, and van der Waals (VDW) radii of neutral atoms can all be derived from crystal structures. However, the VDW radii of ions are a challenge to determine because the atomic distances in crystal structures were determined by a combination of VDW interactions and electrostatic interactions, making it unclear how to define the VDW sphere of ions in such an environment. In the present study, we found that VDW radii, which were determined based on the 0.0015 au electron density contour through a wavefunction analysis on atoms, have excellent agreement with the VDW radii of noble-gas atoms determined experimentally. Based on this criterion, we calculated the VDW radii for various atomic ions across the periodic table, providing a systematic set of VDW radii of ions. Previously we have shown that the 12-6 Lennard-Jones nonbonded model could not simultaneously reproduce the hydration free energy (HFE) and ion−oxygen distance (IOD) for an atomic ion when its charge is +2 or higher. Because of this, we developed the 12-6-4 model to reproduce both properties at the same time by explicitly considering the ion-induced dipole interactions. However, recent studies showed it was possible to use the 12-6 model to simulate both properties simultaneously when an ion has the R min /2 parameter (i.e., the VDW radius) close to the Shannon ionic radius. In the present study, we show that such a "success" is due to an unphysical overfitting, as the VDW radius of an ion should be significantly larger than its ionic radius. Through molecular dynamics simulations, we show that such overfitting causes significant issues when transferring the parameters from ion−water systems to ion− ligand and metalloprotein systems. In comparison, the 12-6-4 model shows significant improvement in comparison to the overfitted 12-6 model, showing excellent transferability across different systems. In summary, although both the 12-6-4 and 12-6 models could reproduce HFE and IOD for an ion, the 12-6-4 model accomplishes such a task based on the consideration of the physics involved, while the 12-6 model accomplishes this through overfitting, which brings significant transferability issues when simulating other systems. Hence, we strongly recommend the use of the 12-6-4 model (or even more sophisticated models) instead of overfitted 12-6 models when simulating complex systems such as metalloproteins.

Atomic radii play important role in scientific research. The covalent radii of atoms, ionic radii of ions, and van der Waals (VDW) radii of neutral atoms can all be derived from crystal structures. However, the VDW radii of ions are a challenge to determine because the atomic distances in crystal structures were determined by a combination of the VDW interactions plus the electrostatic interactions, making it unclear how to define the VDW sphere of ions in such an environment. In the present study, we found that the VDW radii, which were determined based on the 0.0015 au electron density contour through a wavefunction analysis on atoms, have excellent agreement with the VDW radii of noble gas atoms determined experimentally. Based on this criterion, we calculated the VDW radii for various atomic ions across the periodic table, providing a systematic set of VDW radii of ions. Previously we have shown that the 12-6 Lennard-Jones nonbonded model could not simultaneously reproduce the hydration free energy (HFE) and ion-oxygen distance (IOD) for an atomic ion when its charge is +2 or higher. Because of this, we developed the 12-6-4 model to reproduce both properties at the same time by explicitly considering the ion-induced dipole interactions. However, recent studies showed it was possible to use the 12-6 model to simulate both properties simultaneously when an ion has the Rmin/2 parameter (i.e., the VDW radius) close to the Shannon ionic radius. In the present study, we show that such a “success” is due to an unphysical overfitting, as the VDW radius of an ion should be significantly larger than its ionic radius. Through molecular dynamics simulations, we show that such overfitting causes significant issues when transferring the parameters from ion-water systems to ion-ligand and metalloprotein systems. In comparison, the 12-6-4 model shows significant improvement in comparison to the overfitted 12-6 model, showing excellent transferability across different systems. In summary, although both the 12-6-4 and 12-6 models could reproduce HFE and IOD for an ion, the 12-6-4 model accomplishes such a task based on the consideration of the physics involved, while the 12-6 model accomplishes this through overfitting, which brings significant transferability issues when simulating other systems. Hence, we strongly recommend the use of the 12-6-4 model (or even more sophisticated models) instead of overfitted 12-6 models when simulating complex systems such as metalloproteins.

No abstract

Calcium-binding proteins play critical roles in various biological processes such as signal transduction, cell growth, and transcription factor regulation. Ion binding and target binding of Ca 2+ -binding proteins are highly related. Therefore, understanding the ion binding mechanism will benefit the relevant inhibitor design toward the Ca 2+ -binding proteins. The EF-hand is the typical ion binding motif in Ca 2+ -binding proteins. Previous studies indicate that the ion binding affinity of the EF-hand increases with the peptide length, but this mechanism has not been fully understood. Herein, using molecular dynamics simulations, thermodynamic integration calculations, and molecular mechanics Poisson−Boltzmann surface area analysis, we systematically investigated four Ca 2+ -binding peptides containing the EF-hand loop in site III of rabbit skeletal troponin C. These four peptides have 13, 21, 26, and 34 residues. Our simulations reproduced the observed trend that the ion binding affinity increases with the peptide length. Our results implied that the E-helix motif preceding the EF-hand loop, likely the Phe99 residue in particular, plays a significant role in this regulation. The E-helix has a significant impact on the backbone and side-chain conformations of the Asp103 residue, rigidifying important hydrogen bonds in the EF-hand and decreasing the solvent exposure of the Ca 2+ ion, hence leading to more favorable Ca 2+ binding in longer peptides. The present study provides molecular insights into the ion binding in the EF-hand and establishes an important step toward elucidating the responses of Ca 2+ -binding proteins toward the ion and target availability.

Atomic radii play important role in scientific research. The covalent radii of atoms, ionic radii of ions, and van der Waals (VDW) radii of neutral atoms can all be derived from crystal structures. However, the VDW radii of ions are a challenge to determine because the atomic distances in crystal structures were determined by a combination of the VDW interactions plus the electrostatic interactions, making it unclear how to define the VDW sphere of ions in such an environment. In the present study, we found that the VDW radii, which were determined based on the 0.0015 au electron density contour through a wavefunction analysis on atoms, have excellent agreement with the VDW radii of noble gas atoms determined experimentally. Based on this criterion, we calculated the VDW radii for various atomic ions across the periodic table, providing a systematic set of VDW radii of ions. Previously we have shown that the 12-6 Lennard-Jones nonbonded model could not simultaneously reproduce the hydration free energy (HFE) and ion-oxygen distance (IOD) for an atomic ion when its charge is +2 or higher. Because of this, we developed the 12-6-4 model to reproduce both properties at the same time by explicitly considering the ion-induced dipole interactions. However, recent studies showed it was possible to use the 12-6 model to simulate both properties simultaneously when an ion has the Rmin/2 parameter (i.e., the VDW radius) close to the Shannon ionic radius. In the present study, we show that such a “success” is due to an unphysical overfitting, as the VDW radius of an ion should be significantly larger than its ionic radius. Through molecular dynamics simulations, we show that such overfitting causes significant issues when transferring the parameters from ion-water systems to ion-ligand and metalloprotein systems. In comparison, the 12-6-4 model shows significant improvement in comparison to the overfitted 12-6 model, showing excellent transferability across different systems. In summary, although both the 12-6-4 and 12-6 models could reproduce HFE and IOD for an ion, the 12-6-4 model accomplishes such a task based on the consideration of the physics involved, while the 12-6 model accomplishes this through overfitting, which brings significant transferability issues when simulating other systems. Hence, we strongly recommend the use of the 12-6-4 model (or even more sophisticated models) instead of overfitted 12-6 models when simulating complex systems such as metalloproteins.