1
|
Cai Z, Liu T, Lin Q, He J, Lei X, Luo F, Huang Y. Basis for Accurate Protein p Ka Prediction with Machine Learning. J Chem Inf Model 2023; 63:2936-2947. [PMID: 37146199 DOI: 10.1021/acs.jcim.3c00254] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/07/2023]
Abstract
pH regulates protein structures and the associated functions in many biological processes via protonation and deprotonation of ionizable side chains where the titration equilibria are determined by pKa's. To accelerate pH-dependent molecular mechanism research in the life sciences or industrial protein and drug designs, fast and accurate pKa prediction is crucial. Here we present a theoretical pKa data set PHMD549, which was successfully applied to four distinct machine learning methods, including DeepKa, which was proposed in our previous work. To reach a valid comparison, EXP67S was selected as the test set. Encouragingly, DeepKa was improved significantly and outperforms other state-of-the-art methods, except for the constant-pH molecular dynamics, which was utilized to create PHMD549. More importantly, DeepKa reproduced experimental pKa orders of acidic dyads in five enzyme catalytic sites. Apart from structural proteins, DeepKa was found applicable to intrinsically disordered peptides. Further, in combination with solvent exposures, it is revealed that DeepKa offers the most accurate prediction under the challenging circumstance that hydrogen bonding or salt bridge interaction is partly compensated by desolvation for a buried side chain. Finally, our benchmark data qualify PHMD549 and EXP67S as the basis for future developments of protein pKa prediction tools driven by artificial intelligence. In addition, DeepKa built on PHMD549 has been proven an efficient protein pKa predictor and thus can be applied immediately to, for example, pKa database construction, protein design, drug discovery, and so on.
Collapse
Affiliation(s)
- Zhitao Cai
- College of Computer Engineering, Jimei University, Xiamen 361021, China
| | - Tengzi Liu
- College of Computer Engineering, Jimei University, Xiamen 361021, China
| | - Qiaoling Lin
- College of Computer Engineering, Jimei University, Xiamen 361021, China
| | - Jiahao He
- College of Computer Engineering, Jimei University, Xiamen 361021, China
| | - Xiaowei Lei
- College of Computer Engineering, Jimei University, Xiamen 361021, China
| | - Fangfang Luo
- College of Computer Engineering, Jimei University, Xiamen 361021, China
| | - Yandong Huang
- College of Computer Engineering, Jimei University, Xiamen 361021, China
| |
Collapse
|
2
|
Huang Y, Yue Z, Tsai CC, Henderson JA, Shen J. Predicting Catalytic Proton Donors and Nucleophiles in Enzymes: How Adding Dynamics Helps Elucidate the Structure-Function Relationships. J Phys Chem Lett 2018; 9:1179-1184. [PMID: 29461836 PMCID: PMC6555141 DOI: 10.1021/acs.jpclett.8b00238] [Citation(s) in RCA: 27] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/23/2023]
Abstract
Despite the relevance of understanding structure-function relationships, robust prediction of proton donors and nucleophiles in enzyme active sites remains challenging. Here we tested three types of state-of-the-art computational methods to calculate the p Ka's of the buried and hydrogen bonded catalytic dyads in five enzymes. We asked the question what determines the p Ka order, i.e., what makes a residue proton donor vs a nucleophile. The continuous constant pH molecular dynamics simulations captured the experimental p Ka orders and revealed that the negative nucleophile is stabilized by increased hydrogen bonding and solvent exposure as compared to the proton donor. Surprisingly, this simple trend is not apparent from crystal structures and the static structure-based calculations. While the generality of the findings awaits further testing via a larger set of data, they underscore the role of dynamics in bridging enzyme structures and functions.
Collapse
Affiliation(s)
- Yandong Huang
- Department of Pharmaceutical Sciences , University of Maryland School of Pharmacy , Baltimore , Maryland 21201 , United States
| | - Zhi Yue
- Department of Pharmaceutical Sciences , University of Maryland School of Pharmacy , Baltimore , Maryland 21201 , United States
| | - Cheng-Chieh Tsai
- Department of Pharmaceutical Sciences , University of Maryland School of Pharmacy , Baltimore , Maryland 21201 , United States
| | - Jack A Henderson
- Department of Pharmaceutical Sciences , University of Maryland School of Pharmacy , Baltimore , Maryland 21201 , United States
| | - Jana Shen
- Department of Pharmaceutical Sciences , University of Maryland School of Pharmacy , Baltimore , Maryland 21201 , United States
| |
Collapse
|
3
|
Hybrid 2D/3D-quantitative structure-activity relationship and modeling studies perspectives of pepstatin A analogs as cathepsin D inhibitors. Future Med Chem 2017; 10:5-26. [PMID: 29235371 DOI: 10.4155/fmc-2017-0077] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/08/2023] Open
Abstract
AIM Cathepsin D, one of the attractive targets in the treatment of breast cancer, has been implicated in HIV neuropathogenesis with potential proteolytic effects on chemokines. Methodology/result: Diverse modeling tools were used to reveal the key structural features affecting the inhibitory activities of 78 pepstatin A analogs. Analyses were performed to investigate the stability, rationality and fluctuation of the analogs. Results showed a clear correlation between the experimental and predicted activities of the analogs as well as the variation in their activities relative to structural modifications. CONCLUSION The insight gained from this study offers theoretical references for understanding the mechanism of action of cathepsin D and will aid in the design of more potent and clinically-relevant drugs. Graphical abstract [Formula: see text].
Collapse
|
4
|
Ellis CR, Tsai CC, Lin FY, Shen J. Conformational dynamics of cathepsin D and binding to a small-molecule BACE1 inhibitor. J Comput Chem 2017; 38:1260-1269. [PMID: 28370344 PMCID: PMC5403572 DOI: 10.1002/jcc.24719] [Citation(s) in RCA: 19] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/10/2016] [Revised: 11/21/2016] [Accepted: 12/06/2016] [Indexed: 02/06/2023]
Abstract
BACE1 is a major therapeutic target for prevention and treatment of Alzheimer's disease. Developing inhibitors that can selectively target BACE1 in favor of other proteases, especially cathepsin D (CatD), has presented significant challenges. Here, we investigate the conformational dynamics and protonation states of BACE1 and CatD using continuous constant pH molecular dynamics with pH replica-exchange sampling protocol. Despite similar structure, BACE1 and CatD exhibit markedly different active site dynamics. BACE1 displays pH-dependent flap dynamics that controls substrate accessibility, while the CatD flap is relatively rigid and remains open in the pH range 2.5-6. Interestingly, although each protease hydrolyzes peptide bonds, the protonation states of the catalytic dyads are different within the active pH range. The acidic and basic components of the BACE1 catalytic dyad are clear, while either aspartic acid of the CatD catalytic dyad could play the role of acid or base. Finally, we investigate binding of the inhibitor LY2811376 developed by Eli Lilly to BACE1 and CatD. Surprisingly, in the enzyme active pH range, LY2811376 forms a stronger salt bridge with the catalytic dyad in CatD than in BACE1, which might explain the retinal toxicity of the inhibitor related to off-target inhibition of CatD. This work highlights the complexity and challenge in structure-based drug design where receptor-ligand binding induces protonation state change in both the protein and the inhibitor. © 2017 Wiley Periodicals, Inc.
Collapse
Affiliation(s)
- Christopher R. Ellis
- Department of Pharmaceutical Sciences, University of Maryland School of Pharmacy, Baltimore, MD 21201
| | - Cheng-Chieh Tsai
- Department of Pharmaceutical Sciences, University of Maryland School of Pharmacy, Baltimore, MD 21201
| | - Fang-Yu Lin
- Department of Pharmaceutical Sciences, University of Maryland School of Pharmacy, Baltimore, MD 21201
| | - Jana Shen
- Department of Pharmaceutical Sciences, University of Maryland School of Pharmacy, Baltimore, MD 21201
| |
Collapse
|
5
|
Arodola OA, Soliman MES. Molecular Dynamics Simulations of Ligand-Induced Flap Conformational Changes in Cathepsin-D-A Comparative Study. J Cell Biochem 2016; 117:2643-57. [PMID: 27038253 DOI: 10.1002/jcb.25564] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/09/2016] [Accepted: 03/31/2016] [Indexed: 12/11/2022]
Abstract
The flap region in aspartic proteases is a unique structural feature to this class of enzymes, and found to have a profound impact on protein overall structure, function, and dynamics. Understanding the structure and dynamic behavior of the flap regions is crucial in the design of selective inhibitors against aspartic proteases. Cathepsin-D, an aspartic protease enzyme, has been implicated in a long list of degenerative diseases as well as breast cancer progression. Presented herein, for the first time, is a comprehensive description of the conformational flap dynamics of cathepsin-D using a comparative 50 ns "multiple" molecular dynamics simulations. Diverse collective metrics were proposed to accurately define flap dynamics. These are distance d1 between the flap tips residues (Gly79 and Met301); dihedral angle ϕ; in addition to TriCα angles Gly79-Asp33-Asp223, θ1 , and Gly79-Asp223-Met301, θ2 . The maximum distance attained throughout the simulation was 17.42 and 11.47 Å for apo and bound cathepsin-D, respectively, while the minimum distance observed was 8.75 and 6.32 Å for apo and bound cathepsin-D, respectively. The movement of the flap as well as the twist of the active pocket can properly be explained by measuring the angle, θ1 , between Gly79-Asp33-Met301 and correlating it with the distance Cα of the flap tip residues. The asymmetrical opening of the binding cavity was best described by the large shift of -6.26° to +20.94° in the dihedral angle, ϕ, corresponding to the full opening of the flap at a range of 31-33 ns. A wide-range of post-dynamic analyses was also applied in this report to supplement our findings. We believe that this report would augment current efforts in designing potent structure-based inhibitors against cathepsin-D in the treatment of breast cancer and other degenerative diseases. J. Cell. Biochem. 117: 2643-2657, 2016. © 2016 Wiley Periodicals, Inc.
Collapse
Affiliation(s)
- Olayide A Arodola
- Molecular Modelling and Drug Design Research Group, School of Health Sciences, University of KwaZulu-Natal, Westville, Durban 4000, South Africa
| | - Mahmoud E S Soliman
- Molecular Modelling and Drug Design Research Group, School of Health Sciences, University of KwaZulu-Natal, Westville, Durban 4000, South Africa.
| |
Collapse
|
6
|
Morrow G, Kim HJ, Pellerito O, Bourrelle-Langlois M, Le Pécheur M, Groebe K, Tanguay RM. Changes in Drosophila mitochondrial proteins following chaperone-mediated lifespan extension confirm a role of Hsp22 in mitochondrial UPR and reveal a mitochondrial localization for cathepsin D. Mech Ageing Dev 2016; 155:36-47. [PMID: 26930296 DOI: 10.1016/j.mad.2016.02.011] [Citation(s) in RCA: 14] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/11/2015] [Revised: 01/20/2016] [Accepted: 02/25/2016] [Indexed: 10/22/2022]
Abstract
Hsp22 is a small mitochondrial heat shock protein (sHSP) preferentially up-regulated during aging in Drosophila melanogaster. Its developmental expression is strictly regulated and it is rapidly induced in conditions of stress. Hsp22 is one of the few sHSP to be localized inside mitochondria, and is the first sHSP to be involved in the mitochondrial unfolding protein response (UPR(MT)) together with Hsp60, mitochondrial Hsp70 and TRAP1. The UPR(MT) is a pro-longevity mechanism, and interestingly Hsp22 over-expression by-itself increases lifespan and resistance to stress. To unveil the effect of Hsp22 on the mitochondrial proteome, comparative IEF/SDS polyacrylamide 2D gels were done on mitochondria from Hsp22+ flies and controls. Among the proteins influenced by Hsp22 expression were proteins from the electron transport chain (ETC), the TCA cycle and mitochondrial Hsp70. Hsp22 co-migrates with ETC components and its over-expression is associated with an increase in mitochondrial protease activity. Interestingly, the only protease that showed significant changes upon Hsp22 over-expression in the comparative IEF/SDS-PAGE analysis was cathepsin D, which is localized in mitochondria in addition to lysosome in D. melanogaster as evidenced by cellular fractionation. Together the results are consistent with a role of Hsp22 in the UPR(MT) and in mitochondrial proteostasis.
Collapse
Affiliation(s)
- Geneviève Morrow
- Laboratoire de Génétique Cellulaire et Développementale, Département de biologie moléculaire, biochimie médicale et pathologie, Institut de Biologie Intégrative et des Systèmes (IBIS) and PROTEO, Université Laval, Québec, G1V 0A6, Canada
| | - Hyun-Ju Kim
- Laboratoire de Génétique Cellulaire et Développementale, Département de biologie moléculaire, biochimie médicale et pathologie, Institut de Biologie Intégrative et des Systèmes (IBIS) and PROTEO, Université Laval, Québec, G1V 0A6, Canada
| | - Ornella Pellerito
- Laboratoire de Génétique Cellulaire et Développementale, Département de biologie moléculaire, biochimie médicale et pathologie, Institut de Biologie Intégrative et des Systèmes (IBIS) and PROTEO, Université Laval, Québec, G1V 0A6, Canada
| | - Maxime Bourrelle-Langlois
- Laboratoire de Génétique Cellulaire et Développementale, Département de biologie moléculaire, biochimie médicale et pathologie, Institut de Biologie Intégrative et des Systèmes (IBIS) and PROTEO, Université Laval, Québec, G1V 0A6, Canada
| | - Marie Le Pécheur
- Laboratoire de Génétique Cellulaire et Développementale, Département de biologie moléculaire, biochimie médicale et pathologie, Institut de Biologie Intégrative et des Systèmes (IBIS) and PROTEO, Université Laval, Québec, G1V 0A6, Canada
| | | | - Robert M Tanguay
- Laboratoire de Génétique Cellulaire et Développementale, Département de biologie moléculaire, biochimie médicale et pathologie, Institut de Biologie Intégrative et des Systèmes (IBIS) and PROTEO, Université Laval, Québec, G1V 0A6, Canada.
| |
Collapse
|
7
|
Büchold C, Hemberger Y, Heindl C, Welker A, Degel B, Pfeuffer T, Staib P, Schneider S, Rosenthal PJ, Gut J, Morschhäuser J, Bringmann G, Schirmeister T. New cis-configured aziridine-2-carboxylates as aspartic acid protease inhibitors. ChemMedChem 2011; 6:141-52. [PMID: 21082722 DOI: 10.1002/cmdc.201000370] [Citation(s) in RCA: 14] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/11/2022]
Abstract
A series of 52 cis-configured 1-alkyl-3-phenylaziridine-2-carboxylates were synthesized as new pseudo-irreversible inhibitors of Candida albicans secreted aspartic acid protease 1 (SAP1), SAP2, SAP3, and SAP8. Some of the compounds, which were obtained as diastereomers with S,S- and R,R-configured aziridine rings by Cromwell synthesis of racemic (2R,3S+2S,3R)-dibromophenylpropionic acid ester with amines, followed by ester hydrolysis and coupling to hydrophobic amino acid esters, were separated by preparative HPLC. The absolute configuration of the aziridine ring was assigned by a combination of experimental circular dichroism (CD) investigations and quantum chemical CD calculations. In agreement with previous docking studies, the diastereomers all exhibit similar activity. The compounds were found to be more active against the related mammalian enzyme cathepsin D, presumably due to productive interactions of the N-alkyl substituent with the highly lipophilic S2 pocket. The most active inhibitors (5, 9, 10, 21, and 28), characterized by benzyl, cyclohexylmethyl, tert-butyl, or 1,4-dimethylpentyl moieties at the aziridine nitrogen atom, exhibit k(2nd) values between 500 and 900×10³ M⁻¹ min⁻¹ and K(i) values near or below 1 μM for cathepsin D.
Collapse
Affiliation(s)
- Christian Büchold
- Institute of Pharmacy and Food Chemistry, University of Würzburg, Germany
| | | | | | | | | | | | | | | | | | | | | | | | | |
Collapse
|
8
|
Shaw BF, Arthanari H, Narovlyansky M, Durazo A, Frueh DP, Pollastri MP, Lee A, Bilgicer B, Gygi SP, Wagner G, Whitesides GM. Neutralizing positive charges at the surface of a protein lowers its rate of amide hydrogen exchange without altering its structure or increasing its thermostability. J Am Chem Soc 2010; 132:17411-25. [PMID: 21090618 DOI: 10.1021/ja9067035] [Citation(s) in RCA: 29] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022]
Abstract
This paper combines two techniques--mass spectrometry and protein charge ladders--to examine the relationship between the surface charge and hydrophobicity of a representative globular protein (bovine carbonic anhydrase II; BCA II) and its rate of amide hydrogen-deuterium (H/D) exchange. Mass spectrometric analysis indicated that the sequential acetylation of surface lysine-ε-NH3(+) groups--a type of modification that increases the net negative charge and hydrophobicity of the surface of BCA II without affecting its secondary or tertiary structure--resulted in a linear decrease in the aggregate rate of amide H/D exchange at pD 7.4, 15 °C. According to analysis with MS, the acetylation of each additional lysine generated between 1.4 and 0.9 additional hydrogens that are protected from H/D exchange during the 2 h exchange experiment at 15 °C, pD 7.4. NMR spectroscopy demonstrated that none of the hydrogen atoms which became protected upon acetylation were located on the side chain of the acetylated lysine residues (i.e., lys-ε-NHCOCH3) but were instead located on amide NHCO moieties in the backbone. The decrease in rate of exchange associated with acetylation paralleled a decrease in thermostability: the most slowly exchanging rungs of the charge ladder were the least thermostable (as measured by differential scanning calorimetry). This observation--that faster rates of exchange are associated with slower rates of denaturation--is contrary to the usual assumptions in protein chemistry. The fact that the rates of H/D exchange were similar for perbutyrated BCA II (e.g., [lys-ε-NHCO(CH2)2CH3]18) and peracetylated BCA II (e.g., [lys-ε-NHCOCH3]18) suggests that the electrostatic charge is more important than the hydrophobicity of surface groups in determining the rate of H/D exchange. These electrostatic effects on the kinetics of H/D exchange could complicate (or aid) the interpretation of experiments in which H/D exchange methods are used to probe the structural effects of non-isoelectric perturbations to proteins (i.e., phosphorylation, acetylation, or the binding of the protein to an oligonucleotide or to another charged ligand or protein).
Collapse
Affiliation(s)
- Bryan F Shaw
- Department of Chemistry and Chemical Biology, Harvard University, Cambridge, Massachusetts 02138, USA.
| | | | | | | | | | | | | | | | | | | | | |
Collapse
|
9
|
Masson O, Bach AS, Derocq D, Prébois C, Laurent-Matha V, Pattingre S, Liaudet-Coopman E. Pathophysiological functions of cathepsin D: Targeting its catalytic activity versus its protein binding activity? Biochimie 2010; 92:1635-43. [DOI: 10.1016/j.biochi.2010.05.009] [Citation(s) in RCA: 74] [Impact Index Per Article: 5.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/05/2010] [Accepted: 05/14/2010] [Indexed: 11/27/2022]
|
10
|
Liu P, Marzahn MR, Robbins AH, Gutiérrez-de-Terán H, Rodríguez D, McClung SH, Stevens SM, Yowell CA, Dame JB, McKenna R, Dunn BM. Recombinant plasmepsin 1 from the human malaria parasite plasmodium falciparum: enzymatic characterization, active site inhibitor design, and structural analysis. Biochemistry 2009; 48:4086-99. [PMID: 19271776 DOI: 10.1021/bi802059r] [Citation(s) in RCA: 18] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
A mutated form of truncated proplasmepsin 1 (proPfPM1) from the human malaria parasite Plasmodium falciparum, proPfPM1 K110pN, was generated and overexpressed in Escherichia coli. The automaturation process was carried out at pH 4.0 and 4.5, and the optimal catalytic pH of the resulting mature PfPM1 was determined to be pH 5.5. This mature PfPM1 showed comparable binding affinity to peptide substrates and inhibitors with the naturally occurring form isolated from parasites. The S3-S3' subsite preferences of the recombinant mature PfPM1 were explored using combinatorial chemistry based peptide libraries. On the basis of the results, a peptidomimetic inhibitor (compound 1) was designed and yielded 5-fold selectivity for binding to PfPM1 versus the homologous human cathepsin D (hcatD). The 2.8 A structure of the PfPM2-compound 1 complex is reported. Modeling studies were conducted using a series of peptidomimetic inhibitors (compounds 1-6, Table 3) and three plasmepsins: the crystal structure of PfPM2, and homology derived models of PfPM1 and PfPM4.
Collapse
Affiliation(s)
- Peng Liu
- Department of Biochemistry and Molecular Biology, University of Florida, College of Medicine, Gainesville, Florida 32610-0245, USA
| | | | | | | | | | | | | | | | | | | | | |
Collapse
|
11
|
Degel B, Staib P, Rohrer S, Scheiber J, Martina E, Büchold C, Baumann K, Morschhäuser J, Schirmeister T. Cis-Configured Aziridines Are New Pseudo-Irreversible Dual-Mode Inhibitors ofCandida albicans Secreted Aspartic Protease 2. ChemMedChem 2008; 3:302-15. [DOI: 10.1002/cmdc.200700101] [Citation(s) in RCA: 17] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/10/2022]
|