Sample preparation constitutes a crucial and limiting step in structural studies of proteins by NMR. The determination of the solubility and stability (SAS) conditions of biomolecules at millimolar concentrations stays today empirical and hence time- and material-consuming. Only few studies have been recently done in this field and they have highlighted the interest of using crystallogenesis tools to optimise sample conditions. In this study, we have adapted a method based on incomplete factorial design and making use of crystallisation plates to quantify the influence of physico-chemical parameters such as buffer pH and salts on protein SAS. A description of the experimental set up and an evaluation of the method are given by case studies on two functional domains from the bacterial regulatory protein LicT as well as two other proteins. Using this method, we could rapidly determine optimised conditions for extracting soluble proteins from bacterial cells and for preparing purified protein samples sufficiently concentrated and stable for NMR characterisation. The drastic reduction in the time and number of experiments required for searching protein SAS conditions makes this method particularly well-adapted for a systematic investigation on a large range of physico-chemical parameters.