In the present work, we investigate scaling rules for the ionization cross sections of multicharged ions on molecules of biological interest. The cross sections are obtained from distorted-wave calculations for atomic targets combined with a stoichiometric model for the molecules proposed in [Mendez et al. J. Phys B (2020)]. We examine ions with charges Z from +1 to +8 in five nucleobases -adenine, cytosine, guanine, thymine, uracil-, tetrahydrofuran, pyrimidine, and also in water. We propose a scaling with the ion charge, which is valid in the intermediate to high energy range, i.e. 0.2-5 MeV/amu for oxygen impact. We extend our work to a universal scaling for any ion and molecule, merging the forty ion-molecule systems analyzed here into a single band. Furthermore, our model proved to be valid for other molecules too.The ionization of biological molecules by multicharged ions has gained increasing interest due to medical and environmental reasons [1], from medical treatments [2-4] to contaminant recognition in biological materials [5,6]. Many semiempirical [7] and theoretical efforts are currently being undertaken [8][9][10][11][12][13] to get reliable values for the ionization cross sections of these molecules.Recently, we combined the continue distorted-wave calculations (CDW) for atoms and the simple stoichiometric model (SSM) to approximate the ionization cross sections of complex molecular targets by charged ions [8]. The CDW-SSM approximation showed reasonable results for over a hundred of ion-molecule systems. As expected, in the high energy range (i.e. above 5 MeV/amu) the ionization cross sections present the Z 2 dependence predicted by the first Born approximation. However, at intermediate energies, the dependence with Z is more complex, and non-perturvative models are mandatory.The intention of this letter is to give a follow up of our previous work [8] by proposing a scaling with the ion charge Z of the ionization cross sections of complex molecules, valid at intermediate energies. In general, scaling rules are used as first-order tests in experimental measurements and multipurpose codes. Based on [8], we propose a universal scaling for any ion-target system.At intermediate impact energies, Janev and Presnyakov [14] suggest σ Z versus E Z as the natural reduced form of the ionization cross section σ and the incident ion energy E. Much more recently, Montenegro and co-workers [15,16] proposed an alternative scaling by taking into account that σ is a function of Z 2 E. Their scaling, given bykeeps the Z 2 E relationship for any value of α. In Ref.[15], the authors propose α = 4 3 for ionization of He and H 2 by different charged ions. Combining our recent CDW-SSM results [8] and Eq. (1), we propose here a Z-scaling and implement it for forty collisional systems. The ion-molecule systems are composed of eight targets: the DNA and RNA nucleobases -adenine, cytosine, guanine, thymine, uracil-, tetrahydrofuran (THF), pyrimidine, and water; and five charged ions: H + , He +2 , Be +4 , C +6 , and O +8 . We ...