Experimental, computational and chemometrics studies of BSA-vitamin B6 interaction by UV–Vis, FT-IR, fluorescence spectroscopy, molecular dynamics simulation and hard-soft modeling methods

Abstract

The interaction of pyridoxine (Vitamin B6) with bovine serum albumin (BSA) is investigated under pseudo-physiological conditions by UV–Vis, fluorescence and FTIR spectroscopy. The intrinsic fluorescence of BSA was quenched by VB6, which was rationalized in terms of the static quenching mechanism. According to fluorescence quenching calculations, the bimolecular quenching constant (kq), dynamic quenching (KSV) and static quenching (KLB) at 310K were obtained. The efficiency of energy transfer and the distance between the donor (BSA) and the acceptor (VB6) were calculated by Foster’s non-radiative energy transfer theory and were equal to 41.1% and 2.11nm.The collected UV–Vis and fluorescence spectra were combined into a row-and column-wise augmented matrix and resolved by multivariate curve resolution-alternating least squares (MCR-ALS). MCR-ALS helped to estimate the stoichiometry of interactions, concentration profiles and pure spectra for three species (BSA, VB6 and VB6-BSA complex) existed in the interaction procedure. Based on the MCR-ALS results, using mass balance equations, a model was developed and binding constant of complex was calculated using non-linear least squares curve fitting. FT-IR spectra showed that the conformation of proteins was altered in presence of VB6. Finally, the combined docking and molecular dynamics (MD) simulations were used to estimate the binding affinity of VB6 to BSA. Five-nanosecond MD simulations were performed on bovine serum albumin (BSA) to study the conformational features of its ligand binding site. From MD results, eleven BSA snapshots were extracted, at every 0.5ns, to explore the binding affinity (GOLD score) of VB6 using a docking procedure. MD simulations indicated that there is a considerable flexibility in the structure of protein that affected ligand recognition. Structural analyses and docking simulations indicated that VB6 binds to site I and GOLD score values depend on the conformations of both BSA and ligand. Molecular modeling results showed that VB6–BSA complex formed not only on the basis of electrostatic forces, but also on the basis of π–π staking and hydrogen bond. There was an excellent agreement between the experimental and computational results. The results presented in this paper, will offer a reference for detailed and systematic studies on the biological effects and action mechanism of small molecules with proteins.