The energetics of rotation around single bonds (torsions) is a key determinant of the three-dimensional shape that druglike molecules adopt in solution, the solid state, and in different biological environments, which in turn defines their unique physical and pharmacological properties. Therefore, accurate characterization of torsion angle preference and energetics is essential for the success of computational drug discovery and design. Here, we analyze torsional strain in crystal structures of druglike molecules in Cambridge structure database (CSD) and bioactive ligand conformations in protein data bank (PDB), expressing the total strain energy as a sum of strain energy from constituent rotatable bonds. We utilized cloud computing to generate torsion scan profiles of a very large collection of chemically diverse neutral fragments at DFT(B3LYP)/6-31G*//6-31G** or DFT(B3LYP)/6-31+G*//6-31+G** (for sulfur-containing molecule). With the data generated from these ab initio calculations, we performed rigorous analysis of strain due to deviation of observed torsion angles relative to their ideal gas-phase geometries. Contrary to the previous studies based on molecular mechanics, we find that in the crystalline state, molecules generally adopt low-strain conformations, with median per-torsion strain energy in CSD and PDB under one-tenth and one-third of a kcal/mol, respectively. However, for a small fraction (<5%) of motifs, external effects such as steric hindrance and hydrogen bonds result in strain penalty exceeding 2.5 kcal/mol. We find that due to poor quality of PDB structures in general, bioactive structures tend to have higher torsional strain compared to small-molecule crystal conformations. However, in the absence of structural fitting artifacts in PDB structures, protein-induced strain in bioactive conformations is quantitatively similar to those due to the packing forces in small-molecule crystal structures. This analysis allows us to establish strain energy thresholds to help identify biologically relevant conformers in a given ensemble. The work presented here is the most comprehensive study to date that demonstrates the utility and feasibility of gas-phase quantum mechanics (QM) calculations to study conformational preference and energetics of drug-size molecules. Potential applications of this study in computational lead discovery and structure-based design are discussed.
CRISPR-Cas RNA-guided endonucleases hold great promise for disrupting or correcting genomic sequences through site-specific DNA cleavage and repair. However, the lack of methods for cell- and tissue-selective delivery currently limits both research and clinical uses of these enzymes. We report the design and in vitro evaluation of S. pyogenes Cas9 proteins harboring asialoglycoprotein receptor ligands (ASGPrL). In particular, we demonstrate that the resulting ribonucleoproteins (Cas9-ASGPrL RNP) can be engineered to be preferentially internalized into cells expressing the corresponding receptor on their surface. Uptake of such fluorescently labeled proteins in liver-derived cell lines HEPG2 (ASGPr+) and SKHEP (control; diminished ASGPr) was studied by live cell imaging and demonstrates increased accumulation of Cas9-ASGPrL RNP in HEPG2 cells as a result of effective ASGPr-mediated endocytosis. When uptake occurred in the presence of a peptide with endosomolytic properties, we observed receptor-facilitated and cell-type specific gene editing that did not rely on electroporation or the use of transfection reagents. Overall, these in vitro results validate the receptor-mediated delivery of genome-editing enzymes as an approach for cell-selective gene editing and provide a framework for future potential applications to hepatoselective gene editing in vivo.
Inducing α-helicity through side-chain cross-linking is a strategy that has been pursued to improve peptide conformational rigidity and bio-availability. Here we describe the preparation of small peptides tethered to chiral sulfoxide-containing macrocyclic rings. Furthermore, a study of structure-activity relationships (SARs) disclosed properties with respect to ring size, sulfur position, oxidation state, and stereochemistry that show a propensity to induce α-helicity. Supporting data include circular dichroism spectroscopy (CD), NMR spectroscopy, and a single crystal X-ray structure for one such stabilized peptide. Finally, theoretical studies are presented to elucidate the effect of chiral sulfoxides in inducing backbone α-helicity.
Increased fructose consumption and its subsequent metabolism have been implicated in hepatic steatosis, dyslipidemia, obesity, and insulin resistance in humans. Since ketohexokinase (KHK) is the principal enzyme responsible for fructose metabolism, identification of a selective KHK inhibitor may help to further elucidate the effect of KHK inhibition on these metabolic disorders. Until now, studies on KHK inhibition with small molecules have been limited due to the lack of viable in vivo pharmacological tools. Herein we report the discovery of 12, a selective KHK inhibitor with potency and properties suitable for evaluating KHK inhibition in rat models. Key structural features interacting with KHK were discovered through fragment-based screening and subsequent optimization using structure-based drug design, and parallel medicinal chemistry led to the identification of pyridine 12.
GK activators represent a promising strategy for treatment of T2DM; however, drug candidates have failed in clinical trials due to narrow TI between efficacy and hypoglycemia.
A compact and stable bicyclic bridged ketal was developed as a ligand for the asialoglycoprotein receptor (ASGPR). This compound showed excellent ligand efficiency, and the molecular details of binding were revealed by the first X-ray crystal structures of ligand-bound ASGPR. This analogue was used to make potent di- and trivalent binders of ASGPR. Extensive characterization of the function of these compounds showed rapid ASGPR-dependent cellular uptake in vitro and high levels of liver/plasma selectivity in vivo. Assessment of the biodistribution in rodents of a prototypical Alexa647-labeled trivalent conjugate showed selective hepatocyte targeting with no detectable distribution in nonparenchymal cells. This molecule also exhibited increased ASGPR-directed hepatocellular uptake and prolonged retention compared to a similar GalNAc derived trimer conjugate. Selective release in the liver of a passively permeable small-molecule cargo was achieved by retro-Diels-Alder cleavage of an oxanorbornadiene linkage, presumably upon encountering intracellular thiol. Therefore, the multicomponent construct described here represents a highly efficient delivery vehicle to hepatocytes.
Background: Epoxide hydrolases comprise an important class of enzymes that have critical functions in the detoxification of xenobiotics and regulation of signaling molecules. In addition to epoxides, oxetanes have recently been identified as novel substrates of microsomal epoxide hydrolase (mEH). Oxetanes are common scaffolds used in medicinal chemistry design to improve potency and drug-like properties. Metabolism of oxetanes by mEH can result in high uncertainties in the prediction of human clearance due to extrahepatic contribution and large interindividual variability. Therefore, reducing mEH-mediated oxetane metabolism is highly desirable to minimize its contribution to clearance. Objective: The aim of the study is to evaluate whether quantum mechanical parameters are able to predict the hydrolytic rate of mEH-mediated oxetane metabolism in order to guide medicinal chemistry design in order to minimize mEH contribution to clearance. Methods: Quantum mechanics modeling was used to evaluate the hydrolytic rate of twenty-three oxetanes by mEH. All modeling studies were performed with the Maestro software package. Results: The results show that LUMO energy is highly correlated with the diol formation rate of oxetane hydrolysis by mEH for compounds that are structurally similar, while other quantum mechanical parameters are less predictive. The data suggest that the intrinsic reactivity determines the hydrolytic rate of oxetanes. This occurs when the orientations of the molecules in the mEH active site are similar. Predictions of mEH substrate metabolic rates using LUMO are most accurate when comparing subtle structural changes without drastic changes in MW and chemotype. Conclusion: The study suggests that LUMO energy can be used to rank-order oxetanes for their hydrolytic rate by mEH for structurally similar compounds. This finding enables the medicinal chemistry design to reduce mEH-mediated oxetane metabolism based on the calculated LUMO energy.
'/%3e%3cuse xlink:href='%23a' transform='rotate(23.036 .093 25.536)'/%3e%3cuse xlink:href='%23a' transform='rotate(45.87 1.273 16.18)'/%3e%3cuse xlink:href='%23a' transform='rotate(69.945 .996 12.078)'/%3e%3cuse xlink:href='%23a' transform='rotate(20.66 -19.689 31.932)'/%3e%3c/svg%3e)
'%3e%3cpath fill='%23012169' d='M0 0v30h60V0z'/%3e%3cpath stroke='%23fff' stroke-width='6' d='m0 0 60 30m0-30L0 30'/%3e%3cpath stroke='%23C8102E' stroke-width='4' d='m0 0 60 30m0-30L0 30' clip-path='url(%23b)'/%3e%3cpath stroke='%23fff' stroke-width='10' d='M30 0v30M0 15h60'/%3e%3cpath stroke='%23C8102E' stroke-width='6' d='M30 0v30M0 15h60'/%3e%3c/g%3e%3c/svg%3e)