Strategies for the design of inhibitors of aldose reductase, an enzyme showing pronounced induced-fit adaptations.

Aldose reductase is involved in the polyol pathway, catalyzing the reduction of glucose to sorbitol. However, due to pronounced binding site adaptations, the enzyme can operate on a broad palette of structurally diverse substrates ranging from small aliphatic and aromatic aldehydes up to steroid-type ligands. A comparative analysis of the presently accessible crystal structures of aldose reductase complexes reveals four binding-competent protein conformations. Additional relevant conformers are detected through molecular dynamics simulations. They indicate an equilibrium of several conformers which is shifted towards the binding-competent geometries upon ligand binding. Such a manifold system with several alternative binding site conformers requires some tailored concepts in virtual screening. We followed two strategies, both successfully suggesting new micromolar inhibitors. In a first attempt, we concentrated on one preferred conformer and performed a virtual screening, assuming that the binding pocket of aldose reductase adopts only this conformation. In a second approach, we followed a ligand superpositioning method. Ligands were extracted in their bound conformations from three different crystal structures, all accommodating the ligands with different active site conformations. After merging these ligands into one supermolecule, mutual alignments were computed, taking candidate ligands from a screening database. The latter strategy also retrieved several structurally new inhibitors of micromolar potency.

Ordered water and ligand mobility in the HIV-1 integrase-5CITEP complex: a molecular dynamics study.

A 2 ns molecular dynamics simulation has been carried out for the HIV-1 integrase-5CITEP complex in order to understand the role of water in defining the ligand's binding mode and to address issues of binding site flexibility and ligand motion. Although the ligand retains considerable mobility within the active site, a structural water molecule bridging 5CITEP with Asp 64 and Asn 155 is identified in the simulation. Consideration of this water molecule could open a route to new HIV-1 integrase inhibitors.

Active site binding modes of HIV-1 integrase inhibitors.

Using the crystal structure of the first complex of the HIV-1 integrase catalytic core domain with an inhibitor bound to the active site, structural models for the interaction of various inhibitors with integrase were generated by computational docking. For the compound of the crystallographic study, binding modes unaffected by crystal packing have recently been proposed. Although a large search region was used for the docking simulations, the ligands investigated here are found to bind preferably in similar ways close to the active site. The binding site is formed by residues 64-67, 116, 148, 151-152, 155-156, and 159, as well as by residue 92 in case of the largest ligand of the series. The coherent picture of possible interactions of small-molecule inhibitors at the active site provides an improved basis for structure-based ligand design. The recurring motif of tight interaction with the two lysine residues 156 and 159 is suggested to be of prime importance.

Elbow flexibility and ligand-induced domain rearrangements in antibody Fab NC6.8: large effects of a small hapten.

Four 700-ps molecular dynamics simulations were carried out to analyze the structural dynamics of the antigen-binding antibody fragment NC6.8, which is known to exhibit large structural changes upon complexation. The first simulation was started from the x-ray structure of the uncomplexed Fab and produced trajectory averages that closely match the crystallographic results. It allowed assessment of the flexibility of the Fab, revealing an elbow motion of the variable domains with respect to the constant domains. The second simulation was started from the uncomplexed x-ray structure after insertion of the ligand into the binding site. This perturbation resulted in a significantly altered trajectory, with quaternary structural changes corresponding in many aspects to the experimental differences between complexed and uncomplexed state. The observed trend toward a smaller elbow angle and a higher flexion of the H-chain could also be seen in the third simulation, which was started from the x-ray structure of the complex. The changes were revealed to be a clear consequence of the complexation with the ligand because in the fourth simulation (started from the experimental complex structure after removal of the hapten) the Fab remained close to its initial structure. Analyses of the quaternary structure and the binding site of Fab NC6.8 are presented for all four simulations, and possible interpretations are discussed.

Automated docking of ligands to antibodies: methods and applications.

Many approaches to studying protein-ligand interactions by computational docking are currently available. Given the structures of a protein and a ligand, the ultimate goal of all docking methods is to predict the structure of the resulting complex. This requires a suitable representation of molecular structures and properties, search algorithms to efficiently scan the configuration space for favorable interaction geometries, and accurate scoring functions to evaluate and rank the generated orientations. For many of the available methods, tests on experimentally known antibody-antigen or antibody-hapten complexes have appeared in the literature. In addition, some of them have been used in predictive studies on antibody-ligand interactions to provide structural insights where adequate experimental information is missing. The AutoDock program is presented as example of a method for flexibly docking ligands to antibodies. Applying parameters of the second-generation AMBER force field, three antibody-hapten complexes (AN02, DB3, NC6.8) are used as new test cases to analyze the ability of the method to reproduce experimental findings. The X-ray structures could be reconstituted and the corresponding solutions were ranked with best energy score in all cases. Docking to the free instead of the complexed NC6.8 structure indicated the limits of the rigid protein treatment, although fairly good guesses about the location of the binding site and the contact residues could still be obtained if conformational flexibility was allowed at least in the ligand.

Application of multivariate data analysis methods to comparative molecular field analysis (CoMFA) data: proton affinities and pKa prediction for nucleic acids components.

Multivariate data analysis methods (Principal Component Analysis (PCA) and Partial Least Squares (PLS)) are applied to the analysis of the CoMFA (Comparative Molecular Field Analysis) data for several nucleic acids components. The data set includes nitrogenated bases, nucleosides, linear nucleotides, 3', 5'-cyclic nucleotides and oligonucleotides. PCA is applied to study the structure of the CoMFA data and to detect possible outliers in the data set. PLS is applied to correlate the CoMFA data with either calculated AM1 proton affinities or with experimental pKa values. The possibility of making a prediction of pKa values directly from 3D structures of the monomers for polynucleotides is also shown. The influence of the superposition criteria and of conformational changes along the glycosidic bond on the pKa prediction are studied as well.

Ligand binding by antibody IgE Lb4: assessment of binding site preferences using microcalorimetry, docking, and free energy simulations.

Antibody IgE Lb4 interacts favorably with a large number of different compounds. To improve the current understanding of the structural basis of this vast cross-reactivity, the binding of three dinitrophenyl (DNP) amino acids (DNP-alanine, DNP-glycine, and DNP-serine) is investigated in detail by means of docking and molecular dynamics free energy simulations. Experimental binding energies obtained by isothermal titration microcalorimetry are used to judge the results of the computational studies. For all three ligands, the docking procedure proposes two plausible subsites within the binding region formed by the antibody CDR loops. By subsequent molecular dynamics simulations and calculations of relative free energies of binding, one of these subsites, a tyrosine-surrounded pocket, is revealed as the preferred point of complexation. For this subsite, results consistent with experimental observations are obtained; DNP-glycine is found to bind better than DNP-serine, and this, in turn, is found to bind better than DNP-alanine. The suggested binding mode makes it possible to explain both the moderate binding affinity and the differences in binding energy among the three ligands.

Are prions a relic of an early stage of peptide evolution?

The rather unique properties of prions and their presence in very different kinds of living species suggest that this type of molecule was created at a very early stage of evolution and may even represent a relic from a time where peptide evolution was ongoing and RNA/DNA did not yet exist. A comparison of the most frequently occurring amino acid sequences in known prions with the sequences preferentially formed in the salt-induced peptide formation reaction, the most simple mechanism enabling the formation of peptides under primitive earth conditions, shows a remarkable coincidence that strongly supports this hypothesis.

Ligand-induced domain movement in an antibody Fab: molecular dynamics studies confirm the unique domain movement observed experimentally for Fab NC6.8 upon complexation and reveal its segmental flexibility.

Two molecular dynamics simulations were carried out for the antibody Fab NC6.8, both with and without the guanidinium sweetener ligand NC174, in order to assess the segmental flexibility as well as the conformational changes upon ligand binding. Trajectory analyses of the simulation of the uncomplexed Fab suggest low-amplitude motions of the Ig domains with respect to each other, most clearly reflected by a periodic alteration of the elbow angle within a range of 11 degrees. Upon insertion of the hapten into the binding site, the quaternary structure of the Fab exhibits considerable rearrangements: the elbow angle changes by almost 30 degrees, the light chain is elongated and the heavy chain becomes more flexed. Comparison with experiment reveals some interesting agreements with X-ray crystallographic results published previously.

Comparative molecular field analysis of haptens docked to the multispecific antibody IgE(Lb4)

Using comparative molecular field analysis (CoMFA), three-dimensional quantitative structure-activity relationships were developed for 27 haptens which bind to the monoclonal antibody IgE(Lb4). In order to obtain an alignment for these structurally very diverse antigens, the compounds were docked to a previously modeled receptor structure using the automated docking program AUTODOCK (Goodsell, D.S.; Olson, A.J. Proteins: Struct., Funct., Genet. 1990, 8, 195-202). Remarkably, this alignment method yielded highly consistent QSAR models, as indicated by the corresponding cross-validated r2 values (0.809 for a model with carbon as probe atom, 0.773 for a model with hydrogen as probe atom). Conventional alignment failed in providing a basis for self-consistent CoMFAs. Amino acids Tyr H 50, Tyr H 52, and Trp H 95 of the receptor appeared to be of crucial importance for binding of various antigens. These findings are consistent with earlier considerations of aromatic residues being responsible for the multispecificity of certain immunoglobulins.

Comparative docking studies on ligand binding to the multispecific antibodies IgE-La2 and IgE-Lb4.

A large comparative study is presented in which the binding of approximately 30 different ligands to two IgE antibodies (La2 and Lb4) is analyzed by means of an automated-docking procedure based on simulated annealing. The method is able to reproduce experimentally verified binding orientations, as shown by application to the Ig-AN02-hapten complex. The main address of the study is to investigate the concept of antibody multispecificity. Problems and usefulness of docking in this context are discussed. The results indicate reasons for multispecific binding properties and how they can be understood from the topology of the binding site. Though similar in general behaviour, the two antibodies show interesting differences in their binding characteristics. The binding sites of both antibodies are described and the main interacting residues revealed.

Prediction of IgE(Lb4)-ligand complex structures by automated docking.

A mouse monoclonal anti-2,4,6-trinitrophenyl IgE (clone Lb4) was screened with a random set of over 2000 compounds, and several ligands were found to bind with affinities comparable to that of the immunizing hapten (KD in the microM range). An automated docking algorithm was used for the prediction of complex structures formed by 2,4-dinitrophenyl (DNP) and non-DNP ligands in the fragment variable region of IgE(Lb4). All ligands were found to dock in an L-shaped cavity of 15 x 16 x 10 A, surrounded by complementary-determining regions L1, L3, H2 and H3. The ligands were found to occupy the same binding site in different orientations. For rigid ligands the most stable orientation could be predicted with high probability, based on the calculated energy of binding and the occurrence frequencies of identical complexes obtained by repeated simulations. The localization of a flexible ligand (cycrimine-R) was more ambiguous, but it still docked in the same site. The results support a model for heteroligating antibody (Ab) binding sites, where different ligands utilize the total set of available contacts in different combinations. It is suggested that although pseudoenergies calculated by the docking algorithm do not correlate with experimentally measured binding energies, the screening-and-docking procedure can be useful for the mapping of Ab and other receptor binding sites ligating small molecules.

Heteroligation of a mouse monoclonal IgE antibody (La2) with small molecules, analysed by computer-aided automated docking.

A mouse monoclonal anti-TNP IgE antibody (IgE-La2) was screened by a competitive-binding ELISA with a random pool of over 2000 small molecules, mostly drugs, drug derivatives and metabolites. Thirteen of these (naproxene, beta-carboxy-chi-naphthol, oxolinic acid, hymecromone, 8-aminoquinoline, beta-naphthylamine, chi-nitrilo-cinnamic acid, 1,5-diaminonaphthaline, prolonium iodide, diaspirin, 3,4,5-trimethoxy-cinnamic acid, cycrimine, hemimellitic acid) were found to bind as strongly, or stronger, to the antibody as the immunizing hapten. We have used a Monte Carlo search technique for simulated docking of the DNP and non-DNP ligands to a model of the Fv region of IgE(La2). The validity of structural predictions made by the AutoDock program were tested on IgG(ANO2), the three-dimensional structure of which had been obtained previously by X-ray crystallography and 2D-NMR. The rms differences between the experimentally determined and auto-docked complexes in the energetically most favored binding modes were 0.31-0.44 A. Evaluation of structures of IgE(La2)-ligand complexes [including 2,4-dinitrophenol (DNP), 16 DNP amino acids, and the 13 non-DNP ligands listed above] obtained by computer-aided automated docking, suggested the existence of two subsites within an approximately 12 x 18 A2 groove extending between the H and L CDRs. Some of the ligands (DNP-Glu, 8-aminoquinoline, prolonium-I, beta-naphthylamine) were found to bind exclusively to subsite 1, others (DNP-Ala, chi-nitrilo-cinnamic acid, hemimellitic acid, beta-carboxy-chi-naphthol) to subsite 2. The majority of DNP amino acids and other ligands (oxolinic acid, 3,4,5-trimethoxy-cinnamic acid, diaspirin, [R]-cycrimine) were found to occupy an overlapping area including subsites 1 and 2, while some of the compounds (DNP-Asn, DNP-Pro, hymecromone, 1,5-naphthylenediamine) were predicted to interact with either of these subsites with comparable probabilities. When all of the docked La2-ligand complexes were taken into account, five tyrosine residues (H33, L32, L91, L92, L96) were found to provide the majority (53.4%) of all observed contact points. Thus, a multitude of interactions with aromatic residues, and a combinatorial type of interaction within the binding region, seem to be the major factors to explain the mechanism of heteroligation by IgE(La2).