Ligand entry into the calyx of β‐lactoglobulin

Although the thermodynamic principles that control the binding of drug molecules to their protein targets are well understood, the detailed process of how a ligand reaches a protein binding site has been an intriguing question over decades. The short time interval between the encounter between a ligand and its receptor to the formation of the stable complex has prevented experimental observations. Bovine β‐lactoglobulin (βlg) is a lipocalin member that carries fatty acids (FAs) and other lipids in the cellular environment. Βlg accommodates a FA molecule in its highly hydrophobic cavity and exhibits the capability of recognizing a wide variety of hydrophobic ligands. To elucidate the ligand entry process on βlg, we report molecular dynamics simulations of the encounter between palmitate (PA) or laurate (LA) and βlg. Our results show that residues localized in loops at the cavity entrance play an important role in the ligand penetration process. Analysis of the short‐term interaction energies show that the forces operating on the systems lead to average conformations very close to the crystallographic holo‐forms. Whereas the binding free energy analysis using the molecular mechanics Generalized Born surface area method shows that these conformations were thermodynamically favorable. © 2013 Wiley Periodicals, Inc. Biopolymers 101: 744–757, 2014.     

Formation and characterization of a transition state complex of Azotobacter vinelandii nitrogenase

A stable complex is formed between the nitrogenase proteins of Azotobacter vinelandii, aluminium fluoride and MgADP. All nitrogenase activities are inhibited. The complex formation was found to be reversible. An incubation at 50°C recovers nitrogenase activity. The complex has been characterized with respect to protein and nucleotide composition and redox state of the metal-sulphur clusters. Based on the inhibition by aluminium fluoride together with MgADP, it is proposed that a stable transition state complex of nitrogenase is isolated.     

The courtship of proteins: Understanding the encounter complex

The formation of protein complexes involves an encounter complex, in which proteins show few specific interactions and assume many orientations. Recent kinetic and structural studies have shed light on this elusive state. It is generally dominated by electrostatic interactions, although hydrophobic interactions can play a role. During the encounter phase the proteins remain largely solvated. In extreme cases, the proteins only form an encounter complex, and in many other complexes, the encounter state constitutes a significant amount (5% or more), indicating that the energy difference between encounter and productive complexes is small. Thus, the encounter complex represents an essential part of protein complexes.     

A mechanistic study of ferrioxamine B reduction by the biological reducing agent ascorbate in the presence of an iron(II) chelator

The iron overload drug desferal (desferrioxamine B) forms the stable iron complex ferrioxamine B. The reduction potential of ferrioxamine B (E^o = −482 mV versus NHE pH 7) prohibits its reduction by biological reducing agents such as ascorbate, but it was found that the iron(II) chelator 2,2′-bipyridine (bipy) facilitates this reduction. Evidence is given to support the formation of a ternary complex between iron, bipy, and desferrioxamine B as the key step in facilitating the reduction. The equilibrium constant for the formation of the ternary complex was found to be 8.9 × 10⁷ and ternary complex formation is explained in terms of a three step mechanism. The mechanism for the reduction of ferrioxamine B is discussed in terms of rapidly established pre-equilibria which include ternary complex formation, ascorbic acid deprotonation, and encounter complex formation between ascorbate and the ternary complex. These equilibria are followed by rate limiting reduction of the ternary complex. Bipy was found to be a similar facilitator to sulfonated bathophenanthroline for the reduction of ferrioxamine B by ascorbate.     

Conformational transitions in protein-protein association: binding of fasciculin-2 to acetylcholinesterase

The neurotoxin fasciculin-2 (FAS2) is a picomolar inhibitor of synaptic acetylcholinesterase (AChE). The dynamics of binding between FAS2 and AChE is influenced by conformational fluctuations both before and after protein encounter. Submicrosecond molecular dynamics trajectories of apo forms of fasciculin, corresponding to different conformational substates, are reported here with reference to the conformational changes of loop I of this three-fingered toxin. This highly flexible loop exhibits an ensemble of conformations within each substate corresponding to its functions. The high energy barrier found between the two major substates leads to transitions that are slow on the timescale of the diffusional encounter of noninteracting FAS2 and AChE. The more stable of the two apo substates may not be the one observed in the complex with AChE. It seems likely that the more stable apo form binds rapidly to AChE and conformational readjustments then occur in the resulting encounter complex.     

Dynamics in transient complexes of redox proteins

Recent studies have provided experimental information about the initial stage of protein complex formation, the encounter complex. This stage is particularly important in the weak and transient complexes formed between electron transfer proteins and their partners. These studies are discussed and the role of the encounter complex is interpreted in terms of the specific requirements that the biological function puts on these complexes.     

An allosteric-feedback mechanism for protein-assisted group I intron splicing

The I-AniI maturase facilitates self-splicing of a mitochondrial group I intron in Aspergillus nidulans. Binding occurs in at least two steps: first, a specific but labile encounter complex rapidly forms and then this intermediate is slowly resolved into a native, catalytically active RNA/protein complex. Here we probe the structure of the RNA throughout the assembly pathway. Although inherently unstable, the intron core, when bound by I-AniI, undergoes rapid folding to a near-native state in the encounter complex. The next transition includes the slow destabilization and docking into the core of the peripheral stacked helix that contains the 5' splice site. Mutational analyses confirm that both transitions are important for native complex formation. We propose that protein-driven destabilization and docking of the peripheral stacked helix lead to subtle changes in the I-AniI binding site that facilitate native complex formation. These results support an allosteric-feedback mechanism of RNA-protein recognition in which proteins engaged in an intermediate complex can influence RNA structure far from their binding sites. The linkage of these changes to stable binding ensures that the protein and RNA do not get sequestered in nonfunctional complexes.     

On the validity of the quasi-steady state approximation of bimolecular reactions in solution

Two-step binding kinetics are extensively used to study the relative importance of diffusion in biochemical reactions. Classical analysis of this problem assumes ad hoc that the encounter complex is at quasi-steady state (QSS). Using scaling arguments we derive a criterion for the validity of this assumption in the limit of irreversible product formation. We find that the QSS approximation (QSSA) of two-step binding is only valid if the total ligand and receptor concentrations are much smaller than (k2+k-1)/k1, where k1 and k-1 are, respectively, the forward and reverse diffusion encounter rate constants and k2 is the chemical association rate constant. This criterion can be shown to imply that the average time between encounters is much longer than the half-life of the encounter complex and also guarantees that the concentration of the encounter complex is negligible compared to the reactant and product concentrations. Numerical examples of irreversible and reversible cases corroborate our analysis and illustrate that the QSS may be invalid even if k-2相似文献   

New insights into the mechanism of protein-protein association

Association of a protein complex follows a two-step mechanism, with the first step being the formation of an encounter complex that evolves into the final complex. Here, we analyze recent experimental data of the association of TEM1-beta-lactamase with BLIP using theoretical calculations and simulation. We show that the calculated Debye-Hückel energy of interaction for a pair of proteins during association resembles an energy funnel, with the final complex at the minima. All attraction is lost at inter-protein distances of 20 A, or rotation angles of >60 degrees from the orientation of the final complex. For faster-associating protein complexes, the energy funnel deepens and its volume increases. Mutations with the largest impact on association (hotspots for association) have the largest effect on the size and depth of the energy funnel. Analyzing existing evidence, we suggest that the transition state along the association pathway is the formation of the final complex from the encounter complex. Consequently, pairs of proteins forming an encounter complex will tend to dissociate more readily than to evolve into the final complex. Increasing directional diffusion by increasing favorable electrostatic attraction results in a faster forming and slower dissociating encounter complex. The possible applicability of electrostatic calculations for protein-protein docking is discussed.     

Visualizing induced fit in early assembly of the human signal recognition particle

Assembly of almost all ribonucleoprotein complexes involves induced fit in the RNA and, thus, formation of one or more intermediate states. In assembly of the human signal recognition particle (SRP), we show that SRP19 binding to SRP RNA involves obligatory intermediates. An apparent discrepancy exists between the ratio of dissociation and association rate constants, determined in a partitioning experiment, and the equilibrium binding constant; this kinetic signature reflects formation of a stable intermediate in assembly of the ribonucleoprotein complex. Assembly intermediates were observed directly by time-resolved footprinting. SRP19 binds rapidly to SRP RNA to form an initial labile, but structurally specific, encounter complex involving both helices III and IV. Two subsequent steps of structural consolidation yield the native RNA-protein interface. SRP19 binding stabilizes helix IV in the region recognized by SRP54, consistent with protein-protein cooperativity mediated in part by mutual recognition of similar RNA structures. This mechanism illustrates principles general to ribonucleoprotein assembly reactions that rely on recruitment of architectural RNA binding proteins.     

Temperature differentially affects encounter and docking thermodynamics of antibody--antigen association

Using BIACORE SPR, we have examined the mechanism of temperature effects on the binding kinetics of two closely related antibody Fabs (H10 and H26) which recognize coincident epitopes on hen egg-white lysozyme (HEL), and whose association and dissociation kinetics are best described by the two-step conformational change model which we interpret as molecular encounter and docking. Time-course series data obtained at a series of six temperatures (6, 10, 15, 25, 30 and 37 degrees C) showed that temperature differentially affects the rate constants of the encounter and docking steps. Docking is more temperature-sensitive than the encounter step, and energetically less favorable at higher temperatures. At elevated temperatures, the time required for docking is longer and the apparent increase in off-rate reflects the greater proportion of the molecules failing to dock and remaining in the less stable encounter state. As a consequence, distribution of free energy change between the encounter and docking steps is altered. At physiological temperature (37 degrees C) the docking step of the H26 complex is energetically unfavorable and most complexes essentially do not dock. There is a significant decrease in total free energy change of the H26 complex at higher temperatures. Elevated temperature changes the rate-limiting step of H26--HEL association from the encounter to the docking step, but not that of H10--HEL. Our results indicate that the mechanism by which elevated temperature reduces the affinities of antigen--antibody complexes is to decrease the net docking rate, and/or stability of the docked complex; at higher temperatures, a smaller proportion of the complexes actually anneal to a more stable docked state. This mechanism may have broad applicability to other receptor--ligand complexes.     

Rapid-reaction kinetic characterization of the pathway of streptokinase-plasmin catalytic complex formation

Binding of the fibrinolytic proteinase plasmin (Pm) to streptokinase (SK) in a tight stoichiometric complex transforms Pm into a potent proteolytic activator of plasminogen. SK binding to the catalytic domain of Pm, with a dissociation constant of 12 pm, is assisted by SK Lys(414) binding to a Pm kringle, which accounts for a 11-20-fold affinity decrease when Pm lysine binding sites are blocked by 6-aminohexanoic acid (6-AHA) or benzamidine. The pathway of SK.Pm catalytic complex formation was characterized by stopped-flow kinetics of SK and the Lys(414) deletion mutant (SKDeltaK414) binding to Pm labeled at the active site with 5-fluorescein ([5F]FFR-Pm) and the reverse reactions by competitive displacement of [5F]FFR-Pm with active site-blocked Pm. The rate constants for the biexponential fluorescence quenching caused by SK and SKDeltaK414 binding to [5F]FFR-Pm were saturable as a function of SK concentration, reporting encounter complex affinities of 62-110 nm in the absence of lysine analogs and 4900-6500 and 1430-2200 nm in the presence of 6-AHA and benzamidine, respectively. The encounter complex with SKDeltaK414 was approximately 10-fold weaker in the absence of lysine analogs but indistinguishable from that of native SK in the presence of 6-AHA and benzamidine. The studies delineate for the first time the sequence of molecular events in the formation of the SK.Pm catalytic complex and its regulation by kringle ligands. Analysis of the forward and reverse reactions supports a binding mechanism in which SK Lys(414) binding to a Pm kringle accompanies near-diffusion-limited encounter complex formation followed by two slower, tightening conformational changes.     

Potential Regulatory Role of Competitive Encounter Complexes in Paralogous Phosphotransferase Systems

There are two paralogous Escherichia coli phosphotransferase systems, one for sugar import (PTS^sugar) and one for nitrogen regulation (PTS^Ntr), that utilize proteins enzyme I^sugar (EI^sugar) and HPr, and enzyme I^Ntr (EI^Ntr) and NPr, respectively. The enzyme I proteins have similar folds, as do their substrates HPr and NPr, yet they show strict specificity for their cognate partner both in stereospecific protein–protein complex formation and in reversible phosphotransfer. Here, we investigate the mechanism of specific EI^Ntr:NPr complex formation by the study of transient encounter complexes. NMR paramagnetic relaxation enhancement experiments demonstrated transient encounter complexes of EI^Ntr not only with the expected partner, NPr, but also with the unexpected partner, HPr. HPr occupies transient sites on EI^Ntr but is unable to complete stereospecific complex formation. By occupying the non-productive transient sites, HPr promotes NPr transient interaction to productive sites closer to the stereospecific binding site and actually enhances specific complex formation between NPr and EI^Ntr. The cellular level of HPr is approximately 150 times higher than that of NPr. Thus, our finding suggests a potential mechanism for cross-regulation of enzyme activity through formation of competitive encounter complexes.     

Forced Dissociation of the Strand Dimer Interface between C-Cadherin Ectodomains

The force-induced dissociation of the strand dimer interface in C-cadherin has been studied using steered molecular dynamics simulations. The dissociation occurred, without domain unraveling, after the extraction of the conserved trypthophans (Trp2) from their respective hydrophobic pockets. The simulations revealed two stable positions for the Trp2 side chain inside the pocket. The most internal stable position involved a hydrogen bond between the ring Ne of Trp2 and the backbone carbonyl of Glu90. In the second stable position, the aromatic ring is located at the pocket entrance. After extracting the two tryptophans from their pockets, the complex exists in an intermediate bound state that involves a close packing of the tryptophans with residues Asp1 and Asp27 from both domains. Dissociation occurred after this residue association was broken. Simulations carried out with a complex formed between W2A mutants showed that the mutant complex dissociates more easily than the wild type complex does. These results correlate closely with the role of the conserved tryptophans suggested previously by site directed mutagenesis.     

Calcium enhances heparin catalysis of the antithrombin-factor Xa reaction by promoting the assembly of an intermediate heparin-antithrombin-factor Xa bridging complex. Demonstration by rapid kinetics studies

Heparin catalyzes the inhibition of factor Xa by antithrombin mainly through an allosteric activation of the serpin inhibitor, but an alternative heparin bridging mechanism has been suggested to enhance the catalysis in the presence of physiologic calcium levels due to calcium interactions with the Gla domain exposing a heparin binding exosite in factor Xa. To provide direct evidence for this bridging mechanism, we studied the heparin-catalyzed reaction of antithrombin with factor Xa, Gla-domainless factor Xa (GDFXa), and a heparin binding exosite mutant of GDFXa in the absence and presence of calcium using rapid kinetic methods. The pseudo-first-order rate constant for factor Xa inhibition by antithrombin complexed with a long-chain approximately 70-saccharide heparin showed a saturable dependence on inhibitor concentration in the presence but not in the absence of 2.5 mM Ca(2+), indicating the formation of an intermediate heparin-serpin-proteinase encounter complex with a dissociation constant of approximately 90 nM prior to formation of the stable serpin-proteinase complex with a rate constant of approximately 20 s(-1). Similar saturation kinetics were observed for the inhibition of GDFXa by the antithrombin-heparin complex, except that Ca(2+) was not required for the effect. By contrast, no Ca(2+)-dependent saturation of the inhibition rate constant was detectable over the same range of inhibitor concentrations for reactions of either a short-chain approximately 26-saccharide high-affinity heparin-antithrombin complex with factor Xa or the long-chain heparin-antithrombin complex with the heparin binding exosite mutant, GDFXa R240A. These findings suggest that binding of full-length heparin chains to an exosite of factor Xa in the presence of Ca(2+) produces a chain-length-dependent lowering of the dissociation constant for assembly of the intermediate heparin-antithrombin-factor Xa encounter complex, resulting in a several 100-fold rate enhancement by a heparin bridging mechanism.     

Binding of azurin to cytochrome c 551 as investigated by surface plasmon resonance and fluorescence

The interaction between azurin (Az) and cytochrome c 551 (CytC551) from Pseudomonas aeruginosa deserves particular interest for both its physiological aspects and their possible applications in bionano devices. Here, the kinetics of the interaction has been studied by surface plasmon resonance and fluorescence quenching. Surface plasmon resonance data have been successfully interpreted by the heterogeneous ligand model, which predicts the existence of two binding sites on the immobilized Az for CytC551 molecules in solution. On the other hand, the fluorescence study indicates the formation of a complex, with the involvement of the lone Az tryptophan (Trp) at position 48. The two different techniques point out the occurrence of an encounter complex between Az and CytC551 that evolves toward the formation of a more stable complex characterized by an equilibrium dissociation constant K_D typical of transient interactions. Copyright © 2014 John Wiley & Sons, Ltd.     

On the origin of the stereoselective affinity of Nutlin-3 geometrical isomers for the MDM2 protein

The stereoselective affinity of small-molecule binding to proteins is typically broadly explained in terms of the thermodynamics of the final bound complex. Using Brownian dynamics simulations, we show that the preferential binding of the MDM2 protein to the geometrical isomers of Nutlin-3, an effective anticancer lead that works by inhibiting the interaction between the proteins p53 and MDM2, can be explained by kinetic arguments related to the formation of the MDM2:Nutlin-3 encounter complex. This is a diffusively bound state that forms prior to the final bound complex. We find that the MDM2 protein stereoselectivity for the Nutlin-3a enantiomer stems largely from the destabilization of the encounter complex of its mirror image enantiomer Nutlin-3b, by the K70 residue that is located away from the binding site. On the other hand, the trans-Nutlin-3a diastereoisomer exhibits a shorter residence time in the vicinity of MDM2 compared with Nutlin-3a due to destabilization of its encounter complex by the collective interaction of pairs of charged residues on either side of the binding site: Glu25 and Lys51 on one side, and Lys94 and Arg97 on the other side. This destabilization is largely due to the electrostatic potential of the trans-Nutlin-3a isomer being largely positive over extended continuous regions around its structure, which are otherwise well-identified into positive and negative regions in the case of the Nutlin-3a isomer. Such rich insight into the binding processes underlying biological selectivity complements the static view derived from the traditional thermodynamic analysis of the final bound complex. This approach, based on an explicit consideration of the dynamics of molecular association, suggests new avenues for kinetics-based anticancer drug development and discovery.     

Fruitful and Futile Encounters along the Association Reaction between Proteins

The association reaction between pairs of proteins proceeds through an encounter complex that develops into the final complex. Here, we combined Brownian dynamics simulations with experimental studies to analyze the structures of the encounter complexes along the association reaction between TEM1-β-lactamase and its inhibitor, β-lactamase-inhibitor protein. The encounter complex can be considered as an ensemble of short-lived low free-energy states that are stabilized primarily by electrostatic forces and desolvation. For the wild-type, the simulation showed two main encounter regions located outside the physical binding site. One of these regions was located near the experimentally determined transition state. To validate whether these encounters are fruitful or futile, we examined three groups of mutations that altered the encounter. The first group consisted of mutations that increased the experimental rate of association through electrostatic optimization. This resulted in an increase in the size of the encounter region located near the experimentally determined transition state, as well as a decrease in the energy of this region and an increase in the number of successful trajectories (i.e., encounters that develop into complex). A second group of mutations was specifically designed to either increase or decrease the size and energy of the second encounter complex, but either way it did not affect k_on. A third group of mutations consisted of residues that increased k_on without significantly affecting the encounter complexes. These results indicate that the size and energy of the encounter regions are only two of several parameters that lead to fruitful association, and that electrostatic optimization is a major driving force in fast association.     

On the pathway of forming enzymatically productive ligand-protein complexes in lactate dehydrogenase

We have carried out a series of studies on the binding of a substrate mimic to the enzyme lactate dehydrogenase (LDH) using advanced kinetic approaches, which begin to provide a molecular picture of the dynamics of ligand binding for this protein. Binding proceeds via a binding-competent subpopulation of the nonligated form of the protein (the LDH/NADH binary complex) to form a protein-ligand encounter complex. The work here describes the collapse of the encounter complex to form the catalytically competent Michaelis complex. Isotope-edited static Fourier transform infrared studies on the bound oxamate protein complex reveal two kinds of oxamate environments: 1), a major populated structure wherein all significant hydrogen-bonding patterns are formed at the active site between protein and bound ligand necessary for the catalytically productive Michaelis complex and 2), a minor structure in a configuration of the active site that is unfavorable to carry out catalyzed chemistry. This latter structure likely simulates a dead-end complex in the reaction mixture. Temperature jump isotope-edited transient infrared studies on the binding of oxamate with LDH/NADH suggest that the evolution of the encounter complex between LDH/NADH and oxamate collapses via a branched reaction pathway to form the major and minor bound species. The production of the catalytically competent protein-substrate complex has strong similarities to kinetic pathways found in two-state protein folding processes. Once the encounter complex is formed between LDH/NADH and substrate, the ternary protein-ligand complex appears to “fold” to form a compact productive complex in an all or nothing like fashion with all the important molecular interactions coming together at the same time.