A simple theoretical model for fluorescence polarization binding assay development

Fluorescence polarization is a screening technology that is radioactivity free, homogeneous, and ratiometric. The signal measured with this technology is a weighted value of free and bound ligand. As a consequence, saturation curves are accessible only after calculation of the corresponding concentrations of free and bound ligand. To make this technology more accessible to assay development, the authors propose a simple mathematical model that predicts fluorescence polarization values from ligand and receptor total concentrations, depending on the corresponding dissociation constant. This model was validated using data of Bodipy-NDP-alphaMSH binding to MC(5), obtained after either ligand saturation of a receptor preparation or, conversely, receptor saturation of a ligand solution. These experimental data were also used to calculate the actual concentration of free and bound ligand and receptor and to obtain pharmacological constants by Scatchard analysis. A general method is proposed, which facilitates the design of fluorescence polarization binding assays by relying on the representation of theoretical polarization values. This approach is illustrated by the application to 2 systems of very different affinities.     

Ligand binding to the voltage-gated Kv1.5 potassium channel in the open state--docking and computer simulations of a homology model

The binding of blockers to the human voltage-gated Kv1.5 potassium ion channel is investigated using a three-step procedure consisting of homology modeling, automated docking, and binding free energy calculations from molecular dynamics simulations, in combination with the linear interaction energy method. A reliable homology model of Kv1.5 is constructed using the recently published crystal structure of the Kv1.2 channel as a template. This model is expected to be significantly more accurate than earlier ones based on less similar templates. Using the three-dimensional homology model, a series of blockers with known affinities are docked into the cavity of the ion channel and their free energies of binding are calculated. The predicted binding free energies are in very good agreement with experimental data and the binding is predicted to be mainly achieved through nonpolar interactions, whereas the relatively small differences in the polar contribution determine the specificity. Apart from confirming the importance of residues V505, I508, V512, and V516 for ligand binding in the cavity, the results also show that A509 and P513 contribute significantly to the nonpolar binding interactions. Furthermore, we find that pharmacophore models based only on optimized free ligand conformations may not necessarily capture the geometric features of ligands bound to the channel cavity. The calculations herein give a detailed structural and energetic picture of blocker binding to Kv1.5 and this model should thus be useful for further ligand design efforts.     

Absolute binding free energy calculations using molecular dynamics simulations with restraining potentials

The absolute (standard) binding free energy of eight FK506-related ligands to FKBP12 is calculated using free energy perturbation molecular dynamics (FEP/MD) simulations with explicit solvent. A number of features are implemented to improve the accuracy and enhance the convergence of the calculations. First, the absolute binding free energy is decomposed into sequential steps during which the ligand-surrounding interactions as well as various biasing potentials restraining the translation, orientation, and conformation of the ligand are turned "on" and "off." Second, sampling of the ligand conformation is enforced by a restraining potential based on the root mean-square deviation relative to the bound state conformation. The effect of all the restraining potentials is rigorously unbiased, and it is shown explicitly that the final results are independent of all artificial restraints. Third, the repulsive and dispersive free energy contribution arising from the Lennard-Jones interactions of the ligand with its surrounding (protein and solvent) is calculated using the Weeks-Chandler-Andersen separation. This separation also improves convergence of the FEP/MD calculations. Fourth, to decrease the computational cost, only a small number of atoms in the vicinity of the binding site are simulated explicitly, while all the influence of the remaining atoms is incorporated implicitly using the generalized solvent boundary potential (GSBP) method. With GSBP, the size of the simulated FKBP12/ligand systems is significantly reduced, from approximately 25,000 to 2500. The computations are very efficient and the statistical error is small ( approximately 1 kcal/mol). The calculated binding free energies are generally in good agreement with available experimental data and previous calculations (within approximately 2 kcal/mol). The present results indicate that a strategy based on FEP/MD simulations of a reduced GSBP atomic model sampled with conformational, translational, and orientational restraining potentials can be computationally inexpensive and accurate.     

Binding of nonsteroidal anti‐inflammatory drugs to Aβ fibril

Nonsteroidal anti‐inflammatory drugs are considered as potential therapeutic agents against Alzheimer's disease. Using replica exchange molecular dynamics and atomistic implicit solvent model, we studied the mechanisms of binding of naproxen and ibuprofen to the Aβ fibril derived from solid‐state NMR measurements. The binding temperature of naproxen is found to be almost 40 K higher than of ibuprofen implicating higher binding affinity of naproxen. The key factor, which enhances naproxen binding, is strong interactions between ligands bound to the surface of the fibril. The naphthalene ring in naproxen appears to provide a dominant contribution to ligand‐ligand interactions. In contrast, ligand‐fibril interactions cannot explain differences in the binding affinities of naproxen and ibuprofen. The concave fibril edge with the groove is identified as the primary binding location for both ligands. We show that confinement of the ligands to the groove facilitates ligand‐ligand interactions that lowers the energy of the ligands bound to the concave edge compared with those bound to the convex edge. Our simulations appear to provide microscopic rationale for the differing binding affinities of naproxen and ibuprofen observed experimentally. Proteins 2010. © 2010 Wiley‐Liss, Inc.     

The use of transferred nuclear Overhauser effects in the study of the conformations of small molecules bound to proteins

A novel method is proposed for the study of the conformation in solution of small molecules bound to proteins. In transfer of saturation experiments, irradiation at the frequency of a proton in the bound ligand can result in an intensity change in the signal from a different proton in the free excess ligand via a nuclear Overhauser effect between the two protons in the bound ligand. Approximatel calculations show that the observation of such effects depends upon the close spatial proximity (within about 4.0 Å) of the two protons involved and thus gives useful conformational information. Two examples of this method are given, for the binding of trimethoprim and NADP⁺, respectively, to Lactobacillus casei dihydrofolate reductase.     

Conformation of 3'CMP bound to RNase A using TrNOESY

The conditions for accurately determining distance constraints from TrNOESY data on a small ligand (3'CMP) bound to a small protein (RNase A, <14 kDa) are described. For small proteins, normal TrNOESY conditions of 10:1 ligand:protein or greater can lead to inaccurate structures for the ligand-bound conformation due to the contribution of the free ligand to the TrNOESY signals. By using two ligand:protein ratios (2:1 and 5:1), which give the same distance constraints, a conformation of 3'CMP bound to RNase A was determined (glycosidic torsion angle, chi=-166 degrees ; pseudorotational phase angle, 0 degrees < or = P < or =36 degrees ). Ligand-protein NOESY cross peaks were also observed and used to dock 3'CMP into the binding pocket of the apo-protein (7rsa). After energy minimization, the conformation of the 3'CMP:RNase A complex was similar to the X-ray structure (1 rpf) except that a C3'-endo conformation for the ribose ring (rather than C2'-exo conformation) was found in the TrNOESY structure.     

Implications of the ligandin binding site on the binding of non-substrate ligands to Schistosoma japonicum-glutathione transferase

The binding interactions between dimeric glutathione transferase from Schistosoma japonicum (Sj26GST) and bromosulfophthalein (BS) or 8-anilino-1-naphthalene sulfonate (ANS) were characterised by fluorescence spectroscopy and isothermal titration calorimetry (ITC). Both ligands inhibit the enzymatic activity of Sj26GST in a non-competitive form. A stoichiometry of 1 molecule of ligand per mole of dimeric enzyme was obtained for the binding of these ligands. The affinity of BS is higher (K(d)=3.2 microM) than that for ANS (K(d)=195 microM). The thermodynamic parameters obtained by calorimetric titrations are pH-independent in the range of 5.5 to 7.5. The interaction process is enthalpically driven at all the studied temperatures. This enthalpic contribution is larger for the ANS anion than for BS. The strongly favourable enthalpic contribution for the binding of ANS to Sj26GST is compensated by a negative entropy change, due to enthalpy-entropy compensation. DeltaG degrees remains almost invariant over the temperature range studied. The free energy change for the binding of BS to Sj26GST is also favoured by entropic contributions at temperatures below 32 degrees C, thus indicating a strong hydrophobic interaction. Heat capacity change obtained for BS (DeltaC(p) degrees =(-580.3+/-54.2) cal x K(-1) mol(-1)) is twofold larger (in absolute value) than for ANS (DeltaC(p) degrees =(-294.8+/-15.8) cal x K(-1) mol(-1)). Taking together the thermodynamic parameters obtained for these inhibitors, it can be argued that the possible hydrophobic interactions in the binding of these inhibitors to L-site must be accompanied by other interactions whose contribution is enthalpic. Therefore, the non-substrate binding site (designed as ligandin) on Sj26GST may not be fully hydrophobic.     

Estrogen receptor-ligand complexes measured by chip-based nanoelectrospray mass spectrometry: an approach for the screening of endocrine disruptors

In the present report, a method based on chip-based nanoelectrospray mass spectrometry (nanoESI-MS) is described to detect noncovalent ligand binding to the human estrogen receptor alpha ligand-binding domain (hERalpha LBD). This system represents an important environmental interest, because a wide variety of molecules, known as endocrine disruptors, can bind to the estrogen receptor (ER) and induce adverse health effects in wildlife and humans. Using proper experimental conditions, the nanoESI-MS approach allowed for the detection of specific ligand interactions with hERalpha LBD. The relative gas-phase stability of selected hERalpha LBD-ligand complexes did not mirror the binding affinity in solution, a result that demonstrates the prominent role of hydrophobic contacts for stabilizing ER-ligand complexes in solution. The best approach to evaluate relative solution-binding affinity by nanoESI-MS was to perform competitive binding experiments with 17beta-estradiol (E2) used as a reference ligand. Among the ligands tested, the relative binding affinity for hERalpha LBD measured by nanoESI-MS was 4-hydroxtamoxifen approximately diethylstilbestrol > E2 > genistein > bisphenol A, consistent with the order of the binding affinities in solution. The limited reproducibility of the bound to free protein ratio measured by nanoESI-MS for this system only allowed the binding constants (K(d)) to be estimated (low nanomolar range for E2). The specificity of nanoESI-MS combined with its speed (1 min/ligand), low sample consumption (90 pmol protein/ligand), and its sensitivity for ligand (30 ng/mL) demonstrates that this technique is a promising method for screening suspected endocrine disrupting compounds and to qualitatively evaluate their binding affinity.     

How maltose influences structural changes to bind to maltose‐binding protein: Results from umbrella sampling simulation

A well‐studied periplasmic‐binding protein involved in the abstraction of maltose is maltose‐binding protein (MBP), which undergoes a ligand‐induced conformational transition from an open (ligand‐free) to a closed (ligand‐bound) state. Umbrella sampling simulations have been us to estimate the free energy of binding of maltose to MBP and to trace the potential of mean force of the unbinding event using the center‐of‐mass distance between the protein and ligand as the reaction coordinate. The free energy thus obtained compares nicely with the experimentally measured value justifying our theoretical basis. Measurement of the domain angle (N‐terminal‐domain – hinge – C‐terminal‐domain) along the unbinding pathway established the existence of three different states. Starting from a closed state, the protein shifts to an open conformation during the initial unbinding event of the ligand then resides in a semi‐open conformation and later resides predominantly in an open‐state. These transitions along the ligand unbinding pathway have been captured in greater depth using principal component analysis. It is proposed that in mixed‐model, both conformational selection and an induced‐fit mechanism combine to the ligand recognition process in MBP. Proteins 2013. © 2012 Wiley Periodicals, Inc.     

An efficient computational method for calculating ligand binding affinities

Virtual compound screening using molecular docking is widely used in the discovery of new lead compounds for drug design. However, the docking scores are not sufficiently precise to represent the protein-ligand binding affinity. Here, we developed an efficient computational method for calculating protein-ligand binding affinity, which is based on molecular mechanics generalized Born/surface area (MM-GBSA) calculations and Jarzynski identity. Jarzynski identity is an exact relation between free energy differences and the work done through non-equilibrium process, and MM-GBSA is a semimacroscopic approach to calculate the potential energy. To calculate the work distribution when a ligand is pulled out of its binding site, multiple protein-ligand conformations are randomly generated as an alternative to performing an explicit single-molecule pulling simulation. We assessed the new method, multiple random conformation/MM-GBSA (MRC-MMGBSA), by evaluating ligand-binding affinities (scores) for four target proteins, and comparing these scores with experimental data. The calculated scores were qualitatively in good agreement with the experimental binding affinities, and the optimal docking structure could be determined by ranking the scores of the multiple docking poses obtained by the molecular docking process. Furthermore, the scores showed a strong linear response to experimental binding free energies, so that the free energy difference of the ligand binding (ΔΔG) could be calculated by linear scaling of the scores. The error of calculated ΔΔG was within ≈±1.5 kcal•mol⁻¹ of the experimental values. Particularly, in the case of flexible target proteins, the MRC-MMGBSA scores were more effective in ranking ligands than those generated by the MM-GBSA method using a single protein-ligand conformation. The results suggest that, owing to its lower computational costs and greater accuracy, the MRC-MMGBSA offers efficient means to rank the ligands, in the post-docking process, according to their binding affinities, and to compare these directly with the experimental values.     

Overall charge and local charge density of pectin determines the enthalpic and entropic contributions to complexation with β-lactoglobulin

The complex formation between β-lactoglobulin and pectins of varying overall charge and local charge density were investigated. Isothermal titration calorimetry experiments were carried out to determine the enthalpic contribution to the complex formation at pH 4.25 and various ionic strengths. Complex formation was found to be an exothermic process for all conditions. Combination with previously published binding constants by Sperber et al. (Sperber, B. L. H. M.; Cohen Stuart, M. A.; Schols, H. A.; Voragen, A. G. J.; Norde, W. Biomacromolecules 2009, 10, 3246-3252) allows for the determination of the changes in the Gibbs energy and the change in entropy of the system upon complex formation between β-lactoglobulin and pectin. The local charge density of pectin is found to determine the balance between enthalpic and entropic contributions. For a high local charge density pectin, the main contribution to the Gibbs energy is of an enthalpic nature, supported by a favorable entropy effect due to the release of small counterions. A pectin with a low local charge density has a more even distribution of the enthalpic and entropic part to the change of the Gibbs energy. The enthalpic part is reduced due to the lower charge density, while the relative increase of the entropic contribution is thought to be caused by a change in the location of the binding place for pectin on the β-lactoglobulin molecule. The association of the hydrophobic methyl esters on pectin with an exposed hydrophobic region on β-lg results in the release of water molecules from the hydrophobic region and surrounding the methyl esters of the pectin molecule. An increase in the ionic strength decreases the enthalpic contribution due to the shielding of electrostatic attraction in favor of the entropic contribution, supporting the idea that the release of water molecules from hydrophobic areas plays a part in the complex formation.     

Contribution of protein conformation to stereochemistry and reactivity of the active center of heme proteins and enzymes. The existence of horseradish peroxidase conformations and their possible role in the catalysis mechanism

In the spectral region 350-800 nm at 4.2 K we measured magnetic circular dichroism (MCD) spectra of the pentacoordinated complex of protcheme with 2-methylimidazole, deoxyleghemoglobin, neutral and alkaline forms of reduced horseradish peroxidase in the equilibrium states, as well as in non-equilibrium states produced by low-temperature photolysis of their carbon monoxide derivatives. Earlier the corresponding results have been obtained for myoglobin, hemoglobin and cytochromes P-450 and P-420. The energies of Fe-N (proximal His) and Fe-N(pyrroles) bonds and their changes upon ligand binding in heme proteins and enzymes were compared with those in the model heme complex thus providing conformational contribution into stereochemistry of the active site. The examples of weak and strong conformational "pressure" on stereochemistry were analysed and observed. If conformational energy contribution into stereochemistry prevails the electronic one the heme stereochemistry remains unchanged on ligand binding as it was observed for leghemoglobin and alkaline horseradish peroxidase. The change of bond energies in myoglobin and hemoglobin on ligand binding are comparable with those in protein free pentacoordinated protoheme, giving an example of weak conformational contribution to heme stereochemistry. The role of protein conformation energy in the modulation of ligand binding properties of heme in leghemoglobin relative to those in myoglobins is discussed. The most striking result were obtained in the study of reduced horseradish peroxidase in the pH region of 6.0-10.2. It was found that such different perturbations as ligand binding and heme-linked ionization of the distal amino acid residue induce identical changes in heme stereochemistry. Neither heme-linked ionization in the carbon monoxide complex nor the geometry of Fe-Co bond affect the heme local structure of photoproducts. These and other findings suggest a very low conformation mobility of horseradish peroxidase whose protein constraints appear to allow only two preferable geometries of specific amino acid residues that form the heme pocket. The role of the two tertiary structure constraints on the heme in the mechanism of horseradish peroxidase function is discussed. It is supposed that one conformation produces a heme environment suitable for two-electron oxidation of the native enzyme to compound I by hydrogen peroxide while another conformation changes the heme stereochemistry in the direction favourable for back reduction of compound I by the substrate to the resting enzyme through two one-electron steps. The switch from one tertiary structure to another is expected to be induced by substrate bind     

Energetics of ligand-induced conformational flexibility in the lactose permease of Escherichia coli

Isothermal titration calorimetry has been applied to characterize the thermodynamics of ligand binding to wild-type lactose permease (LacY) and a mutant (C154G) that strongly favors an inward facing conformation. The affinity of wild-type or mutant LacY for ligand and the change in free energy (DeltaG) upon binding are similar. However, with the wild type, the change in free energy upon binding is due primarily to an increase in the entropic free energy component (TDeltaS), whereas in marked contrast, an increase in enthalpy (DeltaH) is responsible for DeltaG in the mutant. Thus, wild-type LacY behaves as if there are multiple ligand-bound conformational states, whereas the mutant is severely restricted. The findings also indicate that the structure of the mutant represents a conformational intermediate in the overall transport cycle.     

DNA bis-intercalation: Theoretical calculation of binding curves

A statistical mechanical calculation of the binding properties of DNA bis-intercalators is presented, based on the sequence-generating function method of Lifson. The effects of binding by intercalation of one or both chromophores of a bifunctional intercalating agent are examined. The secular equation for a general model that includes the effects of neighbor (nearest and non-nearest) exclusion and/or cooperativity in the binding of both singly and doubly intercalated ligands is derived. Numerical results for binding curves are presented for a more restricted model in which each type of bound ligand rigorously excludes its nearest neighbor and the total number of sites covered by a doubly intercalated ligand is variable. At low values of free ligand concentration bis-intercalation dominates the binding process, while at high value of free ligand concentration, intercalation of only one chromophore per ligand becomes significant due to the unavailability of contiguous free sites required for bis-intercalation. Also, depending on the binding parameters, the free energy of the system can be lowered by a loss of doubly intercalated ligands in favor of singly intercalated ones. Corresponding to this transition in binding mode, the average number of sites occupied by a bound ligand decreases from that characteristic of bis-intercalation to that characteristic of mono-intercalation as free ligand concentration increases. An analysis of Scatchard plots describing bis-intercalation is presented.     

Evidence that membrane insertion of the cytosolic domain of Bcl-xL is governed by an electrostatic mechanism

Signals from different cellular networks are integrated at the mitochondria in the regulation of apoptosis. This integration is controlled by the Bcl-2 proteins, many of which change localization from the cytosol to the mitochondrial outer membrane in this regulation. For Bcl-xL, this change in localization reflects the ability to undergo a conformational change from a solution to integral membrane conformation. To characterize this conformational change, structural and thermodynamic measurements were performed in the absence and presence of lipid vesicles with Bcl-xL. A pH-dependent model is proposed for the solution to membrane conformational change that consists of three stable conformations: a solution conformation, a conformation similar to the solution conformation but anchored to the membrane by its C-terminal transmembrane domain, and a membrane conformation that is fully associated with the membrane. This model predicts that the solution to membrane conformational change is independent of the C-terminal transmembrane domain, which is experimentally demonstrated. The conformational change is associated with changes in secondary and, especially, tertiary structure of the protein, as measured by far and near-UV circular dichroism spectroscopy, respectively. Membrane insertion was distinguished from peripheral association with the membrane by quenching of intrinsic tryptophan fluorescence by acrylamide and brominated lipids. For the cytosolic domain, the free energy of insertion (DeltaG degrees x) into lipid vesicles was determined to be -6.5 kcal mol(-1) at pH 4.9 by vesicle binding experiments. To test whether electrostatic interactions were significant to this process, the salt dependence of this conformational change was measured and analyzed in terms of Gouy-Chapman theory to estimate an electrostatic contribution of DeltaG degrees el approximately -2.5 kcal mol(-1) and a non-electrostatic contribution of DeltaG degrees nel approximately -4.0 kcal mol(-1) to the free energy of insertion, DeltaG degrees x. Calcium, which blocks ion channel activity of Bcl-xL, did not affect the solution to membrane conformational change more than predicted by these electrostatic considerations. The lipid cardiolipin, that is enriched at mitochondrial contact sites and reported to be important for the localization of Bcl-2 proteins, did not affect the solution to membrane conformational change of the cytosolic domain, suggesting that this lipid is not involved in the localization of Bcl-xL in vivo. Collectively, these data suggest the solution to membrane conformational change is controlled by an electrostatic mechanism. Given the distinct biological activities of these conformations, the possibility that this conformational change might be a regulatory checkpoint for apoptosis is discussed.     

Energetics of subunit assembly and ligand binding in human hemoglobin.

An extensive and self-consistent set of thermodynamic properties has recently been established for the coupled processes of subunit assembly and ligand binding (oxygen and protons) in human hemoglobin. The resulting thermodynamic values permit a consideration of the possible sources of energetic terms accounting for stability of the tetrameric quaternary structures at different stages of ligation, and of the possible sources of cooperative energy. The analysis indicates that: (a) The change in buried surface ara upon oxygenation (i.e., hydrophobic stabilization) does not play a dominant role in stabilizing the unliganded tetramer relative to the liganded tetramer. (b) The pattern of enthalpic and entropic contributions to the free energies of dimer-tetramer. (c) The thermodynamic results are consistent with a dominant role of increased hydrogen bond formation in the deoxy quaternary structure. (d) Within tetramers the variation in free energy for successive oxygenation steps arises from both enthalpic and entropic contributions and the enthalpic contributions are almost entirely attributable to the heats of Bohr proton release. At pH 7.4 the pattern of thermodynamic values suggests that a large contribution to the free energy of cooperativity may arise from the energetics of Bohr proton release. It is suggested that a combination of proton ionization and hydrogen bonding may account for the main energetic features of cooperativity. Possible contributions from fluctuation behavior cannot presently be evaluated.     

Insights from the energetics of water binding at the domain-ligand interface of the Src SH2 domain

SH2 domains play important roles in signal transduction by binding phosphorylated tyrosine residues on cell surface receptors. In an effort to understand the mechanism of ligand binding and more specifically the role of water, we have designed a general computational protocol based on the potential of mean force to compute the thermodynamics of water molecules at the protein-ligand interface for two SH2 domain complexes of the Src kinase, those bound to the two peptides Ac-PQpYEpYI-NH2 and Ac-PQpYIpYV-NH2 where pY indicates a phosphotyrosine. These two peptides were chosen because they have similar binding affinities but very different entropic/enthalpic thermodynamic binding signatures, indicating different interactions with solvent. We find that the isoleucine to valine mutation at position +3 (the third amino acid C-terminal to pY) in the ligand has only limited impact on the water structure. By contrast, the glutamic acid to isoleucine mutation at position +1 has a significant impact by not only abrogating a local hydrophilic binding site but, more importantly and surprisingly, inducing a favorable nonlocal entropic contribution from the water molecules around the phosphorylated tyrosine at the +2 position. Our study demonstrates the validity of the method reported here for exploring the thermodynamic solvation landscape of protein-protein interactions.