Showing papers on "Docking (molecular) published in 1995"

Flexible ligand docking using a genetic algorithm.

Abstract: Two computational techniques have been developed to explore the orientational and conformational space of a flexible ligand within an enzyme. Both methods use the Genetic Algorithm (GA) to generate conformationally flexible ligands in conjunction with algorithms from the DOCK suite of programs to characterize the receptor site. The methods are applied to three enzyme-ligand complexes: dihydrofolate reductase-methotrexate, thymidylate synthase-phenolpthalein and HIV protease-thioketal haloperidol. Conformations and orientations close to the crystallographically determined structures are obtained, as well as alternative structures with low energy. The potential for the GA method to screen a database of compounds is also examined. A collection of ligands is evaluated simultaneously, rather than docking the ligands individually into the enzyme.

Protein docking for low-resolution structures

Abstract: A typical problem for a docking procedure is how to match two molecules with known 3-D structure so as to predict the configuration of their complex. A very serious obstacle to docking is an inherent inaccuracy in the 3-D structures of the molecules. In general, existing molecular recognition techniques are not designed for cases where (i) conformational changes upon macromolecular complex formation are substantial or (ii) the X-ray data on one or both (macro) molecules are not available, and the structures, based on alternative sources (NMR, modeling), are not well defined. We designed a direct computer experiment using molecules totally deprived of any structural features smaller than 7 A. This was performed on the basis of a previously developed docking algorithm. The modified procedure was applied to a number of known protein complexes taken from the Brookhaven Protein Data Bank. In most cases, a pronounced trend towards the correct structure of the molecular complex was clearly indicated and the real binding sites were predicted. The distinction between the prediction of the antigen-antibody complex and other molecular pairs may reflect important differences in the principles of complex formation. The results strongly suggest the use of our recognition procedure for docking studies where the detailed structures of the molecules are lacking.

Flexible ligand docking without parameter adjustment across four ligand–receptor complexes

Abstract: Understanding molecular recognition is one of the fundamental problems in molecular biology. Computationally, molecular recognition is formulated as a docking problem. Ideally, a molecular docking algorithm should be computationally efficient, provide reasonably thorough search of conformational space, obtain solutions with reasonable consistency, and not require parameter adjustments. With these goals in mind, we developed DIVALI (Docking wIth eVolutionary AlgorIthms), a program which efficiently and reliably searches for the possible binding modes of a ligand within a fixed receptor. We use an AMBER-type potential function and search for good ligand conformations using a genetic algorithm (GA). We apply our system to study the docking of both rigid and flexible ligands in four different complexes. Our results indicate that it is possible to find diverse binding modes, including structures like the crystal structure, all with comparable potential function values. To achieve this, certain modifications to the standard GA recipe are essential. © 1995 John Wiley & Sons, Inc.

Pseudospecific docking surfaces on electron transfer proteins as illustrated by pseudoazurin, cytochrome c550 and cytochrome cd1 nitrite reductase.

Abstract: Pseudospecific docking surfaces on electron transfer proteins as illustrated by pseudoazurin, cytochrome c550 and cytochrome cd1 nitrite reductase

A computer modeling postulated mechanism for angiotensin II receptor activation.

Abstract: The angiotensin II receptor of the AT1-type has been modeled starting from the experimentally determined three-dimensional structure of bacteriorhodopsin as the template. Intermediate 3D structures of rhodopsin andβ 2-adrenergic receptors were built because no direct sequence alignment is possible between the AT1 receptor and bacteriorhodopsin. Docking calculations were carried out on the complex of the modeled receptor with AII, and the results were used to analyze the binding possibilities of DuP753-type antagonistic non-peptide ligands. We confirm that the positively charged Lys199 on helix 5 is crucial for ligand binding, as in our model; the charged side chain of this amino acid interacts strongly with the C-terminal carboxyl group of peptide agonists or with the acidic group at the 2′-position of the biphenyl moiety of DuP753-type antagonists. Several other receptor residues which are implicated in the binding of ligands and the activation of receptor by agonists are identified, and their functional role is discussed. Therefore, a plausible mechanism of receptor activation is proposed. The three-dimensional docking model integrates most of the available experimental observations and helps to plan pertinent site-directed mutagenesis experiments which in turn may validate or modify the present model and the proposed mechanism of receptor activation.

Docking flexible molecules: A case study of three proteins

Abstract: A genetic algorithm (GA) conformation search method is used to dock a series of flexible molecules into one of three proteins. The proteins examined are thermolysin (tmn), carboxypeptidase A (cpa), and dihydrofolate reductase (dfr). In the latter two proteins, the crystal ligand was redocked. For thermolysin, we docked eight ligands into a protein conformation derived from a single crystal structure. The bound conformations of the other ligands in tmn are known. In the cpa and dfr cases, and in seven of the eight tmn ligands, the GA docking method found conformations within 1.6 A root mean square (rms) of the relaxed crystal conformation. © 1995 John Wiley & Sons, Inc.

Time-efficient docking of flexible ligands into active sites of proteins

Abstract: We present an algorithm for placing flexible molecules in active sites of proteins. The two major goals in the development of our docking program, called FLEXX, are the explicit exploitation of molecular flexibility of the ligand and the development of a model of the docking process that includes the physico-chemical properties of the molecules. The algorithm consists of three phases: The selection of a base fragment, the placement of the base fragment in the active site, and the incremental construction of the ligand inside the active site. Except for the selection of the base fragment, the algorithm runs without manual intervention. The algorithm is tested by reproducing 11 receptor-ligand complexes known from X-ray crystallography. In all cases, the algorithm predicts a placement of the ligand which is similar to the crystal structure (about 1.5 A RMS deviation or less) in a few minutes on a workstation, assuming that the receptor is given in the bound conformation.

An automated computer vision and robotics-based technique for 3-D flexible biomolecular docking and matching.

Abstract: The generation of binding modes between two molecules, also known as molecular docking, is a key problem in rational drug design and biomolecular recognition. Docking a ligand, e.g., a drug molecule or a protein molecule, to a protein receptor, involves recognition of molecular surfaces as molecules interact at their surface. Recent studies report that the activity of many molecules induces conformational transitions by 'hinge-bending', which involves movements of relatively rigid parts with respect to each other. In ligand-receptor binding, relative rotational movements of molecular substructures about their common hinges have been observed. For automatically predicting flexible molecular interactions, we adapt a new technique developed in Computer Vision and Robotics for the efficient recognition of partially occluded articulated objects. These type of objects consist of rigid parts which are connected by rotary joints (hinges). Our approach is based on an extension and generalization of the Geometric Hashing and Generalized Hough Transform paradigm for rigid object recognition. Unlike other techniques which match each part individually, our approach exploits forcefully and efficiently enough the fact that the different rigid parts do belong to the same flexible molecule. We show experimental results obtained by an implementation of the algorithm for rigid and flexible docking. While the 'correct', crystal-bound complex is obtained with a small RMSD, additional, predictive 'high scoring' binding modes are generated as well. The diverse applications and implications of this general, powerful tool are discussed.

Solution structure of a band 3 peptide inhibitor bound to aldolase: a proposed mechanism for regulating binding by tyrosine phosphorylation.

Abstract: Human erythrocyte band 3 inhibits glycolytic enzymes, including aldolase, by binding these cytoplasmic enzymes at its N-terminus. Phosphorylation of Y8 disrupts inhibition, and there is evidence that in vivo glycolysis levels in erythrocytes are regulated in part by a phosphorylation/dephosphorylation signaling pathway. The structural basis for control by phosphorylation has been investigated by NMR studies on a complex between aldolase and a synthetic peptide corresponding to the first 15 residues of band 3 (MEELQDDYEDMMEEN-NH2). The structure of this band 3 peptide (B3P) when it is bound to rabbit muscle aldolase was determined using the exchange-transferred nuclear Overhauser effect (ETNOE). Two hundred NMR structures for B3P were generated by simulated annealing molecular dynamics with NMR-derived distance restraints and excluding electrostatic terms. Twenty structures were further refined against a force field including full partial charges. The important conformational feature of B3P in the bound state is a folded loop structure involving residues 4-9 and M12 that surrounds Y8 and is stabilized by a hydrophobic cluster with the ring of Y8 sandwiched between the methyl groups of L4 and M12. Differential line broadening indicates that this loop structure binds aldolase in a relatively specific manner, while terminal regions are structurally heterogeneous. To better understand B3P inhibition of aldolase and the mechanism of phosphorylation control, a complex was modeled by docking B3P into the active site of aldolase and optimizing the fit using restrained molecular dynamics and energy minimization. The B3P loop is complementary in conformation to the beta-barrel central core containing the aldolase active site residues. Binding is electrostatic in nature with numerous ionic and hydrogen-bonding interactions involving several conserved lysine and arginine residues of aldolase. How phosphorylation of band 3 could disrupt inhibition was considered by modeling a phosphoryl moiety onto Y8 of B3P. An energetic analysis with respect to rigid phosphate rotation suggests that aldolase inhibition is reversed primarily because of electrostatic repulsion between B3P residues that destabilizes the B3P loop formed in the complex. This proposed intramolecular mechanism for blocking protein--protein association by electrostatic repulsion with the phosphoryl group may be applicable to other protein--protein signaling complexes.

Monte Carlo algorithms for docking to proteins

Abstract: The goal of docking is to predict binding interactions between molecules. We are primarily interested in docking as a tool for the structure-based design of new ligands that could serve as lead compounds for drug development. The program BOXSEARCH uses a Monte Carlo algorithm to explore the relative orientation and position of two molecules. Multiple runs are carried out from different random starting positions and orientations, and the temperature of the system is gradually reduced. An unbiased sampling of low energy states is the result. BOXSEARCH has been tested on a number of known complexes, involving both protein and small molecule ligands. Although a better treatment of solvent effects and of flexibility would improve the ranking of results, the complexes can be reconstructed successfully, even using uncomplexed conformations of the molecules. We are currently implementing two major enhancements. First, the code is being rewritten in a more general and adaptable form, using the object-orie...

Protein Docking in the Absence of Detailed Molecular Structures

Abstract: The design of new computational procedures to predict molecular complexes is a fast developing area stimulated by the growing demands of researchers working in various fields of molecular biology and looking for more powerful tools for their investigations. The problem for molecular recognition (docking) approaches may be shortly formulated as following: how to match two molecules with known 3D structures in order to predict the configuration of their complex? In the general case, no additional prior knowledge on binding sites is assumed to be available.

Drug design based on receptor modeling using a system “BIOCES [E]”

Abstract: A computer system BIOCES [E] enables medicinal chemists to study the interaction between biologically active molecules and their receptors. It is capable of building receptor protein models, docking drug molecules into the receptor model, and analyzing their interacting modes. By using BIOCES [E], the three-dimensional structures of human and rat cathepsin H were predicted based on the crystal structure of papain. The binding sites of both enzymes are highly conserved and the only amino acid replacement found was residue 61 (papain numbering scheme). A well-known inhibitor of thiol proteases, benzyloxy-carbonyl-L-phenylalanyl-L-alanyl-methylene, was fitted into the binding site by a Monte Carlo method. Steric, electrostatic and hydrophobic aspects of the interactions of the inhibitor and the enzymes were analyzed.