Search

Entropy is an important energetic quantity determining the progression of chemical processes. We propose a new approach to obtain hydration entropy directly from probability density functions in state space. We demonstrate the validity of our approach for a series of cations in aqueous solution. Extensive validation of simulation results was performed. Our approach does not make prior assumptions about the shape of the potential energy landscape and is capable of calculating accurate hydration entropy values. Sampling times in the low nanosecond range are sufficient for the investigated ionic systems. Although the presented strategy is at the moment limited to systems for which a scalar order parameter can be derived, this is not a principal limitation of the method. The strategy presented is applicable to any chemical system where sufficient sampling of conformational space is accessible, for example, by computer simulations.

A large-scale evaluation and comparison of four cavity detection algorithms was carried out. The algorithms SiteFinder, fpocket, PocketFinder, and SiteMap were evaluated on a protein test set containing 5416 protein-ligand complexes and 9900 apo forms, corresponding to a subset of the set used earlier for benchmarking the PocketFinder algorithm. For the holo structures, all four algorithms correctly identified a similar amount of pockets (around 95%). SiteFinder, using optimized parameters, SiteMap, and fpocket showed similar pocket ranking performance, which was defined by ranking the correct binding site on rank 1 of the predictions or within the first 5 ranks of the predictions. On the apo structures, PocketFinder especially and also SiteFinder (optimized parameters) performed best, identifying 96% and 84% of all binding sites, respectively. The fpocket program predicts binding sites most accurately among the algorithms evaluated here. SiteFinder needed an average calculation time of 1.6 s compared with 2 min for SiteMap and around 2 s for fpocket.

The intrinsic conformational preferences and structures of the branched trimannoside, alpha-phenyl 3,6-di-O-(alpha-D-mannopyranosyl)-alpha-D-mannopyranoside (which contains the same carbohydrates found in a key subunit of the core pentasaccharide in N-glycans) and its singly hydrated complex, have been investigated in the gas phase isolated at low temperature in a molecular beam expansion. Conformational assignments of their infrared ion dip spectra, based on comparisons between experiment and ONIOM (B3LYP/6-31+G(d):HF/6-31G(d)) and single-point MP2 calculations have identified their preferred structures and relative energies. The unhydrated trimannoside populates a unique structure supported by two strong, central hydrogen bonds linking the central mannose unit (CM), and its two branches (3M and 6M) closely together, through a cooperative hydrogen-bonding network: OH4(CM)-->OH6(3M)-->OH6(6M). A closely bound structure is also retained in the singly hydrated oligosaccharide, with the water molecule bridging across the 3M and 6M branches to provide additional bonding. This structure contrasts sharply with the more open, entropically favored trimannoside structure determined in aqueous solution at 298 K. In principle this structure can be accessed from the isolated trimannoside structure by a simple conformational change, a twist about the alpha(1,3) glycosidic linkage, increasing the dihedral angle psi[C1(3M)-O3(3M)-C3(CM)-C2(CM)] from approximately 74 degrees to approximately 146 degrees to enable accommodation of a water molecule at the centrally bound site occupied by the hydroxymethyl group on the 3M ring and mediation of the water-linked hydrogen-bonded network: OH4(CM) -->OH(W)-->OH6(6M). The creation of a "water pocket" motif localized at the bisecting axis of the trimannoside is strikingly similar to the structure of more complex N-glycans in water, suggesting perhaps a general role for the "bisecting" OH4 group in the central (CM) mannose unit.

This review gives an introduction into ligand - receptor docking and illustrates the basic underlying concepts. An overview of different approaches and algorithms is provided. Although the application of docking and scoring has led to some remarkable successes, there are still some major challenges ahead, which are outlined here as well. Approaches to address some of these challenges and the latest developments in the area are presented. Some aspects of the assessment of docking program performance are discussed. A number of successful applications of structure-based virtual screening are described.

Apoptosis-inducing factor (AIF) is a mitochondrial flavoprotein, which upon apoptosis induction translocates to the nucleus where it interacts with DNA by virtue of positive charges clustered on the AIF surface. Here we show that the AIF interactome, as determined by mass spectroscopy, contains a large panel of ribonucleoproteins, which apparently bind to AIF through the RNA moiety. However, AIF is devoid of any detectable RNAse activity both in vitro and in vivo. Recombinant AIF can directly bind to DNA as well as to RNA. This binding can be visualized by electron microscopy, revealing that AIF can condense DNA, showing a preferential binding to single-stranded over double-stranded DNA. AIF also binds and aggregates single-stranded and structured RNA in vitro. Single-stranded poly A, poly G and poly C, as well double-stranded A/T and G/C RNA competed with DNA for AIF binding with a similar efficiency, thus corroborating a computer-calculated molecular model in which the binding site within AIF is the same for distinct nucleic acid species, without a clear sequence specificity. Among the preferred electron donors and acceptors of AIF, nicotine adenine dinucleotide phosphate (NADP) was particularly efficient in enhancing the generation of higher-order AIF/DNA and AIF/RNA complexes. Altogether, these data support a model in which a direct interaction of AIF contributes to the compaction of nucleic acids within apoptotic cells.

The intrinsic conformer specific vibrational spectra of two important subunits of the core pentasaccharide of N-linked glycans, the alpha(1,3) and alpha(1,6) dimannosides, have been recorded in the gas phase. Coupling these measurements with a computational exploration of their conformational landscapes has enabled their conformational assignment and has identified characteristic vibrational signatures associated with particular conformational families-including those that do or do not display inter-ring hydrogen bonding across the glycosidic linkage. In addition, the IR spectra of the monosaccharide moieties provide benchmarks, through which the conformational assignments can be refined. This introduces a general concept of modularity and secondary structure in oligosaccharides--essential for the success of similar studies on larger oligosaccharides in the future.

Docking programs can generate subsets of a compound collection with an increased percentage of actives against a target (enrichment) by predicting their binding mode (pose) and affinity (score), and retrieving those with the highest scores. Using the QXP and GOLD programs, we compared the ability of six single scoring functions (PLP, Ligscore, Ludi, Jain, ChemScore, PMF) and four composite scoring models (Mean Rank: MR, Rank-by-Vote: Vt, Bayesian Statistics: BS and PLS Discriminant Analysis: DA) to separate compounds that are active against CDK2 from inactives. We determined the enrichment for the entire set of actives (IC50 < 10 microM) and for three activity subsets. In all cases, the enrichment for each subset was lower than for the entire set of actives. QXP outperformed GOLD at pose prediction, but yielded only moderately better enrichments. Five to six scoring functions yielded good enrichments with GOLD poses, while typically only two worked well with QXP poses. For each program, two scoring functions generally performed better than the others (Ligscore2 and Ludi for GOLD; QXP and Jain for QXP). Composite scoring functions yielded better results than single scoring functions. The consensus approaches MR and Vt worked best when separating micromolar inhibitors from inactives. The statistical approaches BS and DA, which require training data, performed best when distinguishing between low and high nanomolar inhibitors. The key observation that all hit rate profiles for all four activity intervals for all scoring schemes for both programs are significantly better than random, is evidence that docking can be successfully applied to enrich compound collections., (Copyright 2005 Wiley-Liss, Inc.)

The conformation of phenyl-substituted monosaccharides (mannose, galactose, and glucose) and their singly hydrated complexes has been investigated in the gas phase by means of a combination of mass selected, conformer specific ultraviolet and infrared double resonance hole burning spectroscopy experiments, and ab initio quantum chemistry calculations. In each case, the water molecule inserts into the carbohydrate at a position where it can replace a weak intramolecular interaction by two stronger intermolecular hydrogen bonds. The insertion can produce significant changes in the conformational preferences of the carbohydrates, and there is a clear preference for structures where cooperative effects enhance the stability of the monosaccharide conformers to which the water molecule chooses to bind. The conclusions drawn from the study of monosaccharide-water complexes are extended to the disaccharide lactose and discussed in the light of the underlying mechanisms that may be involved in the binding of carbohydrate assemblies to proteins and the involvement, or not, of key structural water molecules.

Inhibitors of HIV protease have been shown to have antiapoptotic effects in vitro, yet whether these effects are seen in vivo remains controversial. In this study, we have evaluated the impact of the HIV protease inhibitor (PI) nelfinavir, boosted with ritonavir, in models of nonviral disease associated with excessive apoptosis. In mice with Fas-induced fatal hepatitis, Staphylococcal enterotoxin B-induced shock, and middle cerebral artery occlusion-induced stroke, we demonstrate that PIs significantly reduce apoptosis and improve histology, function, and/or behavioral recovery in each of these models. Further, we demonstrate that both in vitro and in vivo, PIs block apoptosis through the preservation of mitochondrial integrity and that in vitro PIs act to prevent pore function of the adenine nucleotide translocator (ANT) subunit of the mitochondrial permeability transition pore complex.

The beta(1-->4) glycosidic linkage found in lactose is a prevalent structural motif in many carbohydrates and glycoconjugates. Using UV and IR ion-dip spectroscopies to probe benzyl lactoside isolated in the gas phase, we find that the disaccharide unit adopts only a single, rigid structure. Its fully resolved infrared ion-dip spectrum is in excellent agreement with that of the global minimum structure computed ab initio. This has glycosidic torsion angles of phi(H) (H1-C1-O-C4') approximately 180 degrees and psi(H) (C1-O-C4'-H4') approximately 0 degrees which correspond to a rotation of approximately 150 degrees about the glycosidic bond compared to the accepted solution-phase conformation. We discuss the biological implications of this discovery and the generality of the strategies employed in making it.

Novel scoring functions that predict the affinity of a ligand for its receptor have been developed. They were built with several statistical tools (partial least squares, genetic algorithms, neural networks) and trained on a data set of 100 crystal structures of receptor-ligand complexes, with affinities spanning 10 log units. The new scoring functions contain both descriptors generated by the QXP docking program and new descriptors that were developed in-house. These new descriptors are based on solvent accessible surface areas and account for conformational entropy changes and desolvation effects of both ligand and receptor upon binding. The predictive r(2) values for a test set of 24 complexes are in the 0.712-0.741 range and RMS prediction errors in the 1.09-1.16 log K(d) range. Inclusion of the new descriptors led to significant improvements in affinity prediction, compared to scoring functions based on QXP descriptors alone. However, the QXP descriptors by themselves perform better in binding mode prediction. The performance of the linear models is comparable to that of the neural networks. The new functions perform very well, but they still need to be validated as universal tools for the prediction of binding affinity.

Six docking programs (FlexX, GOLD, ICM, LigandFit, the Northwestern University version of DOCK, and QXP) were evaluated in terms of their ability to reproduce experimentally observed binding modes (poses) of small-molecule ligands to macromolecular targets. The accuracy of a pose was assessed in two ways: First, the RMS deviation of the predicted pose from the crystal structure was calculated. Second, the predicted pose was compared to the experimentally observed one regarding the presence of key interactions with the protein. The latter assessment is referred to as interactions-based accuracy classification (IBAC). In a number of cases significant discrepancies were found between IBAC and RMSD-based classifications. Despite being more subjective, the IBAC proved to be a more meaningful measure of docking accuracy in all these cases.

Cyclophilin A (CypA) was determined to interact with apoptosis-inducing factor (AIF) by mass spectroscopy, coimmunoprecipitation, pull-down assays, and molecular modeling. During the initial, caspase-independent stage of chromatin condensation that accompanies apoptosis, AIF and CypA were found to coimmunolocalize in the nucleus. Recombinant AIF and CypA proteins synergized in vitro in the degradation of plasmid DNA, as well as in the capacity to induce DNA loss in purified nuclei. The apoptogenic cooperation between AIF and CypA did not rely on the CypA peptidyl-prolyl cis-trans isomerase activity. In Cyp-expressing cells, AIF overexpression augmented apoptotic chromatinolysis. The AIF-dependent large-scale DNA fragmentation was less pronounced in CypA knockout cells as compared to controls. AIF mutants lacking the CypA-binding domain were inefficient apoptosis sensitizers in transfection experiments. Moreover, AIF failed to sensitize CypA knockout cells to apoptosis induction, and this defect in the AIF response was reversed by reintroduction of the CypA gene into CypA-deficient cells. In summary, AIF and CypA collaborate in chromatinolysis.

Heat shock protein 70 (HSP70) inhibits apoptosis and thereby increases the survival of cells exposed to a wide range of lethal stimuli. HSP70 has also been shown to increase the tumorigenicity of cancer cells in rodent models. The protective function of this chaperone involves interaction and neutralization of the caspase activator apoptotic protease activation factor-1 and the mitochondrial flavoprotein apoptosis-inducing factor (AIF). In this work, we determined by deletional mutagenesis that a domain of AIF comprised between amino acids 150 and 228 is engaged in a molecular interaction with the substrate-binding domain of HSP70. Computer calculations favored this conclusion. On the basis of this information, we constructed an AIF-derived protein, which is cytosolic, noncytotoxic, yet maintains its capacity to interact with HSP70. This protein, designated ADD70, sensitized different human cancer cells to apoptosis induced by a variety of death stimuli by its capacity to interact with HSP70 and therefore to sequester HSP70. Thus, its chemosensitizing effect was lost in cells in which inducible HSP70 genes had been deleted. These data delineate a novel strategy for the selective neutralization of HSP70.

Investigations into the conformational behaviour of macrocyclic ligands 5 and 6 derived from (R,R)-1,2-diaminocyclohexane have been undertaken using molecular modelling, single crystal X-ray diffraction and variable temperature 1H NMR spectroscopy. These have revealed that the lowest energy conformers in both cases do not possess the expected C2-element of symmetry, which can only be accessed at higher temperatures. Instead both molecules exist as C1-conformers at room temperature and in the solid state. In solution a range of dynamic exchange processes is observed which result, in part from the inherent strain in these fused bicyclic systems. An unexpected but characteristic feature of the C1-symmetric conformers is highlighted by the presence of a signal at unexpectedly low field in their 1H NMR spectra due to the interaction of two of the sulfonyl oxygen atoms with one of the bridgehead hydrogen atoms.

A general introduction to molecular modelling techniques in the area of protein-ligand interactions is given. Methods covered range from binding-site analysis to statistical treatment of sets of ligands. The main focus of this paper is on docking and scoring. After an outline of the main concepts, two specific application examples are given.

Suppressor Of Cytokine Signalling 1 (SOCS-1) is one of the proteins responsible for the negative regulation of the JAK-STAT pathway triggered by many cytokines. This important inhibition involves complex formation between SOCS-1 and JAK2, which requires particular structural domains (KIR, ESS and SH2) on SOCS-1. A three-dimensional theoretical model of SOCS-1 is presented here. The model was generated by the application of different modelling techniques, including threading, structure-based modelling, surface analysis and protein docking. The structure accounts for the interactions between SOCS-1 and two other key proteins in the JAK-STAT pathway, namely JAK2 and Elongin BC. The proposed model for the interaction between SOCS-1 and JAK2 suggests that the SOCS-1 suppress the kinase activity of JAK2 by obstructing the catalytic groove of the tyrosine kinase. Subsequent interaction of the JAK-SOCS complex with Elongin BC was also modelled. A sequence and structural comparison between the SH2 domain of SOCS-1 and the SH2 domains of other proteins highlights key residues that could be responsible for SOCS-1 specificity. Currently available mutational data are evaluated. The results are consistent with the experimental data and they provide deeper insights into the inhibitory function of SOCS-1 at a molecular level.

Aiming at new drugs to efficiently treat diseases, in which either increased or decreased levels of active vitamin D are desirable, we have designed some 400 structurally different azole-type inhibitors and examined their capacity to selectively block vitamin D metabolism by CYP24 or synthesis by CYP27B, in human keratinocytes. Based on resulting data, we built pharmacophore models of the active sites using commercial software. The overlay of potent selective compounds indicated similar docking modes in the two-substrate pockets and allowed for identification of bioactive conformations. Superimposing these bioactive conformations with low energy conformers of 25(OH)D(3) suggested that the substrate-mimicked by strong inhibitors in size, shape and lipophilic character-binds to both enzymes in 6s-trans configuration. Pharmacophoric models implied a similar geometry of the substrate sites, nevertheless specific features of CYP24 and CYP27B could be defined. Bulky substituents in alpha-position to the azole caused selectivity for CYP24, whereas bulky substituents in beta-position could result in selectivity for CYP27B. Moreover, studies with small sterically restricted inhibitors revealed a probable location of the 3-OH-group of 25(OH)D(3) in CYP27B. In the absence of crystal structures, our inhibitors are valuable tools to model and understand the active sites of vitamin D hydroxylases, resulting in the design of powerful, selective therapeutics., (Copyright 2002 Wiley-Liss, Inc.)

A theoretical model of human Janus kinase 2 (JAK2) comprising all seven Janus homology domains is presented. The model was generated by application of homology modelling approaches. The three-dimensional structure contains, starting from the N-terminus, FERM (4.1, ezrin, radixin, moesin), SH2 (Src homology region 2), tyrosine kinase-like, and tyrosine kinase domains. The predicted inter-domain orientation in JAK2 is discussed and the currently existing mutational data for Janus kinases are evaluated. Structural details of the SH2 and the FERM domains are presented. The predictions indicate that the SH2 domain is not fully functional. A number of hydrophobic amino acids of the FERM domain that are predicted to be involved in the constitutive association with the cytokine receptors are highlighted. The model gives new insights into the structure-function relationship of this important protein, and areas that could be investigated by mutation studies are highlighted.

The structure of human Janus kinase 2 (JAK2) comprising the two C-terminal domains (JH1 and JH2) was predicted by application of homology modelling techniques. JH1 and JH2 represent the tyrosine kinase and tyrosine kinase-like domains, respectively, and are crucial for function and regulation of the protein. A comparison between the structures of the two domains is made and structural differences are highlighted. Prediction of the relative orientation of JH1 and JH2 was aided by a newly developed method for the detection of correlated amino acid mutations. Analysis of the interactions between the two domains led to a model for the regulatory effect of JH2 on JH1. The predictions are consistent with available experimental data on JAK2 or related proteins and provide an explanation for inhibition of JH1 tyrosine kinase activity by the adjacent JH2 domain.

Blockade of the N-methyl-D-aspartate receptor by uncompetitive antagonists has implications for symptomatic and neuroprotective therapy of various neuropsychiatric diseases. Since the three-dimensional (3D) structure of this ion channel is unknown, the structural requirements for uncompetitive inhibition were investigated by application of a three-step strategy: At first, Ki values were measured for a number of structurally diverse compounds at the phencyclidine (PCP) binding site in postmortem human frontal cortex. Second, a pharmacophore model was developed for this set of molecules employing a mathematical method called graph theory. The resulting pharmacophore provided a very good explanation for the ability of structurally diverse compounds to bind to the same binding site. Using the experimental data and the pharmacophore as a basis for the third step, a three-dimensional quantitative structure-activity relationship (3D-QSAR) applying comparative molecular field analysis (CoMFA) was performed. The QSAR proved to be highly consistent and showed good predictiveness for several additional molecules. The results give a deeper insight into the structural requirements for effective NMDA receptor antagonism and offer the opportunity for improved drug design. The study represents the first quantitative 3D-QSAR model for NMDA receptor blockade, and it comprises structurally very different molecules. That the alignment for a highly consistent CoMFA is based on an automated 3D pharmacophore analysis has important methodological implications.

Thrombopoietin (TPO) is a glycoprotein hormone that regulates red blood cell production. Presented here is a modeling study of the extracellular region of the human thrombopoietin receptor complex, in particular the TPO-receptor interface. The models were developed from structural homology to other cytokines and their receptors. Experimental evidence suggests that the receptor is homodimeric and it was modeled accordingly. Key interactions are shown that correlate with previous cytokine receptor complexes, and the pattern of cysteine bonding (Cys7-Cys151 and Cys29-Cys85) agrees with that experimentally determined for thrombopoietin. These models pave the way for possible mutagenesis experimentation and the design of (ant)agonists.

Following a recent model of human interleukin-7 (IL-7), we present here a modeling study of the extracellular part of the human IL-7 receptor complex, including the IL-7 specific (IL-7R) and the common gamma (gamma c) chains. The investigation is based on structural homology to the complex of human growth hormone (hGH) bound to its receptor (hGHR). For domain 1 of IL-7R two different models are presented which differ in the alignment to hGHR in three regions. However, these differences affect binding to IL-7 in only one region, at the interface between loop EF of domain 1 of IL-7R and helix C of IL-7. The disulfide pattern in domain 1 of IL-7R is predicted to deviate from that observed in hGHR in that the C'-E disulfide (hGHR) is replaced by a C-C' cross-link. The prediction for the gamma c chain is compared with two previous studies. The models of the complex provide insight into the binding of IL-7 to its receptor and have implications for the suggestion of mutagenesis experiments and the design of (ant)agonists.

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 molecular field analysis (CoMFA) is a three-dimensional quantitative structure-activity relationship (3D-QSAR) method which correlates precalculated fields surrounding a set of molecules with some target property. Among others, the electrostatic fields are commonly used. These are usually generated by calculating the Coulomb potential between a probe and the molecules bearing atom-centered point charges. The present study was performed in order to investigate up to which extent different methods to calculate electrostatic potentials can influence the results of a CoMFA. Therefore, a variety of charge calculation methods was applied to a data set consisting of 37 ligands of the benzodiazepine receptor inverse agonist-antagonist active site. These methods included Gasteiger-Marsili, semiempirical (MNDO, AM1, and PM3), and ab initio (HF/STO-3G, HF/3-21G*, and HF/6-31G*) charges. Semiempirical as well as ab initio electron populations were derived both from the Mulliken population analysis (MPA) or from fitting the charges to the molecular electrostatic potential (ESPFIT charges). In addition, the molecular electrostatic potentials (MEPs) resulting from ab initio calculations were mapped directly onto the CoMFA grid. With regard to the cross-validated r(2) values (r(2) cv ) of the resulting QSAR models, the ESPFIT-derived potentials yielded generally higher r(2) cv values than those resulting from MPA charges. For example, at the HF/3-21G* level the r(2) cv rose from 0.61 (MPA-derived potentials) to 0.76 (ESPFIT fields). The MEPs mapped directly onto the CoMFA grid were not superior to the corresponding ESPFIT-derived potentials. Semiempirical ESPFIT charges appeared to be of similar quality compared with ab initio ESPFIT electron populations in the CoMFAs. When no scaling between the steric and electrostatic descriptor matrices was applied, the electrostatic contributions were influenced to a high degree by the magnitude of the corresponding field values. © 1996 by John Wiley & Sons, Inc., (Copyright © 1996 John Wiley & Sons, Inc.)

The three-dimensional structure of human interleukin (IL)-7 has been predicted based on homology to human IL-2, IL-4, granulocyte-macrophage colony stimulating factor and growth hormone. The model has a topology common to other cytokines and displays a unique disulfide pattern. Knowledge of the tertiary structure of IL-7 has implications for analysis of key binding regions, suggestions for mutagenesis experiments and design of (ant)agonists. In this context, the model is discussed and compared with other cytokine structures.

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.

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).

A set of 100 novel 2-heterosubstituted statine derivatives inhibiting human immunodeficiency virus type-1 proteinase has been investigated by comparative molecular field analysis. In order to combine the structural information available from X-ray analyses with a predictive quantitative structure-activity relationship (QSAR) model, docking experiments of a prototype compound into the receptor were performed, and the 'active conformation' was determined. The structure of the receptor was taken from the published X-ray analysis of the proteinase with bound MVT-101, the latter compound exhibiting high structural similarity with the inhibitors investigated. The validity of the resulting QSARs was confirmed in four different ways. (1) The common parameters, namely, the cross-validated r2 values obtained by the leave-one-out (LOO) method (r2cv = 0.572-0.593), and (2) the accurate prediction of a test set of 67 compounds (q2 = 0.552-0.569) indicated a high consistency of the models. (3) Repeated analyses with two randomly selected cross-validation groups were performed and the cross-validated r2 values monitored. The resulting average r2 values were of similar magnitudes compared to those obtained by the LOO method. (4) The coefficient fields were compared with the steric and electrostatic properties of the receptor and showed a high level of compatibility. Further analysis of the results led to the design of a novel class of highly active compounds containing an additional linkage between P1' and P3'. The predicted activities of these inhibitors were also in good agreement with the experimentally determined values.

A new automated procedure to improve the predictive quality of CoMFA models for both training and test sets is described. A model of greater consistency is generated by performing small reorientations of the underlying molecules for which too low activities are calculated. In order to predict activities of test compounds, the most similar molecules in the previously optimized model are identified and used as a basis for the prediction. This method has been applied to two independent sets of dihydrofolate reductase inhibitors (80 compounds each, serving as training sets), resulting in a significant increase of the cross-validated r2 value. For both models, the predictive r2 value for a test set consisting of 70 compounds was improved substantially.

The steric descriptors commonly used in CoMFA--Lennard-Jones 6-12 potential-derived interaction energies calculated between a probe atom and the molecules under investigation--have been replaced by variables indicating the presence of an atom of a particular molecule in predefined volume elements (cubes) within the region enclosing the ensemble of superimposed molecules. The resulting 'atom indicator vectors' were used as steric fields in the subsequent PLS analyses, with and without inclusion of electrostatic Coulomb interaction-derived fields. Application of this method to five training sets (80 compounds each) and five test sets (60 compounds each), randomly selected from an ensemble of 256 dihydrofolate reductase inhibitors, leads to models of significantly higher consistency, as indicated by the cross-validated r2 values for the training sets and the predictive r2 values for the test sets.