Cycloaddition of Phenyltriazolinedione with Carbazole-Alkynes and Yne-Carbamates to Access Diazacyclobutenes.

Previously, we described the synthesis of stable, bicyclic examples of the rather rare diazacyclobutene (DCB) motif by means of a cycloaddition between triazolinediones and electron-rich thiolated alkynes. Here, we report the investigation of the cycloaddition of triazolinediones with related electron-rich yne-carbamates and carbazole-alkynes. Bicyclic DCBs arising from yne-carbamates were isolated in 8-65% yield, while those arising from carbazole-alkynes were isolated in 28-59% yield. Mechanistic studies and characterization of isolable byproducts shed light on the underlying issues leading to poor to moderate yields.

Enolase inhibitors as therapeutic leads for Naegleria fowleri infection.

Infections with the pathogenic free-living amoebae Naegleria fowleri can lead to life-threatening illnesses including catastrophic primary amebic meningoencephalitis (PAM). Efficacious treatment options for these infections are lacking and the mortality rate remains >95% in the US. Glycolysis is very important for the infectious trophozoite lifecycle stage and inhibitors of glucose metabolism have been found to be toxic to the pathogen. Recently, human enolase 2 (ENO2) phosphonate inhibitors have been developed as lead agents to treat glioblastoma multiforme (GBM). These compounds, which cure GBM in a rodent model, are well-tolerated in mammals because enolase 1 (ENO1) is the predominant isoform used systemically. Here, we describe findings that demonstrate that these agents are potent inhibitors of N. fowleri ENO ( Nf ENO) and are lethal to amoebae. In particular, (1-hydroxy-2-oxopiperidin-3-yl) phosphonic acid (HEX) was a potent enzyme inhibitor (IC 50 value of 0.14 ± 0.04 µM) that was toxic to trophozoites (EC 50 value of 0.21 ± 0.02 µM) while the reported CC 50 was >300 µM. Molecular docking simulation revealed that HEX binds strongly to the active site of Nf ENO with a binding affinity of -8.6 kcal/mol. Metabolomic studies of parasites treated with HEX revealed a 4.5 to 78-fold accumulation of glycolytic intermediates upstream of Nf ENO. Last, nasal instillation of HEX increased longevity of amoebae-infected rodents. Two days after infection, animals were treated for 10 days with 3 mg/kg HEX, followed by one week of observation. At the conclusion of the experiment, eight of 12 HEX-treated animals remained alive (resulting in an indeterminable median survival time) while one of 12 vehicle-treated rodents remained, yielding a median survival time of 10.9 days. Brains of six of the eight survivors were positive for amoebae, suggesting the agent at the tested dose suppressed, but did not eliminate, infection. These findings suggest that HEX is a promising lead for the treatment of PAM.

Glucose metabolism in the pathogenic free-living amoebae: Tempting targets for treatment development.

Pathogenic free-living amoebae (pFLA) are single-celled eukaryotes responsible for causing intractable infections with high morbidity and mortality in humans and animals. Current therapeutic approaches include cocktails of antibiotic, antifungal, and antimicrobial compounds. Unfortunately, the efficacy of these can be limited, driving the need for the discovery of new treatments. Pan anti-amebic agents would be ideal; however, identifying these agents has been a challenge, likely due to the limited evolutionary relatedness of the different pFLA. Here, we discuss the potential of targeting amoebae glucose metabolic pathways as the differences between pFLA and humans suggest specific inhibitors could be developed as leads for new therapeutics.

A cathepsin C-like protease mediates the post-translation modification of Toxoplasma gondii secretory proteins for optimal invasion and egress.

Microbial pathogens use proteases for their infections, such as digestion of proteins for nutrients and activation of their virulence factors. As an obligate intracellular parasite, Toxoplasma gondii must invade host cells to establish its intracellular propagation. To facilitate invasion, the parasites secrete invasion effectors from microneme and rhoptry, two unique organelles in apicomplexans. Previous work has shown that some micronemal invasion effectors experience a series of proteolytic cleavages within the parasite's secretion pathway for maturation, such as the aspartyl protease (TgASP3) and the cathepsin L-like protease (TgCPL), localized within the post-Golgi compartment and the endolysosomal system, respectively. Furthermore, it has been shown that the precise maturation of micronemal effectors is critical for Toxoplasma invasion and egress. Here, we show that an endosome-like compartment (ELC)-residing cathepsin C-like protease (TgCPC1) mediates the final trimming of some micronemal effectors, and its loss further results in defects in the steps of invasion, egress, and migration throughout the parasite's lytic cycle. Notably, the deletion of TgCPC1 completely blocks the activation of subtilisin-like protease 1 (TgSUB1) in the parasites, which globally impairs the surface-trimming of many key micronemal invasion and egress effectors. Additionally, we found that Toxoplasma is not efficiently inhibited by the chemical inhibitor targeting the malarial CPC ortholog, suggesting that these cathepsin C-like orthologs are structurally different within the apicomplexan phylum. Collectively, our findings identify a novel function of TgCPC1 in processing micronemal proteins within the Toxoplasma parasite's secretory pathway and expand the understanding of the roles of cathepsin C protease. IMPORTANCE Toxoplasma gondii is a microbial pathogen that is well adapted for disseminating infections. It can infect virtually all warm-blooded animals. Approximately one-third of the human population carries toxoplasmosis. During infection, the parasites sequentially secrete protein effectors from the microneme, rhoptry, and dense granule, three organelles exclusively found in apicomplexan parasites, to help establish their lytic cycle. Proteolytic cleavage of these secretory proteins is required for the parasite's optimal function. Previous work has revealed that two proteases residing within the parasite's secretory pathway cleave micronemal and rhoptry proteins, which mediate parasite invasion and egress. Here, we demonstrate that a cathepsin C-like protease (TgCPC1) is involved in processing several invasion and egress effectors. The genetic deletion of TgCPC1 prevented the complete maturation of some effectors in the parasites. Strikingly, the deletion led to a full inactivation of one surface-anchored protease, which globally impaired the trimming of some key micronemal proteins before secretion. Therefore, this finding represents a novel post-translational mechanism for the processing of virulence factors within microbial pathogens.

Exploring the interactive mechanism of acarbose with the amylase SusG in the starch utilization system of the human gut symbiont Bacteroides thetaiotaomicron through molecular modeling.

The α-amylase, SusG, is a principal component of the Bacteroides thetaiotaomicron (Bt) starch utilization system (Sus) used to metabolize complex starch molecules in the human gastrointestinal (GI) tract. We previously reported the non-microbicidal growth inhibition of Bt by the acarbose-mediated arrest of the Sus as a potential therapeutic strategy. Herein, we report a computational approach using density functional theory (DFT), molecular docking, and molecular dynamics (MD) simulation to explore the interactive mechanism between acarbose and SusG at the atomic level in an effort to understand how acarbose shuts down the Bt Sus. The docking analysis reveals that acarbose binds orthosterically to SusG with a binding affinity of -8.3 kcal/mol. The MD simulation provides evidence of conformational variability of acarbose at the active site of SusG and also suggests that acarbose interacts with the main catalytic residues via a general acid-base double-displacement catalytic mechanism. These results suggest that small molecule competitive inhibition against the SusG protein could impact the entire Bt Sus and eliminate or reduce the system's ability to metabolize starch. This computational strategy could serve as a potential avenue for structure-based drug design to discover other small molecules capable of inhibiting the Sus of Bt with high potency, thus providing a holistic approach for selective modulation of the GI microbiota.

A cathepsin C-like protease post-translationally modifies Toxoplasma gondii secretory proteins for optimal invasion and egress.

Microbial pathogens use proteases for their infections, such as digestion of proteins for nutrients and activation of their virulence factors. As an obligate intracellular parasite, Toxoplasma gondii must invade host cells to establish its intracellular propagation. To facilitate invasion, the parasites secrete invasion effectors from microneme and rhoptry, two unique organelles in apicomplexans. Previous work has shown that some micronemal invasion effectors experience a series of proteolytic cleavages within the parasite's secretion pathway for maturation, such as the aspartyl protease (TgASP3) and the cathepsin L-like protease (TgCPL), localized within the post-Golgi compartment (1) and the endolysosomal system (2), respectively. Furthermore, it has been shown that the precise maturation of micronemal effectors is critical for Toxoplasma invasion and egress (1). Here, we show that an endosome-like compartment (ELC)-residing cathepsin C-like protease (TgCPC1) mediates the final trimming of some micronemal effectors, and its loss further results in defects in the steps of invasion, egress, and migration throughout the parasite's lytic cycle. Notably, the deletion of TgCPC1 completely blocks the activation of subtilisin-like protease 1 (TgSUB1) in the parasites, which globally impairs the surface-trimming of many key micronemal invasion and egress effectors. Additionally, we found that TgCPC1 was not efficiently inhibited by the chemical inhibitor targeting its malarial ortholog, suggesting that these cathepsin C-like orthologs are structurally different within the apicomplexan phylum. Taken together, our findings identify a novel function of TgCPC1 in the processing of micronemal proteins within the secretory pathway of Toxoplasma parasites and expand the understanding of the roles of cathepsin C protease. IMPORTANCE: Toxoplasma gondii is a microbial pathogen that is well adapted for disseminating infections. It can infect virtually all warm-blooded animals. Approximately one-third of the human population carries toxoplasmosis. During infection, the parasites sequentially secrete protein effectors from the microneme, rhoptry, and dense granule, three organelles exclusively found in apicomplexan parasites, to help establish their lytic cycle. Proteolytic cleavage of these secretory proteins is required for the parasite's optimal function. Previous work has revealed that two proteases residing within the parasite's secretory pathway cleave micronemal and rhoptry proteins, which mediate parasite invasion and egress. Here, we demonstrate that a cathepsin C-like protease (TgCPC1) is involved in processing several invasion and egress effectors. The genetic deletion of TgCPC1 prevented the complete maturation of some effectors in the parasites. Strikingly, the deletion led to a full inactivation of one surface-anchored protease, which globally impaired the trimming of some key micronemal proteins before secretion. Therefore, this finding represents a novel post-translational mechanism for the processing of virulence factors within microbial pathogens.

Correlated Mutation in the Evolution of Catalysis in Uracil DNA Glycosylase Superfamily.

Enzymes in Uracil DNA glycosylase (UDG) superfamily are essential for the removal of uracil. Family 4 UDGa is a robust uracil DNA glycosylase that only acts on double-stranded and single-stranded uracil-containing DNA. Based on mutational, kinetic and modeling analyses, a catalytic mechanism involving leaving group stabilization by H155 in motif 2 and water coordination by N89 in motif 3 is proposed. Mutual Information analysis identifies a complexed correlated mutation network including a strong correlation in the EG doublet in motif 1 of family 4 UDGa and in the QD doublet in motif 1 of family 1 UNG. Conversion of EG doublet in family 4 Thermus thermophilus UDGa to QD doublet increases the catalytic efficiency by over one hundred-fold and seventeen-fold over the E41Q and G42D single mutation, respectively, rectifying the strong correlation in the doublet. Molecular dynamics simulations suggest that the correlated mutations in the doublet in motif 1 position the catalytic H155 in motif 2 to stabilize the leaving uracilate anion. The integrated approach has important implications in studying enzyme evolution and protein structure and function.

A structural determinant in the uracil DNA glycosylase superfamily for the removal of uracil from adenine/uracil base pairs.

The uracil DNA glycosylase superfamily consists of several distinct families. Family 2 mismatch-specific uracil DNA glycosylase (MUG) from Escherichia coli is known to exhibit glycosylase activity on three mismatched base pairs, T/U, G/U and C/U. Family 1 uracil N-glycosylase (UNG) from E. coli is an extremely efficient enzyme that can remove uracil from any uracil-containing base pairs including the A/U base pair. Here, we report the identification of an important structural determinant that underlies the functional difference between MUG and UNG. Substitution of a Lys residue at position 68 with Asn in MUG not only accelerates the removal of uracil from mismatched base pairs but also enables the enzyme to gain catalytic activity on A/U base pairs. Binding and kinetic analysis demonstrate that the MUG-K68N substitution results in enhanced ground state binding and transition state interactions. Molecular modeling reveals that MUG-K68N, UNG-N123 and family 5 Thermus thermophiles UDGb-A111N can form bidentate hydrogen bonds with the N3 and O4 moieties of the uracil base. Genetic analysis indicates the gain of function for A/U base pairs allows the MUG-K68N mutant to remove uracil incorporated into the genome during DNA replication. The implications of this study in the origin of life are discussed.

Specificity and catalytic mechanism in family 5 uracil DNA glycosylase.

UDGb belongs to family 5 of the uracil DNA glycosylase (UDG) superfamily. Here, we report that family 5 UDGb from Thermus thermophilus HB8 is not only a uracil DNA glycosyase acting on G/U, T/U, C/U, and A/U base pairs, but also a hypoxanthine DNA glycosylase acting on G/I, T/I, and A/I base pairs and a xanthine DNA glycosylase acting on all double-stranded and single-stranded xanthine-containing DNA. Analysis of potentials of mean force indicates that the tendency of hypoxanthine base flipping follows the order of G/I > T/I, A/I > C/I, matching the trend of hypoxanthine DNA glycosylase activity observed in vitro. Genetic analysis indicates that family 5 UDGb can also act as an enzyme to remove uracil incorporated into DNA through the existence of dUTP in the nucleotide pool. Mutational analysis coupled with molecular modeling and molecular dynamics analysis reveals that although hydrogen bonding to O2 of uracil underlies the UDG activity in a dissociative fashion, Tth UDGb relies on multiple catalytic residues to facilitate its excision of hypoxanthine and xanthine. This study underscores the structural and functional diversity in the UDG superfamily.

The evolution of cefotaximase activity in the TEM ß-lactamase.

The development of a molecular-level understanding of drug resistance through ß-lactamase is critical not only in designing newer-generation antibacterial agents but also in providing insight into the evolutionary mechanisms of enzymes in general. In the present study, we have evaluated the effect of four drug resistance mutations (A42G, E104K, G238S, and M182T) on the cefotaximase activity of the TEM-1 ß-lactamase. Using computational methods, including docking and molecular mechanics calculations, we have been able to correctly identify the relative order of catalytic activities associated with these four single point mutants. Further analyses suggest that the changes in catalytic efficiency for mutant enzymes are correlated to structural changes within the binding site. Based on the energetic and structural analyses of the wild-type and mutant enzymes, structural rearrangement is suggested as a mechanism of evolution of drug resistance through TEM ß-lactamase. The present study not only provides molecular-level insight into the effect of four drug resistance mutations on the structure and function of the TEM ß-lactamase but also establishes a foundation for a future molecular-level analysis of complete evolutionary trajectory for this class of enzymes.

New family of deamination repair enzymes in uracil-DNA glycosylase superfamily.

DNA glycosylases play a major role in the repair of deaminated DNA damage. Previous investigations identified five families within the uracil-DNA glycosylase (UDG) superfamily. All enzymes within the superfamily studied thus far exhibit uracil-DNA glycosylase activity. Here we identify a new class of DNA glycosylases in the UDG superfamily that lacks UDG activity. Instead, these enzymes act as hypoxanthine-DNA glycosylases in vitro and in vivo. Molecular modeling and structure-guided mutational analysis allowed us to identify a unique catalytic center in this class of DNA glycosylases. Based on unprecedented biochemical properties and phylogenetic analysis, we propose this new class of DNA repair glycosylases that exists in bacteria, archaea, and eukaryotes as family 6 and designate it as the hypoxanthine-DNA glycosylase family. This study demonstrates the structural evolvability that underlies substrate specificity and catalytic flexibility in the evolution of enzymatic function.

The evolution of catalytic function in the HIV-1 protease.

The evolution of species is a complex phenomenon based on the optimization of a multidimensional function referred to as fitness. At the level of biomolecular evolution, the fitness function can be reduced to include physiochemical properties relevant to the biological function of a particular molecule. In this work, questions involving the physical-chemical mechanisms underlying the evolution of HIV-1 protease are addressed through molecular simulation and subsequent analysis of thermodynamic properties related to the activity of the enzyme. Specifically, the impact of 40 single amino acid mutations on the binding affinity toward the matrix/capsid (MA/CA) substrate and corresponding transition state intermediate has been characterized using a molecular mechanics Poisson-Boltzmann surface area approach. We demonstrate that this approach is capable of extracting statistically significant information relevant to experimentally determined catalytic activity. Further, no correlation was observed between the effect of mutations on substrate and transition state binding, suggesting independent evolutionary pathways toward optimizing substrate specificity and catalytic activity. In addition, a detailed analysis of calculated binding affinity data suggests that ground-state destabilization (reduced binding affinity for the substrate) could be a contributing factor in the evolutionary optimization of HIV-1 protease. A numerical model is developed to demonstrate that ground-state destabilization is a valid mechanism for activity optimization given the high concentrations of substrate experienced by the functional enzyme in vivo.

Analyzing the robustness of the MM/PBSA free energy calculation method: application to DNA conformational transitions.

The ability to predict and characterize free energy differences associated with conformational equilibria or the binding of biomolecules is vital to understanding the molecular basis of many important biological functions. As biological studies focus on larger molecular complexes and properties of the genome, proteome, and interactome, the development and characterization of efficient methods for calculating free energy becomes increasingly essential. The aim of this study is to examine the robustness of the end-point free energy method termed the molecular mechanics Poisson-Boltzmann solvent accessible surface area (MM/PBSA) method. Specifically, applications of MM/PBSA to the conformational equilibria of nucleic acid (NA) systems are explored. This is achieved by comparing A to B form DNA conformational free energy differences calculated using MM/PBSA with corresponding free energy differences determined with a more rigorous and time-consuming umbrella sampling algorithm. In addition, the robustness of NA MM/PBSA calculations is also evaluated in terms of the sensitivity towards the choice of force field and the choice of solvent model used during conformational sampling. MM/PBSA calculations of the free energy difference between A-form and B-form DNA are shown to be in very close agreement with the PMF result determined using an umbrella sampling approach. Further, it is found that the MM/PBSA conformational free energy differences were also in agreement using either the CHARMM or AMBER force field. The influence of ionic strength on conformational stability was particularly insensitive to the choice of force field. Finally, it is also shown that the use of a generalized Born implicit solvent during conformational sampling results in free energy estimates that deviate slightly from those obtained using explicitly solvated MD simulations in these NA systems.

Identification of Escherichia coli mismatch-specific uracil DNA glycosylase as a robust xanthine DNA glycosylase.

The gene for the mismatch-specific uracil DNA glycosylase (MUG) was identified in the Escherichia coli genome as a sequence homolog of the human thymine DNA glycosylase with activity against mismatched uracil base pairs. Examination of cell extracts led us to detect a previously unknown xanthine DNA glycosylase (XDG) activity in E. coli. DNA glycosylase assays with purified enzymes indicated the novel XDG activity is attributable to MUG. Here, we report a biochemical characterization of xanthine DNA glycosylase activity in MUG. The wild type MUG possesses more robust activity against xanthine than uracil and is active against all xanthine-containing DNA (C/X, T/X, G/X, A/X and single-stranded X). Analysis of potentials of mean force indicates that the double-stranded xanthine base pairs have a relatively narrow energetic difference in base flipping, whereas the tendency for uracil base flipping follows the order of C/U > G/U > T/U > A/U. Site-directed mutagenesis performed on conserved motifs revealed that Asn-140 and Ser-23 are important determinants for XDG activity in E. coli MUG. Molecular modeling and molecular dynamics simulations reveal distinct hydrogen-bonding patterns in the active site of E. coli MUG that account for the specificity differences between E. coli MUG and human thymine DNA glycosylase as well as that between the wild type MUG and the Asn-140 and Ser-23 mutants. This study underscores the role of the favorable binding interactions in modulating the specificity of DNA glycosylases.

Thermodynamic resolution: how do errors in modeled protein structures affect binding affinity predictions?

The present study addresses the effect of structural distortion, caused by protein modeling errors, on calculated binding affinities toward small molecules. The binding affinities to a total of 300 distorted structures based on five different protein-ligand complexes were evaluated to establish a broadly applicable relationship between errors in protein structure and errors in calculated binding affinities. Relatively accurate protein models (less than 2 A RMSD within the binding site) demonstrate a 14.78 (+/-7.5)% deviation in binding affinity from that calculated by using the corresponding crystal structure. For structures of 2-3 A, 3-4 A, and >4 A RMSD within the binding site, the error in calculated binding affinity increases to 20.8 (+/-5.98), 22.79 (+/-11.3), and 29.43 (+/-11.47)%, respectively. The results described here may be used in combination with other tools to evaluate the utility of modeled protein structures for drug development or other ligand-binding studies.

Insights from xanthine and uracil DNA glycosylase activities of bacterial and human SMUG1: switching SMUG1 to UDG.

Single-strand-selective monofunctional uracil DNA glycosylase (SMUG1) belongs to Family 3 of the uracil DNA glycosylase (UDG) superfamily. Here, we report that a bacterial SMUG1 ortholog in Geobacter metallireducens (Gme) and the human SMUG1 enzyme are not only UDGs but also xanthine DNA glycosylases (XDGs). In addition, mutational analysis and molecular dynamics (MD) simulations of Gme SMUG1 identify important structural determinants in conserved motifs 1 and 2 for XDG and UDG activities. Mutations at M57 (M57L) and H210 (H210G, H210M, and H210N), both of which are involved in interactions with the C2 carbonyl oxygen in uracil or xanthine, cause substantial reductions in XDG and UDG activities. Increased selectivity is achieved in the A214R mutant of Gme SMUG1, which corresponds to a position involved in base flipping. This mutation results in an activity profile resembling a human SMUG1-like enzyme as exemplified by the retention of UDG activity on mismatched base pairs and weak XDG activity. MD simulations indicate that M57L increases the flexibility of the motif 2 loop region and specifically A214, which may account for the reduced catalytic activity. G60Y completely abolishes XDG and UDG activity, which is consistent with a modeled structure in which G60Y blocks the entry of either xanthine or uracil to the base binding pocket. Most interestingly, a proline substitution at the G63 position switches the Gme SMUG1 enzyme to an exclusive UDG as demonstrated by the uniform excision of uracil in both double-stranded and single-stranded DNA and the complete loss of XDG activity. MD simulations indicate that a combination of a reduced free volume and altered flexibility in the active-site loops may underlie the dramatic effects of the G63P mutation on the activity profile of SMUG1. This study offers insights on the important role that modulation of conformational flexibility may play in defining specificity and catalytic efficiency.

Comparison of solvation-effect methods for the simulation of peptide interactions with a hydrophobic surface.

In this study we investigated the interaction behavior between thirteen different small peptides and a hydrophobic surface using three progressively more complex methods of representing solvation effects: a united-atom implicit solvation method [CHARMM 19 force field (C19) with Analytical Continuum Electrostatics (ACE)], an all-atom implicit solvation method (C22 with GBMV), and an all-atom explicit solvation method (C22 with TIP3P). The adsorption behavior of each peptide was characterized by the calculation of the potential of mean force as a function of peptide-surface separation distance. The results from the C22/TIP3P model suggest that hydrophobic peptides exhibit relatively strong adsorption behavior, polar and positively-charged peptides exhibit negligible to relatively weak favorable interactions with the surface, and negatively-charged peptides strongly resist adsorption. Compared to the TIP3P model, the ACE and GBMV implicit solvent models predict much stronger attractions for the hydrophobic peptides as well as stronger repulsions for the negatively-charged peptides on the CH(3)-SAM surface. These comparisons provide a basis from which each of these implicit solvation methods may be reparameterized to provide closer agreement with explicitly represented solvation in simulations of peptide and protein adsorption to functionalized surfaces.

An electrostatic basis for the stability of thermophilic proteins.

Two factors provide key contributions to the stability of thermophilic proteins relative to their mesophilic homologues: electrostatic interactions of charged residues in the folded state and the dielectric response of the folded protein. The dielectric response for proteins in a "thermophilic series" globally modulates the thermal stability of its members, with the calculated dielectric constant for the protein increasing from mesophiles to hyperthermophiles. This variability results from differences in the distribution of charged residues on the surface of the protein, in agreement with structural and genetic observations. Furthermore, the contribution of electrostatic interactions to the stability of the folded state is more favorable for thermophilic proteins than for their mesophilic homologues. This leads to the conclusion that electrostatic interactions play an important role in determining the stability of proteins at high temperatures. The interplay between electrostatic interactions and dielectric response also provides further rationalization for the enhanced stability of thermophilic proteins with respect to cold-denaturation. Taken together, the distribution of charged residues and their fluctuations have been shown to be factors in modulating protein stability over the entire range of biologically relevant temperatures.

Native atom types for knowledge-based potentials: application to binding energy prediction.

Knowledge-based potentials have been found useful in a variety of biophysical studies of macromolecules. Recently, it has also been shown in self-consistent studies that it is possible to extract quantities consistent with pair potentials from model structural databases. In this study, we attempt to extend the results obtained from these self-consistent studies toward the extraction of realistic pair potentials from the Protein Data Bank (PDB). The new method utilizes a clustering approach to define atom types within the PDB consistent with the optimal effective pairwise potential. The method has been integrated into the SMoG drug design package, resulting in an improved approach for the rapid and accurate estimation of binding affinities from structural information. Using this approach, it is possible to generate simple knowledge-based potentials that correlate (R = 0.61) with experimental binding affinities in a database of 118 diverse complexes. Furthermore, predictions performed on a random 1/3 of the database consistently show an average unsigned error of 1.5 log Ki units. It is also possible to generate specialized knowledge-based potentials, targeted to specific protein families. This approach is capable of generating potentials that correlate strongly with experimental binding affinities within these families (R = 0.8-0.9). Predictions on 1/3 of these family databases yield average unsigned errors ranging from 1.1 to 1.3 log Ki units. In summary, we describe a physically motivated approach to optimizing knowledge-based potentials for binding energy prediction that can be integrated into a variety of stages within a lead discovery protocol.

Structural mining: self-consistent design on flexible protein-peptide docking and transferable binding affinity potential.

A flexible protein-peptide docking method has been designed to consider not only ligand flexibility but also the flexibility of the protein. The method is based on a Monte Carlo annealing process. Simulations with a distance root-mean-square (dRMS) virtual energy function revealed that the flexibility of protein side chains was as important as ligand flexibility for successful protein-peptide docking. On the basis of mean field theory, a transferable potential was designed to evaluate distance-dependent protein-ligand interactions and atomic solvation energies. The potential parameters were developed using a self-consistent process based on only 10 known complex structures. The effectiveness of each intermediate potential was judged on the basis of a Z score, approximating the gap between the energy of the native complex and the average energy of a decoy set. The Z score was determined using experimentally determined native structures and decoys generated by docking with the intermediate potentials. Using 6600 generated decoys and the Z score optimization criterion proposed in this work, the developed potential yielded an acceptable correlation of R(2) = 0.77, with binding free energies determined for known MHC I complexes (Class I Major Histocompatibility protein HLA-A(*)0201) which were not present in the training set. Test docking on 25 complexes further revealed a significant correlation between energy and dRMS, important for identifying native-like conformations. The near-native structures always belonged to one of the conformational classes with lower predicted binding energy. The lowest energy docked conformations are generally associated with near-native conformations, less than 3.0 Angstrom dRMS (and in many cases less than 1.0 Angstrom) from the experimentally determined structures.