1
|
Chaudhuri D, Datta J, Majumder S, Giri K. Repurposing of drug molecules from FDA database against Hepatitis C virus E2 protein using ensemble docking approach. Mol Divers 2024; 28:1175-1188. [PMID: 37061608 DOI: 10.1007/s11030-023-10646-2] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/22/2023] [Accepted: 03/31/2023] [Indexed: 04/17/2023]
Abstract
Hepatitis C virus, a member of the Flaviviridae family and genus Hepacivirus, is an enveloped, positively single stranded RNA virus. Its surface consists of a heterodimer of E1 and E2 proteins which play a crucial role in receptor binding and membrane fusion. In this study we have used in silico virtual screening by utilizing ensemble docking on the approved drugs. These drugs can bind with high efficiency to the 36 prominent conformations of the CD81 binding site clustered from a total of 3 µs MD simulation data on the E2 protein. We started with 9213 compounds from the FDA list of drugs and progressively came down to 5 compounds which have been seen to bind with very high efficiency to not only all the conformations but also the two predicted druggable pockets that encompass the CD81 binding site. MM/PBSA binding energy calculations also point to the highly stable interaction of the compounds to the E2 protein. This study may in future broaden the arsenal of therapeutics for use against HCV infection and lead to more effective care against the virus.
Collapse
Affiliation(s)
- Dwaipayan Chaudhuri
- Department of Life Sciences, Presidency University, 86/1 College Street, Kolkata, 700073, India
| | - Joyeeta Datta
- Department of Life Sciences, Presidency University, 86/1 College Street, Kolkata, 700073, India
| | - Satyabrata Majumder
- Department of Life Sciences, Presidency University, 86/1 College Street, Kolkata, 700073, India
| | - Kalyan Giri
- Department of Life Sciences, Presidency University, 86/1 College Street, Kolkata, 700073, India.
| |
Collapse
|
2
|
Xu SY, Zhou L, Xu Y, Hong HY, Dai C, Wang YJ, Zheng YG. Recent advances in structure-based enzyme engineering for functional reconstruction. Biotechnol Bioeng 2023; 120:3427-3445. [PMID: 37638646 DOI: 10.1002/bit.28540] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/04/2023] [Revised: 07/27/2023] [Accepted: 08/15/2023] [Indexed: 08/29/2023]
Abstract
Structural information can help engineer enzymes. Usually, specific amino acids in particular regions are targeted for functional reconstruction to enhance the catalytic performance, including activity, stereoselectivity, and thermostability. Appropriate selection of target sites is the key to structure-based design, which requires elucidation of the structure-function relationships. Here, we summarize the mutations of residues in different specific regions, including active center, access tunnels, and flexible loops, on fine-tuning the catalytic performance of enzymes, and discuss the effects of altering the local structural environment on the functions. In addition, we keep up with the recent progress of structure-based approaches for enzyme engineering, aiming to provide some guidance on how to take advantage of the structural information.
Collapse
Affiliation(s)
- Shen-Yuan Xu
- Key Laboratory of Bioorganic Synthesis of Zhejiang Province, College of Biotechnology and Bioengineering, Zhejiang University of Technology, Hangzhou, Zhejiang, People's Republic of China
- Engineering Research Center of Bioconversion and Biopurification of the Ministry of Education, Zhejiang University of Technology, Hangzhou, Zhejiang, People's Republic of China
- The National and Local Joint Engineering Research Center for Biomanufacturing of Chiral Chemicals, Zhejiang University of Technology, Hangzhou, Zhejiang, People's Republic of China
| | - Lei Zhou
- Key Laboratory of Bioorganic Synthesis of Zhejiang Province, College of Biotechnology and Bioengineering, Zhejiang University of Technology, Hangzhou, Zhejiang, People's Republic of China
- Engineering Research Center of Bioconversion and Biopurification of the Ministry of Education, Zhejiang University of Technology, Hangzhou, Zhejiang, People's Republic of China
- The National and Local Joint Engineering Research Center for Biomanufacturing of Chiral Chemicals, Zhejiang University of Technology, Hangzhou, Zhejiang, People's Republic of China
| | - Ying Xu
- Key Laboratory of Bioorganic Synthesis of Zhejiang Province, College of Biotechnology and Bioengineering, Zhejiang University of Technology, Hangzhou, Zhejiang, People's Republic of China
- Engineering Research Center of Bioconversion and Biopurification of the Ministry of Education, Zhejiang University of Technology, Hangzhou, Zhejiang, People's Republic of China
- The National and Local Joint Engineering Research Center for Biomanufacturing of Chiral Chemicals, Zhejiang University of Technology, Hangzhou, Zhejiang, People's Republic of China
| | - Han-Yue Hong
- Key Laboratory of Bioorganic Synthesis of Zhejiang Province, College of Biotechnology and Bioengineering, Zhejiang University of Technology, Hangzhou, Zhejiang, People's Republic of China
- Engineering Research Center of Bioconversion and Biopurification of the Ministry of Education, Zhejiang University of Technology, Hangzhou, Zhejiang, People's Republic of China
- The National and Local Joint Engineering Research Center for Biomanufacturing of Chiral Chemicals, Zhejiang University of Technology, Hangzhou, Zhejiang, People's Republic of China
| | - Chen Dai
- Key Laboratory of Bioorganic Synthesis of Zhejiang Province, College of Biotechnology and Bioengineering, Zhejiang University of Technology, Hangzhou, Zhejiang, People's Republic of China
- Engineering Research Center of Bioconversion and Biopurification of the Ministry of Education, Zhejiang University of Technology, Hangzhou, Zhejiang, People's Republic of China
- The National and Local Joint Engineering Research Center for Biomanufacturing of Chiral Chemicals, Zhejiang University of Technology, Hangzhou, Zhejiang, People's Republic of China
| | - Ya-Jun Wang
- Key Laboratory of Bioorganic Synthesis of Zhejiang Province, College of Biotechnology and Bioengineering, Zhejiang University of Technology, Hangzhou, Zhejiang, People's Republic of China
- Engineering Research Center of Bioconversion and Biopurification of the Ministry of Education, Zhejiang University of Technology, Hangzhou, Zhejiang, People's Republic of China
- The National and Local Joint Engineering Research Center for Biomanufacturing of Chiral Chemicals, Zhejiang University of Technology, Hangzhou, Zhejiang, People's Republic of China
| | - Yu-Guo Zheng
- Key Laboratory of Bioorganic Synthesis of Zhejiang Province, College of Biotechnology and Bioengineering, Zhejiang University of Technology, Hangzhou, Zhejiang, People's Republic of China
- Engineering Research Center of Bioconversion and Biopurification of the Ministry of Education, Zhejiang University of Technology, Hangzhou, Zhejiang, People's Republic of China
- The National and Local Joint Engineering Research Center for Biomanufacturing of Chiral Chemicals, Zhejiang University of Technology, Hangzhou, Zhejiang, People's Republic of China
| |
Collapse
|
3
|
Bhattarai A, Priyadharshini A, Emerson IA. Investigating the binding affinity of andrographolide against human SARS-CoV-2 spike receptor-binding domain through docking and molecular dynamics simulations. J Biomol Struct Dyn 2023; 41:13438-13453. [PMID: 36764825 DOI: 10.1080/07391102.2023.2174596] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/19/2022] [Accepted: 01/24/2023] [Indexed: 02/12/2023]
Abstract
SARS-CoV-2 is a positive-sense single-stranded RNA virus that causes a deadly coronavirus disease (COVID-19) in humans. The infection of SARS-CoV-2 in humans involves a viral surface spike glycoprotein containing the receptor-binding domain (RBD). The interactions of SARS-CoV-2 with the host angiotensin-converting enzyme 2 (ACE2) receptor are mediated by RBD. It binds to the host ACE2 and influences viral replication and disease pathogenesis. Therefore, targeting the RBD to prevent SARS-CoV-2 infections is of utmost importance. In this study, we used docking and molecular dynamics simulations to understand the binding effect of andrographolide on the SARS-CoV-2 spike protein. During docking, a strong binding affinity was observed between the ligand and the target receptor protein. MD results demonstrated higher conformational fluctuations in the ligand-free protein compared to the bound form. Several residues in the active sites make conformational rearrangements for the S protein to interact with the ligand. While RBD experiences conformational transition to gain more stability upon binding with the ligand. This binding is strengthened via several non-covalent interactions that make the complex structure more stable with higher binding affinity. Overall findings of the study may shed some valuable insights concerning the development of potential therapeutics in the strategies for COVID-19 prevention.
Collapse
Affiliation(s)
- Anil Bhattarai
- Bioinformatics Programming Laboratory, Department of Biotechnology, School of Bio-Sciences and Technology, Vellore Institute of Technology, VIT, Vellore, Tamil Nadu, India
| | - Annadurai Priyadharshini
- Bioinformatics Programming Laboratory, Department of Biotechnology, School of Bio-Sciences and Technology, Vellore Institute of Technology, VIT, Vellore, Tamil Nadu, India
| | - Isaac Arnold Emerson
- Bioinformatics Programming Laboratory, Department of Biotechnology, School of Bio-Sciences and Technology, Vellore Institute of Technology, VIT, Vellore, Tamil Nadu, India
| |
Collapse
|
4
|
Huang S, Mei H, Lu L, Kuang Z, Heng Y, Xu L, Liang X, Qiu M, Pan X. Conformational transitions of caspase-6 in substrate-induced activation process explored by perturbation-response scanning combined with targeted molecular dynamics. Comput Struct Biotechnol J 2021; 19:4156-4164. [PMID: 34527189 PMCID: PMC8342898 DOI: 10.1016/j.csbj.2021.07.017] [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: 03/18/2021] [Revised: 07/17/2021] [Accepted: 07/23/2021] [Indexed: 11/23/2022] Open
Abstract
Caspase-6 participates in a series of neurodegenerative pathways, and has aroused widespread attentions as a promising molecular target for the treatment of neurodegeneration. Caspase-6 is a homodimer with 6 central-stranded β-sheets and 5 α-helices in each monomer. Previous crystallographic studies suggested that the 60′s, 90′s and 130′s helices of caspase-6 undergo a distinctive conformational transition upon substrate binding. Although the caspase-6 structures in apo and active states have been determined, the conformational transition process between the two states remains poorly understood. In this work, perturbation-response scanning (PRS) combined with targeted molecular dynamics (TMD) simulations was employed to unravel the atomistic mechanism of the dynamic conformational transitions underlying the substrate-induced activation process of caspase-6. The results showed that the conformational transition of caspase-6 from apo to active states is mainly characterized by structural rearrangements of the substrate-binding site as well as the conformational changes of 60′s and 130′s extended helices. The H-bond interactions between L1, 130′s helix and 90′s helix are proved to be key determinant factors for substrate-induced conformational transition. These findings provide valuable insights into the activation mechanism of caspase-6 as well as the molecular design of caspase-6 inhibitors.
Collapse
Affiliation(s)
- Shuheng Huang
- Key Laboratory of Biorheological Science and Technology (Ministry of Education), College of Bioengineering, Chongqing University, Chongqing 400045, China
| | - Hu Mei
- Key Laboratory of Biorheological Science and Technology (Ministry of Education), College of Bioengineering, Chongqing University, Chongqing 400045, China
| | - Laichun Lu
- Key Laboratory of Biorheological Science and Technology (Ministry of Education), College of Bioengineering, Chongqing University, Chongqing 400045, China
| | - Zuyin Kuang
- Key Laboratory of Biorheological Science and Technology (Ministry of Education), College of Bioengineering, Chongqing University, Chongqing 400045, China
| | - Yu Heng
- Key Laboratory of Biorheological Science and Technology (Ministry of Education), College of Bioengineering, Chongqing University, Chongqing 400045, China
| | - Lei Xu
- Key Laboratory of Biorheological Science and Technology (Ministry of Education), College of Bioengineering, Chongqing University, Chongqing 400045, China
| | - Xiaoqi Liang
- Key Laboratory of Biorheological Science and Technology (Ministry of Education), College of Bioengineering, Chongqing University, Chongqing 400045, China
| | - Minyao Qiu
- Key Laboratory of Biorheological Science and Technology (Ministry of Education), College of Bioengineering, Chongqing University, Chongqing 400045, China
| | - Xianchao Pan
- Department of Medicinal Chemistry, College of Pharmacy, Southwest Medical University, Luzhou, Sichuan 646000, China
| |
Collapse
|
5
|
Bhattarai A, Emerson IA. Exploring the conformational dynamics and flexibility of intrinsically disordered HIV-1 Nef protein using molecular dynamic network approaches. 3 Biotech 2021; 11:156. [PMID: 33747706 DOI: 10.1007/s13205-021-02698-8] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/30/2020] [Accepted: 02/19/2021] [Indexed: 10/22/2022] Open
Abstract
Intrinsically disordered proteins represent a class of proteins that lack fixed and well-defined three-dimensional structures in solution. HIV-1 Nef is an intrinsically disordered peripheral membrane protein involved in the replication and pathogenesis of HIV-1 infection. Nef controls expression levels of cell surface CD4 molecules that are essential for adaptive immunity. Despite the lack of fixed and stable structures, Nef physically interacts with the host cellular proteins (AP-1/MHC-I) and modulates intracellular trafficking pathways. Therefore, it is essential to understand how this dynamic conformational flexibility affects Nef structures and function. In this study, we combined all-atom molecular dynamics (MD) simulations and dynamic network approaches to better understand the structure and dynamics of Nef in two different forms, the free unbound and the bound state. Using the MD simulation approach, we show that the intrinsically disordered Nef exhibit a large dynamic field with more atomic fluctuations and lesser thermodynamic stability in the unbound conditions. The conformations of Nef change over time, and this protein remains more compact, folded, and stable in the bound form. The dynamic network analysis revealed regions of the protein capable of modulating the conformational behavior of the disordered Nef. The average betweenness centrality (BC) unveiled residues that are critical for mediating protein-protein interactions. The average shortest path length (L) and the perturbation response scanning exposed residues that are likely to be important in steering protein conformational changes. Overall, the study demonstrates how all-atom MD simulations combined with the dynamic network approach can be used to gain further insights into the structure and dynamics-function relationship of intrinsically disordered HIV-1 Nef. SUPPLEMENTARY INFORMATION The online version contains supplementary material available at 10.1007/s13205-021-02698-8.
Collapse
|
6
|
Abdizadeh H, Jalalypour F, Atilgan AR, Atilgan C. A Coarse-Grained Methodology Identifies Intrinsic Mechanisms That Dissociate Interacting Protein Pairs. Front Mol Biosci 2020; 7:210. [PMID: 33195399 PMCID: PMC7477071 DOI: 10.3389/fmolb.2020.00210] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/24/2020] [Accepted: 08/03/2020] [Indexed: 11/13/2022] Open
Abstract
We address the problem of triggering dissociation events between proteins that have formed a complex. We have collected a set of 25 non-redundant, functionally diverse protein complexes having high-resolution three-dimensional structures in both the unbound and bound forms. We unify elastic network models with perturbation response scanning (PRS) methodology as an efficient approach for predicting residues that have the propensity to trigger dissociation of an interacting protein pair, using the three-dimensional structures of the bound and unbound proteins as input. PRS reveals that while for a group of protein pairs, residues involved in the conformational shifts are confined to regions with large motions, there are others where they originate from parts of the protein unaffected structurally by binding. Strikingly, only a few of the complexes have interface residues responsible for dissociation. We find two main modes of response: In one mode, remote control of dissociation in which disruption of the electrostatic potential distribution along protein surfaces play the major role; in the alternative mode, mechanical control of dissociation by remote residues prevail. In the former, dissociation is triggered by changes in the local environment of the protein, e.g., pH or ionic strength, while in the latter, specific perturbations arriving at the controlling residues, e.g., via binding to a third interacting partner is required for decomplexation. We resolve the observations by relying on an electromechanical coupling model which reduces to the usual elastic network result in the limit of the lack of coupling. We validate the approach by illustrating the biological significance of top residues selected by PRS on select cases where we show that the residues whose perturbation leads to the observed conformational changes correspond to either functionally important or highly conserved residues in the complex.
Collapse
Affiliation(s)
- Haleh Abdizadeh
- Groningen Biomolecular Sciences and Biotechnology Institute, University of Groningen, Groningen, Netherlands
| | - Farzaneh Jalalypour
- Faculty of Engineering and Natural Sciences, Sabanci University, Istanbul, Turkey
| | - Ali Rana Atilgan
- Faculty of Engineering and Natural Sciences, Sabanci University, Istanbul, Turkey
| | - Canan Atilgan
- Faculty of Engineering and Natural Sciences, Sabanci University, Istanbul, Turkey
| |
Collapse
|
7
|
Bhattarai A, Emerson IA. Computational investigations on the dynamic binding effect of molecular tweezer CLR01 toward intrinsically disordered HIV-1 Nef. Biotechnol Appl Biochem 2020; 68:513-530. [PMID: 32447788 DOI: 10.1002/bab.1957] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/09/2019] [Accepted: 05/12/2020] [Indexed: 01/08/2023]
Abstract
Intrinsically disordered proteins (IDPs) are highly flexible molecules that undergo disorder to order transition through their interaction with other molecules. IDPs play a vital role in several biological processes ranging from molecular recognition to several human diseases through the protein-protein interaction. The dynamic flexibility of IDPs and their implications in several human diseases enable these molecules to serve as novel therapeutic targets. However, the challenging task is to develop novel drugs against IDPs because of their lack of stable structures and the nature of high conformational flexibility. In this study, we have calculated the dynamic binding effect of the supramolecular tweezer CLR01 against the intrinsically disordered HIV-1 Nef by employing molecular docking and dynamics simulation approaches. From docking results, we predicted the strong binding affinity of the tweezer with the target residues of Nef. The docking results were further validated from the molecular dynamics simulation studies confirming the conformational stability of Nef upon tweezer binding. These findings provide useful insights into the development of potent inhibitors for targeting Nef protein functions.
Collapse
Affiliation(s)
- Anil Bhattarai
- Bioinformatics Programming Laboratory, Department of Biotechnology, School of Bio-Sciences and Technology, Vellore Institute of Technology, Vellore, India
| | - Isaac Arnold Emerson
- Bioinformatics Programming Laboratory, Department of Biotechnology, School of Bio-Sciences and Technology, Vellore Institute of Technology, Vellore, India
| |
Collapse
|
8
|
Jalalypour F, Sensoy O, Atilgan C. Perturb-Scan-Pull: A Novel Method Facilitating Conformational Transitions in Proteins. J Chem Theory Comput 2020; 16:3825-3841. [PMID: 32324386 DOI: 10.1021/acs.jctc.9b01222] [Citation(s) in RCA: 17] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/31/2022]
Abstract
Conformational transitions in proteins facilitate precise physiological functions. Therefore, it is crucial to understand the mechanisms underlying these processes to modulate protein function. Yet, studying structural and dynamical properties of proteins is notoriously challenging due to the complexity of the underlying potential energy surfaces (PES). We have previously developed the perturbation-response scanning (PRS) method to identify key residues that participate in the communication network responsible for specific conformational transitions. PRS is based on a residue-by-residue scan of the protein to determine the subset of residues/forces which provide the closest conformational change leading to a target conformational state, inasmuch as linear response theory applies to these motions. Here, we develop a novel method to further evaluate if conformational transitions may be triggered on the PES. We aim to study functionally relevant conformational transitions in proteins by using results obtained from PRS and feeding them as inputs to steered molecular dynamics simulations. The success and the transferability of the method are evaluated on three protein systems having different complexities of motion on the PES: calmodulin, adenylate kinase, and bacterial ferric binding protein. We find that the method captures the target conformation, while providing key residues and the optimum paths with relatively low free energy profiles.
Collapse
Affiliation(s)
- Farzaneh Jalalypour
- Faculty of Engineering and Natural Sciences, Sabanci University, 34956, Istanbul, Turkey
| | - Ozge Sensoy
- School of Engineering and Natural Sciences, Istanbul Medipol University, 34810, Istanbul, Turkey
| | - Canan Atilgan
- Faculty of Engineering and Natural Sciences, Sabanci University, 34956, Istanbul, Turkey.,Sabanci University Nanotechnology Research and Application Center, SUNUM, 34956, Istanbul, Turkey
| |
Collapse
|
9
|
Sheik Amamuddy O, Veldman W, Manyumwa C, Khairallah A, Agajanian S, Oluyemi O, Verkhivker GM, Tastan Bishop Ö. Integrated Computational Approaches and Tools forAllosteric Drug Discovery. Int J Mol Sci 2020; 21:E847. [PMID: 32013012 PMCID: PMC7036869 DOI: 10.3390/ijms21030847] [Citation(s) in RCA: 48] [Impact Index Per Article: 12.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/25/2019] [Revised: 01/20/2020] [Accepted: 01/21/2020] [Indexed: 12/16/2022] Open
Abstract
Understanding molecular mechanisms underlying the complexity of allosteric regulationin proteins has attracted considerable attention in drug discovery due to the benefits and versatilityof allosteric modulators in providing desirable selectivity against protein targets while minimizingtoxicity and other side effects. The proliferation of novel computational approaches for predictingligand-protein interactions and binding using dynamic and network-centric perspectives has ledto new insights into allosteric mechanisms and facilitated computer-based discovery of allostericdrugs. Although no absolute method of experimental and in silico allosteric drug/site discoveryexists, current methods are still being improved. As such, the critical analysis and integration ofestablished approaches into robust, reproducible, and customizable computational pipelines withexperimental feedback could make allosteric drug discovery more efficient and reliable. In this article,we review computational approaches for allosteric drug discovery and discuss how these tools can beutilized to develop consensus workflows for in silico identification of allosteric sites and modulatorswith some applications to pathogen resistance and precision medicine. The emerging realization thatallosteric modulators can exploit distinct regulatory mechanisms and can provide access to targetedmodulation of protein activities could open opportunities for probing biological processes and insilico design of drug combinations with improved therapeutic indices and a broad range of activities.
Collapse
Affiliation(s)
- Olivier Sheik Amamuddy
- Research Unit in Bioinformatics (RUBi), Department of Biochemistry and Microbiology, Rhodes University, Grahamstown 6140, South Africa; (O.S.A.); (W.V.); (C.M.); (A.K.)
| | - Wayde Veldman
- Research Unit in Bioinformatics (RUBi), Department of Biochemistry and Microbiology, Rhodes University, Grahamstown 6140, South Africa; (O.S.A.); (W.V.); (C.M.); (A.K.)
| | - Colleen Manyumwa
- Research Unit in Bioinformatics (RUBi), Department of Biochemistry and Microbiology, Rhodes University, Grahamstown 6140, South Africa; (O.S.A.); (W.V.); (C.M.); (A.K.)
| | - Afrah Khairallah
- Research Unit in Bioinformatics (RUBi), Department of Biochemistry and Microbiology, Rhodes University, Grahamstown 6140, South Africa; (O.S.A.); (W.V.); (C.M.); (A.K.)
| | - Steve Agajanian
- Graduate Program in Computational and Data Sciences, Keck Center for Science and Engineering, Schmid College of Science and Technology, Chapman University, One University Drive, Orange, CA 92866, USA; (S.A.); (O.O.)
| | - Odeyemi Oluyemi
- Graduate Program in Computational and Data Sciences, Keck Center for Science and Engineering, Schmid College of Science and Technology, Chapman University, One University Drive, Orange, CA 92866, USA; (S.A.); (O.O.)
| | - Gennady M. Verkhivker
- Graduate Program in Computational and Data Sciences, Keck Center for Science and Engineering, Schmid College of Science and Technology, Chapman University, One University Drive, Orange, CA 92866, USA; (S.A.); (O.O.)
- Department of Biomedical and Pharmaceutical Sciences, Chapman University School of Pharmacy, Irvine, CA 92618, USA
| | - Özlem Tastan Bishop
- Research Unit in Bioinformatics (RUBi), Department of Biochemistry and Microbiology, Rhodes University, Grahamstown 6140, South Africa; (O.S.A.); (W.V.); (C.M.); (A.K.)
| |
Collapse
|
10
|
Ross CJ, Atilgan AR, Tastan Bishop Ö, Atilgan C. Unraveling the Motions behind Enterovirus 71 Uncoating. Biophys J 2019; 114:822-838. [PMID: 29490244 DOI: 10.1016/j.bpj.2017.12.021] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/17/2017] [Revised: 11/10/2017] [Accepted: 12/12/2017] [Indexed: 01/24/2023] Open
Abstract
Enterovirus 71 can be a severe pathogen in small children and immunocompromised adults. Virus uncoating is a critical step in the infection of the host cell; however, the mechanisms that control this process remain poorly understood. We applied normal mode analysis and perturbation response scanning to several complexes of the virus capsid and present a coarse-graining approach to analyze the full capsid. We show that our method offers an alternative to expressing the system as a set of rigid blocks and accounts for the interconnection between nodes within each subunit and protein interfaces across the capsid. In our coarse-grained approach, the modes associated with capsid expansion are captured in the first three nondegenerate modes and correspond to the changes observed in structural studies of the virus. We show that the resolution of the analysis may be modified without losing information on the global motions leading to uncoating. Perturbation response scanning revealed that a protomer cannot serve as a functional unit to explain deformations of the capsid. Instead, we define a pentamer as the minimum functional unit to investigate changes within the capsid. From the modal analysis and perturbation response scanning, we locate a hotspot region surrounding the fivefold axis. The range of the effect of these single, hotspot residues extend to 140 Å. The perturbation of internal capsid residues in this region displayed greatest propensity to capsid expansion, thus indicating the significant role that the RNA genome may play in triggering uncoating.
Collapse
Affiliation(s)
- Caroline J Ross
- Research Unit in Bioinformatics (RUBi), Department of Biochemistry and Microbiology, Rhodes University, Grahamstown, South Africa
| | - Ali Rana Atilgan
- Faculty of Engineering and Natural Sciences, Sabanci University, Istanbul, Turkey
| | - Özlem Tastan Bishop
- Research Unit in Bioinformatics (RUBi), Department of Biochemistry and Microbiology, Rhodes University, Grahamstown, South Africa
| | - Canan Atilgan
- Faculty of Engineering and Natural Sciences, Sabanci University, Istanbul, Turkey.
| |
Collapse
|
11
|
Abdizadeh H, Atilgan AR, Atilgan C, Dedeoglu B. Computational approaches for deciphering the equilibrium and kinetic properties of iron transport proteins. Metallomics 2018; 9:1513-1533. [PMID: 28967944 DOI: 10.1039/c7mt00216e] [Citation(s) in RCA: 16] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/27/2023]
Abstract
With the advances in three-dimensional structure determination techniques, high quality structures of the iron transport proteins transferrin and the bacterial ferric binding protein (FbpA) have been deposited in the past decade. These are proteins of relatively large size, and developments in hardware and software have only recently made it possible to study their dynamics using standard computational resources. We review computational techniques towards understanding the equilibrium and kinetic properties of iron transport proteins under different environmental conditions. At the level of detail that requires quantum chemical treatments, the octahedral geometry around iron has been scrutinized and it has been established that the iron coordinating tyrosines are in an unusual deprotonated state. At the atomistic level, both the N-lobe and the full bilobal structure of transferrin have been studied under varying conditions of pH, ionic strength and binding of other metal ions by molecular dynamics (MD) simulations. These studies have allowed questions to be answered, among others, on the function of second shell residues in iron release, the role of synergistic anions in preparing the active site for iron binding, and the differences between the kinetics of the N- and the C-lobe. MD simulations on FbpA have led to the detailed observation of the binding kinetics of phosphate to the apo form, and to the conformational preferences of the holo form under conditions mimicking the environmental niches provided by the periplasmic space. To study the dynamics of these proteins with their receptors, one must resort to coarse-grained methodologies, since these systems are prohibitively large for atomistic simulations. A study of the complex of human transferrin (hTf) with its pathogenic receptor by such methods has revealed a potential mechanistic explanation for the defense mechanism that arises in evolutionary warfare. Meanwhile, the motions in the transferrin receptor bound hTf have been shown to disfavor apo hTf dissociation, explaining why the two proteins remain in complex during the recycling process from the endosome to the cell surface. Open problems and possible technological applications related to metal ion binding-release in iron transport proteins that may be handled by hybrid use of quantum mechanical, MD and coarse-grained approaches are discussed.
Collapse
Affiliation(s)
- H Abdizadeh
- Faculty of Engineering and Natural Sciences, Sabancı University, Orhanlı 34956, Tuzla, Istanbul, Turkey.
| | | | | | | |
Collapse
|
12
|
Abdizadeh H, Tamer YT, Acar O, Toprak E, Atilgan AR, Atilgan C. Increased substrate affinity in the Escherichia coli L28R dihydrofolate reductase mutant causes trimethoprim resistance. Phys Chem Chem Phys 2018; 19:11416-11428. [PMID: 28422217 DOI: 10.1039/c7cp01458a] [Citation(s) in RCA: 21] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Dihydrofolate reductase (DHFR) is a ubiquitous enzyme with an essential role in cell metabolism. DHFR catalyzes the reduction of dihydrofolate to tetrahydrofolate, which is a precursor for purine and thymidylate synthesis. Several DHFR targeting antifolate drugs including trimethoprim, a competitive antibacterial inhibitor, have therefore been developed and are clinically used. Evolution of resistance against antifolates is a common public health problem rendering these drugs ineffective. To combat the resistance problem, it is important to understand resistance-conferring changes in the DHFR structure and accordingly develop alternative strategies. Here, we structurally and dynamically characterize Escherichia coli DHFR in its wild type (WT) and trimethoprim resistant L28R mutant forms in the presence of the substrate and its inhibitor trimethoprim. We use molecular dynamics simulations to determine the conformational space, loop dynamics and hydrogen bond distributions at the active site of DHFR for the WT and the L28R mutant. We also report their experimental kcat, Km, and Ki values, accompanied by isothermal titration calorimetry measurements of DHFR that distinguish enthalpic and entropic contributions to trimethoprim binding. Although mutations that confer resistance to competitive inhibitors typically make enzymes more promiscuous and decrease affinity to both the substrate and the inhibitor, strikingly, we find that the L28R mutant has a unique resistance mechanism. While the binding affinity differences between the WT and the mutant for the inhibitor and the substrate are small, the newly formed extra hydrogen bonds with the aminobenzoyl glutamate tail of DHF in the L28R mutant leads to increased barriers for the dissociation of the substrate and the product. Therefore, the L28R mutant indirectly gains resistance by enjoying prolonged binding times in the enzyme-substrate complex. While this also leads to slower product release and decreases the catalytic rate of the L28R mutant, the overall effect is the maintenance of a sufficient product formation rate. Finally, the experimental and computational analyses together reveal the changes that occur in the energetic landscape of DHFR upon the resistance-conferring L28R mutation. We show that the negative entropy associated with the binding of trimethoprim in WT DHFR is due to water organization at the binding interface. Our study lays the framework to study structural changes in other trimethoprim resistant DHFR mutants.
Collapse
Affiliation(s)
- Haleh Abdizadeh
- Faculty of Engineering and Natural Sciences, Sabanci University, Istanbul, Turkey.
| | | | | | | | | | | |
Collapse
|
13
|
Penkler DL, Atilgan C, Tastan Bishop Ö. Allosteric Modulation of Human Hsp90α Conformational Dynamics. J Chem Inf Model 2018; 58:383-404. [PMID: 29378140 DOI: 10.1021/acs.jcim.7b00630] [Citation(s) in RCA: 55] [Impact Index Per Article: 9.2] [Reference Citation Analysis] [Abstract] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
Central to Hsp90's biological function is its ability to interconvert between various conformational states. Drug targeting of Hsp90's regulatory mechanisms, including its modulation by cochaperone association, presents as an attractive therapeutic strategy for Hsp90 associated pathologies. In this study, we utilized homology modeling techniques to calculate full-length structures of human Hsp90α in closed and partially open conformations and used these structures as a basis for several molecular dynamics based analyses aimed at elucidating allosteric mechanisms and modulation sites in human Hsp90α. Atomistic simulations demonstrated that bound adenosine triphosphate (ATP) stabilizes the dimer by "tensing" each protomer, while adenosine diphosphate (ADP) and apo configurations "relax" the complex by increasing global flexibility, the former case resulting in a fully open "v-like" conformation. Dynamic residue network analysis revealed regions of the protein involved in intraprotein communication and identified several key communication hubs that correlate with known functional sites. Pairwise comparison of betweenness centrality, shortest path, and residue fluctuations revealed that a proportional relationship exists between the latter two measurables and an inverse relationship between these two and betweenness centrality. This analysis showed how protein flexibility, degree of compactness, and the distance cutoff used for network construction influence the correlations between these metrics. These findings are novel and suggest shortest path and betweenness centrality to be more relevant quantities to follow for detecting functional residues in proteins compared to residue fluctuations. Perturbation response scanning analysis identified several potential residue sites capable of modulating conformational change in favor of interstate conversion. For the ATP-bound open conformation, these sites were found to overlap with known Aha1 and client binding sites, demonstrating how naturally occurring forces associated with cofactor binding could allosterically modulate conformational dynamics.
Collapse
Affiliation(s)
- David L Penkler
- Research Unit in Bioinformatics (RUBi), Department of Biochemistry and Microbiology, Rhodes University , Grahamstown, 6140, South Africa
| | - Canan Atilgan
- Faculty of Engineering and Natural Sciences, Sabanci University , Tuzla 34956, Istanbul, Turkey
| | - Özlem Tastan Bishop
- Research Unit in Bioinformatics (RUBi), Department of Biochemistry and Microbiology, Rhodes University , Grahamstown, 6140, South Africa
| |
Collapse
|
14
|
Penkler D, Sensoy Ö, Atilgan C, Tastan Bishop Ö. Perturbation-Response Scanning Reveals Key Residues for Allosteric Control in Hsp70. J Chem Inf Model 2017; 57:1359-1374. [PMID: 28505454 DOI: 10.1021/acs.jcim.6b00775] [Citation(s) in RCA: 51] [Impact Index Per Article: 7.3] [Reference Citation Analysis] [Abstract] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/08/2023]
Abstract
Hsp70 molecular chaperones play an important role in maintaining cellular homeostasis, and are implicated in a wide array of cellular processes, including protein recovery from aggregates, cross-membrane protein translocation, and protein biogenesis. Hsp70 consists of two domains, a nucleotide binding domain (NBD) and a substrate binding domain (SBD), each of which communicates via an allosteric mechanism such that the protein interconverts between two functional states, an ATP-bound open conformation and an ADP-bound closed conformation. The exact mechanism for interstate conversion is not as yet fully understood. However, the ligand-bound states of the NBD and SBD as well as interactions with cochaperones such as DnaJ and nucleotide exchange factor are thought to play crucial regulatory roles. In this study, we apply the perturbation-response scanning (PRS) method in combination with molecular dynamics simulations as a computational tool for the identification of allosteric hot residues in the large multidomain Hsp70 protein. We find evidence in support of the hypothesis that substrate binding triggers ATP hydrolysis and that the ADP-substrate complex favors interstate conversion to the closed state. Furthermore, our data are in agreement with the proposal that there is an allosterically active intermediate state between the open and closed states and vice versa, as we find evidence that ATP binding to the closed structure and peptide binding to the open structure allosterically "activate" the respective complexes. We conclude our analysis by showing how our PRS data fit the current opinion on the Hsp70 conformational cycle and present several allosteric hot residues that may provide a platform for further studies to gain additional insight into Hsp70 allostery.
Collapse
Affiliation(s)
- David Penkler
- Research Unit in Bioinformatics (RUBi), Department of Biochemistry and Microbiology, Rhodes University , Grahamstown 6140, South Africa
| | - Özge Sensoy
- School of Engineering and Natural Sciences, Istanbul Medipol University , Beykoz 34810, Istanbul, Turkey
| | - Canan Atilgan
- Faculty of Engineering and Natural Sciences, Sabanci University , Tuzla 34956, Istanbul, Turkey
| | - Özlem Tastan Bishop
- Research Unit in Bioinformatics (RUBi), Department of Biochemistry and Microbiology, Rhodes University , Grahamstown 6140, South Africa
| |
Collapse
|
15
|
Abdizadeh H, Atilgan AR, Atilgan C. Mechanisms by Which Salt Concentration Moderates the Dynamics of Human Serum Transferrin. J Phys Chem B 2017; 121:4778-4789. [DOI: 10.1021/acs.jpcb.7b02380] [Citation(s) in RCA: 16] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/26/2023]
Affiliation(s)
- Haleh Abdizadeh
- Sabanci University, Faculty of Engineering and Natural Sciences, Tuzla, 34956 Istanbul, Turkey
| | - Ali Rana Atilgan
- Sabanci University, Faculty of Engineering and Natural Sciences, Tuzla, 34956 Istanbul, Turkey
| | - Canan Atilgan
- Sabanci University, Faculty of Engineering and Natural Sciences, Tuzla, 34956 Istanbul, Turkey
| |
Collapse
|
16
|
Sang P, Du X, Yang LQ, Meng ZH, Liu SQ. Molecular motions and free-energy landscape of serine proteinase K in relation to its cold-adaptation: a comparative molecular dynamics simulation study and the underlying mechanisms. RSC Adv 2017. [DOI: 10.1039/c6ra23230b] [Citation(s) in RCA: 24] [Impact Index Per Article: 3.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022] Open
Abstract
The physicochemical bases for enzyme cold-adaptation remain elusive.
Collapse
Affiliation(s)
- Peng Sang
- Laboratory of Molecular Cardiology
- Department of Cardiology
- The First Affiliated Hospital of Kunming Medical University
- Kunming
- P. R. China
| | - Xing Du
- Laboratory for Conservation and Utilization of Bio-Resources
- Yunnan University
- Kunming
- P. R. China
- Department of Biochemistry and Molecular Biology
| | - Li-Quan Yang
- College of Agriculture and Biological Science
- Dali University
- Dali
- P. R. China
| | - Zhao-Hui Meng
- Laboratory of Molecular Cardiology
- Department of Cardiology
- The First Affiliated Hospital of Kunming Medical University
- Kunming
- P. R. China
| | - Shu-Qun Liu
- Laboratory of Molecular Cardiology
- Department of Cardiology
- The First Affiliated Hospital of Kunming Medical University
- Kunming
- P. R. China
| |
Collapse
|
17
|
Abdizadeh H, Atilgan C. Predicting long term cooperativity and specific modulators of receptor interactions in human transferrin from dynamics within a single microstate. Phys Chem Chem Phys 2016; 18:7916-26. [DOI: 10.1039/c5cp05107j] [Citation(s) in RCA: 18] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
PRS identifies regions contacting rapidly evolving residues that mechanically manipulate dissociation from the pathogen in the human transferrin–bacterial receptor complex.
Collapse
Affiliation(s)
- Haleh Abdizadeh
- Faculty of Engineering and Natural Sciences
- Sabanci University
- Tuzla
- Turkey
| | - Canan Atilgan
- Faculty of Engineering and Natural Sciences
- Sabanci University
- Tuzla
- Turkey
| |
Collapse
|