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-d
6) 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.