for prot in cell.species_types.get(__type=wc_kb.prokaryote.ProteinSpeciesType):
            prot_specie = prot.species.get_or_create(compartment=cytosol)
            conc = round(abs(random.normal(loc=mean_protein_copy_number,scale=15))) / scipy.constants.Avogadro / mean_volume
            cell.concentrations.get_or_create(id='CONC({})'.format(prot_specie.id), species=prot_specie, value=conc, units=unit_registry.parse_units('M'))