In this study we investigate π-stacking interactions of a variety of aromatic heterocycles with benzene using dispersion corrected density functional theory. We calculate extensive potential energy surfaces for parallel-displaced interaction geometries. We find that dispersion contributes significantly to the interaction energy and is complemented by a varying degree of electrostatic interactions. We identify geometric preferences and minimum interaction energies for a set of 13 5- and 6-membered aromatic heterocycles frequently encountered in small drug-like molecules. We demonstrate that the electrostatic properties of these systems are a key determinant for their orientational preferences. The results of this study can be applied in lead optimization for the improvement of stacking interactions, as it provides detailed energy landscapes for a wide range of coplanar heteroaromatic geometries. These energy landscapes can serve as a guide for ring replacement in structure-based drug design.

Protein dynamics is essential to understand protein function and stability, even though is rarely investigated as the origin of loss-of-function due to genetic variations. Here, we use biochemical, biophysical, cell and computational biology tools to study two loss-of-function and cancer-associated polymorphisms (p.R139W and p.P187S) in human NAD(P)H quinone oxidoreductase 1 (NQO1), a FAD-dependent enzyme which activates cancer pro-drugs and stabilizes several oncosuppressors. We show that p.P187S strongly destabilizes the NQO1 dimer in vitro and increases the flexibility of the C-terminal domain, while a combination of FAD and the inhibitor dicoumarol overcome these alterations. Additionally, changes in global stability due to polymorphisms and ligand binding are linked to the dynamics of the dimer interface, whereas the low activity and affinity for FAD in p.P187S is caused by increased fluctuations at the FAD binding site. Importantly, NQO1 steady-state protein levels in cell cultures correlate primarily with the dynamics of the C-terminal domain, supporting a directional preference in NQO1 proteasomal degradation and the use of ligands binding to this domain to stabilize p.P187S in vivo. In conclusion, protein dynamics are fundamental to understanding loss-of-function in p.P187S, and to develop new pharmacological therapies to rescue this function.The study of NAD(P)H quinone oxidoreductase 1 (NQO1, EC 1.6.5.2) polymorphisms is particularly interesting due to its enhanced expression in several types of cancer 1-3 . NQO1 is a dimeric, two-domain FAD-dependent enzyme (Fig. 1A) which catalyses the two electron reduction of quinones and related substrates through an enzyme-substituted mechanism in which NAD(P)H enters the active site, reduces the FAD and exits as the oxidised form, allowing the subsequent substrate binding and reduction by the FADH 2 4 . Primarily, NQO1 avoids the formation of reactive semiquinones, maintains antioxidants such as α -tocopherol and ubiquinone in their reduced state, and also activates some anticancer bioreductive drugs (e.g. mitomycin C (MMC) and a MMC analogue, EO9) 4,5 . Additionally, NQO1 interacts with tumour suppressors such as p53 and p73 and stabilizes them towards proteasomal degradation 6-9 , while its interaction with the 20S proteasome prevents the degradation of a plethora of proteins with intrinsically disordered regions, including several cell cycle regulators, tumor suppressors and apoptotic proteins 10 .

BackgroundThe search for intrinsic factors, which account for a protein's capability to act as an allergen, is ongoing. Fold stability has been identified as a molecular feature that affects processing and presentation, thereby influencing an antigen's immunologic properties.ObjectiveWe assessed how changes in fold stability modulate the immunogenicity and sensitization capacity of the major birch pollen allergen Bet v 1.MethodsBy exploiting an exhaustive virtual mutation screening, we generated mutants of the prototype allergen Bet v 1 with enhanced thermal and chemical stability and rigidity. Structural changes were analyzed by means of x-ray crystallography, nuclear magnetic resonance, and molecular dynamics simulations. Stability was monitored by using differential scanning calorimetry, circular dichroism, and Fourier transform infrared spectroscopy. Endolysosomal degradation was simulated in vitro by using the microsomal fraction of JAWS II cells, followed by liquid chromatography coupled to mass spectrometry. Immunologic properties were characterized in vitro by using a human T-cell line specific for the immunodominant epitope of Bet v 1 and in vivo in an adjuvant-free BALB/c mouse model.ResultsFold stabilization of Bet v 1 was pH dependent and resulted in resistance to endosomal degradation at a pH of 5 or greater, affecting presentation of the immunodominant T-cell epitope in vitro. These properties translated in vivo into a strong allergy-promoting TH2-type immune response. Efficient TH2 cell activation required both an increased stability at the pH of the early endosome and efficient degradation at lower pH in the late endosomal/lysosomal compartment.ConclusionsOur data indicate that differential pH-dependent fold stability along endosomal maturation is an essential protein-inherent determinant of allergenicity.

Bromodomain containing proteins PB1, SMARCA4, and SMARCA2 are important components of SWI/SNF chromatin remodeling complexes. We identified bromodomain inhibitors that target these proteins and display unusual binding modes involving water displacement from the KAc binding site. The best compound binds the fifth bromodomain of PB1 with a KD of 124 nM, SMARCA2B and SMARCA4 with KD values of 262 and 417 nM, respectively, and displays excellent selectivity over bromodomains other than PB1, SMARCA2, and SMARCA4.

Background: Constitutive activation of STAT3 is associated with the progression of hepatocellular carcinoma (HCC), and abrogation of STAT3 signaling is a potential target for HCC treatment. Results: A novel azaspirane modulates the JAK-STAT pathway in HCC. Conclusion:The lead compound induces apoptosis by down-regulating STAT3 signaling. Significance: This investigation reports a novel inhibitor of the JAK-STAT pathway with the potential to target various cancers.

Cholesterol-phospholipid bilayers continue to be the current state of the art in membrane models and serve as representative systems for studying the effect of cholesterol on the cell membrane. As the mixing of different lipid species requires long spatio-temporal scales, coarse-grained models have gained increasing popularity in modeling such membrane systems. In this paper, a systematic study of the MARTINI coarse-grained model for the DPPC-cholesterol binary system has been performed. We construct the phase diagram of DPPC lipid bilayers in the presence of different cholesterol concentrations and at different temperatures using coarse-grained Molecular Dynamics (MD) simulations with the MARTINI force field. The phase diagram based on the condensation effect is directly comparable to available experimental data and demonstrates qualitative agreement over all cholesterol concentrations. Self-assembled bilayers quantitatively reproduce experimental observables, such as lateral diffusion of lipids, electron density, area per lipid and lipid order parameters. The phase diagram of the DPPC-cholesterol binary system also reveals the profound effect of cholesterol on the physical properties of phospholipid bilayers such lipid order, diffusion, and fluidity. Cholesterol induces the liquid-ordered phase, which increases the fluidity of the phospholipid hydrocarbon chains above the gel to liquid-crystalline phase transition temperature and decreases it below the phase transition. The present study suggests that the MARTINI force field can be successfully used to obtain molecular level insights into cholesterol-DPPC model membranes.

Human proteins are vulnerable towards disease-associated single amino acid replacements affecting protein stability and function. Interestingly, a few studies have shown that consensus amino acids from mammals or vertebrates can enhance protein stability when incorporated into human proteins. Here, we investigate yet unexplored relationships between the high vulnerability of human proteins towards disease-associated inactivation and recent evolutionary site-specific divergence of stabilizing amino acids. Using phylogenetic, structural and experimental analyses, we show that divergence from the consensus amino acids at several sites during mammalian evolution has caused local protein destabilization in two human proteins linked to disease: cancer-associated NQO1 and alanine:glyoxylate aminotransferase, mutated in primary hyperoxaluria type I. We demonstrate that a single consensus mutation (H80R) acts as a disease suppressor on the most common cancer-associated polymorphism in NQO1 (P187S). The H80R mutation reactivates P187S by enhancing FAD binding affinity through local and dynamic stabilization of its binding site. Furthermore, we show how a second suppressor mutation (E247Q) cooperates with H80R in protecting the P187S polymorphism towards inactivation through long-range allosteric communication within the structural ensemble of the protein. Our results support that recent divergence of consensus amino acids may have occurred with neutral effects on many functional and regulatory traits of wild-type human proteins. However, divergence at certain sites may have increased the propensity of some human proteins towards inactivation due to disease-associated mutations and polymorphisms. Consensus mutations also emerge as a potential strategy to identify structural hot-spots in proteins as targets for pharmacological rescue in loss-of-function genetic diseases.

The epidermal growth factor receptor (EGFR), when carrying an activating mutation like del19 or L858R, acts as an oncogenic driver in a subset of lung tumors. While tumor responses to tyrosine kinase inhibitors (TKIs) are accompanied by marked tumor shrinkage, the response is usually not durable. Most patients relapse within two years of therapy often due to acquisition of an additional mutation in EGFR kinase domain that confers resistance to TKIs. Crucially, oncogenic EGFR harboring both resistance mutations, T790M and C797S, can no longer be inhibited by currently approved EGFR TKIs. Here, we describe the discovery of BI-4020, which is a noncovalent, wild-type EGFR sparing, macrocyclic TKI. BI-4020 potently inhibits the above-described EGFR variants and induces tumor regressions in a cross-resistant EGFRdel19 T790M C797S xenograft model. Key was the identification of a highly selective but moderately potent benzimidazole followed by complete rigidification of the molecule through macrocyclization.