Binding Site Flexibility Research Articles

A Specific Two-pore Domain Potassium Channel Blocker Defines the Structure of the TASK-1 Open Pore

Two-pore domain potassium (K2P) channels play a key role in setting the membrane potential of excitable cells. Despite their role as putative targets for drugs and general anesthetics, little is known about the structure and the drug binding site of K2P channels. We describe A1899 as a potent and highly selective blocker of the K2P channel TASK-1. As A1899 acts as an open-channel blocker and binds to residues forming the wall of the central cavity, the drug was used to further our understanding of the channel pore. Using alanine mutagenesis screens, we have identified residues in both pore loops, the M2 and M4 segments, and the halothane response element to form the drug binding site of TASK-1. Our experimental data were used to validate a K2P open-pore homology model of TASK-1, providing structural insights for future rational design of drugs targeting K2P channels.

From Crystal Packing to Molecular Recognition: Prediction and Discovery of a Binding Site on the Surface of Polo‐Like Kinase 1

Picking the pocket: Crystal-packing interactions are present in every macromolecular structure solved with X-ray methods; however, they have not been extensively used to gain insight into molecular recognition. A novel approach is described to utilize these packing interactions for the discovery of a flexible binding site and ligands recognized by it (see picture).

QM/MM study of the structure, energy storage, and origin of the bathochromic shift in vertebrate and invertebrate bathorhodopsins.

By comparing the results from a hybrid quantum mechanics/molecular mechanics method (SORCI+Q//B3LYP/6-31G*:Amber) between vertebrate (bovine) and invertebrate (squid) visual pigments, the mechanism of molecular rearrangements, energy storage, and origin of the bathochromic shift accompanying the transformation of rhodopsin to bathorhodopsin have been evaluated. The analysis reveals that, in the presence of an unrelaxed binding site, bathorhodopsin was found to carry almost 27 kcal/mol energy in both visual pigments and absorb (λ(max)) at 528 nm in bovine and 554 nm in squid. However, when the residues within 4.0 Å radius of the retinal are relaxed during the isomerization event, almost ∼16 kcal/mol energy is lost in squid compared to only ∼8 kcal/mol in bovine. Loss of a larger amount of energy in squid is attributed to the presence of a flexible binding site compared to a rigid binding site in bovine. Structure of the squid bathorhodopsin is characterized by formation of a direct H-bond between the Schiff base and Asn87.

Exploring RNA Drug Binding Using Consistent Charge Models

RNA has emerged as a prospective drug target for a wide range of diseases. However, existing drug discovery tools optimized for protein targets have been largely unsuccessful in producing novel small molecules that target RNA due to RNAs distinct chemical and structural characteristics. Prevalent screening methods assay catalytic activity and are therefore unsuitable for RNA, as few RNAs are catalytically active. Furthermore, RNA has flexible binding sites, precluding the causal relationship between strong binding and inhibition of activity. Computational methods, including molecular docking, can overcome some of these limitations; however, RNA can adopt radically different conformations upon binding small molecules which can be difficult to model computationally. An additional challenge is that current force fields remain underdeveloped in modeling polyanionic nucleic acids with complex electrostatic interactions. Accurately capturing these interactions is crucial to determining precise free energy calculations of binding. We have developed a force field extension model that allows the same charges and force field parameters to be used for both the receptor and ligand, significantly improving the accuracy of capturing interactions between receptor RNAs and potential drug-like small molecules. Here we present details on our force field extension model, predicting the hydration free energy of 503 organic molecules using free energy perturbation demonstrating that our model reproduces experimentally confirmed results. Furthermore, we present docking results of predicted binding affinities and ligand-bound poses for a set of 60 known RNA-ligand complexes from the Protein Databank and compare results with other docking programs. Finally, we predict the binding affinities of 65 small molecules to the HIV-1 TAR RNA and compare results with a high throughput screen.

Roles of Residues Cys69, Asn104, Phe160, Gly232, Ser237, and Asp240 in Extended-Spectrum β-Lactamase Toho-1

Toho-1, which is also designated CTX-M-44, is an extended-spectrum class A β-lactamase that has high activity toward cefotaxime. In this study, we investigated the roles of residues suggested to be critical for the substrate specificity expansion of Toho-1 in previous structural analyses. Six amino acid residues were replaced one by one with amino acids that are often observed in the corresponding position of non-extended-spectrum β-lactamases. The mutants produced in Escherichia coli strains were analyzed both for their kinetic properties and their effect on drug susceptibilities. The results indicate that the substitutions of Asn104 and Ser237 have certain effects on expansion of substrate specificity, while those of Cys69 and Phe160 have less effect, and that of Asp240 has no effect on the hydrolysis of any substrates tested. Gly232, which had been assumed to increase the flexibility of the substrate binding site, was revealed not to be critical for the expansion of substrate specificity of this enzyme, although this substitution resulted in deleterious effects on expression and stability of the enzyme.

A Novel p53 Phosphorylation Site within the MDM2 Ubiquitination Signal: II. A MODEL IN WHICH PHOSPHORYLATION AT SER269 INDUCES A MUTANT CONFORMATION TO p53

The p53 DNA-binding domain harbors a conformationally flexible multiprotein binding site that regulates p53 ubiquitination. A novel phosphorylation site exists within this region at Ser(269), whose phosphomimetic mutation inactivates p53. The phosphomimetic p53 (S269D) exhibits characteristics of mutant p53: stable binding to Hsp70 in vivo, elevated ubiquitination in vivo, inactivity in DNA binding and transcription, increased thermoinstability using thermal shift assays, and λ(max) of intrinsic tryptophan fluorescence at 403 nm rather than 346 nm, characteristic of wild type p53. These data indicate that p53 conformational stability is regulated by a phosphoacceptor site within an exposed flexible surface loop and that this can be destabilized by phosphorylation. To test whether other motifs within p53 have similarly evolved, we analyzed the effect of Ser(215) mutation on p53 function because Ser(215) is another inactivating phosphorylation site in the conformationally flexible PAb240 epitope. The p53(S215D) protein is inactive like p53(S269D), whereas p53(S215A) is as active as p53(S269A). However, the double mutant p53(S215A/S269A) was transcriptionally inactive and more thermally unstable than either individual Ser-Ala loop mutant. Molecular dynamics simulations suggest that (i) solvation of phospho-Ser(215) and phospho-Ser(269) by positive charged residues or solvent water leads to local unfolding, which is accompanied by local destabilization of the N-terminal loop and global destabilization of p53, and (ii) the double alanine 215/269 mutation disrupts hydrogen bonding normally stabilized by both Ser(215) and Ser(269). These data indicate that p53 has evolved two serine phosphoacceptor residues within conformationally flexible epitopes that normally stabilize the p53 DNA-binding domain but whose phosphorylation induces a mutant conformation to wild type p53.

A Novel p53 Phosphorylation Site within the MDM2 Ubiquitination Signal: I. PHOSPHORYLATION AT SER269 IN VIVO IS LINKED TO INACTIVATION OF p53 FUNCTION

p53 is a thermodynamically unstable protein containing a conformationally flexible multiprotein docking site within the DNA-binding domain. A combinatorial peptide chip used to identify the novel kinase consensus site RXSΦ(K/D) led to the discovery of a homologous phosphorylation site in the S10 β-strand of p53 at Ser(269). Overlapping peptide libraries confirmed that Ser(269) was a phosphoacceptor site in vitro, and immunochemical approaches evaluated whether p53 is phosphorylated in vivo at Ser(269). Mutation or phosphorylation of p53 at Ser(269) attenuates binding of the p53-specific monoclonal antibody DO-12, identifying an assay for measuring Ser(269) phosphorylation of p53 in vivo. The mAb DO-12 epitope of p53 is masked via phosphorylation in a range of human tumor cells with WT p53 status, as defined by increased mAb DO-12 binding to endogenous p53 after phosphatase treatment. Phospho-Ser(269)-specific monoclonal antibodies were generated and used to demonstrate that p53 phosphorylation is induced at Ser(269) after irradiation with kinetics similar to those of p53 protein induction. Phosphomimetic mutation at Ser(269) inactivated the transcription activation function and clonogenic suppressor activity of p53. These data suggest that the dynamic equilibrium between native and unfolded states of WT p53 can be modulated by phosphorylation of the conformationally flexible multiprotein binding site in the p53 DNA-binding domain.

Binding site and ligand flexibility revealed by high resolution crystal structures of GluK1 competitive antagonists

The availability of crystal structures for the ligand binding domains of ionotropic glutamate receptors, combined with their key role in synaptic function in the normal and diseased brain, offers a unique selection of targets for pharmaceutical research compared to other drug targets for which the atomic structure of the ligand binding site is not known. Currently only a few antagonist structures have been solved, and these reveal ligand specific conformational changes that hinder rational drug design. Here we report high resolution crystal structures for three kainate receptor GluK1 antagonist complexes which reveal new and unexpected modes of binding, highlighting the continued need for experimentally determined receptor–ligand complexes.

Improved Ligand-Protein Binding Affinity Predictions Using Multiple Binding Modes

Accurate ligand-protein binding affinity prediction, for a set of similar binders, is a major challenge in the lead optimization stage in drug development. In general, docking and scoring functions perform unsatisfactorily in this application. Docking calculations, followed by molecular dynamics simulations and free energy calculations can be applied to improve the predictions. However, for targets with large, flexible binding sites, with no experimentally determined binding modes for a set of ligands, insufficient sampling can decrease the accuracy of the free energy calculations. Cytochrome P450s, a protein family of major importance for drug metabolism, is an example of a challenging target for binding affinity predictions. As a result, the choice of starting structure from the docking solutions becomes crucial. In this study, an iterative scheme is introduced that includes multiple independent molecular dynamics simulations to obtain weighted ensemble averages to be used in the linear interaction energy method. The proposed scheme makes the initial pose selection less crucial for further simulation, as it automatically calculates the relative weights of the various poses. It also properly takes into account the possibility that multiple binding modes contribute similarly to the overall affinity, or of similar compounds occupying very different poses. The method was applied to a set of 12 compounds binding to cytochrome P450 2C9 and it displayed a root mean-square error of 2.9 kJ/mol.

Exploring the Ion Selectivity Properties of a Large Number of Simplified Binding Site Models

The ability to discriminate between different cations efficiently is essential for the proper physiological functioning of many membrane transport proteins. One obvious mechanism of ion selectivity is when a binding site is structurally constrained by the protein architecture and its geometry is precisely adapted to fit an ion of a given size. This mechanism is not effective in the case of flexible protein binding sites that are able to deform structurally or to adapt to a bound ion. In this study, the concept of nontrivial ion selectivity arising in a highly flexible protein binding site conceptually represented as a microdroplet of ligands confined to a small volume is explored. The environment imposed by the spatial confinement is a critical feature of the reduced models. A large number of reduced binding site models (1077) comprising typical ion-coordinating ligands (carbonyl, hydroxyl, carboxylate, water) are constructed and characterized for Na+/K+ and Ca2+/Ba2+ size selectivity using free energy perturbation molecular dynamics simulations. Free energies are highly correlated with the sum of ion-ligand and ligand-ligand mean interactions, but the relative balance of those two contributions is different for K+-selective and Na+-selective binding sites. The analysis indicates that both the number and the type of ligands are important factors in ion selectivity.

Efficient Virtual Screening Using Multiple Protein Conformations Described as Negative Images of the Ligand-Binding Site

The protein structure-based virtual screening is typically accomplished using a molecular docking procedure. However, docking is a fairly slow process that is limited by the available scoring functions that cannot reliably distinguish between active and inactive ligands. In contrast, the ligand-based screening methods that are based on shape similarity identify the active ligands with high accuracy. Here, we show that the usage of negative images of the ligand-binding site, together with shape comparison tools, which are typically used in ligand-based virtual screening, improve the discrimination of active molecules from inactives. In contrast to ligand-based shape comparison, the negative image of the binding site allows identification of compounds whose shape complements the shape of the ligand-binding cavity as closely as possible. Furthermore, the use of several target protein conformations allows the identification of active ligands whose shape is not optimal for crystallized protein conformation. Accordingly, the presented virtual screening method improves the identification of novel lead molecules by concentrating on the optimally shaped molecules for the flexible ligand binding site.

A flexible cation binding site in the multidrug major facilitator superfamily transporter LmrP is associated with variable proton coupling.

The multidrug major facilitator superfamily transporter LmrP from Lactococcus lactis mediates protonmotive-force dependent efflux of amphiphilic ligands from the cell. We compared the role of membrane-embedded carboxylates in transport and binding of divalent cationic propidium and monovalent cationic ethidium. D235N, E327Q, and D142N replacements each resulted in loss of electrogenicity in the propidium efflux reaction, pointing to electrogenic 3H(+)/propidium(2+) antiport. During ethidium efflux, single D142N and D235N replacements resulted in apparent loss of electrogenicity, whereas the E327Q substitution did not affect the energetics, consistent with electrogenic 2H(+)/ethidium(+) antiport. Different roles of carboxylates were also observed in fluorescence anisotropy-based ligand-binding assays. Whereas D235 and E327 were both involved in propidium binding, the loss of one of these carboxylates could be compensated for by the other in ethidium binding. The D142N replacement did not affect the binding of either ligand. These data point to the presence of a dedicated proton binding site containing D142, and a flexible proton/ligand binding site containing D235 and E327, the contributions to proton and ligand binding of which depend on the chemical structure of the bound ligand. Our findings provide the first evidence that multidrug transport by secondary-active transporters can be associated with variable ion coupling.

Methionine motifs of copper transport proteins provide general and flexible thioether-only binding sites for Cu(I) and Ag(I)

Cellular acquisition of copper in eukaryotic organisms is primarily accomplished through high-affinity copper transport proteins (Ctr). The extracellular N-terminal regions of both human and yeast Ctr1 contain multiple methionine residues organized in copper-binding Mets motifs. These motifs comprise combinations of methionine residues arranged in clusters of MXM and MXXM, where X can be one of several amino acids. Model peptides corresponding to 15 different Mets motifs were synthesized and determined to selectively bind Cu(I) and Ag(I), with no discernible affinity for divalent metal ions. These are rare examples of biological thioether-only metal binding sites. Effective dissociation constant (KD) values for the model Mets peptides and Cu(I) were determined by an ascorbic acid oxidation assay and validated through electrospray ionization mass spectrometry and range between 2 and 11 microM. Affinity appears to be independent of pH, the arrangement of the motif, and the composition of intervening amino acids, all of which reveal the generality and flexibility of the MX1-2MX1-2M domain. Circular dichroism spectroscopy, 1H-NMR spectroscopy, and X-ray absorption spectroscopy were also used to characterize the binding event. These results are intended to aid the development of the still unknown mechanism of copper transport across the cell membrane.

Binding-Site Assessment by Virtual Fragment Screening

The accurate prediction of protein druggability (propensity to bind high-affinity drug-like small molecules) would greatly benefit the fields of chemical genomics and drug discovery. We have developed a novel approach to quantitatively assess protein druggability by computationally screening a fragment-like compound library. In analogy to NMR-based fragment screening, we dock ∼11000 fragments against a given binding site and compute a computational hit rate based on the fraction of molecules that exceed an empirically chosen score cutoff. We perform a large-scale evaluation of the approach on four datasets, totaling 152 binding sites. We demonstrate that computed hit rates correlate with hit rates measured experimentally in a previously published NMR-based screening method. Secondly, we show that the in silico fragment screening method can be used to distinguish known druggable and non-druggable targets, including both enzymes and protein-protein interaction sites. Finally, we explore the sensitivity of the results to different receptor conformations, including flexible protein-protein interaction sites. Besides its original aim to assess druggability of different protein targets, this method could be used to identifying druggable conformations of flexible binding site for lead discovery, and suggesting strategies for growing or joining initial fragment hits to obtain more potent inhibitors.

A flexible docking scheme to explore the binding selectivity of PDZ domains

Modeling of protein binding site flexibility in molecular docking is still a challenging problem due to the large conformational space that needs sampling. Here, we propose a flexible receptor docking scheme: A dihedral restrained replica exchange molecular dynamics (REMD), where we incorporate the normal modes obtained by the Elastic Network Model (ENM) as dihedral restraints to speed up the search towards correct binding site conformations. To our knowledge, this is the first approach that uses ENM modes to bias REMD simulations towards binding induced fluctuations in docking studies. In our docking scheme, we first obtain the deformed structures of the unbound protein as initial conformations by moving along the binding fluctuation mode, and perform REMD using the ENM modes as dihedral restraints. Then, we generate an ensemble of multiple receptor conformations (MRCs) by clustering the lowest replica trajectory. Using ROSETTALIGAND, we dock ligands to the clustered conformations to predict the binding pose and affinity. We apply this method to postsynaptic density-95/Dlg/ZO-1 (PDZ) domains; whose dynamics govern their binding specificity. Our approach produces the lowest energy bound complexes with an average ligand root mean square deviation of 0.36 A. We further test our method on (i) homologs and (ii) mutant structures of PDZ where mutations alter the binding selectivity. In both cases, our approach succeeds to predict the correct pose and the affinity of binding peptides. Overall, with this approach, we generate an ensemble of MRCs that leads to predict the binding poses and specificities of a protein complex accurately.

Comparison of current docking tools for the simulation of inhibitor binding by the transmembrane domain of the sarco/endoplasmic reticulum calcium ATPase

Inhibitors of the transmembrane protein sarco/endoplasmic reticulum calcium ATPase (SERCA) are invaluable tools for the study of the enzyme's physiological functions and they have been recognized as a promising new class of anticancer agents. For the discovery of novel enzyme inhibitors, small molecule docking for virtual screens of large compound libraries has become increasingly important. Since the performance of various docking routines varies considerably, depending on the target and the chemical nature of the ligand, we critically evaluated the performance of four frequently used programs — GOLD, AutoDock, Surflex-Dock, and FRED — for the docking of SERCA inhibitors based on the structures of thapsigargin, di-tert-butylhydroquinone, and cyclopiazonic acid. Evaluation criteria were docking accuracy using crystal structures as references, docking reproducibility, and correlation between docking scores and known bioactivities. The best overall results were obtained by GOLD and FRED. Docking runs with conformationally flexible binding sites produced no significant improvement of the results.

Mechanism of Selective Cation Binding to Sodium-Coupled Transporters: Insights from Free Energy Simulations and QM/MM Simulations

Ion-coupled transport of neurotransmitter molecules by secondary amino-acids transporters plays pivotal role in the regulation of neuronal signaling. One of the major events in the transport cycle is ion-substrate coupling and formation of the high-affinity occluded state with bound ions and substrate. Molecular mechanisms of ion-substrate coupling, specificity for a particular cation and the corresponding ion-substrate stoichiometry in secondary transporters has yet to be understood. We have studied Li+/K+/Tl+/Na+ binding and/or selectivity to several transporters with available crystal structures such as the bacterial aspartate transporter GltPh, leucing transporter LeuT and maltose transporter vSGLT using free energy simulations and QM/MM minimization to evaluate the role of different factors in the observed selectivity and ion binding to the protein. Two different mechanisms were found to co-exist for crystallographically characterized binding sites Na1 and Na2 in LeuT and Glt. Furthermore, site Na1 appeared to be well conserved amongst members of different families. To evaluate the role of Na+ binding in the transporter function, we have performed free energy simulations to determine actual cation selectivity as well as binding affinity for sites Na1 and Na2 in the protein. QM/MM minimization was used to characterize the role of the electronic effects of the stabilization of non-native cations such as Li+ and Tl+ in the Na+-selective sites of different transporters. In the case of Tl+ binding to Glt transporter, neighboring residues from a second solvation shell provide a necessary stabilization to the larger cation due to polarization and charge transfer effects implying a rather large flexibility of the metal binding sites.

Dissecting interactions between nucleosides and germination receptors in Bacillus cereus 569 spores

Bacillus cereus 569 spores germinate either with inosine as a sole germinant or with a combination of nucleosides and L-alanine. Whereas the inosine-only germination pathway requires the presence of two different germination receptors (GerI and GerQ) to be activated, the nucleoside/alanine germination pathway only needs one of the two receptors. To differentiate how nucleoside recognition varies between the inosine-only germination pathway and the nucleoside/alanine germination pathway, we tested 61 purine analogues as agonists and antagonists of the two pathways in wild-type, DeltagerI and DeltagerQ spores. The structure-activity relationships of germination agonists and antagonists suggest that the inosine-only germination pathway is restricted to recognize a single germinant (inosine), but can be inhibited in predictable patterns by structurally distinct purine nucleosides. B. cereus spores encoding GerI as the only nucleoside receptor (DeltagerQ mutant) showed a germination inhibition profile similar to wild-type spores treated with inosine only. Thus, GerI seems to have a well-organized binding site that recognizes inosine and inhibitors through specific substrate-protein interactions. Structure-activity analysis also showed that the nucleoside/alanine germination pathway is more promiscuous toward purine nucleoside agonists, and is only inhibited by hydrophobic analogues. B. cereus spores encoding GerQ as the only nucleoside receptor (DeltagerI mutant) behaved like wild-type spores treated with inosine and L-alanine. Thus, the GerQ receptor seems to recognize substrates in a more flexible binding site through non-specific interactions. We propose that the GerI receptor is responsible for germinant detection in the inosine-only germination pathway. On the other hand, supplementing inosine with l-alanine allows bypassing of the GerI receptor to activate the more flexible GerQ receptor.

Deciphering ligand dependent degree of binding site closure and its implication in inhibitor design: A modeling study on human adenosine kinase

Protein flexibility plays a significant role in drug research due to its effect on accurate prediction of ligand binding mode and activity. Adenosine kinase (AK) represents a highly flexible binding site and is known to exhibit large conformational changes as a result of substrate or inhibitor binding. Here we propose a semi-open conformation for ligand binding in human AK, in addition to the known closed and open forms. The modeling study illustrates the necessity of thorough understanding of the conformational states of protein for docking and binding mode prediction. It has been shown that predicting activity in the context of correct binding mode can improve the insight into conserved interactions and mechanism of action for inhibition of AK. Integrating the knowledge about the binding modes of ligands in different conformational states of the protein, separate pharmacophore models were generated and used for virtual screening to explore potential novel hits. In addition, 2D descriptor based clustering was done to differentiate the ligands, binding to closed, semi-open and open conformations of human AK. The results indicated that binding of all AK inhibitors cannot be described by same rules, instead, they represent a rule based preference for inhibition. This inference about tubercidins binding to semi-open conformation of human AK may facilitate in finding much extensive space for AK inhibitors.

Substrate binding site flexibility of the small heat shock protein molecular chaperones.

Small heat shock proteins (sHSPs) serve as a first line of defense against stress-induced cell damage by binding and maintaining denaturing proteins in a folding-competent state. In contrast to the well-defined substrate binding regions of ATP-dependent chaperones, interactions between sHSPs and substrates are poorly understood. Defining substrate-binding sites of sHSPs is key to understanding their cellular functions and to harnessing their aggregation-prevention properties for controlling damage due to stress and disease. We incorporated a photoactivatable cross-linker at 32 positions throughout a well-characterized sHSP, dodecameric PsHsp18.1 from pea, and identified direct interaction sites between sHSPs and substrates. Model substrates firefly luciferase and malate dehydrogenase form strong contacts with multiple residues in the sHSP N-terminal arm, demonstrating the importance of this flexible and evolutionary variable region in substrate binding. Within the conserved alpha-crystallin domain both substrates also bind the beta-strand (beta7) where mutations in human homologs result in inherited disease. Notably, these binding sites are poorly accessible in the sHSP atomic structure, consistent with major structural rearrangements being required for substrate binding. Detectable differences in the pattern of cross-linking intensity of the two substrates and the fact that substrates make contacts throughout the sHSP indicate that there is not a discrete substrate binding surface. Our results support a model in which the intrinsically-disordered N-terminal arm can present diverse geometries of interaction sites, which is likely critical for the ability of sHSPs to protect efficiently many different substrates.