see Re: Hexamer combination from 3mers. about Algorithm::Combinatorics.

Understood. But that doesn't solve the "array of arrays, 1 from each" problem. A shorter version of the OP glob spec would be

@x = (1..2);
@y = qw(one two);
@z = qw(Amy Bob);

@xyz = (\@x, \@y, \@z);
$glob = join( ',', map { '{' . join( ',', @$_ ) . '}' } @xyz );
print "$glob\n";

# prints '{1,2},{one,two},{Amy,Bob}'
[download]

I don't know of a module that can emulate glob on that.

Actually I think your iterator solution is not half bad. It could work very well where nested foreach loops are not an option.

I'd drop the redundant array_limits array, however.