Predicting stability constants of transition metals; Y3+, La3+, and UO2 2+ with organic ligands using the 3D-QSPR methodology

Abstract

Stability constants prediction plays a critical role in the identification and optimization of ligand design for selective complexation of metal ions in solution. Thus, it is important to assess the potential of metal-binding ligand organic in the complex formation process. However, quantitative structure-activity/property relationships (QSAR/QSPR) provide a time-and cost-efficient approach to predict the stability constants of compounds. To this end, we applied a free alignment three-dimensional QSPR technique by generating GRid-INdependent Descriptors (GRINDs) to rationalize the underlying factors effecting on stability constants of transition metals; 105 (Y3+), 186 (La3+), and 66 (UO2 2+) with diverse organic ligands in aqueous solutions at 298 K and an ionic strength of 0.1 M. Kennard– Stone algorithm was employed to split data set to a training set of 75% molecules and a test set of 25% molecules. Fractional factorial design (FFD) and genetic algorithm (GA) applied to derive the most relevant and optimal 3 D molecular descriptors. The selected descriptors using various feature selection were correlated with stability constants by partial least squares (PLS). GA-PLS models were statistically validated ( = 0.96, q2 = 0.82 and R2 pred=0.81 for Y3+; = 0.90, q2 = 0.73 and R2 pred=0.82 for La3+ and = 0.95, q2 = 0.81 and R2 pred=0.88 for UO2 2+), and from the information derived from the graphical results confirmed that hydrogen bonding properties, shape, size, and steric effects are the main parameters influencing stability constant of metal complexation. The provided information in this research can predict the stability constant of the new organic ligand with the transition metals without experimental processes.