Chloroquine (CQ) and hydroxychloroquine (HCQ) have been standard antimalarial drugs since the early 1950s, and very recently, the possibility of their use for the treatment of COVID-19 patients has been considered. To understand the drug mode of action at the submicroscopic level (atoms and molecules), molecular modeling studies with the aid of computational chemistry methods have been of great help. A fundamental step in such theoretical investigations is the knowledge of the predominant drug molecular structure in solution, which is the real environment for the interaction with biological targets. Our strategy to access this valuable information is to perform density functional theory (DFT) calculations of 1H NMR chemical shifts for several plausible molecular conformers and then find the best match with experimental NMR profile in solution (since it is extremely sensitive to conformational changes). Through this procedure, after optimizing 30 trial distinct molecular structures (ωB97x-D/6-31G(d,p)-PCM level of calculation), which may be considered representative conformations, we concluded that the global minimum (named M24), stabilized by an intramolecular N-H hydrogen bond, is not likely to be observed in water, chloroform, and dimethyl sulfoxide (DMSO) solution. Among fully optimized conformations (named M1 to M30, and MD1 and MD2), we found M12 (having no intramolecular H-bond) as the most probable structure of CQ and HCQ in water solution, which is a good approximate starting geometry in drug-receptor interaction simulations. On the other hand, the preferred CQ and HCQ structure in chloroform (and CQ in DMSO-d6) solution was assigned as M8, showing the solvent effects on conformational preferences. We believe that the analysis of 1H NMR data in solution can establish the connection between the macro level (experimental) and the sub-micro level (theoretical), which is not so apparent to us and appears to be more appropriate than the thermodynamic stability criterion in conformational analysis studies.