Cocrystals (or multicomponent crystals) have physico-chemical properties that are different from crystals of pure components. This is significant in drug development, since the desired properties, e.g. solubility, stability and bioavailability, can be tailored by binding two substances into a single crystal without chemical modification of an active component. Here, the FLEXCRYST program suite, implemented with a data mining force field, was used to estimate the relative stability and, consequently, the relative solubility of cocrystals of flavonoids vs their pure crystals, stored in the Cambridge Structural Database. The considerable potency of this approach for in silico screening of cocrystals, as well as their relative solubility, was demonstrated. (C) 2013 Elsevier B.V. All rights reserved.