Predicting absolute ligand binding free energies to a simple model site

A central challenge in structure-based ligand design is the accurate prediction of binding free energies. Here we apply alchemical free energy calculations in explicit solvent to predict ligand binding in a model cavity in T4 lysozyme. Even in this simple site, there are challenges. We made systematic improvements, beginning with single poses from docking, then including multiple poses, additional protein conformational changes, and using an improved charge model. Computed absolute binding free energies had an RMS error of 1.9 kcal/mol relative to previously determined experimental values. In blind prospective tests, the methods correctly discriminated between several true ligands and decoys in a set of putative binders identified by docking. In these prospective tests, the RMS error in predicted binding free energies relative to those subsequently determined experimentally was only 0.6 kcal/mol. X-ray crystal structures of the new ligands bound in the cavity corresponded closely to predictions from the free energy calculations, but sometimes differed from those predicted by docking. Finally, we examined the impact of holding the protein rigid, as in docking, with a view to learning how approximations made in docking affect accuracy and how they may be improved.     

Structural and binding studies of the three-metal center in two mycobacterial PPM Ser/Thr protein phosphatases

Phospho-Ser/Thr protein phosphatases (PPs) are dinuclear metalloenzymes classed into two large families, PPP and PPM, on the basis of sequence similarity and metal ion dependence. The archetype of the PPM family is the α isoform of human PP2C (PP2Cα), which folds into an α/β domain similar to those of PPP enzymes. The recent structural studies of three bacterial PPM phosphatases, Mycobacterium tuberculosis MtPstP, Mycobacterium smegmatis MspP, and Streptococcus agalactiae STP, confirmed the conservation of the overall fold and dinuclear metal center in the family, but surprisingly revealed the presence of a third conserved metal-binding site in the active site. To gain insight into the roles of the three-metal center in bacterial enzymes, we report structural and metal-binding studies of MtPstP and MspP. The structure of MtPstP in a new trigonal crystal form revealed a fully active enzyme with the canonical dinuclear metal center but without the third metal ion bound to the catalytic site. The absence of metal correlates with a partially unstructured flap segment, indicating that the third manganese ion contributes to reposition the flap, but is dispensable for catalysis. Studies of metal binding to MspP using isothermal titration calorimetry revealed that the three Mn²⁺-binding sites display distinct affinities, with dissociation constants in the nano- and micromolar range for the two catalytic metal ions and a significantly lower affinity for the third metal-binding site. In agreement, the structure of inactive MspP at acidic pH was determined at atomic resolution and shown to lack the third metal ion in the active site. Structural comparisons of all bacterial phosphatases revealed positional variations in the third metal-binding site that are correlated with the presence of bound substrate and the conformation of the flap segment, supporting a role of this metal ion in assisting enzyme-substrate interactions.     

Development of prolactin receptor antagonists with reduced pH‐dependence of receptor binding

The cytokine hormone prolactin has a vast number of diverse functions. Unfortunately, it also exhibits tumor growth promoting properties, which makes the development of prolactin receptor antagonists a priority. Prolactin binds to its cognate receptor with much lower affinity at low pH than at physiological pH and since the extracellular environment around solid tumors often is acidic, it is desirable to develop antagonists that have improved binding affinity at low pH. The pK_a value of a histidine side chain is ～6.8 making histidine residues obvious candidates for examination. From evaluation of known molecular structures of human prolactin, of the prolactin receptor and of different complexes of the two, three histidine residues in the hormone–receptor binding site 1 were selected for mutational studies. We analyzed 10 variants by circular dichroism spectroscopy, affinity and thermodynamic characterization of receptor binding by isothermal titration calorimetry combined with in vitro bioactivity in living cells. Histidine residue 27 was recognized as a central hot spot for pH sensitivity and conservative substitutions at this site resulted in strong receptor binding at low pH. Pure antagonists were developed earlier and the histidine mutations were introduced within such background. The antagonistic properties were maintained and the high affinity at low pH conserved. The implications of these findings may open new areas of research in the field of prolactin cancer biology. Copyright © 2010 John Wiley & Sons, Ltd.     

Thermodynamics of protein‐cation interaction: Ca+2 and Mg+2 binding to the fifth binding module of the LDL receptor

The ligand binding domain of the LDL receptor (LDLR) contains seven structurally homologous repeats. The fifth repeat (LR5) is considered to be the main module responsible for the binding of lipoproteins LDL and β‐VLDL. LR5, like the other homologous repeats, is around 40‐residue long and contains three disulfide bonds and a conserved cluster of negatively charged residues surrounding a hexacoordinated calcium ion. The calcium coordinating cage is formed by the backbone oxygens of W193 and D198, and side‐chain atoms of D196, D200, D206, and E207. The functionality of LDLR is closely associated with the presence of calcium. Magnesium ions are to some extent similar to calcium ions. However, they appear to be involved in different physiological events and their concentrations in extracellular and intracellular compartments are regulated by different mechanisms. Whether magnesium ions can play a role in the complex cycle of LDLR internalization and recycling is not known. We report here a detailed study of the interaction between LR5 and these two cations combining ITC, emission fluorescence, high resolution NMR, and MD simulations, at extracellular and endosomal pHs. Our results indicate that the conformational stability and internal dynamics of LR5 are strongly modulated by the specific bound cation. It appears that the difference in binding affinity for these cations is somewhat compensated by their different concentrations in late LDL‐associated endosomes. While the mildly acidic and calcium‐depleted environment in late endosomes has been proposed to contribute significantly to LDL release, the presence of magnesium might assist in efficient LDLR recycling. Proteins 2010. © 2009 Wiley‐Liss, Inc.     

Thermodynamics of DNA binding and distortion by the hyperthermophile chromatin protein Sac7d

Sac7d is a hyperthermophile chromatin protein which binds non-specifically to the minor groove of duplex DNA and induces a sharp kink of 66 degrees with intercalation of valine and methionine side-chains. We have utilized the thermal stability of Sac7d and the lack of sequence specificity to define the thermodynamics of DNA binding over a wide temperature range. The binding affinity for poly(dGdC) was moderate at 25 degrees C (Ka = 3.5(+/-1.6) x 10(6) M(-1)) and increased by nearly an order of magnitude from 10 degrees C to 80 degrees C. The enthalpy of binding was unfavorable at 25 degrees C, and decreased linearly from 5 degrees C to 60 degrees C. A positive binding heat at 25 degrees C is attributed in part to the energy of distorting DNA, and ensures that the temperature of maximal binding affinity (75.1+/-5.6 degrees C) is near the growth temperature of Sulfolobus acidocaldarius. Truncation of the two intercalating residues to alanine led to a decreased ability to bend and unwind DNA at 25 degrees C with a small decrease in binding affinity. The energy gained from intercalation is slightly greater than the free energy penalty of bending duplex DNA. Surprisingly, reduced distortion from the double alanine substitution did not lead to a significant decrease in the heat of binding at 25 degrees C. In addition, an anomalous positive DeltaCp of binding was observed for the double alanine mutant protein which could not be explained by the change in polar and apolar accessible surface areas. Both the larger than expected binding enthalpy and the positive heat capacity can be explained by a temperature dependent structural transition in the protein-DNA complex with a Tm of 15-20 degrees C and a DeltaH of 15 kcal/mol. Data are discussed which indicate that the endothermic transition in the complex is consistent with DNA distortion.     

Thermodynamic and kinetic characterization of ligand binding to the purine riboswitch aptamer domain

Riboswitches are cis-acting genetic regulatory elements found commonly in bacterial mRNAs that consist of a metabolite-responsive aptamer domain coupled to a regulatory switch. Purine riboswitches respond to intracellular concentrations of either adenine or guanine/hypoxanthine to control gene expression. The aptamer domain of the purine riboswitch contains a pyrimidine residue (Y74) that forms a Watson-Crick base-pairing interaction with the bound purine nucleobase ligand that discriminates between adenine and guanine. We sought to understand the structural basis of this specificity and the mechanism of ligand recognition by the purine riboswitch. Here, we present the 2,6-diaminopurine-bound structure of a C74U mutant of the xpt-pbuX guanine riboswitch, along with a detailed thermodynamic and kinetic analysis of nucleobase recognition by both the native and mutant riboswitches. These studies demonstrate clearly that the pyrimidine at position 74 is the sole determinant of purine riboswitch specificity. In addition, the mutant riboswitch binds adenine and adenine derivatives well compared with the guanine-responsive riboswitch. Under our experimental conditions, 2,6-diaminopurine binds the RNA with DeltaH=-40.3 kcal mol(-1), DeltaS=-97.6 cal mol(-1)K(-1), and DeltaG=-10.73 kcal mol(-1). A kinetic determination of the slow rate (0.15 x 10(5)M(-1)s(-1) and 2.1 x 10(5)mM(-1)s(-1) for 2-aminopurine binding the adenine-responsive mutant riboswitch and 7-deazaguanine-binding guanine riboswitch, respectively) of association under varying experimental conditions allowed us to propose a mechanism for ligand recognition by the purine riboswitch. A conformationally dynamic unliganded state for the binding pocket is stabilized first by the Watson-Crick base pairing between the ligand and Y74, and by the subsequent ordering of the J2/3 loop, enclosing the ligand within the three-way junction.     

Production of soluble alphabeta T-cell receptor heterodimers suitable for biophysical analysis of ligand binding

A method to produce alphabeta T-cell receptors (TCRs) in a soluble form suitable for biophysical analysis was devised involving in vitro refolding of a TCR fusion protein. Polypeptides corresponding to the variable and constant domains of each chain of a human and a murine receptor, fused to a coiled coil heterodimerization motif from either c-Jun (alpha) or v-Fos (beta), were overexpressed separately in Escherichia coli. Following recovery from inclusion bodies, the two chains of each receptor were denatured, and then refolded together in the presence of denaturants. For the human receptor, which is specific for the immunodominant influenza A HLA-A2-restricted matrix epitope (M58-66), a heterodimeric protein was purified in milligram yields and found to be homogeneous, monomeric, antibody-reactive, and stable at concentrations lower than 1 microM. Using similar procedures, analogous results were obtained with a murine receptor specific for an influenza nucleoprotein epitope (366-374) restricted by H2-Db. Production of these receptors has facilitated a detailed analysis of viral peptide-Major Histocompatibility Complex (peptide-MHC) engagement by the TCR using both surface plasmon resonance (SPR) and, in the case of the human TCR, isothermal titration calorimetry (ITC) (Willcox et al., 1999). The recombinant methods described should enable a wide range of TCR-peptide-MHC interactions to be studied and may also have implications for the production of other heterodimeric receptor molecules.     

Structural energetics of protein-carbohydrate interactions: Insights derived from the study of lysozyme binding to its natural saccharide inhibitors

High-sensitivity isothermal titration calorimetry was used to characterize the binding of the glycohydrolitic enzyme hen egg-white lysozyme to its natural saccharide inhibitors, chitobiose and chitrotriose. Measurements were done at a pH of 4.7, in the 15 degrees C -45 degrees C temperature range. Using a structural-energetic parameterization derived previously for lectin-carbohydrate associations, both binding enthalpies and entropies for the present systems and for the complex of chitobiose with turkey egg-white lysozyme from the literature were correctly accounted for. These observations suggest that both lysozymes and lectins follow the same structural-energetic behavior in the binding to their ligands. From the analysis of lysozyme data in conjunction with other binding data reported in the literature, an ad hoc parameterization of DeltaCp for protein-carbohydrate complexes was derived for the first time. The novel parameters for both polar and apolar surface areas differed significantly from correlations obtained previously from model compounds and protein-folding data. As DeltaCp is extremely sensitive to changes in solvent structure, this finding indicates that protein-carbohydrate complexes have distinctive hydration properties. According to our analysis, the dehydration of polar groups is the major cause for the observed decrease in DeltaCp, which implies that these groups behave hydrophobically. The contribution of apolar surface areas was found of the expected sign, but their specific weight is much smaller than those obtained in other correlations. This small contribution to DeltaCp is consistent with Lemieux's hypothesis of a low degree of hydration of apolar surfaces on carbohydrates.     

Kinetics and energetics of ligand binding determined by microcalorimetry: insights into active site mobility in a psychrophilic alpha-amylase

A new microcalorimetric method for recording the kinetic parameters k(cat), K(m) and K(i) of alpha-amylases using polysaccharides and oligosaccharides as substrates is described. This method is based on the heat released by glycosidic bond hydrolysis. The method has been developed to study the active site properties of the cold-active alpha-amylase produced by an Antarctic psychrophilic bacterium in comparison with its closest structural homolog from pig pancreas. It is shown that the psychrophilic alpha-amylase is more active on large macromolecular substrates and that the higher rate constants k(cat) are gained at the expense of a lower affinity for the substrate. The active site is able to accommodate larger inhibitory complexes, resulting in a mixed-type inhibition of starch hydrolysis by maltose. A method for recording the binding enthalpies by isothermal titration calorimetry in a low-affinity system has been developed, allowing analysis of the energetics of weak ligand binding using the allosteric activator chloride. It is shown that the low affinity of the psychrophilic alpha-amylase for chloride is entropically driven. The high enthalpic and entropic contributions of activator binding suggest large structural fluctuations between the free and the bound states of the cold-active enzyme. The kinetic and thermodynamic data for the psychrophilic alpha-amylase indicate that the strictly conserved side-chains involved in substrate binding and catalysis possess an improved mobility, responsible for activity in the cold, and resulting from the disappearance of stabilizing interactions far from the active site.     

Iron binding and oxidation kinetics in frataxin CyaY of Escherichia coli

Friedreich's ataxia is associated with a deficiency in frataxin, a conserved mitochondrial protein of unknown function. Here, we investigate the iron binding and oxidation chemistry of Escherichia coli frataxin (CyaY), a homologue of human frataxin, with the aim of better understanding the functional properties of this protein. Anaerobic isothermal titration calorimetry (ITC) demonstrates that at least two ferrous ions bind specifically but relatively weakly per CyaY monomer (K(d) approximately 4 microM). Such weak binding is consistent with the hypothesis that the protein functions as an iron chaperone. The bound Fe(II) is oxidized slowly by O(2). However, oxidation occurs rapidly and completely with H(2)O(2) through a non-enzymatic process with a stoichiometry of two Fe(II)/H(2)O(2), indicating complete reduction of H(2)O(2) to H(2)O. In accord with this stoichiometry, electron paramagnetic resonance (EPR) spin trapping experiments indicate that iron catalyzed production of hydroxyl radical from Fenton chemistry is greatly attenuated in the presence of CyaY. The Fe(III) produced from oxidation of Fe(II) by H(2)O(2) binds to the protein with a stoichiometry of six Fe(III)/CyaY monomer as independently measured by kinetic, UV-visible, fluorescence, iron analysis and pH-stat titrations. However, as many as 25-26 Fe(III)/monomer can bind to the protein, exhibiting UV absorption properties similar to those of hydrolyzed polynuclear Fe(III) species. Analytical ultracentrifugation measurements indicate that a tetramer is formed when Fe(II) is added anaerobically to the protein; multiple protein aggregates are formed upon oxidation of the bound Fe(II). The observed iron oxidation and binding properties of frataxin CyaY may afford the mitochondria protection against iron-induced oxidative damage.     

Studies of lysozyme binding to histamine as a ligand for hydrophobic charge induction chromatography

Histamine was immobilized on Sepharose CL‐6B (Sepharose) for use as a ligand of hydrophobic charge induction chromatography (HCIC) of proteins. Lysozyme adsorption onto Histamine‐Sepharose (HA‐S) was studied by adsorption equilibrium and calorimetry to uncover the thermodynamic mechanism of the protein binding. In both the experiments, the influence of salt (ammonium sulfate and sodium sulfate) was examined. Adsorption isotherms showed that HA‐S exhibited a high salt tolerance in lysozyme adsorption. This property was well explained by the combined contributions of hydrophobic interaction and aromatic stacking. The isotherms were well fitted to the Langmuir equation, and the equilibrium parameters for lysozyme adsorption were obtained. In addition, thermodynamic parameters (ΔH_ads, ΔS_ads, and ΔG_ads) for the adsorption were obtained by isothermal titration calorimetry by titrating lysozyme solutions into the adsorbent suspension. Furthermore, free histamine was titrated into lysozyme solution in the same salt‐buffers. Compared with the binding of lysozyme to free histamine, lysozyme adsorption onto HA‐S was characterized by a less favorable ΔG_ads and an unfavorable ΔS_ads because histamine was covalently attached to Sepharose via a three‐carbon‐chain spacer. Consequently, the immobilized histamine could only associate with the residues on the protein surface rather than those in the hydrophobic pocket, causing a less favorable orientation between histamine and lysozyme. Further comparison of thermodynamic parameters indicated that the unfavorable ΔS_ads was offset by a favorable ΔH_ads, thus exhibiting typical enthalpy‐entropy compensation. Moreover, thermodynamic analyses indicated the importance of the dehydration of lysozyme molecule and HA‐S during the adsorption and a substantial conformational change of the protein during adsorption. The results have provided clear insights into the adsorption mechanisms of lysozyme onto the new HCIC material. © 2009 American Institute of Chemical Engineers Biotechnol. Prog., 2010     

Exploring variation in binding of Protein A and Protein G to immunoglobulin type G by isothermal titration calorimetry

Bacterial Protein A (PrtA) and Protein G (PrtG) are widely used for affinity purification of antibodies. An understanding of how PrtA and PrtG bind to different isotypes of immunoglobulin type G (IgG) and to their corresponding Fc fragments is essential for the development of PrtA and PrtG mimetic ligands and for the establishment of generic processes for the purification of various antibodies. In this paper, the interactions between the two IgG-binding proteins and IgG of two different subclasses, IgG1 and IgG4, as well as their analogous Fc fragments have been studied by isothermal titration calorimetry. The results indicate that both protein ligands bind IgG and Fc fragments strongly with Ka values in the range of 10(7) -10(8) M(-1) and for both ligands, the interaction with both IgG isotypes is enthalpically driven though entropically unfavorable. Moreover, variation in the standard entropic and standard enthalpic contribution to binding between the two isotypes as well as between IgG and Fc fragment implies that the specific interaction with PrtA varies according to IgG isotype. In contrast to PrtA, PrtG bound to F(ab')(2) fragment with a Ka value of 5.1 × 10(5) M(-1) ; thus underscoring the usefulness of PrtA as a preferred ligand for generic antibody purification processes.     

Insights into binding of fatty acids by fatty acid binding proteins

Members of the phylogenetically related intracellular lipid binding protein (iLBP) are characterized by a highly conserved tertiary structure, but reveal distinct binding preferences with regard to ligand structure and conformation, when binding is assessed by the Lipidex method (removal of unbound ligand by hydrophobic polymer) or by isothermal titration calorimetry, a true equilibrium method. Subfamily proteins bind retinoids, subfamily II proteins bind bulky ligands, examples are intestinal bile acid binding protein (I-BABP) and liver fatty acid binding protein (L-FABP) which binds 2 ligand molecules, preferably monounsaturated and n-3 fatty acids. Subfamily III intestinal fatty acid binding protein (I-FABP) binds fatty acid in a bent conformation. The fatty acid bound by subfamily IV FABPs has a U-shaped conformation; here heart (H-) FABP preferably binds n-6, brain (B-) FABP n-3 fatty acids. The ADIFAB-method is a fluorescent test for fatty acid in equilibrium with iLBP and reveals some correlation of binding affinity to fatty acid solubility in the aqueous phase; these data are often at variance with those obtained by the other methods. Thus, in this review published binding data are critically discussed, taking into account on the one hand binding increments calculated for fatty acid double bonds on the basis of the solubility hypothesis, on the other hand the interpretation of calorimetric data on the basis of crystallographic and solution structures of iLBPs.     

Quantification of Ca2+ binding to melanin supports the hypothesis that melanosomes serve a functional role in regulating calcium homeostasis

Calcium regulation in melanocytes affects numerous biological pathways including protecting the redox balance in the cell and regulating the supply of substrate, l ‐tyrosine, for melanogenesis. The pigment contained in the melanocytes, melanin, has been implicated in maintaining calcium homeostasis in the cell and is known to be involved with calcium ion regulation in the inner ear. Herein, the association constant for Ca²⁺ binding to Sepia melanin is determined by isothermal titration calorimetry to be 3.3 (±0.2) × 10³/M. This value is comparable with other well‐established intracellular calcium‐binding proteins that serve to buffer calcium concentrations, lending further support to the hypothesis that melanosomes serve as intracellular mediators of calcium homeostasis in melanocytes. Using this binding constant and the data from a fluorescent Ca²⁺ displacement assay, the pK_a of the carboxyl group coordinated to Ca²⁺ is determined to be 3.1 ± 0.1.     

On the molecular basis of the high affinity binding of basic amino acids to LAOBP,a periplasmic binding protein from Salmonella typhimurium

The rational designing of binding abilities in proteins requires an understanding of the relationship between structure and thermodynamics. However, our knowledge of the molecular origin of high‐affinity binding of ligands to proteins is still limited; such is the case for l ‐lysine–l ‐arginine–l ‐ornithine periplasmic binding protein (LAOBP), a periplasmic binding protein from Salmonella typhimurium that binds to l ‐arginine, l ‐lysine, and l ‐ornithine with nanomolar affinity and to l ‐histidine with micromolar affinity. Structural studies indicate that ligand binding induces a large conformational change in LAOBP. In this work, we studied the thermodynamics of l ‐histidine and l ‐arginine binding to LAOBP by isothermal titration calorimetry. For both ligands, the affinity is enthalpically driven, with a binding ΔCp of ~?300 cal mol^?1 K^?1, most of which arises from the burial of protein nonpolar surfaces that accompanies the conformational change. Osmotic stress measurements revealed that several water molecules become sequestered upon complex formation. In addition, LAOBP prefers positively charged ligands in their side chain. An energetic analysis shows that the protein acquires a thermodynamically equivalent state with both ligands. The 1000‐fold higher affinity of LAOBP for l ‐arginine as compared with l ‐histidine is mainly of enthalpic origin and can be ascribed to the formation of an extra pair of hydrogen bonds. Periplasmic binding proteins have evolved diverse energetic strategies for ligand recognition. STM4351, another arginine binding protein from Salmonella, shows an entropy‐driven micromolar affinity toward l ‐arginine. In contrast, our data show that LAOBP achieves nanomolar affinity for the same ligand through enthalpy optimization. Copyright © 2015 John Wiley & Sons, Ltd.     

Alternative modes of binding of proteins with tandem SH2 domains

The issue of specificity in tyrosine kinase intracellular signaling mediated by src homology 2 (SH2) domains has great importance in the understanding how individual signals maintain their mutual exclusivity and affect downstream responses. Several proteins contain tandem SH2 domains that, on interacting with their ligand, provide a higher level of specificity than can be afforded by the interaction of a single SH2 domain. In this study, we focus on the comparison of two proteins ZAP70 and the p85 subunit of PI 3-kinase, which although distinctly different in function and general structure, possess tandem SH2 domains separated by a linker region and which bind to phosphorylated receptor molecules localized to the cell membrane. Binding studies using isothermal titration calorimetry show that these two proteins interact with peptides mimicking their physiological ligands in very different ways. In the case of the SH2 domains from ZAP70, they interact with a stoichiometry of unity, while p85 is able to make two distinct interactions, one with a stoichiometry of 1:1 and the other with two p85 molecules interacting with one receptor. The observation of two different modes of binding of p85 might be important in providing different cellular responses based on fluctuating intracellular concentration regimes of this protein. Thermodynamic data on both proteins suggest that a conformational change occurs on binding. On investigation of this structural change using a truncated form of p85 (including just the two SH2 domains and the inter-SH2 region), both NMR and circular dichroism spectroscopic studies failed to show significant changes in secondary structure. This suggests that any conformational change associated with binding is small and potentially limited to loop regions of the protein.     

Thermodynamics of zinc binding to hepatitis C virus NS3 protease: A folding by binding event

The hepatitis C virus (HCV) nonstructural protein 3 (NS3) protease is responsible for the processing of the non‐structural region of the viral precursor polyprotein in infected hepatic cells. HCV NS3 is a zinc‐dependent serine protease. The zinc ion, which is bound far away from the active site and considered to have a structural role, is essential for the structural integrity of the protein; furthermore, the ion is required for the hydrolytic activity. Consequently, the NS3 zinc binding site has been considered for a long time as a possible target for drug discovery. As a first step towards this goal, the energetics of the NS3‐zinc interaction and its effect on the NS3 conformation must be established and discussed. The thermodynamic characterization of zinc binding to NS3 protease by isothermal titration calorimetry and spectroscopy is presented here. Spectroscopic and calorimetric results suggest that a considerable conformational change in the protein is coupled to zinc binding. The energetics of the conformational change is comparable to that of the folding of a protein of similar size. Therefore, zinc binding to NS3 protease can be considered as a “folding by binding” event. Proteins 2009. © 2009 Wiley‐Liss, Inc.     

Effect of inorganic phosphate on FMN binding and loop flexibility in Desulfovibrio desulfuricans apo-flavodoxin

The complex between flavin mononucleotide (FMN) and apo-flavodoxin is dominated by isoalloxazine-stacking interactions and 5'-phosphate hydrogen bonds. We show here that FMN binding to Desulfovibrio desulfuricans apo-flavodoxin is faster and the affinity is higher in the presence of inorganic phosphate as compared to in its absence (I=110 mM, pH 7, 20 degrees C). The transition-state of complex formation was investigated by phi-value analysis using Trp60Ala and Tyr98Ala apo-flavodoxin variants. We find that Tyr98 is highly involved in the FMN/protein transition state independent of inorganic phosphate, whereas the participation of Trp60 is modulated by inorganic phosphate. The phi-value for Trp60 is higher without phosphate, implying that at this condition stronger binding of Trp60 is required in the transition state to assure successful complex formation. Consistent with the experimental data, all-atom molecular dynamic simulations reveal that the presence of an anion in the phosphate subsite restricts the mobility of the Trp60-containing loop in terms of both backbone and side-chain movements, but has no effect on the Tyr98-containing loop. The overall thermodynamic stability of apo-flavodoxin is higher in the presence of inorganic phosphate as compared to in its absence (I=110 mM, pH 7, 20 degrees C). Kinetic experiments reveal that the additional stability originates in slower unfolding. The combined experimental and computational observations demonstrate that phosphate has an ordering effect on the Trp60-containing loop, which positions Trp60 favorably for FMN binding and increases the barrier for protein unfolding.     

Exploring the binding dynamics of etoricoxib with human hemoglobin: A spectroscopic,calorimetric, and molecular modeling approach

Etoricoxib, widely used for the treatment of osteoarthritis, rheumatoid arthritis, ankylosing spondylitis, and related conditions has ample affinity to bind with globular proteins. Here, the molecular interaction between purified human hemoglobin (HHb), a major heme protein and etoricoxib, a cyclooxygenase-2 inhibitor was studied by various spectroscopic, calorimetric, and molecular modeling techniques. The binding affected hypochromic changes in the Soret band of hemoglobin (Hb) and induced remarkable quenching of the intrinsic fluorescence property of protein molecules. Synchronous fluorescence studies revealed alterations in tryptophan (Trp) and tyrosine (Tyr) microenvironments of HHb molecule in presence of etoricoxib. Flouremetric and isothermal titration calorimetric studies suggested two binding sites in HHb for etoricoxib at three different temperatures (298.15, 303.15, and 310.15 K). Negative values of Gibbs energy change (ΔG⁰) and enthalpy change (ΔH⁰) strongly suggest that it is spontaneous and exothermic reaction, mainly stabilized by hydrogen bonding as evidenced by sucrose binding assay. These findings support our in silico molecular docking study, which specified the binding site and the non-covalent interactions involved in the association. Moreover, the interaction impacts on structural integrity and functional aspects of HHb as confirmed by decreased α helicity, increased free iron release, increased rate of co-oxidation, and decreased rate of esterase activity. Overall, these studies conclude that etoricoxib leads to a remarkable alteration in the conformational aspects of binding to HHb.

Communicated by Ramaswamy H. Sarma.