Hydroxyethylene isostere inhibitors of human immunodeficiency virus-1 protease: structure-activity analysis using enzyme kinetics, X-ray crystallography, and infected T-cell assays.

Analogues of peptides ranging in size from three to six amino acids and containing the hydroxyethylene dipeptide isosteres Phe psi Gly, Phe psi Ala, Phe psi NorVal, Phe psi Leu, and Phe psi Phe, where psi denotes replacement of CONH by (S)-CH(OH)CH2, were synthesized and studied as HIV-1 protease inhibitors. Inhibition constants (Ki) with purified HIV-1 protease depend strongly on the isostere in the order Phe psi Gly greater than Phe psi Ala greater than Phe psi NorVal greater than Phe psi Leu greater than Phe psi Phe and decrease with increasing length of the peptide analogue, converging to a value of 0.4 nM. Ki values are progressively less dependent on inhibitor length as the size of the P1' side chain within the isostere increases. The structures of HIV-1 protease complexed with the inhibitors Ala-Ala-X-Val-Val-OMe, where X is Phe psi Gly, Phe psi Ala, Phe psi NorVal, and Phe psi Phe, have been determined by X-ray crystallography (resolution 2.3-3.2 A). The crystals exhibit symmetry consistent with space group P6(1) with strong noncrystallographic 2-fold symmetry, and the inhibitors all exhibit 2-fold disorder. The inhibitors bind in similar conformations, forming conserved hydrogen bonds with the enzyme. The Phe psi Gly inhibitor adopts an altered conformation that places its P3' valine side chain partially in the hydrophobic S1' pocket, thus suggesting an explanation for the greater dependence of the Ki value on inhibitor length in the Phe psi Gly series. From the kinetic and crystallographic data, a minimal inhibitor model for tight-binding inhibition is derived in which the enzyme subsites S2-S2' are optimally occupied. The Ki values for several compounds are compared with their potencies as inhibitors of proteolytic processing in T-cell cultures chronically infected with HIV-1 (MIC values) and as inhibitors of acute infectivity (IC50 values). There is a rank-order correspondence, but a 20-1000-fold difference, between the values of Ki and those of MIC or IC50. IC50 values can approach those of Ki but are highly dependent on the conditions of the acute infectivity assay and are influenced by physiochemical properties of the inhibitors such as solubility.     

β-turn tendency in N-methylated peptides with dehydrophenylalanine residue: DFT study

The tendency to adopt β‐turn conformation by model dipeptides with α,β‐dehydrophenylalanine (ΔPhe) residue in the gas phase and in solution is investigated by theoretical methods. We pay special attention to a dependence of conformational properties on the side‐chain configuration of dehydro residue and the influence of N‐methylation on β‐turn stability. An extensive computational study of the conformational preferences of Z and E isomers of dipeptides Ac‐Gly‐(E/Z)‐ΔPhe‐NHMe ( 1a / 1b ) and Ac‐Gly‐(E/Z)‐ΔPhe‐NMe₂ ( 2a / 2b ) by B3LYP/6‐311++G(d,p) and MP2/6‐311++G(d,p) methods is reported. It is shown that, in agreement with experimental data, Ac‐Gly‐(Z)‐ΔPhe‐NHMe has a great tendency to adopt β‐turn conformation. In the gas phase the type II β‐turn is preferred, whereas in the polar environment, the type I. On the other hand, dehydro residue in Ac‐Gly‐(E)‐ΔPhe‐NHMe has a preference to adopt extended conformations in all environments. N‐methylation of C‐terminal amide group, which prevents the formation of 1←4 intramolecular hydrogen bond, change dramatically the conformational properties of studied dehydropeptides. Especially, the tendency to adopt β‐turn conformations is much weaker for the N‐methylated Z isomer (Ac‐Gly‐(Z)‐ΔPhe‐NMe₂), both in vacuo and in the polar environment. On the contrary, N‐methylated E isomer (Ac‐Gly‐(E)‐ΔPhe‐NMe₂) can easier adopt β‐turn conformation, but the backbone torsion angles (?₁, ψ₁, ?₂, ψ₂) are off the limits for common β‐turn types. © 2012 Wiley Periodicals, Inc. Biopolymers 97:518–528, 2012.     

Crystal structures of multidrug-resistant HIV-1 protease in complex with two potent anti-malarial compounds

Two potent inhibitors (compounds 1 and 2) of malarial aspartyl protease, plasmepsin-II, were evaluated against wild type (NL4-3) and multidrug-resistant clinical isolate 769 (MDR) variants of human immunodeficiency virus type-1 (HIV-1) aspartyl protease. Enzyme inhibition assays showed that both 1 and 2 have better potency against NL4-3 than against MDR protease. Crystal structures of MDR protease in complex with 1 and 2 were solved and analyzed. Crystallographic analysis revealed that the MDR protease exhibits a typical wide-open conformation of the flaps (Gly48 to Gly52) causing an overall expansion in the active site cavity, which, in turn caused unstable binding of the inhibitors. Due to the expansion of the active site cavity, both compounds showed loss of direct contacts with the MDR protease compared to the docking models of NL4-3. Multiple water molecules showed a rich network of hydrogen bonds contributing to the stability of the ligand binding in the distorted binding pockets of the MDR protease in both crystal structures. Docking analysis of 1 and 2 showed a decrease in the binding affinity for both compounds against MDR supporting our structure-function studies. Thus, compounds 1 and 2 show promising inhibitory activity against HIV-1 protease variants and hence are good candidates for further development to enhance their potency against NL4-3 as well as MDR HIV-1 protease variants.     

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.     

Effects of PRE and POST therapy drug-pressure selected mutations on HIV-1 protease conformational sampling

Conformational sampling of pre- and post-therapy subtype B HIV-1 protease sequences derived from a pediatric subject infected via maternal transmission with HIV-1 were characterized by double electron–electron resonance spectroscopy. The conformational ensemble of the PRE construct resembles native-like inhibitor bound states. In contrast, the POST construct, which contains accumulated drug-pressure selected mutations, has a predominantly semi-open conformational ensemble, with increased populations of open-like states. The single point mutant L63P, which is contained in PRE and POST, has decreased dynamics, particularly in the flap region, and also displays a closed-like conformation of inhibitor-bound states. These findings support our hypothesis that secondary mutations accumulate in HIV-1 protease to shift conformational sampling to stabilize open-like conformations, while maintaining the predominant semi-open conformation for activity.     

Bioactive conformations of two seminal delta opioid receptor penta‐peptides inferred from free‐energy profiles

Delta‐opioid (DOP) receptors are members of the G protein‐coupled receptor (GPCR) sub‐family of opioid receptors, and are evolutionarily related, with homology exceeding 70%, to cognate mu‐opioid (MOP), kappa‐opioid (KOP), and nociceptin opioid (NOP) receptors. DOP receptors are considered attractive drug targets for pain management because agonists at these receptors are reported to exhibit strong antinociceptive activity with relatively few side effects. Among the most potent analgesics targeting the DOP receptor are the linear and cyclic enkephalin analogs known as DADLE (Tyr‐D ‐Ala‐Gly‐Phe‐D ‐Leu) and DPDPE (Tyr‐D ‐Pen‐Gly‐Phe‐D ‐Pen), respectively. Several computational and experimental studies have been carried out over the years to characterize the conformational profile of these penta‐peptides with the ultimate goal of designing potent peptidomimetic agonists for the DOP receptor. The computational studies published to date, however, have investigated only a limited range of timescales and used over‐simplified representations of the solvent environment. We provide here a thorough exploration of the conformational space of DADLE and DPDPE in an explicit solvent, using microsecond‐scale molecular dynamics and bias‐exchange metadynamics simulations. Free‐energy profiles derived from these simulations point to a small number of DADLE and DPDPE conformational minima in solution, which are separated by relatively small energy barriers. Candidate bioactive forms of these peptides are selected from identified common spatial arrangements of key pharmacophoric points within all sampled conformations. © 2013 Wiley Periodicals, Inc. Biopolymers 101: 21–27, 2014.     

How conformational changes can affect catalysis,inhibition and drug resistance of enzymes with induced-fit binding mechanism such as the HIV-1 protease

A central question is how the conformational changes of proteins affect their function and the inhibition of this function by drug molecules. Many enzymes change from an open to a closed conformation upon binding of substrate or inhibitor molecules. These conformational changes have been suggested to follow an induced-fit mechanism in which the molecules first bind in the open conformation in those cases where binding in the closed conformation appears to be sterically obstructed such as for the HIV-1 protease. In this article, we present a general model for the catalysis and inhibition of enzymes with induced-fit binding mechanism. We derive general expressions that specify how the overall catalytic rate of the enzymes depends on the rates for binding, for the conformational changes, and for the chemical reaction. Based on these expressions, we analyze the effect of mutations that mainly shift the conformational equilibrium on catalysis and inhibition. If the overall catalytic rate is limited by product unbinding, we find that mutations that destabilize the closed conformation relative to the open conformation increase the catalytic rate in the presence of inhibitors by a factor exp(ΔΔG_C/RT) where ΔΔG_C is the mutation-induced shift of the free-energy difference between the conformations. This increase in the catalytic rate due to changes in the conformational equilibrium is independent of the inhibitor molecule and, thus, may help to understand how non-active-site mutations can contribute to the multi-drug-resistance that has been observed for the HIV-1 protease. A comparison to experimental data for the non-active-site mutation L90M of the HIV-1 protease indicates that the mutation slightly destabilizes the closed conformation of the enzyme. This article is part of a Special Issue entitled: The emerging dynamic view of proteins: Protein plasticity in allostery, evolution and self-assembly.     

Evaluation of triazolamers as active site inhibitors of HIV-1 protease

Proteases typically recognize their peptide substrates in extended conformations. General approaches for designing protease inhibitors often consist of peptidomimetics that feature this conformation. Herein we discuss a combination of computational and experimental studies to evaluate the potential of triazole-linked β-strand mimetics as inhibitors of HIV-1 protease activity.     

Reduced-bond tight-binding inhibitors of HIV-1 protease. Fine tuning of the enzyme subsite specificity.

Truncation of a peptide substrate in the N-terminus and replacement of its scissile amide bond with a non-cleavable reduced bond results in a potent inhibitor of HIV-1 protease. A series of such inhibitors has been synthesized, and S2-S3' subsites of the protease binding cleft mapped. The S2 pocket requires bulky Boc or PIV groups, large aromatic Phe residues are preferred in P1 and P1' and Glu in P2'. The S3' pocket prefers Phe over small Ala or Val. Introduction of a Glu residue into the P2' position yields a tight-binding inhibitor of HIV-1 protease, Boc-Phe-[CH2-NH]-Phe-Glu-Phe-OMe, with a subnanomolar inhibition constant. The relevant peptide derived from the same amino acid sequence binds to the protease with a Ki of 110 nM, thus still demonstrating a good fit of the amino acid residues into the protease binding pockets and also the importance of the flexibility of P1-P1' linkage for proper binding. A new type of peptide bond mimetic, N-hydroxylamine -CH2-N(OH)-, has been synthesized. Binding of hydroxylamino inhibitor of HIV-1 protease is further improved with respect to reduced-bond inhibitor.     

Cooperative fluctuations of unliganded and substrate-bound HIV-1 protease: a structure-based analysis on a variety of conformations from crystallography and molecular dynamics simulations

The dynamics of HIV-1 protease, both in unliganded and substrate-bound forms have been analyzed by using an analytical method, Gaussian network model (GNM). The method is applied to different conformations accessible to the protein backbone in the native state, observed in crystal structures and snapshots from fully atomistic molecular dynamics (MD) simulation trajectories. The modes of motion obtained from GNM on different conformations of HIV-1 protease are conserved throughout the MD simulations. The flaps and 40's loop of the unliganded HIV-1 protease structure are identified as the most mobile regions. However, in the liganded structure these flaps lose mobility, and terminal regions of the monomers become more flexible. Analysis of the fast modes shows that residues important for stability are in the same regions of all the structures examined. Among these, Gly86 appears to be a key residue for stability. The contribution of residues in the active site region and flaps to the stability is more pronounced in the substrate-bound form than in the unliganded form. The convergence of modes in GNM to similar regions of HIV-1 protease, regardless of the conformation of the protein, supports the robustness of GNM as a potentially useful and predictive tool.     

Role of hydroxyl group and R/S configuration of isostere in binding properties of HIV-1 protease inhibitors.

The crystal structure of the complex between human immunodeficiency virus type 1 (HIV-1) protease and a peptidomimetic inhibitor of ethyleneamine type has been refined to R factor of 0.178 with diffraction limit 2.5 A. The peptidomimetic inhibitor Boc-Phe-Psi[CH2CH2NH]-Phe-Glu-Phe-NH2 (denoted here as OE) contains the ethyleneamine replacement of the scissile peptide bond. The inhibitor lacks the hydroxyl group which is believed to mimic tetrahedral transition state of proteolytic reaction and thus is suspected to be necessary for good properties of peptidomimetic HIV-1 protease inhibitors. Despite the missing hydroxyl group the inhibition constant of OE is 1.53 nm and it remains in the nanomolar range also towards several available mutants of HIV-1 protease. The inhibitor was found in the active site of protease in an extended conformation with a unique hydrogen bond pattern different from hydroxyethylene and hydroxyethylamine inhibitors. The isostere nitrogen forms a hydrogen bond to one catalytic aspartate only. The other aspartate forms two weak hydrogen bridges to the ethylene group of the isostere. A comparison with other inhibitors of this series containing isostere hydroxyl group in R or S configuration shows different ways of accommodation of inhibitor in the active site. Special attention is devoted to intermolecular contacts between neighbouring dimers responsible for mutual protein adhesion and for a special conformation of Met46 and Phe53 side chains not expected for free protein in water solution.     

Conformational analysis of bacitracin A, a naturally occurring lariat

The proton nmr spectra of bacitracin A in H2O and DMSO-d6 have been assigned and the conformational behavior of the peptide in the two solvents has been compared. Although bacitracin A shows a conformational equilibrium between at least two conformations differing in the relative position of the cyclic and linear domains of the molecule, the spectra in water can be interpreted in terms of a preferred conformation in which the linear part is folded over the cyclic moiety and a turn is present around Ile(8)-DPhe(9).     

Insight into the structural similarity between HIV protease and secreted aspartic protease-2 and binding mode analysis of HIV-Candida albicans inhibitors

The analysis of the structural similarity between Candida albicans Sap2 and HIV-1 aspartic proteases by molecular modeling gave insight into the common requirements for inhibition of both targets. Structure superimposition of Sap2 and HIV-1 protease confirmed the similarity between their active sites and flap regions. HIV-1 protease inhibitors herein investigated can fit the active site of Sap2, adopting very similar ligand-backbone conformations. In particular, key anchoring sites consisting of Gly85 in Sap2 and Ile50 in HIV-1 protease, both belonging to their corresponding flap regions, were found as elements of a similar binding-mode interaction. The knowledge of the molecular basis for binding to both Sap2 and HIV-1 proteases may ultimately lead to the development of single inhibitor acting on both targets.     

Preferred conformations of cyclic Ac-Cys-Pro-Xaa-Cys-NHMe peptides: a model for chain reversal and active site of disulfide oxidoreductase

The conformational study on cyclic Ac-Cys-Pro-Xaa-Cys-NHMe (Ac-CPXC-NHMe; X=Ala, Val, Leu, Aib, Gly, His, Phe, Tyr, Asn and Ser) peptides has been carried out using the Empirical Conformational Energy Program for Peptides, version 3 (ECEPP/3) force field and the hydration shell model in the unhydrated and hydrated states. This work has been undertaken to investigate structural implications of the CPXC sequence as the chain reversal for the initiation of protein folding and as the motif for active site of disulfide oxidoreductases. The backbone conformation DAAA is commonly the most feasible for cyclic CPXC peptides in the hydrated state, which has a type I beta-turn at the Pro-Xaa sequence. The proline residue and the hydrogen bond between backbones of two cystines as well as the formation of disulfide bond appear to play a role in stabilizing this preferred conformation of cyclic CPXC peptides. However, the distributions of backbone conformations and beta-turns may indicate that the cyclic CPXC peptide seems to exist as an ensemble of beta-turns and coiled conformations in aqueous solution. The intrinsic stability of the cyclic CPXC motif itself for the active conformation seems to play a role in determining electrochemical properties of disulfide oxidoreductases.     

Influence of local interactions on protein structure. I. Conformational energy studies of N-acetyl-N′-methylamides of pro-X and X-pro dipeptides

Conformational energy calculations using an Empirical Conformational Energy Program for Peptides (ECEPP) were carried out on the N-acetyl-N′-methylamides of Pro-X, where X = Ala, Asn, Asp, Gly, Leu, Phe, Ser, and Val, and of X-Pro, where X = Ala, Asn, Gly, and Pro. The conformational energy was minimized from starting conformations which included all combinations of low-energy single-residue minima and several standard bend structures. It was found that almost all resulting minima are combinations of low-energy single-residue minima, suggesting that intra residue interactions predominate in determining conformation. The calculations also indicate, however, that inter residue interactions can be important. In addition, librational entropy was found to influence the relative stabilities of some minima. Because of the existence of 10–100 low-energy minima for each dipeptide, the normalized statistical weight of an individual minimum rarely exceeds 0.3, suggesting that these dipeptides have considerable conformational flexibility and exist as statistical ensembles of low-energy structures. The propensity of each dipeptide to form bend conformations was calculated, and the results were compared with available experimental data. It was found that bends are favored in Pro-X dipeptides because ?_Pro is fixed by the pyrrolidine ring in a conformation which is frequently found in bends, but that bends are not favored in X-Pro dipeptides because interactions between the X residue and the pyrrolidine ring restrict the X residue to conformations which are not usually found in bends.     

Mechanism of action of aspartic proteinases. V. Conformational characteristics of fragments of substrate-binding sites in rhizopuspepsin and HIV-1 proteinase

The conformational states of side chains of catalytic Asp residues in active sites of HIV-1 protease and rhizopuspepsin in the potential field of free enzymes were studied by using theoretical conformational analysis. Structural factors that stabilize the conformation of these residues in free enzymes were revealed. Methods of molecular mechanics were used to estimate the stabilization energy of the Met46-Phe53 labile fragments of HIV-1 protease in the potential field of their nearest surrounding amino acid residues for the conformations characteristic of the free protein and similar to that of the protein in enzyme-inhibitor complexes. In solution, the conformational state of the fragments of the free enzyme was concluded to be similar to that observed in the enzyme complex with the ligand and different from that determined by X-ray diffraction analysis. This difference was ascribed to the effect of crystal packing.     

Mechanism of action of aspartic proteases. III. Conformational characteristics of HIV-1 protease inhibitor JG-365]

A set of conformations was shown to be characteristic of the free-state spatial structure of substrate-like inhibitor JG-365 for aspartic protease from HIV-1. Among them, the lowest-energy conformations have a folded form of the peptide backbone. The inhibitor has a noncleavable hydroxyethylamine group with an additional chiral center in its structure. Our calculations showed that only the S-isomer of the inhibitor displays conformational characteristics that practically coincide with those of the native substrate for HIV-1 protease. One of the calculated conformations with a completely extended main chain and a relative energy of 9.5 kcal/mol very closely mimics the experimentally observed structure of the inhibitor in the enzyme-inhibitor complex. The realization of this structure is unlikely for a free inhibitor, because it has only a small number of interresidual noncovalent interactions in the extended conformation; these are presumably compensated for by intermolecular interactions at the active site of the enzyme.     

NMR and molecular dynamics studies of tachykinins: conformation of substance P fragment 4-11

The conformation of the C-terminal octapeptide fragment of Substance P (SP4-11, Pro-Gln-Gln-Phe-Phe-Gly-Leu-Met-NH2) has been investigated by 2D-NMR and MD methods. The octapeptide exists in a blend of conformations. The molecule seems to shuttle between conformations with gamma-bends either at Phe5 or Gly6 or Gln3 or Leu7 and between a nearly extended structure.     

Design, synthesis and structure-activity study of shorter hexa peptide analogues as HIV-1 protease inhibitors

Inhibition of HIV-1 protease enzyme can render the Human Immunodeficiency Virus (HIV-1) non-infectious in vitro. Previous studies have shown that several shorter peptides were discovered as HIV-1 protease inhibitors. In this context, a series of shorter synthetic hexapeptides, Leu-Leu-Glu-Tyr-Val-Xaa (Xaa=Phe, Met, Tyr and Trp), were designed. The synthesized hexa peptides were screened for their HIV-1 protease inhibition. These peptides showed moderately good HIV-1 protease inhibition when compared to acetyl pepstatin.     

Molecular mechanics (MM3(pi)) conformational analysis of molecules containing conjugated pi-electron fragments: Leucomycin-V

The conformations of the 16-membered macrolide antibiotic leucomycin-V (1) were studied with molecular mechanics. Leucomycin-V contains a conjugated pi-electron fragment and necessitates special treatment with the MM3(pi) modeling protocol. Comparison was made with results from the standard MM3 scheme. The CONFLEX conformational search procedure was used for finding low-energy conformations. The computed data are indicative for the existence of mainly one conformation of the macro-ring of 1 and minor participation of several others. Intramolecular hydrogen bonds play important roles for the preferred geometry of the macro-ring and the conformations of the side chains. The most probable macro-ring conformation of 1 is very similar to the preferred conformation of another 16-ring macrolide antibiotic, tylosin. The same order of conformational preference for 1 was estimated with the MM3 and the MM3(pi) methods. Surprisingly, when changing the chirality of the C(9) macro-ring atom of 1, the two methods produced different order of conformational preferences for the 9-epi form (2), as well as enhanced population of several clusters of conformations.