A model for ligand binding to hexacoordinate hemoglobins

Hexacoordinate hemoglobins are heme proteins capable of reversible intramolecular coordination of the ligand binding site by an amino acid side chain from within the heme pocket. Examples of these proteins are found in many living organisms ranging from prokaryotes to humans. The nonsymbiotic hemoglobins (nsHbs) are a class of hexacoordinate heme proteins present in all plants. The nsHb from rice (rHb1) has been used as a model system to develop methods for determining rate constants characterizing binding and dissociation of the His residue responsible for hexacoordination. Measurement of these reactions exploits laser flash photolysis to initiate the reaction from the unligated, pentacoordinate form of the heme protein. A model for ligand binding is presented that incorporates the reaction following rapid mixing with the reaction starting from the pentacoordinate hemoglobin (Hb). This model is based on results indicating that ligand binding to hexacoordinate Hbs is not a simple combination of competing first order (hexacoordination) and second order (exogenous ligand binding) reactions. Ligand binding following rapid mixing is a multiphasic reaction displaying time courses ranging from milliseconds to minutes. The new model incorporates a "closed", slow reacting form of the protein that is not at rapid equilibrium with the reactive conformation. It is also demonstrated that formation of the closed protein species is not dependent on hexacoordination.     

A simple method for the determination of affinity and binding site concentration in receptor binding studies.

In ligand binding studies, it is often difficult to apply kinetic analyses because of an uncertainty in experimental data obtained at high ligand concentrations. Under such circumstances, Kd value (an index of the affinity) and the binding site concentration may be estimated more accurately from the binding of a fixed concentration of labelled ligand observed in the presence of various concentrations of the non-labelled ligand, if the fraction of both labelled and non-labelled ligand bound is small. When there is no cooperative effect of the ligand binding, the Kd value may be calculated by subtracting the concentration of the labelled drug from the concentration of the non-labelled drug to cause a 50% reduction of the saturable binding of the labelled drug. From above values, the binding site concentration may be calculated. The proposed method is capable of examining the cooperativity of the ligand binding, the labelled drug concentration and the specific radioactivity of the labelled drug and does not require large amounts of the labelled drug.     

Ligand binding to proteins: the binding landscape model.

Models of ligand binding are often based on four assumptions: (1) steric fit: that binding is determined mainly by shape complementarity; (2) native binding: that ligands mainly bind to native states; (3) locality: that ligands perturb protein structures mainly at the binding site; and (4) continuity: that small changes in ligand or protein structure lead to small changes in binding affinity. Using a generalization of the 2D HP lattice model, we study ligand binding and explore these assumptions. We first validate the model by showing that it reproduces typical binding behaviors. We observe ligand-induced denaturation, ANS and heme-like binding, and "lock-and-key" and "induced-fit" specific binding behaviors characterized by Michaelis-Menten or more cooperative types of binding isotherms. We then explore cases where the model predicts violations of the standard assumptions. For example, very different binding modes can result from two ligands of identical shape. Ligands can sometimes bind highly denatured states more tightly than native states and yet have Michaelis-Menten isotherms. Even low-population binding to denatured states can cause changes in global stability, hydrogen-exchange rates, and thermal B-factors, contrary to expectations, but in agreement with experiments. We conclude that ligand binding, similar to protein folding, may be better described in terms of energy landscapes than in terms of simpler mass-action models.     

Membrane cholesterol oxidation inhibits ligand binding function of hippocampal serotonin(1A) receptors

We have monitored the ligand binding function of the bovine hippocampal 5-HT(1A) receptor following treatment of native membranes with cholesterol oxidase. Cholesterol oxidase is a water soluble enzyme that acts on the membrane interface to catalyze the conversion of cholesterol to cholestenone. Oxidation of membrane cholesterol significantly inhibits the specific binding of the agonist and antagonist to 5-HT(1A) receptors. Fluorescence polarization measurements of membrane probes incorporated at different locations in the membrane revealed no appreciable effect on membrane order due to the oxidation of cholesterol to cholestenone. These results therefore suggest that the ligand binding function of the 5-HT(1A) receptor is a cholesterol-dependent phenomenon that is not related to the ability of cholesterol to modulate membrane order. Importantly, these results represent the first report on the effect of a cholesterol-modifying agent on the ligand binding function of this important neurotransmitter receptor.     

Mixed mode of ligand-DNA binding results in S-shaped binding curves

S-shaped binding curves often characterize interactions of ligands with nucleic acid molecules as analyzed by different physico-chemical and biophysical techniques. S-shaped experimental binding curves are usually interpreted as indicative of the positive cooperative interactions between the bound ligand molecules. This paper demonstrates that S-shaped binding curves may occur as a result of the "mixed mode" of DNA binding by the same ligand molecule. Mixed mode of the ligand-DNA binding can occur, for example, due to 1) isomerization or dimerization of the ligands in solution or on the DNA lattice, 2) their ability to intercalate the DNA and to bind it within the minor groove in different orientations. DNA-ligand complexes are characterized by the length of the ligand binding site on the DNA lattice (so-called "multiple-contact" model). We show here that if two or more complexes with different lengths of the ligand binding sites could be produced by the same ligand, the dependence of the concentration of the complex with the shorter length of binding site on the total concentration of ligand should be S-shaped. Our theoretical model is confirmed by comparison of the calculated and experimental CD binding curves for bis-netropsin binding to poly(dA-dT) poly(dA-dT). Bis-netropsin forms two types of DNA complexes due to its ability to interact with the DNA as monomers and trimers. Experimental S-shaped bis-netropsin-DNA binding curve is shown to be in good correlation with those calculated on the basis of our theoretical model. The present work provides new insight into the analysis of ligand-DNA binding curves.     

Mutational analysis of the ligand binding domain of the low density lipoprotein receptor

The ligand binding domain of the low density lipoprotein (LDL) receptor contains seven imperfect repeats of a 40-amino acid cysteine-rich sequence. Each repeat contains clustered negative charges that have been postulated as ligand-binding sites. The adjacent region of the protein, the growth factor homology region, contains three cysteine-rich repeats (A-C) whose sequence differs from those in the ligand binding domain. To dissect the contribution of these different cysteine-rich repeats to ligand binding, we used oligonucleotide-directed mutagenesis to alter expressible cDNAs for the human LDL receptor which were then introduced into monkey COS cells by transfection. We measured the ability of the mutant receptors to bind LDL, which contains a single protein ligand for the receptor (apoB-100), and beta-migrating very low density lipoprotein (beta-VLDL), which contains apoB-100 plus multiple copies of another ligand (apoE). The results show that repeat 1 is not required for binding of either ligand. Repeats 2 plus 3 and repeats 6 plus 7 are required for maximal binding of LDL, but not beta-VLDL. Repeat 5 is required for binding of both ligands. Repeat A in the growth factor homology region is required for binding of LDL, but not beta-VLDL. Repeat B is not required for ligand binding. These results support a model for the LDL receptor in which various repeats play additive roles in ligand binding, each repeat making a separate contribution to the binding event.     

Landscape of protein–small ligand binding modes

Elucidating the mechanisms of specific small‐molecule (ligand) recognition by proteins is a long‐standing conundrum. While the structures of these molecules, proteins and ligands, have been extensively studied, protein–ligand interactions, or binding modes, have not been comprehensively analyzed. Although methods for assessing similarities of binding site structures have been extensively developed, the methods for the computational treatment of binding modes have not been well established. Here, we developed a computational method for encoding the information about binding modes as graphs, and assessing their similarities. An all‐against‐all comparison of 20,040 protein–ligand complexes provided the landscape of the protein–ligand binding modes and its relationships with protein‐ and chemical spaces. While similar proteins in the same SCOP Family tend to bind relatively similar ligands with similar binding modes, the correlation between ligand and binding similarities was not very high (R² = 0.443). We found many pairs with novel relationships, in which two evolutionally distant proteins recognize dissimilar ligands by similar binding modes (757,474 pairs out of 200,790,780 pairs were categorized into this relationship, in our dataset). In addition, there were an abundance of pairs of homologous proteins binding to similar ligands with different binding modes (68,217 pairs). Our results showed that many interesting relationships between protein–ligand complexes are still hidden in the structure database, and our new method for assessing binding mode similarities is effective to find them.     

On the benzodiazepine binding pocket in GABAA receptors

Benzodiazepines are used for their sedative/hypnotic, anxiolytic, muscle relaxant, and anticonvulsive effects. They exert their actions through a specific high affinity binding site on the major inhibitory neurotransmitter receptor, the gamma-aminobutyric acid, type A (GABA(A)) receptor channel, where they act as positive allosteric modulators. To start to elucidate the relative positioning of benzodiazepine binding site ligands in their binding pocket, GABA(A) receptor residues thought to reside in the site were individually mutated to cysteine and combined with benzodiazepine analogs carrying substituents reactive to cysteine. Direct apposition of such reactive partners is expected to lead to an irreversible site-directed reaction. We describe here the covalent interaction of alpha(1)H101C with a reactive group attached to the C-7 position of diazepam. This interaction was studied at the level of radioactive ligand binding and at the functional level using electrophysiological methods. Covalent reaction occurs concomitantly with occupancy of the binding pocket. It stabilizes the receptor in its allosterically stimulated conformation. Covalent modification is not observed in wild type receptors or when using mutated alpha(1)H101C-containing receptors in combination with the reactive ligand pre-reacted with a sulfhydryl group, and the modification rate is reduced by the binding site ligand Ro15-1788. We present in addition evidence that gamma(2)Ala-79 is probably located in the access pathway of the ligand to its binding pocket.     

A general method of analysis of ligand-macromolecule equilibria using a spectroscopic signal from the ligand to monitor binding. Application to Escherichia coli single-strand binding protein-nucleic acid interactions

We describe a general method for the analysis of ligand-macromolecule binding equilibria for cases in which the interaction is monitored by a change in a signal originating from the ligand. This method allows the absolute determination of the average degree of ligand binding per macromolecule without any assumptions concerning the number of modes or states for ligand binding or the relationship between the fractional signal change and the fraction of bound ligand. Although this method is generally applicable to any type of signal, we discuss the details of the method as it applies to the analysis of binding data monitored by a change in fluorescence of a ligand upon binding to a nucleic acid. We apply the analysis to the equilibrium binding of Escherichia coli single-strand binding (SSB) protein to single-stranded nucleic acids, which is monitored by the quenching of the intrinsic tryptophan fluorescence of the SSB protein. With this method, one can quantitatively determine the relationship between the fractional signal change of the ligand and the fraction of bound ligand, LB/LT, and rigorously test whether the signal change is directly proportional to LB/LT. For E. coli SSB protein binding to single-stranded nucleic acids in its (SSB)65 binding mode [Lohman, T. M., & Overman, L. B. (1985) J. Biol. Chem. 260, 3594; Chrysogelos, S., & Griffith, J. (1982) Proc. Natl. Acad. Sci. U.S.A. 79, 5803], we show that the fractional quenching of the SSB fluorescence is equal to the fraction of bound SSB.     

A salt-bridge motif involved in ligand binding and large-scale domain motions of the maltose-binding protein

The uptake of nutrients is essential for the survival of bacterial cells. Many specialized systems have evolved, such as the maltose-dependent ABC transport system that transfers oligosaccharides through the cytoplasmic membrane. The maltose/maltodextrin-binding protein (MBP) serves as an initial high-affinity binding component in the periplasm that delivers the bound sugar into the cognate ABC transporter MalFGK(2). We have investigated the domain motions induced by the binding of the ligand maltotriose into the binding cleft using molecular dynamics simulations. We find that MBP is predominantly in the open state without ligand and in the closed state with ligand bound. Oligosaccharide binding induces a closure motion (30.0 degrees rotation), whereas ligand removal leads to domain opening (32.6 degrees rotation) around a well-defined hinge affecting key areas relevant for chemotaxis and transport. Our simulations suggest that a "hook-and-eye" motif is involved in the binding. A salt bridge between Glu-111 and Lys-15 forms that effectively locks the protein-ligand complex in a semiclosed conformation inhibiting any further opening and promoting complete closure. This previously unrecognized feature seems to secure the ligand in the binding site and keeps MBP in the closed conformation and suggests a role in the initial steps of substrate transport.     

A simple method for the determination of affinity and binding site concentration in receptor binding studies

In ligand binding studies, it is often difficult to apply kinetic analyses because of an uncertainty in experimental data obtained at high ligand concentrations. Under such circumstances, K_d value (an index of the affinity) and the binding site concentration may be estimated more accurately from the binding of a fixed concentration of labelled ligand observed in the presence of various concentrations of the non-labelled ligand, if the fraction of both labelled and non-labelled ligand bound is small. When there is no cooperative effect of the ligand binding, the K_d value may be calculated by subtracting the concentration of the labelled drug from the concentration of the non-labelled drug to cause a 50% reduction of the saturable binding of the labelled drug. From above values, the binding site concentration may be calculated. The proposed method is capable of examining the cooperativity of the ligand binding, the labelled drug concentration and the specific radioactivity of the labelled drug and does not require large amounts of the labelled drug.     

Ligand binding by TPR domains

Tetratricopeptide repeat (TPR) domains bind specific peptide ligands and are thought to mediate protein-protein interactions in a variety of biological systems. Here we compare peptide ligand-binding by several different TPR domains. We present specific examples that demonstrate that TPR domains typically undergo little or no structural rearrangement upon ligand binding. Our data suggest that, contrary to a recent proposal, coupled folding and binding is not the common mechanism of ligand recognition by TPR domains.     

A comparison of lac repressor binding to operator and to nonoperator DNA

We have compared the operator and nonoperator DNA binding activities of the lac repressor with respect to inactivation or inhibition by trypsin, heat, actinomycin, and isopropylthiogal-actoside. The two DNA binding activities were found to differ only in their sensitivity to the inducing ligand isopropylthiogal-actoside. Repressor binding to poly(dT-dT-dG)·poly(dC-dA-dA) was shown not to be affected by isopropylthiogalactoside.     

Stabilization of the FK506 binding protein by ligand binding.

Although the rotamase activity of the FK506 binding protein is inhibited by ligand binding, it is hypothesized that the ligand/protein complex itself may be responsible for the immunosuppressive effects of FK506. We have therefore examined the structure of the FK506 binding protein in the presence of an analog of FK506 (FK520) by a combination of fluorescence, CD, FTIR and calorimetry. While only small changes in the overall structure of the protein may be induced by ligand, a large change in thermal stability of the binding protein is observed.     

A membrane-embedded glutamate is required for ligand binding to the multidrug transporter EmrE

EmrE is an Escherichia coli multidrug transporter that confers resistance to a variety of toxins by removing them in exchange for hydrogen ions. The detergent-solubilized protein binds tetraphenylphosphonium (TPP(+)) with a K(D) of 10 nM. One mole of ligand is bound per approximately 3 mol of EmrE, suggesting that there is one binding site per trimer. The steep pH dependence of binding suggests that one or more residues, with an apparent pK of approximately 7.5, release protons prior to ligand binding. A conservative Asp replacement (E14D) at position 14 of the only membrane-embedded charged residue shows little transport activity, but binds TPP(+) at levels similar to those of the wild-type protein. The apparent pK of the Asp shifts to <5.0. The data are consistent with a mechanism requiring Glu14 for both substrate and proton recognition. We propose a model in which two of the three Glu14s in the postulated trimeric EmrE homooligomer deprotonate upon ligand binding. The ligand is released on the other face of the membrane after binding of protons to Glu14.     

Advantages of two- or polyvalent binding of a receptor to the corresponding ligand in comparison to univalent binding

The features of monovalent and bivalent binding of receptors (or antibodies) with a polyvalent ligand (or with an antigen) are considered. It is shown that the rigid connection of the binding sites of the receptor brings to high increase of binding affinity for the corresponding ligand, but only in case if its epitopes are fully complementary to both sites of the receptor binding. If not, then there is no advantage of the binding of bivalent receptor before univalent binding. If the binding sites of the receptor are connected by a flexible linker, then regardless of location of epitopes of the corresponding ligand there is the successful fastening of receptor and ligand. Exactly the connection by a flexible linker is used by Nature in most cases at constructing of polyvalent receptors.