Molecular dynamics simulation of the Pb(II) coordination in biological media via cationic dummy atom models

The coordination of Pb(II) in aqueous solutions containing thiols is a pivotal topic to the understanding of the pollutant potential of this cation. Based on its hard/soft borderline nature, Pb(II) forms stable hydrated ions as well as stable complexes with the thiol groups of proteins. In this paper, the modeling of Pb(II) coordination via classical molecular dynamics simulations was investigated to assess the possible use of non-bonded potentials for the description of the metal–ligand interaction. In particular, this study aimed at testing the capability of cationic dummy atom schemes—in which part of the mass and charge of the Pb(II) is fractioned in three or four sites anchored to the metal center—in reproducing the correct coordination geometry and, also, in describing the hard/soft borderline character of this cation. Preliminary DFT calculations were used to design two topological schemes, PB3 and PB4, that were subsequently implemented in the Amber force field and employed in molecular dynamics simulation of either pure water or thiol/thiolate-containing aqueous solutions. The PB3 scheme was then tested to model the binding of Pb(II) to the lead-sensing protein pbrR. The potential use of CDA topological schemes in the modeling of Pb(II) coordination was here critically discussed.

Dissection of the difference between the group I metal ions in inhibiting GSK3β: a computational study

Glycogen synthase kinase 3β (GSK3β) is a serine/threonine kinase that requires two cofactor Mg(2+) ions for catalysis in regulating many important cellular signals. Experimentally, Li(+) is a competitive inhibitor of GSK3β relative to Mg(2+), while this mechanism is not experienced with other group I metal ions. Herein, we use native Mg(2)(2+)-Mg(1)(2+) GSK3β and its Mg(2)(2+)-M(1)(+) (M = Li, Na, K, and Rb) derivatives to investigate the effect of metal ion substitution on the mechanism of inhibition through two-layer ONIOM-based quantum mechanics/molecular mechanics (QM/MM) calculations and molecular dynamics (MD) simulations. The results of ONIOM calculations elucidate that the interaction of Na(+), K(+), and Rb(+) with ATP is weaker compared to that of Mg(2+) and Li(+) with ATP, and the critical triphosphate moiety of ATP undergoes a large conformational change in the Na(+), K(+), and Rb(+) substituted systems. As a result, the three metal ions (Na(+), K(+), and Rb(+)) are not stable and depart from the active site, while Mg(2+) and Li(+) can stabilize in the active site, evident in MD simulations. Comparisons of Mg(2)(2+)-Mg(1)(2+) and Mg(2)(2+)-Li(1)(+) systems reveal that the inline phosphor-transfer of ATP and the two conserved hydrogen bonds between Lys85 and ATP, together with the electrostatic potential at the Li(1)(+) site, are disrupted in the Mg(2)(2+)-Li(1)(+) system. These computational results highlight the possible mechanism why Li(+) inhibits GSK3β.

Validation of the Existence of Tetrameric Species of Potassium Trimethylsilanolate in the Gas Phase with a Theoretical Cluster Model: Role of the Counterion as Charge Localizer in the Structure

We have investigated, theoretically, the structural properties of potassium trimethylsilanolate in the gas phase at a B3LYP/6-31+G* level. For this purpose, a simplified ionic cluster model based on potassium trimethylsilanolate tetramers, proposed in the literature as the structural units of this compound in the solid state, was developed. Furthermore, we compared the validity of the model with two simpler ones: a monomer of potassium trimethylsilanolate and a trimethylsilanolate anion in the gas phase. The developed ionic cluster model was found to be best in reproducing the experimental structure of potassium trimethylsilanolate, supporting, at the same time, the existence of such tetrameric species (previously identified experimentally from mass spectrometry data by Weiss et al.) in the gas phase. Finally, NBO calculations highlighted the important role of the potassium counterion as a charge localizer in the structure of these chemical species.