Conformational selection and induced changes along the catalytic cycle of Escherichia coli dihydrofolate reductase

Protein function often involves changes between different conformations. Central questions are how these conformational changes are coupled to the binding or catalytic processes during which they occur, and how they affect the catalytic rates of enzymes. An important model system is the enzyme dihydrofolate reductase (DHFR) from Escherichia coli, which exhibits characteristic conformational changes of the active‐site loop during the catalytic step and during unbinding of the product. In this article, we present a general kinetic framework that can be used (1) to identify the ordering of events in the coupling of conformational changes, binding, and catalysis and (2) to determine the rates of the substeps of coupled processes from a combined analysis of nuclear magnetic resonance R₂ relaxation dispersion experiments and traditional enzyme kinetics measurements. We apply this framework to E. coli DHFR and find that the conformational change during product unbinding follows a conformational‐selection mechanism, that is, the conformational change occurs predominantly prior to unbinding. The conformational change during the catalytic step, in contrast, is an induced change, that is, the change occurs after the chemical reaction. We propose that the reason for these conformational changes, which are absent in human and other vertebrate DHFRs, is robustness of the catalytic rate against large pH variations and changes to substrate/product concentrations in E. coli. Proteins 2012;. © 2012 Wiley Periodicals, Inc.     

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.     

On the molecular basis of the high affinity binding of basic amino acids to LAOBP,a periplasmic binding protein from Salmonella typhimurium

The rational designing of binding abilities in proteins requires an understanding of the relationship between structure and thermodynamics. However, our knowledge of the molecular origin of high‐affinity binding of ligands to proteins is still limited; such is the case for l ‐lysine–l ‐arginine–l ‐ornithine periplasmic binding protein (LAOBP), a periplasmic binding protein from Salmonella typhimurium that binds to l ‐arginine, l ‐lysine, and l ‐ornithine with nanomolar affinity and to l ‐histidine with micromolar affinity. Structural studies indicate that ligand binding induces a large conformational change in LAOBP. In this work, we studied the thermodynamics of l ‐histidine and l ‐arginine binding to LAOBP by isothermal titration calorimetry. For both ligands, the affinity is enthalpically driven, with a binding ΔCp of ~?300 cal mol^?1 K^?1, most of which arises from the burial of protein nonpolar surfaces that accompanies the conformational change. Osmotic stress measurements revealed that several water molecules become sequestered upon complex formation. In addition, LAOBP prefers positively charged ligands in their side chain. An energetic analysis shows that the protein acquires a thermodynamically equivalent state with both ligands. The 1000‐fold higher affinity of LAOBP for l ‐arginine as compared with l ‐histidine is mainly of enthalpic origin and can be ascribed to the formation of an extra pair of hydrogen bonds. Periplasmic binding proteins have evolved diverse energetic strategies for ligand recognition. STM4351, another arginine binding protein from Salmonella, shows an entropy‐driven micromolar affinity toward l ‐arginine. In contrast, our data show that LAOBP achieves nanomolar affinity for the same ligand through enthalpy optimization. Copyright © 2015 John Wiley & Sons, Ltd.     

Opening mechanism of a cyclic nucleotide-gated channel based on analysis of single channels locked in each liganded state.

Cyclic nucleotide-gated channels contain four subunits, each with a binding site for cGMP or cAMP in the cytoplasmic COOH-terminal domain. Previous studies of the kinetic mechanism of activation have been hampered by the complication that ligands are continuously binding and unbinding at each of these sites. Thus, even at the single channel level, it has been difficult to distinguish changes in behavior that arise from a channel with a fixed number of ligands bound from those that occur upon the binding and unbinding of ligands. For example, it is often assumed that complex behaviors like multiple conductance levels and bursting occur only as a consequence of changes in the number of bound ligands. We have overcome these ambiguities by covalently tethering one ligand at a time to single rod cyclic nucleotide-gated channels (Ruiz, ML., and J.W. Karpen. 1997. Nature. 389:389-392). We find that with a fixed number of ligands locked in place the channel freely moves between three conductance states and undergoes bursting behavior. Furthermore, a thorough kinetic analysis of channels locked in doubly, triply, and fully liganded states reveals more than one kinetically distinguishable state at each conductance level. Thus, even when the channel contains a fixed number of bound ligands, it can assume at least nine distinct states. Such complex behavior is inconsistent with simple concerted or sequential allosteric models. The data at each level of liganding can be successfully described by the same connected state model (with different rate constants), suggesting that the channel undergoes the same set of conformational changes regardless of the number of bound ligands. A general allosteric model, which postulates one conformational change per subunit in both the absence and presence of ligand, comes close to providing enough kinetically distinct states. We propose an extension of this model, in which more than one conformational change per subunit can occur during the process of channel activation.     

Structures of reduced and ligand‐bound nitric oxide reductase provide insights into functional differences in respiratory enzymes

Nitric oxide reductase (NOR) catalyzes the generation of nitrous oxide (N₂O) via the reductive coupling of two nitric oxide (NO) molecules at a heme/non‐heme Fe center. We report herein on the structures of the reduced and ligand‐bound forms of cytochrome c‐dependent NOR (cNOR) from Pseudomonas aeruginosa at a resolution of 2.3–2.7 Å, to elucidate structure‐function relationships in NOR, and compare them to those of cytochrome c oxidase (CCO) that is evolutionarily related to NOR. Comprehensive crystallographic refinement of the CO‐bound form of cNOR suggested that a total of four atoms can be accommodated at the binuclear center. Consistent with this, binding of bulky acetaldoxime (CH₃‐CH=N‐OH) to the binuclear center of cNOR was confirmed by the structural analysis. Active site reduction and ligand binding in cNOR induced only ～0.5 Å increase in the heme/non‐heme Fe distance, but no significant structural change in the protein. The highly localized structural change is consistent with the lack of proton‐pumping activity in cNOR, because redox‐coupled conformational changes are thought to be crucial for proton pumping in CCO. It also permits the rapid decomposition of cytotoxic NO in denitrification. In addition, the shorter heme/non‐heme Fe distance even in the bulky ligand‐bound form of cNOR (～4.5 Å) than the heme/Cu distance in CCO (～5 Å) suggests the ability of NOR to maintain two NO molecules within a short distance in the confined space of the active site, thereby facilitating N‐N coupling to produce a hyponitrite intermediate for the generation of N₂O. Proteins 2014; 82:1258–1271. © 2013 Wiley Periodicals, Inc.     

Molecular dynamics simulation study of the binding of purine bases to the aptamer domain of the guanine sensing riboswitch

Riboswitches are a novel class of genetic control elements that function through the direct interaction of small metabolite molecules with structured RNA elements. The ligand is bound with high specificity and affinity to its RNA target and induces conformational changes of the RNA''s secondary and tertiary structure upon binding. To elucidate the molecular basis of the remarkable ligand selectivity and affinity of one of these riboswitches, extensive all-atom molecular dynamics simulations in explicit solvent (≈1 μs total simulation length) of the aptamer domain of the guanine sensing riboswitch are performed. The conformational dynamics is studied when the system is bound to its cognate ligand guanine as well as bound to the non-cognate ligand adenine and in its free form. The simulations indicate that residue U51 in the aptamer domain functions as a general docking platform for purine bases, whereas the interactions between C74 and the ligand are crucial for ligand selectivity. These findings either suggest a two-step ligand recognition process, including a general purine binding step and a subsequent selection of the cognate ligand, or hint at different initial interactions of cognate and noncognate ligands with residues of the ligand binding pocket. To explore possible pathways of complex dissociation, various nonequilibrium simulations are performed which account for the first steps of ligand unbinding. The results delineate the minimal set of conformational changes needed for ligand release, suggest two possible pathways for the dissociation reaction, and underline the importance of long-range tertiary contacts for locking the ligand in the complex.     

A novel mechanism of ligand binding and release in the odorant binding protein 20 from the malaria mosquito Anopheles gambiae

Anopheles gambiae mosquitoes that transmit malaria are attracted to humans by the odor molecules that emanate from skin and sweat. Odorant binding proteins (OBPs) are the first component of the olfactory apparatus to interact with odorant molecules, and so present potential targets for preventing transmission of malaria by disrupting the normal olfactory responses of the insect. AgamOBP20 is one of a limited subset of OBPs that it is preferentially expressed in female mosquitoes and its expression is regulated by blood feeding and by the day/night light cycles that correlate with blood‐feeding behavior. Analysis of AgamOBP20 in solution reveals that the apo‐protein exhibits significant conformational heterogeneity but the binding of odorant molecules results in a significant conformational change, which is accompanied by a reduction in the conformational flexibility present in the protein. Crystal structures of the free and bound states reveal a novel pathway for entrance and exit of odorant molecules into the central‐binding pocket, and that the conformational changes associated with ligand binding are a result of rigid body domain motions in α‐helices 1, 4, and 5, which act as lids to the binding pocket. These structures provide new insights into the specific residues involved in the conformational adaptation to different odorants and have important implications in the selection and development of reagents targeted at disrupting normal OBP function.     

Operating Mechanism and Molecular Dynamics of Pheromone-Binding Protein ASP1 as Influenced by pH

Odorant binding protein (OBP) is a vital component of the olfactory sensation system. It performs the specific role of ferrying odorant molecules to odorant receptors. OBP helps insects and types of animal to sense and transport stimuli molecules. However, the molecular details about how OBPs bind or release its odorant ligands are unclear. For some OBPs, the systems'' pH level is reported to impact on the ligands'' binding or unbinding capability. In this work we investigated the operating mechanism and molecular dynamics in bee antennal pheromone-binding protein ASP1 under varying pH conditions. We found that conformational flexibility is the key factor for regulating the interaction of ASP1 and its ligands, and the odorant binds to ASP1 at low pH conditions. Dynamics, once triggered by pH changes, play the key roles in coupling the global conformational changes with the odorant release. In ASP1, the C-terminus, the N-terminus, helix α2 and the region ranging from helices α4 to α5 form a cavity with a novel ‘entrance’ of binding. These are the major regions that respond to pH change and regulate the ligand release. Clearly there are processes of dynamics and hydrogen bond network propagation in ASP1 in response to pH stimuli. These findings lead to an understanding of the mechanism and dynamics of odorant-OBP interaction in OBP, and will benefit chemsensory-related biotech and agriculture research and development.     

On the link between conformational changes,ligand binding and heat capacity

BackgroundConformational changes coupled to ligand binding constitute the structural and energetics basis underlying cooperativity, allostery and, in general, protein regulation. These conformational rearrangements are associated with heat capacity changes. ITC is a unique technique for studying binding interactions because of the simultaneous determination of the binding affinity and enthalpy, and for providing the best estimates of binding heat capacity changes.Scope of reviewStill controversial issues in ligand binding are the discrimination between the “conformational selection model” and the “induced fit model”, and whether or not conformational changes lead to temperature dependent apparent binding heat capacities. The assessment of conformational changes associated with ligand binding by ITC is discussed. In addition, the “conformational selection” and “induced fit” models are reconciled, and discussed within the context of intrinsically (partially) unstructured proteins.Major conclusionsConformational equilibrium is a major contribution to binding heat capacity changes. A simple model may explain both conformational selection and induced fit scenarios. A temperature-independent binding heat capacity does not necessarily indicate absence of conformational changes upon ligand binding. ITC provides information on the energetics of conformational changes associated with ligand binding (and other possible additional coupled equilibria).General significancePreferential ligand binding to certain protein states leads to an equilibrium shift that is reflected in the coupling between ligand binding and additional equilibria. This represents the structural/energetic basis of the widespread dependence of ligand binding parameters on temperature, as well as pH, ionic strength and the concentration of other chemical species. This article is part of a Special Issue entitled Microcalorimetry in the BioSciences — Principles and Applications, edited by Fadi Bou-Abdallah.     

From signal perception to signal transduction: ligand‐induced dimeric switch of DctB sensory domain in solution

Sinorhizobium meliloti DctB is a typical transmembrane sensory histidine kinase, which senses C₄‐dicarboxylic acids (DCA) and regulates the expression of DctA, the DCA transporter. We previously reported the crystal structures of its periplasmic sensory domain (DctBp) in apo and succinate‐bound states, and these structures showed dramatic conformational changes at dimeric level. Here we show a ligand‐induced dimeric switch in solution and a strong correlation between DctBp's dimerization states and the in vivo activities of DctB. Using site‐directed mutagenesis, we identify important determinants for signal perception and transduction. Specifically, we show that the ligand‐binding pocket is essential for DCA‐induced ‘on’ activity of DctB. Mutations at different sections of DctBp's dimerization interface can lock full‐length DctB at either ‘on’ or ‘off’ state, independent of ligand binding. Taken together, these results suggest that DctBp's signal perception and transduction occur through a ‘ligand‐induced dimeric switch’, in which the changes in the dimeric conformations upon ligand binding are responsible for the signal transduction in DctB.     

Principal component analysis of chemical shift perturbation data of a multiple‐ligand‐binding system for elucidation of respective binding mechanism

Chemical shift perturbations (CSPs) in NMR spectra provide useful information about the interaction of a protein with its ligands. However, in a multiple‐ligand‐binding system, determining quantitative parameters such as a dissociation constant (K_d) is difficult. Here, we used a method we named CS‐PCA, a principal component analysis (PCA) of chemical shift (CS) data, to analyze the interaction between bovine β‐lactoglobulin (βLG) and 1‐anilinonaphthalene‐8‐sulfonate (ANS), which is a multiple‐ligand‐binding system. The CSP on the binding of ANS involved contributions from two distinct binding sites. PCA of the titration data successfully separated the CSP pattern into contributions from each site. Docking simulations based on the separated CSP patterns provided the structures of βLG–ANS complexes for each binding site. In addition, we determined the K_d values as 3.42 × 10⁻⁴M² and 2.51 × 10⁻³M for Sites 1 and 2, respectively. In contrast, it was difficult to obtain reliable K_d values for respective sites from the isothermal titration calorimetry experiments. Two ANS molecules were found to bind at Site 1 simultaneously, suggesting that the binding occurs cooperatively with a partial unfolding of the βLG structure. On the other hand, the binding of ANS to Site 2 was a simple attachment without a significant conformational change. From the present results, CS‐PCA was confirmed to provide not only the positions and the K_d values of binding sites but also information about the binding mechanism. Thus, it is anticipated to be a general method to investigate protein–ligand interactions. Proteins 2013. © 2012 Wiley Periodicals, Inc.     

Extending Bell's model: how force transducer stiffness alters measured unbinding forces and kinetics of molecular complexes

Forced unbinding of complementary macromolecules such as ligand-receptor complexes can reveal energetic and kinetic details governing physiological processes ranging from cellular adhesion to drug metabolism. Although molecular-level experiments have enabled sampling of individual ligand-receptor complex dissociation events, disparities in measured unbinding force F_R among these methods lead to marked variation in inferred binding energetics and kinetics at equilibrium. These discrepancies are documented for even the ubiquitous ligand-receptor pair, biotin-streptavidin. We investigated these disparities and examined atomic-level unbinding trajectories via steered molecular dynamics simulations, as well as via molecular force spectroscopy experiments on biotin-streptavidin. In addition to the well-known loading rate dependence of F_R predicted by Bell's model, we find that experimentally accessible parameters such as the effective stiffness of the force transducer k can significantly perturb the energy landscape and the apparent unbinding force of the complex for sufficiently stiff force transducers. Additionally, at least 20% variation in unbinding force can be attributed to minute differences in initial atomic positions among energetically and structurally comparable complexes. For force transducers typical of molecular force spectroscopy experiments and atomistic simulations, this energy barrier perturbation results in extrapolated energetic and kinetic parameters of the complex that depend strongly on k. We present a model that explicitly includes the effect of k on apparent unbinding force of the ligand-receptor complex, and demonstrate that this correction enables prediction of unbinding distances and dissociation rates that are decoupled from the stiffness of actual or simulated molecular linkers.     

Structural parameterization of the binding enthalpy of small ligands

A major goal in ligand and drug design is the optimization of the binding affinity of selected lead molecules. However, the binding affinity is defined by the free energy of binding, which, in turn, is determined by the enthalpy and entropy changes. Because the binding enthalpy is the term that predominantly reflects the strength of the interactions of the ligand with its target relative to those with the solvent, it is desirable to develop ways of predicting enthalpy changes from structural considerations. The application of structure/enthalpy correlations derived from protein stability data has yielded inconsistent results when applied to small ligands of pharmaceutical interest (MW < 800). Here we present a first attempt at an empirical parameterization of the binding enthalpy for small ligands in terms of structural information. We find that at least three terms need to be considered: (1) the intrinsic enthalpy change that reflects the nature of the interactions between ligand, target, and solvent; (2) the enthalpy associated with any possible conformational change in the protein or ligand upon binding; and, (3) the enthalpy associated with protonation/deprotonation events, if present. As in the case of protein stability, the intrinsic binding enthalpy scales with changes in solvent accessible surface areas. However, an accurate estimation of the intrinsic binding enthalpy requires explicit consideration of long-lived water molecules at the binding interface. The best statistical structure/enthalpy correlation is obtained when buried water molecules within 5-7 A of the ligand are included in the calculations. For all seven protein systems considered (HIV-1 protease, dihydrodipicolinate reductase, Rnase T1, streptavidin, pp60c-Src SH2 domain, Hsp90 molecular chaperone, and bovine beta-trypsin) the binding enthalpy of 25 small molecular weight peptide and nonpeptide ligands can be accounted for with a standard error of +/-0.3 kcal x mol(-1).     

Applying conformational selection theory to improve crossdocking efficiency in 3‐phosphoinositide dependent protein kinase‐1

The emerging picture of biomolecular recognition is that of conformational selection followed by induced‐fit. Conformational selection theory states that binding partners exist in various conformations in solution, with binding involving a “selection” between complementary conformers. In this study, we devise a docking protocol that mimics conformational selection in protein–ligand binding and demonstrate that it significantly enhances crossdocking accuracy over Glide's flexible docking protocol, which is widely used in the pharmaceutical industry. Our protocol uses a pregenerated conformational ensemble to simulate ligand flexibility. The ensemble was generated by thorough conformational sampling coupled with conformer minimization. The generated conformers were then rigidly docked in the active site of the protein along with a postdocking minimization step that allows limited induced fit effects to be modeled for the ligand. We illustrate the improved performance of our protocol through crossdocking of 31 ligands to cocomplexed proteins of the kinase 3‐phosphoinositide dependent protein kinase‐1 extracted from the crystal structures 1H1W (ATP bound), 1OKY (staurosporine bound) and 3QD0 (bound to a potent inhibitor). Consistent with conformational selection theory, the performance of our protocol was the best for crossdocking to the cognate protein bound to the natural ligand, ATP. Proteins 2014; 82:436–451. © 2013 Wiley Periodicals, Inc.     

Binding of fluphenazine with human serum albumin in the presence of rutin and quercetin: An evaluation of food‐drug interaction by spectroscopic techniques

The interactions between human serum albumin (HSA) and fluphenazine (FPZ) in the presence or absence of rutin or quercetin were studied by fluorescence, absorption and circular dichroism (CD) spectroscopy and molecular modeling. The results showed that the fluorescence quenching mechanism was static quenching by the formation of an HSA–FPZ complex. Entropy change (ΔS ⁰) and enthalpy change (ΔH ⁰) values were 68.42 J/(mol? K) and ?4.637 kJ/mol, respectively, which indicated that hydrophobic interactions and hydrogen bonds played major roles in the acting forces. The interaction process was spontaneous because the Gibbs free energy change (ΔG ⁰) values were negative. The results of competitive experiments demonstrated that FPZ was mainly located within HSA site I (sub‐domain IIA). Molecular docking results were in agreement with the experimental conclusions of the thermodynamic parameters and competition experiments. Competitive binding to HSA between flavonoids and FPZ decreased the association constants and increased the binding distances of FPZ binding to HSA. The results of absorption, synchronous fluorescence, three‐dimensional fluorescence, and CD spectra showed that the binding of FPZ to HSA caused conformational changes in HSA and simultaneous effects of FPZ and flavonoids induced further HSA conformational changes.     

Escape of a Small Molecule from Inside T4 Lysozyme by Multiple Pathways

The T4 lysozyme L99A mutant is often used as a model system to study small-molecule binding to proteins, but pathways for ligand entry and exit from the buried binding site and the associated protein conformational changes have not been fully resolved. Here, molecular dynamics simulations were employed to model benzene exit from its binding cavity using the weighted ensemble (WE) approach to enhance sampling of low-probability unbinding trajectories. Independent WE simulations revealed four pathways for benzene exit, which correspond to transient tunnels spontaneously formed in previous simulations of apo T4 lysozyme. Thus, benzene unbinding occurs through multiple pathways partially created by intrinsic protein structural fluctuations. Motions of several α-helices and side chains were involved in ligand escape from metastable microstates. WE simulations also provided preliminary estimates of rate constants for each exit pathway. These results complement previous works and provide a semiquantitative characterization of pathway heterogeneity for binding of small molecules to proteins.     

Anchorage‐dependent binding of integrin I‐domain to adhesion ligands

The dynamic interactions between leukocyte integrin receptors and ligands in the vascular endothelium, extracellular matrix, or invading pathogens result in leukocyte adhesion, extravasation, and phagocytosis. This work examined the mechanical strength of the connection between iC3b, a complement component that stimulates phagocytosis, and the ligand‐binding domain, the I‐domain, of integrin α_Mβ₂. Single‐molecule force measurements of α_M I‐domain–iC3b complexes were conducted by atomic force microscope. Strikingly, depending on loading rates, immobilization of the I‐domain via its C‐terminus resulted in a 1.3‐fold to 1.5‐fold increase in unbinding force compared with I‐domains immobilized via the N‐terminus. The force spectra (unbinding force versus loading rate) of the I‐domain–iC3b complexes revealed that the enhanced mechanical strength is due to a 2.4‐fold increase in the lifetime of the I‐domain–iC3b bond. Given the structural and functional similarity of all integrin I‐domains, our result supports the existing allosteric regulatory model by which the ligand binding strength of integrin can be increased rapidly when a force is allowed to stretch the C‐terminus of the I‐domain. This type of mechanism may account for the rapid ligand affinity adjustment during leukocyte migration. Copyright © 2015 John Wiley & Sons, Ltd.     

First structural evidence for the mode of diffusion of aromatic ligands and ligand-induced closure of the hydrophobic channel in heme peroxidases

The mode of binding of aromatic ligands in the substrate binding site on the distal heme side in heme peroxidases is well understood. However, the mode of diffusion through the extended hydrophobic channel and the regulatory role of the channel are not yet clear. To provide answers to these questions, the crystal structure of the complex of lactoperoxidase and 3-amino-1,2,4-triazole (amitrole) has been determined, which revealed the presence of two ligand molecules, one in the substrate binding site and the second in the hydrophobic channel. The binding of ligand in the channel induced a remarkable conformational change in the side chain of Phe254, which flips from its original distant position to interact with the trapped ligand in the hydrophobic channel. As a result, the channel is completely blocked so that no ligand can diffuse through it to the substrate binding site. Another amitrole molecule is bound to lactoperoxidase in the substrate binding site by replacing three water molecules, including the crucial iron-bound water molecule, W1. In this arrangement, the amino nitrogen atom of amitrole occupies the position of W1 and interacts directly with ferric iron. As a consequence, it prevents the binding of H₂O₂ to heme iron. Thus, the interactions of amitrole with lactoperoxidase obstruct both the passage of ligands through the hydrophobic channel as well as the binding of H₂O₂. This explains the amitrole toxicity. From binding studies, the dissociation constant (K _d) for amitrole with lactoperoxidase was found to be approximately 5.5 × 10⁻⁷ M, indicating high affinity.