Electrostatic switches that mediate the pH-dependent conformational change of "short" recombinant human pseudocathepsin D

Basis for Accurate Protein pKa Prediction with Machine Learning

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 pK_a's. To accelerate pH-dependent molecular mechanism research in the life sciences or industrial protein and drug designs, fast and accurate pK_a prediction is crucial. Here we present a theoretical pK_a 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 pK_a 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 pK_a prediction tools driven by artificial intelligence. In addition, DeepKa built on PHMD549 has been proven an efficient protein pK_a predictor and thus can be applied immediately to, for example, pK_a database construction, protein design, drug discovery, and so on.

Predicting Catalytic Proton Donors and Nucleophiles in Enzymes: How Adding Dynamics Helps Elucidate the Structure-Function Relationships

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 K_a's of the buried and hydrogen bonded catalytic dyads in five enzymes. We asked the question what determines the p K_a order, i.e., what makes a residue proton donor vs a nucleophile. The continuous constant pH molecular dynamics simulations captured the experimental p K_a 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.

Hybrid 2D/3D-quantitative structure-activity relationship and modeling studies perspectives of pepstatin A analogs as cathepsin D inhibitors

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].

Conformational dynamics of cathepsin D and binding to a small-molecule BACE1 inhibitor

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.

Molecular Dynamics Simulations of Ligand-Induced Flap Conformational Changes in Cathepsin-D-A Comparative Study

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.

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

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.

New cis-configured aziridine-2-carboxylates as aspartic acid protease inhibitors

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.

Neutralizing positive charges at the surface of a protein lowers its rate of amide hydrogen exchange without altering its structure or increasing its thermostability

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).

Recombinant plasmepsin 1 from the human malaria parasite plasmodium falciparum: enzymatic characterization, active site inhibitor design, and structural analysis

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.