Exploiting high-energy hydration sites for the discovery of potent peptide aldehyde inhibitors of the SARS-CoV-2 main protease with cellular antiviral activity.

Small-molecule antivirals that prevent the replication of the SARS-CoV-2 virus by blocking the enzymatic activity of its main protease (Mpro) are and will be a tenet of pandemic preparedness. However, the peptidic nature of such compounds often precludes the design of compounds within favorable physical property ranges, limiting cellular activity. Here we describe the discovery of peptide aldehyde Mpro inhibitors with potent enzymatic and cellular antiviral activity. This structure-activity relationship (SAR) exploration was guided by the use of calculated hydration site thermodynamic maps (WaterMap) to drive potency via displacement of waters from high-energy sites. Thousands of diverse compounds were designed to target these high-energy hydration sites and then prioritized for synthesis by physics- and structure-based Free-Energy Perturbation (FEP+) simulations, which accurately predicted biochemical potencies. This approach ultimately led to the rapid discovery of lead compounds with unique SAR that exhibited potent enzymatic and cellular activity with excellent pan-coronavirus coverage.

Structure-based identification and neutralization mechanism of tyrosine sulfate mimetics that inhibit HIV-1 entry.

Tyrosine sulfate-mediated interactions play an important role in HIV-1 entry. After engaging the CD4 receptor at the cell surface, the HIV-1 gp120 glycoprotein binds to the CCR5 co-receptor via an interaction that requires two tyrosine sulfates, at positions 10 and 14 in the CCR5-N terminus. Building on previous structure determinations of this interaction, here we report the targeting of these tyrosine sulfate binding sites for drug design through in silico screening of small molecule libraries, identification of lead compounds, and characterization of biological activity. A class of tyrosine sulfate-mimicking small molecules containing a "phenyl sulfonate-linker-aromatic" motif was identified that specifically inhibited binding of gp120 to the CCR5-N terminus as well as to sulfated antibodies that recognize the co-receptor binding region on gp120. The most potent of these compounds bound gp120 with low micromolar affinity and its CD4-induced conformation with K(D)'s as tight as ∼50 nM. Neutralization experiments suggested the targeted site to be conformationally inaccessible prior to CD4 engagement. Primary HIV-1 isolates were weakly neutralized, preincubation with soluble CD4 enhanced neutralization, and engineered isolates with increased dependence on the N terminus of CCR5 or with reduced conformational barriers were neutralized with IC(50) values as low as ∼1 µM. These results reveal the potential of targeting the tyrosine sulfate interactions of HIV-1 and provide insight into how mechanistic barriers, evolved by HIV-1 to evade antibody recognition, also restrict small-molecule-mediated neutralization.

Epik: a software program for pK( a ) prediction and protonation state generation for drug-like molecules.

Epik is a computer program for predicting pK(a) values for drug-like molecules. Epik can use this capability in combination with technology for tautomerization to adjust the protonation state of small drug-like molecules to automatically generate one or more of the most probable forms for use in further molecular modeling studies. Many medicinal chemicals can exchange protons with their environment, resulting in various ionization and tautomeric states, collectively known as protonation states. The protonation state of a drug can affect its solubility and membrane permeability. In modeling, the protonation state of a ligand will also affect which conformations are predicted for the molecule, as well as predictions for binding modes and ligand affinities based upon protein-ligand interactions. Despite the importance of the protonation state, many databases of candidate molecules used in drug development do not store reliable information on the most probable protonation states. Epik is sufficiently rapid and accurate to process large databases of drug-like molecules to provide this information. Several new technologies are employed. Extensions to the well-established Hammett and Taft approaches are used for pK(a) prediction, namely, mesomer standardization, charge cancellation, and charge spreading to make the predicted results reflect the nature of the molecule itself rather just for the particular Lewis structure used on input. In addition, a new iterative technology for generating, ranking and culling the generated protonation states is employed.

Extra precision glide: docking and scoring incorporating a model of hydrophobic enclosure for protein-ligand complexes.

A novel scoring function to estimate protein-ligand binding affinities has been developed and implemented as the Glide 4.0 XP scoring function and docking protocol. In addition to unique water desolvation energy terms, protein-ligand structural motifs leading to enhanced binding affinity are included: (1) hydrophobic enclosure where groups of lipophilic ligand atoms are enclosed on opposite faces by lipophilic protein atoms, (2) neutral-neutral single or correlated hydrogen bonds in a hydrophobically enclosed environment, and (3) five categories of charged-charged hydrogen bonds. The XP scoring function and docking protocol have been developed to reproduce experimental binding affinities for a set of 198 complexes (RMSDs of 2.26 and 1.73 kcal/mol over all and well-docked ligands, respectively) and to yield quality enrichments for a set of fifteen screens of pharmaceutical importance. Enrichment results demonstrate the importance of the novel XP molecular recognition and water scoring in separating active and inactive ligands and avoiding false positives.

Glide: a new approach for rapid, accurate docking and scoring. 2. Enrichment factors in database screening.

Glide's ability to identify active compounds in a database screen is characterized by applying Glide to a diverse set of nine protein receptors. In many cases, two, or even three, protein sites are employed to probe the sensitivity of the results to the site geometry. To make the database screens as realistic as possible, the screens use sets of "druglike" decoy ligands that have been selected to be representative of what we believe is likely to be found in the compound collection of a pharmaceutical or biotechnology company. Results are presented for releases 1.8, 2.0, and 2.5 of Glide. The comparisons show that average measures for both "early" and "global" enrichment for Glide 2.5 are 3 times higher than for Glide 1.8 and more than 2 times higher than for Glide 2.0 because of better results for the least well-handled screens. This improvement in enrichment stems largely from the better balance of the more widely parametrized GlideScore 2.5 function and the inclusion of terms that penalize ligand-protein interactions that violate established principles of physical chemistry, particularly as it concerns the exposure to solvent of charged protein and ligand groups. Comparisons to results for the thymidine kinase and estrogen receptors published by Rognan and co-workers (J. Med. Chem. 2000, 43, 4759-4767) show that Glide 2.5 performs better than GOLD 1.1, FlexX 1.8, or DOCK 4.01.

Design and synthesis of dipeptide nitriles as reversible and potent Cathepsin S inhibitors.

The specificity of the immune response relies on processing of foreign proteins and presentation of antigenic peptides at the cell surface. Inhibition of antigen presentation, and the subsequent activation of T-cells, should, in theory, modulate the immune response. The cysteine protease Cathepsin S performs a fundamental step in antigen presentation and therefore represents an attractive target for inhibition. Herein, we report a series of potent and reversible Cathepsin S inhibitors based on dipeptide nitriles. These inhibitors show nanomolar inhibition of the target enzyme as well as cellular potency in a human B cell line. The first X-ray crystal structure of a reversible inhibitor cocrystallized with Cathepsin S is also reported.