Optimized small‐molecule pull‐downs define MLBP1 as an acyl‐lipid‐binding protein

Abscisic acid (ABA) receptors belong to the START domain superfamily, which encompasses ligand‐binding proteins present in all kingdoms of life. START domain proteins contain a central binding pocket that, depending on the protein, can couple ligand binding to catalytic, transport or signaling functions. In Arabidopsis, the best characterized START domain proteins are the 14 PYR/PYL/RCAR ABA receptors, while the other members of the superfamily do not have assigned ligands. To address this, we used affinity purification of biotinylated proteins expressed transiently in Nicotiana benthamiana coupled to untargeted LC‐MS to identify candidate binding ligands. We optimized this method using ABA–PYL interactions and show that ABA co‐purifies with wild‐type PYL5 but not a binding site mutant. The K_d of PYL5 for ABA is 1.1 μm , which suggests that the method has sufficient sensitivity for many ligand–protein interactions. Using this method, we surveyed a set of 37 START domain‐related proteins, which resulted in the identification of ligands that co‐purified with MLBP1 (At4G01883) or MLP165 (At1G35260). Metabolite identification and the use of authentic standards revealed that MLBP1 binds to monolinolenin, which we confirmed using recombinant MLBP1. Monolinolenin also co‐purified with MLBP1 purified from transgenic Arabidopsis, demonstrating that the interaction occurs in a native context. Thus, deployment of this relatively simple method allowed us to define a protein–metabolite interaction and better understand protein–ligand interactions in plants.     

Q-SiteFinder: an energy-based method for the prediction of protein-ligand binding sites

MOTIVATION: Identifying the location of ligand binding sites on a protein is of fundamental importance for a range of applications including molecular docking, de novo drug design and structural identification and comparison of functional sites. Here, we describe a new method of ligand binding site prediction called Q-SiteFinder. It uses the interaction energy between the protein and a simple van der Waals probe to locate energetically favourable binding sites. Energetically favourable probe sites are clustered according to their spatial proximity and clusters are then ranked according to the sum of interaction energies for sites within each cluster. RESULTS: There is at least one successful prediction in the top three predicted sites in 90% of proteins tested when using Q-SiteFinder. This success rate is higher than that of a commonly used pocket detection algorithm (Pocket-Finder) which uses geometric criteria. Additionally, Q-SiteFinder is twice as effective as Pocket-Finder in generating predicted sites that map accurately onto ligand coordinates. It also generates predicted sites with the lowest average volumes of the methods examined in this study. Unlike pocket detection, the volumes of the predicted sites appear to show relatively low dependence on protein volume and are similar in volume to the ligands they contain. Restricting the size of the pocket is important for reducing the search space required for docking and de novo drug design or site comparison. The method can be applied in structural genomics studies where protein binding sites remain uncharacterized since the 86% success rate for unbound proteins appears to be only slightly lower than that of ligand-bound proteins. AVAILABILITY: Both Q-SiteFinder and Pocket-Finder have been made available online at http://www.bioinformatics.leeds.ac.uk/qsitefinder and http://www.bioinformatics.leeds.ac.uk/pocketfinder     

Comparison analysis of primary ligand‐binding sites in seven‐helix membrane proteins

Seven‐helix transmembrane proteins, including the G‐protein‐coupled receptors (GPCRs), mediate a broad range of fundamental cellular activities through binding to a wide range of ligands. Understanding the structural basis for the ligand‐binding selectivity of these proteins is of significance to their structure‐based drug design. Comparison analysis of proteins' ligand‐binding sites provides a useful way to study their structure‐activity relationships. Various computational methods have been developed for the binding‐site comparison of soluble proteins. In this work, we applied this approach to the analysis of the primary ligand‐binding sites of 92 seven‐helix transmembrane proteins. Results of the studies confirmed that the binding site of bacterial rhodopsins is indeed different from all GPCRs. In the latter group, further comparison of the binding sites indicated a group of residues that could be responsible for ligand‐binding selectivity and important for structure‐based drug design. Furthermore, unexpected binding‐site dissimilarities were observed among adrenergic and adenosine receptors, suggesting that the percentage of the overall sequence identity between a target protein and a template protein alone is not sufficient for selecting the best template for homology modeling of seven‐helix membrane proteins. These results provided novel insight into the structural basis of ligand‐binding selectivity of seven‐helix membrane proteins and are of practical use to the computational modeling of these proteins. © 2010 Wiley Periodicals, Inc. Biopolymers 95: 31–38, 2011.     

Surface‐based protein binding pocket similarity

Protein similarity comparisons may be made on a local or global basis and may consider sequence information or differing levels of structural information. We present a local three‐dimensional method that compares protein binding site surfaces in full atomic detail. The approach is based on the morphological similarity method which has been widely applied for global comparison of small molecules. We apply the method to all‐by‐all comparisons two sets of human protein kinases, a very diverse set of ATP‐bound proteins from multiple species, and three heterogeneous benchmark protein binding site data sets. Cases of disagreement between sequence‐based similarity and binding site similarity yield informative examples. Where sequence similarity is very low, high pocket similarity can reliably identify important binding motifs. Where sequence similarity is very high, significant differences in pocket similarity are related to ligand binding specificity and similarity. Local protein binding pocket similarity provides qualitatively complementary information to other approaches, and it can yield quantitative information in support of functional annotation. Proteins 2011; © 2011 Wiley‐Liss, Inc.     

Theoretical investigation on the glycan‐binding specificity of Agrocybe cylindracea galectin using molecular modeling and molecular dynamics simulation studies

Galectins are β‐galactoside binding proteins which have the ability to serve as potent antitumor, cancer biomarker, and induce tumor cell apoptosis. Agrocybe cylindracea galectin (ACG) is a fungal galectin which specifically recognizes α(2,3)‐linked sialyllactose at the cell surface that plays extensive roles in the biological recognition processes. To investigate the change in glycan‐binding specificity upon mutations, single point and double point site‐directed in silico mutations are performed at the binding pocket of ACG. Molecular dynamics (MD) simulation studies are carried out for the wild‐type (ACG) and single point (ACG1) and double point (ACG2) mutated ACGs to investigate the dynamics of substituted mutants and their interactions with the receptor sialyllactose. Plausible binding modes are proposed for galectin–sialylglycan complexes based on the analysis of hydrogen bonding interactions, total pair‐wise interaction energy between the interacting binding site residues and sialyllactose and binding free energy of the complexes using molecular mechanics–Poisson–Boltzmann surface area. Our result shows that high contribution to the binding in different modes is due to the direct and water‐mediated hydrogen bonds. The binding specificity of double point mutant Y59R/N140Q of ACG2 is found to be high, and it has 26 direct and water‐mediated hydrogen bonds with a relatively low‐binding free energy of −47.52 ± 5.2 kcal/mol. We also observe that the substituted mutant Arg59 is crucial for glycan‐binding and for the preference of α(2,3)‐linked sialyllactose at the binding pocket of ACG2 galectin. When compared with the wild‐type and single point mutant, the double point mutant exhibits enhanced affinity towards α(2,3)‐linked sialyllactose, which can be effectively used as a model for biological cell marker in cancer therapeutics. Copyright © 2015 John Wiley & Sons, Ltd.     

Exploring the role of the phospholipid ligand in endothelial protein C receptor: A molecular dynamics study

Endothelial protein C receptor (EPCR) is a CD1‐like transmembrane glycoprotein with important regulatory roles in protein C (PC) pathway, enhancing PC's anticoagulant, anti‐inflammatory, and antiapoptotic activities. Similarly to homologous CD1d, EPCR binds a phospholipid [phosphatidylethanolamine (PTY)] in a groove corresponding to the antigen‐presenting site, although it is not clear if lipid exchange can occur in EPCR as in CD1d. The presence of PTY seems essential for PC γ‐carboxyglutamic acid (Gla) domain binding. However, the lipid‐free form of the EPCR has not been characterized. We have investigated the structural role of PTY on EPCR, by multiple molecular dynamics (MD) simulations of ligand bound and unbound forms of the protein. Structural changes, subsequent to ligand removal, led to identification of two stable and folded ligand‐free conformations. Compared with the bound form, unbound structures showed a narrowing of the A′ pocket and a high flexibility of the helices around it, in agreement with CD1d simulation. Thus, a lipid exchange with a mechanism similar to CD1d is proposed. In addition, unbound conformations presented a reduced interaction surface for Gla domain, confirming the role of PTY in establishing the proper EPCR conformation for the interaction with its partner protein. Single MD simulations were also obtained for 29 mutant models with predicted structural stability and impaired binding ability. Ligand affinity calculations, based on linear interaction energy method, showed that substitution‐induced conformational changes affecting helices around the A′ pocket were associated to a reduced binding affinity. Mutants responsible for this effect may represent useful reagents for experimental tests. Proteins 2010. © 2010 Wiley‐Liss, Inc.     

Tuning surface‐cross‐linking of molecularly imprinted cross‐linked micelles for molecular recognition in water

Molecular recognition in water is an important challenge in supramolecular chemistry. Surface‐core double cross‐linking of template‐containing surfactant micelles by the click reaction and free radical polymerization yields molecularly imprinted nanoparticles (MINPs) with guest‐complementary binding sites. An important property of MINP‐based receptors is the surface‐cross‐linking between the propargyl groups of the surfactants and a diazide cross‐linker. Decreasing the number of carbons in between the two azides enhanced the binding affinity of the MINPs, possibly by keeping the imprinted binding site more open prior to the guest binding. The depth of the binding pocket can be controlled by the distribution of the hydrophilic/hydrophobic groups of the template and was found to influence the binding in addition to electrostatic interactions between oppositely charged MINPs and guests. Cross‐linkers with an alkoxyamine group enabled two‐stage double surface‐cross‐linking that strengthened the binding constants by an order of magnitude, possibly by expanding the binding pocket of the MINP into the polar region. The binding selectivity among very similar isomeric structures also improved.     

Molecular recognition of Fc‐specific ligands binding onto the consensus binding site of IgG: insights from molecular simulation

Immunoglobulin G (IgG) plays an important role in clinical diagnosis and therapeutics. Meanwhile, the consensus binding site (CBS) on the Fc domain of IgG is responsible for ligand recognition, especially for Fc‐specific ligands. In this study, molecular simulation methods were used to investigate molecular interactions between the CBS of the Fc domain and seven natural Fc‐specific ligands. The analysis on the binding energy of the Fc–ligand complex indicated that hydrophobic interactions provide the main driving force for the Fc–ligand binding processes. The hot spots on the ligands and Fc were identified with the computational alanine scanning approach. It was found that the residues of tryptophan and tyrosine on the ligands have significant contributions for the Fc–ligand binding, while Met252, Ile253, Asn434, His435, and Tyr436 are the key residues of Fc. Moreover, two binding modes based on tryptophan or tyrosine were summarized and constructed according to the pairwise interaction analysis. Guidelines for the rational design of CBS‐specific ligands with high affinity and specificity were proposed. Copyright © 2014 John Wiley & Sons, Ltd.     

Molecular dynamics simulation of the induced‐fit binding process of DNA aptamer and L‐argininamide

Aptamers are rare functional nucleic acids with binding affinity to and specificity for target ligands. Recent experiments have lead to the proposal of an induced‐fit binding mechanism for L ‐argininamide (Arm) and its binding aptamer. However, at the molecular level, this mechanism between the aptamer and its coupled ligand is still poorly understood. The present study used explicit solvent molecular dynamics (MD) simulations to examine the critical bases involved in aptamer‐Arm binding and the induced‐fit binding process at atomic resolution. The simulation results revealed that the Watson‐Crick pair (G10‐C16), C9, A12, and C17 bases play important roles in aptamer‐Arm binding, and that binding of Arm results in an aptamer conformation optimized through a general induced‐fit process. In an aqueous solution, the mechanism has the following characteristic stages: (a) adsorption stage, the Arm anchors to the binding site of aptamer with strong electrostatic interaction; (b) binding stage, the Arm fits into the binding site of aptamer by hydrogen‐bond formation; and (c) complex stabilization stage, the hydrogen bonding and electrostatic interactions cooperatively stabilize the complex structure. This study provides dynamics information on the aptamer‐ligand induced‐fit binding mechanism. The critical bases in aptamer‐ligand binding may provide a guideline in aptamer design for molecular recognition engineering.     

Steric and thermodynamic limits of design for the incorporation of large unnatural amino acids in aminoacyl‐tRNA synthetase enzymes

Orthogonal aminoacyl‐tRNA synthetase/tRNA pairs from archaea have been evolved to facilitate site specific in vivo incorporation of unnatural amino acids into proteins in Escherichia coli. Using this approach, unnatural amino acids have been successfully incorporated with high translational efficiency and fidelity. In this study, CHARMM‐based molecular docking and free energy calculations were used to evaluate rational design of specific protein–ligand interactions for aminoacyl‐tRNA synthetases. A series of novel unnatural amino acid ligands were docked into the p‐benzoyl‐L ‐phenylalanine tRNA synthetase, which revealed that the binding pocket of the enzyme does not provide sufficient space for significantly larger ligands. Specific binding site residues were mutated to alanine to create additional space to accommodate larger target ligands, and then mutations were introduced to improve binding free energy. This approach was used to redesign binding sites for several different target ligands, which were then tested against the standard 20 amino acids to verify target specificity. Only the synthetase designed to bind Man‐α‐O‐Tyr was predicted to be sufficiently selective for the target ligand and also thermodynamically stable. Our study suggests that extensive redesign of the tRNA synthatase binding pocket for large bulky ligands may be quite thermodynamically unfavorable. Proteins 2010. © 2010 Wiley‐Liss, Inc.     

Identification of protein functional surfaces by the concept of a split pocket

The function of a protein is often fulfilled via molecular interactions on its surfaces, so identifying the functional surface(s) of a protein is helpful for understanding its function. Here, we introduce the concept of a split pocket, which is a pocket that is split by a cognate ligand. We use a geometric approach that is site‐specific. Specifically, we first compute a set of all pockets in the protein with its ligand(s) and a set of all pockets with the ligand(s) removed and then compare the two sets of pockets to identify the split pocket(s) of the protein. To reduce the search space and expedite the process of surface partitioning, we design probe radii according to the physicochemical textures of molecules. Our method achieves a success rate of 96% on a benchmark test set. We conduct a large‐scale computation to identify ～19,000 split pockets from 11,328 structures (1.16 million potential pockets); for each pocket, we obtain residue composition, solvent‐accessible area, and molecular volume. With this database of split pockets, our method can be used to predict the functional surfaces of unbound structures. Indeed, the functional surface of an unbound protein may often be found from its similarity to remotely related bound forms that belong to distinct folds. Finally, we apply our method to identify glucose‐binding proteins, including unbound structures. Our study demonstrates the power of geometric and evolutionary matching for studying protein functional evolution and provides a framework for classifying protein functions by local spatial patterns of functional surfaces. Proteins 2009. © 2009 Wiley‐Liss, Inc.     

Structural studies of the C‐terminal 19‐peptide of serum amyloid A and its Pro→Ala variants interacting with human cystatin C

Serum amyloid A (SAA) is a multifunctional acute‐phase protein whose concentration in serum increases markedly following a number of chronic inflammatory and neoplastic diseases. Prolonged high SAA level may give rise to reactive systemic amyloid A (AA) amyloidosis, where the N‐terminal segment of SAA is deposited as amyloid fibrils. Besides, recently, well‐documented association of SAA with high‐density lipoprotein or glycosaminoglycans, in particular heparin/heparin sulfate (HS), and specific interaction between SAA and human cystatin C (hCC), the ubiquitous inhibitor of cysteine proteases, was proved. Using a combination of selective proteolytic excision and high‐resolution mass spectrometry, a hCC binding site in the SAA sequence was determined as SAA(86–104). The role of this SAA C‐terminal fragment as a ligand‐binding locus is still not clear. It was postulated important in native SAA folding and in pathogenesis of AA amyloidosis. In the search of conformational details of this SAA fragment, we did its structure and affinity studies, including its selected double/triple Pro→Ala variants. Our results clearly show that the SAA(86–104) 19‐peptide has rather unordered structure with bends in its C‐terminal part, which is consistent with the previous results relating to the whole protein. The results of affinity chromatography, fluorescent ELISA‐like test, CD and NMR studies point to an importance of proline residues on structure of SAA(86–104). Conformational details of SAA fragment, responsible for hCC binding, may help to understand the objective of hCC–SAA complex formation and its importance for pathogenesis of reactive amyloid A amyloidosis. Copyright © 2015 John Wiley & Sons, Ltd.     

Profiling the interaction of 1‐phenylbenzimidazoles to cyclooxygenases

In the design of 1‐phenylbenzimidazoles as model cyclooxygenase (COX) inhibitors, docking to a series of crystallographic COX structures was performed to evaluate their potential for high‐affinity binding and to reproduce the interaction profile of well‐known COX inhibitors. The effect of ligand‐specific induced fit on the calculations was also studied. To quantitatively compare the pattern of interactions of model compounds to the profile of several cocrystallized COX inhibitors, a geometric parameter, denominated ligand‐receptor contact distance (LRCD), was developed. The interaction profile of several model complexes showed similarity to the profile of COX complexes with inhibitors such as iodosuprofen, iodoindomethacin, diclofenac, and flurbiprofen. Shaping of high‐affinity binding sites upon ligand‐specific induced fit mostly determined both the affinity and the binding mode of the ligands in the docking calculations. The results suggest potential of 1‐phenylbenzimidazole derivatives as COX inhibitors on the basis of their predicted affinity and interaction profile to COX enzymes. The analyses also provided insights into the role of induced fit in COX enzymes. While inhibitors produce different local structural changes at the COX ligand binding site, induced fit allows inhibitors in diverse chemical classes to share characteristic interaction patterns that ensure key contacts to be achieved. Different interaction patterns may also be associated with different inhibitory mechanisms.     

From signal perception to signal transduction: ligand‐induced dimeric switch of DctB sensory domain in solution

Sinorhizobium meliloti DctB is a typical transmembrane sensory histidine kinase, which senses C₄‐dicarboxylic acids (DCA) and regulates the expression of DctA, the DCA transporter. We previously reported the crystal structures of its periplasmic sensory domain (DctBp) in apo and succinate‐bound states, and these structures showed dramatic conformational changes at dimeric level. Here we show a ligand‐induced dimeric switch in solution and a strong correlation between DctBp's dimerization states and the in vivo activities of DctB. Using site‐directed mutagenesis, we identify important determinants for signal perception and transduction. Specifically, we show that the ligand‐binding pocket is essential for DCA‐induced ‘on’ activity of DctB. Mutations at different sections of DctBp's dimerization interface can lock full‐length DctB at either ‘on’ or ‘off’ state, independent of ligand binding. Taken together, these results suggest that DctBp's signal perception and transduction occur through a ‘ligand‐induced dimeric switch’, in which the changes in the dimeric conformations upon ligand binding are responsible for the signal transduction in DctB.     

Theoretical investigation on the diatomic ligand migration process and ligand binding properties in non‐O2‐binding H‐NOX domain

The Nostoc sp (Ns) H‐NOX (heme‐nitric oxide or OXygen‐binding) domain shares 35% sequence identity with soluble guanylate cyclase (sGC) and exhibits similar ligand binding property with the sGC. Previously, our molecular dynamic (MD) simulation work identified that there exists a Y‐shaped tunnel system hosted in the Ns H‐NOX interior, which servers for ligand migration. The tunnels were then confirmed by Winter et al. [PNAS 2011;108(43):E 881–889] recently using x‐ray crystallography with xenon pressured conditions. In this work, to further investigate how the protein matrix of Ns H‐NOX modulates the ligand migration process and how the distal residue composition affects the ligand binding prosperities, the free energy profiles for nitric oxide (NO), carbon monooxide (CO), and O₂ migration are explored using the steered MDs simulation and the ligand binding energies are calculated using QM/MM schemes. The potential of mean force profiles suggest that the longer branch of the tunnel would be the most favorable route for NO migration and a second NO trapping site other than the distal heme pocket along this route in the Ns H‐NOX was identified. On the contrary, CO and O₂ would prefer to diffuse via the shorter branch of the tunnel. The QM/MM (quantum mechanics/molecular mechanics) calculations suggest that the hydrophobic distal pocket of Ns H‐NOX would provide an approximately vacuum environment and the ligand discrimination would be determined by the intrinsic binding properties of the diatomic gas ligand to the heme group. Proteins 2013; 81:1363–1376. © 2013 Wiley Periodicals, Inc.     

Mapping of ligand‐binding cavities in proteins

The complex interactions between proteins and small organic molecules (ligands) are intensively studied because they play key roles in biological processes and drug activities. Here, we present a novel approach to characterize and map the ligand‐binding cavities of proteins without direct geometric comparison of structures, based on Principal Component Analysis of cavity properties (related mainly to size, polarity, and charge). This approach can provide valuable information on the similarities and dissimilarities, of binding cavities due to mutations, between‐species differences and flexibility upon ligand‐binding. The presented results show that information on ligand‐binding cavity variations can complement information on protein similarity obtained from sequence comparisons. The predictive aspect of the method is exemplified by successful predictions of serine proteases that were not included in the model construction. The presented strategy to compare ligand‐binding cavities of related and unrelated proteins has many potential applications within protein and medicinal chemistry, for example in the characterization and mapping of “orphan structures”, selection of protein structures for docking studies in structure‐based design, and identification of proteins for selectivity screens in drug design programs. Proteins 2010. © 2009 Wiley‐Liss, Inc.