Crystallographic studies of the complexes of antiviral protein griffithsin with glucose and N-acetylglucosamine

Crystal structures of complexes of an antiviral lectin griffithsin (GRFT) with glucose and N-acetylglucosamine were solved and refined at high resolution. In both complexes, all six monosaccharide-binding sites of GRFT were occupied and the mode of binding was similar to that of mannose. In our previous attempts to obtain a complex with N-acetylglucosamine by soaking, only a single site was occupied; thus, cocrystallization was clearly superior despite lower concentration of the ligand. Isothermal titration calorimetric experiments with N-acetylglucosamine, glucose, and mannose provided enthalpic evidence of distinct binding differences between the three monosaccharides. A comparison of the mode of binding of different monosaccharides is discussed in the context of the antiviral activity of GRFT, based on specific binding to high-mannose-containing complex carbohydrates found on viral envelopes.     

Domain-swapped structure of the potent antiviral protein griffithsin and its mode of carbohydrate binding

The crystal structure of griffithsin, an antiviral lectin from the red alga Griffithsia sp., was solved and refined at 1.3 A resolution for the free protein and 0.94 A for a complex with mannose. Griffithsin molecules form a domain-swapped dimer, in which two beta strands of one molecule complete a beta prism consisting of three four-stranded sheets, with an approximate 3-fold axis, of another molecule. The structure of each monomer bears close resemblance to jacalin-related lectins, but its dimeric structure is unique. The structures of complexes of griffithsin with mannose and N-acetylglucosamine defined the locations of three almost identical carbohydrate binding sites on each monomer. We have also shown that griffithsin is a potent inhibitor of the coronavirus responsible for severe acute respiratory syndrome (SARS). Antiviral potency of griffithsin is likely due to the presence of multiple, similar sugar binding sites that provide redundant attachment points for complex carbohydrate molecules present on viral envelopes.     

Theoretical investigation on the glycan‐binding specificity of Agrocybe cylindracea galectin using molecular modeling and molecular dynamics simulation studies

Galectins are β‐galactoside binding proteins which have the ability to serve as potent antitumor, cancer biomarker, and induce tumor cell apoptosis. Agrocybe cylindracea galectin (ACG) is a fungal galectin which specifically recognizes α(2,3)‐linked sialyllactose at the cell surface that plays extensive roles in the biological recognition processes. To investigate the change in glycan‐binding specificity upon mutations, single point and double point site‐directed in silico mutations are performed at the binding pocket of ACG. Molecular dynamics (MD) simulation studies are carried out for the wild‐type (ACG) and single point (ACG1) and double point (ACG2) mutated ACGs to investigate the dynamics of substituted mutants and their interactions with the receptor sialyllactose. Plausible binding modes are proposed for galectin–sialylglycan complexes based on the analysis of hydrogen bonding interactions, total pair‐wise interaction energy between the interacting binding site residues and sialyllactose and binding free energy of the complexes using molecular mechanics–Poisson–Boltzmann surface area. Our result shows that high contribution to the binding in different modes is due to the direct and water‐mediated hydrogen bonds. The binding specificity of double point mutant Y59R/N140Q of ACG2 is found to be high, and it has 26 direct and water‐mediated hydrogen bonds with a relatively low‐binding free energy of −47.52 ± 5.2 kcal/mol. We also observe that the substituted mutant Arg59 is crucial for glycan‐binding and for the preference of α(2,3)‐linked sialyllactose at the binding pocket of ACG2 galectin. When compared with the wild‐type and single point mutant, the double point mutant exhibits enhanced affinity towards α(2,3)‐linked sialyllactose, which can be effectively used as a model for biological cell marker in cancer therapeutics. Copyright © 2015 John Wiley & Sons, Ltd.     

Allosteric HIV‐1 integrase inhibitors promote aberrant protein multimerization by directly mediating inter‐subunit interactions: Structural and thermodynamic modeling studies

Allosteric HIV‐1 integrase (IN) inhibitors (ALLINIs) bind at the dimer interface of the IN catalytic core domain (CCD), and potently inhibit HIV‐1 by promoting aberrant, higher‐order IN multimerization. Little is known about the structural organization of the inhibitor‐induced IN multimers and important questions regarding how ALLINIs promote aberrant IN multimerization remain to be answered. On the basis of physical chemistry principles and from our analysis of experimental information, we propose that inhibitor‐induced multimerization is mediated by ALLINIs directly promoting inter‐subunit interactions between the CCD dimer and a C‐terminal domain (CTD) of another IN dimer. Guided by this hypothesis, we have built atomic models of inter‐subunit interfaces in IN multimers by incorporating information from hydrogen‐deuterium exchange (HDX) measurements to drive protein‐protein docking. We have also developed a novel free energy simulation method to estimate the effects of ALLINI binding on the association of the CCD and CTD. Using this structural and thermodynamic modeling approach, we show that multimer inter‐subunit interface models can account for several experimental observations about ALLINI‐induced multimerization, including large differences in the potencies of various ALLINIs, the mechanisms of resistance mutations, and the crucial role of solvent exposed R‐groups in the high potency of certain ALLINIs. Our study predicts that CTD residues Tyr226, Trp235 and Lys266 are involved in the aberrant multimer interfaces. The key finding of the study is that it suggests the possibility of ALLINIs facilitating inter‐subunit interactions between an external CTD and the CCD‐CCD dimer interface.     

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     

Thermodynamic dissection of the binding energetics of KNI-272, a potent HIV-1 protease inhibitor

KNI-272 is a powerful HIV-1 protease inhibitor with a reported inhibition constant in the picomolar range. In this paper, a complete experimental dissection of the thermodynamic forces that define the binding affinity of this inhibitor to the wild-type and drug-resistant mutant V82F/184V is presented. Unlike other protease inhibitors, KNI-272 binds to the protease with a favorable binding enthalpy. The origin of the favorable binding enthalpy has been traced to the coupling of the binding reaction to the burial of six water molecules. These bound water molecules, previously identified by NMR studies, optimize the atomic packing at the inhibitor/protein interface enhancing van der Waals and other favorable interactions. These interactions offset the unfavorable enthalpy usually associated with the binding of hydrophobic molecules. The association constant to the drug resistant mutant is 100-500 times weaker. The decrease in binding affinity corresponds to an increase in the Gibbs energy of binding of 3-3.5 kcal/mol, which originates from less favorable enthalpy (1.7 kcal/mol more positive) and entropy changes. Calorimetric binding experiments performed as a function of pH and utilizing buffers with different ionization enthalpies have permitted the dissection of proton linkage effects. According to these experiments, the binding of the inhibitor is linked to the protonation/deprotonation of two groups. In the uncomplexed form these groups have pKs of 6.0 and 4.8, and become 6.6 and 2.9 in the complex. These groups have been identified as one of the aspartates in the catalytic aspartyl dyad in the protease and the isoquinoline nitrogen in the inhibitor molecule. The binding affinity is maximal between pH 5 and pH 6. At those pH values the affinity is close to 6 x 10(10) M(-1) (Kd = 16 pM). Global analysis of the data yield a buffer- and pH-independent binding enthalpy of -6.3 kcal/mol. Under conditions in which the exchange of protons is zero, the Gibbs energy of binding is -14.7 kcal/mol from which a binding entropy of 28 cal/K mol is obtained. Thus, the binding of KNI-272 is both enthalpically and entropically favorable. The structure-based thermodynamic analysis indicates that the allophenylnorstatine nucleus of KNI-272 provides an important scaffold for the design of inhibitors that are less susceptible to resistant mutations.     

Study on the inhibitory mechanism and binding mode of the hydroxycoumarin compound NSC158393 to HIV‐1 integrase by molecular modeling

Human immunodeficiency virus type 1 integrase (IN) is an essential enzyme in the life cycle of this virus and also an important target for the study of anti‐HIV drugs. In this work, the binding modes of the wild type IN core domain and the two mutants, that is, W132G and C130S, with the 4‐hydroxycoumarin compound NSC158393 were evaluated by using the “relaxed complex” molecular docking approach combined with molecular dynamics (MD) simulations. Based on the monomer MD simulations, both of the two substitutions affect not only the stability of the 128–136 peptides, but also the flexibility of the functional 140s loop. In principle, NSC158393 binds the 128–136 peptides of IN; however, the specific binding modes for the three systems are various. According to the binding mode of NSC158393 with WT, NSC158393 can effectively interfere with the stability of the IN dimer by causing a steric hindrance around the monomer interface. Additionally, through the comparative analysis of the MD trajectories of the wild type IN and the IN‐NSC158393 complex, we found that NSC15893 may also exert its inhibitory function by diminishing the mobility of the function loop of IN. Three key binding residues, that is, W131, K136, and G134, were discovered by energy decomposition calculated with the Molecular Mechanics Generalized Born Surface Area method. Characterized by the largest binding affinity, W131 is likely to be indispensable for the ligand binding. All the above results are consistent with experiment data, providing us some helpful information for understanding the mechanism of the coumarin‐based inhibitors. © 2009 Wiley Periodicals, Inc. Biopolymers 91: 700–709, 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     

Differential binding of SARS-CoV-2 Spike protein variants to its cognate receptor hACE2 using molecular modeling based binding analysis

The current emergence of novel coronavirus, SARS-CoV-2 and its ceaseless expansion worldwide has posed a global health emergency that has adversely affected the humans. With the entire world striving to understand the newly emerged virus, differences in morbidity and infection rate of SARS-CoV-2 have been observed across varied geographic areas, which have been ascribed to viral mutation and evolution over time. The homotrimeric Spike (S) glycoprotein on the viral envelope surface is responsible for binding, priming, and initiating infection in the host. Our phylogeny analysis of 1947 sequences of S proteins indicated there is a change in amino acid (aa) from aspartate (Group-A) to glycine (Group-B) at position 614, near the receptor- binding domain (RBD; aa positions 331-524). The two variants are reported to be in circulation, disproportionately across the world, with Group-A dominant in Asia and Group-B in North America. The trimeric, monomeric, and RBD of S protein of both the variant groups (A & B) were modeled using the Swiss-Model server and were docked with the human receptor angiotensin-converting enzyme 2 (hACE2) employing the PatchDock server and visualized in PyMol. Group-A S protein''s RBD bound imperceptibly to the two binding clefts of the hACE2 protein, on the other hand, Group-B S protein''s RBD perfectly interacted inside the binding clefts of hACE2, with higher number of hydrogen and hydrophobic interactions. This implies that the S protein''s amino acid at position 614 near the core RBD influences its interaction with the cognate hACE2 receptor, which may induce its infectivity that should be explored further with molecular and biochemical studies.     

Spectroscopic and molecular modeling studies of the interaction between morin and polyamidoamine dendrimer

Interactions between the polyamidoamine (PAMAM) dendrimer and drug molecules are of interest for their potential biomedical applications. The goal of this work is to examine the interaction of PAMAM‐C₁₂ 25% dendrimer with morin. The ultraviolet–visible, fluorescence spectroscopic methods as well as molecular modeling were used to analyze drug‐binding mode, binding constants and binding sites, etc. The experimental data showed that the binding constant of morin‐PAMAM‐C₁₂ 25% is about 10⁵ L/mol. The interaction of morin with PAMAM‐C₁₂ 25% is mainly driven by the hydrophobic, electrostatic, hydrogen bonds and van der Waals forces. There are mainly three classes of binding site of morin at the interface of PAMAM‐C₁₂ 25%. These results provided some useful information for self‐assembling and disassembling the PAMAM dendrimer as well as efficient drug delivery and therapeutic applications. Copyright © 2013 John Wiley & Sons, Ltd.     

Design, synthesis, binding, and molecular modeling studies of new potent ligands of cannabinoid receptors

In our ongoing program aimed at the design, synthesis, and biological evaluation of novel cannabinoid receptor ligands derived from olivetol and hexyl-resorcinol, we have designed a structural model for new derivatives on the basis of a previous study. Here we report the synthesis, binding, and molecular modeling studies of new potent compounds with high affinity toward CB(1) and CB(2) receptors. Compounds with amidic 'heads' with alkyloxy chains varying in length from 8 to 12 carbon atoms showed nanomolar affinity for both receptors, depending on the type of aromatic backbone. Two of the new compounds, although not very potent, exhibit selectivity for CB(1) receptors (CB(1)/CB(2)=0.07 and 0.08, respectively). Molecular modeling studies fitted this new class of cannabinoid ligands into a CB(1) receptor model, and the qualitative analysis of the results was in general agreement with the CB(1) affinity constants observed experimentally for these derivatives.     

The active site and substrates binding mode of malonyl-CoA synthetase determined by transferred nuclear Overhauser effect spectroscopy, site-directed mutagenesis, and comparative modeling studies

The active sites and substrate bindings of Rhizobium trifolii molonyl-CoA synthetase (MCS) catalyzing the malonyl-CoA formation from malonate and CoA have been determined based on NMR spectroscopy, site-directed mutagenesis, and comparative modeling methods. The MCS-bound conformation of malonyl-CoA was determined from two-dimensional-transferred nuclear Overhauser effect spectroscopy data. MCS protein folds into two structural domains and consists of 16 alpha-helices, 24 beta-strands, and several long loops. The core active site was determined as a wide cleft close to the end of the small C-terminal domain. The catalytic substrate malonate is placed between ATP and His206 in the MCS enzyme, supporting His206 in its catalytic role as it generates reaction intermediate, malonyl-AMP. These findings are strongly supported by previous biochemical data, as well as by the site-directed mutagenesis data reported here. This structure reveals the biochemical role as well as the substrate specificity that conservative residues of adenylate-forming enzymes have.     

Predicting peptide binding to MHC pockets via molecular modeling, implicit solvation, and global optimization

Development of a computational prediction method based on molecular modeling, global optimization, and implicit solvation has produced accurate structure and relative binding affinity predictions for peptide amino acids binding to five pockets of the MHC molecule HLA-DRB1*0101. Because peptide binding to MHC molecules is essential to many immune responses, development of such a method for understanding and predicting the forces that drive binding is crucial for pharmaceutical design and disease treatment. Underlying the development of this prediction method are two hypotheses. The first is that pockets formed by the peptide binding groove of MHC molecules are independent, separating the prediction of peptide amino acids that bind within individual pockets from those that bind between pockets. The second hypothesis is that the native state of a system composed of an amino acid bound to a protein pocket corresponds to the system's lowest free energy. The prediction method developed from these hypotheses uses atomistic-level modeling, deterministic global optimization, and three methods of implicit solvation: solvent-accessible area, solvent-accessible volume, and Poisson-Boltzmann electrostatics. The method predicts relative binding affinities of peptide amino acids for pockets of HLA-DRB1*0101 by determining computationally an amino acid's global minimum energy conformation. Prediction results from the method are in agreement with X-ray crystallography data and experimental binding assays.     

Understanding the binding mode and function of BMS-488043 against HIV-1 viral entry

A recently discovered small-molecule inhibitor, BMS-488043 (BMS-488), for the invasion of Human immunodeficiency virus Type 1 (HIV-1), shows a high activity against the entry of diversified HIV-1. Docking and molecular dynamic studies have been carried out to understand the binding mode of BMS-488 to gp120 as well as the effect of the small molecule on the conformational change of gp120 induced by CD4 binding. The results indicate that BMS-488 can accommodate in the CD4 binding pocket and interfere the CD4 binding in a noncompetitive mode. The piperazine group of BMS-488 prevents the bridging sheet formation of gp120 induced by the CD4 binding mainly through blocking the rotation of the Trp112 located on the α1 helix of gp120. The bridging sheet formation cannot be blocked for the W112A mutant of gp120 due to the reduced steric hindrance, in agreement with its significant resistance to the BMS inhibitor. The aza-indole ring is likely to interfere the exposure of gp41 by stacking within the β3-β5 and LB loops to disrupt the close packing of Pro212-His66-Phe210. The mode of action of BMS-488 also accommodates many mutagenesis results related to BMS-488 activity.     

Homology modeling and molecular dynamics simulations of the glycine receptor ligand binding domain

We present a homology based model of the ligand binding domain (LBD) of the homopentameric alpha1 glycine receptor (GlyR). The model is based on multiple sequence alignment with other members of the nicotinicoid ligand gated ion channel superfamily and two homologous acetylcholine binding proteins (AChBP) from the freshwater (Lymnaea stagnalis) and saltwater (Aplysia californica) snails with known high resolution structure. Using two template proteins with known structure to model three dimensional structure of a target protein is especially advantageous for sequences with low homology as in the case presented in this paper. The final model was cross-validated by critical evaluation of experimental and published mutagenesis, functional and other biochemical studies. In addition, a complex structure with strychnine antagonist in the putative binding site is proposed based on docking simulation using Autodock program. Molecular dynamics (MD) simulations with simulated annealing protocol are reported on the proposed LBD of GlyR, which is stable in 5 ns simulation in water, as well as for a deformed LBD structure modeled on the corresponding domain determined in low-resolution cryomicroscopy structure of the alpha subunit of the full-length acetylcholine receptor (AChR). Our simulations demonstrate that the beta-sandwich central core of the protein monomer is fairly rigid in the simulations and resistant to deformations in water.     

Structural investigation of selective binding dynamics for the pheromone‐binding protein 1 of the grapevine moth,Lobesia botrana

The European grapevine moth, Lobesia botrana (Denis & Schiffermüller), is a serious pest in vineyards in North and South America. Mating disruption techniques have been used to control and monitor L. botrana on the basis of its sexual communication. This needs a well‐tuned olfactory system, in which it is believed that pheromone‐binding proteins (PBPs) are key players that transport pheromones in the antennae of moths. In this study, the selectivity of a PBP, named as LbotPBP1, was tested by fluorescence binding assays against 11 sex pheromone components and 6 host plant volatiles. In addition, its binding mechanism was predicted on the basis of structural analyses by molecular docking and complex and steered molecular dynamics (SMD). Our results indicate that LbotPBP1 binds selectively to sex pheromone components over certain host plant volatiles, according to both in vitro and in silico tests. Thus, chain length (14 carbon atoms) and functional groups (i.e., alcohol and ester) appear to be key features for stable binding. Likewise, residues such as Phe12, Phe36, and Phe118 could participate in unspecific binding processes, whilst Ser9, Ser56, and Trp114 could participate in the specific recognition and stabilization of sex pheromones instead of host plant volatiles. Moreover, our SMD approach supported 11‐dodecenyl acetate as the best ligand for LbotPBP1. Overall, the dynamics simulations, contact frequency analysis and SMD shed light on the binding mechanism of LbotPBP1 and could overcome the imprecision of molecular docking, supporting the in vitro binding assays. Finally, the role of LbotPBP1 in the chemical ecology of L. botrana is discussed.     

Three‐dimensional structure,binding, and spectroscopic characteristics of the monoclonal antibody 43.1 directed to the carboxyphenyl moiety of fluorescein

Unlike other known anti‐fluorescein antibodies, the monoclonal antibody 43.1 is directed toward the fluorescein's carboxyl phenyl moiety. It demonstrates a very high affinity (K_D ～ 70 pM) and a fast association rate (k_on ～ 2 × 10⁷ M^?1 s^?1). The three‐dimensional structure of the Fab 43.1—fluorescein complex was resolved at 2.4 Å resolution. The antibody binding site is exclusively assembled by the CDR loops. It is comprised of a 14 Å groove‐shaped entrance leading to a 9 Å by 7 Å binding pocket. The highly polar binding pocket complementary encloses the fluorescein's carboxyphenyl moiety and tightly fixes it by multiple hydrogen bonds. The fluorescein's xanthene ring is embedded in the more hydrophobic groove and stacked between the side chains of Tyr37_L and of Arg99_H providing conditions for an excited state electron transfer process. In comparison to fluorescein, the absorption spectrum of the complex in the visible region is shifted to the “red” by 23 nm. The complex demonstrates a very weak fluorescence (Φ_c = 0.0018) with two short lifetime components: 0.03 ns (47%) and 0.8 ns (24%), which reflects a 99.8% fluorescein emission quenching effect upon complex formation. The antibody 43.1 binds fluorescein with remarkable affinity, fast association rate, and strongly quenches its emission. Therefore, it may present a practical interest in applications such as molecular sensors and switches. © 2016 Wiley Periodicals, Inc. Biopolymers 105: 234–243, 2016.     

Structure modeling, ligand binding, and binding affinity calculation (LR-MM-PBSA) of human heparanase for inhibition and drug design

Human heparanase is an endo-beta-D-glycosidase that cleaves heparan sulphate (HS) chains in the extracellular matrix and basement membrane. It is known that the cleavage of HS by heparanase results in cell invasion and metastasis of cancer. Therefore, heparanase is considered an important target for cancer drug development. The three-dimensional structure of heparanase would be useful in the rational design of inhibitors targeted to the enzyme; however, the three-dimensional structure has not yet been determined. In our effort to design inhibitors, we developed a three-dimensional structure of heparanase using a homology-modeling approach. The homology-built structure is consistent to previous bioinformatics and site-mutation experimental results. The heparanase features a (alpha/beta)(8) TIM-barrel fold with two glutamate residues (Glu225 and Glu343) located in the active-site cleft. This feature supports the putative mechanism of proton donor and nucleophilic sites. Docking simulations yielded 41 complex structures, which indicate that the bound inhibitor could block ligand binding into the catalytic site. A free energy of binding model was established for 25 heparanase inhibitors with a training set of 25 heparanase inhibitors using the linear response MM-PBSA approach (LR-MM-PBSA). The correlation between calculated and experimental activity was 0.79 and the reliability of the model was validated with leave-one-out cross-validation method. Its predictive capability was further validated using a test set of 16 inhibitors similar to the training set of inhibitors. The correlation between the predicted and observed activities is significantly improved by the protein "induced-fit" that accounts for the flexibility of the receptor. These interaction and pharmacophore elements provide a unique insight to the rational design of new ligands targeted to the enzyme.     

Ligand binding and thermodynamic stability of a multidomain protein, calmodulin

Chemical and thermal denaturation of calmodulin has been monitored spectroscopically to determine the stability for the intact protein and its two isolated domains as a function of binding of Ca2+ or Mg2+. The reversible urea unfolding of either isolated apo-domain follows a two-state mechanism with relatively low deltaG(o)20 values of approximately 2.7 (N-domain) and approximately 1.9 kcal/mol (C-domain). The apo-C-domain is significantly unfolded at normal temperatures (20-25 degrees C). The greater affinity of the C-domain for Ca2+ causes it to be more stable than the N-domain at [Ca2+] > or = 0.3 mM. By contrast, Mg2+ causes a greater stabilization of the N- rather than the C-domain, consistent with measured Mg2+ affinities. For the intact protein (+/-Ca2+), the bimodal denaturation profiles can be analyzed to give two deltaG(o)20 values, which differ significantly from those of the isolated domains, with one domain being less stable and one domain more stable. The observed stability of the domains is strongly dependent on solution conditions such as ionic strength, as well as specific effects due to metal ion binding. In the intact protein, different folding intermediates are observed, depending on the ionic composition. The results illustrate that a protein of low intrinsic stability is liable to major perturbation of its unfolding properties by environmental conditions and liganding processes and, by extension, mutation. Hence, the observed stability of an isolated domain may differ significantly from the stability of the same structure in a multidomain protein. These results address questions involved in manipulating the stability of a protein or its domains by site directed mutagenesis and protein engineering.     

NMR-based modeling and binding studies of a ternary complex between chicken liver bile acid binding protein and bile acids

Chicken liver bile acid binding protein (cL-BABP) is involved in bile acid transport in the liver cytosol. A detailed study of the mechanism of binding and selectivity of bile acids binding proteins towards the physiological pool of bile salts is a key issue for the complete understanding of the role of these proteins and their involvement in cholesterol homeostasis. In the present study, we modeled the ternary complex of cL-BABP with two molecules of bile salts using the data driven docking program HADDOCK on the basis of NMR and mass spectrometry data. Docking resulted in good 3D models, satisfying the majority of experimental restraints. The docking procedure represents a necessary step to help in the structure determination and in functional analysis of such systems, in view of the high complexity of the 3D structure determination of a ternary complex with two identical ligands. HADDOCK models show that residues involved in binding are mainly located in the C-terminal end of the protein, with two loops, CD and EF, playing a major role in ligand binding. A spine, comprising polarresidues pointing toward the protein interior and involved in motion communication, has a prominent role in ligand interaction. The modeling approach has been complemented with NMR interaction and competition studies of cL-BABP with chenodeoxycholic and cholic acids. A higher affinity for chenodeoxycholic acid was observed and a Kd upper limit estimate was obtained. The binding is highly cooperative and no site selectivity was detected for the different bile salts, thus indicating that site selectivity and cooperativity are not correlated. Differences in physiological pathways and bile salt pools in different species is discussed in light of the binding results thus enlarging the body of knowledge of BABPs biological functions.     

Targeting hepatocellular carcinoma: Synthesis of new pyrazole-based derivatives,biological evaluation,DNA binding,and molecular modeling studies

A new series of pyrazole derivatives was prepared in this work, including pyrazolopyrimidines, pyrazolotriazines, pyrazolylthienopyridines, and 2-(pyrazolylamino)thiazol-4-ones, utilizing 3-amino-5-methyl-1H-pyrazole as a synthetic precursor. Their in vitro anticancer activity was tested on hepatocellular carcinoma cell line, HepG2. The results revealed that the pyrazolylhydrazonoyl cyanide 8, the pyrazolopyrimidine 3, and the pyrazolylaminothiazolone 17 were the most active with IC₅₀ values of 2, 7, and 7 µM respectively in comparison with 5.5 µM for cisplatin as a reference drug. Interestingly, all the synthesized compounds showed higher selectivity index than cisplatin. DNA binding assay was also carried out for the synthesized compounds to rationalize their mechanism of action. Molecular modeling studies, including docking into DNA minor groove, flexible alignment, and surface mapping, were conducted. Results obtained proved the superior DNA-binding affinity of the most active anticancer compounds.