Inhibition of human immunodeficiency virus-1 protease by a C2-symmetric phosphinate. Synthesis and crystallographic analysis

Multimolecular Self-Organization of 1-Acetyl-1,3-bis(haloarylamines) in KOH/DMSO System: From Acetylene Gas and o-Halo Arylamines toward a Higher Molecular Complexity and Diversity

An unprecedented self-organization of three molecules of acetylene, two molecules of o-halo arylamine, and one molecule of water to 1-acetyl-1,3-bis(haloarylamines) in the KOH/DMSO system has been discovered. These high-value 1,3-bis(arylamines) are capable of acting as efficient "platform molecules" to provide one-pot access to hexahydropyrrolo[3,2-b]indoles, 3-haloarylamino-pyrrolidines, benzyliden piperidinols, and other functionalized diamines.

Non-symmetrical bis(aminoalkyl)phosphinates: new ligands with enhanced binding of Cu(ii) ions

Novel, non-symmetrical bis(aminoalkyl)phosphinic acids exhibit enhanced efficiency in Cu(ii) ion binding.

Inhibitors or toxins? Large library target-specific screening of fullerene-based nanoparticles for drug design purpose

Fullerene-based nanoparticles have been the subject of vital interest due to their unique properties and potential application in many areas, including medicine. Here we explore their characteristics that could make them prospective leads for known disease-related proteins. High-throughput virtual screening supported by comprehensive multi-software protein-ligand docking simulation and cheminformatics approaches has been applied in investigation of interactions of 1117 proteins with a 169 fullerene nanoparticles decorated with different small molecules. Moreover, obtained docking results were confirmed by the series of unrestricted all-atom molecular dynamics (MD) simulations. Hydrophobicity of fullerene core along with hydrophilic interaction of side chains plays a key role in binding with the studied proteins. We identified a series of nanoparticles that can lead to development of robust drugs for target proteins and another series that can behave as a highly toxic agent. The structure-activity relationship analysis revealed two significant molecular properties responsible for the binding score values. The application of carefully selected computational techniques and described outcome of the study facilitate development of functional fullerene nanoparticles for drug-like and drug delivery applications.

Virtual drug screen schema based on multiview similarity integration and ranking aggregation

The current drug virtual screen (VS) methods mainly include two categories. i.e., ligand/target structure-based virtual screen and that, utilizing protein-ligand interaction fingerprint information based on the large number of complex structures. Since the former one focuses on the one-side information while the later one focuses on the whole complex structure, they are thus complementary and can be boosted by each other. However, a common problem faced here is how to present a comprehensive understanding and evaluation of the various virtual screen results derived from various VS methods. Furthermore, there is still an urgent need for developing an efficient approach to fully integrate various VS methods from a comprehensive multiview perspective. In this study, our virtual screen schema based on multiview similarity integration and ranking aggregation was tested comprehensively with statistical evaluations, providing several novel and useful clues on how to perform drug VS from multiple heterogeneous data sources. (1) 18 complex structures of HIV-1 protease with ligands from the PDB were curated as a test data set and the VS was performed with five different drug representations. Ritonavir ( 1HXW ) was selected as the query in VS and the weighted ranks of the query results were aggregated from multiple views through four similarity integration approaches. (2) Further, one of the ranking aggregation methods was used to integrate the similarity ranks calculated by gene ontology (GO) fingerprint and structural fingerprint on the data set from connectivity map, and two typical HDAC and HSP90 inhibitors were chosen as the queries. The results show that rank aggregation can enhance the result of similarity searching in VS when two or more descriptions are involved and provide a more reasonable similarity rank result. Our study shows that integrated VS based on multiple data fusion can achieve a remarkable better performance compared to that from individual ones and, thus, serves as a promising way for efficient drug screening, taking advantages of the rapidly accumulated molecule representations and heterogeneous data in the pharmacological area.

Discovery of potent HIV-1 protease inhibitors incorporating sulfoximine functionality

Based on the unique property of sulfoximine and the homodimeric C(2) structural symmetry of HIV-1 protease, a novel class of sulfoximine-based pseudosymmetric HIV-1 protease inhibitors was designed and synthesized. The sulfoximine moiety was demonstrated to be important for HIV-1 protease inhibitor potency. The most active stereoisomer (2S,2'S) displays a potency of 2.5 nM (IC(50)) against HIV-1 protease and an anti-HIV-1 activity of 408 nM (IC(50)). A possible mode of action is proposed.

How inaccuracies in protein structure models affect estimates of protein-ligand interactions: computational analysis of HIV-I protease inhibitor binding

The influence of possible inaccuracies that can arise during homology modeling of protein structures used for ligand binding studies were investigated with the molecular mechanics generalized Born surface area (MM-GBSA) method. For this, a family of well-characterized HIV-I protease-inhibitor complexes was used. Validation of MM-GBSA led to a correlation coefficient ranging from 0.72 to 0.93 between calculated and experimental binding free energies DeltaG. All calculated DeltaG values were based on molecular dynamics simulations with explicit solvent. Errors introduced into the protein structure through misplacement of side-chains during rotamer modeling led to a correlation coefficient between DeltaG(calc) and DeltaG(exp) of 0.75 compared with 0.90 for the correctly placed side chains. This is in contrast to homology models for members of the retroviral protease family with template structures ranging in sequence identity between 32% and 51%. For these protein models, the correlation coefficients vary between 0.84 and 0.87, which is considerably closer to the original protein (0.90). It is concluded that HIV-I low sequence identity with the template structure still allows creating sufficiently reliable homology models to be used for ligand-binding studies, although placement of the rotamers is a critical step during the modeling.

Synthesis and activity of phosphinic tripeptide inhibitors of cathepsin C

Phosphinic tripeptide analogues Gly-Xaapsi[P(O)(OH)CH(2)]-Gly have been developed as inhibitors of cathepsin C (DPP I), a lysosomal, papain-like cysteine protease. The target compounds were synthesised by addition of methyl acrylate to the appropriate phosphinic acids followed by the N-terminus elongation using mixed anhydride procedure. The latter step has been demonstrated to be a suitable method for N-terminal extension of the phosphinic pseudopeptide analogues without requirement of hydroxyphosphinyl protection. The title compounds appeared to be moderate inhibitors of the cathepsin C. However, although designed as transition state analogues, they surprisingly exhibited noncompetitive mode of binding to cathepsin C. Differences in kinetics of C-terminal acids and esters have been additionally observed.

Analysis and optimization of structure-based virtual screening protocols. (3). New methods and old problems in scoring function design

Scoring function research remains a primary focus of current structure-based virtual screening (SVS) technology development. Here, we present an alternative method for scoring function design that attempts to combine crystallographic structural information with data derived from directly within SVS calculations. The technique utilizes a genetic algorithm (GA) to optimize functions based on binding property data derived from multiple virtual screening calculations. These calculations are undertaken on protein data bank (PDB) complex active sites using ligands of known binding mode in conjunction with "noise" compounds. The advantages of such an approach are that the function does not rely on assay data and that it can potentially use the "noise" binding data to recognize the sub-optimal docking interactions inherent in SVS calculations. Initial efforts in technique exploration using DOCK are presented, with comparisons made to existing DOCK scoring functions. An analysis of the problems inherent to scoring function development is also made, including issues in dataset creation and limitations in descriptor utility when viewed from the perspective of docking mode resolution. The future directions such studies might take are also discussed in detail.

Relation between sequence and structure of HIV-1 protease inhibitor complexes: a model system for the analysis of protein flexibility

The flexibility of different regions of HIV-1 protease was examined by using a database consisting of 73 X-ray structures that differ in terms of sequence, ligands or both. The root-mean-square differences of the backbone for the set of structures were shown to have the same variation with residue number as those obtained from molecular dynamics simulations, normal mode analyses and X-ray B-factors. This supports the idea that observed structural changes provide a measure of the inherent flexibility of the protein, although specific interactions between the protease and the ligand play a secondary role. The results suggest that the potential energy surface of the HIV-1 protease is characterized by many local minima with small energetic differences, some of which are sampled by the different X-ray structures of the HIV-1 protease complexes. Interdomain correlated motions were calculated from the structural fluctuations and the results were also in agreement with molecular dynamics simulations and normal mode analyses. Implications of the results for the drug-resistance engendered by mutations are discussed briefly.

Synthesis and comparative study on the reactivity of peptidyl-type phosphinic esters: intramolecular effects in the alkaline and acidic cleavage of methyl beta-carboxyphosphinates

Using the phosphinic analogue of Cbz-Phe-Gly-OEt 1a as a template for this study, several phosphinic esters (2a-g) were prepared, employing an efficient method for each case. The reactivity of these derivatives under conventional deprotection conditions was studied, and the results are listed comparatively. The effect of steric hindrance as well as the contribution of neighboring groups in the rate of hydrolysis of suitably selected beta-carboxyphosphinates under acidic and alkaline deprotection conditions was examined. The results clearly demonstrate that a significant acceleration of phosphinate cleavage occurs due to the intermediacy of a five-membered, mixed anhydride-type species. This was supported by the observation that similar interactions were not observed in the case of hindered alpha-carboxyphosphinate homologous derivatives.

How does a symmetric dimer recognize an asymmetric substrate? A substrate complex of HIV-1 protease

The crystal structure of an actual HIV-1 protease-substrate complex is presented at 2.0 A resolution (R-value of 19.7 % (R(free) 23.3 %)) between an inactive variant (D25N) of HIV-1 protease and a long substrate peptide, Lys-Ala-Arg-Val-Leu-Ala-Glu-Ala-Met-Ser, which covers a full binding epitope of capsid(CA)-p2, cleavage site. The substrate peptide is asymmetric in both size and charge distribution. To accommodate this asymmetry the two protease monomers adopt different conformations burying a total of 1038 A(2) of surface area at the protease-substrate interface. The specificity for the CA-p2 substrate peptide is mainly hydrophobic, as most of the hydrogen bonds are made with the backbone of the peptide substrate. Two water molecules bridge the two monomers through the loops Gly49-Gly52 (Gly49'-Gly52') and Pro79'-Val82' (Pro79-Val82). When other complexes are compared, the mobility of these loops is correlated with the content of the P1 and P1' sites. Interdependence of the conformational changes allows the protease to exhibit its wide range of substrate specificity.

DREAM++: flexible docking program for virtual combinatorial libraries

We present a set of programs, DREAM+2 (Docking and Reaction programs using Efficient seArch Methods written in C++), for docking computationally generated ligands into macromolecular binding sites. DREAM++ is composed of three programs: ORIENT++, REACT++ and SEARCH++. The program ORIENT++ positions molecules in a binding site with the DOCK algorithm. Its output can be used as input to REACT++ and SEARCH+2. The program REACT++ performs user-specific chemical reactions on a docked molecule, so that reaction products can be evaluated for three dimensional complementarity with the macromolecular site. The program SEARCH++ performs an efficient conformation search on the reaction products using a hybrid backtrack and incremental construction algorithm. We have applied the programs to HIV protease-inhibitor complexes as test systems. We found that we can differentiate high-affinity ligands based on several measures: interaction energies, occupancy of protein subsites and the number of successfully docked conformations for each product. Encouraged by the results in the test case, we applied the programs to propose novel inhibitors of HIV protease. These inhibitors can be generated by organic reactions using commercially available reagents. They are alternatives to the inhibitors synthesized by Glaxo.

Solvation studies of DMP323 and A76928 bound to HIV protease: analysis of water sites using grand canonical Monte Carlo simulations

We examine the water solvation of the complex of the inhibitors DMP323 and A76928 bound to HIV-1 protease through grand canonical Monte Carlo simulations, and demonstrate the ability of this method to reproduce crystal waters and effectively predict water positions not seen in the DMP323 or A76928 structures. The simulation method is useful for identifying structurally important waters that may not be resolved in the crystal structures. It can also be used to identify water positions around a putative drug candidate docked into a binding pocket. Knowledge of these water positions may be useful in designing drugs to utilize them as bridging groups or displace them in the binding pocket. In addition, the method should be useful in finding water sites in homology models of enzymes for which crystal structures are unavailable.

Design and synthesis of chiral and racemic phosphonate-based haptens for the induction of aldolase catalytic antibodies

A novel strategy for the generation of aldolase catalytic antibodies, based on the use of antibody-catalyzed enol ester hydrolysis as a 'trigger' to generate a reactive enolate intermediate, is described. A model system to test this strategy was developed and substrate 8 was synthesized. However, the targeted bifunctional haptens 11 and 33 were synthetically inaccessible, and therefore the alternative phosphonate hapten 39 was prepared. The key step in the synthesis of 39 was the direct generation of an unprotected phosphonate precursor via coupling of the secondary alcohol 37 with CH3P(O)Cl2. The chiral counterpart of hapten 39 was also synthesized from alcohol 46, prepared by Corey's asymmetric reduction method. One polyclonal antibody preparation generated from 39 appeared to catalyze the hydrolysis of the secondary acetate 49, but not the desired aldol cyclization of 8. Possible rationales for this finding are discussed.

Molecular docking to ensembles of protein structures

Until recently, applications of molecular docking assumed that the macromolecular receptor exists in a single, rigid conformation. However, structural studies involving different ligands bound to the same target biomolecule frequently reveal modest but significant conformational changes in the target. In this paper, two related methods for molecular docking are described that utilize information on conformational variability from ensembles of experimental receptor structures. One method combines the information into an "energy-weighted average" of the interaction energy between a ligand and each receptor structure. The other method performs the averaging on a structural level, producing a "geometry-weighted average" of the inter-molecular force field score used in DOCK 3.5. Both methods have been applied in docking small molecules to ensembles of crystal and solution structures, and we show that experimentally determined binding orientations and computed energies of known ligands can be reproduced accurately. The use of composite grids, when conformationally different protein structures are available, yields an improvement in computational speed for database searches in proportion to the number of structures.

Engineering human immunodeficiency virus 1 protease heterodimers as macromolecular inhibitors of viral maturation

Dimerization of human immunodeficiency virus type 1 protease (HIV-1 PR) monomers is an essential prerequisite for viral proteolytic activity and the subsequent generation of infectious virus particles. Disruption of the dimer interface inhibits this activity as does formation of heterodimers between wild-type and defective monomers. A structure-based approach was used to identify amino acid substitutions at the dimer interface of HIV-1 PR that facilitate preferential association of heterodimers and inhibit self-association of the defective monomers. Expression of the designed PR monomers inhibits activity of wild-type HIV-1 PR and viral infectivity when assayed in an ex vivo model system. These results show that it is possible to design PR monomers as macromolecular inhibitors that may provide an alternative to small molecule inhibitors for the treatment of HIV infection.

A new class of HIV-1 protease inhibitor: the crystallographic structure, inhibition and chemical synthesis of an aminimide peptide isostere

The essential role of HIV-1 protease (HIV-1 PR) in the viral life cycle makes it an attractive target for the development of substrate-based inhibitors that may find efficacy as anti-AIDS drugs. However, resistance has arisen to potent peptidomimetic drugs necessitating the further development of novel chemical backbones for diversity based chemistry focused on probing the active site for inhibitor interactions and binding modes that evade protease resistance. AQ148 is a potent inhibitor of HIV-1 PR and represents a new class of transition state analogues incorporating an aminimide peptide isostere. A 3-D crystallographic structure of AQ148, a tetrapeptide isostere, has been determined in complex with its target HIV-1 PR to a resolution of 2.5 A and used to evaluate the specific structural determinants of AQ148 potency and to correlate structure-activity relationships within the class of related compounds. AQ148 is a competitive inhibitor of HIV-1 PR with a Ki value of 137 nM. Twenty-nine derivatives have been synthesized and chemical modifications have been made at the P1, P2, P1', and P2' sites. The atomic resolution structure of AQ148 bound to HIV-1 PR reveals both an inhibitor binding mode that closely resembles that of other peptidomimetic inhibitors and specific protein/inhibitor interactions that correlate with structure-activity relationships. The structure provides the basis for the design, synthesis and evaluation of the next generation of hydroxyethyl aminimide inhibitors. The aminimide peptide isostere is a scaffold with favorable biological properties well suited to both the combinatorial methods of peptidomimesis and the rational design of potent and specific substrate-based analogues.

Synthesis of peptide analogues containing phosphonamidate methyl ester functionality: HIV-1 proteinase inhibitors possessing unique cell uptake properties

Stereochemically defined and epimeric phosphonamidate methyl ester-containing peptide analogues were synthesised and were found to be moderate inhibitors of the HIV-1 proteinase. All of the analogues containing the phosphonamidate ester grouping showed a marked ability to enter cells, as highlighted by the approximate equivalence of the IC50 values for enzyme inhibition in solution and inhibition of HIV-1 replication in virus infected cells.

An orally bioavailable HIV-1 protease inhibitor containing an imidazole-derived peptide bond replacement: crystallographic and pharmacokinetic analysis

(2R,4S,5S,1'S)-2-Phenylmethyl-4-hydroxy-5-(tert-butoxycarbonyl) amino-6-phenylhexanoyl-N-(1'-imidazo-2-yl)-2'-methylpropanamide (compound 2) is a tripeptide analogue inhibitor of HIV-1 protease in which a C-terminal imidazole substituent constitutes an isoelectronic, structural mimic of a carboxamide group. Compound 2 is a potent inhibitor of the protease (K(i) = 18 nM) and inhibits HIV-1 acute infectivity of CD4+ T-lymphocytes (IC50 = 570 nM). Crystallographic analysis of an HIV-1 protease-compound 2 complex demonstrates that the nitrogen atoms of the imidazole ring assume the same hydrogen-bonding interactions with the protease as amide linkages in other peptide analogue inhibitors. The sole substitution of the C-terminal carboxamide of a hydroxyethylene-containing tripeptide analogue with an imidazole group imparts greatly improved pharmacokinetic and oral bioavailability properties on the compound compared to its carboxamide-containing homologue (compound 1). While the oral bioavailability of compound 1 in rats was negligible, compound 2 displayed oral bioavailabilities of 30% and 14%, respectively, in rats and monkeys.

Electrostatic complementarity between proteins and ligands. 2. Ligand moieties

Drug design strategies consider factors governing intermolecular interactions to build up putative ligands. In many strategies, the ligand is constructed using fragments which are placed in the site sequentially. The optimization is then performed with each fragment. We would like to examine if this optimization strategy could generate ligands with optimal electrostatic interactions. The electrostatic complementarities between constituent moieties and the receptor site have been calculated. The whole-ligand complementarity does not appear to be the mathematical mean of the individual complementarities, nor have we found a simple relationship between the moiety and whole-ligand complementarities. The results demonstrate clearly that, using a simple model, it is very difficult to predict the electrostatic potential complementarity of the whole ligand from the complementarities of its constituent chemical moieties. This means that ligand design strategies must optimize the electrostatic complementarity globally, and not moiety by moiety.

Electrostatic complementarity between proteins and ligands. 1. Charge disposition, dielectric and interface effects

Electrostatic interactions have always been considered an important factor governing ligand-receptor interactions. Previous work in this field has established the existence of electrostatic complementarity between the ligand and its receptor site. However, this property has not been treated rigorously, and the description remains largely qualitative. In this work, 34 data sets of high quality were chosen from the Brookhaven Protein Databank. The electrostatic complementary has been calculated between the surface potentials; complementarity is absent between adjacent or neighbouring atoms of the ligand and the receptor. There is little difference between complementarities on the total ligand surface and the interfacial region. Altering the homogeneous dielectric to distance-dependent dielectrics reduces the complementarity slightly, but does not affect the pattern of complementarity.

Design of inhibitors of glycogen phosphorylase: a study of alpha- and beta-C-glucosides and 1-thio-beta-D-glucose compounds

alpha-D-Glucose is a weak inhibitor of glycogen phosphorylase b (Ki = 1.7 mM) and acts as a physiological regulator of hepatic glycogen metabolism. Glucose binds to phosphorylase at the catalytic site and results in a conformational change that stabilizes the inactive T state of the enzyme, promoting the action of protein phosphatase 1 and stimulating glycogen synthase. It has been suggested that, in the liver, glucose analogues with greater affinity for glycogen phosphorylase may result in a more effective regulatory agent. Several alpha- and beta-anhydroglucoheptonic acid derivatives and 1-deoxy-1-thio-beta-D-glucose analogues have been synthesized and tested in a series of crystallographic and kinetic binding studies with glycogen phosphorylase. The structural results of the bound enzyme-ligand complexes have been analyzed, together with the resulting affinities, in an effort to understand and exploit the molecular interactions that might give rise to a better inhibitor. This work has shown the following: (i) Similar affinities may be obtained through different sets of interactions. Specifically, in the case of the alpha- and beta-glucose-C-amides, similar Ki's (0.37 and 0.44 mM, respectively) are obtained with the alpha-anomer through interactions from the ligand via water molecules to the protein and with the beta-anomer through direct interaction from the ligand to the protein. Thus, hydrogen bonds through water can contribute binding energy similar to that of hydrogen bonds directly to the protein. (ii) Attempts to improve the inhibition by additional groups did not always lead to the expected result. The addition of nonpolar groups to the alpha-carboxamide resulted in a change in conformation of the pyranose ring from a chair to a skew boat and the consequent loss of favorable hydrogen bonds and increase in the Ki. (iii) The addition of polar groups to the alpha-carboxamide led to compounds with the chair conformation, and in the examples studied, it appears that hydration by a water molecule may provide sufficient stabilization to retain the chair conformation. (iv) The best inhibitor was N-methyl-beta-glucose-C-carboxamide (Ki = 0.16 mM), which showed a 46-fold improvement in Ki from the parent beta-D-glucose. The decrease in Ki may be accounted for by a single hydrogen bond from the amide nitrogen to a main-chain carbonyl oxygen, an increase in entropy through displacement of a water molecule, and favorable van der Waals contacts between the methyl substituent and nonpolar protein residues.(ABSTRACT TRUNCATED AT 250 WORDS)

Use of steady state kinetic methods to elucidate the kinetic and chemical mechanisms of retroviral proteases

Despite the current plethora of structural data of HIV-1 protease and the availability of potent inhibitors, whose structures are based in part on the presumed mechanism of action of this enzyme, our actual understanding of its chemical mechanism has been until now based largely on the precedents of the mammalian and fungal aspartic proteases and static three-dimensional data. The available steady state kinetic data of the protease, as reviewed here, constitute a first step in a detailed description of the mechanism of the enzyme to complement the structural data.

Use of nitrogen-15 kinetic isotope effects to elucidate details of the chemical mechanism of human immunodeficiency virus 1 protease

We have used 15N kinetic isotope effects of the HIV-1 protease-catalyzed peptidolysis of Ac-Ser-Gln-Asn-Tyr-Pro-Val-Val-NH2 to characterize the chemical mechanism of this enzyme. In addition, the multiple isotope effects have been determined by measuring the 15N kinetic isotope effects in both H2O and D2O. The isotope effects, measured on values of V/K, were determined by the incorporation of a radiolabel (tritium and 14C in peptides bearing the heavy and light isotopes, respectively) at a position remote from the isotopically labeled scissile peptide bond, such that the isotope effect was determined by measurement of the change in the 14C/3H ratio in recovered substrates at various fractions of reaction. At pH = 6.0 (37 degrees C), the nitrogen isotope effects were slightly, but significantly, inverse in both solvents: 15(V/K)H2O = 0.995 +/- 0.002, and 15(V/K)D2O = 0.992 +/- 0.003. The observation of an inverse nitrogen kinetic isotope effect implies that bonding to the nitrogen atom is becoming stiffened in a reaction transition state, and since this inverse isotope effect is enhanced in D2O, this isotope effect likely arises from protonation of the proline nitrogen atom.(ABSTRACT TRUNCATED AT 250 WORDS)