Structure-based discovery of novel cruzain inhibitors with distinct trypanocidal activity profiles.

Over 110 years after the first formal description of Chagas disease, the trypanocidal drugs thus far available have limited efficacy and several side effects. This encourages the search for novel treatments that inhibit T. cruzi targets. One of the most studied anti-T. cruzi targets is the cysteine protease cruzain; it is associated with metacyclogenesis, replication, and invasion of the host cells. We used computational techniques to identify novel molecular scaffolds that act as cruzain inhibitors. First, with a docking-based virtual screening, we identified compound 8, a competitive cruzain inhibitor with a Ki of 4.6 µM. Then, aided by molecular dynamics simulations, cheminformatics, and docking, we identified the analog compound 22 with a Ki of 27 µM. Surprisingly, despite sharing the same isoquinoline scaffold, compound 8 presented higher trypanocidal activity against the epimastigote forms, while compound 22, against the trypomastigotes and amastigotes. Taken together, compounds 8 and 22 represent a promising scaffold for further development of trypanocidal compounds as drug candidates for treating Chagas disease.

Staphylococcus aureus Exfoliative Toxin E, Oligomeric State and Flip of P186: Implications for Its Action Mechanism.

Staphylococcal exfoliative toxins (ETs) are glutamyl endopeptidases that specifically cleave the Glu381-Gly382 bond in the ectodomains of desmoglein 1 (Dsg1) via complex action mechanisms. To date, four ETs have been identified in different Staphylococcus aureus strains and ETE is the most recently characterized. The unusual properties of ETs have been attributed to a unique structural feature, i.e., the 180° flip of the carbonyl oxygen (O) of the nonconserved residue 192/186 (ETA/ETE numbering), not conducive to the oxyanion hole formation. We report the crystal structure of ETE determined at 1.61 Å resolution, in which P186(O) adopts two conformations displaying a 180° rotation. This finding, together with free energy calculations, supports the existence of a dynamic transition between the conformations under the tested conditions. Moreover, enzymatic assays showed no significant differences in the esterolytic efficiency of ETE and ETE/P186G, a mutant predicted to possess a functional oxyanion hole, thus downplaying the influence of the flip on the activity. Finally, we observed the formation of ETE homodimers in solution and the predicted homodimeric structure revealed the participation of a characteristic nonconserved loop in the interface and the partial occlusion of the protein active site, suggesting that monomerization is required for enzymatic activity.

In silico identification of noncompetitive inhibitors targeting an uncharacterized allosteric site of falcipain-2.

Falcipain-2 (FP-2) is a Plasmodium falciparum hemoglobinase widely targeted in the search for antimalarials. FP-2 can be allosterically modulated by various noncompetitive inhibitors that have been serendipitously identified. Moreover, the crystal structures of two inhibitors bound to an allosteric site, termed site 6, of the homolog enzyme human cathepsin K (hCatK) suggest that the equivalent region in FP-2 might play a similar role. Here, we conduct the rational identification of FP-2 inhibitors through virtual screenings (VS) of compounds into several pocket-like conformations of site 6, sampled during molecular dynamics (MD) simulations of the free enzyme. Two noncompetitive inhibitors, ZINC03225317 and ZINC72290660, were confirmed using in vitro enzymatic assays and their poses into site 6 led to calculated binding free energies matching the experimental ones. Our results provide strong evidence about the allosteric inhibition of FP-2 through binding of small molecules to site 6, thus opening the way toward the discovery of new inhibitors against this enzyme.

Computational study on the allosteric mechanism of Leishmania major IF4E-1 by 4E-interacting protein-1: Unravelling the determinants of m7GTP cap recognition.

During their life cycle, Leishmania parasites display a fine-tuned regulation of the mRNA translation through the differential expression of isoforms of eukaryotic translation initiation factor 4E (LeishIF4Es). The interaction between allosteric modulators such as 4E-interacting proteins (4E-IPs) and LeishIF4E affects the affinity of this initiation factor for the mRNA cap. Here, several computational approaches were employed to elucidate the molecular bases of the previously-reported allosteric modulation in L. major exerted by 4E-IP1 (Lm4E-IP1) on eukaryotic translation initiation factor 4E 1 (LmIF4E-1). Molecular dynamics (MD) simulations and accurate binding free energy calculations (ΔGbind ) were combined with network-based modeling of residue-residue correlations. We also describe the differences in internal motions of LmIF4E-1 apo form, cap-bound, and Lm4E-IP1-bound systems. Through community network calculations, the differences in the allosteric pathways of allosterically-inhibited and active forms of LmIF4E-1 were revealed. The ΔGbind values show significant differences between the active and inhibited systems, which are in agreement with the available experimental data. Our study thoroughly describes the dynamical perturbations of LmIF4E-1 cap-binding site triggered by Lm4E-IP1. These findings are not only essential for the understanding of a critical process of trypanosomatids' gene expression but also for gaining insight into the allostery of eukaryotic IF4Es, which could be useful for structure-based design of drugs against this protein family.

Water Bridges Play a Key Role in Affinity and Selectivity for Malarial Protease Falcipain-2.

Falcipain-2 (FP-2) is hemoglobinase considered an attractive drug target of Plasmodium falciparum. Recently, it has been shown that peptidomimetic nitriles containing a 3-pyridyl (3Pyr) moiety at P2 display high affinity and selectivity for FP-2 with respect to human cysteine cathepsins (hCats), outperforming other P2-Pyr isomers and analogs. Further characterization demonstrated that certain P3 variants of these compounds possess micromolar inhibition of parasite growth in vitro and no cytotoxicity against human cell lines. However, the structural determinants underlying the selectivity of the 3Pyr-containing nitriles for FP-2 remain unknown. In this work, we conduct a thorough computational study combining MD simulations and free energy calculations to decipher the bases of the selectivity of the aforementioned nitriles. Our results reveal that water bridges involving the nitrogen and one carboxyl oxygen of I85 and D234 of FP-2, respectively, and the nitrogen of the neutral 3Pyr moiety, which are either less prevalent or nonexistent in the other complexes, explain the experimental activity profiles. The presence of crystallographic waters close to the bridging water positions in the studied proteases strongly supports the occurrence of such interactions. Overall, our findings suggest that selective FP-2 inhibitors can be designed by promoting water bridge formation at the bottom of the S2 subsite and/or by introducing complementary groups that displace the bridging water.

Physicochemical and structural properties of lunasin revealed by spectroscopic, chromatographic and molecular dynamics approaches.

Lunasin is a 43-amino acid peptide from seeds and grains with bioavailability in humans and potent chemotherapeutic action against several cancer cell lines. Here, we investigate new information about the physicochemical and structural properties of lunasin using circular dichroism (CD), fluorescence spectroscopy, electrospray ionization-ion mobility spectrometry-mass spectrometry (ESI-IMS-MS), size exclusion chromatography (SEC), molecular dynamics (MD), and bioinformatics. CD analysis and disorder prediction obtained by PONDR indicate that lunasin has a mostly unordered structure. Double wavelength [θ]222nm x [θ]200nm plot data suggests that lunasin is an intrinsically disordered peptide (IDP) in a pre-molten globule-like (PMG-like) state, while CD spectrum deconvolution and MD simulation indicate small ß-strand content. The presence of residual structure was supported by loss of CD signal at 222 nm after treatment with urea and by increasing fluorescence emission upon bis-ANS binding. Lunasin also demonstrated stability to heating up to the temperature of 100 °C, as verified by CD. MD and CD analyses in the presence of TFE and MoRFpred prediction indicated the helix propensity of lunasin. ESI-IMS-MS data revealed that lunasin shows a propensity to form disulfide bonds at the conditions used. MD data also indicated that disulfide bond formation affects the adopted structure, showing a possible role of aspartyl-end in structure stabilization and compaction. In conclusion, our data support a characterization of lunasin as a peptide with an intrinsic disorder in a PMG-like state and reveal new aspects about its structural stability and plasticity, as well as the effects of disulfide bond formation and electrostatic attractions.

Pisum sativum Defensin 1 Eradicates Mouse Metastatic Lung Nodules from B16F10 Melanoma Cells.

Psd1 is a pea plant defensin which can be actively expressed in Pichia pastoris and shows broad antifungal activity. This activity is dependent on fungal membrane glucosylceramide (GlcCer), which is also important for its internalization, nuclear localization, and endoreduplication. Certain cancer cells present a lipid metabolism imbalance resulting in the overexpression of GlcCer in their membrane. In this work, in vitroassays using B16F10 cells showed that labeled fluorescein isothiocyanate FITC-Psd1 internalized into live cultured cells and targeted the nucleus, which underwent fragmentation, exhibiting approximately 60% of cells in the sub-G0/G1 stage. This phenomenon was dependent on GlcCer, and the participation of cyclin-F was suggested. In a murine lung metastatic melanoma model, intravenous injection of Psd1 together with B16F10 cells drastically reduced the number of nodules at concentrations above 0.5 mg/kg. Additionally, the administration of 1 mg/kg Psd1 decreased the number of lung inflammatory cells to near zero without weight loss, unlike animals that received melanoma cells only. It is worth noting that 1 mg/kg Psd1 alone did not provoke inflammation in lung tissue or weight or vital signal losses over 21 days, inferring no whole animal cytotoxicity. These results suggest that Psd1 could be a promising prototype for human lung anti-metastatic melanoma therapy.

Prediction of Noncompetitive Inhibitor Binding Mode Reveals Promising Site for Allosteric Modulation of Falcipain-2.

Falcipain-2 (FP-2) is a Plasmodium falciparum cysteine protease that has been extensively targeted to identify novel antimalarials. Remarkably, previous reports have shown that FP-2 can be allosterically modulated and, for a particular noncompetitive chalcone inhibitor, the existing lines of experimental evidence can guide the prediction of its unknown binding mode to the enzyme in a reliable fashion. In this work, we propose a structure of FP-2 in complex with the aforementioned compound that fulfills all of the experimental data, by employing a combination of molecular modeling tools, such as pocket volume measurements, docking, molecular dynamics (MD) simulations, and free energy calculations. Our results show that the studied inhibitor binds a transient pocket occluded in all of the available FP-2 crystal structures and lying in a region previously characterized as a potential allosteric site in related cysteine proteases. In addition, we detected in silico the occurrence of significant community reorganization in FP-2, increased signal transmission between the allosteric pocket and the active site, and change in loop motions and residue pKa values upon the compound binding, thus providing insight into the uncharacterized allosteric mechanism. Overall, this study yields valuable predictions for the design of novel allosteric inhibitors against FP-2 and other cysteine proteases.

Dissecting a novel allosteric mechanism of cruzain: A computer-aided approach.

Trypanosoma cruzi is the causative agent of Chagas disease, a neglected infection affecting millions of people in tropical regions. There are several chemotherapeutic agents for the treatment of this disease, but most of them are highly toxic and generate resistance. Currently, the development of allosteric inhibitors constitutes a promising research field, since it can improve the accessibility to more selective and less toxic medicines. To date, the allosteric drugs prediction is a state-of-the-art topic in rational structure-based computational design. In this work, a simulation strategy was developed for computational discovery of allosteric inhibitors, and it was applied to cruzain, a promising target and the major cysteine protease of T. cruzi. Molecular dynamics simulations, binding free energy calculations and network-based modelling of residue interactions were combined to characterize and compare molecular distinctive features of the apo form and the cruzain-allosteric inhibitor complexes. By using geometry-based criteria on trajectory snapshots, we predicted two main allosteric sites suitable for drug targeting. The results suggest dissimilar mechanisms exerted by the same allosteric site when binding different potential allosteric inhibitors. Finally, we identified the residues involved in suboptimal paths linking the identified site and the orthosteric site. The present study constitutes the first approximation to the design of cruzain allosteric inhibitors and may serve for future pharmacological intervention. Here, no major effects on active site structure were observed due to compound binding (modification of distance and angles between catalytic residues), which indicates that allosteric regulation in cruzain might be mediated via alterations of its dynamical properties similarly to allosteric inhibition of human cathepsin K (HCatK). The current findings are particularly relevant for the design of allosteric modulators of papain-like cysteine proteases.

Investigation of the binding mode of a novel cruzain inhibitor by docking, molecular dynamics, ab initio and MM/PBSA calculations.

Chagas disease remains a major health problem in South America, and throughout the world. The two drugs clinically available for its treatment have limited efficacy and cause serious adverse effects. Cruzain is an established therapeutic target of Trypanosoma cruzi, the protozoan that causes Chagas disease. Our group recently identified a competitive cruzain inhibitor (compound 1) with an IC50 = 15 µM that is also more synthetically accessible than the previously reported lead, compound 2. Prior studies, however, did not propose a binding mode for compound 1, hindering understanding of the structure-activity relationship and optimization. Here, the cruzain binding mode of compound 1 was investigated using docking, molecular dynamics (MD) simulations with ab initio derived parameters, ab initio calculations, and MM/PBSA. Two ligand protonation states and four binding poses were evaluated. A careful ligand parameterization method was employed to derive more physically meaningful parameters than those obtained by automated tools. The poses of unprotonated 1 were unstable in MD, showing large conformational changes and diffusing away from the binding site, whereas the protonated form showed higher stability and interaction with negatively charged residues Asp161 and Cys25. MM/PBSA also suggested that these two residues contribute favorably to binding of compound 1. By combining results from MD, ab initio calculations, and MM/PBSA, a binding mode of 1 is proposed. The results also provide insights for further optimization of 1, an interesting lead compound for the development of new cruzain inhibitors.

Predicting binding modes of reversible peptide-based inhibitors of falcipain-2 consistent with structure-activity relationships.

Falcipain-2 (FP-2) is a major hemoglobinase of Plasmodium falciparum, considered an important drug target for the development of antimalarials. A previous study reported a novel series of 20 reversible peptide-based inhibitors of FP-2. However, the lack of tridimensional structures of the complexes hinders further optimization strategies to enhance the inhibitory activity of the compounds. Here we report the prediction of the binding modes of the aforementioned inhibitors to FP-2. A computational approach combining previous knowledge on the determinants of binding to the enzyme, docking, and postdocking refinement steps, is employed. The latter steps comprise molecular dynamics simulations and free energy calculations. Remarkably, this approach leads to the identification of near-native ligand conformations when applied to a validation set of protein-ligand structures. Overall, we proposed substrate-like binding modes of the studied compounds fulfilling the structural requirements for FP-2 binding and yielding free energy values that correlated well with the experimental data. Proteins 2017; 85:1666-1683. © 2017 Wiley Periodicals, Inc.

Insights into cytochrome bc1 complex binding mode of antimalarial 2-hydroxy-1, 4-naphthoquinones through molecular modelling

BACKGROUND Malaria persists as a major public health problem. Atovaquone is a drug that inhibits the respiratory chain of Plasmodium falciparum, but with serious limitations like known resistance, low bioavailability and high plasma protein binding. OBJECTIVES The aim of this work was to perform molecular modelling studies of 2-hydroxy-1,4-naphthoquinones analogues of atovaquone on the Qo site of P. falciparum cytochrome bc1 complex (Pfbc1) to suggest structural modifications that could improve their antimalarial activity. METHODS We have built the homology model of the cytochrome b (CYB) and Rieske iron-sulfur protein (ISP) subunits from Pfbc1 and performed the molecular docking of 41 2-hydroxy-1,4-naphthoquinones with known in vitro antimalarial activity and predicted to act on this target. FINDINGS Results suggest that large hydrophobic R2 substituents may be important for filling the deep hydrophobic Qo site pocket. Moreover, our analysis indicates that the H-donor 2-hydroxyl group may not be crucial for efficient binding and inhibition of Pfbc1 by these atovaquone analogues. The C1 carbonyl group (H-acceptor) is more frequently involved in the important hydrogen bonding interaction with His152 of the Rieske ISP subunit. MAIN CONCLUSIONS Additional interactions involving residues such as Ile258 and residues required for efficient catalysis (e.g., Glu261) could be explored in drug design to avoid development of drug resistance by the parasite.

Insights into cytochrome bc1 complex binding mode of antimalarial 2-hydroxy-1,4-naphthoquinones through molecular modelling.

BACKGROUND: Malaria persists as a major public health problem. Atovaquone is a drug that inhibits the respiratory chain of Plasmodium falciparum, but with serious limitations like known resistance, low bioavailability and high plasma protein binding. OBJECTIVES: The aim of this work was to perform molecular modelling studies of 2-hydroxy-1,4-naphthoquinones analogues of atovaquone on the Qo site of P. falciparum cytochrome bc1 complex (Pfbc1) to suggest structural modifications that could improve their antimalarial activity. METHODS: We have built the homology model of the cytochrome b (CYB) and Rieske iron-sulfur protein (ISP) subunits from Pfbc1 and performed the molecular docking of 41 2-hydroxy-1,4-naphthoquinones with known in vitro antimalarial activity and predicted to act on this target. FINDINGS: Results suggest that large hydrophobic R2 substituents may be important for filling the deep hydrophobic Qo site pocket. Moreover, our analysis indicates that the H-donor 2-hydroxyl group may not be crucial for efficient binding and inhibition of Pfbc1 by these atovaquone analogues. The C1 carbonyl group (H-acceptor) is more frequently involved in the important hydrogen bonding interaction with His152 of the Rieske ISP subunit. MAIN CONCLUSIONS: Additional interactions involving residues such as Ile258 and residues required for efficient catalysis (e.g., Glu261) could be explored in drug design to avoid development of drug resistance by the parasite.

Insights into the Interactions of Fasciola hepatica Cathepsin L3 with a Substrate and Potential Novel Inhibitors through In Silico Approaches.

BACKGROUND: Fasciola hepatica is the causative agent of fascioliasis, a disease affecting grazing animals, causing economic losses in global agriculture and currently being an important human zoonosis. Overuse of chemotherapeutics against fascioliasis has increased the populations of drug resistant parasites. F. hepatica cathepsin L3 is a protease that plays important roles during the life cycle of fluke. Due to its particular collagenolytic activity it is considered an attractive target against the infective phase of F. hepatica. METHODOLOGY/PRINCIPAL FINDINGS: Starting with a three dimensional model of FhCL3 we performed a structure-based design of novel inhibitors through a computational study that combined virtual screening, molecular dynamics simulations, and binding free energy (ΔGbind) calculations. Virtual screening was carried out by docking inhibitors obtained from the MYBRIDGE-HitFinder database inside FhCL3 and human cathepsin L substrate-binding sites. On the basis of dock-scores, five compounds were predicted as selective inhibitors of FhCL3. Molecular dynamic simulations were performed and, subsequently, an end-point method was employed to predict ΔGbind values. Two compounds with the best ΔGbind values (-10.68 kcal/mol and -7.16 kcal/mol), comparable to that of the positive control (-10.55 kcal/mol), were identified. A similar approach was followed to structurally and energetically characterize the interface of FhCL3 in complex with a peptidic substrate. Finally, through pair-wise and per-residue free energy decomposition we identified residues that are critical for the substrate/ligand binding and for the enzyme specificity. CONCLUSIONS/SIGNIFICANCE: The present study is the first computer-aided drug design approach against F. hepatica cathepsins. Here we predict the principal determinants of binding of FhCL3 in complex with a natural substrate by detailed energetic characterization of protease interaction surface. We also propose novel compounds as FhCL3 inhibitors. Overall, these results will foster the future rational design of new inhibitors against FhCL3, as well as other F. hepatica cathepsins.

Differential effects of CSF-1R D802V and KIT D816V homologous mutations on receptor tertiary structure and allosteric communication.

The colony stimulating factor-1 receptor (CSF-1R) and the stem cell factor receptor KIT, type III receptor tyrosine kinases (RTKs), are important mediators of signal transduction. The normal functions of these receptors can be compromised by gain-of-function mutations associated with different physiopatological impacts. Whereas KIT D816V/H mutation is a well-characterized oncogenic event and principal cause of systemic mastocytosis, the homologous CSF-1R D802V has not been identified in human cancers. The KIT D816V oncogenic mutation triggers resistance to the RTK inhibitor Imatinib used as first line treatment against chronic myeloid leukemia and gastrointestinal tumors. CSF-1R is also sensitive to Imatinib and this sensitivity is altered by mutation D802V. Previous in silico characterization of the D816V mutation in KIT evidenced that the mutation caused a structure reorganization of the juxtamembrane region (JMR) and facilitated its departure from the kinase domain (KD). In this study, we showed that the equivalent CSF-1R D802V mutation does not promote such structural effects on the JMR despite of a reduction on some key H-bonds interactions controlling the JMR binding to the KD. In addition, this mutation disrupts the allosteric communication between two essential regulatory fragments of the receptors, the JMR and the A-loop. Nevertheless, the mutation-induced shift towards an active conformation observed in KIT D816V is not observed in CSF-1R D802V. The distinct impact of equivalent mutation in two homologous RTKs could be associated with the sequence difference between both receptors in the native form, particularly in the JMR region. A local mutation-induced perturbation on the A-loop structure observed in both receptors indicates the stabilization of an inactive non-inhibited form, which Imatinib cannot bind.

Structural analysis of Pla protein from the biological warfare agent Yersinia pestis: docking and molecular dynamics of interactions with the mammalian plasminogen system.

Yersinia pestis protein Pla is a plasmid-coded outer membrane protein with aspartic-protease activity. Pla exhibits a plasminogen (Plg) activator activity (PAA) that promotes the cleavage of Plg to the active serine-protease form called plasmin. Exactly how Pla activates Plg into plasmin remains unclear. To investigate this event, we performed the interactions between the predicted Plg and Pla protein structures by rigid-body docking with the HEX program and evaluated the complex stability by molecular dynamics (MD) using the GROMACS package programs. The predicted docked complex of Plg-Pla shows the same interaction site predicted by experimental site-direct mutagenesis in other studies. After a total of 8 ns of MD simulation, we observed the relaxation of the beta-barrel structure of Pla and the progressive approximation and stabilization between the cleavage site of Plg into the extracellular loops of Pla, followed by the increase in the number of H bonds. We also report here the aminoacids that participate in the active site and the sub sites of interaction. The total understanding of these interactions can be an important tool for drug design against bacterial proteases.

Molecular dynamics simulations of peptide inhibitors complexed with Trypanosoma cruzi trypanothione reductase.

The drugs against tropical neglected diseases, especially Chagas' Disease, were launched more than 30 years ago, and the development of resistance requires the discovery of new and more effective chemotherapeutic agents. Trypanosoma cruzi has a redox enzyme called trypanothione reductase which was successfully inhibited for peptide derivatives (McKie et al., Amino Acids, 2001, 20: 145). This work aims at studying the mechanism of inhibition of this enzyme through molecular dynamics simulations and evaluating the behavior of some derivatives when inhibiting this protein. We should affirm that any particular molecular dynamics analysis tools (Hbond pattern, 3-D root-mean-square deviation, solvent accessible surface area, etc.) cannot be used apart from the others to justify completely these peptides inhibitory patterns. Based on our results, we reproduced the experimental data and, moreover, we discriminated against a new site in enzyme aperture, which can assist the development of powerful inhibitors against trypanothione reductase enzyme.

Structural analysis of the N-terminal fragment of the antiangiogenic protein endostatin: a molecular dynamics study.

Endostatin is a potent antiangiogenic protein derived from the noncollagenous domain 1 (NC1) of collagen XVIII. The mechanism by which endostatin exerts its antiangiogenic effect is still incompletely understood. It has been shown that the 27 amino acid N-terminal fragment of murine endostatin has antitumor, antimigration, and antipermeability activities comparable to the full soluble protein. To understand how this peptide can exert such elaborate function, we performed structural analysis using molecular dynamics to evaluate the behavior of this fragment in aqueous environment. Here, we show that the N-terminal peptide of murine endostatin is able to assume a well-defined structure, folding into a zinc-dependent ß-hairpin conformation. Analyzing the folding mechanism, we were able to understand why the N-terminal peptide of human endostatin with the same length failed to acquire a stable conformation. Conversely, we were able to predict the successful folding of the R4Q mutant and of a shorter form of the human peptide with 25 residues. Finally, we show that the ß-hairpin conformation assumed by the zinc-bound peptide of murine endostatin has a high structural similarity with fragments of another family of angiogenesis inhibitors: the integrin-binding portion of the NC1 domain of collagen IV. Indeed, our docking simulations show that arresten, canstatin, and the endostatin peptide bind to the same spot of αVß3 integrin, suggesting similar interactions via a common binding site on this receptor.

High temperatures enhance cooperative motions between CBM and catalytic domains of a thermostable cellulase: mechanism insights from essential dynamics.

Cellulases from thermophiles are capable of cleaving sugar chains from cellulose efficiently at high temperatures. The thermo-resistant Cel9A-68 cellulase possesses two important domains: CBM and a catalytic domain connected by a Pro/Ser/Thr rich linker. These domains act cooperatively to allow efficient catalysis. Despite exhaustive efforts to characterize cellulase binding and mechanism of action, a detailed description of the cellulose intrinsic flexibility is still lacking. From computational simulations we studied the temperature influence on the enzyme plasticity, prior to substrate binding. Interestingly, we observed an enhancement of collective motions at high temperatures. These motions are the most representative and describe an intrinsic hinge bending transition. A detailed analysis of these motions revealed an interdomain approximation where D459 and G460, located at the linker region, are the hinge residues. Therefore, we propose a new putative site for mutagenesis targeting the modulation of such conformational transition that may be crucial for activity.

Conformational selection, dynamic restriction and the hydrophobic effect coupled to stabilization of the BIR3 domain of the human X-linked inhibitor of apoptosis protein by the tetrapeptide AVPI.

The XIAP-BIR3 domain blocks a substantial portion of the apoptosis pathway and is an attractive target for novel anticancer agents. The tetrapeptide AVPI, from the protein Smac/DIABLO, binds to the XIAP-BIR3 domain, allowing the cancer cells to die. Here we characterize the binding parameters of AVPI to XIAP-BIR3 and analyze its effects on the thermodynamic stability of this domain. XIAP-BIR3 was exceptionally stable against physical and chemical treatments and became even more stable by interaction with AVPI. Nuclear magnetic resonance experiments demonstrated that conformational selection is taking place upon AVPI interaction with XIAP-BIR3. Molecular dynamics simulations corroborate that the flexibility of XIAP-BIR3 is significantly reduced. The positive binding entropy associated with a loss of conformational entropy involved in the binding indicates that hydrophobic interactions play an important role in the interaction and domain stabilization. The mechanism of XIAP-BIR3 stabilization and its implications for drug affinity optimization are discussed.