Automatic superposition of drug molecules based on their common receptor site

Summary We have prevously developed a new rational method for superposing molecules in terms of submolecular physical and chemical properties, but not in terms of atom positions or chemical structures as has been done in the conventional methods. The program was originally developed for interactive use on a three-dimensional graphic display, providing goodness-of-fit indices on molecular shape, hydrogen bonds, electrostatic interactions and others.Here, we report a new unbiased searching method for the best superposition of molecules, covering all the superposing modes and conformational freedom, as an additional function of the program. The function is based on a novel least-squares method which superposes the expected positions and orientations of hydrogen bonding partners in the receptor that are deduced from both molecules. The method not only gives reliability and reproducibility to the result of the superposition, but also allows us to save labor and time. It is demonstrated that this method is very efficient for finding the correct superposing mode in such systems where hydrogen bonds play important roles.     

Cavity search: An algorithm for the isolation and display of cavity-like binding regions

Summary A set of algorithms designed to enhance the display of protein binding cavities is presented. These algorithms, collectively entitled CAVITY SEARCH, allow the user to isolate and fully define the extent of a particular cavity. Solid modeling techniques are employed to produce a detailed cast of the active site region, which can then be color-coded to show electrostatic and steric interactions between the protein cavity and a bound ligand.     

An automated method for predicting the positions of hydrogen-bonding atoms in binding sites

Hydrogen bonds are the most specific, and therefore predictable of the intermolecular interactions involved in ligand–protein binding. Given the structure of a molecule, it is possible to estimate the positions at which complementary hydrogen-bonding atoms could be found. Crystal-survey data are used in the design of a program, HBMAP, that generates a hydrogen-bond map for any given ligand, which contains all the feasible positions at which a complementary atom could be found. On superposition of ligands, the overlapping regions of their maps represent positions of receptor atoms to which each molecule can bind. The certainty of these positions is increased by the incorporation of a larger number and diversity of molecules. In this work, superposition is achieved using the program HBMATCH, which uses simulated annealing to generate the correspondence between points from the hydrogen-bonding maps of the two molecules. Equivalent matches are distinguished on the basis of their steric similarity. The strategy is tested on a number of ligands for which ligand–protein complexes have been solved crystallographically, which allows validation of the techniques. The receptor atom positions of thermolysin are successfully predicted when the correct superposition is obtained.     

Methods for the prediction of protein-ligand binding sites for structure-based drug design and virtual ligand screening

Structure Based Drug Design (SBDD) is a computational approach to lead discovery that uses the three-dimensional structure of a protein to fit drug-like molecules into a ligand binding site to modulate function. Identifying the location of the binding site is therefore a vital first step in this process, restricting the search space for SBDD or virtual screening studies. The detection and characterisation of functional sites on proteins has increasingly become an area of interest. Structural genomics projects are increasingly yielding protein structures with unknown functions and binding sites. Binding site prediction was pioneered by pocket detection, since the binding site is often found in the largest pocket. More recent methods involve phylogenetic analysis, identifying structural similarity with proteins of known function and identifying regions on the protein surface with a potential for high binding affinity. Binding site prediction has been used in several SBDD projects and has been incorporated into several docking tools. We discuss different methods of ligand binding site prediction, their strengths and weaknesses, and how they have been used in SBDD.     

An analysis of fractal geometry of macromolecular gelation

With fractal geometry theory and based on experiments, an analysis of fractal geometry behavior of gelation of macromolecules was carried out. Using the cross-linking copolymerization of styrene-divinylbenzene (DVB) as an example, through the determinations of the evolution of the molecular weight, size and the dependence of scattering intensity on the angle of macromolecules by employing laser and synchrotron small angle X-ray scattering, respectively, this chemical reaction was described quantitatively, its fractal behavior was analyzed and the fractal dimension was also measured. By avoiding the complex theories on gelation, this approach is based on modern physical techniques and theories to perform the analysis of the behavior of fractal geometry of macromolecular gelation and thus is able to reveal the rules of this kind of complicated gelation more essentially and profoundly.     

Gaussian mapping of chemical fragments in ligand binding sites

We present a new approach to automatically define a quasi-optimal minimal set of pharmacophoric points mapping the interaction properties of a user-defined ligand binding site. The method is based on a fitting algorithm where a grid of sampled interaction energies of the target protein with small chemical fragments in the binding site is approximated by a linear expansion of Gaussian functions. A heuristic approximation selects from this expansion the smallest possible set of Gaussians required to describe the interaction properties of the binding site within a prespecified accuracy. We have evaluated the performance of the approach by comparing the computed Gaussians with the positions of aromatic sites found in experimental protein-ligand complexes. For a set of 53 complexes, good correspondence is found in general. At a 95% significance level, approximately 65% of the predicted interaction points have an aromatic binding site within 1.5 A. We then studied the utility of these points in docking using the program DOCK. Short docking times, with an average of approximately 0.18 s per conformer, are obtained, while retaining, both for rigid and flexible docking, the ability to sample native-like binding modes for the ligand. An average 4-5-fold speed-up in docking times and a similar success rate is estimated with respect to the standard DOCK protocol.     

A comparison of progestin and androgen receptor binding using the CoMFA technique

Summary A series of 48 steroids has been studied with the SYBYL QSAR module using Relative Binding Affinities (RBAs) to progesterone and androgen receptors obtained from the literature. Models for the progesterone and androgen data were developed. Both models show regions where sterics and electrostatics correlate to binding affinity but are different for androgen and progesterone which suggests differences possibly important for receptor selectivity. The progesterone model is more predictive than the androgen (predictive r² of 0.725 vs. 0.545 for progesterone and androgen, respectively).     

Theory and simulation of diffusion-influenced, stochastically gated ligand binding to buried sites

We consider the diffusion-influenced rate coefficient of ligand binding to a site located in a deep pocket on a protein; the binding pocket is flexible and can reorganize in response to ligand entrance. We extend to this flexible protein-ligand system a formalism developed previously [A. M. Berezhkovskii, A, Szabo, and H.-X. Zhou, J. Chem. Phys. 135, 075103 (2011)] for breaking the ligand-binding problem into an exterior problem and an interior problem. Conformational fluctuations of a bottleneck or a lid and the binding site are modeled as stochastic gating. We present analytical and Brownian dynamics simulation results for the case of a cylindrical pocket containing a binding site at the bottom. Induced switch, whereby the conformation of the protein adapts to the incoming ligand, leads to considerable rate enhancement.     

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     

Ligand accessibility to receptor binding sites enhanced by movable polyrotaxanes

Functionalized polyrotaxanes are utilized to investigate the relation to multivalent interactions between the mannose moiety and Con A immobilized surfaces. According to the results of SPR spectroscopy, the mannose-conjugated polyrotaxanes show a higher response than any other mannose conjugate on both surfaces of high- and low-density Con A. Moreover, the results of the FRET analysis suggest that the mobility of α-cyclodextrins in the polyrotaxane more efficiently contributes to their binding interactions in a multivalent manner. This well-defined polyrotaxane system provides control over ligand density, ligand mobility, and gives an efficient response to the biological interaction receptor, which has not been easy to achieve in covalently bound polymeric systems.     

Alloxazine as a ligand for selective binding to adenine opposite AP sites in DNA duplexes and analysis of single-nucleotide polymorphisms

Alloxazine can bind to adenine selectively over other nucleobases opposite an abasic site in DNA duplexes (5'-TCC AGX GCA AC-3'/3'-AGG TCN CGT TG-5', X=AP site, N=A, T, C, G) with a dissociation constant of 0.82 microM (pH 7.0, I=0.11 M, at 5 degrees C), and it is applicable to SNPs typing of PCR amplification products based on the binding-induced fluorescence response.     

Improved mapping of protein binding sites

Computational mapping methods place molecular probes – small molecules or functional groups – on a protein surface in order to identify the most favorable binding positions by calculating an interaction potential. Mapping is an important step in a number of flexible docking and drug design algorithms. We have developed improved algorithms for mapping protein surfaces using small organic molecules as molecular probes. The calculations reproduce the binding of eight organic solvents to lysozyme as observed by NMR, as well as the binding of four solvents to thermolysin, in good agreement with x-ray data. Application to protein tyrosine phosphatase 1B shows that the information provided by the mapping can be very useful for drug design. We also studied why the organic solvents bind in the active site of proteins, in spite of the availability of alternative pockets that can very tightly accommodate some of the probes. A possible explanation is that the binding in the relatively large active site retains a number of rotational states, and hence leads to smaller entropy loss than the binding elsewhere else. Indeed, the mapping reveals that the clusters of the ligand molecules in the protein's active site contain different rotational-translational conformers, which represent different local minima of the free energy surface. In order to study the transitions between different conformers, reaction path and molecular dynamics calculations were performed. Results show that most of the rotational states are separated by low free energy barriers at the experimental temperature, and hence the entropy of binding in the active site is expected to be high.     

The structure-activity relationship of inhibitors of serotonin uptake and receptor binding

Summary An analysis of five different datasets of inhibitors of serotonin uptake has yielded quantitative structure/ activity relationships (QSARs) which delineate the role of steric and hydrophobic properties essential for inhibition by phenylethylamine-type analogues.     

Kinetics of ligand binding to membrane receptors from equilibrium fluctuation analysis of single binding events

Equilibrium fluctuation analysis of single binding events has been used to extract binding kinetics of ligand interactions with cell-membrane bound receptors. Time-dependent total internal reflection fluorescence (TIRF) imaging was used to extract residence-time statistics of fluorescently stained liposomes derived directly from cell membranes upon their binding to surface-immobilized antibody fragments. The dissociation rate constants for two pharmaceutical relevant antibodies directed against different B-cell expressed membrane proteins was clearly discriminated, and the affinity of the interaction could be determined by inhibiting the interaction with increasing concentrations of soluble antibodies. The single-molecule sensitivity made the analysis possible without overexpressed membrane proteins, which makes the assay attractive in early drug-screening applications.     

Extended-ensemble docking to probe dynamic variation of ligand binding sites during large-scale structural changes of proteins

Proteins can sample a broad landscape as they undergo conformational transition between different functional states. At the same time, as key players in almost all cellular processes, proteins are important drug targets. Considering the different conformational states of a protein is therefore central for a successful drug-design strategy. Here we introduce a novel docking protocol, termed extended-ensemble docking, pertaining to proteins that undergo large-scale (global) conformational changes during their function. In its application to multidrug ABC-transporter P-glycoprotein (Pgp), extensive non-equilibrium molecular dynamics simulations employing system-specific collective variables are first used to describe the transition cycle of the transporter. An extended set of conformations (extended ensemble) representing the full transition cycle between the inward- and the outward-facing states is then used to seed high-throughput docking calculations of known substrates, non-substrates, and modulators of the transporter. Large differences are predicted in the binding affinities to different conformations, with compounds showing stronger binding affinities to intermediate conformations compared to the starting crystal structure. Hierarchical clustering of the binding modes shows all ligands preferably bind to the large central cavity of the protein, formed at the apex of the transmembrane domain (TMD), whereas only small binding populations are observed in the previously described R and H sites present within the individual TMD leaflets. Based on the results, the central cavity is further divided into two major subsites, first preferably binding smaller substrates and high-affinity inhibitors, whereas the second one shows preference for larger substrates and low-affinity modulators. These central subsites along with the low-affinity interaction sites present within the individual TMD leaflets may respectively correspond to the proposed high- and low-affinity binding sites in Pgp. We propose further an optimization strategy for developing more potent inhibitors of Pgp, based on increasing its specificity to the extended ensemble of the protein, instead of using a single protein structure, as well as its selectivity for the high-affinity binding site. In contrast to earlier in silico studies using single static structures of Pgp, our results show better agreement with experimental studies, pointing to the importance of incorporating the global conformational flexibility of proteins in future drug-discovery endeavors.

Functional states of P-glycoprotein formed during its full transition cycle (red to blue), captured by molecular dynamics simulations, form a structural framework for extended-ensemble docking of small-molecule ligands of diverse activities.     

Theoretical studies on opioid receptors and ligands. I. Molecular modeling and QSAR studies on the interaction mechanism of fentanyl analogs binding to μ‐opioid receptor

Based on our previous result of the three‐dimensional model of the μ‐opioid receptor, binding conformations of 13 fentanyl analogs and three‐dimensional structures for the complexs of these analogs with μ‐opioid receptor were constructed employing the molecular modeling method and our binding conformation search program for ligands (BCSPL). Energetic calculation and quantitative structure–activity relationship (QSAR) analysis indicated a good correlation between the calculated binding energies of fentanyl analogs and their binding affinities, pK_i's and pK's, and analgesic activities, − log ED₅₀'s. Based on the three‐dimensional models, the possible interaction mechanism of fentanyl analogs with μ‐opioid receptor can be illustrated and the available structure–activity relationship of these analgesic agents can be explained reasonably. © 2000 John Wiley & Sons, Inc. Int J Quant Chem 78: 285–293, 2000