Atomistic simulations of competition between substrates binding to an enzyme

Although the idea that electrostatic potentials generated by enzymes can guide substrates to active sites is well established, it is not always appreciated that the same potentials can also promote the binding of molecules other than the intended substrate, with the result that such enzymes might be sensitive to the presence of competing molecules. To provide a novel means of studying such "electrostatic competition" effects, computer simulation methodology has been developed to allow the diffusion and association of many solute molecules around a single enzyme to be simulated. To demonstrate the power of the methodology, simulations have been conducted on an artificial fusion protein of citrate synthase (CS) and malate dehydrogenase (MDH) to assess the chances of oxaloacetate being channeled between the MDH and CS active sites. The simulations demonstrate that the probability of channeling is strongly dependent on the concentration of the initial substrate (malate) in the solution. In fact, the high concentrations of malate used in experiments appear high enough to abolish any channeling of oxaloacetate. The simulations provide a resolution of a serious discrepancy between previous simulations and experiments and raise important questions relating to the observability of electrostatically mediated substrate channeling in vitro and in vivo.     

Porphyrin and metalloporphyrin binding to DNA polymers: rate and equilibrium binding studies

Interactions of meso-tetrakis(4-N-methylpyridiniumyl)porphyrin [TMpyP(4)] with poly[d(G-C)].poly[d(G-C)] [poly[d(G-C)2] and poly[d(A-T)].poly[d(A-T)] [poly[d(A-T)2] were studied by equilibrium dialysis and stopped-flow dissociation kinetics as a function of [Na+]. Metalloderivatives of TMpyP(4), NiTMpyP(4), and ZnTMpyP(4) were also investigated. The apparent equilibrium binding constants (Kobs) were approximately the same for TMpyP(4) binding to either poly[d(G-C)2] or poly[d(A-T)2] and decreased with increasing [Na+]. The slopes of the plots of log Kobs vs log [Na+] were similar, with values close to -2.7. Contrary to implications in previously reported studies, these data do not indicate that TMpyP(4) prefers to bind to GC sites at low ionic strength and to AT sites at high ionic strength. In contrast, binding of ZnTMpyP(4) to these two polymers is very different. Comparisons of Kobs values at 0.065 M [Na+] indicate that ZnTMpyP(4) binding to AT sites is approximately 200 times more favorable than binding to GC sites, a finding in agreement with previous qualitative observations. Although the binding of the Zn species to the GC polymer was too weak for us to assess the salt effect, the plot of log Kobs vs log [Na+] gave a slope of -2.0 for ZnTMpyP(4) binding to poly[d(A-T)2]. Application of condensation theory for polyelectrolytes suggests similar charge interactions for ZnTMpyP(4) and for TMpyP(4) binding to poly[d(A-T)2]. Likewise, the rates of dissociation from poly[d(A-T)2] were similar for TMpyP(4) and ZnTMpyP(4) [and also NiTMpyP(4)]. However, whereas TMpyP(4) [and NiTMpyP(4)] dissociation from poly[d(G-C)2] was measurable, that for ZnTMpyP(4) was too fast to measure.(ABSTRACT TRUNCATED AT 250 WORDS)     

Coarse-grained models for simulations of multiprotein complexes: application to ubiquitin binding

We develop coarse-grained models and effective energy functions for simulating thermodynamic and structural properties of multiprotein complexes with relatively low binding affinity (K_d > 1 μM) and apply them to binding of Vps27 to membrane-tethered ubiquitin. Folded protein domains are represented as rigid bodies. The interactions between the domains are treated at the residue level with amino-acid-dependent pair potentials and Debye-Hückel-type electrostatic interactions. Flexible linker peptides connecting rigid protein domains are represented as amino acid beads on a polymer with appropriate stretching, bending, and torsion-angle potentials. In simulations of membrane-attached protein complexes, interactions between amino acids and the membrane are described by residue-dependent short-range potentials and long-range electrostatics. We parameterize the energy functions by fitting the osmotic second virial coefficient of lysozyme and the binding affinity of the ubiquitin-CUE complex. For validation, extensive replica-exchange Monte Carlo simulations are performed of various protein complexes. Binding affinities for these complexes are in good agreement with the experimental data. The simulated structures are clustered on the basis of distance matrices between two proteins and ranked according to cluster population. In ∼ 70% of the complexes, the distance root-mean-square is less than 5 Å from the experimental structures. In ∼ 90% of the complexes, the binding interfaces on both proteins are predicted correctly, and in all other cases at least one interface is correct. Transient and nonspecifically bound structures are also observed. With the validated model, we simulate the interaction between the Vps27 multiprotein complex and a membrane-tethered ubiquitin. Ubiquitin is found to bind preferentially to the two UIM domains of Vps27, but transient interactions between ubiquitin and the VHS and FYVE domains are observed as well. These specific and nonspecific interactions are found to be positively cooperative, resulting in a substantial enhancement of the overall binding affinity beyond the ∼ 300 μM of the specific domains. We also find that the interactions between ubiquitin and Vps27 are highly dynamic, with conformational rearrangements enabling binding of Vps27 to diverse targets as part of the multivesicular-body protein-sorting pathway.     

Molecular dynamics simulations and free energy calculations of netropsin and distamycin binding to an AAAAA DNA binding site

Molecular dynamics simulations have been performed on netropsin in two different charge states and on distamycin binding to the minor groove of the DNA duplex d(CGCGAAAAACGCG)·d(CGCGTTTTTCGCG). The relative free energy of binding of the two non-covalently interacting ligands was calculated using the thermodynamic integration method and reflects the experimental result. From 2 ns simulations of the ligands free in solution and when bound to DNA, the mobility and the hydrogen-bonding patterns of the ligands were studied, as well as their hydration. It is shown that even though distamycin is less hydrated than netropsin, the loss of ligand–solvent interactions is very similar for both ligands. The relative mobilities of the ligands in their bound and free forms indicate a larger entropic penalty for distamycin when binding to the minor groove compared with netropsin, partially explaining the lower binding affinity of the distamycin molecule. The detailed structural and energetic insights obtained from the molecular dynamics simulations allow for a better understanding of the factors determining ligand–DNA binding.     

Analyzing thioflavin T binding to amyloid fibrils by an equilibrium microdialysis-based technique

A new approach for the determination of the amyloid fibril - thioflavin T (ThT) binding parameters (the number of binding modes, stoichiometry, and binding constants of each mode) is proposed. This approach is based on the absorption spectroscopy determination of the concentration of free and bound to fibril dye in solutions, which are prepared by equilibrium microdialysis. Furthermore, the proposed approach allowed us, for the first time, to determine the absorption spectrum, molar extinction coefficient, and fluorescence quantum yield of the ThT bound to fibril by each binding modes. This approach is universal and can be used for determining the binding parameters of any dye interaction with a receptor, such as ANS binding to proteins in the molten globule state or to protein amorphous aggregates.     

Equilibrium binding of myristoyllysophosphatidylcholine to bovine myelin basic protein: an example of ligand-mediated acceptor association

The interaction of myristoyllysophosphatidylcholine with bovine myelin basic protein at pH 7.4 and 4.5, I = 0.48, has been investigated by a recycling partition equilibrium technique with Bio-Gel P-2 as the gel phase. Important points to emerge from this direct binding study are that it is a monomeric (not micellar) amphiphile that binds to myelin basic protein, that the amphiphile binds preferentially to the monomeric form of myelin basic protein, that this binding to monomer is highly cooperative, that the similarity of binding behavior in the two environments tested is consistent with the dominance of a hydrophobic contribution to the protein-amphiphile interaction, and that the self-association of myelin basic protein in the presence of phospholipid [Smith, R. (1982) Biochemistry 21, 2697-2701] must reflect the aggregation of a protein-amphiphile complex(es) coupled with concomitant release of some lipid. These findings are then related to earlier nuclear magnetic resonance and circular dichroism studies in which the results were interpreted on the basis that myelin basic protein bound preferentially to micellar phospholipid.     

Equilibrium constant for binding of an actin filament capping protein to the barbed end of actin filaments

Depolymerization of treadmilling actin filaments by a capping protein isolated from bovine brain was used for determination of the equilibrium constant for binding of the capping protein to the barbed ends of actin filaments. When the capping protein blocks monomer consumption at the lengthening barbed ends, monomers continue to be produced at the shortening pointed ends until a new steady state is reached in which monomer production at the pointed ends is balanced by monomer consumption at the uncapped barbed ends. In this way the ratio of capped to uncapped filaments could be determined as a function of the capping protein concentration. Under the experimental conditions (100 mM KCl and 2 mM MgCl2, pH 7.5, 37 degrees C) the binding constant was found to be about 2 X 10(9) M-1. Capping proteins effect the actin monomer concentration only at capping protein concentrations far above the reciprocal of their binding constant. Half-maximal increase of the monomer concentration requires capping of about 99% of the actin filaments. A low proportion of uncapped filaments has a great weight in determining the monomer concentration because association and dissociation reactions occur at the dynamic barbed ends with higher frequencies than at the pointed ends.     

AlexSys: a knowledge-based expert system for multiple sequence alignment construction and analysis

Multiple sequence alignment (MSA) is a cornerstone of modern molecular biology and represents a unique means of investigating the patterns of conservation and diversity in complex biological systems. Many different algorithms have been developed to construct MSAs, but previous studies have shown that no single aligner consistently outperforms the rest. This has led to the development of a number of ‘meta-methods’ that systematically run several aligners and merge the output into one single solution. Although these methods generally produce more accurate alignments, they are inefficient because all the aligners need to be run first and the choice of the best solution is made a posteriori. Here, we describe the development of a new expert system, AlexSys, for the multiple alignment of protein sequences. AlexSys incorporates an intelligent inference engine to automatically select an appropriate aligner a priori, depending only on the nature of the input sequences. The inference engine was trained on a large set of reference multiple alignments, using a novel machine learning approach. Applying AlexSys to a test set of 178 alignments, we show that the expert system represents a good compromise between alignment quality and running time, making it suitable for high throughput projects. AlexSys is freely available from http://alnitak.u-strasbg.fr/∼aniba/alexsys.     

A quantitative analysis of apolipoprotein binding to SR-BI: multiple binding sites for lipid-free and lipid-associated apolipoproteins

Competitive binding experiments were performed using Y1-BS1 adrenal cells to provide information about the interaction of HDL apolipoproteins with scavenger receptor class B, type I (SR-BI). Exchangeable apolipoproteins apolipoprotein A-I (apoA-I), apoA-II, apoE-2, apoE-3, and apoE-4 as phospholipid complexes bind like HDL3 to SR-BI via their multiple amphipathic alpha-helices; the concentrations required to reduce the binding of HDL3 to SR-BI by 50% (IC50) were similar and in the range of 35-50 microgram protein/ml. In the case of apoA-I, peptides corresponding to segments 1-85, 44-65, 44-87, 149-243, and 209-241 all had the same IC50 as each other (P = 0.86), showing that a specific amino acid sequence in apoA-I is not responsible for the interaction with SR-BI. The distribution of charged residues in the amphipathic alpha-helix affects the interaction, with class A and Y helices binding better than class G* helices. Synthetic alpha-helical peptides composed of either l or d amino acids can bind equally to the receptor. Association with phospholipid increases the amount of apolipoprotein binding to SR-BI without altering the affinity of binding. Lipid-free apolipoproteins compete only partially with the binding of HDL to SR-BI, whereas lipidated apolipoproteins compete fully. These results are consistent with the existence of more than one type of apolipoprotein binding site on SR-BI.     

HU binding to DNA: evidence for multiple complex formation and DNA bending

HU, a nonspecific histone-like DNA binding protein, participates in a number of genomic events as an accessory protein and forms multiple complexes with DNA. The HU-DNA binding interaction was characterized by fluorescence, generated with the guanosine analogue 3-methyl-8-(2-deoxy-beta-D-ribofuranosyl)isoxanthopterin (3-MI) directly incorporated into DNA duplexes. The stoichiometry and equilibrium binding constants of complexes formed between HU and 13 and 34 bp DNA duplexes were determined using fluorescence anisotropy and analytical ultracentrifugation. These measurements reveal that three HU molecules bind to the 34 bp duplexes, while two HU molecules bind to the 13 bp duplex. The data are well described by an independent binding site model, and the association constants for the first binding event for both duplexes are similar (approximately 1 x 10(6) M(-1)), indicating that HU binding affinity is independent of duplex length. Further analysis of the binding curves in terms of a nonspecific binding model is indicative that HU binding to DNA exhibits little to no cooperativity. The fluorescence intensity also increases upon HU binding, consistent with decreased base stacking and increased solvent exposure of the 3-MI fluorescence probe. These results are suggestive of a local bending or unwinding of the DNA. On the basis of these results we propose a model in which bending of DNA accompanies HU binding. Up to five complex bands are observed in gel mobility shift assays of HU binding to the 34 bp duplexes. We suggest that protein-induced bending of the DNA leads to the observation of complexes in the gel, which have the same molecular weight but different relative mobilities.     

Free energy simulations for protein ligand binding and stability

Abstract

We summarize several computational techniques to determine relative free energies for condensed-phase systems. The focus is on practical considerations which are capable of making direct contact with experiments. Particular applications include the thermodynamic stability of apo- and holo-myoglobin, insulin dimerization free energy, ligand binding in lysozyme, and ligand diffusion in globular proteins. In addition to provide differential free energies between neighboring states, converged umbrella sampling simulations provide insight into migration barriers and ligand dissociation barriers and analysis of the trajectories yield additional insight into the structural dynamics of fundamental processes. Also, such simulations are useful tools to quantify relative stability changes for situations where experiments are difficult. This is illustrated for NO-bound myoglobin. For the dissociation of benzonitrile from lysozyme it is found that long umbrella sampling simulations are required to approximately converge the free energy profile. Then, however, the resulting differential free energy between the bound and unbound state is in good agreement with estimates from molecular mechanics with generalized Born surface area simulations. Furthermore, comparing the barrier height for ligand escape suggests that ligand dissociation contains a non-equilibrium component.     

Interdomain dynamics and ligand binding: molecular dynamics simulations of glutamine binding protein

Periplasmic binding proteins from Gram-negative bacteria possess a common architecture, comprised of two domains linked by a hinge region, a fold which they share with the neurotransmitter-binding domains of ionotropic glutamate receptors (GluRs). Glutamine-binding protein (GlnBP) is one such protein, whose crystal structure has been solved in both open and closed forms. Multi-nanosecond molecular dynamics simulations have been used to explore motions about the hinge region and how they are altered by ligand binding. Glutamine binding is seen to significantly reduce inter-domain motions about the hinge region. Essential dynamics analysis of inter-domain motion revealed the presence of both hinge-bending and twisting motions, as has been reported for a related sugar-binding protein. Significantly, the influence of the ligand on GlnBP dynamics is similar to that previously observed in simulations of rat glutamate receptor (GluR2) ligand-binding domain. The essential dynamics analysis of GlnBP also revealed a third class of motion which suggests a mechanism for signal transmission in GluRs.     

Computer simulations and experimental studies of gel mobility patterns for weak and strong non-cooperative protein binding to two targets on the same DNA: application to binding of tet repressor variants to multiple and single tet operator sites

A series of computer simulations of gel patterns assuming non-cooperative binding of a protein to two targets on the same DNA fragment was performed and applied to interprete gel mobility shift experiments of Tet repressor-tet operator binding. While a high binding affinity leads to the expected distribution of free DNA, DNA bound by one repressor dimer and DNA bound by two repressor dimers, a lower affinity or an increased electrophoresis time results in the loss of the band corresponding to the singly occupied complex. The doubly occupied complex remains stable under these conditions. This phenomenon is typical for protein binding to DNA fragments with two identical sites. It results from statistical disproportionation of the singly occupied complex in the gel. The lack of the singly occupied complex is commonly taken to indicate cooperative binding, however, our analysis shows clearly, that cooperativity is not needed to interprete these results. Tet repressor proteins and small DNA fragments with two tet operator sites have been prepared from four classes of tetracycline resistance determinants. The results of gel mobility shift analyses of various complexes of these compounds confirm the predictions. Furthermore, calculated gel patterns assuming different gel mobilities of the two singly occupied complexes show discrete bands only if the electrophoresis time is shorter than the inverse of the microscopic dissociation rate constant. Simulations assuming increasing dissociation rates predict that the two bands first merge into one, which then disappears. This behavior was verified by gel mobility analyses of Tet repressor-tet operator titrations at increased salt concentrations as well as by direct footprinting of the complexes in the gel. It is concluded that comparison of the intensities of the single and the double occupation bands allow a rough estimation of the dissociation rate constant. On this basis the sixteen possible Tet repressor-tet operator combinations can be ordered with decreasing binding affinities by a simple gel shift experiment. The implications of these results for gel mobility analyses of other protein-DNA complexes are discussed.     

Rate constants and equilibrium constants for binding of actin to the 1:1 gelsolin-actin complex.

The rate constant and equilibrium constant of association of an actin monomer with 1:1 gelsolin-actin complex isolated from chicken were measured by using fluorescently labeled actin. According to fluorescence stopped-flow experiments, the rate constant of formation of the 1:2 gelsolin-actin complex from 1:1 gelsolin-actin complex and actin was found to be about 2 x 10(7) M-1 s-1 under conditions where gelsolin binds Ca2+. The rate of dissociation of one actin molecule from the 1:2 gelsolin-actin complex was determined by exchange of actin for fluorescently labeled actin. The rate constant of dissociation was about 0.02 s-1. Thus, the equilibrium constant for association of actin with 1:1 gelsolin-actin complex can be calculated to be in the range of 10(9) M-1. The rate of dissociation of actin from 1:2 gelsolin-actin complex was independent of the Ca2+ concentration. Ca2+ affects only the rate of association of actin with 1:1 gelsolin-actin complex.     

A modified equilibrium dialysis method for studying fatty acid binding to proteins

A modified equilibrium dialysis method is described which is suitable for investigating the binding of fatty acids in the form of aqueous micellar dispersions to proteins. The method uses a permeant chromophore which complexes reversibly with free fatty acid within the dialysis bag. The concentration outside the dialysis bag is determined spectrophotometrically. Binding of oleic acid to bovine serum albumin is given as an example. A simplified analysis of fatty acid binding is given and used to indicate the potential of the method.     

Human alpha-thrombin binding to nonpolymerized fibrin-Sepharose: evidence for an anionic binding region

In order to investigate ligand binding sites in alpha-thrombin that interact with nonpolymerized fibrin, fibrinogen was conjugated (with CNBr) to Sepharose 4B and converted to the nonpolymerized fibrin resin with alpha-thrombin. Human alpha-thrombin was bound to the resin at 22 degrees C and eluted with a linear NaCl gradient [50-300 mM in 50 mM tris(hydroxymethyl)aminomethane hydrochloride, pH 7.6] with midpeak elution occurring at an ionic strength that corresponds to 170 +/- 5 mM NaCl. Among various ligands examined, ATP and its analogues caused alpha-thrombin to elute with 125 mM or less salt. Apparent dissociation constants were estimated by the dependence of elution volume on ligand concentration. The most potent ligands for desorption from the column were anionic (e.g., adenine nucleotides), which also inhibit thrombin esterolytic/amidolytic and clotting activity [Conery, B. G., & Berliner, L. J. (1983) Biochemistry 22, 369-375]. The desorption series was at 10 mM concentrations: ATP = ADP greater than pyrophosphate greater than citrate greater than oxalate greater than PO4(3-). Contrastingly, serotonin and related apolar compounds did not cause dissociation of alpha-thrombin from the fibrin resin, even though several of these substances inhibit fibrinogen clotting and esterolytic/amidolytic activities of the enzyme. These data imply that independent sites for apolar and anionic binding in alpha-thrombin are required for converting fibrinogen into clottable fibrin and that alpha-thrombin-fibrin binding involves an anionic site.