Identification of novel drug scaffolds for inhibition of SARS-CoV 3-Chymotrypsin-like protease using virtual and high-throughput screenings

We have used a combination of virtual screening (VS) and high-throughput screening (HTS) techniques to identify novel, non-peptidic small molecule inhibitors against human SARS-CoV 3CLpro. A structure-based VS approach integrating docking and pharmacophore based methods was employed to computationally screen 621,000 compounds from the ZINC library. The screening protocol was validated using known 3CLpro inhibitors and was optimized for speed, improved selectivity, and for accommodating receptor flexibility. Subsequently, a fluorescence-based enzymatic HTS assay was developed and optimized to experimentally screen approximately 41,000 compounds from four structurally diverse libraries chosen mainly based on the VS results. False positives from initial HTS hits were eliminated by a secondary orthogonal binding analysis using surface plasmon resonance (SPR). The campaign identified a reversible small molecule inhibitor exhibiting mixed-type inhibition with a K_i value of 11.1 μM. Together, these results validate our protocols as suitable approaches to screen virtual and chemical libraries, and the newly identified compound reported in our study represents a promising structural scaffold to pursue for further SARS-CoV 3CLpro inhibitor development.     

Proligands with protease‐regulated binding activity identified from cell‐displayed prodomain libraries

A general method was developed for the discovery of protease‐activated binding ligands, or proligands, from combinatorial prodomain libraries displayed on the surface of E. coli. Peptide libraries of candidate prodomains were fused with a matrix metalloprotease‐2 substrate linker to a vascular endothelial growth factor‐binding peptide and sorted using a two‐stage flow cytometry screening procedure to isolate proligands that required protease treatment for binding activity. Prodomains that imparted protease‐mediated switching activity were identified after three sorting cycles using two unique library design strategies. The best performing proligand exhibited a 100‐fold improvement in apparent binding affinity after exposure to protease. This method may prove useful for developing therapeutic and diagnostic ligands with improved systemic targeting specificity.     

In silico design of potent EGFR kinase inhibitors using combinatorial libraries

This paper is an attempt to design 4-anilinoquinazoline compounds having promising anticancer activities against epidermal growth factor (EGFR) kinase inhibition, using virtual combinatorial library approach. Partial least squares method has been applied for the development of a quantitative structure–activity relationship (QSAR) model based on training and test set approaches. The partial least squares model showed some interesting results in terms of internal and external predictability against EGFR kinase inhibition for such type of anilinoquinazoline derivatives. In virtual screening study, out of 4860 compounds in chemical library, 158 compounds were screened and finally, 10 compounds were selected as promising EGFR kinase inhibitors based on their predicted activities from the QSAR model. These derivatives were subjected to molecular docking study to investigate the mode of binding with the EGFR kinase, and the two compounds (ID 3639 and 3399) showing similar type of docking score and binding patterns with that of the existing drug molecules like erlotinib were finally reported.     

Teaching an old scaffold new tricks: monobodies constructed using alternative surfaces of the FN3 scaffold

The fibronectin type III domain (FN3) has become one of the most widely used non-antibody scaffolds for generating new binding proteins. Because of its structural homology to the immunoglobulin domain, combinatorial libraries of FN3 designed to date have primarily focused on introducing amino acid diversity into three loops that are equivalent to antibody complementarity-determining regions. Here, we report an FN3 library that utilizes alternative positions for presenting amino acid diversity. We diversified positions on a β-sheet and surface loops that together form a concave surface. The new library produced binding proteins (termed "monobodies") to multiple target proteins, generally with similar efficacy as the original, loop-focused library. The crystal structure of a monobody generated from the new library in complex with its target, the Abl SH2 domain, revealed that a concave surface of the monobody, as intended in our design, bound to a convex surface of the target with the interface area being among the largest of published structures of monobody-target complexes. This mode of interaction differs from a common binding mode for single-domain antibodies and antibody mimics in which recognition loops recognize clefts in targets. Together, this work illustrates the utilization of different surfaces of a single immunoglobulin-like scaffold to generate binding proteins with distinct characteristics.     

Docking and scoring with alternative side-chain conformations

We describe a scoring and modeling procedure for docking ligands into protein models that have either modeled or flexible side-chain conformations. Our methodical contribution comprises a procedure for generating new potentials of mean force for the ROTA scoring function which we have introduced previously for optimizing side-chain conformations with the tool IRECS. The ROTA potentials are specially trained to tolerate small-scale positional errors of atoms that are characteristic of (i) side-chain conformations that are modeled using a sparse rotamer library and (ii) ligand conformations that are generated using a docking program. We generated both rigid and flexible protein models with our side-chain prediction tool IRECS and docked ligands to proteins using the scoring function ROTA and the docking programs FlexX (for rigid side chains) and FlexE (for flexible side chains). We validated our approach on the forty screening targets of the DUD database. The validation shows that the ROTA potentials are especially well suited for estimating the binding affinity of ligands to proteins. The results also show that our procedure can compensate for the performance decrease in screening that occurs when using protein models with side chains modeled with a rotamer library instead of using X-ray structures. The average runtime per ligand of our method is 168 seconds on an Opteron V20z, which is fast enough to allow virtual screening of compound libraries for drug candidates.     

Novel high-affinity binders of human interferon gamma derived from albumin-binding domain of protein G

Recombinant ligands derived from small protein scaffolds show promise as robust research and diagnostic reagents and next generation protein therapeutics. Here, we derived high-affinity binders of human interferon gamma (hIFNγ) from the three helix bundle scaffold of the albumin-binding domain (ABD) of protein G from Streptococcus G148. Computational interaction energy mapping, solvent accessibility assessment, and in silico alanine scanning identified 11 residues from the albumin-binding surface of ABD as suitable for randomization. A corresponding combinatorial ABD scaffold library was synthesized and screened for hIFNγ binders using in vitro ribosome display selection, to yield recombinant ligands that exhibited K(d) values for hIFNγ from 0.2 to 10 nM. Molecular modeling, computational docking onto hIFNγ, and in vitro competition for hIFNγ binding revealed that four of the best ABD-derived ligands shared a common binding surface on hIFNγ, which differed from the site of human IFNγ receptor 1 binding. Thus, these hIFNγ ligands provide a proof of concept for design of novel recombinant binding proteins derived from the ABD scaffold.     

Combinatorially-generated library of 6-fluoroquinolone analogs as potential novel antitubercular agents: a chemometric and molecular modeling assessment

The virtual combinatorial chemistry approach as a methodology for generating chemical libraries of structurally-similar analogs in a virtual environment was employed for building a general mixed virtual combinatorial library with a total of 53.871 6-FQ structural analogs, introducing the real synthetic pathways of three well known 6-FQ inhibitors. The druggability properties of the generated combinatorial 6-FQs were assessed using an in-house developed drug-likeness filter integrating the Lipinski/Veber rule-sets. The compounds recognized as drug-like were used as an external set for prediction of the biological activity values using a neural-networks (NN) model based on an experimentally-determined set of active 6-FQs. Furthermore, a subset of compounds was extracted from the pool of drug-like 6-FQs, with predicted biological activity, and subsequently used in virtual screening (VS) campaign combining pharmacophore modeling and molecular docking studies. This complex scheme, a powerful combination of chemometric and molecular modeling approaches provided novel QSAR guidelines that could aid in the further lead development of 6-FQs agents.     

Parallel solution- and solid-phase synthesis of spirohydantoin derivatives as neurokinin-1 receptor ligands

The combination of the 3,5-bis(trifluoromethyl)phenyl group with a spirohydantoin motive as a central scaffold was the basis for the design of a combinatorial library targeted towards the neurokinin-1 receptor. A solution- and solid-phase procedure is described and binding affinities of representative compounds presented.     

Virtual screening against highly charged active sites: identifying substrates of alpha-beta barrel enzymes

We have developed a virtual ligand screening method designed to help assign enzymatic function for alpha-beta barrel proteins. We dock a library of approximately 19,000 known metabolites against the active site and attempt to identify the relevant substrate based on predicted relative binding free energies. These energies are computed using a physics-based energy function based on an all-atom force field (OPLS-AA) and a generalized Born implicit solvent model. We evaluate the ability of this method to identify the known substrates of several members of the enolase superfamily of enzymes, including both holo and apo structures (11 total). The active sites of these enzymes contain numerous charged groups (lysines, carboxylates, histidines, and one or more metal ions) and thus provide a challenge for most docking scoring functions, which treat electrostatics and solvation in a highly approximate manner. Using the physics-based scoring procedure, the known substrate is ranked within the top 6% of the database in all cases, and in 8 of 11 cases, it is ranked within the top 1%. Moreover, the top-ranked ligands are strongly enriched in compounds with high chemical similarity to the substrate (e.g., different substitution patterns on a similar scaffold). These results suggest that our method can be used, in conjunction with other information including genomic context and known metabolic pathways, to suggest possible substrates or classes of substrates for experimental testing. More broadly, the physics-based scoring method performs well on highly charged binding sites and is likely to be useful in inhibitor docking against polar binding sites as well. The method is fast (<1 min per ligand), due largely to an efficient minimization algorithm based on the truncated Newton method, and thus, it can be applied to thousands of ligands within a few hours on a small Linux cluster.     

Evaluation of an inverse molecular design algorithm in a model binding site

Computational molecular design is a useful tool in modern drug discovery. Virtual screening is an approach that docks and then scores individual members of compound libraries. In contrast to this forward approach, inverse approaches construct compounds from fragments, such that the computed affinity, or a combination of relevant properties, is optimized. We have recently developed a new inverse approach to drug design based on the dead-end elimination and A* algorithms employing a physical potential function. This approach has been applied to combinatorially constructed libraries of small-molecule ligands to design high-affinity HIV-1 protease inhibitors (Altman et al., J Am Chem Soc 2008;130:6099-6013). Here we have evaluated the new method using the well-studied W191G mutant of cytochrome c peroxidase. This mutant possesses a charged binding pocket and has been used to evaluate other design approaches. The results show that overall the new inverse approach does an excellent job of separating binders from nonbinders. For a few individual cases, scoring inaccuracies led to false positives. The majority of these involve erroneous solvation energy estimation for charged amines, anilinium ions, and phenols, which has been observed previously for a variety of scoring algorithms. Interestingly, although inverse approaches are generally expected to identify some but not all binders in a library, due to limited conformational searching, these results show excellent coverage of the known binders while still showing strong discrimination of the nonbinders.     

Peptidic 1-cyanopyrrolidines: synthesis and SAR of a series of potent,selective cathepsin inhibitors

1-Cyanopyrrolidines have previously been reported to inhibit cysteinyl cathepsins (Falgueyret, J.-P. et al., J. Med. Chem. 2001, 44, 94). In order to optimize binding interactions for a given cathepsin and simultaneously reduce interactions with the other closely related enzymes, small peptidic substituents were introduced to the 1-cyanopyrrolidine scaffold, either at the 2-position starting with proline or at the 3-position of aminopyrrolidines. The resulting novel compounds proved to be micromolar inhibitors of cathepsin B (Cat B) but nanomolar to picomolar inhibitors of cathepsins K, L, and S (Cat K, Cat L, Cat S). Several of the compounds were >20-fold selective versus the other three cathepsins. SAR trends were observed, most notably the remarkable potency of Cat L inhibitors based on the 1-cyano-D-proline scaffold. The selectivity of one such compound, the 94 picomolar Cat L inhibitor 12, was demonstrated at higher concentrations in DLD-1 cells. Although none of the compounds in the proline series that was tested proved to be submicromolar in the in vitro bone resorption assay, two Cat K inhibitors in the 3-substituted pyrrolidine series, 24 and 25 were relatively potent in that assay.     

Optimization of combinatorial mutagenesis

Protein engineering by combinatorial site-directed mutagenesis evaluates a portion of the sequence space near a target protein, seeking variants with improved properties (e.g., stability, activity, immunogenicity). In order to improve the hit-rate of beneficial variants in such mutagenesis libraries, we develop methods to select optimal positions and corresponding sets of the mutations that will be used, in all combinations, in constructing a library for experimental evaluation. Our approach, OCoM (Optimization of Combinatorial Mutagenesis), encompasses both degenerate oligonucleotides and specified point mutations, and can be directed accordingly by requirements of experimental cost and library size. It evaluates the quality of the resulting library by one- and two-body sequence potentials, averaged over the variants. To ensure that it is not simply recapitulating extant sequences, it balances the quality of a library with an explicit evaluation of the novelty of its members. We show that, despite dealing with a combinatorial set of variants, in our approach the resulting library optimization problem is actually isomorphic to single-variant optimization. By the same token, this means that the two-body sequence potential results in an NP-hard optimization problem. We present an efficient dynamic programming algorithm for the one-body case and a practically-efficient integer programming approach for the general two-body case. We demonstrate the effectiveness of our approach in designing libraries for three different case study proteins targeted by previous combinatorial libraries--a green fluorescent protein, a cytochrome P450, and a beta lactamase. We found that OCoM worked quite efficiently in practice, requiring only 1 hour even for the massive design problem of selecting 18 mutations to generate 10? variants of a 443-residue P450. We demonstrate the general ability of OCoM in enabling the protein engineer to explore and evaluate trade-offs between quality and novelty as well as library construction technique, and identify optimal libraries for experimental evaluation.     

A new method for the screening of solid-phase combinatorial libraries for affinity chromatography

A new methodology for the rapid assessment of affinity ligands synthesized by combinatorial solid-phase chemistry is reported. This screening strategy utilizes the target protein conjugated to FITC, and represents an almost universal technique for the preliminary screening of solid-phase combinatorial libraries. The assessment of a triazine-scaffolded solid-phase combinatorial library of ligands, designed to bind to human IgG, was performed with FITC-human IgG, and the results compared with those obtained by conventional affinity chromatographic screening assays. The effect of different molar conjugation ratios of FITC-IgG (F/P) was evaluated. Independently of the F/P ratio, no false negative results were observed, although lower F/P ratios diminished non-specific interactions and the number of false positives. The nature of the substituents on the triazine scaffold was not related to the number of false positive IgG-binding ligands. The reproducibility of the FITC technique, using FITC-human IgG conjugates with low F/P ratio (F/P=2), was also evaluated. The FITC-based technique proved to be efficient and accurate in the identification of strongly binding ligands (binding >50% of loaded protein, by standard affinity chromatographic assays), and is envisaged as a versatile and cost-effective method to screen other systems, and evaluate several binding/elution conditions at small-scale, prior to scale-up to standard affinity chromatography.     

Expression,purification, and characterization of proteins from high-quality combinatorial libraries of the mammalian calmodulin central linker

Combinatorial libraries offer an attractive approach towards exploring protein sequence, structure and function. Although several strategies introduce sequence diversity, the likelihood of identifying proteins with novel functions is increased when the library of genes encodes for folded and soluble structures. Here we present the first application of the binary patterning approach of combinatorial protein library design to the unique central linker region of the highly-conserved protein, calmodulin (CaM). We show that this high-quality approach translates very well to the CaM protein scaffold: All library members over-express and are functionally diverse, having a range of conformations in the presence and absence of calcium as determined by circular dichroism spectroscopy. Collectively, these data support that the binary patterning approach, when applied to the highly-conserved protein fold, can yield large collections of folded, soluble and highly-expressible proteins.     

High-resolution solution structure of two members of a conformationally homogeneous combinatorial peptide library based on the classical zinc-finger motif

Summary We describe the high-resolution structure by NMR of two peptides that belong to a combinatorial library based on the zinc-finger motif. The library represents, to the best of our knowledge, the first example of a conformationally homogeneous peptide library and was obtained by introducing random residues in five positions of the -helical portion of a 26-residue consensus peptide (CP1) belonging to the Cys²-Hys² zinc-finger family. The result was shown to be a highly homogeneous -helical library (Bianchi et al., 1995). The structures of the parent compound (CP1) and of a representative member (CP1m) that was selected by screening the library with a monoclonal antibody are compared in detail as an example of the very high stability of the zinc-finger scaffold upon sequence variability. The two peptides exhibit an extremely high degree of structural similarity. The use of this type of conformationally constrained combinatorial library might represent a step forward in the design of peptidomimetics, as it considerably accelerates the process of the identification of the spatial relationship among the pharmacophoric groups.Abbreviations t-Bu tert-butyloxycarbonyl - Fmoc 9-fluorenylmethoxycarbonyl     

分子对接在基于结构药物设计中的应用

分子对接是研究分子间（如配体如受体）相互作用,并预测其结合模式和亲合力的一种理论模拟方法。近年来,分子对接方法已成为计算机辅助药物研究领域的一项重要技术,在数据库搜寻,组合库设计及蛋白作用研究方面得到了广泛发展。     

Green fluorescent protein as a scaffold for intracellular presentation of peptides.

Peptide aptamers provide probes for biological processes and adjuncts for development of novel pharmaceutical molecules. Such aptamers are analogous to compounds derived from combinatorial chemical libraries which have specific binding or inhibitory activities. Much as it is generally difficult to determine the composition of combinatorial chemical libraries in a quantitative manner, determining the quality and characteristics of peptide libraries displayed in vivo is problematical. To help address these issues we have adapted green fluorescent protein (GFP) as a scaffold for display of conformationally constrained peptides. The GFP-peptide libraries permit analysis of library diversity and expression levels in cells and allow enrichment of the libraries for sequences with predetermined characteristics, such as high expression of correctly folded protein, by selection for high fluorescence.     

Organization of a human UDP-GalNAc:polypeptide, N-acetylgalactosaminyltransferase gene and a related processed pseudogene

We have previously characterized a cDNA that encodes a fulllength human UDP-GalNAc:polypeptide, N-acetyl galactosaminyltransferase(GalNAc-transferase) (J.A. Meurer et al., J. Biochem., 118,568–574, 1995). The present report describes the characterizationof the corresponding human GalNAc-transferase gene and a relatedpseudogene. Two human genomic libraries,     

The inverse problems for some topological indices in combinatorial chemistry.

In the original paper, Goldman et al. (2000) launched the study of the inverse problems in combinatorial chemistry, which is closely related to the design of combinatorial libraries for drug discovery. Following their ideas, we investigate four other topological indices, i.e., the sigma-index, the c-index, the Z-index, and the M(1)-index, with a special emphasis on the sigma-index. Like the Wiener index, these four indices are very popular in combinatorial chemistry and reflect many chemical and physical properties. We give algorithmic and analytical solutions for the inverse problems of the four indices. We also show that the SUBTREEVALUE reconstruction problem for the sigma-index is NP-hard.