Electrostatic complementarity between proteins and ligands. 2. Ligand moieties

Summary Drug design strategies consider factors governing intermolecular interactions to build up putative ligands. In many strategies, the ligand is constructed using fragments which are placed in the site sequentially. The optimization is then performed with each fragment. We would like to examine if this optimization strategy could generate ligands with optimal electrostatic interactions. The electrostatic complementarities between constituent moieties and the receptor site have been calculated. The whole-ligand complementarity does not appear to be the mathematical mean of the individual complementarities, nor have we found a simple relationship between the moiety and whole-ligand complementarities. The results demonstrate clearly that, using a simple model, it is very difficult to predict the electrostatic potential complementarity of the whole ligand from the complementarities of its constituent chemical moieties. This means that ligand design strategies must optimize the electrostatic complementarity globally, and not moiety by moiety.     

Electrostatic complementarity between proteins and ligands. 3. Structural basis

Summary Electrostatic potential complementarity between ligands and their receptor sites is evaluated by the superposition of the electrostatic potential, generated by the receptor, onto the ligand potential over the ligand van der Waals surface. We would like to examine which structural factors generate this pattern of superposition. Example studies suggest that in many ligand-protein pairs, there exist principal formal charges on each molecule, largely responsible for the electrostatic potential complementarity observed. Electrostatic potential complementarity depends on the relative disposition of these principal charges and the ligand van der Waals surface. Simple mathematical models were constructed to predict the complementarity solely from structural considerations. The essential conditions for electrostatic potential complementarity were elucidated. These can be used in ligand design strategies to obtain an electrostatically optimal ligand.     

Improvement of Ligand Affinity and Thermodynamic Properties by NMR‐Based Evaluation of Local Dynamics and Surface Complementarity in the Receptor‐Bound State

The thermodynamic properties of a ligand in the bound state affect its binding specificity. Strict binding specificity can be achieved by introducing multiple spatially defined interactions, such as hydrogen bonds and van der Waals interactions, into the ligand–receptor interface. These introduced interactions are characterized by restricted local dynamics and improved surface complementarity in the bound state. In this study, we experimentally evaluated the local dynamics and the surface complementarity of weak‐affinity ligands in the receptor‐bound state by forbidden coherence transfer analysis in free‐bound exchange systems (Ex‐FCT), using the interaction between a ligand, a myocyte‐enhancer factor 2A (MEF2A) docking peptide, and a receptor, p38α, as a model system. The Ex‐FCT analyses successfully provided information for the rational design of a ligand with higher affinity and preferable thermodynamic properties for p38α.     

Ligand atom partial charges assignment for complementary electrostatic potentials

Summary The design of molecules to fit into the active site of receptors is a rapidly developing area of pharmacology and medicinal chemistry. A good ligand needs a suitable geometry and also appropriate electrostatic properties. The electrostatic properties of the ligand should complement those of the receptor. We present a method for the assignment of atom-centred point charges for a ligand, based on the electrostatic potential of the receptor. These point charges are chosen to give the best possible complementarity to the receptor electrostatic potential over the van der Waals surface of the ligand. We demonstrate that point charges can be chosen to give good electrostatic complementarity, and suggest that a molecule with similar electrostatic properties should bind well to the receptor.     

Molecular Simulations using Spherical Harmonics

Computer-aided drug design is to develop a chemical that binds to a target macromolecule known to play a key role in a disease state. In recognition of ligands by their protein receptors, molecular surfaces are often used because they represent the in-teracting part of molecules and they should reflex the comple-mentarity between ligand and receptor. However, assessing the surface complementarity by searching all relative position of two surfaces is often computationally expensive. The comple-mentarity of lobe-hole is very important in protein-ligand inter-actions. Spherical harmonic models based on expansions of spherical harmonic functions were used as a f‘mgerprint to ap-proximate the binding cavity and the ligand, respectively. This defines a new way to identify the complementarity between lobes and holes. The advantage of this method is that two spherical harmonic surfaces to be compared can be defined sep-arately. This method can be used as a filter to eliminate candi-dates among a large number of conformations, and it will speed up the docking procedure. Therefore, it is possible to select complementary ligands or complementary conformations of a ligand and the macromolecules, by comparing their fingerprints previously stored in a database.     

A shape- and chemistry-based docking method and its use in the design of HIV-1 protease inhibitors

Summary The program DOCK [1,2] has been used successfully to identify molecules which will bind to a specified receptor [3]. The original method ranks molecules based on their shape complementarity to the receptor site and relies on the chemist to bring the appropriate electrostatic or hydrogen bond properties into the molecular skeletons obtained in the search. This is useful when screening a small database of compounds, where it is not likely that molecules with both the correct shape and electrostatic properties will be found. As large databases are more likely to have redundant molecular shapes with a variety of functionality (e.g., members of a congeneric series), it would be useful to have a method which identifies molecules with both the correct shape and functionality. To this end we have modified the DOCK 1.0 method to target user-specified atom types to selected positions in the receptor site. The target sites can be chosen based on structural evidence, calculation or inspection. Targeted-DOCK improves the ability of the DOCK method to find the crystallographically determined binding mode of a ligand. Additionally, targeted-DOCK searches a database of small molecules at 100–1000 times the rate of DOCK 1.0, allowing more molecules to be screened and more sophisticated scoring schemes to be employed. Targeted-DOCK has been used successfully in the design of a novel non-peptide inhibitor of HIV-1 protease.     

Experimental and theoretical evidence of overcharging of calcium silicate hydrate

Electrokinetic measurements such as electrophoresis may show an inversion of the effective surface charge of colloidal particle called overcharging. This phenomenon has been studied by various theoretical approaches but up to now very few attempts of confrontation between theory and experiment have been conducted. In this work we report electrophoretic measurements as well as Monte Carlo simulations of the electrokinetic potential for the surface of calcium silicate hydrate (CSH), which is the major constituent of hydrated cement. In the simulations, the surface charge of CSH nanoparticles in equilibrium with the ionic solution is determined by a single site characteristic and electrostatic interactions between all explicit charges at the surface and in the electric double layer. We will show that ordinary electrostatic interactions are enough to describe all experimental observations. Actually, an excellent agreement is found between experimental and simulated results without any fitting parameter, both with respect to surface titration and electrokinetic behaviour. The agreement extends over a wide range of electrostatic coupling, from a weakly charged surface with mainly monovalent counter-ions to a highly charged one with divalent counter-ions.     

Molecular surface-volume and property matching to superpose flexible dissimilar molecules

Summary Steric complementarity is a prerequisite for ligand-receptor recognition; this implies that drugs with a common receptor binding site should possess sterically similar binding surfaces. This principle is used as the basis for an automatic and unbiased method that superposes molecules. One molecule is rotated and translated to maximize the overlap between the two molecular surface volumes. A fast grid-based method is used to determine the extent of this overlap, and this is optimized using simulated annealing. Matches with high steric similarity scores are then sorted on the basis of both hydrogen-bond and electrostatic similarity between the matched molecules. Flexible molecules are treated as a set of rigid representative conformers. The algorithm has correctly predicted superpositions between a number of pairs of molecules, according to crystallographic data from ligands that have been co-crystallized at common enzyme binding sites.     

SARS冠状病毒主要蛋白酶二聚体静电和疏水效应的研究

从三种冠状病毒主要蛋白酶SARS 3CL, HCoV 3CL和TGEC 3CL蛋白酶结构出发,着重研究了三种蛋白酶二聚体单体之间的静电和疏水相互作用.用连续介质模型有限差分方法计算得到三种蛋白二聚体界面处的静电势,发现三种蛋白酶单体和单体之间静电势分布具有明显的互补性,三种蛋白酶二聚体单体之间具有相同的静电相互作用能.用溶剂可及表面积模型分析了分子表面积及疏水性,发现三种蛋白酶具有相同的疏水分布,其中SARS 3CL蛋白酶疏水率为74%,驱动其单体聚合成二聚体.对三种蛋白酶的去溶剂化能疏水项的计算表明,三种蛋白酶二聚体单体之间具有相似的疏水相互作用能.     

Monitoring the dynamics of ligand-receptor complexes on model membranes

Ligand-induced cross-linking of cell surface receptors is a basic paradigm of signal activation by many transmembrane receptors. After ligand binding, the receptor complexes formed on the membrane are dynamically maintained by two-dimensional protein-protein interactions on the membrane. The biophysical principles governing the dynamics of such interactions have not been understood, mainly because the measurement of lateral interactions on membranes so far has not been experimentally addressed. Here, we describe a generic approach for measuring two-dimensional dissociation rate constants in vitro using a novel high-affinity chelator lipid for reconstituting a ternary cytokine-receptor complex on solid-supported membranes. While monitoring the interaction between the ligand and one of the receptor subunits on the membrane by fluorescence resonance energy transfer, the equilibrium on the surface was perturbed by rapidly tethering a large excess of the unlabeled receptor subunit. Displacement of labeled by unlabeled protein in the ternary complex was detected as a recovery of the donor quenching. Since the dissociation of the ligand-receptor complex in plane of the membrane was the rate-limiting step under these conditions, the two-dimensional rate constant of this process was determined. Strikingly, the two-dimensional dissociation was much slower than ligand dissociation into solution, suggesting that membrane tethering significantly affects the dissociation process. This result highlights the importance of studying ligand-receptor complexes tethered to membranes for understanding the principles governing signal activation by ligand-induced receptor assembling.     

A fast and efficient method to generate biologically relevant conformations

Summary Mutual binding between a ligand of low molecular weight and its macromolecular receptor demands structural complementarity of both species at the recognition site. To predict binding properties of new molecules before synthesis, information about possible conformations of drug molecules at the active site is required, especially if the 3D structure of the receptor is not known. The statistical analysis of small-molecule crystal data allows one to elucidate conformational preferences of molecular fragments and accordingly to compile libraries of putative ligand conformations. A comparison of geometries adopted by corresponding fragments in ligands bound to proteins shows similar distributions in conformation space. We have developed an automatic procedure that generates different conformers of a given ligand. The entire molecule is decomposed into its individual ring and open-chain torsional fragments, each used in a variety of favorable conformations. The latter ones are produced according to the library information about conformational preferences. During this building process, an extensive energy ranking is applied. Conformers ranked as energetically favorable are subjected to an optimization in torsion angle space. During minimization, unfavorable van der Waals interactions are removed while keeping the open-chain torsion angles as close as possible to the experimentally most frequently observed values. In order to assess how well the generated conformers map conformation space, a comparison with experimental data has been performed. This comparison gives some confidence in the efficiency and completeness of this approach. For some ligands that had been structurally characterized by protein crystallography, the program was used to generate sets of some 10 to 100 conformers. Among these, geometries are found that fall convincingly close to the conformations actually adopted by these ligands at the binding site.     

囊泡表面分子间通讯交流体系的构建: 利用受体的分子识别行为调控酶的活性

在人工双层膜囊泡表面, 构建了一个通过人工受体的分子识别行为控制酶反应活性的超分子体系. 体系以生物体细胞信号转导系统为模拟原型, 由作为受体的烷基胺、被受体识别的信号分子吡哆醛衍生物、乳酸脱氢酶、受体和酶之间的媒介物Cu^2＋以及作为体系载体的合成肽脂囊泡五个成分构成．通过UV-vis光谱法及动态光散射测定对体系进行了评价, 结果表明: 随着受体疏水参数增大, 其对信号分子的识别能力增强, 二者呈良好的线性关系; 通过信号分子与囊泡表面静电相互作用的研究表明信号分子具有选择性; 媒介物与信号分子–受体可形成化学计量比为1∶2的配合物, 其形成能力比媒介物与酶的结合能力更强．作为结论, 体系中烷基胺受体对磷酸吡哆醛信号分子的识别有效控制了处于囊泡表面的乳酸脱氢酶的活性．     

Electrostatics of ligand binding: parametrization of the generalized Born model and comparison with the Poisson-Boltzmann approach

An accurate and fast evaluation of the electrostatics in ligand-protein interactions is crucial for computer-aided drug design. The pairwise generalized Born (GB) model, a fast analytical method originally developed for studying the solvation of organic molecules, has been widely applied to macromolecular systems, including ligand-protein complexes. However, this model involves several empirical scaling parameters, which have been optimized for the solvation of organic molecules, peptides, and nucleic acids but not for energetics of ligand binding. Studies have shown that a good solvation energy does not guarantee a correct model of solvent-mediated interactions. Thus, in this study, we have used the Poisson-Boltzmann (PB) approach as a reference to optimize the GB model for studies of ligand-protein interactions. Specifically, we have employed the pairwise descreening approximation proposed by Hawkins et al.(1) for GB calculations and DelPhi for PB calculations. The AMBER all-atom force field parameters have been used in this work. Seventeen protein-ligand complexes have been used as a training database, and a set of atomic descreening parameters has been selected with which the pairwise GB model and the PB model yield comparable results on atomic Born radii, the electrostatic component of free energies of ligand binding, and desolvation energies of the ligands and proteins. The energetics of the 15 test complexes calculated with the GB model using this set of parameters also agrees well with the energetics calculated with the PB method. This is the first time that the GB model has been parametrized and thoroughly compared with the PB model for the electrostatics of ligand binding.     

Charge optimization of the interface between protein kinases and their ligands

Examining the potential for electrostatic complementarity between a ligand and a receptor is a useful technique for rational drug design, and can demonstrate how a system prioritizes interactions when allowed to optimize its charge distribution. In this computational study, we implemented the previously developed, continuum solvent-based charge optimization theory with a simple, quadratic programming algorithm and the UHBD Poisson-Boltzmann solver. This method allows one to compute the best set of point charges for a ligand or ligand region based on the ligand and receptor shape, and the receptor partial charges, by optimizing the binding free energy obtained from a continuum-solvent model. We applied charge optimization to a fragment of the heat-stable protein kinase inhibitor (PKI) of protein kinase A (PKA), to three flavopiridol inhibitors of CDK2, and to cyclin A which interacts with CDK2 to regulate the cell cycle. We found that a combination of global (involving every charge) and local (involving only charges in a local region) optimization can give useful hints for designing better inhibitors. Although some parts of an inhibitor may already contribute significantly to binding, we found that they could still be the most important targets for modifications to obtain stronger binders. In studying the binding of flavopiridol inhibitors to CDK2, comparable binding affinity could be obtained regardless of whether the net charges of the inhibitors were constrained to -2, -1, 0, 1, or 2 during the optimization. This provides flexibility in inhibitor design when a certain net charge of the inhibitor is desired in addition to strong binding affinity. For the study of the PKA-PKI and CDK2-cyclin A interfaces, we identified residues whose charge distributions are already close to optimal and those whose charge distributions could be refined to further improve binding.     

Molecular electrostatic potentials and hydrogen bonding in alpha-, beta-, and gamma-cyclodextrins

Cyclodextrins (CDs) are cyclic oligomers of glucose having the toroid of sugars elaborating a central cavity of varying size depending on the number of glucoses. The central hydrophobic cavity of CD shows a binding affinity toward different guest molecules, which include small substituted benzenes to long chain surfactant molecules leading to a variety of inclusion complexes when the size and shape complementarity of host and guest are compatible. Further, interaction of guest molecules with the outer surface of alpha-CD has also been observed. Primarily it is the electrostatic interactions that essentially constitute a driving force for the formation of inclusion complexes. To gain insights for these interactions, the electronic structure and the molecular electrostatic potentials in alpha-, beta-, and gamma-CDs are derived using the hybrid density functional theory employing the three-parameter exchange correlation functional due to Becke, Lee, Yang, and Parr (B3LYP). The present work demonstrates how the topography of the molecular electrostatic potential (MESP) provides a measure of the cavity dimensions and understanding of the hydrogen-bonded interactions involving primary and secondary hydroxyl groups. In alpha-CD, hydrogen-bonded interactions between primary -OH groups engender a "cone-like" structure, while in beta- or gamma-CD the interactions from the primary -OH with ether oxygen in glucose ring facilitates a "barrel-like" structure. Further, the strength of hydrogen-bonded interactions of primary -OH groups follows the rank order alpha-CD > beta-CD > gamma-CD, while the secondary hydrogen-bonded interactions exhibit a reverse trend. Thus weak hydrogen-bonded interactions prevalent in gamma-CD manifest in shallow MESP minima near hydroxyl oxygens compared to those in alpha- or beta-CD. Furthermore, electrostatic potential topography reveals that the guest molecule tends to penetrate inside the cavity forming the inclusion complex in beta- or gamma-CD.     

A computational model of the 5-HT3 receptor extracellular domain: search for ligand binding sites

A three-dimensional model of the 5-HT₃ receptor extarcellular domain has been derived on the basis of the nicotinic acetylcholine receptor model recently published by Tsigelny et al. Maximum complementarity between the position and characteristics of mutated residues putatively involved in ligand interaction and the pharmacophoric elements derived by the indirect approach applied on several series of 5-HT₃ ligands have been exploited to gain insights into the ligand binding modalities and to speculate on the mechanistic role of the structural components. The analysis of the three-dimensional model allows one to distinguish among amino acids that exert key roles in ligand interactions, subunit architecture, receptor assembly and receptor dynamics. For some of these, alternative roles with respect to the ones hypothesized by experimentalists are assigned. Different binding modalities for agonists and antagonists are highlighted, and residues which probably play a role in the transduction of binding into a change in conformational state of the receptor are suggested. Received: 27 July 2000 / Accepted: 15 September 2000 / Published online: 21 December 2000     

Colloidally Stable CdS Quantum Dots in Water with Electrostatically Stabilized Weak-Binding,Sulfur-Free Ligands

Colloidal quantum dot (QD) photocatalysts have the electrochemical and optical properties to be highly effective for a range of redox reactions. QDs are proven photo-redox catalysts for a variety of reactions in organic solvents but are less prominent for aqueous reactions. Aqueous QD photocatalysts require hydrophilic ligand shells that provide long-term colloidal stability but are not so tight-binding as to prevent catalytic substrates from accessing the QD surface. Common thiolate ligands, which also poison many co-catalysts and undergo photo-oxidative desorption, are therefore often not an option. This paper describes a framework for the design of water-solubilizing ligands that are in dynamic exchange on and off the QD surface, but still provide long-term colloidal stability to CdS QDs. The binding affinity and inter-ligand electrostatic interactions of a bifunctional ligand, aminoethyl phosphonic acid (AEP), are tuned with the pH of the dispersion. The key to colloidal stability is electrostatic stabilization of the monolayer. This work demonstrates a means of mimicking the stabilizing power of a thiolate-bound ligand with a zwitterionic tail group, but without the thiolate binding group.     

Docking of small ligands to low-resolution and theoretically predicted receptor structures

We have developed a simple docking procedure that is able to utilize low-resolution models of proteins created by structure prediction algorithms such as threading or ab initio folding to predict the conformation of receptor-small ligand complexes. In our approach, using only approximate, discretized models of both molecules, we search for the steric and quasi-chemical complementarity between a ligand and the receptor molecules. This averaging procedure allows for the compensation of numerous structural inaccuracies resulting from the theoretical predictions of the receptor structure. The best relative orientation of these two models is obtained by an exhaustive scan over the rigid body's six-dimensional translational and rotational degrees of freedom. The search method is based on a real space grid-searching algorithm, unlike docking methods based on the fast Fourier Transform algorithm. We have applied this algorithm to rebuild structures of several complexes available in the Protein Data Bank. The structures of the receptors are produced by means of our threading algorithm PROSPECTOR, subsequently refined, and then utilized in the docking experiment. In many cases, not only is the localization of the binding site on the receptor surface correctly identified, but the proper orientation of the bounded ligand is also reasonably well reproduced within the level of accuracy of the modeled receptor itself.     

Immunomodulatory peptides from IgSF proteins: a review

The immunoglobulin (Ig) domain is a highly conserved domain predominantly observed in cell surface proteins due to its ability to resist proteolysis. By mutation and selection the Ig domain has evolved to serve diverse biological functions including growth and development, signaling, adhesion and protein-carbohydrate interactions. Collectively, proteins with Ig-like domain constitute the immunoglobulin superfamily (IgSF). The IgSF proteins make up over 2% of human genes constituting the largest gene family in the human genome. Analogous to the complementarity determining regions (CDR)s that form the antigen combining sites of the antibody, the high specificity of the IgSF receptor-ligand interaction is attributed to the sequence and structure of the CDR-like regions unique to each IgSF protein. Hence, CDR-like regions provide ideal templates for the design of mimetics that can potentially perturb specific IgSF receptor/ligand interactions. The determinants of binding are localized near the CDR-like regions, conformation is determined locally and is unique for each loop. In structure based drug design one of the approaches to identify lead agents is to map the receptor/ligand binding epitope onto a small peptide. Data from theoretical, structural and functional studies have been adopted in the design of novel peptide antagonists of the IgSF protein-protein interactions. Many peptide antagonists have shown significant therapeutic potential in multiple animal models. The design of the IgSF peptide analogs, rationale as therapeutic targets, functional efficacy and the clinical benefits are reviewed here.     

Hydroxamate类抑制剂与MMP-2的分子对接研究

通过分子动力学模拟研究了MMP-2和hydroxamate抑制剂之间的作用模式。在分子动力学模拟中,对于催化区的锌离子和其共价结合的配体(包括抑制剂和组氨酸)采用了键合的模型。从模拟的结果可以看到,R^1取代基和MMP-2的S1疏水口袋中的部分残基能形成很好的几何匹配,从而可以产生很强的范德华和疏水相互作用。模拟结果也表明,两个抑制剂和MMP-2之间分别能形成5个和8个氢键,抑制剂B比A活性更高的原因就是能够形成更加有利氢键作用模式。在整个模拟过程中,催化锌都能保持好的五配位形式,配位键的长度也处于稳定的状态,预测得到的MMP-2和其抑制剂的相互作用模式对于全新抑制剂的设计提供了非常重要的结构信息。