Surface free energy of a chemically heterogeneous surface is often treated as an approximately additive quantity through the Cassie equation [Cassie ABD (1948) Discuss Faraday Soc 3:11-16]. However, deviations from additivity are common, and molecular interpretations are still lacking. We use molecular simulations to measure the microscopic analogue of contact angle, θ c , of aqueous nanodrops on heterogeneous synthetic and natural surfaces as a function of surface composition. The synthetic surfaces are layers of graphene functionalized with prototypical nonpolar and polar head group: methyl, amino, and nitrile. We demonstrate positive as well as negative deviations from the linear additivity. We show the deviations reflect the uneven exposure of mixture components to the solvent and the linear relation is recovered if fractions of solvent-accessible surface are used as the measure of composition. As the spatial variations in polarity become of larger amplitude, the linear relation can no longer be obtained. Protein surfaces represent such natural patterned surfaces, also characterized by larger patches and roughness. Our calculations reveal strong deviations from linear additivity on a prototypical surface comprising surface fragments of melittin dimer. The deviations reflect the disproportionately strong influence of isolated polar patches, preferential wetting, and changes in the position of the liquid interface above hydrophobic patches. Because solvent-induced contribution to the free energy of surface association grows as cos θ c , deviations of cos θ c from the linear relation directly reflect nonadditive adhesive energies of biosurfaces. wetting free energy | surface functionalization | nanopatterning | Cassie relation | biointeractions W etting phenomena on chemically heterogeneous surfaces are important in material sciences and biology, in examples ranging from inkjet printing to protein hydration (1, 2). Conventional metrics of surface interactions, designed for homogeneous systems, can often be applied to mixed surfaces characterized by averaged properties of multiple ingredients. The design of composite surface materials and characterization of biosurfaces benefit from combining rules predicting the interfacial free-energy change of wetting, Δγ, from the knowledge about individual constituents and surface composition. In view of the Young equation, Δγ ¼ −γ cos θ c , the strength of intersolute adhesion, W a ∼ −2ðΔγ þ γ sg Þ, relates to contact angle θ c ; here, γ and γ sg denote surface tensions of the solvent and dry solute, respectively (3). Contact angles on macroscopic heterogeneous surfaces are commonly estimated by the Cassie equation (4, 5),[1]developed by assuming linear additivity of the wetting free energies, Δγ. Here, f α is the projected fractional area occupied by component α, θ α is the contact angle on a homogeneous surface of type α, and r α is the Wenzel roughness factor (6), which can be defined as the ratio of solvent-exposed areas of patch of type α in the mixture and that of a ...