A predictive group-contribution statistical associating fluid theory (SAFT-gamma) is developed by extending the molecular-based SAFT-VR equation of state [A. Gil-Villegas et al. J. Chem. Phys. 106, 4168 (1997)] to treat heteronuclear molecules which are formed from fused segments of different types. Our models are thus a heteronuclear generalization of the standard models used within SAFT, comparable to the optimized potentials for the liquid state OPLS models commonly used in molecular simulation; an advantage of our SAFT-gamma over simulation is that an algebraic description for the thermodynamic properties of the model molecules can be developed. In our SAFT-gamma approach, each functional group in the molecule is modeled as a united-atom spherical (square-well) segment. The different groups are thus characterized by size (diameter), energy (well depth) and range parameters representing the dispersive interaction, and by shape factor parameters (which denote the extent to which each group contributes to the overall molecular properties). For associating groups a number of bonding sites are included on the segment: in this case the site types, the number of sites of each type, and the appropriate association energy and range parameters also have to be specified. A number of chemical families (n-alkanes, branched alkanes, n-alkylbenzenes, mono- and diunsaturated hydrocarbons, and n-alkan-1-ols) are treated in order to assess the quality of the SAFT-gamma description of the vapor-liquid equilibria and to estimate the parameters of various functional groups. The group parameters for the functional groups present in these compounds (CH(3), CH(2), CH(3)CH, ACH, ACCH(2), CH(2)=, CH=, and OH) together with the unlike energy parameters between groups of different types are obtained from an optimal description of the pure component phase equilibria. The approach is found to describe accurately the vapor-liquid equilibria with an overall %AAD of 3.60% for the vapor pressure and 0.86% for the saturated liquid density. The fluid phase equilibria of some larger compounds comprising these groups, which are not included in the optimization database and some binary mixtures are examined to confirm the predictive capability of the SAFT-gamma approach. A key advantage of our method is that the binary interaction parameters between groups can be estimated directly from an examination of pure components alone. This means that as a first approximation the fluid-phase equilibria of mixtures of compounds comprising the groups considered can be predicted without the need for any adjustment of the binary interaction parameters (which is common in other approaches). The special case of molecular models comprising tangentially bonded (all-atom and united-atom) segments is considered separately; we comment on the adequacy of such models in representing the properties of real molecules.