Spin alignments in heterospin chains are examined from numerical calculations of model spin Hamiltonians. The Hamiltonians of the heterospin chains mimic an open-shell molecular assemblage composed of an organic biradical in a singlet (S = 0) ground state and a doublet (S = 1/2) monoradical, which are coupled by intermolecular ferromagnetic exchange interactions. It is found from numerical calculations of the spin Hamiltonians that the spin value &2rattach; S-2&2lattach; of the ground-state singlet biradical embedded in the exchange-coupled assemblage deviates from zero and contributes to the bulk magnetization. The alternating chain is found to have two kinds of ground spin states, a high- and a low-spin state. All the spins are parallel to each other in the high-spin state, which is characterized by the spin correlation function of &2rattach; S(i)center dot S-j&2lattach; = 0.25. On the other hand, the spin alignment in the low-spin state is found to be dependent on the topology of the intermolecular exchange interactions. The energy preference of the two states depends on the relative amplitude of the exchange interactions in the chain. The intermolecular ferromagnetic couplings are competing in the ground-state singlet biradical with the intramolecular antiferromagnetic interaction. The appearance of the two kinds of ground states is attributed to a quantum spin frustration effect inherent in the triangular motif of the competing interactions. Magnetic properties of a zigzag chain complex composed of a nitronyl nitroxide biradical with a singlet ground state and Cu(hfac)(2) are examined on the basis of the theoretical calculations. The vanishing magnetic moments, or the product of susceptibility and temperature chi T, at low temperatures observed for the complex are consistent with those of the low-spin state predicted in the theoretical calculations.