Discussion of the catalytic pathway of cysteine proteases based on AM1 calculations

Unrevealing the Proteolytic Activity of RgpB Gingipain from Computational Simulations

Alzheimer's disease represents one of the greatest medical concerns for today's population and health services. Its multifactorial inherent nature represents a challenge for its treatment and requires the development of a broad spectrum of drugs. Recently, the cysteine protease gingipain RgpB has been related to neurodegenerative diseases, including Alzheimer's disease, and its inhibition appears to be a promising neuroprotective strategy. Given these features, a computational study that integrates molecular dynamics (MD) simulations with classical and hybrid quantum mechanics/molecular mechanics (QM/MM) potentials was carried out to unravel the atomistic details of RgpB activity. First, a preliminary study based on principal component analysis (PCA), determined the protonation state of the Cys/His catalytic dyad, as well as the crucial role of a flexible loop that favors reactive interactions of the catalytic residues and the peptide in the precatalytic state in its closed conformation. Then, different mechanisms were explored by means of QM/MM MD simulations. The most favorable mechanism consists of two stages. First is an acylation stage that takes place in two steps where, initially, the sulfur atom of the C244 residue attacks the carbonylic carbon of the peptide and the proton of the C244 residue is transferred to the amino group of the peptide in a concerted manner. Subsequently, the peptide bond is broken, and a fragment of the peptide is released. After that, the deacylation stage takes place in a single step where a water molecule attacks the carbonylic carbon of the peptide and a proton of the water is transferred to the C244 residue. The free energy barrier of the rate limiting step is in very good agreement with available experimental data. The mechanism exhibits an unusual role of H211 residue compared with other cysteine proteases but a crucial role of the peptide in triggering the catalysis. Notably, the atomic and energetic particularities found represent a significant contribution to the comprehension of the reaction mechanism and a great opportunity for the design of efficient inhibitors of gingipain RgpB.

Modeling the archetype cysteine protease reaction using dispersion corrected density functional methods in ONIOM-type hybrid QM/MM calculations; the proteolytic reaction of papain

A proteolytic reaction of papain with a simple peptide model substrate N-methylacetamide has been studied. Our aim was twofold: (i) we proposed a plausible reaction mechanism with the aid of potential energy surface scans and second geometrical derivatives calculated at the stationary points, and (ii) we investigated the applicability of the dispersion corrected density functional methods in comparison with the popular hybrid generalized gradient approximations (GGA) method (B3LYP) without such a correction in the QM/MM calculations for this particular problem. In the resting state of papain the ion pair and neutral forms of the Cys-His catalytic dyad have approximately the same energy and they are separated by only a small barrier. Zero point vibrational energy correction shifted this equilibrium slightly to the neutral form. On the other hand, the electrostatic solvation free energy corrections, calculated using the Poisson-Boltzmann method for the structures sampled from molecular dynamics simulation trajectories, resulted in a more stable ion-pair form. All methods we applied predicted at least a two elementary step acylation process via a zwitterionic tetrahedral intermediate. Using dispersion corrected DFT methods the thioester S-C bond formation and the proton transfer from histidine occur in the same elementary step, although not synchronously. The proton transfer lags behind (or at least does not precede) the S-C bond formation. The predicted transition state corresponds mainly to the S-C bond formation while the proton is still on the histidine Nδ atom. In contrast, the B3LYP method using larger basis sets predicts a transition state in which the S-C bond is almost fully formed and the transition state can be mainly featured by the Nδ(histidine) to N(amid) proton transfer. Considerably lower activation energy was predicted (especially by the B3LYP method) for the next amide bond breaking elementary step of acyl-enzyme formation. Deacylation appeared to be a single elementary step process in all the methods we applied.

Distinct Roles of Catalytic Cysteine and Histidine in the Protease and Ligase Mechanisms of Human Legumain As Revealed by DFT-Based QM/MM Simulations

The cysteine protease enzyme legumain hydrolyzes peptide bonds with high specificity after asparagine and under more acidic conditions after aspartic acid [BakerE. N.J. Mol. Biol.1980, 141, 441−4847003158; BakerE. N.; J. Mol. Biol.1977, 111, 207–210859183; DrenthJ.; Biochemistry1976, 15, 3731–3738952885; MenardR.; J. Cell. Biochem.1994, 137; PolgarL.Eur. J. Biochem.1978, 88, 513–521689035; StorerA. C.; Methods Enzymol.1994, 244, 486–5007845227. Remarkably, legumain additionally exhibits ligase activity that prevails at pH > 5.5. The atomic reaction mechanisms including their pH dependence are only partly understood. Here we present a density functional theory (DFT)-based quantum mechanics/molecular mechanics (QM/MM) study of the detailed reaction mechanism of both activities for human legumain in solution. Contrasting the situation in other papain-like proteases, our calculations reveal that the active site Cys189 must be present in the protonated state for a productive nucleophilic attack and simultaneous rupture of the scissile peptide bond, consistent with the experimental pH profile of legumain-catalyzed cleavages. The resulting thioester intermediate (INT1) is converted by water attack on the thioester into a second intermediate, a diol (INT2), which is released by proton abstraction by Cys189. Surprisingly, we found that ligation is not the exact reverse of the proteolysis but can proceed via two distinct routes. Whereas the transpeptidation route involves aminolysis of the thioester (INT1), at pH 6 a cysteine-independent, histidine-assisted ligation route was found. Given legumain’s important roles in immunity, cancer, and neurodegenerative diseases, our findings open up possibilities for targeted drug design in these fields.

Conformational interconversions in peptide beta-turns: analysis of turns in proteins and computational estimates of barriers

The two most important beta-turn features in peptides and proteins are the type I and type II turns, which differ mainly in the orientation of the central peptide unit. Facile conformational interconversion is possible, in principle, by a flip of the central peptide unit. Homologous crystal structures afford an opportunity to structurally characterize both possible conformational states, thus allowing identification of sites that are potentially stereochemically mobile. A representative data set of 250 high-resolution (</=2.0 A), non-homologous protein crystal structures and corresponding variant and homologous entries, obtained from the Brookhaven Protein Data Bank, was examined to identify turns that are assigned different conformational types (type I/type II) in related structures. A total of 55 examples of beta-turns were identified as possible candidates for a stereochemically mobile site. Of the 55 examples, 45 could be classified as a potential site for interconversion between type I and type II beta-turns, while ten correspond to flips from type I' to type II' structures. As a further check, the temperature factors of the central peptide unit carbonyl oxygen atom of the 55 examples were examined. The analysis reveals that the turn assignments are indeed reliable. Examination of the secondary structures at the flanking positions of the flippable beta-turns reveals that seven examples occur in the loop region of beta-hairpins, indicating that the formation of ordered secondary structures on either side of the beta-turn does not preclude local conformational variations. In these beta-turns, Pro (11 examples), Lys (nine examples) and Ser (seven examples) were most often found at the i+1 position. Glycine was found to occur overwhelmingly at position i+2 (28 examples), while Ser (seven examples) and Asn (six examples) were amongst the most frequent residues. Activation energy barriers for the interconversion between type I and type II beta-turns were computed using the peptide models Ac-Pro-Aib-NHMe and Ac-Pro-Gly-NHMe within the framework of the AM1 semi-empirical molecular orbital procedure. In order to have a uniform basis for comparison and to eliminate the distracting influence of the deviation of backbone dihedral angles from that expected for ideal beta-turns, the dihedral angles phii+1 and psii+2 were fixed at the ideal values (phii+1=-60 degrees and psii+2=0 degrees). The other two angles (psii+1 and phii+2) were varied systematically to go from type II to type I beta-turn structures. The computational results suggest that there exists one stereospecific, concerted flip of the central peptide unit involving correlated single bond rotation that can occur with an activation barrier of the order of 3 kcal/mol. The results presented here suggest that conformational variations in beta-turns are observed in protein crystal structures and such changes may be an important dynamic feature in solution.