Binding modes of a new epoxysuccinyl-peptide inhibitor of cysteine proteases. Where and how do cysteine proteases express their selectivity?

An enormously active and selective azapeptide inhibitor of cathepsin B

The peptidomimetic Z-Arg-Leu-Arg-Agly-Ile-Val-OMe (where Agly means alpha-aza-glycyl, -NHNHCO-) is the strongest (K(i) = 480 pM) and the most selective inhibitor of cathepsin B to date, being approximately 2310 times as active to cathepsin B as to cathepsin K. In this paper we introduce the peptide and seek to rationalize its structure-activity relationships using molecular dynamics (MD) and NMR. It is shown that the -Agly-moiety restrains the peptide backbone to a bent shape, contrary to its parent peptide (with Gly in position 4), having its backbone extended and flexible. This fold is maintained in the plug covalently bound to the cathepsin B Cys29, in compliance with similar bends already observed in two other azapeptides attached to the active sites of cathepsin B. The MD simulation of the Z-Arg-Leu-Arg-Agly approximately cathepsin B complex suggests that, contrary to other potent inhibitors of cathepsin B, the current double Arg(1)/Arg(3) inhibitor, while maintaining the fold is able to form a unique ion cluster involving both Arg residues on the inhibitor part and two acidic Glu171 and Glu245 on the cathepsin B part, thus enhancing the affinity and subsequently the inhibiting power and selectivity of Z-Arg-Leu-Arg-Agly-Ile-Val-OMe to the observed extreme extent.

Proregion of Acanthoscelides obtectus cysteine proteinase: a novel peptide with enhanced selectivity toward endogenous enzymes

Acanthoscelides obtectus is a devastating storage insect pest capable of causing severe bean crop losses. In order to maintain their own development, insect pest larvae feed continuously, synthesizing efficient digestive enzymes. Among them, cysteine proteinases (CPs) are commonly produced as inactive precursors (procysteines), requiring a cleavage of the peptide proregion to become active. The proregion fits tightly into the active site of procysteines, efficiently preventing their activity. In this report, a CP cDNA (cpao) was isolated from A. obtectus midgut larvae. In silico studies indicated that the complete CP sequence contains a hydrophobic signal peptide, a prodomain and a conserved catalytic region. Moreover, the encoding cDNA contains 963bp translating into a 321 residue protein, CPAo, which was expressed in E. coli, fused with thioredoxin. Enzymatic assays using the recombinant protein revealed that the enzyme was catalytically active, being able to cleave the synthetic substrate Z-Phe-Arg-7-AMC. Additionally, this report also focuses the cpao propeptide (PCPAo) subcloning and expression. The expressed propeptide efficiently inhibited CPAo, as well as digestive CP of other bean bruchids. Little or no activity was found against proteolytic enzymes of two other coleopterans: Rhyzopertha dominica and Anthonomus grandis. The data reported here indicate the possibility of endogenous propeptides as a novel strategy on bruchids control, which could be applicable to bean improvement programs.

Prediction of Peptide Retention at Different HPLC Conditions from Multiple Linear Regression Models

To quantitatively characterize the structure of a peptide and to predict its gradient retention time at given HPLC conditions three structural descriptors are used: (i) logarithm of the sum of retention times of the amino acids composing the peptide, log SumAA, (ii) logarithm of the van der Waals volume of the peptide, log VDW(Vol), (iii) and the logarithm of the peptide's calculated n-octanol-water partition coefficient, clog P. The log SumAA descriptor is obtained from empirical data for 20 natural amino acids, determined in a given HPLC system. The two other descriptors are calculated from the peptides' structural formulas using molecular modeling methods. The quantitative structure-retention relationships (QSRR), build by multiple linear regression, describe HPLC retention of peptide on a given chromatographic system on which the retention of the 20 amino acids was predetermined. A structurally diversified series of 98 peptides was employed. The predicted gradient retention times on several chromatographic systems were in good agreement with the experimental data. The QSRR equations, derived for a given system operated at variable gradient times and temperatures allowed for the prediction of peptide retention in that system. Matching the experimental HPLC retention to the theoretically predicted for a presumed peptide could facilitate original protein identification in proteomics. In conjunction with MS data, prediction of the retention time for a given peptide might be used to improve the confidence of peptide identifications and to increase the number of correctly identified peptides.

Azapeptides structurally based upon inhibitory sites of cystatins as potent and selective inhibitors of cysteine proteases

A series of azapeptides as potential inhibitors of cysteine proteases were synthesized. Their structures, based on the binding center of cystatins, contain an azaglycine residue (Agly) in place of the evolutionarily conserved glycine residue in the N-terminal part of the enzyme binding region of cystatins. Incorporation of Agly should lead to deactivation of the acyl-enzyme complex formed against nucleophilic attack by water molecules in the final step of peptide bond hydrolysis. The majority of synthesized azapeptides shows high inhibitory potency toward the investigated cysteine proteases, papain, cathepsin B, and cathepsin K. One of them, Z-Arg-Leu-Val-Agly-Ile-Val-OMe (compound 17), which contains in its sequence the amino acid residues from the N-terminal binding segment as well as the hydrophobic residues from the first binding loop of human cystatin C, proved to be a highly potent and selective inhibitor of cathepsin B. It inhibits cathepsin B with a K(i) value of 0.088 nM. To investigate the influence of the structure of compound 17 for its inhibitory properties, we determined its conformation by means of NMR studies and theoretical calculations. The Z-Arg-Leu-Val-Agly fragment, covalently linked to Cys29 of cathepsin B, was also developed and modeled, in the catalytic pocket of the enzyme, through a molecular dynamics approach, to analyze ligand-protein interactions in detail. Analysis of the simulation trajectories generated using the AMBER force field provided us with atomic-level understanding of the conformational variability of this inhibitor, which is discussed in the context of other experimental and theoretical data.

The investigation of the effects of counterions in protein dynamics simulations

Molecular simulations able to exactly represent solvated charged proteins are helpful in understanding protein dynamics, structure and function. In the present study we have used two different starting structures of papain (a typical, stable, globular protein of intermediate net charge) and different modeling procedures to evaluate some effects of counterions in simulations. A number of configurations have been generated and relaxed for each system by various combinations of constrained simulated annealing and molecular dynamics procedures, using the AMBER force field. The analysis of trajectories shows that the simulations of solvated proteins are moderately sensitive to the presence of counterions. However, this sensitivity is highly dependent on the starting model and different procedures of equilibration used. The neutralized systems tend to evince smaller root mean square deviations regardless of the system investigated and the simulation procedure used. The results of parameterized fitting of the simulated structures to the crystallographic data, giving quantitative measure of the total charge influence on the stability of various elements of the secondary structure, revealed a clear scatter of different reactions of various systems' secondary structures to counterions addition: some systems apparently were stabilized when neutralized, while the others were not. Thus, one cannot unequivocally state, despite consideration of specific simulation conditions, whether protein secondary structures are more stable when they have neutralized charges. This suggests that caution should be taken when claiming the stabilizing effect of counterions in simulations other than those involving small, unstable polypeptides or highly charged proteins.