Protein Conformational Changes Are Detected and Resolved Site Specifically by Second-Harmonic Generation

We present here a straightforward, broadly applicable technique for real-time detection and measurement of protein conformational changes in solution. This method is based on tethering proteins labeled with a second-harmonic generation (SHG) active dye to supported lipid bilayers. We demonstrate our method by measuring the conformational changes that occur upon ligand binding with three well-characterized proteins labeled at lysine residues: calmodulin (CaM), maltose-binding protein (MBP), and dihydrofolate reductase (DHFR). We also create a single-site cysteine mutant of DHFR engineered within the Met20 catalytic loop region and study the protein’s structural motion at this site. Using published x-ray crystal structures, we show that the changes in the SHG signals upon ligand binding are the result of structural motions that occur at the labeled sites between the apo and ligand-bound forms of the proteins, which are easily distinguished from each other. In addition, we demonstrate that different magnitudes of the SHG signal changes are due to different and specific ligand-induced conformational changes. Taken together, these data illustrate the potential of the SHG approach for detecting and measuring protein conformational changes for a wide range of biological applications.     

Probing mechanisms of resistance to the tuberculosis drug isoniazid: Conformational changes caused by inhibition of InhA, the enoyl reductase from Mycobacterium tuberculosis

The frontline tuberculosis drug isoniazid (INH) inhibits InhA, the NADH-dependent fatty acid biosynthesis (FAS-II) enoyl reductase from Mycobacterium tuberculosis (MTB), via formation of a covalent adduct with NAD(+) (the INH-NAD adduct). Resistance to INH can be correlated with many mutations in MTB, some of which are localized in the InhA cofactor binding site. While the InhA mutations cause a substantial decrease in the affinity of InhA for NADH, surprisingly the same mutations result in only a small impact on binding of the INH-NAD adduct. Based on the knowledge that InhA interacts in vivo with other components of the FAS-II pathway, we have initiated experiments to determine whether enzyme inhibition results in structural changes that could affect protein-protein interactions involving InhA and how these ligand-induced conformational changes are modulated in the InhA mutants. Significantly, while NADH binding to wild-type InhA is hyperbolic, the InhA mutants bind the cofactor with positive cooperativity, suggesting that the mutations permit access to a second conformational state of the protein. While cross-linking studies indicate that enzyme inhibition causes dissociation of the InhA tetramer into dimers, analytical ultracentrifugation and size exclusion chromatography reveal that ligand binding causes a conformational change in the protein that prevents cross-linking across one of the dimer-dimer interfaces in the InhA tetramer. Interestingly, a similar ligand-induced conformational change is also observed for the InhA mutants, indicating that the mutations modulate communication between the subunits without affecting the two conformational states of the protein that are present.     

Inhibitors of different structure induce distinguishing conformations in the omega loop,Cys69-Cys96, of mouse acetylcholinesterase

We have shown previously that association of reversible active site ligands induces a conformational change in an omega loop (Omega loop), Cys(69)-Cys(96), of acetylcholinesterase. The fluorophore acrylodan, site-specifically incorporated at positions 76, 81, and 84, on the external portion of the loop not lining the active site gorge, shows changes in its fluorescence spectrum that reflect the fluorescent side chain moving from a hydrophobic environment to become more solvent-exposed. This appears to result from a movement of the Omega loop accompanying ligand binding. We show here that the loop is indeed flexible and responds to conformational changes induced by both active center and peripheral site inhibitors (gallamine and fasciculin). Moreover, phosphorylation and carbamoylation of the active center serine shows distinctive changes in acrylodan fluorescence spectra at the Omega loop sites, depending on the chirality and steric dimensions of the covalently conjugated ligand. Capping of the gorge with fasciculin, although it does not displace the bound ligand, dominates in inducing a conformational change in the loop. Hence, the ligand-induced conformational changes are distinctive and suggest multiple loop conformations accompany conjugation at the active center serine. The fluorescence changes induced by the modified enzyme may prove useful in the detection of organophosphates or exposure to cholinesterase inhibitors.     

Multistate approaches in computational protein design

Computational protein design (CPD) is a useful tool for protein engineers. It has been successfully applied towards the creation of proteins with increased thermostability, improved binding affinity, novel enzymatic activity, and altered ligand specificity. Traditionally, CPD calculations search and rank sequences using a single fixed protein backbone template in an approach referred to as single-state design (SSD). While SSD has enjoyed considerable success, certain design objectives require the explicit consideration of multiple conformational and/or chemical states. Cases where a "multistate" approach may be advantageous over the SSD approach include designing conformational changes into proteins, using native ensembles to mimic backbone flexibility, and designing ligand or oligomeric association specificities. These design objectives can be efficiently tackled using multistate design (MSD), an emerging methodology in CPD that considers any number of protein conformational or chemical states as inputs instead of a single protein backbone template, as in SSD. In this review article, recent examples of the successful design of a desired property into proteins using MSD are described. These studies employing MSD are divided into two categories-those that utilized multiple conformational states, and those that utilized multiple chemical states. In addition, the scoring of competing states during negative design is discussed as a current challenge for MSD.     

Conformational selection in silico: loop latching motions and ligand binding in enzymes

Ligand binding frequently induces significant conformational changes in a protein receptor. Understanding and predicting such conformational changes represent an important challenge for computational biology, including applications to structure-based drug design. We describe an approach to this problem based on the assumption that the holo state is at least transiently populated in the absence of a ligand; this hypothesis has been referred to as "conformational selection." Here, we apply a method that tests this hypothesis on a challenging class of ligand-induced conformational changes, which we refer to as loop latching: the closing of a loop around an active site that sequesters the ligand from solvent. The method uses a combination of replica exchange molecular dynamics and a loop prediction algorithm to generate low-energy loop structures, and docking to select the conformation appropriate for binding a particular ligand. On a test set of six proteins, it yields loop structures including hololike conformations, generally below 2 A RMSD from the liganded structure, for loops that span up to 15 residues. Docking serves as a stringent test of the predictions. In five of the six cases, the predicted loop conformations improve the ranks of cognate ligands relative to using the apo structure, although the results remain, in most cases, significantly worse than using a holo structure. The poses of the cognate ligands are correct in four of the six test cases, while they are correct for five of the six using a holo structure.     

The translocation mechanism of P-glycoprotein

Multidrug transporters are involved in mediating the failure of chemotherapy in treating several serious diseases. The archetypal multidrug transporter P-glycoprotein (P-gp) confers resistance to a large number of chemically and functionally unrelated anti-cancer drugs by mediating efflux from cancer cells. The ability to efflux such a large number of drugs remains a biological enigma and the lack of mechanistic understanding of the translocation pathway used by P-gp prevents rational design of compounds to inhibit its function. The translocation pathway is critically dependent on ATP hydrolysis and drug interaction with P-gp is possible at one of a multitude of allosterically linked binding sites. However, aspects such as coupling stoichiometry, molecular properties of binding sites and the nature of conformational changes remain unresolved or the centre of considerable controversy. The present review attempts to utilise the available data to generate a detailed sequence of events in the translocation pathway for this dexterous protein.     

Profiling the Protein Targets of Unmodified Bio‐Active Molecules with Drug Affinity Responsive Target Stability and Liquid Chromatography/Tandem Mass Spectrometry

Identifying the target proteins of bioactive small molecules is a key step in understanding mode‐of‐action of the drug and addressing the underlying mechanisms responsible for a particular phenotype. Proteomics has been successfully used to elucidate the target protein profiles of unmodified and ligand‐modified bioactive small molecules. In the latter approach, compounds can be modified via click chemistry and combined with activity‐based protein profiling. Target proteins are then enriched by performing a pull‐down with the modified ligand. Methods that utilize unmodified bioactive small molecules include the cellular thermal shift assay, thermal proteome profiling, stability of proteins from rates of oxidation, and the drug affinity responsive target stability (DARTS) determination (or read‐out). This review highlights recent proteomic approaches utilizing data‐dependent analysis and data‐independent analysis to identify target proteins by DARTS. When combined with liquid chromatography/tandem mass spectrometry, DARTS enables the identification of proteins that bind to drug molecules that leads to a conformational change in the target protein(s). In addition, an effective strategy is proposed for selecting the target protein(s) from within the pool of analyzed candidates. With additional complementary methods, the biologically relevant target proteins that bind to the small bio‐active molecules can be further validated.     

Analysis of Binding Interaction of Curcumin and Diacetylcurcumin with Human and Bovine Serum Albumin Using Fluorescence and Circular Dichroism Spectroscopy

The current study reports the binding of curcumin (CUR) as the main pharmacologically active ingredient of turmeric and diacetylcurcumin (DAC) as a bioactive derivative of curcumin to human serum albumin (HSA) and bovine serum albumin (BSA). The apparent binding constants and number of substantive binding sites have been evaluated by fluorescence quenching method. The distance (r) between donor (HSA and BSA) and acceptor (CUR and DAC) was obtained on the basis of the Förster’s theory of non-radiative energy transfer. The minor changes on the far-UV circular dichroism spectra resulted in partial changes in the calculated secondary structure contents of HSA and BSA. The negligible alteration in the secondary structure of both albumin proteins indicated that ligand-induced conformational changes are localized to the binding site and do not involve considerable changes in protein folding. The visible CD spectra indicated that the optical activity observed during the ligand binding due to induced-protein chirality. All of the achieved results suggested the important role of the phenolic OH group of CUR in the binding process.     

Unusually stable and long-lived ligand-induced conformations of integrins

Integrins are a large family of cell surface receptors that are involved in a wide range of biological processes. The integrin alpha(IIb)beta3 (glycoprotein IIb-IIIa) is a major platelet glycoprotein heterodimeric receptor that mediates platelet aggregation and is currently a target for pharmaceutical intervention. Ligand binding to the receptor has been shown to induce conformational changes by physical methods and the exposure of neoepitopes (the ligand-induced binding sites). Here we show that the antagonist XP280 induces a conformation that is stable to treatment with SDS and that the protein retains this conformation for several days even after dissociation of the inhibitor. These ligand-induced conformational changes take place with purified protein and on intact platelets. They are competable with an RGDS peptide and are stable to reduction but not boiling or treatment with EDTA. The retention of an altered conformation in the absence of the ligand implies the possibility of ligand-induced alteration of biological function even in the absence of ligand. Finally, similar behavior is observed with the integrin alpha(v)beta3, suggesting that access to SDS stable conformations may be conserved throughout the integrin superfamily. The unusual stability, long-lived nature, and potential generality of these conformations could have profound implications for integrin biology.     

Ligand-induced conformational changes in a thermophilic ribose-binding protein

Background

Members of the periplasmic binding protein (PBP) superfamily are involved in transport and signaling processes in both prokaryotes and eukaryotes. Biological responses are typically mediated by ligand-induced conformational changes in which the binding event is coupled to a hinge-bending motion that brings together two domains in a closed form. In all PBP-mediated biological processes, downstream partners recognize the closed form of the protein. This motion has also been exploited in protein engineering experiments to construct biosensors that transduce ligand binding to a variety of physical signals. Understanding the mechanistic details of PBP conformational changes, both global (hinge bending, twisting, shear movements) and local (rotamer changes, backbone motion), therefore is not only important for understanding their biological function but also for protein engineering experiments.

Results

Here we present biochemical characterization and crystal structure determination of the periplasmic ribose-binding protein (RBP) from the hyperthermophile Thermotoga maritima in its ribose-bound and unliganded state. The T. maritima RBP (tmRBP) has 39% sequence identity and is considerably more resistant to thermal denaturation ( ^app T _m value is 108°C) than the mesophilic Escherichia coli homolog (ecRBP) ( ^app T _m value is 56°C). Polar ligand interactions and ligand-induced global conformational changes are conserved among ecRBP and tmRBP; however local structural rearrangements involving side-chain motions in the ligand-binding site are not conserved.

Conclusion

Although the large-scale ligand-induced changes are mediated through similar regions, and are produced by similar backbone movements in tmRBP and ecRBP, the small-scale ligand-induced structural rearrangements differentiate the mesophile and thermophile. This suggests there are mechanistic differences in the manner by which these two proteins bind their ligands and are an example of how two structurally similar proteins utilize different mechanisms to form a ligand-bound state.     

Ligand-induced conformational changes within a hexameric Acyl-CoA thioesterase

Acyl-coenzyme A (acyl-CoA) thioesterases play a crucial role in the metabolism of activated fatty acids, coenzyme A, and other metabolic precursor molecules including arachidonic acid and palmitic acid. These enzymes hydrolyze coenzyme A from acyl-CoA esters to mediate a range of cellular functions including β-oxidation, lipid biosynthesis, and signal transduction. Here, we present the crystal structure of a hexameric hot-dog domain-containing acyl-CoA thioesterase from Bacillus halodurans in the apo-form and provide structural and comparative analyses to the coenzyme A-bound form to identify key conformational changes induced upon ligand binding. We observed dramatic ligand-induced changes at both the hot-dog dimer and the trimer-of-dimer interfaces; the dimer interfaces in the apo-structure differ by over 20% and decrease to about half the size in the ligand-bound state. We also assessed the specificity of the enzyme against a range of fatty acyl-CoA substrates and have identified a preference for short-chain fatty acyl-CoAs. Coenzyme A was shown both to negatively regulate enzyme activity, representing a direct inhibitory feedback, and consistent with the structural data, to destabilize the quaternary structure of the enzyme. Coenzyme A-induced conformational changes in the C-terminal helices of enzyme were assessed through mutational analysis and shown to play a role in regulating enzyme activity. The conformational changes are likely to be conserved from bacteria through to humans and provide a greater understanding, particularly at a structural level, of thioesterase function and regulation.     

Allosteric sites can be identified based on the residue–residue interaction energy difference

Allosteric drugs act at a distance to regulate protein functions. They have several advantages over conventional orthosteric drugs, including diverse regulation types and fewer side effects. However, the rational design of allosteric ligands remains a challenge, especially when it comes to the identification allosteric binding sites. As the binding of allosteric ligands may induce changes in the pattern of residue–residue interactions, we calculated the residue–residue interaction energies within the allosteric site based on the molecular mechanics generalized Born surface area energy decomposition scheme. Using a dataset of 17 allosteric proteins with structural data for both the apo and the ligand‐bound state available, we used conformational ensembles generated by molecular dynamics simulations to compute the differences in the residue–residue interaction energies in known allosteric sites from both states. For all the known sites, distinct interaction energy differences (>25%) were observed. We then used CAVITY, a binding site detection program to identify novel putative allosteric sites in the same proteins. This yielded a total of 31 “druggable binding sites,” of which 21 exhibited >25% difference in residue interaction energies, and were hence predicted as novel allosteric sites. Three of the predicted allosteric sites were supported by recent experimental studies. All the predicted sites may serve as novel allosteric sites for allosteric ligand design. Our study provides a computational method for identifying novel allosteric sites for allosteric drug design. Proteins 2015; 83:1375–1384. © 2014 Wiley Periodicals, Inc.     

Induced fit or conformational selection? The role of the semi-closed state in the maltose binding protein

A full characterization of the thermodynamic forces underlying ligand-associated conformational changes in proteins is essential for understanding and manipulating diverse biological processes, including transport, signaling, and enzymatic activity. Recent experiments on the maltose binding protein (MBP) have provided valuable data about the different conformational states implicated in the ligand recognition process; however, a complete picture of the accessible pathways and the associated changes in free energy remains elusive. Here we describe results from advanced accelerated molecular dynamics (aMD) simulations, coupled with adaptively biased force (ABF) and thermodynamic integration (TI) free energy methods. The combination of approaches allows us to track the ligand recognition process on the microsecond time scale and provides a detailed characterization of the protein's dynamic and the relative energy of stable states. We find that an induced-fit (IF) mechanism is most likely and that a mechanism involving both a conformational selection (CS) step and an IF step is also possible. The complete recognition process is best viewed as a "Pac Man" type action where the ligand is initially localized to one domain and naturally occurring hinge-bending vibrations in the protein are able to assist the recognition process by increasing the chances of a favorable encounter with side chains on the other domain, leading to a population shift. This interpretation is consistent with experiments and provides new insight into the complex recognition mechanism. The methods employed here are able to describe IF and CS effects and provide formally rigorous means of computing free energy changes. As such, they are superior to conventional MD and flexible docking alone and hold great promise for future development and applications to drug discovery.     

Autoinhibited proteins as promising drug targets

Current drug discovery efforts generally focus on a limited number of protein classes, typically including proteins with well-defined catalytic active sites (e.g., kinases) or ligand binding sites (e.g., G protein-coupled receptors). Nevertheless, many clinically important pathways are mediated by proteins with no such obvious targets for small molecule inhibitors. Allosteric inhibitors offer an alternative approach to inhibition of protein activities, particularly for proteins that undergo conformational changes as part of their activity cycle. Proteins regulated by autoinhibitory domains represent one broad class of proteins that meets this criterion. In this article, we discuss the potential of autoinhibited proteins as targets for allosteric inhibitors and describe two examples of small molecules that act by stabilizing native autoinhibited conformations of their targets. We propose that proteins regulated by autoinhibition may be generally amenable to allosteric inhibition by small molecules that stabilize the native, autoinhibited fold.     

Construction and optimization of a family of genetically encoded metabolite sensors by semirational protein engineering

A family of genetically-encoded metabolite sensors has been constructed using bacterial periplasmic binding proteins (PBPs) linearly fused to protein fluorophores. The ligand-induced conformational change in a PBP allosterically regulates the relative distance and orientation of a fluorescence resonance energy transfer (FRET)-compatible protein pair. Ligand binding is transduced into a macroscopic FRET observable, providing a reagent for in vitro and in vivo ligand-measurement and visualization. Sensors with a higher FRET signal change are required to expand the dynamic range and allow visualization of subtle analyte changes under high noise conditions. Various observations suggest that factors other than inter-fluorophore separation contribute to FRET transfer efficiency and the resulting ligand-dependent spectral changes. Empirical and rational protein engineering leads to enhanced allosteric linkage between ligand binding and chromophore rearrangement; modifications predicted to decrease chromophore rotational averaging enhance the signal change, emphasizing the importance of the rotational freedom parameter kappa2 to FRET efficiency. Tighter allosteric linkage of the PBP and the fluorophores by linker truncation or by insertion of chromophores into the binding protein at rationally designed sites gave rise to sensors with improved signal change. High-response sensors were obtained with fluorescent proteins attached to the same binding PBP lobe, suggesting that indirect allosteric regulation during the hinge-bending motion is sufficient to give rise to a FRET response. The optimization of sensors for glucose and glutamate, ligands of great clinical interest, provides a general framework for the manipulation of ligand-dependent allosteric signal transduction mechanisms.     

DynaBiS: A hierarchical sampling algorithm to identify flexible binding sites for large ligands and peptides

Knowing the ligand or peptide binding site in proteins is highly important to guide drug discovery, but experimental elucidation of the binding site is difficult. Therefore, various computational approaches have been developed to identify potential binding sites in protein structures. However, protein and ligand flexibility are often neglected in these methods due to efficiency considerations despite the recognition that protein–ligand interactions can be strongly affected by mutual structural adaptations. This is particularly true if the binding site is unknown, as the screening will typically be performed based on an unbound protein structure. Herein we present DynaBiS, a hierarchical sampling algorithm to identify flexible binding sites for a target ligand with explicit consideration of protein and ligand flexibility, inspired by our previously presented flexible docking algorithm DynaDock. DynaBiS applies soft-core potentials between the ligand and the protein, thereby allowing a certain protein–ligand overlap resulting in efficient sampling of conformational adaptation effects. We evaluated DynaBiS and other commonly used binding site identification algorithms against a diverse evaluation set consisting of 26 proteins featuring peptide as well as small ligand binding sites. We show that DynaBiS outperforms the other evaluated methods for the identification of protein binding sites for large and highly flexible ligands such as peptides, both with a holo or apo structure used as input.     

High yield bacterial expression of active c-Abl and c-Src tyrosine kinases

The Abl and Src tyrosine kinases are key signaling proteins that are of considerable interest as drug targets in cancer and many other diseases. The regulatory mechanisms that control the activity of these proteins are complex, and involve large-scale conformational changes in response to phosphorylation and other modulatory signals. The success of the Abl inhibitor imatinib in the treatment of chronic myelogenous leukemia has shown the potential of kinase inhibitors, but the rise of drug resistance in patients has also shown that drugs with alternative modes of binding to the kinase are needed. The detailed understanding of mechanisms of protein-drug interaction and drug resistance through biophysical methods demands a method for the production of active protein on the milligram scale. We have developed a bacterial expression system for the kinase domains of c-Abl and c-Src, which allows for the quick expression and purification of active wild-type and mutant kinase domains by coexpression with the YopH tyrosine phosphatase. This method makes practical the use of isotopic labeling of c-Abl and c-Src for NMR studies, and is also applicable for constructs containing the SH2 and SH3 domains of the kinases.     

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.     

Computational Study of Synthetic Agonist Ligands of Ionotropic Glutamate Receptors

Neurological glutamate receptors are among the most important and intensely studied protein ligand binding systems in humans. They are crucial for the functioning of the central nervous system and involved in a variety of pathologies. Apart from the neurotransmitter glutamate, several artificial, agonistic and antagonistic ligands are known. Of particular interest here are novel photoswitchable agonists that would open the field of optogenetics to glutamate receptors. The receptor proteins are complex, membrane-bound multidomain oligomers that undergo large scale functional conformational changes, making detailed studies of their atomic structure challenging. Therefore, a thorough understanding of the microscopic details of ligand binding and receptor activation remains elusive in many cases. This topic has been successfully addressed by theoretical studies in the past and in this paper, we present extensive molecular dynamics simulation and free energy calculation results on the binding of AMPA and an AMPA derivative, which is the basis for designing light-sensitive ligands. We provide a two-step model for ligand binding domain activation and predict binding free energies for novel compounds in good agreement to experimental observations.