Elucidating TH301's influence on the torsion angle of CRY1 W399 using replica exchange with solute tempering (REST) molecular dynamics (MD) simulations

Abstract

Cryptochrome 1 (CRY1) is a protein involved in the circadian clock and associated with various diseases. Targeting CRY1 for drug development requires the discovery of competitive inhibitors that target its FAD binding site through ubiquitination. During the development of compounds to regulate CRY1, an intriguing compound called TH301 was identified. Despite binding to CRY1, TH301 does not induce the expected reaction and is considered an inactive compound. However, it has been observed that TH301 affects the torsion angle of CRY1's W399 residue, which plays a crucial role in the regulation of ubiquitination by influencing the movement of the lid loop. In our research, we aimed to understand how TH301 induces the torsion angle of CRY1's W399 to shift to an “out-form” by performing REST-based MD simulations. The cyclopentane of TH301 tends to align parallel with W292, creating a repulsive force when W399 is in the “in-form”, leading to a flip. In the “out-form”, W399's side chain interacts with TH301's chlorobenzene through a π–π interaction, stabilizing this pose. This analysis helps identify compounds binding to CRY1 and filter out inactive ones. We found that assessing the interaction energy between TH301 and W399 is crucial to evaluate whether W399 flips or not. These findings contribute to the development of drugs targeting CRY1 and enhance our understanding of its regulatory mechanisms.