Search

Gap junction channels (GJCs) mediate intercellular communication by connecting two neighboring cells and enabling direct exchange of ions and small molecules. Cell coupling via connexin-43 (Cx43) GJCs is important in a wide range of cellular processes in health and disease1-3, yet the structural basis of Cx43 function and regulation has not been determined until now. Here we describe the structure of a human Cx43 GJC solved by cryo-EM and single particle analysis at 2.26 Å resolution. The pore region of Cx43 GJC features several lipid-like densities per Cx43 monomer, located close to a putative lateral access site at the monomer boundary. We found a previously undescribed conformation on the cytosolic side of the pore, formed by the N-terminal domain and the transmembrane helix 2 of Cx43 and stabilized by a small molecule. Structures of the Cx43 GJC and hemichannels in nanodiscs reveal a similar gate arrangement. The features of the Cx43 GJC and hemichannel cryo-EM maps and the channel properties revealed by molecular dynamics simulations suggest that the captured states of Cx43 are consistent with a closed state.

The cationic porphyrin TMPyP4 is a well-established DNA G-quadruplex (G4) binding ligand that can stabilize different topologies via multiple binding modes. However, TMPyP4 can have both a stabilizing and destabilizing effect on RNA G4 structures. The structural mechanisms that mediate RNA G4 unfolding remain unknown. Here, we report on the TMPyP4-induced RNA G4 unfolding mechanism studied by well-tempered metadynamics (WT-MetaD) with supporting biophysical experiments. The simulations predict a two-state mechanism of TMPyP4 interaction via a groove-bound and a top-face-bound conformation. The dynamics of TMPyP4 stacking on the top tetrad disrupts Hoogsteen H-bonds between guanine bases, resulting in the consecutive TMPyP4 intercalation from top-to-bottom G-tetrads. The results reveal a striking correlation between computational and experimental approaches and validate WT-MetaD simulations as a powerful tool for studying RNA G4-ligand interactions.

[Image: see text] Whether recent updates and new releases of atomistic force fields can model the structural and dynamical properties of proteins containing both folded and partially disordered domains is still unclear. To address this fundamental question, we tested eight recently released force fields against our set of nuclear magnetic resonance (NMR) observables for a complex and medically relevant system, the major factor VIII binding region on the von Willebrand factor. This biomedically important region comprises both a folded and a partially structured domain. By using an enhanced sampling technique (temperature replica-exchange molecular dynamics simulations), we find that some force fields indeed rise to the challenge and capture the structural and dynamical features of the NMR ensemble and, therefore, are the appropriate choice for simulations of proteins with partially structured domains. What is more, we show that only such force fields can qualitatively capture the effects of a pathogenic mutation on the structural ensemble.

The rise of multi-drug resistance in bacterial pathogens is one of the grand challenges facing medical science. A major concern is the speed of development of β-lactamase-mediated resistance in Gram-negative species, thus putting at risk the efficacy of the most recently approved antibiotics and inhibitors, including carbapenems and avibactam, respectively. New strategies to overcome resistance are urgently required, which will ultimately be facilitated by a deeper understanding of the mechanisms that regulate the function of β-lactamases such as the Klebsiella Pneumoniae carbapenemases (KPCs). Using enhanced sampling computational methods together with site-directed mutagenesis, we report the identification of two “hydrophobic networks” in the KPC-2 enzyme, the integrity of which has been found to be essential for protein stability and corresponding resistance. Present throughout the structure, these networks are responsible for the structural integrity and allosteric signaling. Disruption of the networks leads to a loss of the KPC-2 mediated resistance phenotype, resulting in restored susceptibility to different classes of β-lactam antibiotics including carbapenems and cephalosporins. The ”hydrophobic networks” were found to be highly conserved among class-A β-lactamases, which implies their suitability for exploitation as a potential target for therapeutic intervention.

In all living organisms, coenzyme A (CoA) is an essential cofactor with a unique design allowing it to function as an acyl group carrier and a carbonyl-activating group in diverse biochemical reactions. It is synthesized in a highly conserved process in prokaryotes and eukaryotes that requires pantothenic acid (vitamin B5), cysteine and ATP. CoA and its thioester derivatives are involved in major metabolic pathways, allosteric interactions and the regulation of gene expression. A novel unconventional function of CoA in redox regulation has been recently discovered in mammalian cells and termed protein CoAlation. Here, we report for the first time that protein CoAlation occurs at a background level in exponentially growing bacteria and is strongly induced in response to oxidizing agents and metabolic stress. Over 12% of Staphylococcus aureus gene products were shown to be CoAlated in response to diamide-induced stress. In vitro CoAlation of S. aureus glyceraldehyde-3-phosphate dehydrogenase was found to inhibit its enzymatic activity and to protect the catalytic cysteine 151 from overoxidation by hydrogen peroxide. These findings suggest that in exponentially growing bacteria, CoA functions to generate metabolically active thioesters, while it also has the potential to act as a low-molecular-weight antioxidant in response to oxidative and metabolic stress.

p38α is a Ser/Thr protein kinase involved in a variety of cellular processes and pathological conditions, which makes it a promising pharmacological target. Although the activity of the enzyme is highly regulated, its molecular mechanism of activation remains largely unexplained, even after decades of research. By using state-of-the-art molecular dynamics simulations, we decipher the key elements of the complex molecular mechanism refined by evolution to allow for a fine tuning of p38α kinase activity. Our study describes for the first time the molecular effects of different regulators of the enzymatic activity, and provides an integrative picture of the activation mechanism that explains the seemingly contradictory X-ray and NMR data. AK was partly supported by the European Union Seventh Framework Programme (FP7/Marie Curie Actions/COFUND, no. 600404). FLG and GS acknowledge financial support from EPSRC (grant EP/M013898/1). ARN was supported by a grant from the European Commission (ERC, project ID: 294665). MO is an ICREA research fellow. MO acknowledges the support from two grants from the European Commission H2020 program (BioExcel project ID: 676559 and ELIXIR- Excellerate project ID: 675728) and the Spanish MINECO (BIO2015/64802-R).

Cryptic pockets, that is, sites on protein targets that only become apparent when drugs bind, provide a promising alternative to classical binding sites for drug development. Here, we investigate the nature and dynamical properties of cryptic sites in four pharmacologically relevant targets, while comparing the efficacy of various simulation-based approaches in discovering them. We find that the studied cryptic sites do not correspond to local minima in the computed conformational free energy landscape of the unliganded proteins. They thus promptly close in all of the molecular dynamics simulations performed, irrespective of the force-field used. Temperature-based enhanced sampling approaches, such as Parallel Tempering, do not improve the situation, as the entropic term does not help in the opening of the sites. The use of fragment probes helps, as in long simulations occasionally it leads to the opening and binding to the cryptic sites. Our observed mechanism of cryptic site formation is suggestive of an interplay between two classical mechanisms: induced-fit and conformational selection. Employing this insight, we developed a novel Hamiltonian Replica Exchange-based method "SWISH" (Sampling Water Interfaces through Scaled Hamiltonians), which combined with probes resulted in a promising general approach for cryptic site discovery. We also addressed the issue of "false-positives" and propose a simple approach to distinguish them from druggable cryptic pockets. Our simulations, whose cumulative sampling time was more than 200 μs, help in clarifying the molecular mechanism of pocket formation, providing a solid basis for the choice of an efficient computational method.

Author(s): Fechner, Uli; de Graaf, Chris; Torda, Andrew E; Gussregen, Stefan; Evers, Andreas; Matter, Hans; Hessler, Gerhard; Richmond, Nicola J; Schmidtke, Peter; Segler, Marwin HS; Waller, Mark P; Pleik, Stefanie; Shea, Joan-Emma; Levine, Zachary; Mullen, Ryan; van den Broek, Karina; Epple, Matthias; Kuhn, Hubert; Truszkowski, Andreas; Zielesny, Achim; Fraaije, Johannes Hans; Gracia, Ruben Serral; Kast, Stefan M; Bulusu, Krishna C; Bender, Andreas; Yosipof, Abraham; Nahum, Oren; Senderowitz, Hanoch; Krotzky, Timo; Schulz, Robert; Wolber, Gerhard; Bietz, Stefan; Rarey, Matthias; Zimmermann, Markus O; Lange, Andreas; Ruff, Manuel; Heidrich, Johannes; Onlia, Ionut; Exner, Thomas E; Boeckler, Frank M; Bermudez, Marcel; Firaha, Dzmitry S; Holloczki, Oldamur; Kirchner, Barbara; Tautermann, Christofer S; Volkamer, Andrea; Eid, Sameh; Turk, Samo; Rippmann, Friedrich; Fulle, Simone; Saleh, Noureldin; Saladino, Giorgio; Gervasio, Francesco L; Haensele, Elke; Banting, Lee; Whitley, David C; Oliveira Santos, Jana Sopkova-de; Bureau, Ronan; Clark, Timothy; Sandmann, Achim; Lanig, Harald; Kibies, Patrick; Heil, Jochen; Hoffgaard, Franziska; Frach, Roland; Engel, Julian; Smith, Steven; Basu, Debjit; Rauh, Daniel; Kohlbacher, Oliver; Boeckler, Frank M; Essex, Jonathan W; Bodnarchuk, Michael S; Ross, Gregory A; Finkelmann, Arndt R; Goller, Andreas H; Schneider, Gisbert; Husch, Tamara; Schutter, Christoph; Balducci, Andrea; Korth, Martin; Ntie-Kang, Fidele; Gunther, Stefan; Sippl, Wolfgang; Mbaze, Luc Meva'a

Fibroblast growth factor receptors (FGFRs) are recognized therapeutic targets in cancer. We here describe insights underpinning the impact of mutations on FGFR1 and FGFR3 kinase activity and drug efficacy, using a combination of computational calculations and experimental approaches including cellular studies, X-ray crystallography and biophysical and biochemical measurements. Our findings reveal that some of the tested compounds, in particular TKI258, could provide therapeutic opportunity not only for patients with primary alterations in FGFR but also for acquired resistance due to the gatekeeper mutation. The accuracy of the computational methodologies applied here shows a potential for their wider application in studies of drug binding and in assessments of functional and mechanistic impacts of mutations, thus assisting efforts in precision medicine., Highlights • Point mutations in FGFR can result in distinct signatures for drug efficacy. • Gatekeeper mutations in FGFR1 and FGFR3 retain sensitivity to TKI258. • Theoretical assessments of drug binding and allosteric effects show high level of accuracy.

Due to its inhibition of the Abl kinase domain in the BCR-ABL fusion protein, imatinib is strikingly effective in the initial stage of chronic myeloid leukemia with more than 90% of the patients showing complete remission. However, as in the case of most targeted anti-cancer therapies, the emergence of drug resistance is a serious concern. Several drug-resistant mutations affecting the catalytic domain of Abl and other tyrosine kinases are now known. But, despite their importance and the adverse effect that they have on the prognosis of the cancer patients harboring them, the molecular mechanism of these mutations is still debated. Here by using long molecular dynamics simulations and large-scale free energy calculations complemented by in vitro mutagenesis and microcalorimetry experiments, we model the effect of several widespread drug-resistant mutations of Abl. By comparing the conformational free energy landscape of the mutants with those of the wild-type tyrosine kinases we clarify their mode of action. It involves significant and complex changes in the inactive-to-active dynamics and entropy/enthalpy balance of two functional elements: the activation-loop and the conserved DFG motif. What is more the T315I gatekeeper mutant has a significant impact on the binding mechanism itself and on the binding kinetics., Author Summary Imatinib remains the most important and studied anti-cancer drug for cancer therapy in its new paradigm. Due to its inhibition of the Abl kinase domain, imatinib is strikingly effective in the initial stage of chronic myeloid leukemia with more than 90% of the patients showing a complete remission. However, the emergence of drug resistance is a serious concern. Here, we investigate the molecular mechanism of drug-resistant mutations which, despite the importance and the adverse effect on cancer patients prognosis, is still debated. Our extensive molecular simulations and free energy calculations are consistent with an allosteric effect of the single-point drug-resistance-causing mutations on the conformational dynamics. Two partially independent conformational changes play a role. Our findings might help the design of anti-cancer therapies to overcome drug resistance and be used to predict the clinical relevance of new drug-resistant mutants found by genetic screenings of tumor samples.

Regulation of the c-Abl (ABL1) tyrosine kinase is important because of its role in cellular signaling, and its relevance in the leukemiogenic counterpart (BCR-ABL). Both auto-inhibition and full activation of c-Abl are regulated by the interaction of the catalytic domain with the Src Homology 2 (SH2) domain. The mechanism by which this interaction enhances catalysis is not known. We combined computational simulations with mutagenesis and functional analysis to find that the SH2 domain conveys both local and global effects on the dynamics of the catalytic domain. Locally, it regulates the flexibility of the αC helix in a fashion reminiscent of cyclins in cyclin-dependent kinases, reorienting catalytically important motifs. At a more global level, SH2 binding redirects the hinge motion of the N and C lobes and changes the conformational equilibrium of the activation loop. The complex network of subtle structural shifts that link the SH2 domain with the activation loop and the active site may be partially conserved with other SH2-domain containing kinases and therefore offer additional parameters for the design of conformation-specific inhibitors.

Regulation of the c-Abl (ABL1) tyrosine kinase is important because of its role in cellular signaling, and its relevance in the leukemiogenic counterpart (BCR-ABL). Both auto-inhibition and full activation of c-Abl are regulated by the interaction of the catalytic domain with the Src Homology 2 (SH2) domain. The mechanism by which this interaction enhances catalysis is not known. We combined computational simulations with mutagenesis and functional analysis to find that the SH2 domain conveys both local and global effects on the dynamics of the catalytic domain. Locally, it regulates the flexibility of the αC helix in a fashion reminiscent of cyclins in cyclin-dependent kinases, reorienting catalytically important motifs. At a more global level, SH2 binding redirects the hinge motion of the N and C lobes and changes the conformational equilibrium of the activation loop. The complex network of subtle structural shifts that link the SH2 domain with the activation loop and the active site may be partially conserved with other SH2-domain containing kinases and therefore offer additional parameters for the design of conformation-specific inhibitors., Author Summary The Abl kinase is a key player in many crucial cellular processes. It is also an important anti-cancer drug target, because a mutation leading to the fusion protein Bcr-Abl is the main cause for chronic myeloid leukemia (CML). Abl inhibitors are currently the only pharmaceutical treatment for CML. There are two main difficulties associated with the development of kinase inhibitors: the high similarity between active sites of different kinases, which makes selectivity a challenge, and mutations leading to resistance, which make it mandatory to search for alternative drugs. One important factor controlling Abl is the interplay between the catalytic domain and an SH2 domain. We used computer simulations to understand how the interactions between the domains modify the dynamic of the kinase and detected both local and global effects. Based on our computer model, we suggested mutations that should alter the domain-domain interplay. Consequently, we tested the mutants experimentally and found that they support our hypothesis. We propose that our findings can be of help for the development of new classes of Abl inhibitors, which would modify the domain-domain interplay instead of interfering directly with the active site.

Large-scale conformational transitions represent both a challenge and an opportunity for computational drug design. Exploring the conformational space of a druggable target with sufficient detail is computationally demanding. However, if it were possible to fully account for target flexibility, one could exploit this knowledge to rationally design more potent and more selective drug candidates. Here, we discuss how molecular dynamics together with free energy algorithms based on Metadynamics and Path Collective Variables can be used to study both large-scale conformational transitions and ligand binding to flexible targets. We show real-life examples of how these methods have been applied in the case of cyclin-dependent kinases, a family of flexible targets that shows promise in cancer therapy.

Cyclin‐dependent kinases (CDKs) are mostly known for their role in the cell cycle regulation. The activation mechanism of all CDKs involves the association with a regulatory protein, generally a cyclin, that binds to the kinase unit and stabilizes a catalytically active conformation. Active and inactive conformations of CDKs are characterized by the different spatial localization of two typical elements, namely the activation loop and an □‐helix, whose amino‐acid composition varies throughout the family.