1
|
Patel A, Sinha S, Arantes P, Palermo G. Unveiling Cas8 Dynamics and Regulation within a transposon-encoded Cascade-TniQ Complex. BIORXIV : THE PREPRINT SERVER FOR BIOLOGY 2024:2024.06.21.600075. [PMID: 38948825 PMCID: PMC11213026 DOI: 10.1101/2024.06.21.600075] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 07/02/2024]
Abstract
Cascade is a class 1, type 1 CRISPR-Cas system with a variety of roles in prokaryote defense, specifically against DNA-based viruses. The Vibrio Cholerae transposon, Tn6677, encodes a variant of the type 1F Cascade known as type 1F-3. This Cascade variant complexes with a homodimer of the transposition protein TniQ and leverages the sequence specificity of Cascade to direct the integration activity of the heteromeric transposase tnsA/B, resulting in site-specific transposition of Tn6677. We desire to uncover the molecular details behind R Loop formation of 'Cascade-TniQ.' Due to the lack of a complete model of Cascade-TniQ available at atom-level resolution, we first build a complete model using AlphaFold V2.1. We then simulate this model via classical molecular dynamics and umbrella sampling to study an important regulatory component within Cascade-TniQ, known as the Cas8 'bundle.' Particularly, we show that this alpha helical bundle experiences a free energy barrier to its large-scale translatory motions and relative free energies of its states primarily dependent on a loop within a Cas7 subunit in Cascade-TniQ. Further, we comment on additional structural and dynamical regulatory points of Cascade-TniQ during R Loop formation, such as Cascade-TniQ backbone rigidity, and the potential role TniQ plays in regulating bundle dynamics. In summary, our outcomes provide the first all-atom dynamic representation of one of the largest CRISPR systems, with information that can contribute to understanding the mechanism of nucleic acid binding and, eventually, to transposase recruitment itself. Such information may prove informative to advance genome engineering efforts.
Collapse
|
2
|
Coppola F, Cimino P, Petrone A, Rega N. Evidence of Excited-State Vibrational Mode Governing the Photorelaxation of a Charge-Transfer Complex. J Phys Chem A 2024; 128:1620-1633. [PMID: 38381887 DOI: 10.1021/acs.jpca.3c08366] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/23/2024]
Abstract
Modern, nonlinear, time-resolved spectroscopic techniques have opened new doors for investigating the intriguing but complex world of photoinduced ultrafast out-of-equilibrium phenomena and charge dynamics. The interaction between light and matter introduces an additional dimension, where the complex interplay between electronic and vibrational dynamics needs the most advanced theoretical-computational protocols to be fully understood on the molecular scale. In this study, we showcase the capabilities of ab initio molecular dynamics simulation integrated with a multiresolution wavelet protocol to carefully investigate the excited-state relaxation dynamics in a noncovalent complex involving tetramethylbenzene (TMB) and tetracyanoquinodimethane (TCNQ) undergoing charge transfer (CT) upon photoexcitation. Our protocol provides an accurate description that facilitates a direct comparison between transient vibrational analysis and time-resolved spectroscopic signals. This molecular level perspective enhances our understanding of photorelaxation processes confined in the adiabatic regime and offers an improved interpretation of vibrational spectra. Furthermore, it enables the quantification of anharmonic vibrational couplings between high- and low-frequency modes, specifically the TCNQ "rocking" and "bending" modes. Additionally, it identifies the primary vibrational mode that governs the adiabaticity between the ground state and the CT state. This comprehensive understanding of photorelaxation processes holds significant importance in the rational design and precise control of more efficient photovoltaic and sensor devices.
Collapse
Affiliation(s)
- Federico Coppola
- Scuola Superiore Meridionale, Largo San Marcellino 10, I-80138 Napoli, Italy
| | - Paola Cimino
- Department of Chemical Sciences, University of Napoli Federico II, Complesso Universitario di M.S. Angelo, 80126 Napoli, Italy
| | - Alessio Petrone
- Scuola Superiore Meridionale, Largo San Marcellino 10, I-80138 Napoli, Italy
- Department of Chemical Sciences, University of Napoli Federico II, Complesso Universitario di M.S. Angelo, 80126 Napoli, Italy
- Istituto Nazionale Di Fisica Nucleare, sezione di Napoli, Complesso Universitario di Monte S. Angelo ed. 6, 80126 Napoli, Italia
| | - Nadia Rega
- Scuola Superiore Meridionale, Largo San Marcellino 10, I-80138 Napoli, Italy
- Department of Chemical Sciences, University of Napoli Federico II, Complesso Universitario di M.S. Angelo, 80126 Napoli, Italy
- Istituto Nazionale Di Fisica Nucleare, sezione di Napoli, Complesso Universitario di Monte S. Angelo ed. 6, 80126 Napoli, Italia
| |
Collapse
|
3
|
Multi-ligand molecular docking, simulation, free energy calculations and wavelet analysis of the synergistic effects between natural compounds baicalein and cubebin for the inhibition of the main protease of SARS-CoV-2. J Mol Liq 2023; 374:121253. [PMID: 36694691 PMCID: PMC9854241 DOI: 10.1016/j.molliq.2023.121253] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/11/2022] [Revised: 12/16/2022] [Accepted: 01/10/2023] [Indexed: 01/18/2023]
Abstract
Combination drugs have been used for several diseases for many years since they produce better therapeutic effects. However, it is still a challenge to discover candidates to form a combination drug. This study aimed to investigate whether using a comprehensive in silico approach to identify novel combination drugs from a Chinese herbal formula is an appropriate and creative strategy. We, therefore, used Toujie Quwen Granules for the main protease (Mpro) of SARS-CoV-2 as an example. We first used molecular docking to identify molecular components of the formula which may inhibit Mpro. Baicalein (HQA004) is the most favorable inhibitory ligand. We also identified a ligand from the other component, cubebin (CHA008), which may act to support the proposed HQA004 inhibitor. Molecular dynamics simulations were then performed to further elucidate the possible mechanism of inhibition by HQA004 and synergistic bioactivity conferred by CHA008. HQA004 bound strongly at the active site and that CHA008 enhanced the contacts between HQA004 and Mpro. However, CHA008 also dynamically interacted at multiple sites, and continued to enhance the stability of HQA004 despite diffusion to a distant site. We proposed that HQA004 acted as a possible inhibitor, and CHA008 served to enhance its effects via allosteric effects at two sites. Additionally, our novel wavelet analysis showed that as a result of CHA008 binding, the dynamics and structure of Mpro were observed to have more subtle changes, demonstrating that the inter-residue contacts within Mpro were disrupted by the synergistic ligand. This work highlighted the molecular mechanism of synergistic effects between different herbs as a result of allosteric crosstalk between two ligands at a protein target, as well as revealed that using the multi-ligand molecular docking, simulation, free energy calculations and wavelet analysis to discover novel combination drugs from a Chinese herbal remedy is an innovative pathway.
Collapse
Key Words
- ADME/T, absorption, distribution, metabolism, excretion and toxicity
- COVID-19
- COVID-19, Coronavirus disease 2019
- Combination drug therapy
- Computer simulation
- Computers molecular
- H-bonds, hydrogen bonds
- LD50, median lethal dose
- MD, molecular dynamics
- MM-PBSA, molecular mechanics Poisson Boltzmann surface area
- Mpro, main protease
- Natural products
- PAINS, Pan-assay interference compounds
- RCO, inter-residue contact order
- RMSF, root-mean-square-fluctuation
- SARS-CoV-2, severe acute respiratory syndrome coronavirus 2
- SMILES, Simplified Molecular-input Line-entry System
- TCMSP, traditional Chinese medicine systems pharmacology database and analysis platform
- TQG, Toujie Quwen Granule
- Virus diseases
Collapse
|
4
|
Wordom update 2: A user-friendly program for the analysis of molecular structures and conformational ensembles. Comput Struct Biotechnol J 2023; 21:1390-1402. [PMID: 36817953 PMCID: PMC9929209 DOI: 10.1016/j.csbj.2023.01.026] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/14/2022] [Revised: 01/19/2023] [Accepted: 01/20/2023] [Indexed: 01/29/2023] Open
Abstract
We present the second update of Wordom, a user-friendly and efficient program for manipulation and analysis of conformational ensembles from molecular simulations. The actual update expands some of the existing modules and adds 21 new modules to the update 1 published in 2011. The new adds can be divided into three sets that: 1) analyze atomic fluctuations and structural communication; 2) explore ion-channel conformational dynamics and ionic translocation; and 3) compute geometrical indices of structural deformation. Set 1 serves to compute correlations of motions, find geometrically stable domains, identify a dynamically invariant core, find changes in domain-domain separation and mutual orientation, perform wavelet analysis of large-scale simulations, process the output of principal component analysis of atomic fluctuations, perform functional mode analysis, infer regions of mechanical rigidity, analyze overall fluctuations, and perform the perturbation response scanning. Set 2 includes modules specific for ion channels, which serve to monitor the pore radius as well as water or ion fluxes, and measure functional collective motions like receptor twisting or tilting angles. Finally, set 3 includes tools to monitor structural deformations by computing angles, perimeter, area, volume, β-sheet curvature, radial distribution function, and center of mass. The ring perception module is also included, helpful to monitor supramolecular self-assemblies. This update places Wordom among the most suitable, complete, user-friendly, and efficient software for the analysis of biomolecular simulations. The source code of Wordom and the relative documentation are available under the GNU general public license at http://wordom.sf.net.
Collapse
|
5
|
Spassov DS, Atanasova M, Doytchinova I. A role of salt bridges in mediating drug potency: A lesson from the N-myristoyltransferase inhibitors. Front Mol Biosci 2023; 9:1066029. [PMID: 36703920 PMCID: PMC9871453 DOI: 10.3389/fmolb.2022.1066029] [Citation(s) in RCA: 7] [Impact Index Per Article: 7.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/10/2022] [Accepted: 12/16/2022] [Indexed: 01/12/2023] Open
Abstract
The salt bridge is the strongest non-covalent interaction in nature and is known to participate in protein folding, protein-protein interactions, and molecular recognition. However, the role of salt bridges in the context of drug design has remained not well understood. Here, we report that a common feature in the mechanism of inhibition of the N-myristoyltransferases (NMT), promising targets for the treatment of protozoan infections and cancer, is the formation of a salt bridge between a positively charged chemical group of the small molecule and the negatively charged C-terminus of the enzyme. Substituting the inhibitor positively charged amine group with a neutral methylene group prevents the formation of the salt bridge and leads to a dramatic activity loss. Molecular dynamics simulations have revealed that salt bridges stabilize the NMT-ligand complexes by functioning as molecular clips that stabilize the conformation of the protein structure. As such, the creation of salt bridges between the ligands and their protein targets may find an application as a valuable tool in rational drug design.
Collapse
|
6
|
Marmesin and Marmelosin Interact with the Heparan Sulfatase-2 Active Site: Potential Mechanism for Phytochemicals from Bael Fruit Extract as Antitumor Therapeutics. OXIDATIVE MEDICINE AND CELLULAR LONGEVITY 2023; 2023:9982194. [PMID: 36644581 PMCID: PMC9836799 DOI: 10.1155/2023/9982194] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 09/16/2022] [Revised: 11/09/2022] [Accepted: 11/10/2022] [Indexed: 01/06/2023]
Abstract
Human heparan sulfatase-2 (HSULF-2) is an oncoprotein overexpressed in the surface of all types of tumor cells and its activity plays a critical role in cancer survival and progression. Our previous studies have shown that bael fruit extract, containing marmesin and marmelosin, inhibits the HSULF-2 activity and kills breast tumor cells, but the mechanism of these processes remains fairly known mainly because the HSULF-2's 3D structure is partially known. Herein, we aimed at providing an in silico molecular mechanism of the inhibition of human HSULF-2 by phytochemicals from bael fruit extract. Pharmacokinetic parameters of the main phytochemicals contained in the bael fruit extract, sequence-based 3D structure of human HSULF-2, and the interaction of bael fruit's phytochemicals with the enzyme active site was modeled, evaluated, and verified. Docking studies revealed marmesin and marmelosin as potential inhibitors with binding score -8.5 and -7.7 Kcal/mol; these results were validated using molecular dynamics simulations, which exhibited higher stability of the protein-ligand complexes. Taking together, with our earlier in vitro data, our computational analyses suggest that marmesin and marmelosin interact at the active site of HSULF-2 providing a potential mechanism for its inhibition and consequent antitumor activity by phytochemicals contained in the bael fruit extract.
Collapse
|
7
|
Exploring 129Xe NMR parameters for structural investigation of biomolecules: relativistic, solvent, and thermal effects. J Mol Model 2022; 28:372. [DOI: 10.1007/s00894-022-05365-8] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/09/2022] [Accepted: 10/24/2022] [Indexed: 11/11/2022]
|
8
|
Gonçalves MA, Gonçalves AS, Franca TCC, Santana MS, da Cunha EFF, Ramalho TC. Improved Protocol for the Selection of Structures from Molecular Dynamics of Organic Systems in Solution: The Value of Investigating Different Wavelet Families. J Chem Theory Comput 2022; 18:5810-5818. [DOI: 10.1021/acs.jctc.2c00593] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Mateus A. Gonçalves
- Chemistry Department, Federal University of Lavras, 37200-000Lavras, MG, Brazil
| | - Arlan S. Gonçalves
- Federal Institute of Education Science and Technology of Espírito Santo, Vila Velha29100-000, Brazil
- Federal Institute of Education Science and Technology of Espírito Santo, Unit Goiabeiras, 29056-264Vitória, Espírito Santo, Brazil
| | - Tanos C. C. Franca
- Laboratory of Molecular Modeling Applied to Chemical and Biological Defense (LMACBD), Military Institute of Engineering (IME), Rio de Janeiro22290-270, Brazil
| | - Mozarte S. Santana
- Chemistry Department, Federal University of Lavras, 37200-000Lavras, MG, Brazil
| | | | - Teodorico C. Ramalho
- Chemistry Department, Federal University of Lavras, 37200-000Lavras, MG, Brazil
- Faculty of Informatics and Management, Center for Basic and Applied Research, University of Hradec Kralove, Hradec Kralove50003, Czech Republic
| |
Collapse
|
9
|
Zhao H, Yu J, Weng G, Yu J, Wang E, Gao J, Liu H, Hou T, Wang Z, Kang Y. Structural view on the role of the TRD loop in regulating DNMT3A activity: a molecular dynamics study. Phys Chem Chem Phys 2022; 24:15791-15801. [PMID: 35758413 DOI: 10.1039/d2cp02031a] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
DNA methyltransferase 3A (DNMT3A) has been regarded as a potential epigenetic target for the development of cancer therapeutics. A number of DNMT3A inhibitors have been reported, but most of them do not have good potency, high selectivity and/or low cytotoxicity. It has been suggested that a non-conserved region around the target recognition domain (TRD) loop is implicated in the DNMT3A activity under the allosteric regulation of the ATRX-DNMT3-DNMT3L (ADD) domain, but the molecular mechanism of the regulation of the TRD loop on the DNMT3A activity needs to be elucidated. In this study, based on the reported crystal structures, the dynamics of the TRD loop in different multimerization with/without the bound guest molecule, namely the ADD domain or the DNA molecule, was investigated using conventional molecular dynamics (MD) and umbrella sampling simulations. The simulation results illustrate that the TRD loop exhibits relatively higher flexibility than the other components in the whole catalytic domain (CD), which could be well stabilized into different local minima through the binding with either the ADD domain or the DNA molecule by forming tight hydrogen-bond and salt-bridge networks involving distinct residues. Moreover, the movement of the TRD loop away from the catalytic loop upon activation could be triggered simply by the detachment of the ADD domain, but not necessarily induced by the ADD domain relocation on the CD. All these dynamic structural details could be a supplement to the previously reported crystal structure, which underlines the importance of the structural flexibility for the critical residues in the TRD loop, arousing more interest in the rational design of novel DNMT3A inhibitors targeting this region.
Collapse
Affiliation(s)
- Hong Zhao
- College of Pharmaceutical Sciences, Zhejiang University, Hangzhou 310058, Zhejiang, China. .,Alibaba-Zhejiang University Joint Research Center of Future Digital Healthcare, China
| | - Jie Yu
- College of Pharmaceutical Sciences, Zhejiang University, Hangzhou 310058, Zhejiang, China.
| | - Gaoqi Weng
- College of Pharmaceutical Sciences, Zhejiang University, Hangzhou 310058, Zhejiang, China.
| | - Jiahui Yu
- College of Pharmaceutical Sciences, Zhejiang University, Hangzhou 310058, Zhejiang, China.
| | - Ercheng Wang
- College of Pharmaceutical Sciences, Zhejiang University, Hangzhou 310058, Zhejiang, China.
| | - Junbo Gao
- College of Pharmaceutical Sciences, Zhejiang University, Hangzhou 310058, Zhejiang, China.
| | - Huanxiang Liu
- Faculty of Applied Science, Macao Polytechnic University, Macao, SAR, China
| | - Tingjun Hou
- College of Pharmaceutical Sciences, Zhejiang University, Hangzhou 310058, Zhejiang, China.
| | - Zhe Wang
- College of Pharmaceutical Sciences, Zhejiang University, Hangzhou 310058, Zhejiang, China.
| | - Yu Kang
- College of Pharmaceutical Sciences, Zhejiang University, Hangzhou 310058, Zhejiang, China.
| |
Collapse
|
10
|
Investigating an efficient and accurate protocol for sampling structures from molecular dynamics simulations: a close look by different wavelet families. Theor Chem Acc 2021. [DOI: 10.1007/s00214-021-02816-y] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/20/2022]
|
11
|
Hassan I, Ferraro F, Imhof P. Effect of the Hydration Shell on the Carbonyl Vibration in the Ala-Leu-Ala-Leu Peptide. Molecules 2021; 26:2148. [PMID: 33917998 PMCID: PMC8068333 DOI: 10.3390/molecules26082148] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/01/2021] [Revised: 03/30/2021] [Accepted: 03/30/2021] [Indexed: 11/16/2022] Open
Abstract
The vibrational spectrum of the Ala-Leu-Ala-Leu peptide in solution, computed from first-principles simulations, shows a prominent band in the amide I region that is assigned to stretching of carbonyl groups. Close inspection reveals combined but slightly different contributions by the three carbonyl groups of the peptide. The shift in their exact vibrational signature is in agreement with the different probabilities of these groups to form hydrogen bonds with the solvent. The central carbonyl group has a hydrogen bond probability intermediate to the other two groups due to interchanges between different hydrogen-bonded states. Analysis of the interaction energies of individual water molecules with that group shows that shifts in its frequency are directly related to the interactions with the water molecules in the first hydration shell. The interaction strength is well correlated with the hydrogen bond distance and hydrogen bond angle, though there is no perfect match, allowing geometrical criteria for hydrogen bonds to be used as long as the sampling is sufficient to consider averages. The hydrogen bond state of a carbonyl group can therefore serve as an indicator of the solvent's effect on the vibrational frequency.
Collapse
Affiliation(s)
- Irtaza Hassan
- Institute for Theoretical Physics, Freie Universtiät Berlin, Arnimallee 14, 14195 Berlin, Germany;
| | - Federica Ferraro
- Computer Chemistry Center, Friedrich-Alexander University (FAU) Erlangen-Nürnberg, Nägelsbachstrasse 25, 91052 Erlangen, Germany;
| | - Petra Imhof
- Institute for Theoretical Physics, Freie Universtiät Berlin, Arnimallee 14, 14195 Berlin, Germany;
- Computer Chemistry Center, Friedrich-Alexander University (FAU) Erlangen-Nürnberg, Nägelsbachstrasse 25, 91052 Erlangen, Germany;
| |
Collapse
|
12
|
Ben Imeddourene A, Esque J, André I. Combining multi-scale modelling methods to decipher molecular motions of a branching sucrase from glycoside-hydrolase family 70. PLoS One 2018; 13:e0201323. [PMID: 30067837 PMCID: PMC6070258 DOI: 10.1371/journal.pone.0201323] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/25/2018] [Accepted: 07/12/2018] [Indexed: 01/23/2023] Open
Abstract
Among α-transglucosylases from Glycoside-Hydrolase family 70, the ΔN123-GB-CD2 enzyme derived from the bifunctional DSR-E from L. citreum NRRL B-1299 is particularly interesting as it was the first described engineered Branching Sucrase, not able to elongate glucan polymers from sucrose substrate. The previously reported overall structural organization of this multi-domain enzyme is an intricate U-shape fold conserved among GH70 enzymes which showed a certain conformational variability of the so-called domain V, assumed to play a role in the control of product structures, in available X-ray structures. Understanding the role of functional dynamics on enzyme reaction and substrate recognition is of utmost interest although it remains a challenge for biophysical methods. By combining long molecular dynamics simulation (1μs) and multiple analyses (NMA, PCA, Morelet Continuous Wavelet Transform and Cross Correlations Dynamics), we investigated here the dynamics of ΔN123-GB-CD2 alone and in interaction with sucrose substrate. Overall, our results provide the detailed picture at atomic level of the hierarchy of motions occurring along different timescales and how they are correlated, in agreement with experimental structural data. In particular, detailed analysis of the different structural domains revealed cooperative dynamic behaviors such as twisting, bending and wobbling through anti- and correlated motions, and also two structural hinge regions, of which one was unreported. Several highly flexible loops surrounding the catalytic pocket were also highlighted, suggesting a potential role in the acceptor promiscuity of ΔN123-GBD-CD2. Normal modes and essential dynamics underlined an interesting two-fold dynamic of the catalytic domain A, pivoting about an axis splitting the catalytic gorge in two parts. The comparison of the conformational free energy landscapes using principal component analysis of the enzyme in absence or in presence of sucrose, also revealed a more harmonic basin when sucrose is bound with a shift population of the bending mode, consistent with the substrate binding event.
Collapse
Affiliation(s)
- Akli Ben Imeddourene
- Laboratoire d’Ingénierie des Systèmes Biologiques et Procédés, LISBP, Université de Toulouse, CNRS, INRA, INSA, Toulouse, France
| | - Jérémy Esque
- Laboratoire d’Ingénierie des Systèmes Biologiques et Procédés, LISBP, Université de Toulouse, CNRS, INRA, INSA, Toulouse, France
| | - Isabelle André
- Laboratoire d’Ingénierie des Systèmes Biologiques et Procédés, LISBP, Université de Toulouse, CNRS, INRA, INSA, Toulouse, France
- * E-mail:
| |
Collapse
|
13
|
Botello-Smith WM, Alsamarah A, Chatterjee P, Xie C, Lacroix JJ, Hao J, Luo Y. Polymodal allosteric regulation of Type 1 Serine/Threonine Kinase Receptors via a conserved electrostatic lock. PLoS Comput Biol 2017; 13:e1005711. [PMID: 28827795 PMCID: PMC5578689 DOI: 10.1371/journal.pcbi.1005711] [Citation(s) in RCA: 15] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/25/2017] [Revised: 08/31/2017] [Accepted: 08/05/2017] [Indexed: 11/19/2022] Open
Abstract
Type 1 Serine/Threonine Kinase Receptors (STKR1) transduce a wide spectrum of biological signals mediated by TGF-β superfamily members. The STKR1 activity is tightly controlled by their regulatory glycine-serine rich (GS) domain adjacent to the kinase domain. Despite decades of studies, it remains unknown how physiological or pathological GS domain modifications are coupled to STKR1 kinase activity. Here, by performing molecular dynamics simulations and free energy calculation of Activin-Like Kinase 2 (ALK2), we found that GS domain phosphorylation, FKBP12 dissociation, and disease mutations all destabilize a D354-R375 salt-bridge, which normally acts as an electrostatic lock to prevent coordination of adenosine triphosphate (ATP) to the catalytic site. We developed a WAFEX-guided principal analysis and unraveled how phosphorylation destabilizes this highly conserved salt-bridge in temporal and physical space. Using current-flow betweenness scores, we identified an allosteric network of residue-residue contacts between the GS domain and the catalytic site that controls the formation and disruption of this salt bridge. Importantly, our novel network analysis approach revealed how certain disease-causing mutations bypass FKBP12-mediated kinase inhibition to produce leaky signaling. We further provide experimental evidence that this salt-bridge lock exists in other STKR1s, and acts as a general safety mechanism in STKR1 to prevent pathological leaky signaling. In summary, our study provides a compelling and unifying allosteric activation mechanism in STKR1 kinases that reconciles a large number of experimental studies and sheds light on a novel therapeutic avenue to target disease-related STKR1 mutants. Kinases play central role in essential physiological process and are attractive therapeutic drug targets. One of the important kinase families is Type 1 Serine/Threonine Kinase Receptors (STKR1), which control gene expression in response to extracellular growth factors. The activities of STKR1 are tightly controlled by their regulatory domain, which is distant from the kinase catalytic site. The underlying molecular mechanism is elucidated here. We identified that formation or disruption of a highly conserved charge-charge interaction located near the ATP binding site, mediates the physiological inhibition or activation of STKR1. We find that the stability of this charge-charge interaction is remotely controlled by interactions propagated from the distant regulatory domain. Several disease-causing mutations are located at the regulatory domain. We demonstrate how those mutations bypass these endogenous STKR1 inhibition mechanisms to produce pathological phenotypes. This study provides a general activation mechanism in STKR1 kinases, thus may benefit understanding the molecular mechanism of diseases and drug development.
Collapse
Affiliation(s)
- Wesley M. Botello-Smith
- College of pharmacy, Western University of Health Sciences, Pomona, CA, United States of America
| | - Abdelaziz Alsamarah
- College of pharmacy, Western University of Health Sciences, Pomona, CA, United States of America
| | - Payal Chatterjee
- College of pharmacy, Western University of Health Sciences, Pomona, CA, United States of America
| | - Chen Xie
- College of Veterinary Sciences, Western University of Health Sciences, Pomona, CA, United States of America
| | - Jerome J. Lacroix
- Graduate College of Biomedical Sciences, Western University of Health Sciences, Pomona, CA, United States of America
| | - Jijun Hao
- College of Veterinary Sciences, Western University of Health Sciences, Pomona, CA, United States of America
- * E-mail: (YL); (JH)
| | - Yun Luo
- College of pharmacy, Western University of Health Sciences, Pomona, CA, United States of America
- * E-mail: (YL); (JH)
| |
Collapse
|