Protein-ligand binding with a missing species

Considerable experimental evidence has been produced recently that shows that in the binding of oxygen or carbon monoxide to certain tetrameric hemoglobins, the triply-ligated species is virtually non-existent. The binding polynomial representing this phenomenon for the general case is P(x) = 1 + beta 1x + ... + beta n-1xn-1 + beta nxn, where beta n-1 is nearly zero. The zeros, factorization and associated Hill plots of such binding polynomials with beta n-1 = 0 are investigated for the general case, and are analyzed in detail for n = 3 and n = 4. These results are then compared with the results obtained from experimental data on a number of tetrameric hemoglobins for which beta 3 is small. One concludes that, apart from the slope of the high-saturation asymptote of the Hill plot, a small perturbation of beta 3 from zero produces small changes in other properties associated with the binding process, such as fractional saturation, maximum Hill slope, and zeros and factorization of the binding polynomial.     

A probabilistic model for fitting MWC polynomials in protein-ligand binding

Given a binding polynomial in Adair form, A(x) = 1 + beta 1 x + ... + beta n x n, beta i greater than or equal to 0, a basic problem is to determine a method of fitting a model polynomial to A(x) and a quantitative measure of the goodness of fit. This paper presents such a method for fitting Monod-Wyman-Changeux (MWC) model polynomials when A(x) is of degree three or four. The method of fitting is based on the property that the zeros of an MWC polynomial of any degree lie on a circle in the complex plane. The parameters in the MWC model are determined so that if possible this circle coincides with the circle on which lie the zeros of A(x). The measure of goodness of fit is provided by a probabilistic model which gives the probability that a binding polynomial has its zeros on a circle on which lie the zeros of an MWC polynomial and if so, the probability that the juxtaposition of the two sets of zeros can occur by chance alone.     

Location of zeros of Wiener and distance polynomials

The geometry of polynomials explores geometrical relationships between the zeros and the coefficients of a polynomial. A classical problem in this theory is to locate the zeros of a given polynomial by determining disks in the complex plane in which all its zeros are situated. In this paper, we infer bounds for general polynomials and apply classical and new results to graph polynomials namely Wiener and distance polynomials whose zeros have not been yet investigated. Also, we examine the quality of such bounds by considering four graph classes and interpret the results.     

A new measure of cooperativity in protein-ligand binding

An allosteric binding system consisting of a single ligand and a nondissociating macromolecule having multiple binding sites can be represented by a binding polynomial. Various properties of the binding process can be obtained by analyzing the coefficients of the binding polynomial and such functions as the binding curve and the Hill plot. The Hill plot has an asymptote of unit slope at each end and the departure of the slope from unity at any point can be used to measure the effective interaction free energy at that point. Of particular interest in detecting and measuring cooperativity are extrema of the Hill slope and its value at the half-saturation point. If the binding polynomial is symmetric, then there is an extremum of the Hill slope at the half-saturation point. This value, the Hill coefficient, is a convenient measure of cooperativity. The purpose of this paper is to express the Hill coefficient for symmetric binding polynomials in terms of the roots of the polynomial and to give an interpretation of cooperativity in terms of the geometric pattern of the roots in the complex plane. This interpretation is then applied to the binding polynomials for the MWC (Monod-Wyman-Changeux) and KNF (Koshland-Nemethy-Filmer) models.     

An analogue of the binding polynomial for the case of ligands binding to an aggregating macromolecule

It is known that when a macromolecule exists as a mixture of forms with different molecular weights in equilibrium in dilute solution with ligands then the saturation functions cannot be derived by differentiating a binding polynomial. However, we demonstrate that in such circumstances a saturation function can be expressed in terms of the binding polynomials for each of the polymeric forms of macromolecule and integration of this using the conservation equation affords a function, Q, which is the analogue of the binding polynomial for a non-aggregating system. The existence of Q leads to heterotropic linkage relationships.     

Ligand binding distributions in nucleic acids

We illustrate a new method for the determination of the complete binding polynomial for nucleic acids based on experimental titration data with respect to ligand concentration. From the binding polynomial, one can then calculate the distribution function for the number of ligands bound at any ligand concentration. The method is based on the use of a finite set of moments of the binding distribution function, which are obtained from the titration curve. Using the maximum-entropy method, the moments are then used to construct good approximations to the binding distribution function. Given the distribution functions at different ligand concentrations, one can calculate all of the coefficients in the binding polynomial no matter how many binding sites a molecule has. Knowledge of the complete binding polynomial in turn yields the thermodynamics of binding. This method gives all of the information that can be obtained from binding isotherms without the assumption of any specific molecular model for the nature of the binding. Examples are given for the binding of Mn(2+) and Mg(2+) to t-RNA and for the binding of Mg(2+) and I(6) to poly-C using literature data.     

Free energy of proton binding in proteins

In this article we use literature data on the titration of denatured ribonuclease to test the accuracy of proton-binding distributions obtained using our recent approach employing moments. We find that using only the local slope of the titration curve at a small number of points (five, for example) we can reproduce the detailed proton-binding distribution at all pH values. Our method gives the complete proton-binding polynomial for a given protein and each coefficient in this polynomial in turn yields the free energy for binding a given number of protons in all ways to the protein. Using these net free energies, we can then compute the average proton-binding free energy per proton as a function of the fraction of protons bound. We find that this function is remarkably similar for different proteins, even for proteins that exhibit quite different titration behavior. For the special case of binding to independent sites, we obtain simple relations for the first and last terms in the free energy per-proton function. For this special case we also can calculate the distribution functions giving the probability that a molecule has a given number of positive or negative charges and the joint distribution that a molecule simultaneously has a given number of positive and negative charge.     

Effects of heterogeneity and cooperativity on the forms of binding curves for multivalent ligands

An exploratory investigation is made of the binding behavior that is likely to be encountered with multivalent ligands under circumstances where a single intrinsic binding constant does not suffice to describe all acceptor-ligand interactions. Numerical simulations of theoretical binding behavior have established that current criteria for recognizing heterogeneity and cooperativity of acceptor sites on the basis of the deviation of the binding curve from rectangular hyperbolic form for univalent ligands also apply to the interpretation of the corresponding binding curves for multivalent ligands. However, for systems in which the source of the departure from equivalence and independence of binding sites resides in the ligand, these criteria are reversed. On the basis of these observations a case is then made for attributing results of an experimental binding study of the interaction between pyruvate kinase and muscle myofibrils to positive cooperativity of enzyme sites rather than to heterogeneity or negative cooperativity of the myofibrillar sites.     

Uncovering Community Structures with Initialized Bayesian Nonnegative Matrix Factorization

Uncovering community structures is important for understanding networks. Currently, several nonnegative matrix factorization algorithms have been proposed for discovering community structure in complex networks. However, these algorithms exhibit some drawbacks, such as unstable results and inefficient running times. In view of the problems, a novel approach that utilizes an initialized Bayesian nonnegative matrix factorization model for determining community membership is proposed. First, based on singular value decomposition, we obtain simple initialized matrix factorizations from approximate decompositions of the complex network’s adjacency matrix. Then, within a few iterations, the final matrix factorizations are achieved by the Bayesian nonnegative matrix factorization method with the initialized matrix factorizations. Thus, the network’s community structure can be determined by judging the classification of nodes with a final matrix factor. Experimental results show that the proposed method is highly accurate and offers competitive performance to that of the state-of-the-art methods even though it is not designed for the purpose of modularity maximization.     

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.     

Description of the topological entanglement of DNA catenanes and knots by a powerful method involving strand passage and recombination

We utilize a recently discovered, powerful method to classify the topological state of knots and catenanes. In this method, each such form is associated with a unique polynomial. These polynomials allow a rigorous determination of whether knotted or catenated DNA molecules that appear distinct actually are, and indicate the structure of related molecules. A tabulation is given of the polynomials for all possible stereoisomers of many of the knotted and catenated forms that are found in DNA. The polynomials for a substrate DNA molecule and the products obtained from it by either recombination or strand passage by a topoisomerase are related by a simple theorem. This theorem affords natural applications of the polynomial method to these processes. Examples are presented involving site-specific recombination by the transposon Tn3-encoded resolvase and the phage lambda integrase, in which product structure is predicted as a function of crossover mechanism.     

Protein-Binding Polynomials

Using a new method recently published for analyzing the binding isotherms of biopolymers (Poland, 2000a), we calculate the complete binding polynomials for lysozyme, insulin, and serum albumin from published titration data. These three proteins have, respectively, 22, 32, and 184 dissociable protons and hence are represented by series in powers of the hydrogen ion concentration with the highest powers in the series being the numbers just indicated. Given the complete binding polynomial, the distribution function giving the concentration of all states of proton binding can then be calculated at any pH.     

Protein-binding polynomials

Using a new method recently published for analyzing the binding isotherms of biopolymers (Poland, 2000a), we calculate the complete binding polynomials for lysozyme, insulin, and serum albumin from published titration data. These three proteins have, respectively, 22, 32, and 184 dissociable protons and hence are represented by series in powers of the hydrogen ion concentration with the highest powers in the series being the numbers just indicated. Given the complete binding polynomial, the distribution function giving the concentration of all states of proton binding can then be calculated at any pH.     

Ligand-induced conformational changes: improved predictions of ligand binding conformations and affinities

A computational docking strategy using multiple conformations of the target protein is discussed and evaluated. A series of low molecular weight, competitive, nonpeptide protein tyrosine phosphatase inhibitors are considered for which the x-ray crystallographic structures in complex with protein tyrosine phosphatase 1B (PTP1B) are known. To obtain a quantitative measure of the impact of conformational changes induced by the inhibitors, these were docked to the active site region of various structures of PTP1B using the docking program FlexX. Firstly, the inhibitors were docked to a PTP1B crystal structure cocrystallized with a hexapeptide. The estimated binding energies for various docking modes as well as the RMS differences between the docked compounds and the crystallographic structure were calculated. In this scenario the estimated binding energies were not predictive inasmuch as docking modes with low estimated binding energies corresponded to relatively large RMS differences when aligned with the corresponding crystal structure. Secondly, the inhibitors were docked to their parent protein structures in which they were cocrystallized. In this case, there was a good correlation between low predicted binding energy and a correct docking mode. Thirdly, to improve the predictability of the docking procedure in the general case, where only a single target protein structure is known, we evaluate an approach which takes possible protein side-chain conformational changes into account. Here, side chains exposed to the active site were considered in their allowed rotamer conformations and protein models containing all possible combinations of side-chain rotamers were generated. To evaluate which of these modeled active sites is the most likely binding site conformation for a certain inhibitor, the inhibitors were docked against all active site models. The receptor rotamer model corresponding to the lowest estimated binding energy is taken as the top candidate. Using this protocol, correct inhibitor binding modes could successfully be discriminated from proposed incorrect binding modes. Moreover, the ranking of the estimated ligand binding energies was in good agreement with experimentally observed binding affinities.     

Osmometric analysis of quasi-indefinite mixed associations of the first kind

Mixed associations of the type A + B----AB, A + AB----A2B, ..., A + Ai-1 B----AiB, ... are readily analyzed by osmometric methods. The equilibrium molar concentration of A, mA, is obtained very simply from mA = meq-m0B; here meq = c/Meqn is the equilibrium molar concentration of all associating species and m0B denotes the stoichiometric or original molar concentration of B. The quantity mB can then be obtained from methods developed by Steiner. The value of the binding polynomial lambda is given by lambda = m0B/mB; lambda is a function of mA only. In principle, one can evaluate the equilibrium constants (kA,B,etc.) by fitting lambda to the appropriate polynomial in mA of degree n (n = 2, 3, ...). The binding polynomial lambda is analogous to polynomials encountered in the analysis of self-associations. By making some simple assumptions one can develop four analogs of two sequential, equal equilibrium constant (SEK) or two attenuated equilibrium constant (AK) models. With the aid of r (the number average degree of binding), g (the osmotic coefficient), lambda, as well as mA and mB, one can evaluate the equilibrium constant or constants. The methods developed here can be extended to the nonideal case.     

Determination of binding constants for cooperative site-specific protein-DNA interactions using the gel mobility-shift assay

We have investigated the question of whether the gel mobility-shift assay can provide data that are useful to the demonstration of cooperativity in the site-specific binding of proteins to DNA. Three common patterns of protein-DNA interaction were considered: (i) the cooperative binding of a protein to two sites (illustrated by the Escherichia coli Gal repressor); (ii) the cooperative binding of a bidentate protein to two sites (illustrated by the E. coli Lac repressor); and (iii) the cooperative binding of a protein to three sites (illustrated by the lambda cI repressor). A simple, rigorous, and easily extendable statistical mechanical approach to the derivation of the binding equations for the different patterns is presented. Both simulated and experimental data for each case are analyzed. The mobility-shift assay provides estimates of the macroscopic binding constants for each step of ligation based on its separation of liganded species by the number of ligands bound. Resolution of the binding constants depends on the precision with which the equilibrium distribution of liganded species is determined over the entire range of titration of each of the sites. However, the evaluation of cooperativity from the macroscopic binding constants is meaningful only for data that are also accurate. Some criteria that are useful in evaluating accuracy are introduced and illustrated. Resolution of cooperative effects is robust only for the simplest case, in which there are two identical protein binding sites. In this case, cooperative effects of up to 1,000-fold are precisely determined. For heterogeneous sites, cooperative effects of greater than 1,000-fold are resolvable, but weak cooperativity is masked by the heterogeneity. For three-site systems, only averaged pair-wise cooperative effects are resolvable.     

Conditions for the various types of co-operativity allowed by the adair,equation and their relation with the algebraic character of binding polynomials

A convenient way to obtain for any number, n, of sites, the functions of the constants of the Adair equation that decide the type of co-operativity of ligand binding to a non-dissociating protein is given and is illustrated by the examples n = 4 and n = 5. These functions are invariants of the binding polynomial and various of its derivatives.Although there are some simple sufficient conditions (inequalities relating successive Adair constants) for some co-operativity types, the full necessary and sufficient conditions even for uniform positive and negative co-operativity depend on very complicated functions of the constants for n > 4.However there are alternative ways of writing binding polynomials known as canonical forms. Up to at least n = 5, and probably beyond, the conditions that are complicated in terms of Adair constants are very simple in terms of the constants of canonical forms. For instance any fourth-degree polynomial can be written in the form p(x - α)⁴ + q(x - β)⁴ + 6μ (x - α)²(x - β)² although in three different ways. For one of these ways, the sign of μ distinguishes between mixed and uniform co-operativity. For any kind of mixed co-operativity μ > 0, while μ < 0 corresponds to uniform co-operativity. Advantages of the use of canonical forms are briefly commented on.