Global optimization constitutes an important and fundamental problem in theoretical studies in many chemical fields, such as catalysis, materials or separations problems. In this paper, a novel algorithm has been developed for the global optimization of large systems including neat and ligated clusters in gas phase, and supported clusters in periodic boundary conditions. The method is based on an updated artificial bee colony (ABC) algorithm method, that allows for adaptive-learning during the search process. The new algorithm is tested against four classes of systems of diverse chemical nature: gas phase Au<sub>55</sub>, ligated Au<sub>8</sub><sup>2+</sup>, Au<sub>8</sub> supported on graphene oxide and defected rutile, and a large cluster assembly [Co<sub>6</sub>Te<sub>8</sub>(PEt<sub>3</sub>)<sub>6</sub>][C<sub>60</sub>]<sub>𝑛</sub>, with sizes ranging between 1 to 3 nm and containing up to 1300 atoms. Reliable global minima (GMs) are obtained for all cases, either confirming published data or reporting new lower energy structures. The algorithm and interface to other codes in the form of an independent program, Northwest Potential Energy Search Engine (NWPEsSe), is freely available and it provides a powerful and efficient approach for global optimization of nanosized cluster systems.