Quantitative structure-activity relationships for a series of inhibitors of cruzain from Trypanosoma cruzi: Molecular modeling, CoMFA and CoMSIA studies

Human parasitic diseases are the foremost threat to human health and welfare around the world. Trypanosomiasis is a very serious infectious disease against which the currently available drugs are limited and not effective. Therefore, there is an urgent need for new chemotherapeutic agents. One attractive drug target is the major cysteine protease from Trypanosoma cruzi, cruzain. In the present work, comparative molecular field analysis (CoMFA) and comparative molecular similarity indices analysis (CoMSIA) studies were conducted on a series of thiosemicarbazone and semicarbazone derivatives as inhibitors of cruzain. Molecular modeling studies were performed in order to identify the preferred binding mode of the inhibitors into the enzyme active site, and to generate structural alignments for the three-dimensional quantitative structure-activity relationship (3D QSAR) investigations. Statistically significant models were obtained (CoMFA. r(2) = 0.96 and q(2) = 0.78; CoMSIA, r(2) = 0.91 and q(2) = 0.73), indicating their predictive ability for untested compounds. The models were externally validated employing a test set, and the predicted values were in good agreement with the experimental results. The final QSAR models and the information gathered from the 3D CoMFA and CoMSIA contour maps provided important insights into the chemical and structural basis involved in the molecular recognition process of this family of cruzain inhibitors, and should be useful for the design of new structurally related analogs with improved potency. (C) 2009 Elsevier Inc. All rights reserved.

Effectiveness of commercial inhibitors against subtype F HIV-1 protease

Subtype F wild type HIV protease has been kinetically characterized using six commercial inhibitors (amprenavir, indinavir, lopinavir, nelfinavir, ritonavir and saquinavir) commonly used for HIV/AIDS treatment, as well as inhibitor TL-3 and acetylpepstatin. We also obtained kinetic parameters for two multi-resistant proteases (one of subtype B and one of subtype F) harboring primary and secondary mutations selected by intensive treatment with ritonavir/nelfinavir. This newly obtained biochemical data shows that all six studied commercially available protease inhibitors are significantly less effective against subtype F HIV proteases than against HIV proteases of subtype B, as judged by increased K(i) and biochemical fitness (vitality) values. Comparison with previously reported kinetic values for subtype A and C HIV proteases show that subtype F wild type proteases are significantly less susceptible to inhibition. These results demonstrate that the accumulation of natural polymorphisms in subtype F proteases yields catalytically more active enzymes with a large degree of cross-resistance, which thus results in strong virus viability.

Comparative molecular field analysis of a series of inhibitors of HIV-1 protease

Several protease inhibitors have reached the world market in the last fifteen years, dramatically improving the quality of life and life expectancy of millions of HIV-infected patients. In spite of the tremendous research efforts in this area, resistant HIV-1 variants are constantly decreasing the ability of the drugs to efficiently inhibit the enzyme. As a consequence, inhibitors with novel frameworks are necessary to circumvent resistance to chemotherapy. In the present work, we have created 3D QSAR models for a series of 82 HIV-1 protease inhibitors employing the comparative molecular field analysis (CoMFA) method. Significant correlation coefficients were obtained (q(2) = 0.82 and r(2) = 0.97), indicating the internal consistency of the best model, which was then used to evaluate an external test set containing 17 compounds. The predicted values were in good agreement with the experimental results, showing the robustness of the model and its substantial predictive power for untested compounds. The final QSAR model and the information gathered from the CoMFA contour maps should be useful for the design of novel anti-HIV agents with improved potency.

Pharmacophore Model of Cruzain Inhibitors

Chagas disease, caused by the protozoan Trypanosoma cruzi, is one of the most serious amongst the so-called neglected diseases in Latin America, specially in Brazil. So far there has been no effective treatment for the chronic phase of this disease. Cruzain is a major cysteine protease of T cruzi and it is recognized as a valid target for Chagas disease chemotherapy. The mechanism of cruzain action is associated with the nucleophilic attack of an activated sulfur atom towards electrophilic groups. In this report, features of a putative pharmacophore model of the enzyme, developed as a virtual screening tool for the selection of potential cruzain inhibitors, are described. The final proposed model was applied to the ZINC v.7 database and afterwards experimentally validated by an enzymatic inhibition assay. One of the compounds selected by the model showed cruzain inhibition in the low micromolar range.

Integration of Ligand- and Target-Based Virtual Screening for the Discovery of Cruzain Inhibitors

A myriad of methods are available for virtual screening of small organic compound databases. In this study we have successfully applied a quantitative model of consensus measurements, using a combination of 3D similarity searches (ROCS and EON), Hologram Quantitative Structure Activity Relationships (HQSAR) and docking (FRED, FlexX, Glide and AutoDock Vina), to retrieve cruzain inhibitors from collected databases. All methods were assessed individually and then combined in a Ligand-Based Virtual Screening (LBVS) and Target-Based Virtual Screening (TBVS) consensus scoring, using Receiving Operating Characteristic (ROC) curves to evaluate their performance. Three consensus strategies were used: scaled-rank-by-number, rank-by-rank and rank-by-vote, with the most thriving the scaled-rank-by-number strategy, considering that the stiff ROC curve appeared to be satisfactory in every way to indicate a higher enrichment power at early retrieval of active compounds from the database. The ligand-based method provided access to a robust and predictive HQSAR model that was developed to show superior discrimination between active and inactive compounds, which was also better than ROCS and EON procedures. Overall, the integration of fast computational techniques based on ligand and target structures resulted in a more efficient retrieval of cruzain inhibitors with desired pharmacological profiles that may be useful to advance the discovery of new trypanocidal agents.

2D QSAR and similarity studies on cruzain inhibitors aimed at improving selectivity over cathepsin L

Hologram quantitative structure-activity relationships (HQSAR) were applied to a data set of 41 cruzain inhibitors. The best HQSAR model (Q(2) = 0.77; R-2 = 0.90) employing Surflex-Sim, as training and test sets generator, was obtained using atoms, bonds, and connections as fragment distinctions and 4-7 as fragment size. This model was then used to predict the potencies of 12 test set compounds, giving satisfactory predictive R-2 value of 0,88. The contribution maps obtained from the best HQSAR model are in agreement with the biological activities of the study compounds. The Trypanosoma cruzi cruzain shares high similarity with the mammalian homolog cathepsin L. The selectivity toward cruzam was checked by a database of 123 compounds, which corresponds to the 41 cruzain inhibitors used in the HQSAR model development plus 82 cathepsin L inhibitors. We screened these compounds by ROCS (Rapid Overlay of Chemical Structures), a Gaussian-shape volume overlap filter that can rapidly identify shapes that match the query molecule. Remarkably, ROCS was able to rank the first 37 hits as being only cruzain inhibitors. In addition, the area under the curve (AUC) obtained with ROCS was 0.96, indicating that the method was very efficient to distinguishing between cruzain and cathepsin L inhibitors. (c) 2007 Elsevier Ltd. All rights reserved.

Recombinant expression, localization and in vitro inhibition of midgut cysteine peptidase (Sl-CathL) from sugarcane weevil, Sphenophorus levis

A cDNA coding for a digestive cathepsin L, denominated Sl-CathL, was isolated from a cDNA library of Sphenophorus levis larvae, representing the most abundant EST (10.49%) responsible for proteolysis in the midgut. The open reading frame of 972 bp encodes a preproenzyme similar to midgut cathepsin L-like enzymes in other coleopterans. Recombinant Sl-CathL was expressed in Pichia pastoris, with molecular mass of about 42 kDa. The recombinant protein was catalytically activated at low pH and the mature enzyme of 39 kDa displayed thermal instability and maximal activity at 37 degrees C and pH 6.0. Immunocytochemical analysis revealed Sl-CathL production in the midgut epithelium and secretion from vesicles containing the enzyme into the gut lumen, confirming an important role for this enzyme in the digestion of the insect larvae. The expression profile identified by RT-PCR through the biological cycle indicates that Sl-CathL is mainly produced in larval stages, with peak expression in 30-day-old larvae. At this stage, the enzyme is 1250-fold more expressed than in the pupal fase, in which the lowest expression level is detected. This enzyme is also produced in the adult stage, albeit in lesser abundance, assuming the presence of a different array of enzymes in the digestive system of adults. Tissue-specific analysis revealed that Sl-CathL mRNA synthesis occurs fundamentally in the larval midgut, thereby confirming its function as a digestive enzyme, as detected in immunolocalization assays. The catalytic efficiency of the purified recombinant enzyme was calculated using different substrates (Z-Leu-Arg-AMC, Z-Arg-Arg-AMC and Z-Phe-Arg-AMC) and rSl-CathL exhibited hydrolysis preference for Z-Leu-Arg-AMC (k(cat)/K-m = 37.53 mM S-1), which is similar to other insect cathepsin L-like enzymes. rSl-CathL activity inhibition assays were performed using four recombinant sugarcane cystatins. rSl-CathL was strongly inhibited by recombinant cystatin CaneCPI-4 (K-i = 0.196 nM), indicating that this protease is a potential target for pest control. (C) 2011 Elsevier Ltd. All rights reserved.

Design of potent and selective human cathepsin K inhibitors that span the active site

Potent and selective active-site-spanning inhibitors have been designed for cathepsin K, a cysteine protease unique to osteoclasts. They act by mechanisms that involve tight binding intermediates, potentially on a hydrolytic pathway. X-ray crystallographic, MS, NMR spectroscopic, and kinetic studies of the mechanisms of inhibition indicate that different intermediates or transition states are being represented that are dependent on the conditions of measurement and the specific groups flanking the carbonyl in the inhibitor. The species observed crystallographically are most consistent with tetrahedral intermediates that may be close approximations of those that occur during substrate hydrolysis. Initial kinetic studies suggest the possibility of irreversible and reversible active-site modification. Representative inhibitors have demonstrated antiresorptive activity both in vitro and in vivo and therefore are promising leads for therapeutic agents for the treatment of osteoporosis. Expansion of these inhibitor concepts can be envisioned for the many other cysteine proteases implicated for therapeutic intervention.

Peptide inhibitors of HIV-1 protease and viral infection of peripheral blood lymphocytes based on HIV-1 Vif

We recently reported that HIV-1 Vif (virion infectivity factor) inhibits HIV-1 protease in vitro and in bacteria, suggesting that it may serve as the basis for the design of new protease inhibitors and treatment for HIV-1 infection. To evaluate this possibility, we synthesized peptide derivatives from the region of Vif, which inhibits protease, and tested their activity on protease. In an assay of cleavage of virion-like particles composed of HIV-1 Gag precursor polyprotein, full-length recombinant Vif, and a peptide consisting of residues 21–65 of Vif, but not a control peptide or BSA, inhibited protease activity. Vif21–65 blocked protease at a molar ratio of two to one. We then tested this peptide and a smaller peptide, Vif41–65, for their effects on HIV-1 infection of peripheral blood lymphocytes. Both Vif peptides inhibited virus expression below the limit of detection, but control peptides had no effect. To investigate its site of action, Vif21–65 was tested for its effect on Gag cleavage by protease during HIV-1 infection. We found that commensurate with its reduction of virus expression, Vif21–65 inhibited the cleavage of the polyprotein p55 to mature p24. These results are similar to those obtained by using Ro 31–8959, a protease inhibitor in clinical use. We conclude that Vif-derived peptides inhibit protease during HIV-1 infection and may be useful for the development of new protease inhibitors.

Molecular ordering of the Fas-apoptotic pathway: The Fas/APO-1 protease Mch5 is a CrmA-inhibitable protease that activates multiple Ced-3/ICE-like cysteine proteases

The Fas/APO-1-receptor associated cysteine protease Mch5 (MACH/FLICE) is believed to be the enzyme responsible for activating a protease cascade after Fas-receptor ligation, leading to cell death. The Fas-apoptotic pathway is potently inhibited by the cowpox serpin CrmA, suggesting that Mch5 could be the target of this serpin. Bacterial expression of proMch5 generated a mature enzyme composed of two subunits, which are derived from the precursor proenzyme by processing at Asp-227, Asp-233, Asp-391, and Asp-401. We demonstrate that recombinant Mch5 is able to process/activate all known ICE/Ced-3-like cysteine proteases and is potently inhibited by CrmA. This contrasts with the observation that Mch4, the second FADD-related cysteine protease that is also able to process/activate all known ICE/Ced-3-like cysteine proteases, is poorly inhibited by CrmA. These data suggest that Mch5 is the most upstream protease that receives the activation signal from the Fas-receptor to initiate the apoptotic protease cascade that leads to activation of ICE-like proteases (TX, ICE, and ICE-relIII), Ced-3-like proteases (CPP32, Mch2, Mch3, Mch4, and Mch6), and the ICH-1 protease. On the other hand, Mch4 could be a second upstream protease that is responsible for activation of the same protease cascade in CrmA-insensitive apoptotic pathways.

Unexpected Activity of a Novel Kunitz-Type Inhibtior: Inhibition of Cysteine Proteases but not Serine Proteases

Kunitz-type (KT) protease inhibitors are low molecular weight proteins classically defined as serine protease inhibitors. We identified a novel secreted KT inhibitor associated with the gut and parenchymal tissues of the infective juvenile stage of Fasciola hepatica, a helminth parasite of medical and veterinary importance. Unexpectedly, recombinant KT inhibitor (rFhKT1) exhibited no inhibitory activity towards serine proteases but was a potent inhibitor of the major secreted cathepsin L cysteine proteases of F. hepatica, FhCL1 and FhCL2, and of human cathepsins L and K (Ki = 0.4 nM - 27 nM). FhKT1 prevented the auto-catalytic activation of FhCL1 and FhCL2 and formed stable complexes with the mature enzymes. Pull-down experiments from adult parasite culture medium showed that rFhKT1 interacts specifically with native secreted FhCL1, FhCL2 and FhCL5. Substitution of the unusual P1 Leu15 within the exposed reactive loop of FhKT1 for the more commonly found Arg (FhKT1Leu15/Arg15) had modest adverse effects on the cysteine protease inhibition but conferred potent activity against the serine protease trypsin (Ki = 1.5 nM). Computational docking and sequence analysis provided hypotheses for the exclusive binding of FhKT1 to cysteine proteases, the importance of the Leu15 in anchoring the inhibitor into the S2 active site pocket, and the inhibitor's selectivity towards FhCL1, FhCL2 and human cathepsins L and K. FhKT1 represents a novel evolutionary adaptation of KT protease inhibitors by F. hepatica, with its prime purpose likely in the regulation of the major parasite-secreted proteases and/or cathepsin L-like proteases of its host.

Analysis of human immunodeficiency virus type 1 drug resistance in children receiving nucleoside analogue reverse-transcriptase inhibitors plus nevirapine, nelfinavir, or ritonavir (Pediatric AIDS Clinical Trials Group 377)

In Pediatric AIDS Clinical Trials Group 377, antiretroviral therapy-experienced children were randomized to 4 treatment arms that included different combinations of stavudine, lamivudine (3TC), nevirapine (Nvp), nelfinavir (Nfv), and ritonavir (Rtv). Previous treatment with zidovudine (Zdv), didanosine (ddI), or zalcitabine (ddC) was acceptable. Drug resistance ((R)) mutations were assessed before study treatment (baseline) and at virologic failure. Zdv(R), ddI(R), and ddC(R) mutations were detected frequently at baseline but were not associated with virologic failure. Children with drug resistance mutations at baseline had greater reductions in virus load over time than did children who did not. Nvp(R) and 3TC(R) mutations were detected frequently at virologic failure, and Nvp(R) mutations were more common among children receiving 3-drug versus 4-drug Nvp-containing regimens. Children who were maintained on their study regimen after virologic failure accumulated additional Nvp(R) and 3TC(R) mutations plus Rtv(R) and Nfv(R) mutations. However, Rtv(R) and Nfv(R) mutations were detected at unexpectedly low rates.

Structural Basis of Drug Resistance by Genetic Variants of HIV Type 1 Clade C Protease from India

Using computer modeling of three-dimensional structures and structural information available on the crystal structures of HIV-1 protease, we investigated the structural effects of mutations, in treatment-naive and treatment-exposed individuals from India and postulated mechanisms of resistance in clade C variants. A large number of models (14) have been generated by computational mutation of the available crystal structures of drug bound proteases. Localized energy minimization was carried out in and around the sites of mutation in order to optimize the geometry of interactions present. Most of the mutations result in structural differences at the flap that favors the semiopen state of the enzyme. Some of the mutations were also found to confer resistance by affecting the geometry of the active site. The E35D mutation affects the flap structure in clade B strains and E35N and E35K mutation, seen in our modeled strains, have a more profound effect. Common polymorphisms at positions 36 and 63 in clade C also affected flap structure. Apart from a few other residues Gln-58, Asn-83, Asn-88, and Gln-92 and their interactions are important for the transition from the closed to the open state. Development of protease inhibitors by structure-based design requires investigation of mechanisms operative for clade C to improve the efficacy of therapy.

Conserved water-mediated H-bonding dynamics of catalytic Asn 175 in plant thiol protease

The role of invariant water molecules in the activity of plant cysteine protease is ubiquitous in nature. On analysing the 11 different Protein DataBank (PDB) structures of plant thiol proteases, the two invariant water molecules W I and W2 (W220 and W222 in the template 1PPN structure) were observed to form H-bonds with the Ob atom of Asn 175. Extensive energy minimization and molecular dynamics simulation studies up to 2 ns on all the PDB and solvated structures clearly revealed the involvement of the H-bonding association of the two water molecules in fixing the orientation of the asparagine residue of the catalytic triad. From this study, it is suggested that H-bonding of the water molecule at the W1 invariant site better stabilizes the Asn residue at the active site of the catalytic triad.