Untargeted metabolomics aims to provide a global picture of the metabolites present in the system under study. To this end, making a careful choice of sample preparation is mandatory to obtain reliable and reproducible biological information. In this study, eight different sample preparation techniques were evaluated using Caulobacter crescentus as a model for Gram-negative bacteria. Two cell retrieval systems, two quenching and extraction solvents, and two cell disruption procedures were combined in a full factorial experimental design. To fully exploit the multivariate structure of the generated data, the ANOVA multiblock orthogonal partial least squares (AMOPLS) algorithm was employed to decompose the contribution of each factor studied and their potential interactions for a set of annotated metabolites. All main effects of the factors studied were found to have a significant contribution on the total observed variability. Cell retrieval, quenching and extraction solvent, and cell disrupting mechanism accounted respectively for 27.6%, 8.4%, and 7.0% of the total variability. The reproducibility and metabolome coverage of the sample preparation procedures were then compared and evaluated in terms of relative standard deviation (RSD) on the area for the detected metabolites. The protocol showing the best performance in terms of recovery, versatility, and variability was centrifugation for cell retrieval, using MeOH:H2O (8:2) as quenching and extraction solvent, and freeze-thaw cycles as the cell disrupting mechanism.