The active sites of 3 types of aspartic proteases are modeled, based on crystallographic coordinates of endothiapepsin and of a model of HIV-l protease. The enthalpies of deprotonation from neutral to mono-anion and to dianion are calculated with semiempirical minimal neglect of differential overlap, hydrogen bonding corrected (MNDO/H). This quantum mechanical study of models for the active sites of pepsins, human renin and retroviral aspartic proteases demonstrates that the replacements of Thr-218 from pepsins by Ala in human renin and of both Ser-35 and Thr-218 by alanines in retroviral proteases increases the proton affinity and modulates the charge distribution of those active sites compared to the pepsins.Protease, Aspartic; HIV-l; Renin; Pepsin; Model, molecular
I~RODUCTIONAPs are a family of endopeptidases which have many common features [l J in addition to structural homology [2]. The catalytic machinery includes two closely positioned aspartic acids with highly conserved neighbor residues. Most APs are maximally active at low pH ('acid proteases'). They have diverse substrate specificities but are inhibited, to varying extents, by a common inhibitor, pepstatin A [3-51. APs are known to be involved in a few disease conditions such as AIDS, essential h~rtension, several cancers and others f6] and are now a major target for designing inhibitors.Recent research establishes that HIV-I protease, which processes the gag and pol polyprotein precursors towards the production of mature viruses, is an AP [7]. Like other retroviral proteases, it has a much shorter polypeptide chain than pepsins and renins, and was postulated to function in dimeric form [B]. Its crystal structure [9,10] confirms those predictions. Other divergences of HIV-l protease are the higher pH for optimal activity f7,11] compared to pepsins, and its much weaker inhibition by pepstatin A [4,7,12] Fig.1 shows the alignment of 8 active site residues which are crucial for maintaining the hydrogen-bonded structure of the catalytic centers. The role of the fourth residue of each partial sequence (which is mutated in h. renin and retroviral proteases compared to pepsins) in regulating the properties of active site aspartates has rarely been addressed and the source of optimal pH variations among the AP is poorly understood [20]. The conservation of spatial relations of those residues in the crystal structures of many native and inhibited enzymes supports the conclusion that it is possible to model the active site of an AP for which only primary structure is known, by using the published coordinates of another. We have taken this approach to study the active sites of h. renin and retroviral proteases based on the published coordinates of endothiapepsin. Models of renin were previously constructed for studying their interactions with inhibitors. However, some details of the structure are missing due to the lack of hydrogen coordinates in X-ray studies of proteins. Some of those proton positions are extremely important for hydrogen bonding.