基于“底物包膜”假说筛选新型HIV-1蛋白酶抑制剂

基于“底物包膜”假说, 以现有HIV-1蛋白酶抑制剂Darunavir为模板构建药效团模型并对中药化学数据库进行搜索; 采用分子对接方法进一步考察化合物与HIV-1蛋白酶结合情况及其与“底物包膜”符合程度, 优先选出两个化合物Annomonicin和去乙酰蟾蜍它灵; 应用分子动力学方法对这两个化合物进行动力学模拟, 观察它们与蛋白酶结合的复合物在动力学过程中的稳定性并计算其结合自由能, 综合评价筛选结果, 最终确定化合物Annomonicin具有更潜在的深入研究价值.     

Energetic basis for drug resistance of HIV-1 protease mutants against amprenavir

Amprenavir (APV) is a high affinity (0.15 nM) HIV-1 protease (PR) inhibitor. However, the affinities of the drug resistant protease variants V32I, I50V, I54V, I54M, I84V and L90M to amprenavir are decreased 3 to 30-fold compared to the wild-type. In this work, the popular molecular mechanics Poisson-Boltzmann surface area method has been used to investigate the effectiveness of amprenavir against the wild-type and these mutated protease variants. Our results reveal that the protonation state of Asp25/Asp25′ strongly affects the dynamics, the overall affinity and the interactions of the inhibitor with individual residues. We emphasize that, in contrast to what is often assumed, the protonation state may not be inferred from the affinities but requires pK_a calculations. At neutral pH, Asp25 and Asp25′ are ionized or protonated, respectively, as suggested from pK_a calculations. This protonation state was thus mainly considered in our study. Mutation induced changes in binding affinities are in agreement with the experimental findings. The decomposition of the binding free energy reveals the mechanisms underlying binding and drug resistance. Drug resistance arises from an increase in the energetic contribution from the van der Waals interactions between APV and PR (V32I, I50V, and I84V mutant) or a rise in the energetic contribution from the electrostatic interactions between the inhibitor and its target (I54M and I54V mutant). For the V32I mutant, also an increased free energy for the polar solvation contributes to the drug resistance. For the L90M mutant, a rise in the van der Waals energy for APV-PR interactions is compensated by a decrease in the polar solvation free energy such that the net binding affinity remains unchanged. Detailed understanding of the molecular forces governing binding and drug resistance might assist in the design of new inhibitors against HIV-1 PR variants that are resistant against current drugs.     

Interpretation of protein/ligand crystal structure using QM/MM calculations: case of HIV-1 protease/metallacarborane complex

Deltahedral metallacarborane compounds have recently been discovered as potent, specific, stable, and nontoxic inhibitors of HIV-1 protease (PR), the major target for AIDS therapy. The 2.15 A-resolution X-ray structure has exhibited a nonsymmetrical binding of the parental compound [Co(3+)-(C2B9H11)2](-) (GB-18) into PR dimer and a symmetrical arrangement in the crystal of two PR dimer complexes into a tetramer. In order to explore structural and energetic details of the inhibitor binding, quantum mechanics coupled with molecular mechanics approach was utilized. Realizing the close positioning of anionic inhibitors in the active site cavity, the possibility of an exchange of structural water molecules Wat50 and Wat128 by Na+ counterions was studied. The energy profiles for the rotation of the GB-18 molecules along their longitudinal axes in complex with PR were calculated. The results show that two Na+ counterions are present in the active site cavity and provide energetically favorable and unfavorable positions for carbon atoms within the carborane cages. Eighty-one rotamer combinations of four molecules of GB-18 bound to PR out of 4 x 10(5) are predicted to be highly populated. These results lay ground for further calculations of interaction energies between GB-18 and amino acids of PR active site and will make it possible to interpret computationally the binding of similar metallacarborane molecules to PR as well as to resistant PR variants. Moreover, this computational tool will allow the design of new, more potent metallacarborane-based HIV-1 protease inhibitors.     

Strategy for discovering chemical inhibitors of human cyclophilin a: focused library design, virtual screening, chemical synthesis and bioassay

The discovery of cyclophilin A (CypA) inhibitor is now of special interest in the treatment of immunological disorders. In this work, using a strategy integrating focused combinatorial library design, virtual screening, chemical synthesis, and bioassay, a series of novel small molecular CypA inhibitors have been discovered. First, using the fragments taken from our previously discovered CypA inhibitors (Bioorg. Med. Chem. 2006, 14, 2209-2224) as building blocks, we designed a focused combinatorial library containing 255 molecules employing the LD1.0 program (J. Comb. Chem. 2005, 7, 398-406) developed by us. Sixteen compounds (1a-e, 2a-b, 3a-b, and 4a-g) were selected by using virtual screening against the X-ray crystal structure of CypA as well as druglike analysis for further synthesis and bioassay. All these sixteen molecules are CypA binders with binding affinities (K(D) values) ranging from 0.076 to 41.0 microM, and five of them (4a, 4c, and 4e-g) are potent CypA inhibitors with PPIase inhibitory activities (IC(50) values) of 0.25-6.43 microM. The hit rates for binders and inhibitors are as high as 100% and 31.25%, respectively. Remarkably, both the binding affinity and inhibitory activity of the most potent compound increase approximately 10 times than that of the most active compound discovered previously. The high hit rate and the high potency of the new CypA inhibitors demonstrated the efficiency of the strategy for focused library design and screening. In addition, the novel chemical entities reported in this study could be leads for discovering new therapies against the CypA pathway.     

Does a diol cyclic urea inhibitor of HIV-1 protease bind tighter than its corresponding alcohol form? A study by free energy perturbation and continuum electrostatics calculations

The cyclic urea inhibitors of HIV-1 protease generally have two hydroxyl groups on the seven-membered ring. In this study, free energy perturbation and continuum electrostatic calculations were used to study the contributions of the two hydroxyl groups to the binding affinity and solubility of a cyclic urea inhibitor DMP323. The results indicated that the inhibitor with one hydroxyl group has better binding affinity and solubility than the inhibitor with two hydroxyl groups. Therefore, removal of one hydroxyl group from DMP323 may help to improve the properties of DMP323. This is also likely to be true for other cyclic urea inhibitors. The study also illustrated the difficulty in accurate modeling of the binding affinities of HIV-1 protease inhibitors, which involves many possible protonation states of the two catalytic aspartic acids in the active site of the enzyme.     

Heterocyclic HIV-1 protease inhibitors

[formula: see text] A series of simple heterocyclic HIV-1 protease inhibitors were developed on the basis of size, shape, and electronic complementarity to the active site of the enzyme. The C2-symmetric heterocycles do not contain a transition-state isostere nor are they active site directed irreversible inhibitors; thus, they represent the success of a new design strategy. The first generation heterocycles inhibit the protease in the micromolar range, whereas control compounds show no bioactivity at the same concentrations.     

Intricate Recognition of Glycolipid-Like Compounds by HIV-1 Envelope Proteins Evaluated with Surface Plasmon Resonance Imaging

Fusion of human immunodeficiency virus (HIV) to the cell membrane occurs by the specific binding of an envelope protein of HIV-1 (gp120 and gp160) and a glycosphingolipid of the cell membrane. In this study, quantitative and array-based affinity evaluation of gp120 and gp160 was performed by surface plasmon resonance (SPR) and the SPR imaging technique using a substrate immobilized with glycolipid-like compounds (Gb3, GM3, and Lac). Quantitative affinity evaluation showed that gp160 specifically bound to Gb3 and Lac compared with GM3, whereas gp120 showed lower binding affinity and specificity. Array-based evaluation showed that gp160 binds to Gb3 more favorably than Lac and GM3.     

SAR by MS

RNAs have recently emerged as an exciting new target for small molecule therapeutics. Conventional HTS discovery strategies measuring disruption of RNAprotein interactions have proven unsuccessful. We describe a ligand-based drug discovery strategy that addresses the inherent difficulties RNA targets. The strategy is based on: 1) using a MS spectrometry (MS)-based assay to measure the affinity of compounds for a target; 2) performing competitive binding experiments and molecular modeling with the motifs to determine the binding site(s) of the ligands; 3) design and synthesis of derivatives of interesting binders to establish the linking sites; 4) identifying the appropriate linker group using MS; 5) fusing motifs into a more complex structure to afford higher affinity compounds. Example of applying this strategy to identify new classes of lead molecules with affinity and specificity for ribosomal RNA targets will be presented.     

Combating susceptibility to drug resistance: lessons from HIV-1 protease

Drug resistance is a major obstacle in modern medicine. However, resistance is rarely considered in drug development and may inadvertently be facilitated, as many designed inhibitors contact residues that can mutate to confer resistance, without significantly impairing function. Contemporary drug design often ignores the detailed atomic basis for function and primarily focuses on disrupting the target's activity, which is necessary but not sufficient for developing a robust drug. In this study, we examine the impact of drug-resistant mutations in HIV-1 protease on substrate recognition and demonstrate that most primary active site mutations do not extensively contact substrates, but are critical to inhibitor binding. We propose a general, structure-based strategy to reduce the probability of drug resistance by designing inhibitors that interact only with those residues that are essential for function.     

The application of thermodynamic methods in drug design

The optimization of lead compounds as viable drug candidates involves the optimization of their binding affinity towards the selected target. The binding affinity, K_a, is determined by the Gibbs energy of binding, ΔG, which in turn is determined by the enthalpy, ΔH, and entropy, ΔS, changes (ΔG=ΔH−TΔS). In principle, many combinations of ΔH and ΔS values can give rise to the same ΔG value and, therefore, elicit the same binding affinity. However, enthalpically dominated ligands do not behave the same as entropically dominated ligands. Current paradigms in drug design usually generate highly hydrophobic and conformationally constrained ligands. The thermodynamic signature of these ligands is an entropically dominated binding affinity often accompanied by an unfavorable binding enthalpy. Conformationally constrained ligands cannot easily adapt to changes in the geometry of the binding site, being therefore highly susceptible to drug resistance mutations or naturally occurring genetic polymorphisms. The design of ligands with the capability to adapt to a changing target requires the introduction of certain elements of flexibility or the relaxation of some conformational constraints. Since these compounds pay a larger conformational entropy penalty upon binding, the optimization of their binding affinity requires the presence of a favorable binding enthalpy. In this paper, experimental and computational strategies aimed at identifying and optimizing enthalpic ligands will be discussed and applied to the case of HIV-1 protease inhibitors. It is shown that a thermodynamic guide to drug design permits the identification of drug candidates with a lower susceptibility to target mutations causing drug resistance.     

Enhancing protein backbone binding--a fruitful concept for combating drug-resistant HIV

The evolution of drug resistance is one of the most fundamental problems in medicine. In HIV/AIDS, the rapid emergence of drug-resistant HIV-1 variants is a major obstacle to current treatments. HIV-1 protease inhibitors are essential components of present antiretroviral therapies. However, with these protease inhibitors, resistance occurs through viral mutations that alter inhibitor binding, resulting in a loss of efficacy. This loss of potency has raised serious questions with regard to effective long-term antiretroviral therapy for HIV/AIDS. In this context, our research has focused on designing inhibitors that form extensive hydrogen-bonding interactions with the enzyme's backbone in the active site. In doing so, we limit the protease's ability to acquire drug resistance as the geometry of the catalytic site must be conserved to maintain functionality. In this Review, we examine the underlying principles of enzyme structure that support our backbone-binding concept as an effective means to combat drug resistance and highlight their application in our recent work on antiviral HIV-1 protease inhibitors.     

Potent and selective peptidyl boronic acid inhibitors of the serine protease prostate-specific antigen

Prostate cancer cells produce high (microgram to milligram/milliliter) levels of the serine protease Prostate-Specific Antigen (PSA). PSA is enzymatically active in the extracellular fluid surrounding prostate cancers but is found at 1,000- to 10,000-fold lower concentrations in the circulation, where it is inactivated due to binding to abundant serum protease inhibitors. The exclusive presence of high levels of active PSA within prostate cancer sites makes PSA an attractive candidate for targeted imaging and therapeutics. A synthetic approach based on a peptide substrate identified first peptide aldehyde and then boronic acid inhibitors of PSA. The best of these had the sequence Cbz-Ser-Ser-Lys-Leu-(boro)Leu, with a K_i for PSA of 65 nM. The inhibitor had a 60-fold higher K_i for chymotrypsin. A validated model of PSA's catalytic site confirmed the critical interactions between the inhibitor and residues within the PSA enzyme.     

Binding of novel fullerene inhibitors to HIV-1 protease: insight through molecular dynamics and molecular mechanics Poisson–Boltzmann surface area calculations

The objectives of this study include the design of a series of novel fullerene-based inhibitors for HIV-1 protease (HIV-1 PR), by employing two strategies that can also be applied to the design of inhibitors for any other target. Additionally, the interactions which contribute to the observed exceptionally high binding free energies were analyzed. In particular, we investigated: (1) hydrogen bonding (H-bond) interactions between specific fullerene derivatives and the protease, (2) the regions of HIV-1 PR that play a significant role in binding, (3) protease changes upon binding and (4) various contributions to the binding free energy, in order to identify the most significant of them. This study has been performed by employing a docking technique, two 3D-QSAR models, molecular dynamics (MD) simulations and the molecular mechanics Poisson–Boltzmann surface area (MM–PBSA) method. Our computed binding free energies are in satisfactory agreement with the experimental results. The suitability of specific fullerene derivatives as drug candidates was further enhanced, after ADMET (absorption, distribution, metabolism, excretion and toxicity) properties have been estimated to be promising. The outcomes of this study revealed important protein–ligand interaction patterns that may lead towards the development of novel, potent HIV-1 PR inhibitors.     

Rational design and synthesis of a novel class of active site-targeted HIV protease inhibitors containing a hydroxymethylcarbonyl isostere. Use of phenylnorstatine or allophenylnorstatine as a transition-state mimic.

A novel class of HIV-1 protease inhibitors containing a hydroxymethylcarbonyl (HMC) isostere were designed from the substrate transition state and synthesized. Phenylnorstatine [Pns; (2R,3S)-3-amino-2-hydroxy-4-phenylbutyric acid] and the 2S diastereomer, (2S,3S)-3-amino-2-hydroxy-4-phenylbutyric acid, named allophenylnorstatine (Apns) were effective transition-state mimics, and incorporation of Pns-Pro or Apns-Pro at the P1-P1' site gave potent and specific HIV-1 protease inhibitors. In the inhibitory assays, the chemically synthesized [Ala67,95] HIV-1 protease was used.     

Optimization of bivalent glutathione S-transferase inhibitors by combinatorial linker design

Dimeric glutathione S-transferases (GSTs) are pharmacological targets for several diseases, including cancer. Isoform specificity has been difficult to achieve due to their overlapping substrate selectivity. Here we demonstrate the utility of bivalent GST inhibitors and their optimization via combinatorial linker design. A combinatorial library with dipeptide linkers emanating symmetrically from a central scaffold (bis-3,5-aminomethyl benzoic acid, AMAB) to connect two ethacrynic acid moieties was prepared and decoded via iterative deconvolution, against the isoforms GSTA1-1 and GSTP1-1. The library yielded high affinity GSTA1-1 selective inhibitors (70-120-fold selectivity) and with stoichiometry of one inhibitor: one GSTA1-1 dimer. Saturation Transfer Difference (STD) NMR with one of these inhibitors, with linker structure (Asp-Gly-AMAB-Gly-Asp) and K(D) = 42 nM for GSTA1-1, demonstrates that the Asp-Gly linker interacts tightly with GSTA1-1, but not P1-1. H/D exchange mass spectrometry was used to map the protein binding site and indicates that peptides within the intersubunit cleft and in the substrate binding site are protected by inhibitor from solvent exchange. A model is proposed for the binding orientation of the inhibitor, which is consistent with electrostatic complementarity between the protein cleft and inhibitor linker as the source of isoform selectivity and high affinity. The results demonstrate the utility of combinatorial, or "irrational", linker design for optimizing bivalent inhibitors.     

Charge optimization of the interface between protein kinases and their ligands

Examining the potential for electrostatic complementarity between a ligand and a receptor is a useful technique for rational drug design, and can demonstrate how a system prioritizes interactions when allowed to optimize its charge distribution. In this computational study, we implemented the previously developed, continuum solvent-based charge optimization theory with a simple, quadratic programming algorithm and the UHBD Poisson-Boltzmann solver. This method allows one to compute the best set of point charges for a ligand or ligand region based on the ligand and receptor shape, and the receptor partial charges, by optimizing the binding free energy obtained from a continuum-solvent model. We applied charge optimization to a fragment of the heat-stable protein kinase inhibitor (PKI) of protein kinase A (PKA), to three flavopiridol inhibitors of CDK2, and to cyclin A which interacts with CDK2 to regulate the cell cycle. We found that a combination of global (involving every charge) and local (involving only charges in a local region) optimization can give useful hints for designing better inhibitors. Although some parts of an inhibitor may already contribute significantly to binding, we found that they could still be the most important targets for modifications to obtain stronger binders. In studying the binding of flavopiridol inhibitors to CDK2, comparable binding affinity could be obtained regardless of whether the net charges of the inhibitors were constrained to -2, -1, 0, 1, or 2 during the optimization. This provides flexibility in inhibitor design when a certain net charge of the inhibitor is desired in addition to strong binding affinity. For the study of the PKA-PKI and CDK2-cyclin A interfaces, we identified residues whose charge distributions are already close to optimal and those whose charge distributions could be refined to further improve binding.     

Drug-resistant molecular mechanism of CRF01_AE HIV-1 protease due to V82F mutation

Human immunodeficiency virus type 1 protease (HIV-1 PR) is one of the major targets of anti-AIDS drug discovery. The circulating recombinant form 01 A/E (CRF01_AE, abbreviated AE) subtype is one of the most common HIV-1 subtypes, which is infecting more humans and is expanding rapidly throughout the world. It is, therefore, necessary to develop inhibitors against subtype AE HIV-1 PR. In this work, we have performed computer simulation of subtype AE HIV-1 PR with the drugs lopinavir (LPV) and nelfinavir (NFV), and examined the mechanism of resistance of the V82F mutation of this protease against LPV both structurally and energetically. The V82F mutation at the active site results in a conformational change of 79′s loop region and displacement of LPV from its proper binding site, and these changes lead to rotation of the side-chains of residues D25 and I50′. Consequently, the conformation of the binding cavity is deformed asymmetrically and some interactions between PR and LPV are destroyed. Additionally, by comparing the interactive mechanisms of LPV and NFV with HIV-1 PR we discovered that the presence of a dodecahydroisoquinoline ring at the P1′ subsite, a [2-(2,6-dimethylphenoxy)acetyl]amino group at the P2′ subsite, and an N2 atom at the P2 subsite could improve the binding affinity of the drug with AE HIV-1 PR. These findings are helpful for promising drug design. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Protein-Templated Hit Identification through an Ugi Four-Component Reaction**

Kinetic target-guided synthesis represents an efficient hit-identification strategy, in which the protein assembles its own inhibitors from a pool of complementary building blocks via an irreversible reaction. Herein, we pioneered an in situ Ugi reaction for the identification of novel inhibitors of a model enzyme and binders for an important drug target, namely, the aspartic protease endothiapepsin and the bacterial β-sliding clamp DnaN, respectively. Highly sensitive mass-spectrometry methods enabled monitoring of the protein-templated reaction of four complementary reaction partners, which occurred in a background-free manner for endothiapepsin or with a clear amplification of two binders in the presence of DnaN. The Ugi products we identified show low micromolar activity on endothiapepsin or moderate affinity for the β-sliding clamp. We succeeded in expanding the portfolio of chemical reactions and biological targets and demonstrated the efficiency and sensitivity of this approach, which can find application on any drug target.