Search

The DNA-repair enzyme MutT homolog 1 (MTH1) is a potential target for a broad range of tumors. Its substrate binding site features a non-catalytical pair of aspartic acids which resembles the catalytic dyad of aspartic proteases. We hypothesized that inhibitors of the latter might be re-targeted for MTH1 despite the two enzyme classes having different substrates and catalyze different reactions. We selected from the crystal structures of holo aspartic proteases a library of nearly 350 inhibitors for in silico screening. Three fragment hits were identified by docking and scoring according to a force field-based energy with continuum dielectric solvation. These fragments showed good ligand efficiency in a colorimetric assay (MW 300 Da and IC50 50 μM). Molecular dynamics simulations were carried out for determining the most favorable interaction patterns. On the basis of the simulation results we evaluated in vitro seven commercially available compounds, two of which showed submicromolar potency for MTH1. To obtain definitive evidence of the predicted binding modes we solved the crystal structures of five of the 10 inhibitors predicted in silico. The final step of hit optimization was guided by protein crystallography and involved the synthesis of a single compound, the lead 11, which shows nanomolar affinity for MTH1 in two orthogonal binding assays, and selectivity higher than 2000-fold against its original target (BACE1). The high rate of fragment-hit identification and the fast optimization suggest that ligand retargeting by binding site analogy is an efficient strategy for drug design.

Eph receptor tyrosine kinases, particularly EphA2 and EphB4, represent promising candidates for molecular imaging due to their essential role in cancer progression and therapy resistance. Xanthine derivatives were identified to be potent Eph receptor inhibitors with IC50 values in the low nanomolar range (1-40 nm). These compounds occupy the hydrophobic pocket of the ATP-binding site in the kinase domain. Based on lead compound 1, we designed two fluorine-18-labelled receptor tyrosine kinase inhibitors ([18F]2/3) as potential tracers for positron emission tomography (PET). Docking into the ATP-binding site allowed us to find the best position for radiolabelling. The replacement of the methyl group at the uracil residue ([18F]3) rather than the methyl group of the phenoxy moiety ([18F]2) by a fluoropropyl group was predicted to preserve the affinity of the lead compound 1. Herein, we point out a synthesis route to [18F]2 and [18F]3 and the respective tosylate precursors as well as a labelling procedure to insert fluorine-18. After radiolabelling, both radiotracers were obtained in approximately 5% radiochemical yield with high radiochemical purity (98%) and a molar activity of10 GBq μmol-1. In line with the docking studies, first cell experiments revealed specific, time-dependent binding and uptake of [18F]3 to EphA2 and EphB4-overexpressing A375 human melanoma cells, whereas [18F]2 did not accumulate at these cells. Since both tracers [18F]3 and [18F]2 are stable in rat blood, the novel radiotracers might be suitable for in vivo molecular imaging of Eph receptors with PET.

We report a crystallographic analysis of small-molecule ligands of the human YTHDC1 domain that recognizes N6-methylated adenine (m

Recent years have witnessed rapid developments of computer-aided drug design methods, which have reached accuracy that allows their routine practical applications in drug discovery campaigns. Protein structure-based methods are useful for the prediction of binding modes of small molecules and their relative affinity. The high-throughput docking of up to 106 small molecules followed by scoring based on implicit-solvent force field can robustly identify micromolar binders using a rigid protein target. Molecular dynamics with explicit solvent is a low-throughput technique for the characterization of flexible binding sites and accurate evaluation of binding pathways, kinetics, and thermodynamics. In this review we highlight recent advancements in applications of ligand docking tools and molecular dynamics simulations to ligand identification and optimization.

T-cell receptor (TCR) recognition of the myelin basic protein (MBP) peptide presented by major histocompatibility complex (MHC) protein HLA-DR2a, one of the MHC class II alleles associated with multiple sclerosis, is highly variable. Interactions in the trimolecular complex between the TCR of the MBP83-99-specific T cell clone 3A6 with the MBP-peptide/HLA-DR2a (abbreviated TCR/pMHC) lead to substantially different proliferative responses when comparing the wild-type decapeptide MBP90-99 and a superagonist peptide, which differs mainly in the residues that point toward the TCR. Here, we investigate the influence of the peptide sequence on the interface and intrinsic plasticity of the TCR/pMHC trimolecular and pMHC bimolecular complexes by molecular dynamics simulations. The intermolecular contacts at the TCR/pMHC interface are similar for the complexes with the superagonist and the MBP self-peptide. The orientation angle between TCR and pMHC fluctuates less in the complex with the superagonist peptide. Thus, the higher structural stability of the TCR/pMHC tripartite complex with the superagonist peptide, rather than a major difference in binding mode with respect to the self-peptide, seems to be responsible for the stronger proliferative response.

The high-throughput docking protocol called ALTA-VS (anchor-based library tailoring approach for virtual screening) was developed in 2005 for the efficient in silico screening of large libraries of compounds by preselection of only those molecules that have optimal fragments (anchors) for the protein target. Here we present an updated version of ALTA-VS with a broader range of potential applications. The evaluation of binding energy makes use of a classical force field with implicit solvent in the continuum dielectric approximation. In about 2 days per protein target on a 96-core compute cluster (equipped with Xeon E3-1280 quad core processors at 2.5 GHz), the screening of a library of nearly 77 000 diverse molecules with the updated ALTA-VS protocol has resulted in the identification of 19, 3, 3, and 2 μM inhibitors of the human bromodomains ATAD2, BAZ2B, BRD4(1), and CREBBP, respectively. The success ratio (i.e., number of actives in a competition binding assay in vitro divided by the number of compounds tested) ranges from 8% to 13% in dose-response measurements. The poses predicted by fragment-based docking for the three ligands of the BAZ2B bromodomain were confirmed by protein X-ray crystallography.

Bromodomains are protein modules that bind to acetylated lysine side chains in histones and other proteins. The bromodomain adjacent to zinc finger domain protein 2B (BAZ2B) has been reported to be poorly druggable. Here, we screened an in-house library of 350 fragments by automatic docking to the BAZ2B bromodomain. The top 12 fragments according to the predicted binding energy were selected for experiments of soaking into apo crystals of BAZ2B which yielded the structure of the complex for four of them, which is a hit rate of 33%. Additional crystal structures were solved for BAZ2B and two scaffolds identified by analogy. For three topologically similar fragments, the crystal structures reveal binding modes with different penetration, i.e., with zero, one, and two water molecules, respectively, located between the fragment and the side chain of a conserved tyrosine (Tyr1901) in the bottom of the acetyl lysine pocket of BAZ2B. Furthermore, a remarkable stereoselectivity of the acetyl lysine pocket emerges from the crystal structures of the bromodomains of BAZ2B and SMARCA4 in complex with the chiral diol MPD (2-methyl-2,4-pentanediol).

BRPF1 plays a scaffolding role in transcription. We report on fragment screening by high-throughput docking to the BRPF1 bromodomain which resulted in six chemotypes with very favorable ligand efficiency (0.45–0.50 kcal/mol per non-hydrogen atom). Twenty crystal structures of BRPF1/ligand complexes show structural conservation in the acetyllysine binding site, common binding motifs, and unusual interactions (e.g., the replacement of a conserved water molecule). The structural information is useful for the design of chemical probes.

Small-molecule hits for the bromodomains of CREBBP and BAZ2B have been identified by scaffold hopping followed by docking of a set of ∼200 compounds containing the acetyl indole scaffold. Chemical synthesis of nearly 30 derivatives has resulted in ligands of representatives of three subfamilies of human bromodomains with favorable ligand efficiency. The X-ray crystal structures of three different bromodomains (CREBBP, BAZ2B, and BRPF1b) in complex with acetyl indole derivatives reveal the influence of the gatekeeper residue on the orientation of small-molecule ligands in the acetyl lysine binding site.

This is the peer-reviewed version of the following article: Peón, A., Robles, A., Blanco, B., Convertino, M., Thompson, P., & Hawkins, A. et al. (2017). Reducing the Flexibility of Type II Dehydroquinase for Inhibition: A Fragment-Based Approach and Molecular Dynamics Study. Chemmedchem, 12(18), 1512-1524, which has been published in final form at https://doi.org/10.1002/cmdc.201700396. This article may be used for non-commercial purposes in accordance with Wiley-VCH Terms and Conditions for Self-Archiving A multidisciplinary approach was used to identify and optimize a quinazolinedione‐based ligand that would decrease the flexibility of the substrate‐covering loop (catalytic loop) of the type II dehydroquinase from Helicobacter pylori. This enzyme, which is essential for the survival of this bacterium, is involved in the biosynthesis of aromatic amino acids. A computer‐aided fragment‐based protocol (ALTA) was first used to identify the aromatic fragments able to block the interface pocket that separates two neighboring enzyme subunits and is located at the active site entrance. Chemical modification of its non‐aromatic moiety through an olefin cross‐metathesis and Seebach's self‐reproduction of chirality synthetic principle allowed the development of a quinazolinedione derivative that disables the catalytic loop plasticity, which is essential for the enzyme′s catalytic cycle. Molecular dynamics simulations revealed that the ligand would force the catalytic loop into an inappropriate arrangement for catalysis by strong interactions with the catalytic tyrosine and by expelling the essential arginine out of the active site Secretaría de Estado de Investigación, Desarrollo e Innovación. Grant Numbers: SAF2013-42899-R, SAF2016-75638-R Consellería de Cultura, Educación e Ordenación Universitaria, Xunta de Galicia. Grant Number: ED431G/09 European Regional Development Fund Ministerio de Educación, Cultura y Deporte SI

Bromodomains are α-helical bundles of approximately 110 residues that recognize acetylated lysine side chains mainly on histone tails. Bromodomains are known to play an important role in cancer and inflammation, and as such, significant efforts are being made to identify small-molecule inhibitors of these epigenetic reader proteins. Here, explicit solvent molecular dynamics (MD) simulations of two bromodomains (BAZ2B and CREBBP) are used to analyze the water molecules that seem to be conserved at the bottom of the acetyl-lysine binding site in most crystal structures of bromodomains. The MD runs suggest that the occupancy of the structured water molecules is influenced by conformational transitions of the loop that connects helices Z and A. Additional simulations in the presence of 50 molecules of cosolvent (i.e., 440 mM of dimethylsulfoxide, methanol, or ethanol) indicate that some of the structured water molecules can be displaced transiently. The residence time in the acetyl-lysine binding site is calculated to be about 1 ns, 2-5 ns, and 10-30 ns for methanol, ethanol, and dimethylsulfoxide, respectively, while the affinity of the three cosolvents for BAZ2B and CREBBP is in the range of 50-500 mM. The results described have implications for ligand design, suggesting that only structured water molecules that do not exchange with cosolvent should be maintained in crystal structures used for docking campaigns, and that hydroxy substituents should be incorporated in the ligand so as to map the structured water molecules replaced by (m)ethanol.

A detailed understanding of the conformational dynamics of biological molecules is difficult to obtain by experimental techniques due to resolution limitations in both time and space. Computer simulations avoid these in theory but are often too short to sample rare events reliably. Here we show that the progress index-guided sampling (PIGS) protocol can be used to enhance the sampling of rare events in selected parts of biomolecules without perturbing the remainder of the system. The method is very easy to use as it only requires as essential input a set of several features representing the parts of interest sufficiently. In this feature space, new states are discovered by spontaneous fluctuations alone and in unsupervised fashion. Because there are no energetic biases acting on phase space variables or projections thereof, the trajectories PIGS generates can be analyzed directly in the framework of transition networks. We demonstrate the possibility and usefulness of such focused explorations of biomolecules with two loops that are part of the binding sites of bromodomains, a family of epigenetic "reader" modules. This real-life application uncovers states that are structurally and kinetically far away from the initial crystallographic structures and are also metastable. Representative conformations are intended to be used in future high-throughput virtual screening campaigns.

Multiple molecular dynamics simulations with explicit solvent at room temperature and at 400 K were carried out to characterize designed ankyrin repeat (AR) proteins with full-consensus repeats. Using proteins with one to five repeats, the stability of the native structure was found to increase with the number of repeats. The C-terminal capping repeat, originating from the natural guanine-adenine-binding protein, was observed to denature first in almost all high-temperature simulations. Notably, a stable intermediate is found in experimental equilibrium unfolding studies of one of the simulated consensus proteins. On the basis of simulation results, this intermediate is interpreted to represent a conformation with a denatured C-terminal repeat. To validate this interpretation, constructs without C-terminal capping repeat were prepared and did not show this intermediate in equilibrium unfolding experiments. Conversely, the capping repeats were found to be essential for efficient folding in the cell and for avoiding aggregation, presumably because of their highly charged surface. To design a capping repeat conferring similar solubility properties yet even higher stability, eight point mutations adapting the C-cap to the consensus AR and adding a three-residue extension at the C-terminus were predicted in silico and validated experimentally. The in vitro full-consensus proteins were also compared with a previously published designed AR protein, E3_5, whose internal repeats show 80% identity in primary sequence. A detailed analysis of the simulations suggests that networks of salt bridges between beta-hairpins, as well as additional interrepeat hydrogen bonds, contribute to the extraordinary stability of the full consensus.

Very few selective inhibitors of the zeta-chain associated protein kinase 70 kDa (ZAP70) have been reported despite its importance in autoimmune diseases. Here, to induce a fit of the so-called gatekeeper residue (Met414) and hydrophobic pocket next to it, a potent Janus kinase 2 (JAK2) inhibitor was first docked into the ATP binding site of ZAP70 by structural alignment of the kinase domains. The resulting model of the complex between ZAP70 and the JAK2 inhibitor was then relaxed by an explicit solvent molecular dynamics simulation with restraints on the backbone atoms. High-throughput docking into the induced-fit conformation of ZAP70 generated by molecular dynamics has revealed 10 low-micromolar inhibitors which correspond to six distinct chemotypes. One of these ZAP70 inhibitors has an IC50 of 110 nM for JAK2.

Bromodomains are protein modules that selectively recognize histones by binding to acetylated lysines. Here, we have carried out multiple molecular dynamics simulations of 20 human bromodomains to investigate the flexibility of their binding site. Some bromodomains show alternative side chain orientations of three evolutionarily conserved residues: the Asn involved in acetyl-lysine binding and two conserved aromatic residues. Furthermore, for the BAZ2B and CREBBP bromodomains we observe occlusion of the binding site which is coupled to the displacement of the two aromatic residues. In contrast to available structures, the simulations reveal large variability of the binding site accessibility. The simulations suggest that the flexibility of the bromodomain binding site and presence of self-occluded metastable states influence the recognition of acetyl-lysine on histone tails.

The ankyrin repeat is one of the most common protein motifs and is involved in protein-protein interactions. It consists of 33 residues that assume a beta-hairpin helix-loop-helix fold. Mutagenesis and kinetic experiments (Phi-value analysis of the folding transition state) have shown that the tumor suppressor p16(INK4a), a four-repeat protein, unfolds sequentially starting from the two N-terminal repeats. Here, the flexibility of p16(INK4a) at room temperature and its unfolding mechanism at high temperature have been investigated by multiple molecular dynamics runs in explicit water for a total simulation time of 0.65 micros. The transition state ensemble (TSE) of p16(INK4a) was identified by monitoring both the deviation from the experimental Phi values and sudden conformational changes along the unfolding trajectories. Conformations in the TSE have a mainly unstructured second repeat whereas the other repeats are almost completely folded. A rigid-body displacement of the first repeat involving both a rotation and translation is observed in all molecular dynamics simulations at high temperature. The Trp(15), Pro(75), and Ala(76) side-chains are more buried in the TSE than the native state. The sequential unfolding starting at the second repeat is in agreement with the mutagenesis studies whereas the displacement of the first repeat and the presence of nonnative interactions at the TSE are simulation results which supplement the experimental data. Furthermore, the unfolding trajectories reveal the presence of two on-pathway intermediates with partial alpha-helical structure. Finally, on the basis of the available experimental and simulation results we suggest that in modular proteins the shift of the folding TSE toward the native structure upon reduction of the number of tandem repeats is consistent with the Hammond effect.

The aspartic protease beta-secretase (BACE) cleaves the amyloid precursor protein into a 42 residue beta-peptide, which is the principal biochemical marker of Alzheimer's disease. Multiple explicit-water molecular dynamics simulations of the apo and inhibitor bound structures of BACE indicate that both open- and closed-flap conformations are accessible at room temperature and should be taken into account for inhibitor design. Correlated motion is observed within each of the two lobes of BACE, as well as for the interfacial region. A self-inhibited conformation with the side chain of Tyr71 occupying the S(1) pocket is present in some of the unbound simulations. The reversible loss of the side chain hydrogen bond between the catalytic Asp32 and Ser35, due to the concomitant reorientation of the Ser35 hydroxyl group and a water molecule conserved in pepsin-like enzymes, provides further evidence for the suggestion that Ser35 assists in proton acceptance and release by Asp32 during catalysis.

As part of the cytoskeleton, actin is essential for the morphology, motility, and division of eukaryotic cells. Recent X-ray fiber diffraction studies have shown that the conformation of monomeric actin is flattened upon incorporation into the filament by a relative rotation of its two major domains. The antiproliferative activity of latrunculin, a macrolide toxin produced by sponges, seems to be related to its binding to monomeric actin and inhibition of polymerization. Yet, the mechanism of inhibition is not known in detail. Here, multiple explicit water molecular dynamics simulations show that latrunculin binding hinders the conformational transition related to actin polymerization. In particular, the presence of latrunculin at the interface of the two major domains of monomeric actin reduces the correlated displacement of Domain 2 with respect to Domain 1. Moreover, higher rotational flexibility between the two major domains is observed in the absence of ATP as compared to ATP-bound actin, offering a possible explanation as to why actin polymerizes more favorably in the absence of nucleotides. Proteins 2012; © 2012 Wiley Periodicals, Inc.

An improved version of the fragment-based flexible ligand docking approach SEED–FFLD is tested on inhibitors of human immunodeficiency virus type 1 protease, human α-thrombin and the estrogen receptor β. The docking results indicate that it is possible to correctly reproduce the binding mode of inhibitors with more than ten rotatable bonds if the strain in their covalent geometry upon binding is not large. A high degree of convergence towards a unique binding mode in multiple runs of the genetic algorithm is proposed as a necessary condition for successful docking. © 2003 Wiley Periodicals, Inc. J Comput Chem 25: 412–422, 2004

A procedure based on semi-empirical quantum mechanical (QM) calculations of interaction energy is proposed for the rapid screening of compound poses generated by high-throughput docking. Small molecules (consisting of 2-10 atoms and termed "probes") are overlapped with polar groups in the binding site of the protein target. The interaction energy values between each compound pose and the probes, calculated by a semi-empirical Hamiltonian, are used as filters. The QM probe method does not require fixed partial charges and takes into account polarization and charge-transfer effects which are not captured by conventional force fields. The procedure is applied to screen approximately 100 million poses (of 2.7 million commercially available compounds) obtained by high-throughput docking in the ATP binding site of the tyrosine kinase erythropoietin-producing human hepatocellular carcinoma receptor B4 (EphB4). Three QM probes on the hinge region and one at the entrance pocket are employed to select for binding affinity, while a QM probe on the side chain of the so-called gatekeeper residue (a hypervariable residue in the kinome) is used to enforce selectivity. The poses with favorable interactions with the five QM probes are filtered further for hydrophobic matching and low ligand strain. In this way, a single-digit micromolar inhibitor of EphB4 with a relatively good selectivity profile is identified in a multimillion-compound library upon experimental tests of only 23 molecules.

High throughput docking (HTD) using high performance computing platforms is a multidisciplinary challenge. To handle HTD data effectively and efficiently, we have developed a distributed virtual screening data management system (DVSDMS) in which the data handling and the distribution of jobs are realized by the open-source structured query language database software MySQL. The essential concept of DVSDMS is the separation of the data management from the docking and ranking applications. DVSDMS can be used to dock millions of molecules effectively, monitor the process in real time, analyze docking results promptly, and process up to 10(8) poses by energy ranking techniques. In an HTD campaign to identify kinase inhibitors a low cost Linux PC has allowed DVSDMS to efficiently assign the workload to more than 500 computing clients. Notably, in a stress test of DVSDMS that emulated a large number of clients, about 60 molecules per second were distributed to the clients for docking, which indicates that DVSDMS can run efficiently on very large compute cluster (up to about 40000 cores).

Plasmepsins (PMs) are essential proteases of the plasmodia parasites and are therefore promising targets for developing drugs against malaria. We have discovered six inhibitors of PM II by high-throughput fragment-based docking of a diversity set of approximately 40,000 molecules, and consensus scoring with force field energy functions. Using the common scaffold of the three most active inhibitors (IC(50)=2-5 microM), another seven inhibitors were identified by substructure search. Furthermore, these 13 inhibitors belong to at least three different classes of compounds. The in silico approach was very effective since a total of 13 active compounds were discovered by testing only 59 molecules in an enzymatic assay. This hit rate is about one to two orders of magnitude higher than those reported for medium- and high-throughput screening techniques in vitro. Interestingly, one of the inhibitors identified by docking was halofantrine, an antimalarial drug of unknown mechanism. Explicit water molecular dynamics simulations were used to discriminate between two putative binding modes of halofantrine in PM II.

Rational ligand design is a complex problem that can be divided into three parts: the search for optimal positions and orientations of functional groups in the binding site, the connection of such positions to form candidate ligands, and the estimation of their binding constants. Approaches for addressing the first two parts of the problem are described in the present work. They are applied to the construction of peptide ligands in the binding site of the human immunodeficiency virus 1 (HIV-1) proteinase. The primary objective is to test the method by comparison of the results with the MVT-101 complex structure for which coordinates are available; the results obtained with the liganded and unliganded proteinase structure are used to examine the utility of the latter for binding studies. A secondary objective is to show how to find new inhibitor candidates. The multiple copy simultaneous search (MCSS) method is utilized to search for optimal positions and orientations of a set of functional groups. For peptide ligands, functional groups corresponding to the protein main chain (N-methylacetamide) and to protein side chains (e.g., methanol, ethyl guanidinium) are used. The resulting N-methylacetamide minima are connected to form hexapeptide main chains with a simple pseudoenergy function that permits a complete search of all possible ways of connecting the minima. Side chains are added to the main-chain candidates by application of the same pseudoenergy function to the appropriate functional group minima. A set of 15 hexapeptides with the sequence of MVT-101 is then minimized by a Monte Carlo scheme, which allows for escape from local minima. Comparison of the MCSS results with the structure of MVT-101 in the HIV-1 binding site showed that all of its functional group positions correspond (within 2.4 A) to some (usually more than one) MCSS minima. There were also many other low-energy MCSS minima which do not appear in any known inhibitors, e.g., methyl ammonium minima in the neighborhood of the catalytic aspartates. Among the 15 lowest minima are seven hexapeptides with the same main-chain orientation as the one found by X-ray crystallography for the inhibitor MVT-101 in the binding site and eight with the main chain oriented in the opposite direction; the latter tend to be more stable. [Addendum: These results are in agreement with recent high-resolution crystallographic data provided after the study was completed.(ABSTRACT TRUNCATED AT 400 WORDS)

Motivation and method: Small-molecule inhibitors targeting the adenosine triphosphate (ATP) binding pocket of the catalytic domain of protein kinases have potential to become drugs devoid of (major) side effects, particularly if they bind selectively. Here, the sequences of the 518 human kinases are first mapped onto the structural alignment of 116 kinases of known three-dimensional structure. The multiple structure alignment is then used to encode the known strategies for developing selective inhibitors into a fingerprint. Finally, a network analysis is used to partition the kinases into clusters according to similarity of their fingerprints, i.e. physico-chemical characteristics of the residues responsible for selective binding. Results: For each kinase the network analysis reveals the likelihood to find selective inhibitors targeting the ATP binding site. Systematic guidelines are proposed to develop selective inhibitors. Importantly, the network analysis suggests that the tyrosine kinase EphB4 has high selectivity potential, which is consistent with the selectivity profile of two novel EphB4 inhibitors. Contact: dhuang@bioc.uzh.ch; caflisch@bioc.uzh.ch Supplementary information: Supplementary data are available at Bioinformatics online.

In drug discovery, the occurrence of false positives is a major hurdle in the search for lead compounds that can be developed into drugs. A small-molecular-weight compound that inhibits dengue virus protease at low micromolar levels was identified in a screening campaign. Binding to the enzyme was confirmed by isothermal titration calorimetry (ITC) and nuclear magnetic resonance (NMR). However, a structure-activity relationship study that ensued did not yield more potent leads. To further characterize the parental compound and its analogues, we developed a high-speed, low-cost, quantitative fluorescence quenching assay. We observed that specific analogues quenched dengue protease fluorescence and showed variation in IC(50) values. In contrast, nonspecifically binding compounds did not quench its fluorescence and showed similar IC(50) values with steep dose-response curves. We validated the assay using single Trp-to-Ala protease mutants and the competitive protease inhibitor aprotinin. Specific compounds detected in the binding assay were further analyzed by competitive ITC, NMR, and surface plasmon resonance, and the assay's utility in comparison with these biophysical methods is discussed. The sensitivity of this assay makes it highly useful for hit finding and validation in drug discovery. Furthermore, the technique can be readily adapted for studying other protein-ligand interactions.

The two-component NS2B-NS3 protease of West Nile virus is essential for its replication and presents an attractive target for drug development. Here, we describe protocols for the high-yield expression of stable isotope-labelled samples in vivo and in vitro. We also describe the use of NMR spectroscopy to determine the binding mode of new low molecular mass inhibitors of the West Nile virus NS2B-NS3 protease which were discovered using high-throughput in vitro screening. Binding to the substrate-binding sites S1 and S3 is confirmed by intermolecular NOEs and comparison with the binding mode of a previously identified low molecular mass inhibitor. Our results show that all these inhibitors act by occupying the substrate-binding site of the protease rather than by an allosteric mechanism. In addition, the NS2B polypeptide chain was found to be positioned near the substrate-binding site, as observed previously in crystal structures of the protease in complex with peptide inhibitors or bovine pancreatic trypsin inhibitor. This indicates that the new low molecular mass compounds, although inhibiting the protease, also promote the proteolytically active conformation of NS2B, which is very different from the crystal structure of the protein without inhibitor.

A fragment-based docking procedure followed by substructure search were used to identify active-site beta-secretase inhibitors from a composite set of about 300 000 and a library of nearly 180 000 small molecules, respectively. EC(50) values less than 10 microM were measured in at least one of two different mammalian cell-based assays for 12 of the 72 purchased compounds. In particular, the phenylureathiadiazole 2 and the diphenylurea derivative 3 are promising lead compounds for beta-secretase inhibition.

The relative importance of amino acid sequence and native topology in the unfolding process of two SH3 domains and two circular permutants was investigated by 120 molecular dynamics runs at 375 K for a total simulation time of 0.72 micros. The alpha-spectrin (aSH3) and src SH3 (sSH3) domains, which have the same topology and a sequence identity of only 34%, show similar unfolding pathways. The disappearance of the three-stranded antiparallel beta-sheet is the last unfolding event, in agreement with a large repertoire of kinetic data derived from point mutations as well as glycine insertions and disulfide crosslinks. Two alternative routes of beta-sheet unfolding have emerged from the analysis of the trajectories. One is statistically preferred in aSH3 (n-src loop breaks before distal hair-pin) and the inverse in sSH3. An elongation of the beta2-beta3 hairpin was observed during the unfolding of sSH3 at 375 K and in 300 K simulations started from the putative transition state of sSH3 in accord with unusual kinetic data for point mutations at the n-src loop. The change of connectivity in the permutants influenced the sequence of unfolding events mainly at the permutation site. Regions where the connectivity remained unaffected showed the same chronology of contact disappearance. Taken together with previous folding simulations of two designed three-stranded antiparallel beta-sheet peptides, these results indicate that, at least for small beta-sheet proteins, the folding mechanism is primarily defined by the native state topology, whilst specific interactions determine the statistically predominant folding route.

A new method is presented to quan- titatively estimate and graphically display the pro- pensity of nonpolar groups to bind at the surface of proteins. It is based on the calculation of the bind- ing energy, i.e., van der Waals interaction plus pro- tein electrostatic desolvation, of a nonpolar probe sphere rolled over the protein surface, and on the color coding of this quantity on a smooth molecular surface (hydrophobicity map). The method is vali- dated on ten protein-ligand complexes and is shown to distinguish precisely where polar and nonpolar groups preferentially bind. Comparisons with exist- ing approaches, like the display of the electrostatic potential or the curvature, illustrate the advantages and the better predictive power of the present method. Hydrophobicity maps will play an impor- tant role in the characterization of binding sites for the large number of proteins emerging from the genome projects and structure modeling ap- proaches. Proteins 1999;37:565-575.

A new method is presented for computer-aided ligand design by combinatorial selection of fragments that bind favorably to a macromolecular target of known three-dimensional structure. Firstly, the multiple-copy simultaneous-search procedure (MCSS) is used to exhaustively search for optimal positions and orientations of functional groups on the surface of the macromolecule (enzyme or receptor fragment). The MCSS minima are then sorted according to an approximated binding free energy, whose solvation component is expressed as a sum of separate electrostatic and nonpolar contributions. The electrostatic solvation energy is calculated by the numerical solution of the linearized Poisson-Boltzmann equation, while the nonpolar contribution to the binding free energy is assumed to be proportional to the loss in solvent-accessible surface area. The program developed for computational combinatorial ligand design (CCLD) allows the fast and automatic generation of a multitude of highly diverse compounds, by connecting in a combinatorial fashion the functional groups in their minimized positions. The fragments are linked as two atoms may be either fused, or connected by a covalent bond or a small linker unit. To avoid the combinatorial explosion problem, pruning of the growing ligand is performed according to the average value of the approximated binding free energy of its fragments. The method is illustrated here by constructing candidate ligands for the active site of human alpha-thrombin. The MCSS minima with favorable binding free energy reproduce the interaction patterns of known inhibitors. Starting from these fragments, CCLD generates a set of compounds that are closely related to high-affinity thrombin inhibitors. In addition, putative ligands with novel binding motifs are suggested. Probable implications of the MCSS-CCLD approach for the evolving scenario of drug discovery are discussed.