1
|
Bisindolylmaleimides New Ligands of CaM Protein. Molecules 2022; 27:molecules27217161. [PMID: 36363988 PMCID: PMC9653884 DOI: 10.3390/molecules27217161] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/11/2022] [Revised: 09/28/2022] [Accepted: 09/28/2022] [Indexed: 11/06/2022] Open
Abstract
In the present study, we reported the interactions at the molecular level of a series of compounds called Bisindolylmaleimide, as potential inhibitors of the calmodulin protein. Bisindolylmaleimide compounds are drug prototypes derived from Staurosporine, an alkaloid with activity for cancer treatment. Bisindolylmaleimide compounds II, IV, VII, X, and XI, are proposed and reported as possible inhibitors of calmodulin protein for the first time. For the above, a biotechnological device was used (fluorescent biosensor hCaM M124C-mBBr) to directly determine binding parameters experimentally (Kd and stoichiometry) of these compounds, and molecular modeling tools (Docking, Molecular Dynamics, and Chemoinformatic Analysis) to carry out the theoretical studies and complement the experimental data. The results indicate that this compound binds to calmodulin with a Kd between 193–248 nM, an order of magnitude lower than most classic inhibitors. On the other hand, the theoretical studies support the experimental results, obtaining an acceptable correlation between the ΔGExperimental and ΔGTheoretical (r2 = 0.703) and providing us with complementary molecular details of the interaction between the calmodulin protein and the Bisindolylmaleimide series. Chemoinformatic analyzes bring certainty to Bisindolylmaleimide compounds to address clinical steps in drug development. Thus, these results make these compounds attractive to be considered as possible prototypes of new calmodulin protein inhibitors.
Collapse
|
2
|
Khan S. Conformational spread drives the evolution of the calcium-calmodulin protein kinase II. Sci Rep 2022; 12:8499. [PMID: 35589775 PMCID: PMC9120016 DOI: 10.1038/s41598-022-12090-y] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/26/2021] [Accepted: 05/03/2022] [Indexed: 11/17/2022] Open
Abstract
The calcium calmodulin (Ca2+/CaM) dependent protein kinase II (CaMKII) decodes Ca2+ frequency oscillations. The CaMKIIα isoform is predominantly expressed in the brain and has a central role in learning. I matched residue and organismal evolution with collective motions deduced from the atomic structure of the human CaMKIIα holoenzyme to learn how its ring architecture abets function. Protein dynamic simulations showed its peripheral kinase domains (KDs) are conformationally coupled via lateral spread along the central hub. The underlying β-sheet motions in the hub or association domain (AD) were deconvolved into dynamic couplings based on mutual information. They mapped onto a coevolved residue network to partition the AD into two distinct sectors. A second, energetically stressed sector was added to ancient bacterial enzyme dimers for assembly of the ringed hub. The continued evolution of the holoenzyme after AD–KD fusion targeted the sector’s ring contacts coupled to the KD. Among isoforms, the α isoform emerged last and, it alone, mutated rapidly after the poikilotherm–homeotherm jump to match the evolution of memory. The correlation between dynamics and evolution of the CaMKII AD argues single residue substitutions fine-tune hub conformational spread. The fine-tuning could increase CaMKIIα Ca2+ frequency response range for complex learning functions.
Collapse
Affiliation(s)
- Shahid Khan
- Molecular Biology Consortium, Lawrence Berkeley National Laboratory, Berkeley, CA, 94720, USA. .,SBA School of Science and Engineering, LUMS, Lahore, Pakistan. .,Laboratory of Cell Biology, NINDS, NIH, Bethesda, MD, 20892, USA.
| |
Collapse
|
3
|
Kazan IC, Sharma P, Rahman MI, Bobkov A, Fromme R, Ghirlanda G, Ozkan SB. Design of novel cyanovirin-N variants by modulation of binding dynamics through distal mutations. eLife 2022; 11:67474. [PMID: 36472898 PMCID: PMC9725752 DOI: 10.7554/elife.67474] [Citation(s) in RCA: 5] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/12/2021] [Accepted: 11/28/2022] [Indexed: 12/07/2022] Open
Abstract
We develop integrated co-evolution and dynamic coupling (ICDC) approach to identify, mutate, and assess distal sites to modulate function. We validate the approach first by analyzing the existing mutational fitness data of TEM-1 β-lactamase and show that allosteric positions co-evolved and dynamically coupled with the active site significantly modulate function. We further apply ICDC approach to identify positions and their mutations that can modulate binding affinity in a lectin, cyanovirin-N (CV-N), that selectively binds to dimannose, and predict binding energies of its variants through Adaptive BP-Dock. Computational and experimental analyses reveal that binding enhancing mutants identified by ICDC impact the dynamics of the binding pocket, and show that rigidification of the binding residues compensates for the entropic cost of binding. This work suggests a mechanism by which distal mutations modulate function through dynamic allostery and provides a blueprint to identify candidates for mutagenesis in order to optimize protein function.
Collapse
Affiliation(s)
- I Can Kazan
- Center for Biological Physics and Department of Physics, Arizona State UniversityTempeUnited States,School of Molecular Sciences, Arizona State UniversityTempeUnited States
| | - Prerna Sharma
- School of Molecular Sciences, Arizona State UniversityTempeUnited States
| | | | - Andrey Bobkov
- Sanford Burnham Prebys Medical Discovery InstituteLa JollaUnited States
| | - Raimund Fromme
- School of Molecular Sciences, Arizona State UniversityTempeUnited States
| | - Giovanna Ghirlanda
- School of Molecular Sciences, Arizona State UniversityTempeUnited States
| | - S Banu Ozkan
- Center for Biological Physics and Department of Physics, Arizona State UniversityTempeUnited States
| |
Collapse
|
4
|
Ullah H, Mahmud S, Hossain MJ, Islam MSB, Kibria KMK. Immunoinformatic identification of the epitope-based vaccine candidates from Maltoporin, FepA and OmpW of Shigella Spp, with molecular docking confirmation. INFECTION GENETICS AND EVOLUTION 2021; 96:105129. [PMID: 34737105 DOI: 10.1016/j.meegid.2021.105129] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 06/01/2021] [Revised: 10/16/2021] [Accepted: 10/28/2021] [Indexed: 11/17/2022]
Abstract
Shigella is a bacterial pathogen that causes shigellosis, fatal bacillary dysentery, responsible for a higher level of mortality worldwide. We adopted a number of computational approaches to predict potential epitope-based vaccine candidates of immunogenic proteins of Shigella spp. We selected three cell surface proteins of the bacterium according to their antigenicity using the VaxiJen server, including, FepA, Maltoporin, and OmpW. The sequence analyses by the IEDB server resulted in three 15-mer peptides of the core epitope, FTAEHTQSV, FLVNQTLTL, and MRAGSATVR from FepA, Maltoporin, and OmpW, respectively, as the most potential epitopes that have an affinity with both cytotoxic and helper T-cells. Moreover, the epitopes showed 73.76%, 99.0%, and 93.07% world population coverage, along with 100% conservancy among the Shigella subspecies. The molecular docking simulation studies were performed to verify the interactions between the peptides and the respective HLAs. Docking analyses showed that the Epitope-MHC complexes had a higher level of global energy score dictating strong binding. We have also predicted B-cell epitopes from the sequences of these three proteins. In vivo study of the proposed epitope might contribute to the development of a functional and efficient vaccine, which might be an effective way to elude dysentery from the world.
Collapse
Affiliation(s)
- Hedayet Ullah
- Department of Biotechnology and Genetic Engineering, Faculty of Life Science, Mawlana Bhashani Science and Technology University, Tangail 1902, Bangladesh
| | - Shahin Mahmud
- Department of Biotechnology and Genetic Engineering, Faculty of Life Science, Mawlana Bhashani Science and Technology University, Tangail 1902, Bangladesh
| | - Md Jakir Hossain
- Department of Biotechnology and Genetic Engineering, Faculty of Life Science, Mawlana Bhashani Science and Technology University, Tangail 1902, Bangladesh
| | - Md Shaid Bin Islam
- Department of Biotechnology and Genetic Engineering, Faculty of Life Science, Mawlana Bhashani Science and Technology University, Tangail 1902, Bangladesh
| | - K M Kaderi Kibria
- Department of Biotechnology and Genetic Engineering, Faculty of Life Science, Mawlana Bhashani Science and Technology University, Tangail 1902, Bangladesh.
| |
Collapse
|
5
|
Nde J, Zhang P, Ezerski JC, Lu W, Knapp K, Wolynes PG, Cheung MS. Coarse-Grained Modeling and Molecular Dynamics Simulations of Ca 2+-Calmodulin. Front Mol Biosci 2021; 8:661322. [PMID: 34504868 PMCID: PMC8421859 DOI: 10.3389/fmolb.2021.661322] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/30/2021] [Accepted: 07/21/2021] [Indexed: 12/21/2022] Open
Abstract
Calmodulin (CaM) is a calcium-binding protein that transduces signals to downstream proteins through target binding upon calcium binding in a time-dependent manner. Understanding the target binding process that tunes CaM’s affinity for the calcium ions (Ca2+), or vice versa, may provide insight into how Ca2+-CaM selects its target binding proteins. However, modeling of Ca2+-CaM in molecular simulations is challenging because of the gross structural changes in its central linker regions while the two lobes are relatively rigid due to tight binding of the Ca2+ to the calcium-binding loops where the loop forms a pentagonal bipyramidal coordination geometry with Ca2+. This feature that underlies the reciprocal relation between Ca2+ binding and target binding of CaM, however, has yet to be considered in the structural modeling. Here, we presented a coarse-grained model based on the Associative memory, Water mediated, Structure, and Energy Model (AWSEM) protein force field, to investigate the salient features of CaM. Particularly, we optimized the force field of CaM and that of Ca2+ ions by using its coordination chemistry in the calcium-binding loops to match with experimental observations. We presented a “community model” of CaM that is capable of sampling various conformations of CaM, incorporating various calcium-binding states, and carrying the memory of binding with various targets, which sets the foundation of the reciprocal relation of target binding and Ca2+ binding in future studies.
Collapse
Affiliation(s)
- Jules Nde
- Department of Physics, University of Houston, Houston, TX, United States.,Center for Theoretical Biological Physics, Rice University, Houston, TX, United States
| | - Pengzhi Zhang
- Department of Physics, University of Houston, Houston, TX, United States
| | - Jacob C Ezerski
- Department of Physics, University of Houston, Houston, TX, United States
| | - Wei Lu
- Center for Theoretical Biological Physics, Rice University, Houston, TX, United States
| | - Kaitlin Knapp
- Center for Theoretical Biological Physics, Rice University, Houston, TX, United States
| | - Peter G Wolynes
- Center for Theoretical Biological Physics, Rice University, Houston, TX, United States
| | - Margaret S Cheung
- Department of Physics, University of Houston, Houston, TX, United States.,Center for Theoretical Biological Physics, Rice University, Houston, TX, United States
| |
Collapse
|
6
|
Serçinoğlu O, Ozbek P. Sequence-structure-function relationships in class I MHC: A local frustration perspective. PLoS One 2020; 15:e0232849. [PMID: 32421728 PMCID: PMC7233585 DOI: 10.1371/journal.pone.0232849] [Citation(s) in RCA: 18] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/18/2019] [Accepted: 04/22/2020] [Indexed: 12/22/2022] Open
Abstract
Class I Major Histocompatibility Complex (MHC) binds short antigenic peptides with the help of Peptide Loading Complex (PLC), and presents them to T-cell Receptors (TCRs) of cytotoxic T-cells and Killer-cell Immunglobulin-like Receptors (KIRs) of Natural Killer (NK) cells. With more than 10000 alleles, human MHC (Human Leukocyte Antigen, HLA) is the most polymorphic protein in humans. This allelic diversity provides a wide coverage of peptide sequence space, yet does not affect the three-dimensional structure of the complex. Moreover, TCRs mostly interact with HLA in a common diagonal binding mode, and KIR-HLA interaction is allele-dependent. With the aim of establishing a framework for understanding the relationships between polymorphism (sequence), structure (conserved fold) and function (protein interactions) of the human MHC, we performed here a local frustration analysis on pMHC homology models covering 1436 HLA I alleles. An analysis of local frustration profiles indicated that (1) variations in MHC fold are unlikely due to minimally-frustrated and relatively conserved residues within the HLA peptide-binding groove, (2) high frustration patches on HLA helices are either involved in or near interaction sites of MHC with the TCR, KIR, or tapasin of the PLC, and (3) peptide ligands mainly stabilize the F-pocket of HLA binding groove.
Collapse
Affiliation(s)
- Onur Serçinoğlu
- Department of Bioengineering, Recep Tayyip Erdogan University, Faculty of Engineering, Fener, Rize, Turkey
| | - Pemra Ozbek
- Department of Bioengineering, Marmara University, Faculty of Engineering, Goztepe, Istanbul, Turkey
- * E-mail:
| |
Collapse
|
7
|
Molecular mechanism of multispecific recognition of Calmodulin through conformational changes. Proc Natl Acad Sci U S A 2017; 114:E3927-E3934. [PMID: 28461506 DOI: 10.1073/pnas.1615949114] [Citation(s) in RCA: 35] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/18/2022] Open
Abstract
Calmodulin (CaM) is found to have the capability to bind multiple targets. Investigations on the association mechanism of CaM to its targets are crucial for understanding protein-protein binding and recognition. Here, we developed a structure-based model to explore the binding process between CaM and skMLCK binding peptide. We found the cooperation between nonnative electrostatic interaction and nonnative hydrophobic interaction plays an important role in nonspecific recognition between CaM and its target. We also found that the conserved hydrophobic anchors of skMLCK and binding patches of CaM are crucial for the transition from high affinity to high specificity. Furthermore, this association process involves simultaneously both local conformational change of CaM and global conformational changes of the skMLCK binding peptide. We found a landscape with a mixture of the atypical "induced fit," the atypical "conformational selection," and "simultaneously binding-folding," depending on the synchronization of folding and binding. Finally, we extend our discussions on multispecific binding between CaM and its targets. These association characteristics proposed for CaM and skMLCK can provide insights into multispecific binding of CaM.
Collapse
|
8
|
Lipstein N, Göth M, Piotrowski C, Pagel K, Sinz A, Jahn O. Presynaptic Calmodulin targets: lessons from structural proteomics. Expert Rev Proteomics 2017; 14:223-242. [DOI: 10.1080/14789450.2017.1275966] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/20/2022]
Affiliation(s)
- Noa Lipstein
- Department of Molecular Neurobiology, Max-Planck-Institute of Experimental Medicine, Göttingen, Germany
| | - Melanie Göth
- Institute of Chemistry and Biochemistry, Free University Berlin, Berlin & Fritz Haber Institute of the Max-Planck-Society, Berlin, Germany
| | - Christine Piotrowski
- Department of Pharmaceutical Chemistry & Bioanalytics, Institute of Pharmacy, Martin Luther University Halle-Wittenberg, Halle (Saale), Germany
| | - Kevin Pagel
- Institute of Chemistry and Biochemistry, Free University Berlin, Berlin & Fritz Haber Institute of the Max-Planck-Society, Berlin, Germany
| | - Andrea Sinz
- Department of Pharmaceutical Chemistry & Bioanalytics, Institute of Pharmacy, Martin Luther University Halle-Wittenberg, Halle (Saale), Germany
| | - Olaf Jahn
- Proteomics Group, Max-Planck-Institute of Experimental Medicine, Göttingen, Germany
| |
Collapse
|
9
|
Parra RG, Schafer NP, Radusky LG, Tsai MY, Guzovsky AB, Wolynes PG, Ferreiro DU. Protein Frustratometer 2: a tool to localize energetic frustration in protein molecules, now with electrostatics. Nucleic Acids Res 2016; 44:W356-60. [PMID: 27131359 PMCID: PMC4987889 DOI: 10.1093/nar/gkw304] [Citation(s) in RCA: 148] [Impact Index Per Article: 18.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/05/2016] [Accepted: 04/12/2016] [Indexed: 11/14/2022] Open
Abstract
The protein frustratometer is an energy landscape theory-inspired algorithm that aims at localizing and quantifying the energetic frustration present in protein molecules. Frustration is a useful concept for analyzing proteins’ biological behavior. It compares the energy distributions of the native state with respect to structural decoys. The network of minimally frustrated interactions encompasses the folding core of the molecule. Sites of high local frustration often correlate with functional regions such as binding sites and regions involved in allosteric transitions. We present here an upgraded version of a webserver that measures local frustration. The new implementation that allows the inclusion of electrostatic energy terms, important to the interactions with nucleic acids, is significantly faster than the previous version enabling the analysis of large macromolecular complexes within a user-friendly interface. The webserver is freely available at URL: http://frustratometer.qb.fcen.uba.ar.
Collapse
Affiliation(s)
- R Gonzalo Parra
- Protein Physiology Lab, Dep de Química Biológica, Facultad de Ciencias Exactas y Naturales, UBA-CONICET-IQUIBICEN, Buenos Aires, Argentina
| | - Nicholas P Schafer
- Interdisciplinary Nanoscience Center, Department of Molecular Biology and Genetics, Aarhus University, DK-8000 Aarhus, Denmark
| | - Leandro G Radusky
- Structural Bioinformatics Group, Dep de Química Biológica, Facultad de Ciencias Exactas y Naturales, UBA-CONICET-IQUIBICEN, Buenos Aires, Argentina
| | - Min-Yeh Tsai
- Center for Theoretical Biological Physics and Department of Chemistry, Rice University, Houston, TX 77005, USA
| | - A Brenda Guzovsky
- Protein Physiology Lab, Dep de Química Biológica, Facultad de Ciencias Exactas y Naturales, UBA-CONICET-IQUIBICEN, Buenos Aires, Argentina
| | - Peter G Wolynes
- Center for Theoretical Biological Physics and Department of Chemistry, Rice University, Houston, TX 77005, USA
| | - Diego U Ferreiro
- Protein Physiology Lab, Dep de Química Biológica, Facultad de Ciencias Exactas y Naturales, UBA-CONICET-IQUIBICEN, Buenos Aires, Argentina
| |
Collapse
|
10
|
Using the folding landscapes of proteins to understand protein function. Curr Opin Struct Biol 2016; 36:67-74. [PMID: 26812092 DOI: 10.1016/j.sbi.2016.01.001] [Citation(s) in RCA: 31] [Impact Index Per Article: 3.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/19/2015] [Revised: 12/31/2015] [Accepted: 01/06/2016] [Indexed: 11/20/2022]
Abstract
Proteins fold on a biologically-relevant timescale because of a funnel-shaped energy landscape. This landscape is sculpted through evolution by selecting amino-acid sequences that stabilize native interactions while suppressing stable non-native interactions that occur during folding. However, there is strong evolutionary selection for functional residues and these cannot be chosen to optimize folding. Their presence impacts the folding energy landscape in a variety of ways. Here, we survey the effects of functional residues on folding by providing several examples. We then review how such effects can be detected computationally and be used as assays for protein function. Overall, an understanding of how functional residues modulate folding should provide insights into the design of natural proteins and their homeostasis.
Collapse
|