By taking advantage of the crystallographic data of 3-hydroxy-3-methylglutaryl coenzyme A reductase (HMGR) complexed with statins, a quantum biochemistry study based on the density functional theory is performed to estimate the interaction energy for each statin when one considers binding pockets of different sizes. Assuming a correlation between statin potency and the strength of the total HMGR–statin binding energy, clinical data as well as IC₅₀ values of these cholesterol-lowering drugs are successfully explained only after stabilization of the calculated total binding energy for a larger size of the ligand-interacting HGMR region, one with a radius of at least 12.0 Å. Actually, the binding pocket radius suggested by classic works, which was based solely on the interpretation of crystallographic data of the HMGR–statin complex, is smaller than that necessary to achieve total binding energy convergence in our simulations. Atorvastatin and rosuvastatin are shown to be the most strongly bound HMGR inhibitors, while simvastatin and fluvastatin are the weakest ones. A binding site, interaction energy between residues and statin atoms, and residues domain (BIRD) panel is constructed, indicating clear quantum biochemistry-based routes for the development of new statin derivatives.

Results of classical force field geometry optimizations for twisted graphene nanoribbons with a number of twists Nt varying from 0 to 7 (the case Nt=1 corresponds to a half-twist Möbius nanoribbon) are presented in this work. Their structural stability was investigated using the Brenner reactive force field. The best classical molecular geometries were used as input for semiempirical calculations, from which the electronic properties (energy levels, HOMO, LUMO orbitals) were computed for each structure. CI wavefunctions were also calculated in the complete active space framework taking into account eigenstates from HOMO−4 to LUMO+4, as well as the oscillator strengths corresponding to the first optical transitions in the UV-VIS range. The lowest energy molecules were found less symmetric than initial configurations, and the HOMO-LUMO energy gaps are larger than the value found for the nanographene used to build them due to electronic localization effects created by the twisting. A high number of twists leads to a sharp increase of the HOMO→LUMO transition energy. We suggest that some twisted nanoribbons could form crystals stabilized by dipolar interactions.