Characterization of surface binding sites in glycoside hydrolases: A computational study

Structural properties of carbohydrate surface binding sites (SBSs) were investigated with computational methods. Eighty‐five SBSs of 44 enzymes in 119 Protein Data Bank (PDB) files were collected as a dataset. On the basis of SBSs shape, they were divided into 3 categories: flat surfaces, clefts, and cavities (types A, B, and C, respectively). Ligand varieties showed the correlation between shape of SBSs and ligands size. To reduce cut‐off differences in each SBSs with different ligand size, molecular docking were performed. Molecular docking results were used to refine SBSs classification and binding sites cut‐off. Docking results predicted putative ligands positions and displayed dependence of the ligands binding mode to the structural type of SBSs. Physicochemical properties of SBSs were calculated for all docking results with YASARA Structure. The results showed that all SBSs are hydrophilic, while their charges could vary and depended to ligand size and defined cut‐off. Surface binding sites type B had highest average values of solvent accessible surface area. Analysis of interactions showed that hydrophobic interactions occur more than hydrogen bonds, which is related to the presence of aromatic residues and carbohydrates interactions.     

Steric effects on multivalent ligand-receptor binding: exclusion of ligand sites by bound cell surface receptors.

Steric effects can influence the binding of a cell surface receptor to a multivalent ligand. To account for steric effects arising from the size of a receptor and from the spacing of binding sites on a ligand, we extend a standard mathematical model for ligand-receptor interactions by introducing a steric hindrance factor. This factor gives the fraction of unbound ligand sites that are accessible to receptors, and thus available for binding, as a function of ligand site occupancy. We derive expressions for the steric hindrance factor for various cases in which the receptor covers a compact region on the ligand surface and the ligand expresses sites that are distributed regularly or randomly in one or two dimensions. These expressions are relevant for ligands such as linear polymers, proteins, and viruses. We also present numerical algorithms that can be used to calculate steric hindrance factors for other cases. These theoretical results allow us to quantify the effects of steric hindrance on ligand-receptor kinetics and equilibria.     

Conformational Transitions in 3D Model of Bovine Testicular Hyaluronidase during Molecular Docking with Glycosaminoglycan Ligands

In silico molecular docking of the trimer repeating unit of chondroitin sulfate (sulfated hexasaccharide) and tetramer repeating unit of heparin (sulfated octasaccharide) to the 3D model of bovine testicular hyaluronidase by the methods of computational chemistry demonstrated the presence of eight significant binding sites for these ligands (cs1–cs8). The interaction of the active site of the enzyme with the heparin ligand, which inactivates the enzyme, and the protective effect of the chondroitin sulfate ligands bound to the surface sites of the biocatalyst molecule were theoretically studied using calculation approaches. We sequentially determined binding sites for the chondroitin sulfate ligands (in positions cs2, cs4, cs7, cs8 or cs1, cs2, cs4, cs7, cs8) critical for the protein structure stabilization, whose occupancy is theoretically sufficient to prevent irreversible deformations of the enzyme molecule when the heparin ligand is introduced into the cavity of its active site. Theoretical detection of these ‘sensibility points’ on the hyaluronidase globule indicates the possibility of regulating its functioning under the binding of the glycosaminoglycan ligands that initiate the fine formation of an effective type of the surface electrostatic potential. The interaction of the glycosaminoglycan ligands with hyaluronidase is mainly determined by electrostatic forces.     

Bromodomains识别并结合配体的机理研究

关于酵母重组蛋白内的Bromodomain识别乙酰化赖氨酸的研究近年来受到广泛的关注,但是其识别配体并与之相紧密结合的机理有待进一步的研究。本文采用2015年开发的结合位点拓扑学方法(FCTM)和分子动力学模拟的方式对Bromodomains识别并结合配体的机理进行了充分研究,其中分子动力学模拟时间达24 ns。通过FCTM方法发现结合位点的几何结构具有高度的凹性,且其alphaspace达到了131。分子动力学模拟的结果显示:在模拟的过程中结合位点表面的脯氨酸(Pro66)始终对配体保持着强的分子间相互作用,同时pocket内的水分子分布对配体的氢键网络也一直存在影响。以上结果表明Bromodomains识别并结合配体有两个重要因素:蛋白结构域自身的几何结构和配体受到来自于结合位点表面的氨基酸分子相互作用和pocket内水分子的氢键网络作用。     

β-Lactoglobulin's Conformational Requirements for Ligand Binding at the Calyx and the Dimer Interphase: a Flexible Docking Study

β-lactoglobulin (BLG) is an abundant milk protein relevant for industry and biotechnology, due significantly to its ability to bind a wide range of polar and apolar ligands. While hydrophobic ligand sites are known, sites for hydrophilic ligands such as the prevalent milk sugar, lactose, remain undetermined. Through the use of molecular docking we first, analyzed the known fatty acid binding sites in order to dissect their atomistic determinants and second, predicted the interaction sites for lactose with monomeric and dimeric BLG. We validated our approach against BLG structures co-crystallized with ligands and report a computational setup with a reduced number of flexible residues that is able to reproduce experimental results with high precision. Blind dockings with and without flexible side chains on BLG showed that: i) 13 experimentally-determined ligands fit the calyx requiring minimal movement of up to 7 residues out of the 23 that constitute this binding site. ii) Lactose does not bind the calyx despite conformational flexibility, but binds the dimer interface and an alternate Site C. iii) Results point to a probable lactolation site in the BLG dimer interface, at K141, consistent with previous biochemical findings. In contrast, no accessible lysines are found near Site C. iv) lactose forms hydrogen bonds with residues from both monomers stabilizing the dimer through a claw-like structure. Overall, these results improve our understanding of BLG''s binding sites, importantly narrowing down the calyx residues that control ligand binding. Moreover, our results emphasize the importance of the dimer interface as an insufficiently explored, biologically relevant binding site of particular importance for hydrophilic ligands. Furthermore our analyses suggest that BLG is a robust scaffold for multiple ligand-binding, suitable for protein design, and advance our molecular understanding of its ligand sites to a point that allows manipulation to control binding.     

Allosteric sites can be identified based on the residue–residue interaction energy difference

Allosteric drugs act at a distance to regulate protein functions. They have several advantages over conventional orthosteric drugs, including diverse regulation types and fewer side effects. However, the rational design of allosteric ligands remains a challenge, especially when it comes to the identification allosteric binding sites. As the binding of allosteric ligands may induce changes in the pattern of residue–residue interactions, we calculated the residue–residue interaction energies within the allosteric site based on the molecular mechanics generalized Born surface area energy decomposition scheme. Using a dataset of 17 allosteric proteins with structural data for both the apo and the ligand‐bound state available, we used conformational ensembles generated by molecular dynamics simulations to compute the differences in the residue–residue interaction energies in known allosteric sites from both states. For all the known sites, distinct interaction energy differences (>25%) were observed. We then used CAVITY, a binding site detection program to identify novel putative allosteric sites in the same proteins. This yielded a total of 31 “druggable binding sites,” of which 21 exhibited >25% difference in residue interaction energies, and were hence predicted as novel allosteric sites. Three of the predicted allosteric sites were supported by recent experimental studies. All the predicted sites may serve as novel allosteric sites for allosteric ligand design. Our study provides a computational method for identifying novel allosteric sites for allosteric drug design. Proteins 2015; 83:1375–1384. © 2014 Wiley Periodicals, Inc.     

Protection of all cleavable sites of DNA with the multiple CGCG or continuous CGG sites from the restriction enzyme,indicative of simultaneous binding of small ligands

The anthracenone ligands (1–12) with a keto-phenol and a hydroxamic acid unit were synthesized and evaluated by a restriction enzyme inhibition assay. DNA substrates composed of multiple CGCG or CGG sites are fully hydrolyzed by a restriction enzyme that is selective for each sequence. Under such conditions, the full-length DNA substrate remains only when the ligand binds to all binding sites and protects it from hydrolysis by the restriction enzymes. In the assay using AccII and the 50-mer DNA substrates containing a different number of CGCG sites at different non-binding AT base pair intervals, the more the CGCG sites, the more the full-length DNA increased. Namely, simultaneous binding of the ligand (5) to the CGCG sites increased in the order of (CGCG)₅>(CGCG)₂>(CGCG)₁. Furthermore, the length of the spacer of the hydroxamic acid to the anthracenone skeleton played an important role in the preference for the number of the d(A/T) base pairs between the CGCG sites. The long spacer-ligand (5) showed a preference to the CGCG sites with five AT pairs, and the short spacer-ligand (10) to that with two AT pairs. The ligand (12) with the shortest spacer showed a preference in simultaneous binding to the 54-mer DNA composed of 16 continuous CGG sites in the assay using the restriction enzyme Fnu4HI that hydrolyzes the d(GCGGC)/d(CGCCG) site. Application of these ligands to biological systems including the repeat DNA sequence should be of significant interest.     

De novo protein fold families expand the designable ligand binding site space

A major challenge in designing proteins de novo to bind user-defined ligands with high affinity is finding backbones structures into which a new binding site geometry can be engineered with high precision. Recent advances in methods to generate protein fold families de novo have expanded the space of accessible protein structures, but it is not clear to what extend de novo proteins with diverse geometries also expand the space of designable ligand binding functions. We constructed a library of 25,806 high-quality ligand binding sites and developed a fast protocol to place (“match”) these binding sites into both naturally occurring and de novo protein families with two fold topologies: Rossman and NTF2. Each matching step involves engineering new binding site residues into each protein “scaffold”, which is distinct from the problem of comparing already existing binding pockets. 5,896 and 7,475 binding sites could be matched to the Rossmann and NTF2 fold families, respectively. De novo designed Rossman and NTF2 protein families can support 1,791 and 678 binding sites that cannot be matched to naturally existing structures with the same topologies, respectively. While the number of protein residues in ligand binding sites is the major determinant of matching success, ligand size and primary sequence separation of binding site residues also play important roles. The number of matched binding sites are power law functions of the number of members in a fold family. Our results suggest that de novo sampling of geometric variations on diverse fold topologies can significantly expand the space of designable ligand binding sites for a wealth of possible new protein functions.     

Localization of the delta-like-1-binding site in human Notch-1 and its modulation by calcium affinity

The Notch signaling pathway plays a key role in a myriad of cellular processes, including cell fate determination. Despite extensive study of the downstream consequences of receptor activation, very little molecular data are available for the initial binding event between the Notch receptor and its ligands. In this study, we have expressed and purified a natively folded wild-type epidermal growth factor-like domain (EGF) 11-14 construct from human Notch-1 and have used flow cytometry and surface plasmon resonance analysis to demonstrate a calcium-dependent interaction with the human ligand Delta-like-1. Site-directed mutagenesis of three of the calcium-binding sites within the Notch-(11-14) fragment indicated that only loss of calcium binding to EGF12, and not EGF11 or EGF13, abrogates ligand binding. Further mapping of the ligand-binding site within this region by limited proteolysis of Notch wild-type and mutant fragments suggested that EGF12 rather than EGF11 contains the major Delta-like-1-binding site. Analysis of an extended fragment EGF-(10-14), where EGF11 is placed in a native context, surprisingly demonstrated a reduction in ligand binding, suggesting that EGF10 modulates binding by limiting access of ligand. This inhibition could be overcome by the introduction of a calcium binding mutation in EGF11, which decouples the EGF-(10-11) module interface. This study therefore demonstrates that long range calcium-dependent structural perturbations can influence the affinity of Notch for its ligand, in the absence of any post-translational modifications.     

Static and dynamic roles of extracellular loops in G-protein-coupled receptors: a mechanism for sequential binding of thyrotropin-releasing hormone to its receptor.

Small ligands generally bind within the seven transmembrane-spanning helices of G-protein-coupled receptors, but their access to the binding pocket through the closely packed loops has not been elucidated. In this work, a model of the extracellular loops of the thyrotropin-releasing hormone (TRH) receptor (TRHR) was constructed, and molecular dynamics simulations and quasi-harmonic analysis have been performed to study the static and dynamic roles of the extracellular domain. The static analysis based on curvature and electrostatic potential on the surface of TRHR suggests the formation of an initial recognition site between TRH and the surface of its receptor. These results are supported by experimental evidence. A quasi-harmonic analysis of the vibrations of the extracellular loops suggest that the low-frequency motions of the loops will aid the ligand to access its transmembrane binding pocket. We suggest that all small ligands may bind sequentially to the transmembrane pocket by first interacting with the surface binding site and then may be guided into the transmembrane binding pocket by fluctuations in the extracellular loops.     

Automated site-directed drug design: a general algorithm for knowledge acquisition about hydrogen-bonding regions at protein surfaces

This is the first of four papers that begin to explore the possibility of automated site-directed drug design. A general outline is given of the logical steps involved in approaching the problem. The starting point is the process of knowledge acquisition about the site. An algorithm is described here for the construction of a map of hydrogen-bonding regions at protein surfaces directly from the Brookhaven Protein Data Bank coordinates. Hydrogen-bonding atoms are located, intramolecular bonds are searched for, hydrogen-bonding atoms at the surface are found and hydrogen-bonding regions are computed at the accessible surface. A grid is placed within each region discovered and the probability of hydrogen bonding at each grid point is computed. The output of the program is a map of hydrogen-bonding regions displayed within a user-defined window. This information can be used as part of a knowledge base for the automatic construction of novel ligands to fit specified binding sites.     

Hydration energy landscape of the active site cavity in cytochrome P450cam

Hydration of protein cavities influences protein stability, dynamics, and function. Protein active sites usually contain water molecules that, upon ligand binding, are either displaced into bulk solvent or retained to mediate protein–ligand interactions. The contribution of water molecules to ligand binding must be accounted for to compute accurate values of binding affinities. This requires estimation of the extent of hydration of the binding site. However, it is often difficult to identify the water molecules involved in the binding process when ligands bind on the surface of a protein. Cytochrome P450cam is, therefore, an ideal model system because its substrate binds in a buried active site, displacing partially disordered solvent, and the protein is well characterized experimentally. We calculated the free energy differences for having five to eight water molecules in the active site cavity of the unliganded enzyme from molecular dynamics simulations by thermodynamic integration employing a three-stage perturbation scheme. The computed free energy differences between the hydration states are small (within 12 kJ mol⁻¹) but distinct. Consistent with the crystallographic determination and studies employing hydrostatic pressure, we calculated that, although ten water molecules could in principle occupy the volume of the active site, occupation by five to six water molecules is thermodynamically most favorable. Proteins 32:381–396, 1998. © 1998 Wiley-Liss, Inc.     

Unbinding pathways of an agonist and an antagonist from the 5-HT3 receptor

The binding sites of 5-HT3 and other Cys-loop receptors have been extensively studied, but there are no data on the entry and exit routes of ligands for these sites. Here we have used molecular dynamics simulations to predict the pathway for agonists and antagonists exiting from the 5-HT3 receptor binding site. The data suggest that the unbinding pathway follows a tunnel at the interface of two subunits, which is approximately 8 A long and terminates approximately 20 A above the membrane. The exit routes for an agonist (5-HT) and an antagonist (granisetron) were similar, with trajectories toward the membrane and outward from the ligand binding site. 5-HT appears to form many hydrogen bonds with residues in the unbinding pathway, and experiments show that mutating these residues significantly affects function. The location of the pathway is also supported by docking studies of granisetron, which show a potential binding site for granisetron on the unbinding route. We propose that leaving the binding pocket along this tunnel places the ligands close to the membrane and prevents their immediate reentry into the binding pocket. We anticipate similar exit pathways for other members of the Cys-loop receptor family.     

Studies on the binding sites of IgG2 monoclonal antibodies recognized by terpyridine‐based affinity ligands

This investigation has examined the origin of the molecular recognition associated with the interaction of monoclonal IgG2's with terpyridine‐based ligands immobilized onto agarose‐derived chromatographic adsorbents. Isothermal titration calorimetric (ITC) methods have been employed to acquire thermodynamic data associated with the IgG2‐ligand binding. These ITC investigations have documented that different enthalpic and entropic processes are involved depending on the nature of the chemical substituents in the core structure of the terpyridinyl moiety. In addition, molecular docking studies have been carried out with IgG2 structures with the objective to identify possible ligand binding sites and key interacting amino acid residues. These molecular docking experiments with the different terpyridine‐based ligands have shown that all of the examined ligands can potentially undergo favorable interactions with a site located within the Fab region of the IgG2. However, another favorable binding site was also identified from the docking poses to exist within the Fc region of the IgG2 for some, but not all, of the ligands studied. These investigations have provided a basis to elucidate the unique binding properties and chromatographic behaviors shown by several substituted terpyridine ligands in their interaction with IgGs of different isotype. Copyright © 2016 John Wiley & Sons, Ltd.     

Identification of a modifier site on acetylcholinesterase with affinity for d-tubocurarine

Decamethonium and d-tubocurarine displace N-methylacridinium ion, a potent fluorescent inhibitor of acetylcholinesterase, from the surface of the enzyme. Decamethonium is competitive with N-methylacridinium which indicates that the binding sites for these ligands overlap. However, the displacement of N-methylacridinium ion by d-tubocurarine requires the existence of a binding site for d-tubocurarine in addition to the active site. Since the affinities for d-tubocurarine at both sites are comparable, two well defined ligand binding sites must exist for each catalytic site that is titratable by 7-dimethylcarbamyl-N-methylquinolinium iodide.     

Analysis of the structure and dynamics of human serum albumin

A computational approach to designing a peptide-based ligand for the purification of human serum albumin (HSA) was undertaken using molecular docking and molecular dynamics (MD) simulation. A three-step procedure was performed to design a specific ligand for HSA. Based on the candidate pocket structure of HSA (warfarin binding site), a peptide library was built. These peptides were then docked into the pocket of HSA using the GOLD program. The GOLDscore values were used to determine the affinity of peptides for HSA. Consequently, the dipeptide Trp–Trp, which shows a high GOLDscore value, was selected and linked to a spacer arm of Lys[CO(CH₂)₅NH] on the surface of ECH-lysine sepharose 4 gel. For further evaluation, the Autodock Vina program was used to dock the linked compound into the pocket of HSA. The docking simulation was performed to obtain a first guess of the binding structure of the spacer–Trp–Trp–HSA complex and subsequently analyzed by MD simulations to assess the reliability of the docking results. These MD simulations indicated that the ligand–HSA complex remains stable, and water molecules can bridge between the ligand and the protein by hydrogen bonds. Finally, absorption spectroscopic studies were performed to illustrate the appropriateness of the binding affinity of the designed ligand toward HSA. These studies demonstrate that the designed dipeptide can bind preferentially to the warfarin binding site. Graphical Abstract

Three-step computational approach to the design of a dipeptide ligand for human serum albumin purification exploiting structure-based docking and molecular dynamics simulation