Inhibition of integrin αvβ6 is a promising approach to the treatment of fibrotic disease such as idiopathic pulmonary fibrosis. Screening a small library combining head groups that stabilize the bent-closed conformation of integrin αIIbβ3 with αv integrin binding motifs resulted in the identification of hit compounds that bind the bent-closed conformation of αvβ6. Crystal structures of these compounds bound to αvβ6 and related integrins revealed opportunities to increase potency and selectivity, and these efforts were accelerated using accurate free energy perturbation (FEP+) calculations. Optimization of PK parameters including permeability, bioavailability, clearance, and half-life resulted in the discovery of development candidate MORF-627, a highly selective inhibitor of αvβ6 that stabilizes the bent-closed conformation and has good oral PK. Unfortunately, the compound showed toxicity in a 28-day NHP safety study, precluding further development. Nevertheless, MORF-627 is a useful tool compound for studying the biology of integrin αvβ6.

D-Serine is a co-agonist of the N-methyl D-aspartate (NMDA) receptor, a key excitatory neurotransmitter receptor. In the brain, D-Serine is synthesized from its L-isomer by serine racemase and is metabolized by the D-amino acid oxidase (DAO, DAAO), a flavoenzyme that catalyzes the oxidative degradation of D-amino acids including D-serine to the corresponding α-keto acids. Many studies have linked decreased D-serine concentration and/or increased DAO expression and enzyme activity to NMDA dysfunction and schizophrenia. Thus, many companies have explored the possibility of employing DAO inhibitors for the treatment of schizophrenia and other indications. Powered by the Schrödinger computational modeling platform, we initiated a research program to identify novel DAO inhibitors with best-in-class properties. The program execution leveraged an hDAO FEP+ model to prospectively predict compound hDAO inhibitory potency and prioritize design ideas from both human design and computer enumeration by our AutoDesigner algorithm. A novel class of DAO inhibitors with desirable pharmacokinetic and brain penetration properties was discovered from this effort. In an in vivo mouse PK/PD model, tool compound 37 demonstrated modulation of D-serine concentrations in the plasma and brain through inhibition of DAO function. Continued SAR work has led to significant potency improvement in both DAO biochemical and cell assays. Our modeling technology on this program has not only enhanced the efficiency of medicinal chemistry execution, it has also helped to identify a previously unexplored subpocket for further SAR development.

There has been a remarkable increase in the number of biologics, especially monoclonal antibodies, in the market over the last decade. In addition to attaining the desired binding to their targets, a crucial aspect is the 'developability' of these drugs, which includes several desirable properties such as high solubility, low viscosity and aggregation, physico-chemical stability, low immunogenicity and low poly-specificity. The lack of any of these desirable properties can lead to significant hurdles in advancing them to the clinic and are often discovered only during late stages of drug development. Hence, in silico methods for early detection of these properties, particularly the ones that affect aggregation and solubility in the earlier stages can be highly beneficial. We have developed a computational framework based on a large and diverse set of protein specific descriptors that is ideal for making liability predictions using a QSPR (quantitative structure-property relationship) approach. This set offers a high degree of feature diversity that may coarsely be classified based on (1) sequence (2) structure and (3) surface patches. We assess the sensitivity and applicability of these descriptors in four dedicated case studies that are believed to be representative of biophysical characterizations commonly employed during the development process of a biologics drug candidate. In addition to data sets obtained from public sources, we have validated the descriptors on novel experimental data sets in order to address antibody developability and to generate prospective predictions on Adnectins. The results show that the descriptors are well suited to assist in the improvement of protein properties of systems that exhibit poor solubility or aggregation.

We present a reliable and accurate solution to the induced fit docking problem for protein-ligand binding by combining ligand-based pharmacophore docking (Phase), rigid receptor docking (Glide), and protein structure prediction (Prime) with explicit solvent molecular dynamics simulations. We provide an in-depth description of our novel methodology and present results for 41 targets consisting of 415 cross-docking cases divided amongst a training and test set. For both the training and test-set, we compute binding modes with a ligand-heavy atom RMSD to within 2.5 Å or better in over 90% of cross-docking cases compared to less than 70% of cross-docking cases using our previously published induced-fit docking algorithm and less than 41% using rigid receptor docking. Applications of the predicted ligand-receptor structure in free energy perturbation calculations is demonstrated for both public data and in active drug discovery projects, both retrospectively and prospectively.

Proteins and ligands sample a conformational ensemble that governs molecular recognition, activity, and dissociation. In structure-based drug design, access to this conformational ensemble is critical to understand the balance between entropy and enthalpy in lead optimization. However, ligand conformational heterogeneity is currently severely underreported in crystal structures in the Protein Data Bank, owing in part to a lack of automated and unbiased procedures to model an ensemble of protein-ligand states into X-ray data. Here, we designed a computational method, qFit-ligand, to automatically resolve conformationally averaged ligand heterogeneity in crystal structures, and applied it to a large set of protein receptor-ligand complexes. In an analysis of the cancer related BRD4 domain, we found that up to 29% of protein crystal structures bound with drug-like molecules present evidence of unmodeled, averaged, relatively isoenergetic conformations in ligand-receptor interactions. In many retrospective cases, these alternate conformations were adventitiously exploited to guide compound design, resulting in improved potency or selectivity. Combining qFit-ligand with high-throughput screening or multitemperature crystallography could therefore augment the structure-based drug design toolbox.

Protein aggregation is a phenomenon that has attracted considerable attention within the pharmaceutical industry from both a developability standpoint (to ensure stability of protein formulations) and from a research perspective for neurodegenerative diseases. Experimental identification of aggregation behavior in proteins can be expensive; and hence, the development of accurate computational approaches is crucial. The existing methods for predicting protein aggregation rely mostly on the primary sequence and are typically trained on amyloid-like proteins. However, the training bias toward beta amyloid peptides may worsen prediction accuracy of such models when applied to larger protein systems. Here, we present a novel algorithm to identify aggregation-prone regions in proteins termed "AggScore" that is based entirely on three-dimensional structure input. The method uses the distribution of hydrophobic and electrostatic patches on the surface of the protein, factoring in the intensity and relative orientation of the respective surface patches into an aggregation propensity function that has been trained on a benchmark set of 31 adnectin proteins. AggScore can accurately identify aggregation-prone regions in several well-studied proteins and also reliably predict changes in aggregation behavior upon residue mutation. The method is agnostic to an amyloid-specific aggregation context and thus may be applied to globular proteins, small peptides and antibodies.

A novel method for exploring macrocycle conformational space, Prime macrocycle conformational sampling (Prime-MCS), is introduced and evaluated in the context of other available algorithms (Molecular Dynamics, LowModeMD in MOE, and MacroModel Baseline Search). The algorithms were benchmarked on a data set of 208 macrocycles which was curated for diversity from the Cambridge Structural Database, the Protein Data Bank, and the Biologically Interesting Molecule Reference Dictionary. The algorithms were evaluated in terms of accuracy (ability to reproduce the crystal structure), diversity (coverage of conformational space), and computational speed. Prime-MCS most reliably reproduced crystallographic structures for RMSD thresholds >1.0 Å, most often produced the most diverse conformational ensemble, and was most often the fastest algorithm. Detailed analysis and examination of both typical and outlier cases were performed to reveal characteristics, shortcomings, expected performance, and complementarity of the methods.

Although many popular docking programs include a facility to account for covalent ligands, large-scale systematic docking validation studies of covalent inhibitors have been sparse. In this paper, we present the development and validation of a novel approach for docking and scoring covalent inhibitors, which consists of conventional noncovalent docking, heuristic formation of the covalent attachment point, and structural refinement of the protein-ligand complex. This approach combines the strengths of the docking program Glide and the protein structure modeling program Prime and does not require any parameter fitting for the study of additional covalent reaction types. We first test this method by predicting the native binding geometry of 38 covalently bound complexes. The average RMSD of the predicted poses is 1.52 Å, and 76% of test set inhibitors have an RMSD of less than 2.0 Å. In addition, the apparent affinity score constructed herein is tested on a virtual screening study and the characterization of the SAR properties of two different series of congeneric compounds with satisfactory success.

ABSTRACTWe present the blinded prediction results in the Second Antibody Modeling Assessment (AMA‐II) using a fully automatic antibody structure prediction method implemented in the programs BioLuminate and Prime. We have developed a novel knowledge based approach to model the CDR loops, using a combination of sequence similarity, geometry matching, and the clustering of database structures. The homology models are further optimized with a physics‐based energy function (VSGB2.0), which improves the model quality significantly. H3 loop modeling remains the most challenging task. Our ab initio loop prediction performs well for the H3 loop in the crystal structure context, and allows improved results when refining the H3 loops in the context of homology models. For the 10 human and mouse derived antibodies in this assessment, the average RMSDs for the homology model Fv and framework regions are 1.19 Å and 0.74 Å, respectively. The average RMSDs for five non‐H3 CDR loops range from 0.61 Å to 1.05 Å, and the H3 loop average RMSD is 2.91 Å using our knowledge‐based loop prediction approach. The ab initio H3 loop predictions yield an average RMSD of 1.28 Å when performed in the context of the crystal structure and 2.67 Å in the context of the homology modeled structure. Notably, our method for predicting the H3 loop in the crystal structure environment ranked first among the seven participating groups in AMA‐II, and our method made the best prediction among all participants for seven of the ten targets. Proteins 2014; 82:1646–1655. © 2014 Wiley Periodicals, Inc.

Structure-based virtual screening plays an important role in drug discovery and complements other screening approaches. In general, protein crystal structures are prepared prior to docking in order to add hydrogen atoms, optimize hydrogen bonds, remove atomic clashes, and perform other operations that are not part of the x-ray crystal structure refinement process. In addition, ligands must be prepared to create 3-dimensional geometries, assign proper bond orders, and generate accessible tautomer and ionization states prior to virtual screening. While the prerequisite for proper system preparation is generally accepted in the field, an extensive study of the preparation steps and their effect on virtual screening enrichments has not been performed. In this work, we systematically explore each of the steps involved in preparing a system for virtual screening. We first explore a large number of parameters using the Glide validation set of 36 crystal structures and 1,000 decoys. We then apply a subset of protocols to the DUD database. We show that database enrichment is improved with proper preparation and that neglecting certain steps of the preparation process produces a systematic degradation in enrichments, which can be large for some targets. We provide examples illustrating the structural changes introduced by the preparation that impact database enrichment. While the work presented here was performed with the Protein Preparation Wizard and Glide, the insights and guidance are expected to be generalizable to structure-based virtual screening with other docking methods.

Antibodies have the capability of binding a wide range of antigens due to the diversity of the six loops constituting the complementarity determining region (CDR). Among the six loops, the H3 loop is the most diverse in structure, length, and sequence identity. Prediction of the three‐dimensional structures of antibodies, especially the CDR loops, is an important step in the computational design and engineering of novel antibodies for improved affinity and specificity. Although it has been demonstrated that the conformation of the five non‐H3 loops can be accurately predicted by comparing their sequences against databases of canonical loop conformations, no such connection has been established for H3 loops. In this work, we present the results for ab initio structure prediction of the H3 loop using conformational sampling and energy calculations with the program Prime on a dataset of 53 loops ranging in length from 4 to 22 residues. When the prediction is performed in the crystal environment and including symmetry mates, the median backbone root mean square deviation (RMSD) is 0.5 Å to the crystal structure, with 91% of cases having an RMSD of less than 2.0 Å. When the prediction is performed in a noncrystallographic environment, where the scaffold is constructed by swapping the H3 loops between homologous antibodies, 70% of cases have an RMSD below 2.0 Å. These results show promise for ab initio loop predictions applied to modeling of antibodies. © 2012 Wiley Periodicals, Inc.

We present a novel protein-ligand docking method that accurately accounts for both ligand and receptor flexibility by iteratively combining rigid receptor docking (Glide) with protein structure prediction (Prime) techniques. While traditional rigid-receptor docking methods are useful when the receptor structure does not change substantially upon ligand binding, success is limited when the protein must be "induced" into the correct binding conformation for a given ligand. We provide an in-depth description of our novel methodology and present results for 21 pharmaceutically relevant examples. Traditional rigid-receptor docking for these 21 cases yields an average RMSD of 5.5 A. The average ligand RMSD for docking to a flexible receptor for the 21 pairs is 1.4 A; the RMSD is < or =1.8 A for 18 of the cases. For the three cases with RMSDs greater than 1.8 A, the core of the ligand is properly docked and all key protein/ligand interactions are captured.