The application of structure-based in silico methods to drug discovery is still considered a major challenge, especially when the x-ray structure of the target protein is unknown. Such is the case with human G protein-coupled receptors (GPCRs), one of the most important families of drug targets, where in the absence of x-ray structures, one has to rely on in silico 3D models. We report repeated success in using ab initio in silico GPCR models, generated by the PREDICT method, for blind in silico screening when applied to a set of five different GPCR drug targets. More than 100,000 compounds were typically screened in silico for each target, leading to a selection of <100 ''virtual hit'' compounds to be tested in the lab. In vitro binding assays of the selected compounds confirm high hit rates, of 12-21% (full dose-response curves, K i < 5 M). In most cases, the best hit was a novel compound (New Chemical Entity) in the 1-to 100-nM range, with very promising pharmacological properties, as measured by a variety of in vitro and in vivo assays. These assays validated the quality of the hits as lead compounds for drug discovery. The results demonstrate the usefulness and robustness of ab initio in silico 3D models and of in silico screening for GPCR drug discovery.modeling ͉ in silico screening ͉ structure-based G protein-coupled receptors (GPCRs) are membraneembedded proteins, responsible for communication between the cell and its environment (1). As a consequence, many major diseases, such as hypertension, cardiac dysfunction, depression, anxiety, obesity, inflammation, and pain, involve malfunction of these receptors (2), making them among the most important drug targets for pharmacological intervention (3-5). Thus, whereas GPCRs are only a small subset of the human genome, they are the targets for Ϸ50% of all recently launched drugs (6). As targets of paramount importance, it is expected that drug discovery for GPCRs would benefit from the introduction of computational methodologies (7), especially as these methods can be used in conjunction with such experimental methods as high-throughput screening (8, 9), NMR, and crystallography (10).Unfortunately, GPCRs, like other membrane-embedded proteins, have characteristics that make their 3D structure extremely difficult to determine experimentally. To date, the only GPCR for which a 3D structure was determined by x-ray crystallography is bovine rhodopsin (11), which is unique among GPCRs in that its ligand, retinal, is covalently bound and that it responds to light rather than to ligand binding. Hence, in the case of GPCRs, the limited availability of structural data has forced the computational design of ligands to heavily rely on ligand-based techniques. Indeed, for many GPCRs, the natural ligand can provide a good starting point, leading to useful pharmacophore models that can be used for identifying lead structures with novel scaffolds (6). These methods have been successfully applied for the discovery of peptide agonists to the somatostatin receptor (12) and for...