Strategies to search and design stabilizers of protein-protein interactions: a feasibility study

Since protein-protein interactions play a pivotal role in the communication on the molecular level in virtually every biological system and process, the search and design for modulators of such interactions is of utmost importance. In recent years many inhibitors for specific protein-protein interactions have been developed, however, in only a few cases, small and druglike molecules are able to interfere in the complex formation of proteins. On the other hand, there are several small molecules known to modulate protein-protein interactions by means of stabilizing an already assembled complex. To achieve this goal, a ligand is binding to a pocket, which is located rim-exposed at the interface of the interacting proteins, for example as the phytotoxin Fusicoccin, which stabilizes the interaction of plant H+-ATPase and 14-3-3 protein by nearly a factor of 100. To suggest alternative leads, we performed a virtual screening campaign to discover new molecules putatively stabilizing this complex. Furthermore, we screen a dataset of 198 transient recognition protein-protein complexes for cavities, which are located rim-exposed at their interfaces. We provide evidence for high similarity between such rim-exposed cavities and usual ligands accommodating active sites of enzymes. This analysis suggests that rim-exposed cavities at protein-protein interfaces are druggable binding sites. Therefore, the principle of stabilizing protein-protein interactions seems to be a promising alternative to the approach of the competitive inhibition of such interactions by small molecules.     

Selectivity by Small-Molecule Inhibitors of Protein Interactions Can Be Driven by Protein Surface Fluctuations

Small-molecules that inhibit interactions between specific pairs of proteins have long represented a promising avenue for therapeutic intervention in a variety of settings. Structural studies have shown that in many cases, the inhibitor-bound protein adopts a conformation that is distinct from its unbound and its protein-bound conformations. This plasticity of the protein surface presents a major challenge in predicting which members of a protein family will be inhibited by a given ligand. Here, we use biased simulations of Bcl-2-family proteins to generate ensembles of low-energy conformations that contain surface pockets suitable for small molecule binding. We find that the resulting conformational ensembles include surface pockets that mimic those observed in inhibitor-bound crystal structures. Next, we find that the ensembles generated using different members of this protein family are overlapping but distinct, and that the activity of a given compound against a particular family member (ligand selectivity) can be predicted from whether the corresponding ensemble samples a complementary surface pocket. Finally, we find that each ensemble includes certain surface pockets that are not shared by any other family member: while no inhibitors have yet been identified to take advantage of these pockets, we expect that chemical scaffolds complementing these “distinct” pockets will prove highly selective for their targets. The opportunity to achieve target selectivity within a protein family by exploiting differences in surface fluctuations represents a new paradigm that may facilitate design of family-selective small-molecule inhibitors of protein-protein interactions.     

Protein stabilization by specific binding of guanidinium to a functional arginine-binding surface on an SH3 domain

Guanidinium hydrochloride (GuHCl) at low concentrations significantly stabilizes the Fyn SH3 domain. In this work, we have demonstrated that this stabilizing effect is manifested through a dramatic (five- to sixfold) decrease in the unfolding rate of the domain with the folding rate being affected minimally. This behavior contrasts to the effect of NaCl, which stabilizes this domain by accelerating the folding rate. These data imply that the stabilizing effect of GuHCl is not predominantly ionic in nature. Through NMR studies, we have identified a specific binding site for guanidinium, and we have determined a dissociation constant of 90 mM for this interaction. The guanidinium-binding site overlaps with a functionally important arginine-binding pocket on the domain surface, and we have shown that GuHCl is a specific inhibitor of the peptide-binding activity of the domain. A different SH3 domain possessing a similar arginine-binding pocket is also thermodynamically stabilized by GuHCl. These data suggest that many proteins that normally interact with arginine-containing ligands may also be able to specifically interact with guanidinium. Thus, some caution should be used when using GuHCl as a denaturant in protein folding studies. Since arginine-mediated interactions are often important in the energetics of protein-protein interactions, our observations could be relevant for the design of small molecule inhibitors of protein-protein interactions.     

Structural propensities of kinase family proteins from a Potts model of residue co‐variation

Understanding the conformational propensities of proteins is key to solving many problems in structural biology and biophysics. The co‐variation of pairs of mutations contained in multiple sequence alignments of protein families can be used to build a Potts Hamiltonian model of the sequence patterns which accurately predicts structural contacts. This observation paves the way to develop deeper connections between evolutionary fitness landscapes of entire protein families and the corresponding free energy landscapes which determine the conformational propensities of individual proteins. Using statistical energies determined from the Potts model and an alignment of 2896 PDB structures, we predict the propensity for particular kinase family proteins to assume a “DFG‐out” conformation implicated in the susceptibility of some kinases to type‐II inhibitors, and validate the predictions by comparison with the observed structural propensities of the corresponding proteins and experimental binding affinity data. We decompose the statistical energies to investigate which interactions contribute the most to the conformational preference for particular sequences and the corresponding proteins. We find that interactions involving the activation loop and the C‐helix and HRD motif are primarily responsible for stabilizing the DFG‐in state. This work illustrates how structural free energy landscapes and fitness landscapes of proteins can be used in an integrated way, and in the context of kinase family proteins, can potentially impact therapeutic design strategies.     

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.     

Anti-apoptotic Bcl-XL protein in complex with BH3 peptides of pro-apoptotic Bak, Bad, and Bim proteins: comparative molecular dynamics simulations

The Bcl-2 family of proteins plays a central role in the regulation of mitochondrial outer-membrane permeabilization, a critical step in apoptosis. Heterodimerization between the pro- and anti-apoptotic members of Bcl-2 family is a key event in this process. Anti-apoptotic proteins have high levels of expression in many cancers and they have different affinities for different pro-apoptotic proteins. Experimentally determined structures of all members of Bcl-2 proteins have remarkably similar helical fold despite poor amino acid sequence identity. Peptides representing BH3 region of pro-apoptotic proteins have been shown to bind the hydrophobic cleft of anti-apoptotic proteins and this segment is responsible in modulating the apoptotic pathways in living cells. Understanding the molecular basis of protein-protein recognition is required to develop inhibitors specific to a particular anti-apoptotic protein. We have carried out molecular dynamics simulations on the anti-apoptotic Bcl-X(L) protein in complex with three different BH3 peptides derived from pro-apoptotic Bak, Bad and Bim proteins. Each complex structure was simulated for a period of 50 ns after 2.5 ns equilibration. Analysis of the simulation results showed that in the Bcl-X(L) protein, the helix containing the BH3 region is more flexible than other helices in all three simulations. A network of strong hydrophobic interactions exists between four of the six helices and they contribute significantly to the stability of this helix bundle protein. Analysis of Bcl-X(L)-BH3 peptide interactions reveals the role of loop residues in the protein-peptide interactions in all three simulations. Bad and Bim peptides maintain strong hydrophobic and hydrophilic interactions with the helix preceding the central hydrophobic helix. Residues from this helix interact with an Arg residue in Bad and Bim peptides. This Arg residue is next to the conserved Leu residue and is replaced by Ala in Bak. Absence of these interactions and the helix propensity are likely to be the cause for Bak peptide's weaker binding affinity with the Bcl-X(L) protein. The results of this study have implications in the design of Bcl-X(L)-specific inhibitors.     

Rho guanine dissociation inhibitors: pivotal molecules in cellular signalling.

The small G proteins of the Ras family act as bimodal relays in the transfer of intracellular signals. This is a dynamic phenomenon involving a cascade of protein-protein interactions modulated by chemical modifications, structural rearrangements and intracellular relocalisations. Most of the small G proteins could be operationally defined as proteins having two conformational states, each of which interacts with different cellular partners. These two states are determined by the nature of the bound nucleotide, GDP or GTP. This capacity to cycle between a GDP-bound conformation and a GTP-bound conformation enables them to filter, to amplify or to temporise the upstream signals that they receive. Thus the control of this cycle is crucial. Membrane anchoring of the proteins in the Ras family is a prerequisite for their activity. Most of the proteins in the Rho/Rac and Rab subfamilies of Ras proteins cycle between cytosol and membranes. Then the control of membrane association/dissociation is an other important regulation level. This review will describe one family of crucial regulators acting on proteins in the Rho/Rac family-the Rho guanine nucleotide dissociation inhibitors, or RhoGDIs. As yet, only three RhoGDIs have been described: RhoGDI-1, RhoGDI-2 (or D4/Ly-GDI) and RhoGDI-3. RhoGDI 1 and 2 are cytosolic and participate in the regulation of both the GDP/GTP cycle and the membrane association/dissociation cycle of Rho/Rac proteins. The non-cytosolic RhoGDI-3 seems to act in a slightly different way.     

Bromodomain and extra-terminal (BET) family proteins: New therapeutic targets in major diseases

The bromodomains and extra-terminal domain (BET) family proteins recognize acetylated chromatin through their bromodomains (BDs) and help in regulating gene expression. BDs are chromatin ‘readers’: by interacting with acetylated lysines on the histone tails, they recruit chromatin-regulating proteins on the promoter region to regulate gene expression and repression. Extensive efforts have been employed by scientific communities worldwide to identify and develop potential inhibitors of BET family BDs to regulate protein expression by inhibiting acetylated histone (H3/H4) interactions. Several small molecule inhibitors have been reported, which not only have high affinity but also have high specificity to BET BDs. These developments make BET family proteins an important therapeutic targets for major diseases such as cancer, neurological disorders, obesity and inflammation. Here, we review and discuss the structural biology of BET family BDs and their applications in major diseases.     

Protein inhibitors of serine proteinases

Serine proteinases and their natural protein inhibitors belong to the most intensively studied models of protein-protein recognition. Protein inhibitors do not form a single group but can be divided into about 20 different families. Global structures of proteins representing different inhibitor families are completely different and comprise alpha-helical proteins, beta-sheet proteins, alpha/beta-proteins and different folds of small disulfide-rich proteins. Three different types of inhibitors can be distinguished: canonical (standard mechanism) inhibitors, non-canonical inhibitors, and serpins. The canonical inhibitor binds to the enzyme through the exposed and convex binding loop, which is complementary to the active site of the enzyme. The mechanism of inhibition in this group is consistently very similar and resembles that of an ideal substrate. Non-canonical inhibitors, originating from blood sucking organisms, specifically block enzymes of the blood clotting cascade. The interaction is mediated through inhibitor N-terminus which binds to the proteinase forming a parallel beta-sheet. There are also extensive secondary interactions which provide an additional buried area and contribute significantly to the strength and specificity of recognition. Serpins are major proteinase inhibitors occurring in plasma. Similarly to canonical inhibitors, serpins interact with their target proteinases in a substrate-like manner. However, in the case of serpins, cleavage of a single peptide bond in a flexible and exposed binding loop leads to dramatic structural changes.     

Structural Basis for Influence of Viral Glycans on Ligand Binding by Influenza Hemagglutinin

Binding of cell surface glycans by influenza hemagglutinin controls viral attachment and infection of host cells. This binding is a three-way interaction between viral proteins, host glycans, and viral glycans; many structural details of this interaction have been difficult to resolve. Here, we use a series of 100-ns molecular dynamics simulations to further analyze available crystallographic data on hemagglutinin-ligand interactions. Based on our simulations, we predict that the viral glycans contact the host glycans within 1-2 residues of the ligand-binding site. We also predict that the glycan-glycan interactions contain both stabilizing and destabilizing components. These predictions suggest a structural means to explain why changes to viral glycosylation alter the efficiency and selectivity of ligand binding. We also predict that the proximity of these interactions to the ligand-binding pocket will impact the binding affinity of small glycomimetic ligands analogous to the influenza neuraminidase inhibitors currently in clinical use.     

Exploring the role of cation-pi interactions in glycoproteins lipid-binding proteins and RNA-binding proteins

We have analyzed and compared the influence of cation-pi interactions in glycoproteins (GPs), lipid-binding proteins (LBPs) and RNA-binding proteins (RBPs) in this study. We observed that all the proteins included in the study had profound cation-pi interactions. There is an average of one energetically significant cation-pi interaction for every 71 residues in GPs, for every 58 residues in LBPs and for every 64 residues in RBPs. Long-range contacts are predominant in all the three types of proteins studied. The pair-wise cation-pi interaction energy between the positively charged and aromatic residues shows that Arg-Trp pair energy was the strongest among all six possible pairs in all the three types of proteins studied. There were considerable differences in the preference of cation-pi interacting residues to different secondary structure elements and ASA and these might contribute to differences in biochemical functions of GPs, LBPs and RBPs. It was interesting to note that all the five residues involved in cation-pi interactions were found to have stabilization centers in GPs, LBPs and RBPs. Majority of the cation-pi interacting residues investigated in the present study had a conservation score of 6, the cutoff value used to identify the stabilizing residues. A small percentage of cation-pi interacting residues were also present as stabilizing residues. The cation-pi interaction-forming residues play an important role in the structural stability of in GPs, LBPs and RBPs. The results obtained in this study will be helpful in further understanding the stability, specificity and differences in the biochemical functions of GPs, LBPs and RBPs.     

Optimal local propensities for model proteins

Lattice models of proteins were used to examine the role of local propensities in stabilizing the native state of a protein, using techniques drawn from spin-glass theory to characterize the free-energy landscapes. In the strong evolutionary limit, optimal conditions for folding are achieved when the contributions from local interactions to the stability of the native state is small. Further increasing the local interactions rapidly decreases the foldability. © 1995 Wiley-Liss, Inc.     

Prokaryote-derived protein inhibitors of peptidases: A sketchy occurrence and mostly unknown function

In metazoan organisms protein inhibitors of peptidases are important factors essential for regulation of proteolytic activity. In vertebrates genes encoding peptidase inhibitors constitute up to 1% of genes reflecting a need for tight and specific control of proteolysis especially in extracellular body fluids. In stark contrast unicellular organisms, both prokaryotic and eukaryotic consistently contain only few, if any, genes coding for putative peptidase inhibitors. This may seem perplexing in the light of the fact that these organisms produce large numbers of proteases of different catalytic classes with the genes constituting up to 6% of the total gene count with the average being about 3%. Apparently, however, a unicellular life-style is fully compatible with other mechanisms of regulation of proteolysis and does not require protein inhibitors to control their intracellular and extracellular proteolytic activity. So in prokaryotes occurrence of genes encoding different types of peptidase inhibitors is infrequent and often scattered among phylogenetically distinct orders or even phyla of microbiota. Genes encoding proteins homologous to alpha-2-macroglobulin (family I39), serine carboxypeptidase Y inhibitor (family I51), alpha-1-peptidase inhibitor (family I4) and ecotin (family I11) are the most frequently represented in Bacteria. Although several of these gene products were shown to possess inhibitory activity, with an exception of ecotin and staphostatins, the biological function of microbial inhibitors is unclear. In this review we present distribution of protein inhibitors from different families among prokaryotes, describe their mode of action and hypothesize on their role in microbial physiology and interactions with hosts and environment.     

Biochemical and pharmacological profiling of the pro-survival protein Bcl-xL

Bcl-2 family proteins are key mediators of programmed cell death. Over-expression of anti-apoptotic Bcl-2 family members (Bcl-2, Bcl-xL, and Mcl-1) has been associated with tumor progression and chemotherapeutic resistance. Pharmacological agents that neutralize the functions of anti-apoptotic Bcl-2 family proteins have emerged as a promising new class of anti-cancer agents. Biochemical analyses have demonstrated that small molecule inhibitors and some pro-apoptotic proteins exhibit distinct binding preferences for anti-apoptotic proteins. While numerous structures of anti-apoptotic proteins bound to ligands have been reported, the source of this selectivity is still unclear. Here, we present a systematic analysis of a series of Bcl-xL variants that contain mutations within the hydrophobic ligand-binding cleft. The ability of these Bcl-xL mutants to interact with both small molecule inhibitors and BH3 peptides was determined. These studies provide information on the contributions of specific residues to small molecule inhibitor binding and shed light on the ligand selectivity of these therapeutically important proteins.     

Protein fusions of BPI with CETP retain functions inherent to each

Cholesteryl ester transfer protein (CETP), bactericidal/permeability inducing protein (BPI), and lipopolysaccharide binding protein (LBP) are members of the lipid transfer/lipopolysaccharide binding protein (LT/LBP) family of proteins that share a common secondary/tertiary structure. Despite this commonality of structure, very different patterns of lipid binding and protein-protein interactions are observed among the family members. BPI was previously shown to retain aspects of its own function when part of it was fused with LBP to form a chimeric protein. We have extended those observations to CETP. Some aspects of cholesteryl ester transfer function can be maintained in a chimeric protein even when over 40% of the sequence is from BPI. Further replacement of an additional 60 amino acids resulted in a complete loss of CETP function even though the chimera was able to retain some BPI-like properties. These artificial fusions retain BPI functions such as lipopolysaccharide (LPS) binding and protein-protein interactions that are not observed with native CETP. BPI-CETP chimeras are inhibited by LPS but cannot be inhibited by small molecule CETP inhibitors as effectively as native CETP. These results localize the site of LPS binding in BPI to a region no larger than the amino terminal 155 amino acids. This region can participate in some protein-protein interactions similar to intact BPI. Chimeras containing the amino terminus of CETP and the carboxy terminus of BPI did not retain any observable CETP function. These results further confirm the modular nature of the LT/LBP family of proteins but also highlight the discrete nature of their individual functions.     

Detecting molecular interactions that stabilize, activate and guide ligand-binding of the sodium/proton antiporter MjNhaP1 from Methanococcus jannaschii

Integral membrane proteins are involved in virtually every cellular process. Precisely regulating these machineries would allow controlling many human and vertebrate diseases. Embedded into cellular membranes, membrane proteins establish molecular interactions that sensitively react to environmental changes and to molecular compounds, such as ligands or inhibitors. We applied atomic force microscopy (AFM) to image the Na(+)/H(+) antiporter MjNhaP1 from Methanococcus jannaschii, and single-molecule force spectroscopy (SMFS) to probe molecular interactions that drive the protein structure-function relationship. High-resolution AFM topographs showed the dimeric assembly of MjNhaP1 being reconstituted into a lipid bilayer. SMFS of MjNhaP1 unraveled molecular interactions stabilizing individual structural domains. Transmembrane domains exhibited certain probabilities to unfold individually or cooperatively with other domains resulting in different unfolding pathways. Helices VIII and X established pH sensitive interactions altering significantly upon MjNhaP1 activation, while removal of the ligand (Na(+)) destabilized the entire antiporter except helix VIII. It is assumed that Asp234/235 of helix VIII are involved in the ligand-binding site and that helix X plays a functional role in the activation of the transporter.     

Differential scanning fluorometry signatures as indicators of enzyme inhibitor mode of action: case study of glutathione S-transferase

Differential scanning fluorometry (DSF), also referred to as fluorescence thermal shift, is emerging as a convenient method to evaluate the stabilizing effect of small molecules on proteins of interest. However, its use in the mechanism of action studies has received far less attention. Herein, the ability of DSF to report on inhibitor mode of action was evaluated using glutathione S-transferase (GST) as a model enzyme that utilizes two distinct substrates and is known to be subject to a range of inhibition modes. Detailed investigation of the propensity of small molecule inhibitors to protect GST from thermal denaturation revealed that compounds with different inhibition modes displayed distinct thermal shift signatures when tested in the presence or absence of the enzyme's native co-substrate glutathione (GSH). Glutathione-competitive inhibitors produced dose-dependent thermal shift trendlines that converged at high compound concentrations. Inhibitors acting via the formation of glutathione conjugates induced a very pronounced stabilizing effect toward the protein only when GSH was present. Lastly, compounds known to act as noncompetitive inhibitors exhibited parallel concentration-dependent trends. Similar effects were observed with human GST isozymes A1-1 and M1-1. The results illustrate the potential of DSF as a tool to differentiate diverse classes of inhibitors based on simple analysis of co-substrate dependency of protein stabilization.     

Electrostatic contributions to the stability of hyperthermophilic proteins.

Electrostatic contributions to the folding free energy of several hyperthermophilic proteins and their mesophilic homologs are calculated. In all the cases studied, electrostatic interactions are more favorable in the hyperthermophilic proteins. The electrostatic free energy is found not to be correlated with the number of ionizable amino acid residues, ion pairs or ion pair networks in a protein, but rather depends on the location of these groups within the protein structure. Moreover, due to the large free energy cost associated with burying charged groups, buried ion pairs are found to be destabilizing unless they undergo favorable interactions with additional polar groups, including other ion pairs. The latter case involves the formation of stabilizing ion pair networks as is observed in a number of proteins. Ion pairs located on the protein surface also provide stabilizing interactions in a number of cases. Taken together, our results suggest that many hyperthermophilic proteins enhance electrostatic interactions through the optimum placement of charged amino acid residues within the protein structure, although different design strategies are used in different cases. Other physical mechanisms are also likely to contribute, however optimizing electrostatic interactions offers a simple means of enhancing stability without disrupting the core residues characteristic of different protein families.