Next Article in Journal
Real-World Analysis of Nivolumab and Atezolizumab Efficacy in Previously Treated Patients with Advanced Non-Small Cell Lung Cancer
Next Article in Special Issue
Substrate-Based Design of Cytosolic Nucleotidase IIIB Inhibitors and Structural Insights into Inhibition Mechanism
Previous Article in Journal
Hypericum perforatum and Its Ingredients Hypericin and Pseudohypericin Demonstrate an Antiviral Activity against SARS-CoV-2
Previous Article in Special Issue
Correction: Alfei, S.; Zuccari, G. Recommendations to Synthetize Old and New β-Lactamases Inhibitors: A Review to Encourage Further Production. Pharmaceuticals 2022, 15, 384
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Covalent Inhibition of the Human 20S Proteasome with Homobelactosin C Inquired by QM/MM Studies

by
Natalia Serrano-Aparicio
,
Silvia Ferrer
and
Katarzyna Świderek
*
BioComp Group, Department of Physical and Analytical Chemistry, Universitat Jaume I, 12071 Castellón de la Plana, Spain
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Pharmaceuticals 2022, 15(5), 531; https://doi.org/10.3390/ph15050531
Submission received: 31 March 2022 / Revised: 21 April 2022 / Accepted: 22 April 2022 / Published: 25 April 2022
(This article belongs to the Special Issue Design of Enzyme Inhibitors as Potential Drugs 2022)

Abstract

:
20S proteasome is a main player in the protein degradation pathway in the cytosol, thus intervening in multiple pivotal cellular processes. Over the years the proteasome has emerged as a crucial target for the treatment of many diseases such as neurodegenerative diseases, cancer, autoimmune diseases, developmental disorders, cystic fibrosis, diabetes, cardiac diseases, atherosclerosis, and aging. In this work, the mechanism of proteasome covalent inhibition with bisbenzyl-protected homobelactosin C (hBelC) was explored using quantum mechanics/molecular mechanics (QM/MM) methods. Molecular dynamic simulations were used to describe key interactions established between the hBelC and its unique binding mode in the primed site of the β5 subunit. The free energy surfaces were computed to characterize the kinetics and thermodynamics of the inhibition process. This study revealed that although the final inhibition product for hBelC is formed according to the same molecular mechanism as one described for hSalA, the free energy profile of the reaction pathway differs significantly from the one previously reported for γ-lactam-β-lactone containing inhibitors in terms of the height of the activation barrier as well as the stabilization of the final product. Moreover, it was proved that high stabilization of the covalent adduct formed between β5-subunit and hBelC, together with the presence of aminocarbonyl side chain in the structure of the inhibitor which prevents the hydrolysis of the ester bond from taking place, determines its irreversible character.

Graphical Abstract

1. Introduction

The ubiquitin-proteasome system (UPS) is responsible for hydrolyzing more than 90% of all cytosolic proteins, where the 26S proteasome, a large ~2.5 MDa molecular complex, is the central player in this non-lysosomal protein degradation pathway [1]. The 26S proteasome is formed by two subcomplexes, the catalytic core particle (CP) or 20S proteasome, and one or more 19S regulatory particles (RP) that serve as activators. This complex system was discovered in 1969 when an ATP-dependent proteolysis process was identified [2]. Later it was connected to a stable polypeptide, named ubiquitin for being present in all tissues of eukaryotic organisms [3,4]. Eventually, the intracellular proteolytic system was identified as a big multimolecular protein initially called “cylindrin” due to its cylindric shape [5], afterward named as the 20S proteasome [6,7,8] and confirmed to be the catalytic core of the ATP dependent proteolysis in the UPS [9]. An overall achievement in the characterization of the ubiquitin conjugation system and its role in specific proteolysis labeling was awarded with the Nobel Prize in Chemistry in 2004 to Avram Hershko, Aaron Ciechanover, and Irwin Rose [10].
It was discovered that the UPS degrades proteins in a multistep process, where target proteins are linked by covalent bonds to the ubiquitin in presence of ATP [11,12]. This small protein binds to the 19S RP of the 26S proteasome, where the protein destined for degradation is unfolded and inserted into the multicatalytic chamber, the CP, to be rapidly discomposed into smaller polypeptides [13,14,15]. Since the 20S proteasome’s primary role is to degrade key regulatory and aberrant proteins [16], it procures the maintenance of cellular homeostasis controlling protein turnover. [15] Being the main protein degradation pathway in the cytosol, UPS has a wide spectrum of polypeptide substrates, and its function intervenes in multiple crucial cellular processes [17]. The UPS participates in cell cycle progression, [18,19,20] cellular signaling [21,22], genome integrity [23,24], regulation of endocrine pathways [25,26], apoptosis [27,28], transcriptional regulation [29,30], metabolism regulation [31,32], immune responses [33,34], and in various human diseases and pathogenesis [35,36]. Therefore, all these functions make the proteasome a crucial target for therapeutic intervention against, for instance, neurodegenerative diseases [37,38], cancer [39,40,41], autoimmune diseases, developmental disorders, cystic fibrosis, diabetes, cardiac diseases, atherosclerosis, or just aging [42]. Some infectious diseases such as Chaga’s disease [43], malaria [44,45,46], and tuberculosis [47] have been also related to proteasomal function.
Except for bacteria (excluding actinomycetes), which harbor complexes of two homohexameric rings with one active subunit each [48], the 20S proteasomes consist of four stacked homoheptameric rings formed by two types of basic subunits known as α and β in the rest of living organisms. As shown in Figure 1 each of the two outer rings is formed by seven α subunits and each of the two inner rings by seven β subunits.
The α proteins have a key role in the complex assembly due to the formation of the entry gate to the interior of the proteasome. The β subunits bear the catalytic activity. [49,50] In eukaryotes, only three active β subunits (i.e., β1, β2, and β5 in constitutive variant) were found, that present different specificity towards the substrates due to the different architecture of S1 pockets located in their binding sites. The β5 subunit, characterized by the presence of Met45 in the hydrophobic S1 pocket, is known as chymotrypsin-like protease. The β2 subunit with hydrophilic S1 cavity with negatively charged Asp55 residue is known as trypsin-like protein. Finally, the β1 owes its caspase-like character to positively charged Arg45 residue located in the bottom of its hydrophilic pocket [51].
Most known proteasome inhibitors are peptide-like molecules bearing electrophilic warheads that reversibly or irreversibly bind to the catalytic OγThr1 forming covalent adducts. Therefore, proteasome inhibitors are commonly classified according to the type of functional groups of this reactive portion of the molecule. Seven different inhibitor classes are recognized at the moment, i.e., aldehydes, boronic acids, α,β-epoxyketones, α-ketoaldehydes (glyoxals), vinyl sulfones, vinyl amides (syrbactins), and β-lactones [52]. So far, only two boronic-acid, Bortezomib (Velcade®, Millennium Pharmaceuticals, Cambridge, MA, USA), and Ixazomib (Ninlaro®, Takeda Pharmaceuticals, MA, USA), and one α,β-epoxyketone, Carfilzomib (Kyprolis®, Amgen Inc., Thousand Oaks, CA, USA), have been approved by the U.S. FDA [53,54,55,56,57]. All three molecules inhibit the β5/β5i subunit.
Herein, the inhibition mechanism of the β5 subunit from human 20S proteasome with the homo-analog of belactosin C (BelC) was explored using a computational approach including molecular dynamic (MD) simulations with classical and quantum mechanics/molecular mechanics (QM/MM) potentials. Originally, isolated from Streptomyces sp. UCK14, the BelC was found to exhibit antitumor activity [58,59]. This activity was shown to increase significantly upon acetylation of the free amino group and esterification or amidation of the carboxyl group, as well as displacement of the ornithine moiety with lysine, resulting in the homo-analog homobelactosin C. Both structures of BelC and its homolog are shown in Figure 2A. Introduced modifications to the structure of BelC resulted in IC50 on the nanomolar level measured against human pancreoma and colon cancer [60,61]. This high antitumor activity was related to proteasome inhibition [62,63,64,65]. The crystal structure obtained for yeast CP [66] in a complex with bis-benzyl-protected homobelactosin C (hBelC) covalently bound in the active site suggested that the β-lactone (equally as other covalent proteasome inhibitors) undergoes nucleophilic attack by oxygen gamma of N-termini threonine (OγThr1) resulting in ester bond formation, as shown in Figure 2B. Nevertheless, it was demonstrated in the same study that this inhibitor presents high selectivity towards the primed site of the chymotrypsin-like β5 active site, adopting a completely different orientation to the one observed previously for β-lactone-containing proteasome inhibitors such as omuralide (Omu) or salinosporamide A (SalA). Moreover, based on the final product orientation, it was suggested that hBelC follows a unique mechanism preventing cleavage of its covalent adduct formed with the active site. The resistance of the ester bond was assigned to the presence of a 4-aminocarbonyl side chain that according to crystallographic studies can block the access of the water molecule due to its specific orientation in the binding pocket.
Therefore, the main goal of this work was to clarify the mechanism of inhibition involving hBelC, including recognition and inactivation steps. Insight into the specific interactions established between the inhibitor and residues from the primed site of the enzyme can be very helpful for several reasons. Mainly, it can provide a guideline that can be used to refine the recognition step by eliminating unfavorable interactions and improving the existing favorable ones. This is especially important in new strategies used to design polypharmacological anticancer drugs [67] where hBelC could serve as a template for the dual inhibition of the proteasome and fatty acid as proposed by Romo and co-workers. [61] Moreover, the information gathered at the molecular level on the ligand interaction with this part of the active site may contribute to a better understanding of the catalytic activity of the β5 subunit, especially since little is known about interactions of the substrate with this region.
In the case of the inactivation step, the mechanism of inhibition was studied, and the structure of the final product complex was analyzed to shed light on the hydrolysis step of the newly formed covalent adduct. Free energy surfaces were explored assuming the most favorable molecular mechanism, which was proposed based on our previous experience. Finally, a comparative analysis of results obtained for the hBelC and γ-lactam-β-lactone, homo-salinosporamide A (hSalA) was carried out to provide information on similarities and differences in the inhibition process offered by molecules belonging to the same warhead family.

2. Results and Discussion

In general, covalent inhibition is a two-step process. [68,69] First, an inhibitor reversibly associates with the target enzyme, and a protein-inhibitor complex (E-I) is formed. The potency of this first step is defined by the binding constant Ki. In the second step, the reaction occurs between the two reactive entities of the inhibitor and the enzyme, respectively, forming a protein-inhibitor covalent complex (E-I). This chemical transformation is characterized by the rate constant of inactivation, kinact as summarized in Equation (1)
E + I K i E I k i n a c t E I
where E stands for the enzyme, I is an inhibitor, EI is a non-covalent enzyme-inhibitor complex, and E-I represents a covalent adduct formed between enzyme and inhibitor. Computational research on covalent inhibitors often focuses solely on the energetics of the kinetic step, as high accuracy insights into the binding process seem to be still beyond the scope of currently available in silico techniques. Nevertheless, in this work, long MD simulations with a classical force field were carried out for the covalent E-hBelC complex to shed light on the recognition step. Subsequently, the mechanism of covalent adduct formation that occurs in the inactivation process was explored using QM/MM methods. The obtained results for both steps are described below.

2.1. Recognition Step

A key part of the optimization of inhibitors in pharmaceutical drug development is to vary the molecular design to enhance the complementarity of the chemical features of the compound with the position of the side chains of the amino acids in the binding pocket of a target enzyme [70]. In the special case of the covalent inhibitors, the recognition step requires an additional condition to be fulfilled. This is related to the optimal relative position of the warhead with respect to the catalytically important residues to facilitate the covalent adduct formation. According to crystallographic data, hBelC seems to fulfill both requirements satisfactorily. Nevertheless, because MD simulations were carried out in the final product of inhibition, the position of the warhead will not be analyzed here in detail, and therefore all attention will be focused on the behavior of the inhibitor tail within the binding cavity.
The shape of the site that binds hBelC to the active site of the β5 chain is, in fact, determined by the orientation of five subunits contained in the two β rings of the proteasome, i.e., β6, β5, and β4 subunits from one ring, and β3′ and β4′ from the neighboring one. As shown in Figure 3, the biggest portion of the hBelC inhibitor is bound to the β5 subunit that involves an active site together with the S1 pocket. Analysis of the hydrophobicity map confirms the existence of well-known high hydrophobicity of this pocket, at the same time indicating the hydrophilic character of the active site due to the presence of Lys33 with the hydrophobicity of −1.5 according to Eisenberg hydrophobicity scale [71] and Thr1 with the hydrophobicity of −0.05. Due to the molecular characterization of the inhibitor structure, the S1′ pocket remains empty and consecutive S2′ and S3′ pockets of the primed site are occupied by O-Bn (benzyl group) and methyl group, respectively.
It was observed that the N-Cbz group (benzyl carbamate) attached at the end of the hBelC tail does not enter the S4′ pocket, even though this pocket is highly hydrophobic. The hydrophobic character of this pocket is achieved due to the presence of Ile25 and Val26 of the β4, Val32 of the β5, and Phe170 and Ile171 of the β4′ subunit with the values of the hydrophobicity ranging between +1.08 and +1.38 (as provided by the Eisenberg scale). It is possible that the behavior of this N-Cbz group originates in its large size or is dictated by the attachment of the benzyl group to the sp2 oxygen instead sp3 carbon atom, as would be expected for the side chain of polypeptide substrate. In consequence, the presence of an oxygen atom in this position may force the aromatic ring to adopt an alternative orientation from desired in the cavity. The lack of geometrical complementarity of the N-Cbz group towards this pocket results in its high mobility as shown in Figure 4A,B where the various conformations explored by the N-Cbz group within the active site are shown together with the evolution of root-mean-square-deviation (RMSD) computed for the heavy atoms of the inhibitor during the simulations. As shown on the plots, the inhibitor experiences the highest deviations from the original position due to the flexibility of its tail. This finding proves that the N-Cbz group does not create specific interactions within the binding site of the 20S proteasome and is not crucial for the recognition step.
On the other hand, the O-Bn group seems to play an important role, because it binds very tightly inside the corresponding S2′ pocket and remains there during the overall MD simulations. The size and shape of this pocket seem to be established by the presence of the side chains of three residues located on its bottom, Asp115, Val114, and Tyr113, by the beta-sheet, that consists of beta-strands of the β5 subunit on one side, and by the backbone of Ser23, and the side chain of Asn24 from the β4 subunit on the other. Favorable interactions, as shown in Figure 4, created between these residues and hBelC ensure the stability of its binding.
Finally, the hydrogen bond (H-bond) contacts established between hBelC and the binding cavity during MD simulations were determined (see Figure 4). Based on the obtained results, we concluded that it is rather unlikely that the recognition step depends on this type of interaction. In addition to Gly47, which is involved in H-bond interaction with the carbonyl oxygen of the β-lactone ring, only one more H-bond was identified, established between the NH2-group of the side chain of Asn24 from the β4 chain and the oxygen of the peptide bond from the backbone of the side-chain of lactone ring of the hBelC. In all, the primed site occupied by this inhibitor was not found to provide more H-bond interactions that could be indicated as key for the recognition step.
It is very likely that the minimum requirement for an inhibitor structure to form a non-covalent adduct with the active site of the β5 subunit of the human 20S proteasome is primarily related to the molecular architecture of the P1, as has been proven for Omu or SalA. However, it seems that in the case of hBelC where a small P1 group is employed, the presence of the O-Bn group in the P2′ position may play an important supporting role.

2.2. Inactivation Step

In our previous studies, the mechanisms of inhibition with promising candidates, such as dihydroeponemycin [72] and SalA [73], i.e., α,β-epoxyketone, and γ-lactam-β-lactone, have been revealed using QM/MM MD methods. Based on computationally delivered results supported by the newest experimental data [74], it was concluded that the originally proposed and widely accepted mechanism of inhibition suggesting that only one residue of the active site, N-termini Thr1 [75], is involved in the process was incomplete. As demonstrated by obtained results, active participation of another residue, i.e., Lys33 was required to accomplish inhibition within a realistic energetic cost. [72] The role of this residue in the inactivation step of inhibition was further confirmed in studies focused on homo-salinosporamide A (hSalA) [76].
Therefore, guided by the previously acquired knowledge, we assume that in the case of the hBelC, the inhibition mechanism should be very close to the one previously described for hSalA. Thus, it was proposed, that the inhibition process takes place in three steps, as illustrated in Figure 5.
Accordingly, it was proposed that in the first step, oxygen gamma (Oγ) of Thr1 is activated by nitrogen zeta (Nζ) of Lys33, which abstracts the hydrogen gamma (Hγ), and consequently, the Oγ atom attacks the carbonyl carbon (C1βlac) of the β-lactone ring generating a covalent adduct in the form of a tetrahedral intermediate (E-TI). In the second step, the bond between carbonyl carbon (C1βlac) and the oxygen (O2βlac) of β-lactone is cleft resulting in a ring-opening. This step is accompanied by the accumulation of negative charge on O2βlac. As demonstrated in our previous work, [76] such a newly formed adduct is expected to be thermodynamically unfavorable, and thus an additional chemical step was proposed to neutralize the separation of charges observed in the active site of intermediate 2 (E-I2), i.e., the negative charge of the oxygen (O2βlac) of the open lactone ring and positive charge on nitrogen zeta (Nζ) of Lys33. A desired final neutral product may be obtained, as proposed previously, [76] by rearranging the positions of hydrogen atoms within the active site. This rearrangement could be achieved by a double proton transfer that would involve the transfer of one proton from Lys33 to the nitrogen (NThr1) of the amino group of Thr1 and the second proton from the same amino group to the oxygen (O2βlac) of the open β-lactone. Formation of the stable product of inhibition (E-PC) is expected once this last process is accomplished.
The existence of E-TI was previously excluded based on experimental studies on β-lactone inhibitors in serine proteases [77] in which it was suggested that the tetrahedral transition state may not be generated in the acylation reaction pathway, due to the high ring strain energy within the β-lactone. Consequently, we decided to explore the evolution of the first and second steps on the same free energy surface to provide evidence on whether the nucleophilic attack and ring-opening take place in a stepwise manner. Therefore, both chemical transformations were explored by controlling the antisymmetric combination of the d(Oγ-Hγ) and d(Hγ-Nζ) distances, which describe the process of proton transfer from Thr1 to Lys33, together with the antisymmetric combination of d(O2βlac-C1βlac) and d(Oγ-C1βlac) distances that define the process of β-lactone ring-opening and nucleophilic attack of Thr1 on its carbonyl carbon, respectively. The computed free energy surface is provided in Figure 6.
The minimum free energy pathway from non-covalent adduct (E-hBelC) leads, as originally assumed, to a metastable tetrahedral intermediate (E-TI) through the first transition state (TS1). Subsequently, the system is crossing the second transition state (TS2) to form the second intermediate (E-I2). The existence of both, TS1 and TS2 was confirmed by optimization and characterization of their structures at the M06-2X/AMBER level of theory as shown in Figure 7. All key distances and charges obtained for optimized stationary structures and determined during reaction progress are provided in Tables S1 and S2 of Supplementary Materials.
Subsequently, the free energy surface for the last step was explored controlling the d(Hγ-NThr1) distance describing the transfer of the proton (Hγ) originally located on Lys33 to the nitrogen atom (NThr1) of NH2-group of Thr1, together with the d(HThr1-O2βlac) distance that describes second proton transfer from the amino group of threonine to the negatively charged oxygen of open β-lactone ring (see Figure 6). The minimum energy path indicates that both protons are transferred in a concerted way, and the NH2- group of Thr1 serves as a proton shuttle in the inactivation step. A detail of the structure of the transition state (TS3) for the last step of the inhibition mechanism, optimized at the M06-2X/AMBER level of theory, is shown in Figure 7.
Both computed free energy surfaces were consequently used to obtain the free energy profile for the complete inactivation step of the inhibition which is presented in Figure 8. Since the mechanism of inhibition with hBelC and hSalA appear to be identical considering the sequence of events occurring during the chemical reaction, the newly obtained energetic profile was further confronted with the same profile computed for the inhibition of β5-subunit of the 20S proteasome with hSalA, deduced from our previous work [76].
According to the computed Gibbs free energy profile, it was found that the rate-determining step corresponds, as in other cases of β-lactone containing inhibitors, to the first chemical step, i.e., the nucleophilic attack of Thr1 on the β-lactone ring. However, and surprisingly, the obtained free energy barrier of 24.9 kcal·mol−1 for hBelC is meaningfully higher than the one computed for molecules with a γ-lactam-β-lactone warhead that were ranging between 20.3 and 20.4 kcal·mol−1 for SalA and hSala, respectively [73,76] This different result for hBelC may be caused by the higher flexibility of the single β-lactone ring due to (i) the presence of the sec-butyl group in the P1 position, and (ii) the absence of γ-lactam ring in the warhead structure that induces rigidness of the ring and ensures the position of the electrophilic center towards the oxygen gamma of Thr1. The sec-butyl can be responsible for weaker recognition characteristics than the hydroxymethyl-cyclohexenyl group of SalA or hSalA. Alternatively, the high barrier could simply originate in the weaker electrophilic center offered by the smaller warhead of hBelC. To explore these possible scenarios, first, the ESP partial charges were computed on each atom belonging to the warheads of hBelC and hSalA inhibitors in gas phase. As shown in Figure 9A, the carbonyl carbon (C1βlac) of the β-lactone ring carries a higher positive charge in the case of hSalA (0.717 e) than in hBelC (0.637 e) due to its intramolecular characteristic. Therefore, a lower positive charge computed on C1βlac in hBelC confirms that this warhead should provide a weaker electrophilic center. Although a relatively small shift of 0.08 e between charge distribution on C1βlac atom in the β-lactone with and without attachment of a γ-lactam ring was observed, nevertheless it can result in meaningful change in the free energy barrier. In fact, such an effect was already observed previously in the case of C–N bond hydrolysis catalyzed by Candida antarctica lipase B where the charge shift of only 0.12 e on the nucleophile acceptor determined the possibility of the reaction taking place [78].
Possible deviation of the β-lactone ring of hBelC from the optimal position was examined by analyzing the electrophile-nucleophile distance together with the Bürgi–Dunitz [79,80,81] (αBD) and Flippin–Lodge [82] (αFL) angles that describe the nucleophilic attack trajectory. The obtained values of αBD and αFL were compared with those determined for the first step of inhibition with hSalA. The evolution of these two angles along the progress of the first step collected for hBelC and hSalA inhibitors is shown in Figure 9B As can be seen, the relative position between the nucleophile (OγThr1) and electrophilic center (C1βlac) slightly changes when comparing both Michaelis complexes. It was found that while in E•hBelC complex the αBD of 102° is closer to the value of 105 ± 5° (value determined by Bürgi et al. as the angle that ensures a reliable position for the nucleophile to attack) than in E•hSalA complex (αBD = 97°), and the estimated value of αFL, that in the ideal case is expected to be 0°, presents higher deviation of the nucleophile from the normal plane to the electrophile plane in E•hBelC complex (αFL = 9°) than in E•hSalA (αFL = 6°). Thus, a higher deviation of αFL angle from the ideal angle equal to 0° together with the lower divergence of αBD versus the ideal angle of 105°, measured at the initiation point of the first step of inhibition, do not allow for the conclusion that the increment of the free energy barrier computed for the first step of inhibition with hBelC is related to the energetic cost of the nucleophile and electrophile position rearrangement. Moreover, it is rather expected that due to the opposite behavior of αBD and αFL angles their overall contribution to energetic change would rather cancel out and the energetic barrier should be closer to the one computed originally for hSalA.
On the contrary, the electrophile-nucleophile distance was found to be longer in E•hBelC complex (2.46 Å) than in E•hSalA (2.35 Å). Therefore, elongation of the electrophile-nucleophile distance together with a weaker electrophilic center offered by hBelC can account as energetically unfavorable factors that can explain the rise in the barrier height observed when the hBelC is used to block the 20S proteasome activity.
The barrier of the rate-limiting step is not the only difference found when comparing the energetic profiles computed for the inhibition process with hBelC and hSalA. As shown in Figure 8, meaningful dissimilarity in the stabilization of intermediate 2 (E-I2) was detected. This complex, as described before, is a covalent adduct of Thr1 with an open β-lactone ring whose creation is accompanied by NζLys(+) and O2βlac(−) ionic pair formation. Interestingly, while in hSalA this complex is c.a. 10 kcal·mol−1 higher than the non-covalent complex (EhSalA), in hBelC it is almost equally stable as EhBelC, with computationally predicted relative free energy of 1.9 kcal·mol−1. Therefore, this unexpected result was further analyzed using both structural and energetic factors.
First, the structure of E-I2 intermediates formed between β5-subunit and hBelC and hSalA were compared by the geometrical overlay of the active sites. The result of this structural alignment is shown in Figure 10A. It was observed that although atoms of the P1 group, three carbons of the β-lactone ring, and the oxygen atom of the carbonyl group (O1βlac) from both inhibitors occupy almost identical positions in the active site, the orientation adopted by the oxygen (O2βlac) of the open ring in the case of hBelC, is meaningfully deviated from the one observed for hSalA. As indicated in Figure 10A, this O2βlac atom of hBelC, besides the translation, experiences additional rotation during C1βlac-O2βlac bond cleavage that resulted in changing its position towards the one originally occupied by a methyl group (–CH3) in the lactam ring in hSalA. This atom rearrangement was confirmed by a comparative analysis of the evolution of dihedral angles defined by the positions of C1βlac-C4βlac-C3βlac-O2βlac atoms in hBelC, and C1βlac-C4βlac-C3βlac-O2βlac and C1βlac-C4βlac-C3βlac-CMeβlac atoms in hSalA, respectively, which were collected along the second step of inhibition, that corresponds to C1βlac-O2βlac bond cleavage. The evolution of dihedrals shown in Figure 10B indicates that in order to accomplish this chemical step, the C1βlac-O2βlac bond scission in hBelC is accompanied by the β-lactone ring puckering inversion. The dramatic change of puckering coordinate from c.a. 11° in E-TI to c.a. −3° in TS2 was observed for the four-member ring of hBelC. Puckering inversion of the β-lactone ring was not observed in the case of hSalA for which values of dihedral defined by C1βlac-C4βlac-C3βlac-O2βlac atoms changes from c.a. 19° in E-TI to 27° in TS2.
Since it is expected that such an evident conformational change as the one observed for hBelC cannot go unnoticed by the environment, the next step of the analysis was to test if the better stabilization of E-I2 complex in the case of this inhibitor can originate in the improvement of its interactions established within the active site, by comparison with hSalA. For this purpose, interaction energy (Eint), as a sum of electrostatic (Eelec) and vdW (EvdW) interactions, was computed for hSalA and hBelC in their E-I2 adducts. The results of these calculations are presented in Figure 10C,D. For clarity, only those residues for which the obtained difference in computed Eint for studied inhibitors was higher than 1 kcal·mol−1 were provided. During the analysis, it was found that despite two residues, i.e., Asp17, and Asp167 that increase meaningfully their unfavorable interactions with hBelC with respect to hSalA (from 9.7 ± 3.3 to 18.8 ± 3.0, and from 14.2 ± 3.4 to 18.8 ± 3.2 kcal·mol−1, respectively), the conformational change in hBelC resulted in gaining much more favorable contacts. Thus, a meaningful rise in the favorable interactions was observed (higher than 5 kcal·mol−1 by comparison to Eint values computed in E-I2 for hSalA) for hBelC interacting with Arg19, Lys33, Val128, Gly129, and Ser130 residues that can be considered as a key for stabilizing newly formed E-I2 complex. In all cases, the dominant contribution of Eelec to the overall value of Eint is observed. Importantly, it was found that the interaction with Gly47 from the oxyanion hole was also affected by the conformational change as reflected by the higher value of Eint computed for hBelC in comparison to hSalA, i.e., −7.8 ± 1.3 versus −2.6 ± 3.3 kcal·mol−1, respectively.
Regardless the E-I2 complex formation was found to be thermodynamically favored for hBelC by comparison with the same intermediate obtained for hSalA inhibitor, the formation of its structure did not ensure overall product stability. Thus, as commented before, the neutralization of the active site by double proton transfer was required to obtain the desired irreversible product of inhibition. As delivered from calculations, to reach the final product complex (E-PC) the system must cross over an additional not negligible free energy barrier of 19.1 kcal·mol−1. However, the newly created E-PC was found to be much more stable in Gibbs free energy, 15.6 kcal·mol−1 lower than the non-covalent E•hBelC complex and 7.1 kcal·mol−1 lower in comparison to the final product of inhibition formed with hSalA (E-hSalA). This higher stability of the final product in the case of hBelC can be explained using the same arguments as discovered in the analysis carried out for the E-I2 complex. Therefore, and contrary to hSalA which was proven to reversibly inhibit the β5-subunit of the 20S proteasome, the inhibition reaction with hBelC seems to present an irreversible character, because the reversible process that would lead to non-covalent adduct from E-PC requires to surmount the barrier of 40.5 kcal·mol−1, what it is rather unfeasible.
The remaining alternative pathway to recover the activity of the enzyme, once the E-PC is formed is through the hydrolysis of the ester bond formed between the Thr1 and C1βlac atom of the lactam ring. The process of hydrolysis would follow a standard reaction mechanism pathway in which in the first step water molecule would attack the C1βlac atom. This step can be achieved within a reasonable energetic cost only if the initial position of the water molecule is geometrically favorable. In other words, the nucleophilic attack trajectory for this molecule should not be disturbed by the presence of another residue from the enzyme or substituents of the inhibitor. In fact, we have previously provided quantitative proof that in the case of hSalA water molecules cannot access the C1βlac atom due to the presence of the hydroxyl group in opened β-lactone ring that resulted in very high free energy barriers for this process. Moreover, it was also demonstrated that the attack on the planar electrophilic center can be achieved only from one side that is opposite to the localization of the highly hydrophobic pocket occupied by the P1 group. [76] Very useful in this analysis appears to be the evaluation of the behavior of the distance established between O2βlac and C1βlac atoms as well as αBD and αFL angles in the E-PC complex. Therefore, in the case of hBelC, results of 100 ns MD simulations obtained for the E-PC model were analyzed in this respect. The distribution of chosen distances together with αBD and αFL angles explored during MD simulations is provided in Figure 11. In this analysis, not only the relative position of O2βlac to C1βlac was examined but also of NhBelC and C1′hBelC atoms of an aminocarbonyl side chain. The atom positions in the structure of hBelC are given in Figure 12.
As commented before, the optimal position for the nucleophilic attack would be defined by a reasonable short distance to the electrophilic center complemented with a suitable location of the nucleophile in the three-dimensional space, which can be judged based on αBD and αFL values. If any of chosen herein atoms for analysis fulfill these requirements it can be assumed that it creates a steric hindrance and prevents water access to the target of its attack. Since this role was proven to be played by O2βlac in previous work, our analysis of hBelC inhibitor started from this atom. Despite, the distance between O2βlac and C1βlac being short (3.3 ± 0.1 Å) in the E-PC complex, as expected, the αBD is far from the ideal value (105°) and it turned out to be much more obtuse (140 ± 1°) than the same angle of 104° determined for inhibition product created with hSalA along the full relaxation of the system along MD simulations. [76] This change must be dictated by the rotation experienced by O2βlac along the second step of the reaction, which, as discussed before, induces higher stabilization of the E-I2 but at the same time, it exposes access to the C1βlac carbon. Nevertheless, it was observed that due to the characteristic of the hBelC structure, a nitrogen atom, NhBelC of an aminocarbonyl side chain can take over the role of O2βlac and prevent hydrolysis. Considering the preserved short NhBelC-C1βlac distance of 3.3 ± 0.1 Å and the very promising αBD and αFL values of 91 ± 1° and 1 ± 1°, it can be concluded that the hydrolysis of E-PC is rather unfeasible. Therefore, the results of this analysis agree with the previous proposal suggesting that hBelC follows a unique mechanism preventing cleavage of its covalent adduct formed with the active site [66].
Finally, the short distance of 2.7 Å established between O2βlac and oxygen of Arg19 (OArg19) of the protein backbone that was observed in the X-ray structure of the E-hBelC [66], was also analyzed and was found to be long (3.68 ± 0.25 Å) in the structure obtained after computational exploration of the last process of the inactivation step. However, as it was observed during longer MD simulations the contact between the hydroxyl group of open lactone ring and Arg17 is oscillating between 2.7 and 7.1 Å (as shown in Figure S1 of Supplementary Materials) indicating that the H-bond interaction established between these atoms in the case of the human 20S proteasome and in the given condition of temperature is very labile.

3. Computational Methods

3.1. System Setup

The molecular model for the inhibition studies of human 20S proteasome with hBelC was prepared using the crystal structure of the CP (PDB ID: 5LF1) [83] with bound SalA molecule. This proteasome-SalA complex was prepared and used in our previous studies. [73] For the purpose of this work, SalA was substituted by the structure of hBelC in its final inhibition product as found in the crystal structure of 20S proteasome from Saccharomyces cerevisiae (PDB ID: 3E47) [66] Missing force field parameters for the hBelC covalent adduct was obtained using the Antechamber software [84] (see Table S3 of Supplementary Materials). The pKa shift of all titratable residues was predicted with PropKa software ver. 3.1 [85,86], providing the same results and consequently the same protonation states at pH 7.0 as described in our previous work. [72] Missing hydrogen atoms were added to the enzyme together with 43 positively charged sodium counterions that were put in the most electrostatically favorable positions in order to neutralize the total negative charge of the system. Hydrogen and counterions were added using the tLEAP [87] module of the AmberTools package. Subsequently, the system was soaked within an orthorhombic box of TIP3P [88] water molecules, with a size of 17.1 × 16.7 × 20.0 nm3 providing a model consisting of 542,592 atoms. To describe the protein and water molecules the AMBER [89] and TIP3P force fields, respectively, were employed and the NAMD [90] software was used as an MD engine. The equilibration protocol for MD simulations involved a preliminary minimization and gradual heating of the system to 310 K with 0.001 K temperature increments, followed by 100 ps of non-biased NPT equilibration and 100 ps of non-biased NVT equilibration. The final model was used as starting structure for long NVT MD simulations. A cut-off for non-bonding interactions was set between 14.5 to 16 Å using a smooth switching function. Although improved algorithms to compute non-electrostatic interactions have been recently developed [91,92] in this work interactions between the QM and MM moieties were treated by means of standard Lennard-Jones potential [93]. During 100 ns of MD simulations, in order to significantly improve time and reduce the cost of calculations, the positions of all residues beyond 40 Å from the inhibitor were fixed. 40 Å is slightly higher than the double size of radius chosen as the external limit for the cut-off function. The temperature during the simulations was controlled using the Langevin thermostat, [94] and the pressure for the NPT equilibration with the Nosé-Hoover Langevin piston [95] pressure control.
All simulations were analyzed according to the evolution of the RMSD of the protein backbone as well as of heavy atoms of hBelC and key distances. The system was considered equilibrated after 50 ns of MD simulations (see Figure 4B) and Figure S2 of Supplementary Materials).

3.2. QM/MM Calculations

A representative structure of the product complex E-hBelC from the MD simulation was selected as the starting point for exploring the inhibition mechanism. The size of the system was reduced, cutting the box of waters to a water sphere of 60 Å centered on residue Thr1 of the β5 active site, leading to a system with 149,520 atoms. The AMBER and TIP3P force fields, as implemented in the fDynamo library [96,97], were used for describing the protein, counterions, and water solvent. The cut-off scheme for non-bonding interactions was the same as the one used in the classical MD simulations protocol. Additionally, the positions of the atoms of residues beyond 20 Å from the inhibitor were fixed. The QM part of the model was selected including only the essential 72 atoms, as shown in Figure 12. This includes part of the inhibitor, where the warhead is contained, full N-termini Thr1 residue, part of residue Thr2, and the backbone of residue Lys33. Three-link atoms [98] were inserted on the QM/MM boundaries, placed in the Cα-Cβ bond of Lys33, the C-Cα bond of Thr2, and C2′-C3′ bond in the aliphatic chain of the inhibitor. The Austin Model, AM1 [99] semiempirical Hamiltonian, and the Minnesota density functional, M06-2X [100] with the standard 6-31+G(d,p) basis set, were employed to treat the QM sub-set of atoms.

3.3. Potential Energy Surfaces (PES)

Adequate distinguished reaction coordinates, ξ, were selected for exploring each chemical step of the inhibition mechanism with hBelC. Each PESs was generated by grid scanning, where the step size for hydrogen transfer was controlled every 0.05 Å, while distances between heavy atoms were changed by 0.1 Å. The resulting PES (shown in Figure S3 of Supplementary Materials) allowed identifying the minimum energy path (MEP). All stationary points observed along the MEP were optimized at a low and high level of theory (AM1/AMBER and M06-2X/AMBER, respectively) employing a micro-macro iteration method using Baker’s algorithm [101] and they were characterized by computing the matrix of the second energy derivatives. All computed TSs at M06-2X/MM (the Cartesian coordinates for QM atoms are provided in Table S4 of Supplementary Materials) were connected to the expected minima by tracing intrinsic reaction coordinate (IRC) paths, to further confirm the explored inhibition mechanism.

3.4. Free Energy Surfaces (FES)

A series of QM/MM MD simulations employing the umbrella sampling (US) method [102] as implemented in fDynamo, were computed using previously generated structures at PES at the AM1/MM level of theory. A force constant of 2500 kJ·mol−1·Å−2 was employed to constrain the reaction coordinate and 310 K was set as the simulation temperature. An initial equilibration of 5 ps was carried out, followed by 20 ps of production at every window. Finally, the weighted histogram analysis method (WHAM) [103] was used to integrate the obtained results in terms of potentials of mean force (PMF). A density tolerance of 10−3 to consider the WHAM calculation converged.
The PMF was obtained as a function of the distinguished reaction coordinate, W(ξ). The PMF is related to the normalized probability of finding the system at a particular value of the chosen coordinate, as shown in Equation (2),
W ξ = C k T l n ρ r N δ ξ r N ξ d r N 1
Then, the free energy of activation can be expressed as follows:
Δ G ξ = W ξ W ξ R + G ξ ξ R
where the superscripts indicate the value of the reaction coordinate at the reactants (R), and at the TS (‡), and GξR) is the free energy associated with setting the reaction coordinate to a specific value at the reactant state. Normally this last term makes a small contribution, and the activation free energy is directly estimated from the PMF change between the maximum of the profile and the reactant’s minimum,
Δ G ξ W ξ W ξ R = Δ W ξ
The selection of the reaction coordinate is usually trivial when the mechanism can be driven by a single internal coordinate or a simple combination (as the antisymmetric combination of two interatomic distances). However, this is not the case for the reaction steps explored in this work. Instead, it was necessary to obtain a much more computationally demanding 2D-PMF using two coordinates: ξ1 and ξ2. The 2D-PMF is related to the probability of finding the system at distinct values of these two coordinates,
W ξ = C k T l n ρ r N δ ξ 1 r N ξ 1 δ ξ 2 r N ξ 2 d r N 2
To estimate the free energy barriers from this quantity, the one-dimensional PMF changes were recovered and a maximum probability reaction path on the 2D-PMF surface was traced integrating over the perpendicular coordinate.

3.5. Spline Corrections (SP)

To improve the quality of obtained results and to reduce the possible errors associated with the semiempirical method used during free energy simulations, high-level corrections were employed using DFT method. Spline corrections were applied through the energy function [104,105,106]:
E = E L L / M M + S Δ E L L H L ξ 1 , ξ 2
where the final energy is obtained from a correction term computed using the single-point energy difference between the high-level (HL) and the low-level (LL) for the QM sub-set of atoms. As mentioned above, the AM1 method was applied as LL method while as the HL the hybrid M06-2X functional with the standard 6-31+G(d,p) basis set was used. The Gaussian09 [107] program combined with fDynamo was employed for DFT/AMBER calculations. In this work, a total of 1798 and 1568 US windows and the respective single-point calculations at the M06-2X/MM level were computed to produce the final QM/MM free energy profiles.

4. Conclusions

In this work, the key molecular factors that are responsible for the irreversible character of hBelC in the covalent inhibition of β5-subunit of human 20S proteasome have been revealed using classical MD and QM/MM MD simulations. The herein studied β-lactone containing inhibitor, which is a close analog to a natural product, BelC, seems to be a good candidate because it reveals promising activity against human pancreoma and colon cancer.
The present study provides insights into the recognition and inactivation steps that must be completed to successfully block the activity of the enzyme according to the covalent inhibition mechanism. It shows that hBelC consists of two key fragments that ensure the recognition step, i.e., the sec-butyl and O-benzyl group in P1 and P2′ positions, respectively. Both substituents bind in hydrophobic cavities that ensure favorable weak nonbonding interactions. Moreover, results suggest that while the presence of sec-butyl is crucial as it is responsible for the positioning of the warhead in the active site, the role of the substituent present at P2′ position is rather secondary.
The study of the inactivation step reveals that hBelC irreversibly inhibits the β5-subunit according to the same molecular mechanism as the one observed on the inhibition with hSalA. [76] The most probable reaction mechanism based on the computed free energy profiles indicates that hBelC blocks the active site of the enzyme according to the same sequence of chemical events as previously described for hSalA. Although all three steps of reaction were found to be identical, the computed free energy profile showed meaningful differences. The activation barrier for the rate-limiting step that corresponds to the nucleophilic attack of Thr1 on C1βlac carbon is higher than the one reported previously for hSalA. Computational results revealed that the observed rise in the barrier is due to the weaker electrophilic center on C1βlac offered by hBelC, as well as the longer distance established between nucleophile and electrophile in Michaelis complex, E•hBelC than in E•hSalA. Additionally, stronger stabilization of the final product complex ensures the irreversible character of hBelC which is, in fact, opposite to the reversible behavior of hSalA. Based on our analysis it was concluded that the higher stability of the final covalent complex E-hBelC is related to the conformational change of the β-lactone ring where the flexible oxygen experiences additional rotation during reaction progress, which causes the appearance of more favorable interactions with the residues of the active site.
The existence of an alternative hydrolysis pathway allowing the recovery of the β5-subunit activity was also excluded. The results of our studies have provided quantitative evidence that hydrolysis of the final product of inhibition with hBelC is unfeasible due to the presence of the aminocarbonyl side chain which most probably blocks the water molecule access to the electrophilic center by occupying the most optimal position for nucleophilic attack in the active site.
We strongly believe that the knowledge gained in this study for the inhibition mechanism of the human 20S proteasome with hBelC explored at the atomistic level may be useful in the future. This can be especially relevant in redesigning structures of the already known candidates or designing new ones with additional properties such as for instance polypharmacological anticancer drugs, where hBelC could serve as a template for the dual inhibition of the proteasome and fatty acid. Finally, the results obtained in this work can additionally serve as a reference in a future study on the 20S immunoproteasome providing evidence confirming or contradicting the hypothesis that small sequence changes observed in the primed site of β5i subunit are responsible for increased hBelC affinity towards its binding pocket.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/ph15050531/s1. Table S1: Key distances for stationary structures optimized and characterized at M06-2X/AMBER level of theory with 6-31+(d,p) basis set, along 20S proteasome inhibition pathway with hBelC. Table S2: ESP charges on key atoms computed for stationary structures optimized and characterized at M06-2X/AMBER level of theory with 6-31+(d,p) basis set. Figure S1: Evolution of the distance established between O2βlac and OArg19 during 100 ns MD simulation in E-PC complex. Table S3: Force field parameters for the hBelC covalent adduct formed with Thr1 of β5 subunit of 20S proteasome. Figure S2: Time-dependent evolution of the total energy, temperature, RMSD of the protein backbone, RMSD of the atoms of hBelC and evolution of key distances during 100 ns of classical MD simulations in E-PC complex. Figure S3: Free energy surfaces computed at AM1/AMBER level for a full three-step inhibition process.; Table S4: Cartesian coordinates of QM atoms of optimized and characterized transition state structures at M06-2X/AMBER level of theory with 6-31+(d,p) basis set, along 20S proteasome inhibition pathway with hBelC.

Author Contributions

N.S.-A. and S.F. performed the calculations. K.Ś. conceived and designed the study, supervised calculations, prepared the figures, and wrote the manuscript. K.Ś., N.S.-A. and S.F. discussed and analyzed the data. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Spanish Ministerio de Ciencia, Innovación y Universidades (ref. PID2019-107098RJ-I00 and ref: RYC2020-030596-I), the Generalitat Valenciana (ref. SEJI/2020/007), and Universitat Jaume I (UJI-A2019-04).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data is contained within the article or Supplementary Materials.

Acknowledgments

This work was supported by the Spanish Ministerio de Ciencia, Innovación y Universidades (ref. PID2019-107098RJ-I00), the Generalitat Valenciana (ref. SEJI/2020/007), Universitat Jaume I (UJI-A2019-04), K.Ś. thanks Ministerio de Ciencia, Innovación y Universidades for a Ramón y Cajal contract (ref: RYC2020-030596-I). The authors acknowledge computational resources from the Servei d’Informàtica of Universitat Jaume I.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Hyams, J.S. Gastroenterology and Nutrition, 2nd ed.; Elsevier Saunders: Philadelphia, PA, USA, 2001; Volume 13, ISBN 978-1-4377-2603-9. [Google Scholar]
  2. Etlinger, J.D.; Goldberg, A.L. A soluble ATP-dependent proteolytic system responsible for the degradation of abnormal proteins in reticulocytes. Proc. Natl. Acad. Sci. USA 1977, 74, 54–58. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Ciehanover, A.; Hod, Y.; Hershko, A. A heat-stable polypeptide component of an ATP-dependent proteolytic system from reticulocytes. Biochem. Biophys. Res. Commun. 1978, 81, 1100–1105. [Google Scholar] [CrossRef]
  4. Goldstein, G.; Scheid, M.; Hammerling, U.; Schlesinger, D.H.; Niall, H.D.; Boyse, E.A. Isolation of a polypeptide that has lymphocyte-differentiating properties and is probably represented universally in living cells. Proc. Natl. Acad. Sci. USA 1975, 72, 11–15. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  5. Harris, J. The isolation and purification of a macromolecular protein component from the human erythrocyte ghost. Biochim. Biophys. Acta (BBA)-Protein Struct. 1969, 188, 31–42. [Google Scholar] [CrossRef]
  6. Hase, J.; Kobashi, K.; Nakai, N.; Mitsui, K.; Iwata, K.; Takadera, T. The quaternary structure of carp muscle alkaline protease. Biochim. Biophys. Acta (BBA)-Enzym. 1980, 611, 205–213. [Google Scholar] [CrossRef]
  7. Arrigo, A.-P.; Tanaka, K.; Goldberg, A.L.; Welch, W.J. Identity of the 19S ‘prosome’ particle with the large multifunctional protease complex of mammalian cells (the proteasome). Nature 1988, 331, 192–194. [Google Scholar] [CrossRef]
  8. Wilk, S.; Orlowski, M. Cation-Sensitive Neutral Endopeptidase: Isolation and Specificity of the Bovine Pituitary Enzyme. J. Neurochem. 1980, 35, 1172–1182. [Google Scholar] [CrossRef]
  9. Eytan, E.; Ganoth, D.; Armon, T.; Hershko, A. ATP-dependent incorporation of 20S protease into the 26S complex that degrades proteins conjugated to ubiquitin. Proc. Natl. Acad. Sci. USA 1989, 86, 7751–7755. [Google Scholar] [CrossRef] [Green Version]
  10. The Nobel Prize in Chemistry 2005 NobelPrize.Org. Nobel Prize Outreach AB 2021. Available online: https://www.nobelprize.org/prizes/chemistry/2005/summary/ (accessed on 4 September 2021).
  11. Hershko, A.; Ciechanover, A.; Heller, H.; Haas, A.L.; Rose, I.A. Proposed role of ATP in protein breakdown: Conjugation of protein with multiple chains of the polypeptide of ATP-dependent proteolysis. Proc. Natl. Acad. Sci. USA 1980, 77, 1783–1786. [Google Scholar] [CrossRef] [Green Version]
  12. Ciechanover, A.; Elias, S.; Heller, H.; Ferber, S.; Hershko, A. Characterization of the heat-stable polypeptide of the ATP-dependent proteolytic system from reticulocytes. J. Biol. Chem. 1980, 255, 7525–7528. [Google Scholar] [CrossRef]
  13. De Martino, G.N.; Proske, R.J.; Moomaw, C.R.; Strong, A.A.; Song, X.; Hisamatsu, H.; Tanaka, K.; Slaughter, C.A. Identification, Purification, and Characterization of a PA700-dependent Activator of the Proteasome. J. Biol. Chem. 1996, 271, 3112–3118. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Adams, G.M.; Crotchett, B.; Slaughter, C.A.; DeMartino, G.N.; Gogol, E.P. Formation of Proteasome−PA700 Complexes Directly Correlates with Activation of Peptidase Activity. Biochemistry 1998, 37, 12927–12932. [Google Scholar] [CrossRef]
  15. Thibaudeau, T.A.; Smith, D.M. A Practical Review of Proteasome Pharmacology. Pharmacol. Rev. 2019, 71, 170–197. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  16. Patrick, G.L. An Introduction to Medicinal Chemistry, 6th ed.; Oxford University Press: New York, NY, USA, 1995; ISBN 9780198749691. [Google Scholar]
  17. Mishra, R.; Upadhyay, A.; Prajapati, V.K.; Mishra, A. Proteasome-mediated proteostasis: Novel medicinal and pharmacological strategies for diseases. Med. Res. Rev. 2018, 38, 1916–1973. [Google Scholar] [CrossRef] [PubMed]
  18. Hershko, A.; Ciechanover, A. The ubiquitin system. Annu. Rev. Biochem. 1998, 67, 425–479. [Google Scholar] [CrossRef] [PubMed]
  19. Peters, J.-M. The Anaphase-Promoting Complex: Proteolysis in Mitosis and Beyond. Mol. Cell 2002, 9, 931–943. [Google Scholar] [CrossRef]
  20. Bassermann, F.; Eichner, R.; Pagano, M. The ubiquitin proteasome system—Implications for cell cycle control and the targeted treatment of cancer. Biochim. Biophys. Acta Bioenerg. 2014, 1843, 150–162. [Google Scholar] [CrossRef] [Green Version]
  21. Bonvini, P.; Rosa, H.D.; Vignes, N.; Rosolen, A. Ubiquitination and proteasomal degradation of nucleophosmin-anaplastic lymphoma kinase induced by 17-allylamino-demethoxygeldanamycin: Role of the co-chaperone carboxyl heat shock protein 70-interacting protein. Cancer Res. 2004, 64, 3256–3264. [Google Scholar] [CrossRef] [Green Version]
  22. Didcock, L.; Young, D.F.; Goodbourn, S.; Randall, R.E. The V Protein of Simian Virus 5 Inhibits Interferon Signalling by Targeting STAT1 for Proteasome-Mediated Degradation. J. Virol. 1999, 73, 9928–9933. [Google Scholar] [CrossRef] [Green Version]
  23. Desai, S.D.; Li, T.-K.; Rodriguez-Bauman, A.; Rubin, E.H.; Liu, L.F. Ubiquitin/26S proteasome-mediated degradation of topoisomerase I as a resistance mechanism to camptothecin in tumor cells. Cancer Res. 2001, 61, 5926–5932. [Google Scholar]
  24. Ratner, J.N.; Balasubramanian, B.; Corden, J.; Warren, S.L.; Bregman, D.B. Ultraviolet Radiation-induced Ubiquitination and Proteasomal Degradation of the Large Subunit of RNA Polymerase II. J. Biol. Chem. 1998, 273, 5184–5189. [Google Scholar] [CrossRef] [Green Version]
  25. Rui, L.; Fisher, T.L.; Thomas, J.; White, M.F. Regulation of Insulin/Insulin-like Growth Factor-1 Signaling by Proteasome-mediated Degradation of Insulin Receptor Substrate-2. J. Biol. Chem. 2001, 276, 40362–40367. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  26. Sun, X.J.; Goldberg, J.L.; Qiao, L.Y.; Mitchell, J.J. Insulin-induced insulin receptor substrate-1 degradation is mediated by the proteasome degradation pathway. Diabetes 1999, 48, 1359–1364. [Google Scholar] [CrossRef] [PubMed]
  27. McDonald, E.R.; El-Deiry, W.S. Suppression of caspase-8- and -10-associated RING proteins results in sensitization to death ligands and inhibition of tumor cell growth. Proc. Natl. Acad. Sci. USA 2004, 101, 6170–6175. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  28. Suzuki, Y.; Nakabayashi, Y.; Takahashi, R. Ubiquitin-protein ligase activity of X-linked inhibitor of apoptosis protein promotes proteasomal degradation of caspase-3 and enhances its anti-apoptotic effect in Fas-induced cell death. Proc. Natl. Acad. Sci. USA 2001, 98, 8662–8667. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  29. Miller, D.M.; Thomas, S.D.; Islam, A.; Muench, D.; Sedoris, K. c-Myc and Cancer Metabolism. Clin. Cancer Res. 2012, 18, 5546–5553. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  30. Gregory, M.A.; Hann, S.R. c-Myc Proteolysis by the Ubiquitin-Proteasome Pathway: Stabilization of c-Myc in Burkitt’s Lymphoma Cells. Mol. Cell. Biol. 2000, 20, 2423–2435. [Google Scholar] [CrossRef] [Green Version]
  31. DeBerardinis, R.J.; Thompson, C.B. Cellular Metabolism and Disease: What Do Metabolic Outliers Teach Us? Cell 2012, 148, 1132–1144. [Google Scholar] [CrossRef] [Green Version]
  32. Kelly, J.M.; Summers, M.; Park, H.S.; Milligan, L.P.; McBride, B.W. Cellular Energy Metabolism and Regulation; Academic Press: New York, NY, USA, 1991; Volume 74, ISBN 978-0-12-066150-3. [Google Scholar]
  33. Bhoj, V.G.; Chen, Z.J. Ubiquitylation in innate and adaptive immunity. Nature 2009, 458, 430–437. [Google Scholar] [CrossRef]
  34. Zhang, M.; Wang, L.; Zhao, X.; Zhao, K.; Meng, H.; Zhao, W.; Gao, C. TRAF-interacting protein (TRIP) negatively regulates IFN-β production and antiviral response by promoting proteasomal degradation of TANK-binding kinase 1. J. Exp. Med. 2012, 209, 1703–1711. [Google Scholar] [CrossRef] [Green Version]
  35. Glickman, M.H.; Ciechanover, A. The Ubiquitin-Proteasome Proteolytic Pathway: Destruction for the Sake of Construction. Physiol. Rev. 2002, 82, 373–428. [Google Scholar] [CrossRef] [PubMed]
  36. Sakamoto, K.M. Ubiquitin-dependent proteolysis: Its role in human diseases and the design of therapeutic strategies. Mol. Genet. Metab. 2002, 77, 44–56. [Google Scholar] [CrossRef]
  37. Layfield, R.; Lowe, J.; Bedford, L. The ubiquitin–proteasome system and neurodegenerative disorders. Essays Biochem. 2005, 41, 157–171. [Google Scholar] [CrossRef] [PubMed]
  38. McKinnon, C.; Tabrizi, S. The Ubiquitin-Proteasome System in Neurodegeneration. Antioxid. Redox Signal. 2014, 21, 2302–2321. [Google Scholar] [CrossRef]
  39. Dou, Q.P.; Smith, D.M.; Daniel, K.G.; Kazi, A. Interruption of tumor cell cycle progression through proteasome inhibition: Implications for cancer therapy. Prog. Cell Cycle Res. 2003, 5, 441–446. [Google Scholar]
  40. Papandreou, C.N. The Proteasome as a Target for Cancer Treatment. Am. J. Cancer 2005, 4, 359–372. [Google Scholar] [CrossRef] [Green Version]
  41. Dick, L.R.; Fleming, P.E. Building on bortezomib: Second-generation proteasome inhibitors as anti-cancer therapy. Drug Discov. Today 2010, 15, 243–249. [Google Scholar] [CrossRef]
  42. Schmidt, M.; Finley, D. Regulation of proteasome activity in health and disease. Biochim. Biophys. Acta Mol. Cell Res. 2014, 1843, 13–25. [Google Scholar] [CrossRef] [Green Version]
  43. Zmuda, F.; Sastry, L.; Shepherd, S.M.; Jones, D.; Scott, A.; Craggs, P.D.; Cortes, A.; Gray, D.W.; Torrie, L.S.; De Rycker, M. Identification of Novel Trypanosoma cruzi Proteasome Inhibitors Using a Luminescence-Based High-Throughput Screening Assay. Antimicrob. Agents Chemother. 2019, 63, e00309-19. [Google Scholar] [CrossRef] [Green Version]
  44. Cromm, P.M.; Crews, C.M. The Proteasome in Modern Drug Discovery: Second Life of a Highly Valuable Drug Target. ACS Central Sci. 2017, 3, 830–838. [Google Scholar] [CrossRef] [Green Version]
  45. Aminake, M.N.; Arndt, H.-D.; Pradel, G. The proteasome of malaria parasites: A multi-stage drug target for chemotherapeutic intervention? Int. J. Parasitol. Drugs Drug Resist. 2012, 2, 1–10. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  46. Le Chapelain, C.; Groll, M. Rational Design of Proteasome Inhibitors as Antimalarial Drugs. Angew. Chem. Int. Ed. 2016, 55, 6370–6372. [Google Scholar] [CrossRef] [PubMed]
  47. Gandotra, S.; Schnappinger, D.; Monteleone, M.; Hillen, W.; Ehrt, S. In vivo gene silencing identifies the Mycobacterium tuberculosis proteasome as essential for the bacteria to persist in mice. Nat. Med. 2007, 13, 1515–1520. [Google Scholar] [CrossRef] [PubMed]
  48. Bochtler, M.; Ditzel, L.; Groll, M.; Huber, R. Crystal structure of heat shock locus V (HslV) from Escherichia coli. Proc. Natl. Acad. Sci. USA 1997, 94, 6070–6074. [Google Scholar] [CrossRef] [Green Version]
  49. Whitby, F.G.; Masters, E.I.; Kramer, L.; Knowlton, J.R.; Yao, Y.; Wang, C.C.; Hill, C.P. Structural basis for the activation of 20S proteasomes by 11S regulators. Nature 2000, 408, 115–120. [Google Scholar] [CrossRef] [PubMed]
  50. Groll, M.; Bajorek, M.; Köhler, A.; Moroder, L.; Rubin, D.M.; Huber, R.; Glickman, M.H.; Finley, D. A gated channel into the proteasome core particle. Nat. Struct. Mol. Biol. 2000, 7, 1062–1067. [Google Scholar] [CrossRef]
  51. Dick, T.P.; Nussbaum, A.K.; Deeg, M.; Heinemeyer, W.; Groll, M.; Schirle, M.; Keilholz, W.; Stevanović, S.; Wolf, D.H.; Huber, R.; et al. Contribution of Proteasomal β-Subunits to the Cleavage of Peptide Substrates Analyzed with Yeast Mutants. J. Biol. Chem. 1998, 273, 25637–25646. [Google Scholar] [CrossRef] [Green Version]
  52. Huber, E.M. Introduction. In Structural and Functional Characterization of the Immunoproteasome; Huber, E.M., Ed.; Springer Theses TS; Springer International Publishing: Cham, Switzerland, 2013; pp. 1–18. ISBN 978-3-319-01555-2. [Google Scholar]
  53. Dou, Q.P.; Goldfarb, R.H. Bortezomib (millennium pharmaceuticals). IDrugs Investig. Drugs J. 2002, 5, 828–834. [Google Scholar]
  54. Kane, R.C.; Bross, P.F.; Farrell, A.T.; Pazdur, R. Velcade®: U.S. FDA Approval for the Treatment of Multiple Myeloma Progressing on Prior Therapy. Oncologist 2003, 8, 508–513. [Google Scholar] [CrossRef]
  55. Siegel, D.S.; Martin, T.; Wang, M.; Vij, R.; Jakubowiak, A.J.; Lonial, S.; Trudel, S.; Kukreti, V.; Bahlis, N.; Alsina, M.; et al. A phase 2 study of single-agent carfilzomib (PX-171-003-A1) in patients with relapsed and refractory multiple myeloma. Blood 2012, 120, 2817–2825. [Google Scholar] [CrossRef] [Green Version]
  56. Parlati, F.; Lee, S.J.; Aujay, M.; Suzuki, E.; Levitsky, K.; Lorens, J.B.; Micklem, D.R.; Ruurs, P.; Sylvain, C.; Lu, Y.; et al. Carfilzomib can induce tumor cell death through selective inhibition of the chymotrypsin-like activity of the proteasome. Blood 2009, 114, 3439–3447. [Google Scholar] [CrossRef] [Green Version]
  57. Moreau, P.; Masszi, T.; Grzasko, N.; Bahlis, N.J.; Hansson, M.; Pour, L.; Sandhu, I.; Ganly, P.; Baker, B.W.; Jackson, S.R.; et al. Oral Ixazomib, Lenalidomide, and Dexamethasone for Multiple Myeloma. N. Engl. J. Med. 2016, 374, 1621–1634. [Google Scholar] [CrossRef] [PubMed]
  58. Mizukami, T.; Asai, A.; Yamashita, Y.; Katahira, R.; Hasegawa, A.; Ochiai, K.; Akinaga, S. UCK 14 Compounds. U.S. Patent 5,663, 298, 1997. Available online: https://patentimages.storage.googleapis.com/53/ba/fc/d9b04aa7e971dd/US5663298.pdf (accessed on 1 March 2022).
  59. Asai, A.; Hasegawa, A.; Ochiai, K.; Yamashita, Y.; Mizukami, T. Belactosin A, a Novel Antitumor Antibiotic Acting on Cyclin/CDK Mediated Cell Cycle Regulation, Produced by Streptomyces sp. J. Antibiot. 2000, 53, 81–83. [Google Scholar] [CrossRef] [Green Version]
  60. Asai, A.; Tsujita, T.; Sharma, S.V.; Yamashita, Y.; Akinaga, S.; Funakoshi, M.; Kobayashi, H.; Mizukami, T. A new structural class of proteasome inhibitors identified by microbial screening using yeast-based assay. Biochem. Pharmacol. 2004, 67, 227–234. [Google Scholar] [CrossRef] [PubMed]
  61. Zhu, M.; Harshbarger, W.D.; Robles, O.; Krysiak, J.; Hull, K.G.; Cho, S.W.; Richardson, R.D.; Yang, Y.; Garcia, A.; Spiegelman, L.; et al. A strategy for dual inhibition of the proteasome and fatty acid synthase with belactosin C-orlistat hybrids. Bioorganic Med. Chem. 2017, 25, 2901–2916. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  62. Gillessen, S.; Groettrup, M.; Cerny, T. The Proteasome, a New Target for Cancer Therapy. Onkologie 2002, 25, 534–539. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  63. Almond, J.B.; Cohen, G.M. The Proteasome: A Target Novel for Cancer Chemotherapy. Leukemia 2002, 16, 433–443. [Google Scholar] [CrossRef] [Green Version]
  64. Elliott, P.J.; Zollner, T.M.; Boehncke, W.-H. Proteasome inhibition: A new anti-inflammatory strategy. Klin. Wochenschr. 2003, 81, 235–245. [Google Scholar] [CrossRef]
  65. Armstrong, A.; Scutt, J.N. Total synthesis of (+)-belactosin A. Chem. Commun. 2004, 4, 510–511. [Google Scholar] [CrossRef] [Green Version]
  66. Groll, M.; Larionov, O.V.; Huber, R.; de Meijere, A. Inhibitor-binding mode of homobelactosin C to proteasomes: New insights into class I MHC ligand generation. Proc. Natl. Acad. Sci. USA 2006, 103, 4576–4579. [Google Scholar] [CrossRef] [Green Version]
  67. Knight, Z.A.; Lin, H.; Shokat, K.M. Targeting the cancer kinome through polypharmacology. Nat. Cancer 2010, 10, 130–137. [Google Scholar] [CrossRef] [PubMed]
  68. Ghosh, A.K.; Samanta, I.; Mondal, A.; Liu, W.R. Covalent Inhibition in Drug Discovery. ChemMedChem 2019, 14, 889–906. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  69. Baillie, T.A. Targeted Covalent Inhibitors for Drug Design. Angew. Chem. Int. Ed. 2016, 55, 13408–13421. [Google Scholar] [CrossRef]
  70. Huggins, D.J.; Sherman, W.; Tidor, B. Rational Approaches to Improving Selectivity in Drug Design. J. Med. Chem. 2012, 55, 1424–1444. [Google Scholar] [CrossRef]
  71. Eisenberg, D.; Schwarz, E.; Komaromy, M.; Wall, R. Analysis of membrane and surface protein sequences with the hydrophobic moment plot. J. Mol. Biol. 1984, 179, 125–142. [Google Scholar] [CrossRef]
  72. Aparicio, N.S.; Świderek, K.; Moliner, V. Theoretical study of the inhibition mechanism of human 20S proteasome by dihydroeponemycin. Eur. J. Med. Chem. 2019, 164, 399–407. [Google Scholar] [CrossRef] [PubMed]
  73. Serrano-Aparicio, N.; Moliner, V.; Świderek, K. Nature of Irreversible Inhibition of Human 20S Proteasome by Salinosporamide A. The Critical Role of Lys–Asp Dyad Revealed from Electrostatic Effects Analysis. ACS Catal. 2021, 11, 3575–3589. [Google Scholar] [CrossRef]
  74. Huber, E.M.; Heinemeyer, W.; Li, X.; Arendt, C.S.; Hochstrasser, M.; Groll, M. A unified mechanism for proteolysis and autocatalytic activation in the 20S proteasome. Nat. Commun. 2016, 7, 10900. [Google Scholar] [CrossRef] [Green Version]
  75. Groll, M.; Huber, R.; Potts, B.C.M. Crystal Structures of Salinosporamide A (NPI-0052) and B (NPI-0047) in Complex with the 20S Proteasome Reveal Important Consequences of β-Lactone Ring Opening and a Mechanism for Irreversible Binding. J. Am. Chem. Soc. 2006, 128, 5136–5141. [Google Scholar] [CrossRef]
  76. Serrano-Aparicio, N.; Moliner, V.; Świderek, K. On the Origin of the Different Reversible Characters of Salinosporamide A and Homosalinosporamide A in the Covalent Inhibition of the Human 20S Proteasome. ACS Catal. 2021, 11, 11806–11819. [Google Scholar] [CrossRef]
  77. Kim, D.H.; Park, J.-I.; Chung, S.J.; Park, J.D.; Park, N.-K.; Han, J.H. Cleavage of β-lactone ring by serine protease. Mechanistic implications. Bioorg. Med. Chem. 2002, 10, 2553–2560. [Google Scholar] [CrossRef]
  78. Galmés, M.A.; García-Junceda, E.; Świderek, K.; Moliner, V. Exploring the Origin of Amidase Substrate Promiscuity in CALB by a Computational Approach. ACS Catal. 2020, 10, 1938–1946. [Google Scholar] [CrossRef] [Green Version]
  79. Burgi, H.B.; Dunitz, J.D.; Shefter, E. Geometrical reaction coordinates. II. Nucleophilic addition to a carbonyl group. J. Am. Chem. Soc. 1973, 95, 5065–5067. [Google Scholar] [CrossRef]
  80. Burgi, H.B.; Dunitz, J.D.; Lehn, J.; Wipff, G. Stereochemistry of reaction paths at carbonyl centres. Tetrahedron 1974, 30, 1563–1572. [Google Scholar] [CrossRef]
  81. Fleming, I. Molecular Orbitals and Organic Chemical Reactions, Student Edition; John Wiley & Sons: Hoboken NJ, USA, 2011; pp. 214–215. ISBN 9780470746608. [Google Scholar]
  82. Heathcock, C.H. Understanding and Controlling Diastereofacial Selectivity in Carbon-Carbon Bond-Forming Reactions. Aldrichimica Acta 1990, 23, 94–111. [Google Scholar]
  83. Schrader, J.; Henneberg, F.; Mata, R.A.; Tittmann, K.; Schneider, T.R.; Stark, H.; Bourenkov, G.; Chari, A. The inhibition mechanism of human 20 S proteasomes enables next-generation inhibitor design. Science 2016, 353, 594–598. [Google Scholar] [CrossRef] [Green Version]
  84. Wang, J.; Wang, W.; Kollman, P.A.; Case, D.A. Automatic atom type and bond type perception in molecular mechanical calculations. J. Mol. Graph. Model. 2006, 25, 247–260. [Google Scholar] [CrossRef]
  85. Søndergaard, C.R.; Olsson, M.H.M.; Rostkowski, M.; Jensen, J.H. Improved Treatment of Ligands and Coupling Effects in Empirical Calculation and Rationalization of pKa Values. J. Chem. Theory Comput. 2011, 7, 2284–2295. [Google Scholar] [CrossRef]
  86. Olsson, M.H.M.; Søndergaard, C.R.; Rostkowski, M.; Jensen, J.H. PROPKA3: Consistent Treatment of Internal and Surface Residues in Empirical pKa Predictions. J. Chem. Theory Comput. 2011, 7, 525–537. [Google Scholar] [CrossRef]
  87. Schafmeister, C.E.A.; Ross, W.S.; Romanovski, V. LEAP; University of California: San Francisco, CA, USA, 1995. [Google Scholar]
  88. Jorgensen, W.L.; Chandrasekhar, J.; Madura, J.D.; Impey, R.W.; Klein, M.L. Comparison of simple potential functions for simulating liquid water. J. Chem. Phys. 1983, 79, 926–935. [Google Scholar] [CrossRef]
  89. Duan, Y.; Wu, C.; Chowdhury, S.; Lee, M.C.; Xiong, G.; Zhang, W.; Yang, R.; Cieplak, P.; Luo, R.; Lee, T.; et al. A point-charge force field for molecular mechanics simulations of proteins based on condensed-phase quantum mechanical calculations. J. Comput. Chem. 2003, 24, 1999–2012. [Google Scholar] [CrossRef]
  90. Phillips, J.C.; Braun, R.; Wang, W.; Gumbart, J.; Tajkhorshid, E.; Villa, E.; Chipot, C.; Skeel, R.D.; Kalé, L.; Schulten, K. Scalable molecular dynamics with NAMD. J. Comput. Chem. 2005, 26, 1781–1802. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  91. Giovannini, T.; Lafiosca, P.; Cappelli, C. A General Route to Include Pauli Repulsion and Quantum Dispersion Effects in QM/MM Approaches. J. Chem. Theory Comput. 2017, 13, 4854–4870. [Google Scholar] [CrossRef] [PubMed]
  92. Giovannini, T.; Lafiosca, P.; Chandramouli, B.; Barone, V.; Cappelli, C. Effective yet reliable computation of hyperfine coupling constants in solution by a QM/MM approach: Interplay between electrostatics and non-electrostatic effects. J. Chem. Phys. 2019, 150, 124102. [Google Scholar] [CrossRef] [PubMed]
  93. Lennard-Jones, J.E. Cohesion. Proc. Phys. Soc. 1931, 43, 461–482. [Google Scholar] [CrossRef]
  94. Grest, G.S.; Kremer, K. Molecular dynamics simulation for polymers in the presence of a heat bath. Phys. Rev. A 1986, 33, 3628–3631. [Google Scholar] [CrossRef]
  95. Martyna, G.J.; Tobias, D.J.; Klein, M.L. Constant pressure molecular dynamics algorithms. J. Chem. Phys. 1994, 101, 4177–4189. [Google Scholar] [CrossRef]
  96. Field, M.J.; Albe, M.; Bret, C.; Martin, F.P.-D.; Thomas, A. The dynamo library for molecular simulations using hybrid quantum mechanical and molecular mechanical potentials. J. Comput. Chem. 2000, 21, 1088–1100. [Google Scholar] [CrossRef]
  97. Krzemińska, A.; Paneth, P.; Moliner, V.; Świderek, K. Binding Isotope Effects as a Tool for Distinguishing Hydrophobic and Hydrophilic Binding Sites of HIV-1 RT. J. Phys. Chem. B 2014, 119, 917–927. [Google Scholar] [CrossRef]
  98. Field, M.J.; Bash, P.A.; Karplus, M. A combined quantum mechanical and molecular mechanical potential for molecular dynamics simulations. J. Comput. Chem. 1990, 11, 700–733. [Google Scholar] [CrossRef]
  99. Dewar, M.J.S.; Zoebisch, E.G.; Healy, E.F.; Stewart, J.J.P. Development and use of quantum mechanical molecular models. 76. AM1: A new general purpose quantum mechanical molecular model. J. Am. Chem. Soc. 1985, 107, 3902–3909. [Google Scholar] [CrossRef]
  100. Zhao, Y.; Truhlar, D.G. The M06 suite of density functionals for main group thermochemistry, thermochemical kinetics, noncovalent interactions, excited states, and transition elements: Two new functionals and systematic testing of four M06-class functionals and 12 other functionals. Theor. Chem. Acc. 2008, 120, 215–241. [Google Scholar] [CrossRef] [Green Version]
  101. Baker, J. An algorithm for the location of transition states. J. Comput. Chem. 1986, 7, 385–395. [Google Scholar] [CrossRef]
  102. Torrie, G.M.; Valleau, J.P. Nonphysical sampling distributions in Monte Carlo free-energy estimation: Umbrella sampling. J. Comput. Phys. 1977, 23, 187–199. [Google Scholar] [CrossRef]
  103. Kumar, S.; Rosenberg, J.M.; Bouzida, D.; Swendsen, R.; Kollman, P.A. THE weighted histogram analysis method for free-energy calculations on biomolecules. I. The method. J. Comput. Chem. 1992, 13, 1011–1021. [Google Scholar] [CrossRef]
  104. Ruiz-Pernía, J.J.; Silla, E.; Tuñón, I.; Martí, S. Hybrid Quantum Mechanics/Molecular Mechanics Simulations with Two-Dimensional Interpolated Corrections: Application to Enzymatic Processes. J. Phys. Chem. B 2006, 110, 17663–17670. [Google Scholar] [CrossRef]
  105. Chuang, Y.-Y.; Corchado, A.J.C.; Truhlar, D.G. Mapped Interpolation Scheme for Single-Point Energy Corrections in Reaction Rate Calculations and a Critical Evaluation of Dual-Level Reaction Path Dynamics Methods. J. Phys. Chem. A 1999, 103, 1140–1149. [Google Scholar] [CrossRef]
  106. Świderek, K.; Tuñón, I.; Martí, S.; Moliner, V. Protein Conformational Landscapes and Catalysis. Influence of Active Site Conformations in the Reaction Catalyzed by L-Lactate Dehydrogenase. ACS Catal. 2015, 5, 1172–1185. [Google Scholar] [CrossRef] [Green Version]
  107. Frisch, M.J.; Trucks, G.W.; Schlegel, H.B.; Scuseria, G.E.; Robb, M.A.; Cheeseman, J.R.; Scalmani, G.; Barone, V.; Petersson, G.A.; Nakatsuji, H.; et al. Gaussian 09; Gaussian, Inc.: Wallingford, CT, USA, 2009; Available online: https://gaussian.com/glossary/g09/ (accessed on 1 March 2022).
Figure 1. Schematic representation of constitutive 20S proteasome with the indicated position of catalytically active subunits in β-rings.
Figure 1. Schematic representation of constitutive 20S proteasome with the indicated position of catalytically active subunits in β-rings.
Pharmaceuticals 15 00531 g001
Figure 2. Structure of (−)-Belactosin C and bisbenzyl-protected N-Cbz-O-Bn-(−)-Homobelactosin C (A) (where Bn stands for a benzyl group and Cbz for benzyl carbamates) together with the general scheme for the covalent inhibition reaction (B).
Figure 2. Structure of (−)-Belactosin C and bisbenzyl-protected N-Cbz-O-Bn-(−)-Homobelactosin C (A) (where Bn stands for a benzyl group and Cbz for benzyl carbamates) together with the general scheme for the covalent inhibition reaction (B).
Pharmaceuticals 15 00531 g002
Figure 3. Structure of the binding site of N-Cbz-O-Bn-(−)-homobelactosin C in human 20S Proteasome with indicated positions of the Thr1 and Lys33, catalytically important residues in the active site and five key subunits β6, β5, β4, β3′and β4′ that participate in its formation (on the (left)) together with its hydrophobicity map generated based on the Eisenberg hydrophobicity scale (on the (right)).
Figure 3. Structure of the binding site of N-Cbz-O-Bn-(−)-homobelactosin C in human 20S Proteasome with indicated positions of the Thr1 and Lys33, catalytically important residues in the active site and five key subunits β6, β5, β4, β3′and β4′ that participate in its formation (on the (left)) together with its hydrophobicity map generated based on the Eisenberg hydrophobicity scale (on the (right)).
Pharmaceuticals 15 00531 g003
Figure 4. (A) Eleven snapshots illustrating a variety of conformations explored by hBelC in the binding cavity of 20S proteasome during the last 50 ns of MD simulations. (B) The root-mean-square-deviation computed for positions of heavy atoms of hBelC in its covalent complex with 20S proteasome during 100 ns of classical MD simulations. (C) H-bond contacts established between hBelC and residues from the binding pocket during the last 50 ns of simulations defined using geometric criteria (3.0 Å for Donor···Acceptor distance, and 135.0° for Donor-H⋯Acceptor angle) (D) Interaction energies computed for hBelC and residues of S2′ pocket as well as residues of the oxyanion hole.
Figure 4. (A) Eleven snapshots illustrating a variety of conformations explored by hBelC in the binding cavity of 20S proteasome during the last 50 ns of MD simulations. (B) The root-mean-square-deviation computed for positions of heavy atoms of hBelC in its covalent complex with 20S proteasome during 100 ns of classical MD simulations. (C) H-bond contacts established between hBelC and residues from the binding pocket during the last 50 ns of simulations defined using geometric criteria (3.0 Å for Donor···Acceptor distance, and 135.0° for Donor-H⋯Acceptor angle) (D) Interaction energies computed for hBelC and residues of S2′ pocket as well as residues of the oxyanion hole.
Pharmaceuticals 15 00531 g004
Figure 5. Schematic representation of the proposed three-step inhibition mechanism of β5-subunits of the 20S proteasome with N-Cbz-O-Bn-(−)-homobelactosin C.
Figure 5. Schematic representation of the proposed three-step inhibition mechanism of β5-subunits of the 20S proteasome with N-Cbz-O-Bn-(−)-homobelactosin C.
Pharmaceuticals 15 00531 g005
Figure 6. Free energy surfaces explored for all three proposed steps of the chemical transformations occurring along the 20S proteasome inhibition mechanism with N-Cbz-O-Bn-(−)-homobelactosin C. PMFs were obtained at AM1/AMBER level and corrected at M06-2X/AMBER level of theory with 6-31+G(d,p) basis set (see Computational Methods for details). Values of distinguished coordinates and energy are given in Å and in kcal·mol−1, respectively. Black stars indicate the positions of optimized and characterized transition state structures detected along the reaction progress at the M06-2X/AMBER level of theory with a 6-31+G(d,p) basis set.
Figure 6. Free energy surfaces explored for all three proposed steps of the chemical transformations occurring along the 20S proteasome inhibition mechanism with N-Cbz-O-Bn-(−)-homobelactosin C. PMFs were obtained at AM1/AMBER level and corrected at M06-2X/AMBER level of theory with 6-31+G(d,p) basis set (see Computational Methods for details). Values of distinguished coordinates and energy are given in Å and in kcal·mol−1, respectively. Black stars indicate the positions of optimized and characterized transition state structures detected along the reaction progress at the M06-2X/AMBER level of theory with a 6-31+G(d,p) basis set.
Pharmaceuticals 15 00531 g006
Figure 7. Optimized and characterized transition state structures along the minimum energy reaction pathway at the M06-2X/AMBER level of theory with a 6-31+(d,p) basis set. Key distances are in Å while values of imaginary wavenumbers are in cm−1.
Figure 7. Optimized and characterized transition state structures along the minimum energy reaction pathway at the M06-2X/AMBER level of theory with a 6-31+(d,p) basis set. Key distances are in Å while values of imaginary wavenumbers are in cm−1.
Pharmaceuticals 15 00531 g007
Figure 8. Free energy profile computed for 20S proteasome inhibition mechanism by N-Cbz-O-Bn-(−)-homobelactosin C (hBelC) (in red) and homo-Salinosporamide A (hSalA) (in gray). Results for hSalA were taken from [76].
Figure 8. Free energy profile computed for 20S proteasome inhibition mechanism by N-Cbz-O-Bn-(−)-homobelactosin C (hBelC) (in red) and homo-Salinosporamide A (hSalA) (in gray). Results for hSalA were taken from [76].
Pharmaceuticals 15 00531 g008
Figure 9. (A) Distribution of ESP charges on the β-lactone ring in hBelC and hSalA inhibitors computed using ChelpG method at M06-2X/6-31G+(d,p) level of theory with 6-31+G(d,p) basis set in gas phase. (B) Evolution of nucleophilic attack trajectory described by Bürgi–Dunitz (αBD) and Flippin–Lodge (αFL) angles along minimum energy pathway computed at M06-2X/6-31G+(d,p) level of theory for the first step of 20S proteasome inhibition with hBelC and hSalA.
Figure 9. (A) Distribution of ESP charges on the β-lactone ring in hBelC and hSalA inhibitors computed using ChelpG method at M06-2X/6-31G+(d,p) level of theory with 6-31+G(d,p) basis set in gas phase. (B) Evolution of nucleophilic attack trajectory described by Bürgi–Dunitz (αBD) and Flippin–Lodge (αFL) angles along minimum energy pathway computed at M06-2X/6-31G+(d,p) level of theory for the first step of 20S proteasome inhibition with hBelC and hSalA.
Pharmaceuticals 15 00531 g009
Figure 10. (A) Overlay of the optimized structures of E-I2 complexes formed between the active site of β5-subunit and hBelC or hSalA. For clarity, only atoms of warheads and Thr1 are shown. (B) The evolution of dihedral angles established between planes defined by positions of C1βlac-C4βlac-C3βlac -O2βlac and C1βlac-C4βlac-C3βlac-CMeβlac atoms of the β-lactone ring along the C1βlac-O2βlac bond cleavage pathway. (C) The active site of 20S proteasome in its E-I2 complex with hBelC. Positions of the residues with favorable interactions are shown as cyan spheres while the residues with unfavorable interactions are shown as pink spheres. G47 participating in the formation of the oxyanion hole is shown as an orange sphere. (D) Interaction energy (Eint) established between hBelC or hSalA inhibitors and residues from the vicinity of both warheads in E-I2 adduct together with electrostatic Eelec and EvdW contribution to the overall Eint values in hBelC. The average values of Eint were obtained based on 1000 structures generated during 100 ps of QM/MM MD simulations at AM1/AMBER level of theory. For clarity, only those residues for which the obtained difference in computed Eint for studies inhibitors was higher than 1 kcal·mol−1 are provided.
Figure 10. (A) Overlay of the optimized structures of E-I2 complexes formed between the active site of β5-subunit and hBelC or hSalA. For clarity, only atoms of warheads and Thr1 are shown. (B) The evolution of dihedral angles established between planes defined by positions of C1βlac-C4βlac-C3βlac -O2βlac and C1βlac-C4βlac-C3βlac-CMeβlac atoms of the β-lactone ring along the C1βlac-O2βlac bond cleavage pathway. (C) The active site of 20S proteasome in its E-I2 complex with hBelC. Positions of the residues with favorable interactions are shown as cyan spheres while the residues with unfavorable interactions are shown as pink spheres. G47 participating in the formation of the oxyanion hole is shown as an orange sphere. (D) Interaction energy (Eint) established between hBelC or hSalA inhibitors and residues from the vicinity of both warheads in E-I2 adduct together with electrostatic Eelec and EvdW contribution to the overall Eint values in hBelC. The average values of Eint were obtained based on 1000 structures generated during 100 ps of QM/MM MD simulations at AM1/AMBER level of theory. For clarity, only those residues for which the obtained difference in computed Eint for studies inhibitors was higher than 1 kcal·mol−1 are provided.
Pharmaceuticals 15 00531 g010
Figure 11. Distribution of distances established between O2βlac, NhBelC, and C1′hBelC to C1βlac together with the distribution of Bürgi–Dunitz (αBD) and Flippin-Lodge (αFL) angles obtained for 10,000 snapshots generated during 100 ns MD simulations on E-PC complex.
Figure 11. Distribution of distances established between O2βlac, NhBelC, and C1′hBelC to C1βlac together with the distribution of Bürgi–Dunitz (αBD) and Flippin-Lodge (αFL) angles obtained for 10,000 snapshots generated during 100 ns MD simulations on E-PC complex.
Pharmaceuticals 15 00531 g011
Figure 12. Schematic representation of the active site of β5 subunit of human 20S proteasome with hBelC. The grey area includes atoms described at the QM level of theory. The position of link atoms is depicted as three black dots.
Figure 12. Schematic representation of the active site of β5 subunit of human 20S proteasome with hBelC. The grey area includes atoms described at the QM level of theory. The position of link atoms is depicted as three black dots.
Pharmaceuticals 15 00531 g012
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Serrano-Aparicio, N.; Ferrer, S.; Świderek, K. Covalent Inhibition of the Human 20S Proteasome with Homobelactosin C Inquired by QM/MM Studies. Pharmaceuticals 2022, 15, 531. https://doi.org/10.3390/ph15050531

AMA Style

Serrano-Aparicio N, Ferrer S, Świderek K. Covalent Inhibition of the Human 20S Proteasome with Homobelactosin C Inquired by QM/MM Studies. Pharmaceuticals. 2022; 15(5):531. https://doi.org/10.3390/ph15050531

Chicago/Turabian Style

Serrano-Aparicio, Natalia, Silvia Ferrer, and Katarzyna Świderek. 2022. "Covalent Inhibition of the Human 20S Proteasome with Homobelactosin C Inquired by QM/MM Studies" Pharmaceuticals 15, no. 5: 531. https://doi.org/10.3390/ph15050531

APA Style

Serrano-Aparicio, N., Ferrer, S., & Świderek, K. (2022). Covalent Inhibition of the Human 20S Proteasome with Homobelactosin C Inquired by QM/MM Studies. Pharmaceuticals, 15(5), 531. https://doi.org/10.3390/ph15050531

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop