Solvated interaction energy: from small-molecule to antibody drug design.

Scoring functions are ubiquitous in structure-based drug design as an aid to predicting binding modes and estimating binding affinities. Ideally, a scoring function should be broadly applicable, obviating the need to recalibrate and refit its parameters for every new target and class of ligands. Traditionally, drugs have been small molecules, but in recent years biologics, particularly antibodies, have become an increasingly important if not dominant class of therapeutics. This makes the goal of having a transferable scoring function, i.e., one that spans the range of small-molecule to protein ligands, even more challenging. One such broadly applicable scoring function is the Solvated Interaction Energy (SIE), which has been developed and applied in our lab for the last 15 years, leading to several important applications. This physics-based method arose from efforts to understand the physics governing binding events, with particular care given to the role played by solvation. SIE has been used by us and many independent labs worldwide for virtual screening and discovery of novel small-molecule binders or optimization of known drugs. Moreover, without any retraining, it is found to be transferrable to predictions of antibody-antigen relative binding affinities and as accurate as functions trained on protein-protein binding affinities. SIE has been incorporated in conjunction with other scoring functions into ADAPT (Assisted Design of Antibody and Protein Therapeutics), our platform for affinity modulation of antibodies. Application of ADAPT resulted in the optimization of several antibodies with 10-to-100-fold improvements in binding affinity. Further applications included broadening the specificity of a single-domain antibody to be cross-reactive with virus variants of both SARS-CoV-1 and SARS-CoV-2, and the design of safer antibodies by engineering of a pH switch to make them more selective towards acidic tumors while sparing normal tissues at physiological pH.

Programmable Attenuation of Antigenic Sensitivity for a Nanobody-Based EGFR Chimeric Antigen Receptor Through Hinge Domain Truncation.

Epidermal growth factor family receptor (EGFR) is commonly overexpressed in many solid tumors and an attractive target for chimeric antigen receptor (CAR)-T therapy, but as EGFR is also expressed at lower levels in healthy tissues a therapeutic strategy must balance antigenic responsiveness against the risk of on-target off-tumor toxicity. Herein, we identify several camelid single-domain antibodies (also known as nanobodies) that are effective EGFR targeting moieties for CARs (EGFR-sdCARs) with very strong reactivity to EGFR-high and EGFR-low target cells. As a strategy to attenuate their potent antigenic sensitivity, we performed progressive truncation of the human CD8 hinge commonly used as a spacer domain in many CAR constructs. Single amino acid hinge-domain truncation progressively decreased both EGFR-sdCAR-Jurkat cell binding to EGFR-expressing targets and expression of the CD69 activation marker. Attenuated signaling in hinge-truncated EGFR-sdCAR constructs increased selectivity for antigen-dense EGFR-overexpressing cells over an EGFR-low tumor cell line or healthy donor derived EGFR-positive fibroblasts. We also provide evidence that epitope location is critical for determining hinge-domain requirement for CARs, as hinge truncation similarly decreased antigenic sensitivity of a membrane-proximal epitope targeting HER2-CAR but not a membrane-distal EGFRvIII-specific CAR. Hinge-modified EGFR-sdCAR cells showed clear functional attenuation in Jurkat-CAR-T cells and primary human CAR-T cells from multiple donors in vitro and in vivo. Overall, these results indicate that hinge length tuning provides a programmable strategy for throttling antigenic sensitivity in CARs targeting membrane-proximal epitopes, and could be employed for CAR-optimization and improved tumor selectivity.

Antibody mutations favoring pH-dependent binding in solid tumor microenvironments: Insights from large-scale structure-based calculations.

Antibody-based therapeutics for treatment of various tumors have grown rapidly in recent years. Unfortunately, safety issues, attributed to off-tumor effects and cytotoxicity, are still a significant concern with the standard of care. Improvements to ensure targeted delivery of antitumor pharmaceuticals are desperately needed. We previously demonstrated that incorporating histidyl pH-switches in an anti-HER2 antibody induced selective antigen binding under acidic pH conditions (MAbs 2020;12:1682866). This led to an improved safety profile due to preferential targeting of the oncoprotein in the acidic solid tumor microenvironment. Following this success, we expanded this approach to a set of over 400 antibody structures complexed with over 100 different human oncoproteins, associated with solid tumors. Calculations suggested that mutations to His of certain residue types, namely Trp, Arg, and Tyr, could be significantly more successful for inducing pH-dependent binding under acidic conditions. Furthermore, 10 positions within the complementarity-determining region were also predicted to exhibit greater successes. Combined, these two accessible metrics could serve as the basis for a sequence-based engineering of pH-selective binding. This approach could be applied to most anticancer antibodies, which lack detailed structural characterization.

Redesigning an antibody H3 loop by virtual screening of a small library of human germline-derived sequences.

The design of superior biologic therapeutics, including antibodies and engineered proteins, involves optimizing their specific ability to bind to disease-related molecular targets. Previously, we developed and applied the Assisted Design of Antibody and Protein Therapeutics (ADAPT) platform for virtual affinity maturation of antibodies (Vivcharuk et al. in PLoS One 12(7):e0181490, https://doi.org/10.1371/journal.pone.0181490 , 2017). However, ADAPT is limited to point mutations of hot-spot residues in existing CDR loops. In this study, we explore the possibility of wholesale replacement of the entire H3 loop with no restriction to maintain the parental loop length. This complements other currently published studies that sample replacements for the CDR loops L1, L2, L3, H1 and H2. Given the immense sequence space theoretically available to H3, we focused on the virtual grafting of over 5000 human germline-derived H3 sequences from the IGMT/LIGM database increasing the diversity of the sequence space when compared to using crystalized H3 loop sequences. H3 loop conformations are generated and scored to identify optimized H3 sequences. Experimental testing of high-ranking H3 sequences grafted into the framework of the bH1 antibody against human VEGF-A led to the discovery of multiple hits, some of which had similar or better affinities relative to the parental antibody. In over 75% of the tested designs, the re-designed H3 loop contributed favorably to overall binding affinity. The hits also demonstrated good developability attributes such as high thermal stability and no aggregation. Crystal structures of select re-designed H3 variants were solved and indicated that although some deviations from predicted structures were seen in the more solvent accessible regions of the H3 loop, they did not significantly affect predicted affinity scores.

Hydrogen-deuterium exchange mass spectrometry reveals three unique binding responses of mAbs directed to the catalytic domain of hCAIX.

Human carbonic anhydrase (hCAIX), an extracellular enzyme that catalyzes the reversible hydration of CO2, is often overexpressed in solid tumors. This enzyme is instrumental in maintaining the survival of cancer cells in a hypoxic and acidic tumor microenvironment. Absent in most normal tissues, hCAIX is a promising therapeutic target for detection and treatment of solid tumors. Screening of a library of anti-hCAIX monoclonal antibodies (mAbs) previously identified three therapeutic candidates (mAb c2C7, m4A2 and m9B6) with distinct biophysical and functional characteristics. Selective binding to the catalytic domain was confirmed by yeast surface display and isothermal calorimetry, and deeper insight into the dynamic binding profiles of these mAbs upon binding were highlighted by bottom-up hydrogen-deuterium exchange mass spectrometry (HDX-MS). Here, a conformational and allosterically silent epitope was identified for the antibody-drug conjugate candidate c2C7. Unique binding profiles are described for both inhibitory antibodies, m4A2 and m9B6. M4A2 reduces the ability of the enzyme to hydrate CO2 by steric gating at the entrance of the catalytic cavity. Conversely, m9B6 disrupts the secondary structure that is necessary for substrate binding and hydration. The synergy of these two inhibitory mechanisms is demonstrated in in vitro activity assays and HDX-MS. Finally, the ability of m4A2 to modulate extracellular pH and intracellular metabolism is reported. By highlighting three unique modes by which hCAIX can be targeted, this study demonstrates both the utility of HDX-MS as an important tool in the characterization of anti-cancer biotherapeutics, and the underlying value of CAIX as a therapeutic target.

Structure-based engineering of pH-dependent antibody binding for selective targeting of solid-tumor microenvironment.

Recent development of monoclonal antibodies as mainstream anticancer agents demands further optimization of their safety for use in humans. Potent targeting and/or effector activities on normal tissues is an obvious toxicity concern. Optimization of specific tumor targeting could be achieved by taking advantage of the extracellular acidity of solid tumors relative to normal tissues. Here, we applied a structure-based computational approach to engineer anti-human epidermal growth factor receptor 2 (Her2) antibodies with selective binding in the acidic tumor microenvironment. We used an affinity maturation platform in which dual-pH histidine-scanning mutagenesis was implemented for pH selectivity optimization. Testing of a small set of designs for binding to the recombinant Her2 ectodomain led to the identification of antigen-binding fragment (Fab) variants with the desired pH-dependent binding behavior. Binding selectivity toward acidic pH was improved by as much as 25-fold relative to the parental bH1-Fab. In vitro experiments on cells expressing intact Her2 confirmed that designed variants formatted as IgG1/k full-size antibodies have high affinity and inhibit the growth of tumor spheroids at a level comparable to that of the benchmark anti-Her2 antibody trastuzumab (Herceptin®) at acidic pH, whereas these effects were significantly reduced at physiological pH. In contrast, both Herceptin and the parental bH1 antibody exhibited strong cell binding and growth inhibition irrespective of pH. This work demonstrates the feasibility of computational optimization of antibodies for selective targeting of the acidic environment such as that found in many solid tumors.

Role of the non-hypervariable FR3 D-E loop in single-domain antibody recognition of haptens and carbohydrates.

Single-domain antibodies (sdAbs), the variable domains of camelid heavy chain-only antibodies, are generally thought to poorly recognize nonproteinaceous small molecules and carbohydrates in comparison with conventional antibodies. However, the structures of anti-methotrexate, anti-triclocarban and anti-cortisol sdAbs revealed unexpected contributions of the non-hypervariable "CDR4" loop, formed between ß-strands D and E of framework region 3, in binding. Here, we investigated the potential role of CDR4 in sdAb binding to a hapten, 15-acetyl-deoxynivalenol (15-AcDON), and to carbohydrates. We constructed and panned a phage-displayed library in which CDR4 of the 15-AcDON-specific sdAb, NAT-267, was extended and randomized. From this library, we identified one sdAb, MA-232, bearing a 14-residue insertion in CDR4 and showing improved binding to 15-AcDON by ELISA and surface plasmon resonance. On the basis of these results, we constructed a second set of phage-displayed libraries in which the CDR4 and other regions of three hapten- or carbohydrate-binding sdAbs were diversified. With the goal of identifying sdAbs with novel glycan-binding specificities, we panned the library against four tumor-associated carbohydrate antigens but were unable to enrich binding phages. Thus, we conclude that while CDR4 may play a role in binding of some rare hapten-specific sdAbs, diversifying this region through molecular engineering is probably not a general solution to sdAb carbohydrate recognition in the absence of a paired VL domain.

The adaptor protein Ste50 directly modulates yeast MAPK signaling specificity through differential connections of its RA domain.

Discriminating among diverse environmental stimuli is critical for organisms to ensure their proper development, homeostasis, and survival. Saccharomyces cerevisiae regulates mating, osmoregulation, and filamentous growth using three different MAPK signaling pathways that share common components and therefore must ensure specificity. The adaptor protein Ste50 activates Ste11p, the MAP3K of all three modules. Its Ras association (RA) domain acts in both hyperosmolar and filamentous growth pathways, but its connection to the mating pathway is unknown. Genetically probing the domain, we found mutants that specifically disrupted mating or HOG-signaling pathways or both. Structurally these residues clustered on the RA domain, forming distinct surfaces with a propensity for protein-protein interactions. GFP fusions of wild-type (WT) and mutant Ste50p show that WT is localized to the shmoo structure and accumulates at the growing shmoo tip. The specifically pheromone response-defective mutants are severely impaired in shmoo formation and fail to localize ste50p, suggesting a failure of association and function of Ste50 mutants in the pheromone-signaling complex. Our results suggest that yeast cells can use differential protein interactions with the Ste50p RA domain to provide specificity of signaling during MAPK pathway activation.

Binding and functional profiling of antibody mutants guides selection of optimal candidates as antibody drug conjugates.

An increasingly appreciated conundrum in the discovery of antibody drug conjugates (ADCs) is that an antibody that was selected primarily for strong binding to its cancer target may not serve as an optimal ADC. In this study, we performed mechanistic cell-based experiments to determine the correlation between antibody affinity, avidity, internalization and ADC efficacy. We used structure-guided design to assemble a panel of antibody mutants with predicted Her2 affinities ranging from higher to lower relative to the parent antibody, Herceptin. These antibodies were ranked for binding via SPR and via flow-cytometry on high-Her2 SKOV3 cells and low-Her2 MCF7 cells, the latter acting as a surrogate for low-Her2 normal cells. A subpanel of variants, representative of different Her2-binding affinities (2 strong, 2 moderate and 3 weak), were further screened via high-content imaging for internalization efficacies in high versus low-Her2 cells. Finally, these antibodies were evaluated in ADC cytotoxicity screening assays (using DM1 and MMAE secondary antibodies) and as antibody-drug conjugates (DM1 and PNU159682). Our results identified specific but weak Her2-binding variants as optimal candidates for developing DM1 and PNU ADCs since they exhibited high potencies (low to sub-nM) in high-Her2 SKOV3 cells and low toxicities in low-Her2 cells. The 2 strong-affinity variants were highly potent in SKOV3 cells but also showed significant toxicities in low-Her2 cells and therefore are predicted to be toxic in normal tissues. Our findings show that pharmacological profiling of an antibody library in multiple binding and functional assays allows for selection of optimal ADCs.

Enhanced anti-metastatic bioactivity of an IGF-TRAP re-engineered to improve physicochemical properties.

The insulin-like growth factor (IGF) axis has been implicated in the progression of malignant disease and identified as a clinically important therapeutic target. Several IGF-1 receptor (IGF-1R) targeting drugs including humanized monoclonal antibodies have advanced to phase II/III clinical trials, but to date, have not progressed to clinical use, due, at least in part, to interference with insulin receptor signalling. We previously reported on the production of a soluble fusion protein consisting of the extracellular domain of human IGF-1R fused to the Fc portion of human IgG1 (first generation IGF-TRAP) that bound human IGF-1 and IGF-2 with a 3 log higher affinity than insulin. We showed that the IGF-TRAP had potent anti-cancer activity in several pre-clinical models of aggressive carcinomas. Here we report on the re-engineering of the IGF-TRAP with the aim of improving physicochemical properties and suitability for clinical applications. We show that cysteine-serine substitutions in the Fc hinge region of IGF-TRAP eliminated high-molecular-weight oligomerized species, while a further addition of a flexible linker, not only improved the pharmacokinetic profile, but also enhanced the therapeutic profile of the IGF-TRAP, as evaluated in an experimental colon carcinoma metastasis model. Dose-response profiles of the modified IGF-TRAPs correlated with their bio-availability profiles, as measured by the IGF kinase-receptor-activation (KIRA) assay, providing a novel, surrogate biomarker for drug efficacy. This study provides a compelling example of structure-based re-engineering of Fc-fusion-based biologics for better manufacturability that also significantly improved pharmacological parameters. It identifies the re-engineered IGF-TRAP as a potent anti-cancer therapeutic.

Assisted Design of Antibody and Protein Therapeutics (ADAPT).

Effective biologic therapeutics require binding affinities that are fine-tuned to their disease-related molecular target. The ADAPT (Assisted Design of Antibody and Protein Therapeutics) platform aids in the selection of mutants that improve/modulate the affinity of antibodies and other biologics. It uses a consensus z-score from three scoring functions and interleaves computational predictions with experimental validation, significantly enhancing the robustness of the design and selection of mutants. The platform was tested on three antibody Fab-antigen systems that spanned a wide range of initial binding affinities: bH1-VEGF-A (44 nM), bH1-HER2 (3.6 nM) and Herceptin-HER2 (0.058 nM). Novel triple mutants were obtained that exhibited 104-, 46- and 32-fold improvements in binding affinity for each system, respectively. Moreover, for all three antibody-antigen systems over 90% of all the intermediate single and double mutants that were designed and tested showed higher affinities than the parent sequence. The contributions of the individual mutants to the change in binding affinity appear to be roughly additive when combined to form double and triple mutants. The new interactions introduced by the affinity-enhancing mutants included long-range electrostatics as well as short-range nonpolar interactions. This diversity in the types of new interactions formed by the mutants was reflected in SPR kinetics that showed that the enhancements in affinities arose from increasing on-rates, decreasing off-rates or a combination of the two effects, depending on the mutation. ADAPT is a very focused search of sequence space and required only 20-30 mutants for each system to be made and tested to achieve the affinity enhancements mentioned above.

An engineered transforming growth factor ß (TGF-ß) monomer that functions as a dominant negative to block TGF-ß signaling.

The transforming growth factor ß isoforms, TGF-ß1, -ß2, and -ß3, are small secreted homodimeric signaling proteins with essential roles in regulating the adaptive immune system and maintaining the extracellular matrix. However, dysregulation of the TGF-ß pathway is responsible for promoting the progression of several human diseases, including cancer and fibrosis. Despite the known importance of TGF-ßs in promoting disease progression, no inhibitors have been approved for use in humans. Herein, we describe an engineered TGF-ß monomer, lacking the heel helix, a structural motif essential for binding the TGF-ß type I receptor (TßRI) but dispensable for binding the other receptor required for TGF-ß signaling, the TGF-ß type II receptor (TßRII), as an alternative therapeutic modality for blocking TGF-ß signaling in humans. As shown through binding studies and crystallography, the engineered monomer retained the same overall structure of native TGF-ß monomers and bound TßRII in an identical manner. Cell-based luciferase assays showed that the engineered monomer functioned as a dominant negative to inhibit TGF-ß signaling with a Ki of 20-70 nm Investigation of the mechanism showed that the high affinity of the engineered monomer for TßRII, coupled with its reduced ability to non-covalently dimerize and its inability to bind and recruit TßRI, enabled it to bind endogenous TßRII but prevented it from binding and recruiting TßRI to form a signaling complex. Such engineered monomers provide a new avenue to probe and manipulate TGF-ß signaling and may inform similar modifications of other TGF-ß family members.

Exhaustive docking and solvated interaction energy scoring: lessons learned from the SAMPL4 challenge.

We continued prospective assessments of the Wilma-solvated interaction energy (SIE) platform for pose prediction, binding affinity prediction, and virtual screening on the challenging SAMPL4 data sets including the HIV-integrase inhibitor and two host-guest systems. New features of the docking algorithm and scoring function are tested here prospectively for the first time. Wilma-SIE provides good correlations with actual binding affinities over a wide range of binding affinities that includes strong binders as in the case of SAMPL4 host-guest systems. Absolute binding affinities are also reproduced with appropriate training of the scoring function on available data sets or from comparative estimation of the change in target's vibrational entropy. Even when binding modes are known, SIE predictions lack correlation with experimental affinities within dynamic ranges below 2 kcal/mol as in the case of HIV-integrase ligands, but they correctly signaled the narrowness of the dynamic range. Using a common protein structure for all ligands can reduce the noise, while incorporating a more sophisticated solvation treatment improves absolute predictions. The HIV-integrase virtual screening data set consists of promiscuous weak binders with relatively high flexibility and thus it falls outside of the applicability domain of the Wilma-SIE docking platform. Despite these difficulties, unbiased docking around three known binding sites of the enzyme resulted in over a third of ligands being docked within 2 Å from their actual poses and over half of the ligands docked in the correct site, leading to better-than-random virtual screening results.

Sulfonation of glycopeptide antibiotics by sulfotransferase StaL depends on conformational flexibility of aglycone scaffold.

Although glycopeptide antibiotics (GPAs), including vancomycin and teicoplanin, represent the most important class of anti-infective agents in the treatment of serious gram-positive bacterial infections, their usefulness is threatened by the emergence of resistant strains. GPAs are complex natural products consisting of a heptapeptide skeleton assembled via nonribosomal peptide synthesis and constrained through multiple crosslinks, with diversity resulting from enzymatic modifications by a variety of tailoring enzymes, which can be used to produce GPA analogues that could overcome antibiotic resistance. GPA-modifying sulfotransferases are promising tools for generating the unique derivatives. Despite significant sequence and structural similarities, these sulfotransferases modify distinct side chains on the GPA scaffold. To provide insight into the spatial diversity of modifications, we have determined the crystal structure of the ternary complex of bacterial sulfotransferase StaL with the cofactor product 3'-phosphoadenosine 5'-phosphate and desulfo-A47934 aglycone substrate. Desulfo-A47934 binds with the hydroxyl group on the 4-hydroxyphenylglycine in residue 1 directed toward the 3'-phosphoadenosine 5'-phosphate and hydrogen-bonded to the catalytic His67. Homodimeric StaL can accommodate GPA substrate in only one of the two active sites because of potential steric clashes. Importantly, the aglycone substrate demonstrates a flattened conformation, in contrast to the cup-shaped structures observed previously. Analysis of the conformations of this scaffold showed that despite the apparent rigidity due to crosslinking between the side chains, the aglycone scaffold displays substantial flexibility, important for enzymatic modifications by the GPA-tailoring enzymes. We also discuss the potential of using the current structural information in generating unique GPA derivatives.

Engineering and therapeutic application of single-chain bivalent TGF-ß family traps.

Deregulation of TGF-ß superfamily signaling is a causative factor in many diseases. Here we describe a protein engineering strategy for the generation of single-chain bivalent receptor traps for TGF-ß superfamily ligands. Traps were assembled using the intrinsically disordered regions flanking the structured binding domain of each receptor as "native linkers" between two binding domains. This yields traps that are approximately threefold smaller than antibodies and consists entirely of native receptor sequences. Two TGF-ß type II receptor-based, single-chain traps were designed, termed (TßRII)2 and (TßRIIb)2, that have native linker lengths of 35 and 60 amino acids, respectively. Both single-chain traps exhibit a 100 to 1,000 fold higher in vitro ligand binding and neutralization activity compared with the monovalent ectodomain (TßRII-ED), and a similar or slightly better potency than pan-TGF-ß-neutralizing antibody 1D11 or an Fc-fused receptor trap (TßRII-Fc). Despite its short in vivo half-life (<1 hour), which is primarily due to kidney clearance, daily injections of the (TßRII)2 trap reduced the growth of 4T1 tumors in BALB/c mice by 50%, an efficacy that is comparable with 1D11 (dosed thrice weekly). In addition, (TßRII)2 treatment of mice with established 4T1 tumors (100 mm(3)) significantly inhibited further tumor growth, whereas the 1D11 antibody did not. Overall, our results indicate that our rationally designed bivalent, single-chain traps have promising therapeutic potential.

Selective inhibition of USP7.

The deubiquitinating enzyme USP7 is an emerging oncology and antiviral target. Reverdy et al., in this issue of Chemistry & Biology, disclose the first small-molecule inhibitor selective for USP7, which recapitulates its knockdown in cancer cells and hence demonstrates the therapeutic feasibility of USP7 inhibitors.

Binding the atypical RA domain of Ste50p to the unfolded Opy2p cytoplasmic tail is essential for the high-osmolarity glycerol pathway.

Activation of the high-osmolarity glycerol (HOG) pathway for osmoregulation in the yeast Saccharomyces cerevisiae involves interaction of the adaptor Ste50p with the cytoplasmic tail of single-transmembrane protein Opy2p. We have determined the solution structure of the Ste50p-RA (Ras association) domain, and it shows an atypical RA fold lacking the beta1 and beta2 strands of the canonical motif. Although the core of the RA domain is fully functional in the pheromone response, an additional region is required for the HOG pathway activation. Two peptide motifs within the intrinsically disordered cytoplasmic tail of Opy2p defined by NMR spectroscopy physically interact with the Step50p-RA domain. These Opy2p-derived peptides bind overlapping regions of the Step50p-RA domain with similarly weak affinities, suggesting a multivalent interaction of these proteins as a crucial point of control of the HOG pathway. As well, overall selection of signaling pathways depends on functionally distinct regions of the Ste50p-RA domain, implicating this element in the control of global regulatory decisions.

Structure and function of the glycopeptide N-methyltransferase MtfA, a tool for the biosynthesis of modified glycopeptide antibiotics.

There is a considerable interest in the modification of existing antibiotics to generate new antimicrobials. Glycopeptide antibiotics (GPAs) are effective against serious Gram-positive bacterial pathogens including methicillin-resistant Staphylococcus aureus. However, resistance to these antibiotics is becoming a serious problem requiring new strategies. We show that the Amycolatopsis orientalis (S)-adenosyl-L-methionine-dependent methyltransferase MtfA, from the vancomycin-class GPA chloroeremomycin biosynthetic pathway, catalyzes in vivo and in vitro methyl transfer to generate methylated GPA derivatives of the teicoplanin class. The crystal structure of MtfA complexed with (S)-adenosyl-L-methionine, (S)-adenosylhomocysteine, or sinefungin inhibitor, coupled with mutagenesis, identified His228 as a likely general base required for methyl transfer to the N terminus of the glycopeptide. Computational docking and molecular dynamics simulations were used to model binding of demethyl-vancomycin aglycone to MtfA. These results demonstrate its utility as a tool for engineering methylated analogs of GPAs.

The propeptide of cruzipain--a potent selective inhibitor of the trypanosomal enzymes cruzipain and brucipain, and of the human enzyme cathepsin F.

Papain-like cysteine proteases of pathogenic protozoa play important roles in parasite growth, differentiation and host cell invasion. The main cysteine proteases of Trypanosoma cruzi (cruzipain) and of Trypanosoma brucei (brucipain) are validated targets for the development of new chemotherapies. These proteases are synthesized as precursors and activated upon removal of the N-terminal prodomain. Here we report potent and selective inhibition of cruzipain and brucipain by the recombinant full-length prodomain of cruzipain. The propeptide did not inhibit human cathepsins S, K or B or papain at the tested concentrations, and moderately inhibited human cathepsin V. Human cathepsin F was very efficiently inhibited (K(i) of 32 pm), an interesting finding indicating that cruzipain propeptide is able to discriminate cathepsin F from other cathepsin L-like enzymes. Comparative structural modeling and analysis identified the interaction between the beta1p-alpha3p loop of the propeptide and the propeptide-binding loop of mature enzymes as a plausible cause of the observed inhibitory selectivity.

Solvated interaction energy (SIE) for scoring protein-ligand binding affinities. 1. Exploring the parameter space.

We present a binding free energy function that consists of force field terms supplemented by solvation terms. We used this function to calibrate the solvation model along with the binding interaction terms in a self-consistent manner. The motivation for this approach was that the solute dielectric-constant dependence of calculated hydration gas-to-water transfer free energies is markedly different from that of binding free energies (J. Comput. Chem. 2003, 24, 954). Hence, we sought to calibrate directly the solvation terms in the context of a binding calculation. The five parameters of the model were systematically scanned to best reproduce the absolute binding free energies for a set of 99 protein-ligand complexes. We obtained a mean unsigned error of 1.29 kcal/mol for the predicted absolute binding affinity in a parameter space that was fairly shallow near the optimum. The lowest errors were obtained with solute dielectric values of Din = 20 or higher and scaling of the intermolecular van der Waals interaction energy by factors ranging from 0.03 to 0.15. The high apparent Din and strong van der Waals scaling may reflect the anticorrelation of the change in solvated potential energy and configurational entropy, that is, enthalpy-entropy compensation in ligand binding (Biophys. J. 2004, 87, 3035-3049). Five variations of preparing the protein-ligand data set were explored in order to examine the effect of energy refinement and the presence of bound water on the calculated results. We find that retaining water in the final protein structure used for calculating the binding free energy is not necessary to obtain good results; that is the continuum solvation model is sufficient. Virtual screening enrichment studies on estrogen receptor and thymidine kinase showed a good ability of the binding free energy function to recover true hits in a collection of decoys.