The present work introduces a novel method for drug research based on the sequential building of linked multivariate statistical models, each one introducing a different level of drug description. The use of multivariate methods allows us to overcome the traditional one-target assumption and to link in vivo endpoints with drug binding profiles, involving multiple receptors. The method starts with a set of drugs, for which in vivo or clinical observations and binding affinities for potentially relevant receptors are known, and allows obtaining predictions of the in vivo endpoints highlighting the most influential receptors. Moreover, provided that the structure of the receptor binding sites is known (experimentally or by homology modeling), the proposed method also highlights receptor regions and ligandreceptor interactions that are more likely to be linked to the in vivo endpoints, which is information of high interest for the design of novel compounds. The method is illustrated by a practical application dealing with the study of the metabolic side effects of antipsychotic drugs. Herein, the method detects related receptors confirmed by experimental results. Moreover, the use of structural models of the receptor binding sites allows identifying regions and ligand-receptor interactions that are involved in the discrimination between antipsychotic drugs that show metabolic side effects and those that do not. The structural results suggest that the topology of a hydrophobic sandwich involving residues in transmembrane helices (TM) 3, 5, and 6, as well as the assembling of polar residues in TM5, are important discriminators between target/antitarget receptors. Ultimately, this will provide useful information for the design of safer compounds inducing fewer side effects.The classic formulation in which a disease is associated to a single biological receptor is clearly an oversimplification, valid only for exceptional situations. Recent advances in genomics, proteomics, and systems biology depict a far more complex scenario, much less comfortable to work with, in which the biochemical mechanisms of diseases and therapies involve many different biological receptors linked in complex interaction pathways, which are usually poorly understood. The best possible therapies for many complex diseases will probably require a deep and detailed understanding of those relationships. However, because of the extreme complexity of the involved mechanisms, we will not achieve such detailed understanding in the immediate future. Hence, there is a need for new methodological approaches that are able to advance in this way. In this work, we propose a novel method for exploiting pharmacologic information obtained for sets of currently available drugs that are known to interact with multiple related receptors. This method is intended to improve our understanding of the biological mechanisms implicated in the pharmacological treatments but also to provide hints about how the chemical structure of the studied drugs could be improve...