Characterizing the conformational ensemble of PROTAC degraders in solutions <i>via</i> atomistic simulations

Targeted protein degradation has emerged as a promising strategy for developing novel therapeutics, particularly for "undruggable" disease-related proteins. One approach is the use of PROteolysis TArgeting Chimeras (PROTACs) degraders, which induce the formation of ternary complexes between the target protein and E3 ligase, leading to ubiquitination and degradation of the target protein. Understanding the conformational behavior of PROTACs in solutions and how it relates to their pharmacokinetic properties and membrane permeability is crucial for optimizing PROTAC design and efficacy. Due to the large size and flexibility of PROTACs and their chameleonic character, it is essential to understand their conformational ensembles, and how they depend on the environment. Here, we introduce a novel methodology for exploring PROTAC conformational behavior using atomistic simulations. We employ the enhanced sampling method parallel bias metadynamics, where we bias generic local collective variables, specifically all rotatable dihedral angles, thereby avoiding the considerable challenge of identifying suitable global collective variables for biasing. The methodology allows for obtaining free energy surfaces of global CVs <i>via</i> reweighting. We apply the method to the prototypical case of the MZ1 PROTAC degrader, which targets bromodomain-containing protein-4 (BRD4) for degradation <i>via</i> the von Hippel-Lindau (VHL) E3 ligase, and elucidate its conformational behavior in different solvents, allowing us to gain insights into the chameleonic property of MZ1. Our results confirm that MZ1 adopts distinct conformations depending on the solvent, exhibiting collapsed conformations in water and chloroform, and extended conformations in DMSO. Collapsed conformations in chloroform have been correlated to increased cell permeability. Thus, our results show that MZ1 takes on conformations suitable for membrane permeation in apolar environments. Our methodological framework is generally applicable to large flexible molecules like PROTACs, and the results demonstrate its efficiency, laying the groundwork for similar investigations for other PROTACs and other "beyond-rule-of-5" drug candidates. This work provides valuable insights into the design and optimization of PROTACs, ultimately contributing to developing novel therapeutics for "undruggable" proteins.