Probing the binding mechanism of capecitabine to human serum albumin using spectrometric methods,molecular modeling,and chemometrics approach

Capecitabine as a prodrug of 5-Fluorouracil plays an important role in the treatment of breast and gastrointestinal cancers. Herein, in view of the importance of this drug in chemotherapy, interaction mechanism between Capecitabine (CAP) and human serum albumin (HSA) as a major transport protein in the blood circulatory system has been investigated by using a combination of spectroscopic and molecular modeling methods. The fluorescence spectroscopic results revealed that capecitabine could effectively quench the intrinsic fluorescence of HSA through a static quenching mechanism. Evaluation of the thermodynamic parameters suggested that the binding process was spontaneous while hydrogen bonds and van der Waals forces played a major role in this interaction. The value of the binding constant (K_b = 1.820 × 10⁴) suggested a moderate binding affinity between CAP and HSA which implies its easy diffusion from the circulatory system to the target tissue. The efficiency of energy transfer and the binding distance between the donor (HSA) and acceptor (CAP) were determined according to forster theory of nonradiation energy transfer as 0.410 and 4.135 nm, respectively. Furthermore, UV–Vis spectroscopic results confirmed that the interaction was occurred between HSA and CAP and caused conformational and micro-environmental changes of HSA during the interaction. Multivariate curve resolution-alternating least square (MCR-ALS) methodology as an efficient chemometric tool was used to separate the overlapped spectra of the species. The MCR-ALS result was exploited to estimate the stoichiometry of interaction and to provide concentration and structural information about HSA-CAP interactions. Molecular docking studies suggested that CAP binds mainly to the subdomain IIA of HSA, which were compatible with those obtained by experimental data. Finally, molecular dynamics simulation (MD) was performed on the best docked complex by considering the permanence and flexibility of HSA-CAP complex in the binding site. MD result showed that CAP could steadily bind to HSA in the site I based on the formation of hydrogen bond and π-π stacking interaction in addition to hydrophobic force.     

Identification of potent inhibitors against snake venom metalloproteinase (SVMP) using molecular docking and molecular dynamics studies

Snake venom metalloproteinase (SVMP) (Echis coloratus (Carpet viper) is a multifunctional enzyme that is involved in producing several symptoms that follow a snakebite, such as severe local hemorrhage, nervous system effects and tissue necrosis. Because the three-dimensional (3D) structure of SVMP is not known, models were constructed, and the best model was selected based on its stereo-chemical quality. The stability of the modeled protein was analyzed through molecular dynamics (MD) simulation studies. Structure-based virtual screening was performed, and 15 potential molecules with the highest binding energies were selected. Further analysis was carried out with induced fit docking, Prime/MM–GBSA (ΔG_Bind calculations), quantum-polarized ligand docking, and density functional theory calculations. Further, the stability of the lead molecules in the SVMP-active site was examined using MD simulation. The results showed that the selected lead molecules were highly stable in the active site of SVMP. Hence, these molecules could potentially be selective inhibitors of SVMP. These lead molecules can be experimentally validated, and their backbone structural scaffold could serve as building blocks in designing drug-like molecules for snake antivenom.     

An explorative study on Staphylococcus aureus MurE inhibitor: induced fit docking,binding free energy calculation,and molecular dynamics

Staphylococcus aureus MurE enzyme catalyzes the addition of l-lysine as third residue of the peptidoglycan peptide moiety. Due to the high substrate specificity and its ubiquitous nature among bacteria, MurE enzyme is considered as one of the potential target for the development of new therapeutic agents. In the present work, induced fit docking (IFD), binding free energy calculation, and molecular dynamics (MD) simulation were carried out to elucidate the inhibition potential of 2-thioxothiazolidin-4-one based inhibitor 1 against S. aureus MurE enzyme. The inhibitor 1 formed majority of hydrogen bonds with the central domain residues Asn151, Thr152, Ser180, Arg187, and Lys219. Binding free-energy calculation by MM-GBSA approach showed that van der Waals (ΔG_vdW, ?57.30?kcal/mol) and electrostatic solvation (ΔG_solv, ?36.86?kcal/mol) energy terms are major contributors for the inhibitor binding. Further, 30-ns MD simulation was performed to validate the stability of ligand–protein complex and also to get structural insight into mode of binding. Based on the IFD and MD simulation results, we designed four new compounds D1–D4 with promising binding affinity for the S. aureus MurE enzyme. The designed compounds were subjected to the extra-precision docking and binding free energy was calculated for complexes. Further, a 30-ns MD simulation was performed for D1/4C13 complex.     

Intestinal transport of TRH analogs through PepT1: the role of in silico and in vitro modeling

The present study involves molecular docking, molecular dynamics (MD) simulation studies, and Caco‐2 cell monolayer permeability assay to investigate the effect of structural modifications on PepT1‐mediated transport of thyrotropin releasing hormone (TRH) analogs. Molecular docking of four TRH analogs was performed using a homology model of human PepT1 followed by subsequent MD simulation studies. Caco‐2 cell monolayer permeability studies of four TRH analogs were performed at apical to basolateral and basolateral to apical directions. Inhibition experiments were carried out using Gly‐Sar, a typical PepT1 substrate, to confirm the PepT1‐mediated transport mechanism of TRH analogs. P_app of the four analogs follows the order: NP‐1894 < NP‐2378 < NP‐1896 < NP‐1895. Higher absorptive transport was observed in the case of TRH analogs, indicating the possibility of a carrier‐mediated transport mechanism. Further, the significant inhibition of the uptake of Gly‐Sar by TRH analogs confirmed the PepT1‐mediated transport mechanism. Glide docking scores of all the four analogues were in good agreement with their transport rates, suggesting the role of substrate binding affinity in the PepT1‐mediated transport of TRH analogs. MD simulation studies revealed that the polar interactions with amino acid residues present in the active site are primarily responsible for substrate binding, and a downward trend was observed with the increase in bulkiness at the N‐histidyl moiety of TRH analogs. Copyright © 2014 John Wiley & Sons, Ltd.     

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

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 (k_q), dynamic quenching (K_SV) and static quenching (K_LB) at 310 K 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.11 nm.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.5 ns, 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.     

Structure and dynamics of the full-length M1 muscarinic acetylcholine receptor studied by molecular dynamics simulations

The three-dimensional structure of full-length structure of the M₁ muscarinic receptor was obtained through the fragmental homology modeling procedure. A 10-ns molecular dynamics (MD) simulation of the protein imbedded in a lipid slab and surrounded by water molecules was further used to relax the model. It was found that the homology model corresponded to the conformation in the ground state, since no significant motions of the backbone of transmembrane domains were observed. Furthermore, the reliability of the model was validated by analyzing key inter-helical contacts, sidechain-sidechain interactions, the formation of stable aromatic microdomains (clusters) and the docking of acetylcholine to its binding site. Moreover, a few conserved interactions observed in the X-ray structure of rhodopsin, such as inter-helical sidechain-sidechain hydrogen bonds were accurately reproduced in the MD simulation. The coupling of ACh to its binding site was found to be dominated by π-cation and salt bridge interactions, while its conformational space was restrained through van der Waals and hydrogen bond interactions. In general, such features were in very good agreement with the available experimental as well as with theoretical data. Considering the above, the structural information obtained in this study can be used a starting point to investigate the activation mechanism of the receptor and the ability to develop selective agonists and allosteric modulators which could be used for the treatment of Alzheimer’s disease.     

3D-QSAR,molecular docking and molecular dynamics studies of a series of RORγt inhibitors

The discovery of clinically relevant inhibitors of retinoic acid receptor-related orphan receptor-gamma-t (RORγt) for autoimmune diseases therapy has proven to be a challenging task. In the present work, to find out the structural features required for the inhibitory activity, we show for the first time a three-dimensional quantitative structure–activity relationship (3D-QSAR), molecular docking and molecular dynamics (MD) simulations for a series of novel thiazole/thiophene ketone amides with inhibitory activity at the RORγt receptor. The optimum CoMFA and CoMSIA models, derived from ligand-based superimposition I, exhibit leave-one-out cross-validated correlation coefficient (R²_cv) of .859 and .805, respectively. Furthermore, the external predictive abilities of the models were evaluated by a test set, producing the predicted correlation coefficient (R²_pred) of .7317 and .7097, respectively. In addition, molecular docking analysis was applied to explore the binding modes between the inhibitors and the receptor. MD simulation and MM/PBSA method were also employed to study the stability and rationality of the derived conformations, and the binding free energies in detail. The QSAR models and the results of molecular docking, MD simulation, binding free energies corroborate well with each other and further provide insights regarding the development of novel RORγt inhibitors with better activity.     

Homology modeling and docking studies of FabH (β-ketoacyl-ACP synthase III) enzyme involved in type II fatty acid biosynthesis of Chlorella variabilis: a potential algal feedstock for biofuel production

The concept of using microalgae as an alternative renewable source of biofuel has gained much importance in recent years. However, its commercial feasibility is still an area of concern for researchers. Unraveling the fatty acid metabolic pathway and understanding structural features of various key enzymes regulating the process will provide valuable insights to target microalgae for augmented oil content. FabH (β-ketoacyl-acyl carrier protein synthase; KAS III) is a condensing enzyme catalyzing the initial elongation step of type II fatty acid biosynthetic process and acyl carrier protein (ACP) facilitates the shuttling of the fatty acyl intermediates to the active site of the respective enzymes in the pathway. In the present study, a reliable three-dimensional structure of FabH from Chlorella variabilis, an oleaginous green microalga was modeled and subsequently the key residues involved in substrate binding were determined by employing protein–protein docking and molecular dynamics (MD) simulation protocols. The FabH-ACP complex having the lowest docking energy score showed the binding of ACP to the electropositive FabH surface with strong hydrogen bond interactions. The MD simulation results indicated that the substrate-complexed FabH adopted a more stable conformation than the free enzyme. Further, the FabH structure retained its stability throughout the simulation although noticeable displacements were observed in the loop regions. Molecular simulation studies suggested the importance of crucial hydrogen bonding of the conserved Arg⁹¹ of FabH with Glu⁵³ and Asp⁵⁶ of ACP for exhibiting high affinity between the enzyme and substrate. The molecular modeling results are consistent with available experimental results on the flexibility of FabH and the present study provides first in silico insights into the structural and dynamical aspect of catalytic mechanism of FabH, which could be used for further site-specific mutagenic experiments to develop engineered high oil-yielding microalgal strains for biofuel production.     

Structure-based approach for the study of thyroid hormone receptor binding affinity and subtype selectivity

Thyroid hormone (TH) possesses the ability to lower cholesterol and improve cardiac performance, which have prompted the efforts to design analogs that can utilize the cholesterol-lowering property without adversely affecting heart function. In order to gain insights into the interaction mechanism for agonists at the active site of thyroid hormone receptor β (TRβ), quantitative structure–activity relationship (QSAR) models have been developed on TRβ agonists, significant statistical coefficients were obtained (CoMFA, R²_cv, .732), (CoMSIA, R²_cv, .853), indicating the internal consistency of the models, the obtained models were further validated using the test set, the acquired R²_pred values .7054 and .7129 were in good agreement with the experimental results. The key amino acids affecting ligand binding were identified by molecular docking, and the detailed binding modes of the compounds with different activities were also determined. Furthermore, molecular dynamics (MD) simulations were conducted to assess the reliability of the derived models and the docking results. Moreover, TH exerts significant physiological effects through modulation of the two human thyroid hormone receptor subtypes. Because TRβ and TRα locate in different target cells, selective TR ligands would target specific tissues regulated by one receptor without affecting the other. Thus, the 3D information was analyzed to reveal the most relevant structural features involved in selectivity. The findings serve as the basis for further investigation into selective TRβ/TRα agonists.     

Molecular modeling studies,synthesis and biological evaluation of dabigatran analogues as thrombin inhibitors

In this work, 48 thrombin inhibitors based on the structural scaffold of dabigatran were analyzed using a combination of molecular modeling techniques. We generated three-dimensional quantitative structure–activity relationship (3D-QSAR) models based on three alignments for both comparative molecular field analysis (CoMFA) and comparative molecular similarity index analysis (CoMSIA) to highlight the structural requirements for thrombin protein inhibition. In addition to the 3D-QSAR study, Topomer CoMFA model also was established with a higher leave-one-out cross-validation q² and a non-cross-validation r², which suggest that the three models have good predictive ability. The results indicated that the steric, hydrophobic and electrostatic fields play key roles in QSAR model. Furthermore, we employed molecular docking and re-docking simulation explored the binding relationship of the ligand and the receptor protein in detail. Molecular docking simulations identified several key interactions that were also indicated through 3D-QSAR analysis. On the basis of the obtained results, two compounds were designed and predicted by three models, the biological evaluation in vitro (IC₅₀) demonstrated that these molecular models were effective for the development of novel potent thrombin inhibitors.     

Assessment of the interaction procedure between Pt(IV) prodrug [Pt(5,5′-dmbpy)Cl4 and human serum albumin: Combination of spectroscopic and molecular modeling technique

In this study, a cytotoxic Pt(IV) complex [Pt(5,5′-dmbpy)Cl₄ (5,5′-dmbpy is 5,5′-dimethyl-2,2′-bipyridine) was selected to investigate its affinity to human serum albumin (HSA) by spectroscopy and molecular docking methods. This complex has a bidentate nitrogen donor ligand with four chloride anions attached to a Pt(IV) metal in a distorted octahedral environment. The ?uorescence data showed this complex quench the intrinsic ?uorescence of HSA through a static quenching mechanism. The binding constant (K_b) and the number of binding sites (n) were obtained based on the results of fluorescence measurements. UV–vis, circular dichroism spectroscopy, and three-dimensional fluorescence spectroscopy proved that the Pt(IV) complex could slightly change the secondary structure of protein. Thermodynamic parameters show that the Pt(IV) complex binds to HSA through electrostatic and Vander Waals interactions with one binding site. The molecular docking results confirmed the spectroscopic results and showed that Pt(IV) complex is embedded into subdomain IIA of protein. The aim of this study is to describe the performance of effective anti-cancer drugs when faced with proteins such as HSA.     

Investigation on the binding mode of benzothiophene analogues as potent factor IXa (FIXa) inhibitors in thrombosis by CoMFA,docking and molecular dynamic studies

Recently, benzothiophenes attract much attention of interest due to its possible inhibitory activity targeting FIXa, a blood coagulation factor that is essential for the amplification or consolidation phase of blood coagulation. To explore this inhibitory mechanism, three-dimensional quantitative structure–activity relationship (3D-QSAR), molecular docking and molecular dynamics (MD) studies on a series of 84 benzothiophene analogues, for the first time, were performed. As a result, a highly predictive CoMFA model was developed with the q²?=?0.52, r²?=?0.97 and r²_pred?=?0.81, respectively. The CoMFA contour maps, the docking analysis, as well as the MD simulation results are all in a good agreement, proving the reliability and robustness of the model. These models and the information, we hoped, would be helpful in screening and development of novel drugs against thrombosis prior to synthesis.     

Induced fit docking,free energy calculation and molecular dynamics studies on Mycobacterium tuberculosis alanine racemase inhibitor

Alar, a Pyridoxal 5′-phosphate (PLP)-dependent bacterial enzyme is responsible for the racemisation of L-alanine into D-alanine which is essential for the peptidoglycan biosynthesis in both Gram-positive and Gram-negative bacteria. In the present study, we performed induced fit docking, binding free energy calculation and molecular dynamics simulation to elucidate the Alar inhibition potential of 1,2,4-thiadiazolidine-3,5-dione-based inhibitor 1. The inhibitor binds to the hydrophobic groove of Alar and the binding was found to be stable throughout 20-ns MD simulation. Induced fit docking result showed that Lys42, Tyr46, Tyr175 and Tyr364 residues are primarily responsible for the stabilisation of inhibitor–protein complex. Further, high negative van der Waals binding free energy value of –38.88 kcal/mol, indicated it as the main driving force for the inhibitor binding. Based on the information obtained from this study, we designed few molecules as potent Alar inhibitor. In order to gain structural insight and to validate the stability of complex, we performed 20-ns MD simulation of the designed molecule D1. Results obtained from this study can be used for the design of M. tuberculosis Alar potent inhibitors lacking affinity for the co-factor PLP.     

Curcumin specifically binds to the human calcium–calmodulin-dependent protein kinase IV: fluorescence and molecular dynamics simulation studies

Calcium–calmodulin-dependent protein kinase IV (CAMK4) plays significant role in the regulation of calcium-dependent gene expression, and thus, it is involved in varieties of cellular functions such as cell signaling and neuronal survival. On the other hand, curcumin, a naturally occurring yellow bioactive component of turmeric possesses wide spectrum of biological actions, and it is widely used to treat atherosclerosis, diabetes, cancer, and inflammation. It also acts as an antioxidant. Here, we studied the interaction of curcumin with human CAMK4 at pH 7.4 using molecular docking, molecular dynamics (MD) simulations, fluorescence binding, and surface plasmon resonance (SPR) methods. We performed MD simulations for both neutral and anionic forms of CAMK4-curcumin complexes for a reasonably long time (150 ns) to see the overall stability of the protein–ligand complex. Molecular docking studies revealed that the curcumin binds in the large hydrophobic cavity of kinase domain of CAMK4 through several hydrophobic and hydrogen-bonded interactions. Additionally, MD simulations studies contributed in understanding the stability of protein–ligand complex system in aqueous solution and conformational changes in the CAMK4 upon binding of curcumin. A significant increase in the fluorescence intensity at 495 nm was observed (λ_exc = 425 nm), suggesting a strong interaction of curcumin to the CAMK4. A high binding affinity (K_D = 3.7 × 10^?8 ± .03 M) of curcumin for the CAMK4 was measured by SPR further indicating curcumin as a potential ligand for the CAMK4. This study will provide insights into designing a new inspired curcumin derivatives as therapeutic agents against many life-threatening diseases.     

Exploring the role of the phospholipid ligand in endothelial protein C receptor: A molecular dynamics study

Endothelial protein C receptor (EPCR) is a CD1‐like transmembrane glycoprotein with important regulatory roles in protein C (PC) pathway, enhancing PC's anticoagulant, anti‐inflammatory, and antiapoptotic activities. Similarly to homologous CD1d, EPCR binds a phospholipid [phosphatidylethanolamine (PTY)] in a groove corresponding to the antigen‐presenting site, although it is not clear if lipid exchange can occur in EPCR as in CD1d. The presence of PTY seems essential for PC γ‐carboxyglutamic acid (Gla) domain binding. However, the lipid‐free form of the EPCR has not been characterized. We have investigated the structural role of PTY on EPCR, by multiple molecular dynamics (MD) simulations of ligand bound and unbound forms of the protein. Structural changes, subsequent to ligand removal, led to identification of two stable and folded ligand‐free conformations. Compared with the bound form, unbound structures showed a narrowing of the A′ pocket and a high flexibility of the helices around it, in agreement with CD1d simulation. Thus, a lipid exchange with a mechanism similar to CD1d is proposed. In addition, unbound conformations presented a reduced interaction surface for Gla domain, confirming the role of PTY in establishing the proper EPCR conformation for the interaction with its partner protein. Single MD simulations were also obtained for 29 mutant models with predicted structural stability and impaired binding ability. Ligand affinity calculations, based on linear interaction energy method, showed that substitution‐induced conformational changes affecting helices around the A′ pocket were associated to a reduced binding affinity. Mutants responsible for this effect may represent useful reagents for experimental tests. Proteins 2010. © 2010 Wiley‐Liss, Inc.     

Comparison of computational model and X-ray crystal structure of human serotonin transporter: potential application for the pharmacology of human monoamine transporters

Abstract

The human serotonin transporter (hSERT) played a significant role in neurological process whose structural basis had been analysed for many years. Recently, the first homology model was constructed for hSERT based on the crystal structure of drosophila melanogaster dopamine transporter was published, and the inhibitory mechanism underlying the binding mode between hSERT and approved antidepressants was substantially investigated by molecular dynamics (MD) simulation. Right after this publication, the X-ray crystallographic structures of hSERT were reported, which provided a good opportunity to reassess the performance of previous simulation. In this study, the analyses of side-chain contact map, stereochemical quality and ligand-binding pocket were firstly conducted, which revealed that the constructed homology model of hSERT could successfully reproduce the reported crystal structure. Secondly, the approved antidepressant escitalopram was docked into the X-ray structure, and its binding pose was consistent with the reported docking pose in the homology model. Finally, MD simulation were performed based on the crystal structure of hSERT, and structural features revealed as critical for escitalopram-hSERT interaction by previous simulation were successfully recaptured. Thus, the newly reported X-ray crystal structure of hSERT was precisely predicted by computational model, which demonstrated its reliability in understanding the pharmacology of other human monoamine transporters whose 3-D structure remained unknown.     

Prediction of HIV-1 entry inhibitors neomycin-arginine conjugates interaction with the CD4-gp120 binding site by molecular modeling and multistep docking procedure

Developing of multi-target HIV-1 entry inhibitors represents an important avenue of drug therapy. Two such inhibitors are hexa-arginine-neomycin-conjugate (NeoR6) and nona-d-arginine-neomycin-conjugate (Neo-r9). Our findings that NeoR6-resistant mutations appear in the gp120 constant regions; and NeoR6 is not CCR5 antagonist, but inhibits CXCR4 and CCR5 HIV-1 using isolates, led us to suggest that NeoR6 may inhibit HIV-1 entry by interfering with the CD4-gp120 binding. To support this notion, we constructed a homology model of unliganded HIV-1_IIIB gp120 and docked NeoR6 and Neo-r9 to it, using a multistep docking procedure: geometric-electrostatic docking by MolFit; flexible ligand docking by Autodock3 and final refinement of the obtained complexes by Discover3. Binding free energies were calculated by MM-PBSA methodology. The model predicts competitive inhibition of CD4-gp120 binding by NeoR6 and Neo-r9. We determined plausible binding sites between constructed CD4-bound gp120 trimer and homology modeled membranal CXCR4, and tested NeoR6 and Neo-r9 interfering with this interaction. These models support our notion that another mechanism of anti-HIV-1 activity of NeoR6 is inhibition of gp120-CXCR4 binding. These structural models and interaction of NeoR6 and Neo-r9 with gp120 and CXCR4 provide a powerful approach for structural based drug design for selective targeting of HIV-1 entry and/or for inhibition of other retroviruses with similar mechanism of entry.     

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