In this work, we present a user-friendly python package, which builds upon AmberTools, to provide an automated process for generating conformers of arbitrary, metal-containing or metalfree complexes, in explicit solvent. The user only needs to provide an input structure, the desired QM method to be used to the geometry optimization and force constant calculations, total system charge and multiplicity, as well as the solvent to be used to the simulation (if applicable). We provide full support for a large number of pre-parametrized solvents, as well as the ability to parametrize any solvent of choice. Additionally, as many metal-containing complexes are associated with a counterion, we provide seven pre-parametrized ions, as well as the ability to parametrize an ion of choice. By interfacing to ORCA 5, the user has access to state-of-the-art QM methods for structure calculations. For system charge assignment, an interface to Multiwfn 3.8, 53 offers the user a plethora of charge calculations schemes, in addition to the Merz-Kollman RESP scheme recommended by MCPB. After the parametrization is complete, the system is solvated and a suggested equilibration procedure is offered to the user, using the Amber MD package. Once a simulation is performed, the resulting trajectory file can be analyzed via a script provided by PyConSolv. The clustering itself is based on the popular cpptraj package and returns a number of conformers, ranked by their QM energy, based on single point calculations.