Binding Site Flexibility Research Articles

Physical signals for protein–DNA recognition

This paper discovers consensus physical signals around eukaryotic splice sites, transcription start sites, and replication origin start and end sites on a genome-wide scale based on their DNA flexibility profiles calculated by three different flexibility models. These salient physical signals are localized highly rigid and flexible DNAs, which may play important roles in protein–DNA recognition by the sliding search mechanism. The found physical signals lead us to a detailed hypothetical view of the search process in which a DNA-binding protein first finds a genomic region close to the target site from an arbitrary starting location by three-dimensional (3D) hopping and intersegment transfer mechanisms for long distances, and subsequently uses the one-dimensional (1D) sliding mechanism facilitated by the localized highly rigid DNAs to accurately locate the target flexible binding site within 30 bp (base pair) short distances. Guided by these physical signals, DNA-binding proteins rapidly search the entire genome to recognize a specific target site from the 3D to 1D pathway. Our findings also show that current promoter prediction programs (PPPs) based on DNA physical properties may suffer from lots of false positives because other functional sites such as splice sites and replication origins have similar physical signals as promoters do.

Flexible Docking in PDZ domains using Elastic Network Model and Replica Exchange Molecular Dynamics

Molecular docking is a key tool for studying protein-ligand or protein-protein interactions and for designing new drugs. Majority of the current binding/docking methods attempt to predict the bound ligand by keeping the protein (receptor) fixed and moving the target ligand around the binding site while performing an energy minimization. Modeling of protein binding site flexibility is still a challenging problem due to the large conformational space that must be sampled and inaccurate energy function.Here, we discuss a new strategy for achieving successful flexible docking of peptides to PDZ domains. Due to promiscuous behavior of PDZ domain proteins, two different sets of conformations are obtained by perturbing the unbound structure along the normal modes of elastic network model (ENM) responsible for Class I and Class II type binding. A restrained replica exchange molecular dynamics (REMD) is applied to these perturbed structures to explore the conformational space of the protein receptor. After restrained-REMD, different peptides are docked to each individual snapshots of the receptor to generate a collection of docked complexes of different stabilities.

An FTIR Investigation of Flanking Sequence Effects on the Structure and Flexibility of DNA Binding Sites

Fourier transform infrared (FTIR) spectroscopy and a library of FTIR marker bands have been used to examine the structure and relative flexibilities conferred by different flanking sequences on the EcoRI binding site. This approach allowed us to examine unique peaks and subtle changes in the spectra of d(AAAGAATTCTTT)(2), d(TTCGAATTCGAA)(2), and d(CGCGAATTCGCG)(2) and thereby identify local changes in base pairing, base stacking, backbone conformation, glycosidic bond rotation, and sugar puckering in the studied sequences. The changes in flanking sequences induce differences in the sugar puckers, glycosidic bond rotation, and backbone conformations. Varying levels of local flexibility are observed within the sequences in agreement with previous biological activity assays. The results also provide supporting evidence for the presence of a splay in the G(4)-C(9) base pair of the EcoRI binding site and a potential pocket of flexibility at the G(4) cleavage site that have been proposed in the literature. In sum, we have demonstrated that FTIR is a powerful methodology for studying the effect of flanking sequences on DNA structure and flexibility, for it can provide information about the local structure of the nucleic acid and the overall relative flexibilities conferred by different flanking sequences.

Protein structure based rational design of ecdysone agonists

We review the impact of protein X-ray crystallography on the rational design of insecticides that act as agonists of the ligand-binding domain of the Ecdysone receptor (EcR). As the EcR is a target specific to insects, these compounds potentially constitute new chemical classes of safe insecticides. The increased insight relative to that from ligand-only based (Quantitative) Structure–Activity Relations (QSARs), classical 2D-Hansch type or 3D-CoMFA/CoMSIA (Comparative Molecular Field/Similarity Analysis), is discussed. The importance of protein X-ray structure determination in support of the discovery process is stressed as the simplistic lock-and-key picture fails due to the remarkable flexibility of the EcR ligand binding site. Several new non-steroidal chemical classes of ecdysone agonists, designed by guidance from protein X-ray studies, are described.

Fully Automated Molecular Mechanics Based Induced Fit Protein−Ligand Docking Method

We describe a method for docking a ligand into a protein receptor while allowing flexibility of the protein binding site. The method employs a multistep procedure that begins with the generation of protein and ligand conformations. An initial placement of the ligand is then performed by computing binding site hotspots. This initial placement is followed by a protein side-chain refinement stage that models protein flexibility. The final step of the process is an energy minimization of the ligand pose in the presence of the rigid receptor. Thus the algorithm models flexibility of the protein at two stages, before and after ligand placement. We validated this method by performing docking and cross docking studies of eight protein systems for which crystal structures were available for at least two bound ligands. The resulting rmsd values of the 21 docked protein-ligand complexes showed values of 2 A or less for all but one of the systems examined. The method has two critical benefits for high throughput virtual screening studies. First, no user intervention is required in the docking once the initial binding site selection has been made in the protein. Second, the initial protein conformation generation needs to be performed only once for a given binding region. Also, the method may be customized in various ways depending on the particular scenario in which dockings are being performed. Each of the individual steps of the method is fully independent making it straightforward to explore different variants of the high level workflow to further improve accuracy and performance.

Active Site Pressurization: A New Tool for Structure-Guided Drug Design and Other Studies of Protein Flexibility

We present a new molecular dynamics methodology to assist in structure-based drug design and other studies that seek to predict protein deformability. Termed Active Site Pressurization (ASP), the new methodology simply injects a resin into the ligand binding-site of a protein during the course of a molecular dynamics simulation such that novel, energetically reasonable protein conformations are generated in an unbiased way that may be better representations of the ligand binding conformation than are currently available. Here we apply two different versions of the ASP methodology to three proteins, cytochrome P450cam, PcrA helicase, and glycogen synthase kinase 3beta (GSK3beta), and show that the method is capable of inducing significant conformational changes when compared to the X-ray crystal structures. Application of the ASP methodology therefore provides a view of binding site flexibility that is a rich source of data for inclusion in a variety of further investigations, including high-throughput virtual screening, lead hopping, revealing alternative modes of deformation, and revealing hidden exit and entrance tunnels.

Structure and functional analysis of the IGF-II/IGF2R interaction

Embryonic development and normal growth require exquisite control of insulin-like growth factors (IGFs). In mammals the extracellular region of the cation-independent mannose-6-phosphate receptor has gained an IGF-II-binding function and is termed type II IGF receptor (IGF2R). IGF2R sequesters IGF-II; imbalances occur in cancers and IGF2R is implicated in tumour suppression. We report crystal structures of IGF2R domains 11-12, 11-12-13-14 and domains 11-12-13/IGF-II complex. A distinctive juxtaposition of these domains provides the IGF-II-binding unit, with domain 11 directly interacting with IGF-II and domain 13 modulating binding site flexibility. Our complex shows that Phe19 and Leu53 of IGF-II lock into a hydrophobic pocket unique to domain 11 of mammalian IGF2Rs. Mutagenesis analyses confirm this IGF-II 'binding-hotspot', revealing that IGF-binding proteins and IGF2R have converged on the same high-affinity site.

Outliers in SAR and QSAR: 2. Is a flexible binding site a possible source of outliers?

Structure-activity relationship (SAR) and/or quantitative structure-activity relationship (QSAR) studies play an important role in a lead optimization of drug discovery research. When there is a lack of ligand-bound protein structural information, one of the assumptions in SAR and QSAR studies is that similar analogs bind to the same binding site in a similar binding mode. In such studies, outliers have often been observed, especially in QSAR. However, most of these studies have focused their attention on the development of QSAR and left outliers unattended. We searched ligand-bound X-ray crystal structures from the protein structure database to find evidences that could indicate a possible source of outliers in SAR or QSAR. Our results showed the possibility of conformational changes in a flexible binding site as one possible source of outliers.

Structure of PICK1 and other PDZ domains obtained with the help of self-binding C-terminal extensions.

PDZ domains are protein-protein interaction modules that generally bind to the C termini of their target proteins. The C-terminal four amino acids of a prospective binding partner of a PDZ domain are typically the determinants of binding specificity. In an effort to determine the structures of a number of PDZ domains we have included appropriate four residue extensions on the C termini of PDZ domain truncation mutants, designed for self-binding. Multiple truncations of each PDZ domain were generated. The four residue extensions, which represent known specificity sequences of the target PDZ domains and cover both class I and II motifs, form intermolecular contacts in the expected manner for the interactions of PDZ domains with protein C termini for both classes. We present the structures of eight unique PDZ domains crystallized using this approach and focus on four which provide information on selectivity (PICK1 and the third PDZ domain of DLG2), binding site flexibility (the third PDZ domain of MPDZ), and peptide-domain interactions (MPDZ 12th PDZ domain). Analysis of our results shows a clear improvement in the chances of obtaining PDZ domain crystals by using this approach compared to similar truncations of the PDZ domains without the C-terminal four residue extensions.

Promiscuous Binding of Ligands by β-Lactoglobulin Involves Hydrophobic Interactions and Plasticity

Bovine β-lactoglobulin (βLG) binds a variety of hydrophobic ligands, though precisely how is not clear. To understand the structural basis of this promiscuous binding, we studied the interaction of βLG with palmitic acid (PA) using heteronuclear NMR spectroscopy. The titration was monitored using tryptophan fluorescence and a HSQC spectrum confirmed a 1:1 stoichiometry for the PA-βLG complex. Upon the binding of PA, signal disappearances and large changes in chemical shifts were observed for the residues located at the entrance and bottom of the cavity, respectively. This observation indicates that the lower region makes a rigid connection with PA whereas the entrance is more flexible. The result is in contrast to the binding of PA to intestinal fatty acid-binding protein, another member of the calycin superfamily, in which structural consolidation occurs upon ligand binding. On the other hand, the ability of βLG to accommodate various hydrophobic ligands resembles that of GroEL, in which a large hydrophobic cavity and flexible binding site confer the ability to bind various hydrophobic substrates. Considering these observations, it is suggested that, in addition to the presence of the hydrophobic cavity, the plasticity of the entrance region makes possible the binding of hydrophobic ligands of various shapes. Thus, in contrast to the specific binding seen for many enzymes, βLG provides an example of binding with low specificity but high affinity, which may play an important role in protein–ligand and protein–protein networks.

Synergistic Use of Compound Properties and Docking Scores in Neural Network Modeling of CYP2D6 Binding: Predicting Affinity and Conformational Sampling

Cytochrome P450 2D6 (CYP2D6) is used to develop an approach for predicting affinity and relevant binding conformation(s) for highly flexible binding sites. The approach combines the use of docking scores and compound properties as attributes in building a neural network (NN) model. It begins by identifying segments of CYP2D6 that are important for binding specificity, based on structural variability among diverse CYP enzymes. A family of distinct, low-energy conformations of CYP2D6 are generated using simulated annealing (SA) and a collection of 82 compounds with known CYP2D6 affinities are docked. Interestingly, docking poses are observed on the backside of the heme as well as in the known active site. Docking scores for the active site binders, along with compound-specific attributes, are used to train a neural network model to properly bin compounds as strong binders, moderate binders, or nonbinders. Attribute selection is used to preselect the most important scores and compound-specific attributes for the model. A prediction accuracy of 85+/-6% is achieved. Dominant attributes include docking scores for three of the 20 conformations in the ensemble as well as the compound's formal charge, number of aromatic rings, and AlogP. Although compound properties were highly predictive attributes (12% improvement over baseline) in the NN-based prediction of CYP2D6 binders, their combined use with docking score attributes is synergistic (net increase of 23% above baseline). Beyond prediction of affinity, attribute selection provides a way to identify the most relevant protein conformation(s), in terms of binding competence. In the case of CYP2D6, three out of the ensemble of 20 SA-generated structures are found to be the most predictive for binding.

How Different are Structurally Flexible and Rigid Binding Sites? Sequence and Structural Features Discriminating Proteins that Do and Do not Undergo Conformational Change upon Ligand Binding

Proteins are dynamic molecules and often undergo conformational change upon ligand binding. It is widely accepted that flexible loop regions have a critical functional role in enzymes. Lack of consideration of binding site flexibility has led to failures in predicting protein functions and in successfully docking ligands with protein receptors. Here we address the question: which sequence and structural features distinguish the structurally flexible and rigid binding sites? We analyze high-resolution crystal structures of ligand bound (holo) and free (apo) forms of 41 proteins where no conformational change takes place upon ligand binding, 35 examples with moderate conformational change, and 22 cases where a large conformational change has been observed. We find that the number of residue–residue contacts observed per-residue (contact density) does not distinguish flexible and rigid binding sites, suggesting a role for specific interactions and amino acids in modulating the conformational changes. Examination of hydrogen bonding and hydrophobic interactions reveals that cases that do not undergo conformational change have high polar interactions constituting the binding pockets. Intriguingly, the large, aromatic amino acid tryptophan has a high propensity to occur at the binding sites of examples where a large conformational change has been noted. Further, in large conformational change examples, hydrophobic–hydrophobic, aromatic–aromatic, and hydrophobic–polar residue pair interactions are dominant. Further analysis of the Ramachandran dihedral angles (φ, ψ) reveals that the residues adopting disallowed conformations are found in both rigid and flexible cases. More importantly, the binding site residues adopting disallowed conformations clustered narrowly into two specific regions of the l-Ala Ramachandran map. Examination of the dihedral angles changes upon ligand binding shows that the magnitude of φ, ψ changes are in general minimal, although some large changes particularly between right-handed α-helical and extended conformations are seen. Our work further provides an account of conformational changes in the dihedral angles space. The findings reported here are expected to assist in providing a framework for predicting protein–ligand complexes and for template-based prediction of protein function.

Twenty-five Years of Delila and Molecular Information Theory

A brief personal history is given about how information theory can be applied to binding sites of genetic control molecules on nucleic acids. The primary example used is ribosome binding sites in Escherichia coli. Once the sites are aligned, the information needed to describe the sites can be computed using Claude Shannon's method. This is displayed by a computer graphic called a sequence logo. The logo represents an average binding site, and the mathematics easily allows one to determine the components of this average. That is, given a set of binding sites, the information for individual binding sites can also be computed. One can go further and predict the information of sites that are not in the original data set. Information theory also allows one to model the flexibility of ribosome binding sites, and this led us to a simple model for ribosome translational initiation in which the molecular components fit together only when the ribosome is at a good ribosome binding site. Since information theory is general, the same mathematics applies to human splice junctions, where we can predict the effect of sequence changes that cause human genetic diseases and cancer. The second example given is the Pribnow 'box' which, when viewed by the information theory method, reveals a mechanism for initiation of both transcription and DNA replication. Replication, transcription, splicing, and translation into protein represent the central dogma, so these examples show how molecular information theory is contributing to our knowledge of basic biology.

Ion selectivity in potassium channels

Potassium channels are tetrameric membrane-spanning proteins that provide a selective pore for the conduction of K + across the cell membranes. One of the main physiological functions of potassium channels is efficient and very selective transport of K + ions through the membrane to the cell. Classical views of ion selectivity are summarized within a historical perspective, and contrasted with the molecular dynamics (MD) simulations free energy perturbation (FEP) performed on the basis of the crystallographic structure of the KcsA phospholipid membrane. The results show that the KcsA channel does not select for K + ions by providing a binding site of an appropriate (fixed) cavity size. Rather, selectivity for K + arises directly from the intrinsic local physical properties of the ligands coordinating the cation in the binding site, and is a robust feature of a pore symmetrically lined by backbone carbonyl groups. Further analysis reveals that it is the interplay between the attractive ion–ligand (favoring smaller cation) and repulsive ligand–ligand interactions (favoring larger cations) that is the basic element governing Na +/K + selectivity in flexible protein binding sites. Because the number and the type of ligands coordinating an ion directly modulate such local interactions, this provides a potent molecular mechanism to achieve and maintain a high selectivity in protein binding sites despite a significant conformational flexibility.

Crystal Structure of Mouse Succinic Semialdehyde Reductase AKR7A5: Structural Basis for Substrate Specificity,

The aldo-keto reductases make up a superfamily of enzymes which can reduce a variety of aldehydes and ketones to their corresponding alcohols. Within each family are distinct preferences for certain substrates, presumably reflecting their role within the cell. The original member of the AKR7A subfamily was purified from liver as an aflatoxin dialdehyde reductase AKR7A1. However, recent additions to the family have revealed that even closely related enzymes have clear substrate preferences with AKR7A2, AKR7A4, and AKR7A5 showing much higher affinities for succinic semialdehyde (SSA) than does AKR7A1. To investigate the structural basis of this specificity, the crystal structure of mouse AKR7A5 has been determined to better than 2.5 A resolution. The structure is of the ternary complex of the enzyme with NADP+ and tartrate as an inhibitor. This structure has the same overall fold as the previously determined structure of AKR7A1; however, there are a number of differences in loops around the active site that contribute to observed differences in the substrate specificity between the AKR7A enzymes. Several differences are the result of bulky hydrophobic residues found in AKR7A5, namely, Met44, Trp77, and Trp224, which significantly restrict the size and modify the architecture of the substrate-binding pocket, producing a tighter or less flexible binding site for SSA than in AKR7A1. Site-directed mutagenesis was used to introduce Met44, Trp77, and Trp224 individually into AKR7A1, to test if they improved the affinity of the enzyme for SSA. Each mutation showed improved affinity for SSA, with Trp77Met having the largest effect. This confirms the role of these amino acids as substrate determinants for SSA.

Binding Site Flexibility: Molecular Simulation of Partial and Full Agonists within a Glutamate Receptor

Ionotropic glutamate receptors mediate fast synaptic transmission in the mammalian central nervous system and play an important role in many different functions, including memory and learning. They have also been implicated in a variety of neuropathologies and as such have generated widespread interest in their structure and function. Molecular Dynamics simulations (5 x 20 ns) of the ligand-binding core of the GluR2 glutamate receptor were performed. Through simulations of both wild type and the L650T mutant, we show that the degree of protein flexibility can be correlated with the extent to which the binding cleft is open. In agreement with recent experiments, the simulations of kainate with the wild-type construct show a slight increase in beta-sheet content that we are able to localize to two specific regions. During one simulation, the protein made a transition from an open-cleft conformation to a closed-cleft conformation. This closed cleft conformation closely resembles the closed-cleft crystal structure, thus indicating a potential pathway for conformational change associated with receptor activation. Analysis of the binding pocket suggests that partial agonists possess a greater degree of flexibility within the pocket that may help to explain why they are less efficient at opening the channel than full agonists. Examination of water molecules surrounding the ligands reveals that mobility in distinct subsites can be a discriminator between full and partial agonism and will be an important consideration in the design of drugs against these receptors.

Receptor Flexibility in de Novo Ligand Design and Docking

One of the major problems in computational drug design is incorporation of the intrinsic flexibility of protein binding sites. This is particularly crucial in ligand binding events, when induced fit can lead to protein structure rearrangements. As a consequence of the huge conformational space available to protein structures, receptor flexibility is rarely considered in ligand design procedures. In this work, we present an algorithm for integrating protein binding-site flexibility into de novo ligand design and docking processes. The approach allows dynamic rearrangement of amino acid side chains during the docking and design simulations. The impact of protein conformational flexibility is investigated in the docking of highly active inhibitors in the binding sites of acetylcholinesterase and human collagenase (matrix metalloproteinase-1) and in the design of ligands in the S1' pocket of MMP-1. The results of corresponding simulations for both rigid and flexible binding sites are compared in order to gauge the influence of receptor flexibility in drug discovery protocols.

Genetic and fluorescence studies of affinity maturation in related antibodies

Affinity maturation, the process by which an organism's response to infection becomes more specific and more effective over time, occurs after somatic hypermutation of antibody genes in B-cells. This increase in affinity might be a result of the evolution of either specific interactions between antigen and antibody over time (enthalpic factors) or antibody binding site rigidification (entropic factors) or both. Here, monoclonal antibodies, derived from antibodies elicited at different points in the murine immune response after inoculation with the same diketone hapten, have been characterized both genetically and functionally. Though this hapten has previously been shown to produce the catalytic aldolase antibody 38C2, antibodies described here are not catalytic and unlike 38C2, form no covalent enzyme–substrate complex. Thus, they provide a system in which to assess contributions to the evolution of binding affinity. The genes for these non-catalytic antibodies have been sequenced and analyzed both with regard to their relationships to germ line genes, to each other, and to two commercially available catalytic aldolase antibodies. Consequences of particular mutations for antigen binding behavior are discussed. The protein products of these genes have been expressed, purified, and binding properties measured by two complementary techniques: the hapten-induced quenching of the native antibody fluorescence and the changes in the anisotropy of Prodan (6-propionyl-2-(dimethylamino)naphthalene), a fluorescent hapten analogue. Differences in binding affinity are related back to differences in the lengths and amino acid sequences of the complementary determining region 3 (CDR3) binding loop. Taken together with our earlier results on binding site heterogeneity from tryptophan lifetime analysis [Mohan, G.S., Chiu, P.T., Southern, C.A., O’Hara, P.B., 2004. Steady-state and multifrequency phase fluorometry studies of binding site flexibility in related antibodies. J. Phys. Chem. A 108, 7871–7877], affinity appears to be modulated by a combination of entropic and enthalpic factors, and not dominated by one or the other. Because these antibodies are not related to the same germ line gene, however, these results do not provide evidence for the dominance of enthalpy or entropy in evolving binding affinity in this system.

Ligand Docking and Design in a Flexible Receptor Site

Intrinsic protein flexibility is often ignored in the drug discovery process, however, there is increasing evidence suggesting that many therapeutic targets are flexible. In this work, we describe a method for incorporating receptor flexibility, based on protein side-chain rearrangements, in ligand docking and design. The approach is applied to the docking of a highly potent inhibitor in the acetylcholinesterase binding site and to the generation of ligands in the Sl' cavity of human collagenase. Simulations are conducted for both static and flexible binding sites and from the results we assess the impact of receptor flexibility in drug discovery procedures.

RAW264.7 cells lack prostaglandin-dependent autoregulation of tumor necrosis factor-α secretion

Studies of the response of RAW264.7 cells (RAW) to lipopolysaccharide (LPS) were carried out to determine why these cells do not demonstrate the prostaglandin (PG)-dependent autocrine regulation of tumor necrosis factor-alpha (TNF-alpha) secretion observed in primary resident peritoneal macrophages (RPMs). The major cyclooxygenase (COX) product of LPS-stimulated RAW was PGD2, with lesser amounts of PGE2. LPS-treated RAW produced PGs more slowly and reached their maximal PG synthetic rate later than did LPS-treated RPMs, as a result of lower constitutive COX-1 expression and a slower rate of COX-2 induction. Cytosolic phospholipase A2 and levels of free arachidonic acid were similar in RAW and RPMs. In contrast to RPMs, LPS-treated RAW produced high quantities of TNF-alpha, which were not altered in the presence of COX inhibitors. This failure of endogenous PGs to suppress TNF-alpha secretion was explained by the absence of the prostaglandin D2 receptor and the low levels of PGE2 produced during the first 2 h of the LPS response. These studies demonstrate that autocrine regulation of TNF-alpha secretion in response to LPS is greatly facilitated by a COX-1-mediated rapid accumulation of PGs as well by a correspondence between the PGs produced and the receptors expressed by the cells.