NMR-based analysis of protein–ligand interactions

Physiological processes are mainly controlled by intermolecular recognition mechanisms involving protein–protein and protein–ligand (low molecular weight molecules) interactions. One of the most important tools for probing these interactions is high-field solution nuclear magnetic resonance (NMR) through protein-observed and ligand-observed experiments, where the protein receptor or the organic compounds are selectively detected. NMR binding experiments rely on comparison of NMR parameters of the free and bound states of the molecules. Ligand-observed methods are not limited by the protein molecular size and therefore have great applicability for analysing protein–ligand interactions. The use of these NMR techniques has considerably expanded in recent years, both in chemical biology and in drug discovery. We review here three major ligand-observed NMR methods that depend on the nuclear Overhauser effect—transferred nuclear Overhauser effect spectroscopy, saturation transfer difference spectroscopy and water–ligand interactions observed via gradient spectroscopy experiments—with the aim of reporting recent developments and applications for the characterization of protein–ligand complexes, including affinity measurements and structural determination.     

An improved adaptive genetic algorithm for protein–ligand docking

A new optimization model of molecular docking is proposed, and a fast flexible docking method based on an improved adaptive genetic algorithm is developed in this paper. The algorithm takes some advanced techniques, such as multi-population genetic strategy, entropy-based searching technique with self-adaptation and the quasi-exact penalty. A new iteration scheme in conjunction with above techniques is employed to speed up the optimization process and to ensure very rapid and steady convergence. The docking accuracy and efficiency of the method are evaluated by docking results from GOLD test data set, which contains 134 protein-ligand complexes. In over 66.2% of the complexes, the docked pose was within 2.0 A root-mean-square deviation (RMSD) of the X-ray structure. Docking time is approximately in proportion to the number of the rotatable bonds of ligands.     

Ranking docking poses by graph matching of protein–ligand interactions: lessons learned from the D3R Grand Challenge 2

A novel docking challenge has been set by the Drug Design Data Resource (D3R) in order to predict the pose and affinity ranking of a set of Farnesoid X receptor (FXR) agonists, prior to the public release of their bound X-ray structures and potencies. In a first phase, 36 agonists were docked to 26 Protein Data Bank (PDB) structures of the FXR receptor, and next rescored using the in-house developed GRIM method. GRIM aligns protein–ligand interaction patterns of docked poses to those of available PDB templates for the target protein, and rescore poses by a graph matching method. In agreement with results obtained during the previous 2015 docking challenge, we clearly show that GRIM rescoring improves the overall quality of top-ranked poses by prioritizing interaction patterns already visited in the PDB. Importantly, this challenge enables us to refine the applicability domain of the method by better defining the conditions of its success. We notably show that rescoring apolar ligands in hydrophobic pockets leads to frequent GRIM failures. In the second phase, 102 FXR agonists were ranked by decreasing affinity according to the Gibbs free energy of the corresponding GRIM-selected poses, computed by the HYDE scoring function. Interestingly, this fast and simple rescoring scheme provided the third most accurate ranking method among 57 contributions. Although the obtained ranking is still unsuitable for hit to lead optimization, the GRIM–HYDE scoring scheme is accurate and fast enough to post-process virtual screening data.     

Studying protein–protein affinity and immobilized ligand–protein affinity interactions using MS-based methods

This review discusses the most important current methods employing mass spectrometry (MS) analysis for the study of protein affinity interactions. The methods are discussed in depth with particular reference to MS-based approaches for analyzing protein–protein and protein–immobilized ligand interactions, analyzed either directly or indirectly. First, we introduce MS methods for the study of intact protein complexes in the gas phase. Next, pull-down methods for affinity-based analysis of protein–protein and protein–immobilized ligand interactions are discussed. Presently, this field of research is often called interactomics or interaction proteomics. A slightly different approach that will be discussed, chemical proteomics, allows one to analyze selectivity profiles of ligands for multiple drug targets and off-targets. Additionally, of particular interest is the use of surface plasmon resonance technologies coupled with MS for the study of protein interactions. The review addresses the principle of each of the methods with a focus on recent developments and the applicability to lead compound generation in drug discovery as well as the elucidation of protein interactions involved in cellular processes. The review focuses on the analysis of bioaffinity interactions of proteins with other proteins and with ligands, where the proteins are considered as the bioactives analyzed by MS.     

Inhibition of TRAF6-Ubc13 interaction in NFkB inflammatory pathway by analyzing the hotspot amino acid residues and protein–protein interactions using molecular docking simulations

Protein–protein interactions (PPIs) are important in most of the biochemical processes. Hotspot amino acid residues in proteins are the most important contributors for proper protein–protein interactions. Hotspot amino acid residues have been looked down upon as important therapeutic targets in inhibiting PPIs. Interaction between TRAF6 and Ubc13 is a crucial point in the NFkB inflammatory pathway. Dysfunction of the NFkB pathway is associated with numerous human diseases including cancer and neurodenegeration disorders. Ubc13 also interacts specifically to TRAF6 and not with other proteins of the TRAF family and this makes the TRAF6-Ubc13 complex an important target for specific inhibition. Hence, interfering with the TRAF6-Ubc13 association may prove effective in suppressing the NFkB disease pathway. In the present study, we searched the TRAF6-Ubc13 interaction interface to analyze their binding hotspot amino acid residues using various computational techniques. Heterocyclic compounds are known for their medicinal properties. We screened for heterocyclic analogues to the known TRAF6 inhibitor PDTC, to predict a better inhibitor using in silico protein–ligand and protein–protein interaction studies. Our in silico prediction results suggest that tetrahydro-2-thiophenecarbothioamide (Chemspider ID 36027528) binds one of the major hot-spot residues of TRAF6-Ubc13 interface and can be a better alternative in suppressing TRA6-Ubc13 complex formation in chronic inflammation than PDTC.     

Temporal resolution of protein–protein interactions in the live-cell plasma membrane

We have recently devised a method to quantify interactions between a membrane protein (“bait”) and a fluorophore-labeled protein (“prey”) directly in the live-cell plasma membrane (Schwarzenbacher et al. Nature Methods 5:1053–1060 2008). The idea is to seed cells on surfaces containing micro-patterned antibodies against the exoplasmic domain of the bait, and monitor the co-patterning of the fluorescent prey via fluorescence microscopy. Here, we characterized the time course of bait and prey micropattern formation upon seeding the cells onto the micro-biochip. Patterns were formed immediately after contact of the cells with the surface. Cells were able to migrate over the chip surface without affecting the micropattern contrast, which remained constant over hours. On single cells, bait contrast may be subject to fluctuations, indicating that the bait can be released from and recaptured on the micropatterns. We conclude that interaction studies can be performed at any time-point ranging from 5 min to several hours post seeding. Monitoring interactions with time opens up the possibility for new assays, which are briefly sketched in the discussion section.     

GalaxyDock BP2 score: a hybrid scoring function for accurate protein–ligand docking

Protein–ligand docking is a useful tool for providing atomic-level understanding of protein functions in nature and design principles for artificial ligands or proteins with desired properties. The ability to identify the true binding pose of a ligand to a target protein among numerous possible candidate poses is an essential requirement for successful protein–ligand docking. Many previously developed docking scoring functions were trained to reproduce experimental binding affinities and were also used for scoring binding poses. However, in this study, we developed a new docking scoring function, called GalaxyDock BP2 Score, by directly training the scoring power of binding poses. This function is a hybrid of physics-based, empirical, and knowledge-based score terms that are balanced to strengthen the advantages of each component. The performance of the new scoring function exhibits significant improvement over existing scoring functions in decoy pose discrimination tests. In addition, when the score is used with the GalaxyDock2 protein–ligand docking program, it outperformed other state-of-the-art docking programs in docking tests on the Astex diverse set, the Cross2009 benchmark set, and the Astex non-native set. GalaxyDock BP2 Score and GalaxyDock2 with this score are freely available at http://galaxy.seoklab.org/softwares/galaxydock.html.     

Analysis of Fas–ligand interactions using a molecular model of the receptor–ligand interface

A molecular model of the complex between Fas and its ligand was generated to better understand the location and putative effects of site-specific mutations, analyze interactions at the Fas–FasL interface, and identify contact residues. The modeling study was conservative in the sense that regions in Fas and its ligand which could not be predicted with confidence were omitted from the model to ensure accuracy of the analysis. Using the model, it was possible to map four of five N-linked glycosylation sites in Fas and FasL and to study 10 of 11 residues previously identified by mutagenesis as important for binding. Interactions involving six of these residues could be analyzed in detail and their importance for binding was rationalized based on the model. The predicted structure of the Fas–FasL interface was consistent with the experimentally established importance of these residues for binding. In addition, five previously not targeted residues were identified and predicted to contribute to binding via electrostatic interactions. Despite its limitations, the study provided a much improved basis to understand the role of Fas and FasL residues for binding compared to previous residue mapping studies using only a molecular model of Fas.     

Insights into the binding mechanism of 2,5-substituted 4-pyrone derivatives as therapeutic agents for fused dimeric interactions: A computational study using QTAIM,dynamics and docking simulations of protein–ligand complexes

The study is focused on examining 2,5-Substituted 4-Pyrone based compounds through quantum chemical and topological analysis techniques, evaluating the properties of these compounds, including their geometrical structure, intermolecular interactions and assess their possible applications. Additionally, the molecular stability, charge delocalization and UV-Visible data was investigated and compared with the calculated energy and oscillator strength using the TD-DFT approach. The researchers observed that charge transfer occurred within the molecule, indicated by the HOMO and LUMO energies. It was also found that the compound exhibited planarity and higher chemical reactivity. The calculated Mulliken charges and molecular electrostatic potential were used to interpret the Fukui index data that help predict reactive sites and understand the reactivity patterns of specific atoms in a compound. The study is aimed to understand the role of NCI in the molecule under investigation using electron localization functions and localized orbit locator methods. Molecular docking and ADMET studies were conducting involving a detailed MD simulation of a protein-ligand complex using the OPLS3e force field and the SPC water model. These findings could prove to be beneficial in developing new therapeutic agents with various pharmacological effects and potential toxicities.     

Protein–ligand interfaces are polarized: discovery of a strong trend for intermolecular hydrogen bonds to favor donors on the protein side with implications for predicting and designing ligand complexes

Understanding how proteins encode ligand specificity is fascinating and similar in importance to deciphering the genetic code. For protein–ligand recognition, the combination of an almost infinite variety of interfacial shapes and patterns of chemical groups makes the problem especially challenging. Here we analyze data across non-homologous proteins in complex with small biological ligands to address observations made in our inhibitor discovery projects: that proteins favor donating H-bonds to ligands and avoid using groups with both H-bond donor and acceptor capacity. The resulting clear and significant chemical group matching preferences elucidate the code for protein-native ligand binding, similar to the dominant patterns found in nucleic acid base-pairing. On average, 90% of the keto and carboxylate oxygens occurring in the biological ligands formed direct H-bonds to the protein. A two-fold preference was found for protein atoms to act as H-bond donors and ligand atoms to act as acceptors, and 76% of all intermolecular H-bonds involved an amine donor. Together, the tight chemical and geometric constraints associated with satisfying donor groups generate a hydrogen-bonding lock that can be matched only by ligands bearing the right acceptor-rich key. Measuring an index of H-bond preference based on the observed chemical trends proved sufficient to predict other protein–ligand complexes and can be used to guide molecular design. The resulting Hbind and Protein Recognition Index software packages are being made available for rigorously defining intermolecular H-bonds and measuring the extent to which H-bonding patterns in a given complex match the preference key.     

Recent developments in protein–ligand affinity mass spectrometry

This review provides an overview of direct and indirect technologies to screen protein–ligand interactions with mass spectrometry. These technologies have as a key feature the selection or affinity purification of ligands in mixtures prior to detection. Specific fields of interest for these technologies are metabolic profiling of bioactive metabolites, natural extract screening, and the screening of libraries for bioactives, such as parallel synthesis libraries and small combichem libraries. The review addresses the principles of each of the methods discussed, with a focus on developments in recent years, and the applicability of the methods to lead generation and development in drug discovery.     

Label-free bioelectronic detection of aptamer–protein interactions

We demonstrate for the first time the utility of nucleic acid aptamers for electrochemical detection of proteins. Highly specific and sensitive label-free detection of the target protein is achieved by combining aptamer-coated magnetic beads and chronopotentiometric stripping measurements of the captured protein (in connection to the intrinsic electroactivity of the protein). Lysozyme has thus been detected selectively in a mixture containing a large excess of six proteins and amino acids (both electroactive and non-electroactive), with a detection limit of 350 fmol (7 nM). While aptamer-based electronic sensors are in their infancy, such devices offer attractive opportunities for electrochemical detection of proteins and for developing proteomic chips.     

Receptor–ligand interactions in a reconstituted bilayer lipid membrane

A simple method is described to reconstitute membrane receptors into bilayer lipid membranes (BLMs). After reconstitution, the receptor still retains its ligand activity. Furthermore, the relationship between receptor–ligand interactions and electrical properties of reconstituted BLMs such as membrane capacitance (C_m) and membrane resistance (R_m) was studied. When glycophorin in erythrocyte and asialoglycoprotein in hepatocyte were taken as examples, it was found that the resistance of reconstituted BLM decreased when adding blood type monoclonal antibody or the solutions of galactose, respectively, and the decrease is ligand-concentration dependent; however, the membrane capacitance was not influenced. This provides a simple, practical approach to determining the interactions between the receptor and its ligand.     

Investigation of protein–protein interactions in living cells by chemical crosslinking and mass spectrometry

The identification of protein–protein interactions within their physiological environment is the key to understanding biological processes at the molecular level. However, the artificial nature of in vitro experiments, with their lack of other cellular components, may obstruct observations of specific cellular processes. In vivo analyses can provide information on the processes within a cell that might not be observed in vitro. Chemical crosslinking combined with mass spectrometric analysis of the covalently connected binding partners allows us to identify interacting proteins and to map their interface regions directly in the cell. In this paper, different in vivo crosslinking strategies for deriving information on protein–protein interactions in their physiological environment are described.     

Development of small molecules designed to modulate protein–protein interactions

Summary Protein–protein interactions are ubiquitous, essential to almost all known biological processes, and offer attractive opportunities for therapeutic intervention. Developing small molecules that modulate protein–protein interactions is challenging, owing to the large size of protein-complex interface, the lack of well-defined binding pockets, etc. We describe a general approach based on the “privileged-structure hypothesis” [Che, Ph.D. Thesis, Washington University, 2003] – that any organic templates capable of mimicking surfaces of protein-recognition motifs are potential privileged scaffolds as protein-complex antagonists – to address the challenges inherent in the discovery of small-molecule inhibitors of protein–protein interactions.This paper is adapted from a presentation at the 230th National Meeting of the American Chemical Society, Washington DC, August 28 – September 1, 2005, Abstract COMP-136.     

Electrochemistry of membrane proteins and protein–lipid assemblies

1. Download : Download high-res image (124KB)
2. Download : Download full-size image