Design, synthesis, and evaluation of nonpeptidic inhibitors of human rhinovirus 3C protease.

The design, synthesis, and biological evaluation of reversible, nonpeptidic inhibitors of human rhinovirus (HRV) 3C protease (3CP) are reported. A novel series of 2,3-dioxindoles (isatins) were designed that utilized a combination of protein structure-based drug design, molecular modeling, and structure-activity relationship (SAR). The C-2 carbonyl of isatin was envisioned to react in the active site of HRV 3CP with the cysteine responsible for catalytic proteolysis, thus forming a stabilized transition state mimic. Molecular-modeling experiments using the apo crystal structure of human rhinovirus-serotype 14 (HRV-14) 3CP and a peptide substrate model allowed us to design recognition features into the P1 and P2 subsites, respectively, from the 5- and 1-positions of isatin. Attempts to optimize recognition properties in the P1 subsite using SAR at the 5-position were performed. In addition, a series of ab initio calculations were carried out on several 5-substituted isatins to investigate the stability of sulfide adducts at C-3. The inhibitors were prepared by general synthetic methods, starting with commercially available 5-substituted isatins in nearly every case. All compounds were tested for inhibition of purified HRV-14 3CP. Compounds 8, 14, and 19 were found to have excellent selectivity for HRV-14 3CP compared to other proteolytic enzymes, including chymotrypsin and cathepsin B. Selected compounds were assayed for antiviral activity against HRV-14-infected HI-HeLa cells. A 2.8 A cocrystal structure of derivative 19 covalently bound to human rhinovirus-serotype 2 (HRV-2) 3CP was solved and revealed that the isatin was situated in essentially the same conformation as modeled.

Substituted benzamide inhibitors of human rhinovirus 3C protease: structure-based design, synthesis, and biological evaluation.

A series of nonpeptide benzamide-containing inhibitors of human rhinovirus (HRV) 3C protease was identified using structure-based design. The design, synthesis, and biological evaluation of these inhibitors are reported. A Michael acceptor was combined with a benzamide core mimicking the P1 recognition element of the natural 3CP substrate. alpha,beta-Unsaturated cinnamate esters irreversibly inhibited the 3CP and displayed antiviral activity (EC(50) 0.60 microM, HRV-16 infected H1-HeLa cells). On the basis of cocrystal structure information, a library of substituted benzamide derivatives was prepared using parallel synthesis on solid support. A 1.9 A cocrystal structure of a benzamide inhibitor in complex with the 3CP revealed a binding mode similar to that initially modeled wherein covalent attachment of the nucleophilic cysteine residue is observed. Unsaturated ketones displayed potent reversible inhibition but were inactive in the cellular antiviral assay and were found to react with nucleophilic thiols such as DTT.

Structure-based design, synthesis, and biological evaluation of irreversible human rhinovirus 3C protease inhibitors. 4. Incorporation of P1 lactam moieties as L-glutamine replacements.

The structure-based design, chemical synthesis, and biological evaluation of various human rhinovirus (HRV) 3C protease (3CP) inhibitors which incorporate P1 lactam moieties in lieu of an L-glutamine residue are described. These compounds are comprised of a tripeptidyl or peptidomimetic binding determinant and an ethyl propenoate Michael acceptor moiety which forms an irreversible covalent adduct with the active site cysteine residue of the 3C enzyme. The P1-lactam-containing inhibitors display significantly increased 3CP inhibition activity along with improved antirhinoviral properties relative to corresponding L-glutamine-derived molecules. In addition, several lactam-containing compounds exhibit excellent selectivity for HRV 3CP over several other serine and cysteine proteases and are not appreciably degraded by a variety of biological agents. One of the most potent inhibitors (AG7088, mean antirhinoviral EC90 approximately 0.10 microM, n = 46 serotypes) is shown to warrant additional preclinical development to explore its potential for use as an antirhinoviral agent.

Tripeptide aldehyde inhibitors of human rhinovirus 3C protease: design, synthesis, biological evaluation, and cocrystal structure solution of P1 glutamine isosteric replacements.

The investigation of tripeptide aldehydes as reversible covalent inhibitors of human rhinovirus (HRV) 3C protease (3CP) is reported. Molecular models based on the apo crystal structure of HRV-14 3CP and other trypsin-like serine proteases were constructed to approximate the binding of peptide substrates, generate transition state models of P1-P1' amide cleavage, and propose novel tripeptide aldehydes. Glutaminal derivatives have limitations since they exist predominantly in the cyclic hemiaminal form. Therefore, several isosteric replacements for the P1 carboxamide side chain were designed and incorporated into the tripeptide aldehydes. These compounds were found to be potent inhibitors of purified HRV-14 3CP with Kis ranging from 0.005 to 0.64 microM. Several have low micromolar antiviral activity when tested against HRV-14-infected H1-HeLa cells. The N-acetyl derivative 3 was also shown to be active against HRV serotypes 2, 16, and 89. High-resolution cocrystal structures of HRV-2 3CP, covalently bound to compounds 3, 15, and 16, were solved. These cocrystal structures were analyzed and compared with our original HRV-14 3CP-substrate and inhibitor models.

Backbone entropy of loops as a measure of their flexibility: application to a ras protein simulated by molecular dynamics.

The flexibility of surface loops plays an important role in protein-protein and protein-peptide recognition; it is commonly studied by Molecular Dynamics or Monte Carlo stimulations. We propose to measure the relative backbone flexibility of loops by the difference in their backbone conformational entropies, which are calculated here with the local states (LS) method of Meirovitch. Thus, one can compare the entropies of loops of the same protein or, under certain simulation conditions, of different proteins. These loops should be equal in size but can differ in their sequence of amino acids residues. This methodology is applied successfully to three segments of 10 residues of a Ras protein simulated by the stochastic boundary molecular dynamics procedure. For the first time estimates of backbone entropy differences are obtained, and their correlation with B factors is pointed out; for example, the segments which consist of residues 60-65 and 112-117 have average B factors of 67 and 18 A2, respectively, and entropy difference T delta S = 5.4 +/- 0.1 kcal/mol at T = 300 K. In a large number of recent publications the entropy due to the fast motions (on the ps-ns time scale) of N-H and C-H vectors has been obtained from their order parameter, measured in nuclear magnetic resonance spin relaxation experiments. This enables one to estimate differences in the entropy of protein segments due to folding-unfolding transitions, for example. However, the vectors are assumed to be independent, and the effect of the neglected correlations is unknown; our method is expected to become an important tool for assessing this approximation. The present calculations, obtained with the LS method, suggest that the errors involved in experimental entropy differences might not be large; however, this should be verified in each case. Potential applications of entropy calculations to rational drug design are discussed.

Structure-assisted design of mechanism-based irreversible inhibitors of human rhinovirus 3C protease with potent antiviral activity against multiple rhinovirus serotypes.

Human rhinoviruses, the most important etiologic agents of the common cold, are messenger-active single-stranded monocistronic RNA viruses that have evolved a highly complex cascade of proteolytic processing events to control viral gene expression and replication. Most maturation cleavages within the precursor polyprotein are mediated by rhinovirus 3C protease (or its immediate precursor, 3CD), a cysteine protease with a trypsin-like polypeptide fold. High-resolution crystal structures of the enzyme from three viral serotypes have been used for the design and elaboration of 3C protease inhibitors representing different structural and chemical classes. Inhibitors having alpha,beta-unsaturated carbonyl groups combined with peptidyl-binding elements specific for 3C protease undergo a Michael reaction mediated by nucleophilic addition of the enzyme's catalytic Cys-147, resulting in covalent-bond formation and irreversible inactivation of the viral protease. Direct inhibition of 3C proteolytic activity in virally infected cells treated with these compounds can be inferred from dose-dependent accumulations of viral precursor polyproteins as determined by SDS/PAGE analysis of radiolabeled proteins. Cocrystal-structure-assisted optimization of 3C-protease-directed Michael acceptors has yielded molecules having extremely rapid in vitro inactivation of the viral protease, potent antiviral activity against multiple rhinovirus serotypes and low cellular toxicity. Recently, one compound in this series, AG7088, has entered clinical trials.

Interaction of a novel GDP exchange inhibitor with the Ras protein.

Mutated, tumorigenic Ras is present in a variety of human tumors. Compounds that inhibit tumorigenic Ras function may be useful in the treatment of Ras-related tumors. The interaction of a novel GDP exchange inhibitor (SCH-54292) with the Ras-GDP protein was studied by NMR spectroscopy. The binding of the inhibitor to the Ras protein was enhanced at low Mg2+ concentrations, which enabled the preparation of a stable complex for NMR study. To understand the enhanced inhibitor binding and the increased GDP dissociation rates of the Ras protein, the conformational changes of the Ras protein at low Mg2+ concentrations was investigated using two-dimensional 1H-15N HSQC experiments. The Ras protein existed in two conformations in slow exchange on the NMR time scale under such conditions. The conformational changes mainly occurred in the GDP binding pocket, in the switch I and the switch II regions, and were reversible. The Ras protein resumed its regular conformation after an excess amount of Mg2+ was added. A model of the inhibitor in complex with the Ras-GDP protein was derived from intra- and intermolecular NOE distance constraints, and revealed that the inhibitor bound to the critical switch II region of the Ras protein.