Molecular dynamics study of conformational changes in human serum albumin by binding of fatty acids

Human serum albumin (HSA) binds with fatty acids under normal physiologic conditions. To date, there is little published information on the tertiary structure of HSA-fatty acid complex in aqueous solution. In the present study, we used molecular dynamics (MD) simulations to elucidate possible structural changes of HSA brought about by the binding of fatty acids. Both unliganded HSA and HSA-fatty acid complex models for MD calculations were constructed based on the X-ray crystal structures. Five myristates (MYRs) were bound in the HSA-fatty acid complex model. In the present MD study, the motion of domains I and III caused by the binding of MYR molecules increased the radius of gyration of HSA. Root-mean-square fluctuations from the MD simulations revealed that the atomic fluctuations of the specific amino acids at drug-binding site I that can regulate the drug-binding affinity were increased by the binding of MYR molecules. Primary internal motions, characterized by the first three principal components, were observed mainly at domains I and III in the principal component analysis for trajectory data. The directional motion projected on the first principal component of unliganded HSA was conserved in HSA-MYR complex as the third principal directional motion with higher frequency. However, the third principal directional motion in unliganded HSA turned into the first principal directional motion with lower frequency in the HSA-MYR complex. Thus, the present MD study provides insights into the possible conformational changes of HSA caused by the binding of fatty acids.     

Exploration of human serum albumin binding sites by docking and molecular dynamics flexible ligand–protein interactions

Five‐nanosecond molecular dynamics (MD) simulations were performed on human serum albumin (HSA) to study the conformational features of its primary ligand binding sites (I and II). Additionally, 11 HSA snapshots were extracted every 0.5 ns to explore the binding affinity (K_d) of 94 known HSA binding drugs using a blind docking procedure. MD simulations indicate that there is considerable flexibility for the protein, including the known sites I and II. Movements at HSA sites I and II were evidenced by structural analyses and docking simulations. The latter enabled the study and analysis of the HSA–ligand interactions of warfarin and ketoprofen (ligands binding to sites I and II, respectively) in greater detail. Our results indicate that the free energy values by docking (K_d observed) depend upon the conformations of both HSA and the ligand. The 94 HSA–ligand binding K_d values, obtained by the docking procedure, were subjected to a quantitative structure‐activity relationship (QSAR) study by multiple regression analysis. The best correlation between the observed and QSAR theoretical (K_d predicted) data was displayed at 2.5 ns. This study provides evidence that HSA binding sites I and II interact specifically with a variety of compounds through conformational adjustments of the protein structure in conjunction with ligand conformational adaptation to these sites. These results serve to explain the high ligand‐promiscuity of HSA. © 2009 Wiley Periodicals, Inc. Biopolymers 93: 161–170, 2010. This article was originally published online as an accepted preprint. The “Published Online” date corresponds to the preprint version. You can request a copy of the preprint by emailing the Biopolymers editorial office at biopolymers@wiley.com     

Unraveling the binding mechanism of asiatic acid with human serum albumin and its biological implications

Asiatic acid (AsA), a naturally occurring pentacyclictriterpenoid found in Centella asiatica, plays a major role in neuroprotection, anticancer, antioxidant, and hepatoprotective activities. Human serum albumin (HSA), a blood plasma protein, participates in the regulation of plasma osmotic pressure and transports endogenous and exogenous substances. The study undertaken to analyze the drug-binding mechanisms of HSA is crucial in understanding the bioavailability of drugs. In this study, we analyzed the cytotoxic activity of AsA on HepG2 (human hepatocellular carcinoma) cell lines and its binding, conformational, docking, molecular simulation studies with HSA under physiological pH 7.2. These studies revealed a clear decrease in the viability of HepG2 cells upon exposure to AsA in a dose-dependent manner with an IC50 of 45?μM. Further studies showed the quenching of intrinsic fluorescence of HSA by AsA with a binding constant of K_AsA?=?3.86?±?0.01?×?10⁴?M^?1, which corresponds to the free energy of (ΔG) ?6.3?kcal?M^?1 at 25?°C. Circular dichroism (CD) studies revealed that there is a clear decrease in the α-helical content from 57.50?±?2.4 to 50%?±?2.3 and an increase in the β-turns from 25?±?0.65 to 29%?±?0.91 and random coils from 17.5%?±?0.95 to 21%?±?1.2, suggesting partial unfolding of HSA. Autodock studies revealed that the AsA is bound to the subdomain IIA with hydrophobic and hydrophilic interactions. From molecular dynamics, simulation data (RMSD, Rg and RMSF) emphasized the local conformational changes and rigidity of the residues of both HSA and HSA–AsA complexes.     

Functional roles of magnesium binding to extracellular signal-regulated kinase 2 explored by molecular dynamics simulations and principal component analysis

Molecular dynamics (MD) simulations coupled with principal component (PC) analysis were carried out to study functional roles of Mg²⁺ binding to extracellular signal-regulated kinase 2 (ERK2). The results suggest that Mg²⁺ binding heavily decreases eigenvalue of the first principal component and totally inhibits motion strength of ERK2, which favors stabilization of ERK2 structure. Binding free energy predictions indicate that Mg²⁺ binding produces an important effect on binding ability of adenosine triphosphate (ATP) to ERK2 and strengthens the ATP binding. The calculations of residue-based free energy decomposition show that lack of Mg²⁺ weakens interactions between the hydrophobic rings of ATP and five residues I29, V37, A50, L105, and L154. Hydrogen bond analyses also prove that Mg²⁺ binding increases occupancies of hydrogen bonds formed between ATP and residues K52, Q103, D104, and M106. We expect that this study can provide a significant theoretical hint for designs of anticancer drugs targeting ERK2.     

Spectroscopic and molecular modeling studies on the binding of the flavonoid luteolin and human serum albumin

Luteolin (LUT) is a polyphenolic compound, found in a variety of fruits, vegetables, and seeds, which has a variety of pharmacological properties. In the present contribution, binding of LUT to human serum albumin (HSA), the most abundant carrier protein in the blood, was investigated with the aim of describing the binding mode and parameters of the interaction. The application of circular dichroism, UV‐Vis absorption, fluorescence, Raman and surface‐enhanced Raman scattering spectroscopy combined with molecular modeling afforded a clear picture of the association mode of LUT to HSA. Specific interactions with protein amino acids were evidenced. LUT was found to be associated in subdomain IIA where an interaction with Trp‐214 is established. Hydrophobic and electrostatic interactions are the major acting forces in the binding of LUT to HSA. The HSA conformations were slightly altered by the drug complexation with reduction of α‐helix and increase of β‐turns structures, suggesting a partial protein unfolding. Also the configuration of at least two disulfide bridges were altered. Furthermore, the study of molecular modeling afforded the binding geometry. © 2009 Wiley Periodicals, Inc. Biopolymers 91: 917–927, 2009. This article was originally published online as an accepted preprint. The “Published Online” date corresponds to the preprint version. You can request a copy of the preprint by emailing the Biopolymers editorial office at biopolymers@wiley.com     

Unveiling the binding mode of perfluorooctanoic acid to human serum albumin

Perfluorooctanoic acid (PFOA) is a toxic compound that is absorbed and distributed throughout the body by noncovalent binding to serum proteins such as human serum albumin (hSA). Though the interaction between PFOA and hSA has been already assessed using various analytical techniques, a high resolution and detailed analysis of the binding mode is still lacking. We report here the crystal structure of hSA in complex with PFOA and a medium‐chain saturated fatty acid (FA). A total of eight distinct binding sites, four occupied by PFOAs and four by FAs, have been identified. In solution binding studies confirmed the 4:1 PFOA‐hSA stoichiometry and revealed the presence of one high and three low affinity binding sites. Competition experiments with known hSA‐binding drugs allowed locating the high affinity binding site in sub‐domain IIIA. The elucidation of the molecular basis of the interaction between PFOA and hSA might provide not only a better assessment of the absorption and elimination mechanisms of these compounds in vivo but also have implications for the development of novel molecular receptors for diagnostic and biotechnological applications.     

Conformational changes and allosteric communications in human serum albumin due to ligand binding

It is well recognized that knowledge of structure alone is not sufficient to understand the fundamental mechanism of biomolecular recognition. Information of dynamics is necessary to describe motions involving relevant conformational states of functional importance. We carried out principal component analysis (PCA) of structural ensemble, derived from 84 crystal structures of human serum albumin (HSA) with different ligands and/or different conditions, to identify the functionally important collective motions, and compared with the motions along the low-frequency modes obtained from normal mode analysis of the elastic network model (ENM) of unliganded HSA. Significant overlap is observed in the collective motions derived from PCA and ENM. PCA and ENM analysis revealed that ligand selects the most favored conformation from accessible equilibrium structures of unliganded HSA. Further, we analyzed dynamic network obtained from molecular dynamics simulations of unliganded HSA and fatty acids- bound HSA. Our results show that fatty acids-bound HSA has more robust community network with several routes to communicate among different parts of the protein. Critical nodes (residues) identified from dynamic network analysis are in good agreement with allosteric residues obtained from sequence-based statistical coupling analysis method. This work underscores the importance of intrinsic structural dynamics of proteins in ligand recognition and can be utilized for the development of novel drugs with optimum activity.     

Probing the binding interaction of AKR with human serum albumin by multiple fluorescence spectroscopy and molecular modeling

Human serum albumin (HSA) is the major transport protein affording endogenous and exogenous substances in plasma. It can affect the behavior and efficacy of chemicals in vivo through the binding interaction. AKR (3-O-α-l-arabinofuranosyl-kaempferol-7-O-α-l-rhamnopyranoside) is a flavonoid diglycoside with modulation of estrogen receptors (ERs). Herein, we investigated the binding interaction between AKR and HSA by multiple fluorescence spectroscopy and molecular modeling. As a result, AKR specifically binds in site I of HSA through hydrogen bonds, van der Waals force, and electrostatic interaction. The formation of AKR–HSA complex in binding process is spontaneously exothermic and leads to the static fluorescence quenching through affecting the microenvironment around the fluorophores. The complex also affects the backbone of HSA and makes AKR access to fluorophores. Molecular modeling gives the visualization of the interaction between AKR and HSA as well as ERs. The affinity of AKR with HSA is higher than the competitive site marker Warfarin. In addition, docking studies reveal the binding interaction of AKR with ERs through hydrogen bonds, van der Waals force, hydrophobic, and electrostatic interactions. And AKR is more favorable to ERβ. These results unravel the binding interaction of AKR with HSA and mechanism as an ERs modulator.     

Fipronil recognition by the FA1 site of human serum albumin

Fipronil is a broad‐spectrum pesticide widely used in agriculture, horticulture, and forestry. Because fipronil can cause a variety of toxic effects in animals and humans, its use is authorized as a pesticide in veterinary medicinal products for pets, but not for the treatment of livestock animals whose products are intended for consumption. Recently, however, the presence of fipronil residues has been detected in the eggs and meat of layer hens from farms located in different European countries. Given the relevance of fipronil toxicity for human health, it is important to gain information concerning its fate in the human body, including its binding mode to human serum albumin (HSA), the most abundant protein in plasma. Here, the inhibition of heme‐Fe(III) binding to the fatty acid site 1 (FA1) of HSA by fipronil is reported. Docking simulations support functional data, indicating that the FA1 site is the preferential cleft for fipronil recognition by HSA. The affinity of fipronil for HSA (K_f = 1.9 × 10^?6 M, at pH 7.3, and 20.0°C) may be relevant in vivo. Indeed, HSA could play a pivotal role in fipronil transport and scavenging, thus reducing the pesticide‐free plasmatic levels, with consequent reduced systemic toxicity. In turn, fipronil binding to the FA1 site of HSA could impair the recognition of endogenous and exogenous molecules.     

Understanding the effects of two bound glucose in Sudlow site I on structure and function of human serum albumin: theoretical studies

Human serum albumin (HSA) is the most abundant protein found in blood serum. It carries essential metabolites and many drugs. The glycation of HSA causes abnormal biological effects. Importantly, glycated HSA (GHSA) is of interest as a biomarker for diabetes. Recently, the first HSA structure with bound pyranose (GLC) and open-chain (GLO) glucose at Sudlow site I has been crystallised. We therefore employed molecular dynamics (MD) simulations and ONIOM calculations to study the dynamic nature of two bound glucose in a pre-glycated HSA (pGHSA) and observe how those sugars alter a protein structure comparing to wild type (Apo) and fatty acid-bound HSA (FA). Our analyses show that the overall structural stability of pGHSA is similar to Apo and FA, except Sudlow site II. Having glucose induces large protein flexibility at Sudlow site II. Besides, the presence of glucose causes W214 to reorient resulting in a change in W214 microenvironment. Considering sugars, both sugars are exposed to water, but GLO is more solvent-accessible. ONIOM results show that glucose binding is favoured for HSA (?115.04 kcal/mol) and GLO (?85.10 kcal/mol) is more preferable for Sudlow site I over GLC (?29.94 kcal/mol). GLO can strongly react with K195 and K199, whereas K195 and K199 provide slightly repulsive forces for GLC. This can confirm that an open-chain GLO is more favourable inside a pocket.     

Cantharidin inhibits competitively heme‐Fe(III) binding to the FA1 site of human serum albumin

Cantharidin, a monoterpene isolated from the insect blister beetle, has long been used as a medicinal agent in the traditional Chinese medicine. Cantharidin inhibits a subgroup of serine/threonine phosphatases, thus inducing cell growth inhibition and cytotoxicity. Cantharidin has anticancer activity in vitro, since it is able of inducing p53‐dependent apoptosis and double‐strand breakage of DNA in cancer cells. Although the toxicity of cantharidin to the gastrointestinal and urinary tracts prevents its medical use, it is a promising lead compound for chemical modification to develop new anticancer therapeutics. In fact, cantharidin does not cause myelosuppression and displays anticancer activity against cells with a multidrug resistance phenotype. Here, the competitive inhibitory effect of cantharidin on heme‐Fe(III) binding to the fatty acid site 1 (FA1) of human serum albumin (HSA) is reported. Docking and molecular dynamics simulations support functional data indicating the preferential binding of cantharidin to the FA1 site of HSA. Present results may be relevant in vivo as HSA could transport cantharidin, which in turn could affect heme‐Fe(III) scavenging by HSA.     

Metal-catalyzed oxidation of human serum albumin does not alter the interactive binding to the two principal drug binding sites

It is well known that various physiological factors such as pH, endogenous substances or post-translational modifications can affect the conformational state of human serum albumin (HSA). In a previous study, we reported that both pH- and long chain fatty acid-induced conformational changes can alter the interactive binding of ligands to the two principal binding sites of HSA, namely, site I and site II. In the present study, the effect of metal-catalyzed oxidation (MCO) caused by ascorbate/oxygen/trace metals on HSA structure and the interactive binding between dansyl-L-asparagine (DNSA; a site I ligand) and ibuprofen (a site II ligand) at pH 6.5 was investigated. MCO was accompanied by a time-dependent increase in carbonyl content in HSA, suggesting that the HSA was being oxidized. In addition, The MCO of HSA was accompanied by a change in net charge to a more negative charge and a decrease in thermal stability. SDS-PAGE patterns and α-helical contents of the oxidized HSAs were similar to those of native HSA, indicating that the HSA had not been extensively structurally modified by MCO. MCO also caused a selective decrease in ibuprofen binding. In spite of the changes in the HSA structure and ligand that bind to site II, no change in the interactive binding between DNSA and ibuprofen was observed. These data indicated that amino acid residues in site II are preferentially oxidized by MCO, whereas the spatial relationship between sites I and II (e.g. the distance between sites), the flexibility or space of each binding site are not altered. The present findings provide insights into the structural characteristics of oxidized HSA, and drug binding and drug-drug interactions on oxidized HSA.     

Molecular dynamics and principal components of potassium binding with human telomeric intra-molecular G-quadruplex

Telomere assumes intra-molecular G-quadruplex that is a significant drug target for inhibiting telomerase maintenance of telomeres in cancer. Metal cations have been recognized as playing important roles in stabilizing G-quadruplex, but their binding processes to human telomeric G-quadruplex remain uncharacterized. To investigate the detailed binding procedures, molecular dynamics simulations were conducted on the hybrid [3+ 1] form-one human telomeric intra-molecular G-quadruplex. We show here that the binding of a potassium ion to a G-tetrad core is mediated by two alternative pathways. Principal component analysis illustrated the dominant concerted motions of G-quadruplex occurred at the loop domains. MM-PBSA calculations revealed that binding was energetically favorable and driven by the electrostatic interactions. The lower binding site was found more constructive favorable for binding. Our data provide useful information on a potassium-mediated stable structure of human telomeric intra-molecular G-quadruplex, implicating in ion disorder associated conformational changes and targeted drug design.     

Determination on the binding of thiadiazole derivative to human serum albumin: a spectroscopy and computational approach

4-[3-acetyl-5-(acetylamino)-2,3-dihydro-1,3,4-thiadiazole-2-yl]phenyl benzoate from the family of thiadiazole derivative has been newly synthesized. It has good anticancer activity as well as antibacterial and less toxic in nature, its binding characteristics are therefore of huge interest for understanding pharmacokinetic mechanism of the drug. The binding of thiadiazole derivative to human serum albumin (HSA) has been investigated by studying its quenching mechanism, binding kinetics and the molecular distance, r between the donor (HSA) and acceptor (thiadiazole derivative) was estimated according to Forster’s theory of non-radiative energy transfer. The Gibbs free energy (ΔG), enthalpy (ΔH) and entropy (ΔS) changes of temperature-dependent K_b was calculated, which explains that the reaction is spontaneous and exothermic. The microenvironment of HSA have also been studied using synchronous fluorescence spectroscopy, and the feature of thiadiazole derivative-induced structural changes of HSA have been carried using Fourier transform infrared spectroscopy and the Molecular modelling simulations explore the hydrophobic and hydrogen bonding interactions.     

Probing the combination of erlotinib hydrochloride,an anticancer drug,and human serum albumin: Spectroscopic,molecular docking,and molecular dynamic analyses

Human serum albumin (HSA) is a globular and monomeric protein in plasma that transports many drugs and compounds. Binding of some drugs to HSA can lead to changes in its stability and biological function. We investigated the binding interactions between erlotinib hydrochloride (Erlo) and HSA. Erlo is used to treat lung, pancreatic, and some other cancers. Experimental data showed that the fluorescence emission of the protein was quenched by Erlo using a static quenching mechanism. The calculation of the binding constant, K_b (1.57 × 10⁵ M⁻¹ at 300 K), confirmed the existence of a moderate binding interaction between Erlo and HSA. The interaction was enthalpy driven, spontaneous, and exothermic. The calculated thermodynamic parameters in agreement with simulation and molecular docking data showed that van der Waals and hydrogen bond forces played an important role in the interaction process. Molecular docking results indicated that Erlo has more affinity to bind to subdomain IIA (site I) of HSA. Molecular dynamics simulation analysis showed that the protein is stable in the presence of Erlo under simulation conditions.     

Study on the molecular recognition action of lamivudine by human serum albumin

In this work, the interaction of an anti‐HIV drug lamivudine and human serum albumin (HSA) was studied by multispectroscopic and molecular modeling methods. The fluorescence emission spectra showed that the fluorescence of HSA was quenched by lamivudine through static mechanism with HSA‐lamivudine complex produced at ground state. According to the binding equilibriums observed at 4 different temperatures, the number of binding site, binding constant, enthalpy change, entropy change, and Gibbs free energy change of the interaction were calculated. The results indicated that there was only 1 main binding site under present concentration condition, and then, the location of this binding site was ascertained by molecular probe experiments using warfarin and ibuprofen as site markers. The interaction was a spontaneous and exothermic process. Hydrogen bonds and van der Waals force were the major power that fixed lamivudine on Sudlow's site I in subdomain IIA of HSA molecule. The distance between donor and acceptor was determined by Förster's nonradiative fluorescence resonance energy transfer theory. Circular dichroism spectra exhibited the alteration of HSA's secondary structures. Molecular modeling investigation revealed the structure of HSA‐lamivudine complex, including the conformation of lamivudine in binding site, the amino residues close to lamivudine, and the interaction forces between receptor and ligand. The study may be beneficial to therapists in understanding the distribution of lamivudine in vivo and explaining its drug‐resistant mechanism in clinical diagnosis.     

Study on the multiple sites binding of human serum albumin and porphyrin by affinity capillary electrophoresis

Equations to describe the two sites binding between proteins and ligands were deduced. According to these equations, not only the binding constants, but also the mole fraction of proteins in different forms could be obtained. Using the published data on the interaction between human serum albumin (HSA) and three kinds of porphyrin (coproporphyrin (CP), uroporphyrin I (UP) and protoporphyrin (PP)), a further study on their binding was carried out. It was concluded that there may exist two binding sites with the binding constants at the first site, proved to be the preferential one, being 6.50 x l0(5), 1.94 x 10(6) and 8.94 x 10(5), respectively. In addition, it was also demonstrated that the two binding sites of HSA with CP and UP might be of different kinds, though those of HSA and PP were of the same kind but at different positions.     

Interaction of pirprofen enantiomers with human serum albumin.

The interaction of pirprofen enantiomers with human serum albumin (HSA) was investigated by means of high-performance liquid chromatography (HPLC), circular dichroism (CD), and 1H NMR spectroscopy. HPLC experiments indicated that both pirprofen enantiomers were bound to one class of high-affinity binding sites (n(+) = 1.91 +/- 0.13, K(+) = (4.09 +/- 0.64) x 10(5) M-1, n(-) = 2.07 +/- 0.13, K(-) = (6.56 +/- 1.35) x 10(5) M-1) together with nonspecific binding (n'K'(+) = (1.51 +/- 0.21) x 10(4) M-1, n'K'(-) = (0.88 +/- 0.13) x 10(-4) M-1). Slight stereoselectivity in specific binding was demonstrated by the difference in product n(+)K(+) = (0.77 +/- 0.08) x 10(6) M-1 vs. n(-)K(-) = (1.30 +/- 0.21) x 10(6) M-1, i.e., the ratio n(-)K(-)/n(+)K(+) = 1.7. CD measurements showed changes in the binding sites located on the aromatic amino acid side chains (a small positive band at 315 nm and a pronounced negative extrinsic Cotton effect in the region 250-280 nm). The protein remains, however, in its predominantly alpha-helical conformation. The 1H NMR difference spectra confirmed that both pirprofen enantiomers interacted with HSA specifically, most probably with site II on the albumin molecule.     

Synthesis,characterization, and binding assessment with human serum albumin of three bipyridine lanthanide(III) complexes

In this work, the terbium(III), dysprosium(III), and ytterbium(III) complexes containing 2, 2′-bipyridine (bpy) ligand have been synthesized and characterized using CHN elemental analysis, FT-IR, UV–Vis and ¹H-NMR techniques and their binding behavior with human serum albumin (HSA) was studied by UV–Vis, fluorescence and molecular docking examinations. The experimental data indicated that all three lanthanide complexes have high binding affinity to HSA with effective quenching of HSA fluorescence via static mechanism. The binding parameters, the type of interaction, the value of resonance energy transfer, and the binding distance between complexes and HSA were estimated from the analysis of fluorescence measurements and Förster theory. The thermodynamic parameters suggested that van der Waals interactions and hydrogen bonds play an important role in the binding mechanism. While, the energy transfer from HSA molecules to all these complexes occurs with high probability, the order of binding constants (BpyTb > BpyDy > BpyYb) represents the importance of radius of Ln³⁺ ion in the complex-HSA interaction. The results of molecular docking calculation and competitive experiments assessed site 3 of HSA, located in subdomain IB, as the most probable binding site for these ligands and also indicated the microenvironment residues around the bound mentioned complexes. The computational results kept in good agreement with experimental data.     

Experimental and computational studies on the binding of diazinon to human serum albumin

In the present research, the binding properties of diazinon (DZN), as an organophosphorus herbicide, to human serum albumin (HSA) were investigated using combination of spectroscopic, electrochemistry, and molecular modeling techniques. Changes in the UV–Vis and FT-IR spectra were observed upon ligand binding along with a significant degree of tryptophan fluorescence quenching on complex formation. The obtained results from spectroscopic and electrochemistry experiments along with the computational studies suggest that DZN binds to residues located in subdomains IIA of HSA with binding constant about 1410.9 M^?1 at 300 K. From the thermodynamic parameters calculated according to the van’t Hoff equation, the enthalpy change ΔH° and entropy change ΔS° were found to be ?16.695 and 0.116 KJ/mol K, respectively. The primary binding pattern is determined by hydrophobic interaction and hydrogen binding occurring in so-called site I of HSA. DZN could slightly alter the secondary structure of HSA. All of experimental results are supported by computational techniques such as docking and molecular dynamics simulation using a HSA crystal model.