Immunophilins and HIV-1 V3 Loop For Structure-Based Anti-AIDS Drug Design

Abstract

The model of the structural complex of cyclophilin A (CycA) belonging to the immunophilins family with the HIV-MN gpl20 V3 loop was generated, and the computer-aided design of the immunophilin-derived peptide able to mask the biologically crucial V3 segments was implemented.

To this end, the following problems were solved: (i) the NMR-based conformational analysis of the HIV-MN V3 loop was put into effect, and its low energy structure fitting the input experimental observations was determined; (ii) molecular docking of this V3 structure with the X-ray conformation of CycA was carried out, and the energy refining the simulated structural complex was performed; (iii) the matrix of inter-atomic distances for the amino acids of the molecules forming part of the built over-molecular ensemble was computed, the types of interactions responsible for its stabilization were analyzed, and the CycA stretch, which accounts for the binding to V3, was identified; (iv) the most probable 3D structure for this stretch in the unbound state was predicted, and its collation with the X-ray structure for the corresponding site of CycA was performed; (v) the potential energy function and its constituents were studied for the structural complex generated by molecular docking of the V3 loop with the CycA peptide offering the virtual molecule that imitates the CycA segment, making a key contribution to the interactions of the native protein with the HIV-1 principal neutralizing determinant; (vi) as a result of the studies above, the designed molecule was shown to be capable of the efficacious blockading the functionally crucial V3 sites; and (vii) based on the joint analysis of the evidence obtained previously and in the present study, the composition of the peptide cocktail presenting the promising anti-AIDS pharmacological substance was developed.

The molecules simulated here by molecular modeling methods may become the first representatives of a new class of the chemical compounds (immunophilin-derived peptides) offering the looking-forward basic structures for the design of efficacious and safe antiviral agents.     

The binding of a glycoprotein 120 V3 loop peptide to HIV-1 neutralizing antibodies. Structural implications

The structural and antigenic properties of a peptide ("CRK") derived from the V3 loop of HIV-1 gp120 protein were studied using NMR and SPR techniques. The sequence of CRK corresponds to the central portion of the V3 loop containing the highly conserved "GPGR" residue sequence. Although the biological significance of this conserved sequence is unknown, the adoption of conserved secondary structure (type II beta-turn) in this region has been proposed. The tendency of CRK (while free or conjugated to protein), to adopt such structure and the influence of such structure upon CRK antigenicity were investigated by NMR and SPR, respectively. Regardless of conjugation, CRK is conformationally averaged in solution but a weak tendency of the CRK "GPGR" residues to adopt a beta-turn conformation was observed after conjugation. The influence of GPGR structure upon CRK antigenicity was investigated by measuring the affinities of two cognate antibodies: "5023A" and "5025A," for CRK, protein-conjugated CRK and gp120 protein. Each antibody bound to all the antigens with nearly the same affinity. From these data, it appears that: (a) antibody binding most likely involves an induced fit of the peptide and (b) the gp120 V3 loop is probably conformationally heterogeneous. Since 5023A and 5025A are HIV-1 neutralizing antibodies, neutralization in these cases appears to be independent of adopted GPGR beta-turn structure.     

Rapid protein-ligand docking using soft modes from molecular dynamics simulations to account for protein deformability: binding of FK506 to FKBP

Most current docking methods to identify possible ligands and putative binding sites on a receptor molecule assume a rigid receptor structure to allow virtual screening of large ligand databases. However, binding of a ligand can lead to changes in the receptor protein conformation that are sterically necessary to accommodate a bound ligand. An approach is presented that allows relaxation of the protein conformation in precalculated soft flexible degrees of freedom during ligand-receptor docking. For the immunosuppressant FK506-binding protein FKBP, the soft flexible modes are extracted as principal components of motion from a molecular dynamics simulation. A simple penalty function for deformations in the soft flexible mode is used to limit receptor protein deformations during docking that avoids a costly recalculation of the receptor energy by summing over all receptor atom pairs at each step. Rigid docking of the FK506 ligand binding to an unbound FKBP conformation failed to identify a geometry close to experiment as favorable binding site. In contrast, inclusion of the flexible soft modes during systematic docking runs selected a binding geometry close to experiment as lowest energy conformation. This has been achieved at a modest increase of computational cost compared to rigid docking. The approach could provide a computationally efficient way to approximately account for receptor flexibility during docking of large numbers of putative ligands and putative docking geometries.     

New general approach for determining the solution structure of a ligand bound weakly to a receptor: structure of a fibrinogen Aalpha-like peptide bound to thrombin (S195A) obtained using NOE distance constraints and an ECEPP/3 flexible docking program.

A new approach incorporating flexible docking simulations and NMR data is presented for calculating the bound conformation of a ligand that interacts weakly with an enzyme. This approach consists of sampling directly the conformation of a flexible ligand inside a receptor active site containing surrounding flexible loops. To make this sampling efficient, a ligand-growing procedure has been adopted. Optimization of the ECEPP/3-plus-NOE constraint function is carried out by using a collective variable Monte Carlo minimization technique. Numerous energy minimizations are made possible for such a large system by using a Bezier splines energy grid technique. This new flexible docking approach was applied to determine the structure of a fibrinogen Aalpha-like peptide (7DFLAEGGGVRGPRV20) bound to an active site mutant of thrombin [thrombin(S195A)]. Structure calculations of the bound ligand, using 2D-transferred NOESY distance constraints in the DIANA program, showed that the N-terminal portion of the peptide (D7-R16) involves a chain reversal, whereas the C-terminal portion (G17-V20) adopts a fold that exists in several different orientations. In addition, the ECEPP/3 flexible docking package was used to assess the conformational variability of the ligand and surrounding 60D-insertion loop of thrombin. Amino acid residues (17-20) of the peptide interact with a region of the enzyme that exhibits broad specificity, with a preferred direction between the 60D-insertion loop and Pro37 of thrombin.     

Correlation between biological activity and binding energy in systems of integrin with cyclic RGD-containing binders: a QM/MM molecular dynamics study

We here report a combined quantum mechanical/molecular mechanical (QM/MM) molecular dynamics (MD) study on the binding interactions between the α(V)β(3) integrin and eight cyclic arginine-glycine-aspartate (RGD) containing peptides. The initial conformation of each peptide within the binding site of the integrin was determined by docking the ligand to the reactive site of the integrin crystal structure with the aid of docking software FRED. The subsequent QM/MM MD simulations of the complex structures show that these eight cyclic RGD-peptides have a generally similar interaction mode with the binding site of the integrin to the cyclo(RGDf-N[M]V) analog found in the crystal structure. Still, there are subtle differences in the interactions of peptide ligands with the integrin, which contribute to the different inhibition activities. The averaged QM/MM protein-ligand interaction energy (IE) is remarkably correlated to the biological activity of the ligand. The IE, as well as a three-variable model which is somewhat interpretable, thus can be used to predict the bioactivity of a new ligand quantitatively, at least within a family of analogs. The present study establishes a helpful protocol for advancing lead compounds to potent inhibitors.     

Determining the invariant structure elements of the HIV-1 variable V3 loops: insight into the HIV-MN and HIV-Haiti isolates

The computer approaches that combined the 3D protein structure modeling with the mathematical statistics methods were used to compute the NMR-based 3D structures of the HIV-1 gp120 V3 loop for the HIV-MN and HIV-Haiti isolates in water as well as to compare their conformational characteristics with the purpose of determining the structure elements common for the two virus modifications. As a result, the variability of the amino acid sequence was found to stimulate the considerable structural rearrangements of the V3 loop. However, despite this fact, one functionally crucial stretch of V3 and a greater portion of its residues were shown to preserve the conformations in the viral strains of interest. To reveal the structural motifs and individual amino acids giving rise to the close conformations in the HIV-MN and HIV-Haiti V3 loops regardless of the sequence and environment variability, the simulated structures were collated with those deciphered previously in terms of NMR data in a water/trifluoroethanol mixed solvent. The structure elements and single residues of V3 residing in its biologically significant sites and keeping the conformations in all of the cases at question are considered to be the promising targets for anti-AIDS drugs studies. In this context, the structurally inflexible motifs of V3 presenting the weak units in the virus protection system may be utilized as the most convenient landing-places for molecular docking of the V3 loop and ligand structures followed by selecting chemical compounds suitable as a basis for the design of safe and effective antiviral agents.     

Structure comparison of native and mutant human recombinant FKBP12 complexes with the immunosuppressant drug FK506 (tacrolimus).

The consequences of site-directed mutagenesis experiments are often anticipated by empirical rules regarding the expected effects of a given amino acid substitution. Here, we examine the effects of "conservative" and "nonconservative" substitutions on the X-ray crystal structures of human recombinant FKBP12 mutants in complex with the immunosuppressant drug FK506 (tacrolimus). R42K and R42I mutant complexes show 110-fold and 180-fold decreased calcineurin (CN) inhibition, respectively, versus the native complex, yet retain full peptidyl prolyl isomerase (PPIase) activity, FK506 binding, and FK506-mediated PPIase inhibition. Interestingly, the structure of the R42I mutant complex is better conserved than that of the R42K mutant complex when compared to the native complex structure, within both the FKBP12 protein and FK506 ligand regions of the complexes, and with respect to temperature factors and RMS coordinate differences. This is due to compensatory interactions mediated by two newly ordered water molecules in the R42I complex structure, molecules that act as surrogates for the missing arginine guanidino nitrogens of R42. The absence of such surrogate solvent interactions in the R42K complex leads to some disorder in the so-called "40s loop" that encompasses the substituent. One rationalization proposed for the observed loss in CN inhibition in these R42 mutant complexes invokes indirect effects leading to a misorientation of FKBP12 and FK506 structural elements that normally interact with calcineurin. Our results with the structure of the R42I complex in particular suggest that the observed loss of CN inhibition might also be explained by the loss of a specific R42-mediated interaction with CN that cannot be mimicked effectively by the solvent molecules that otherwise stabilize the conformation of the 40s loop in that structure.     

Solution structure of a neurotrophic ligand bound to FKBP12 and its effects on protein dynamics.

The structure of a recently reported neurotrophic ligand, 3-(3-pyridyl)-1-propyl(2S)-1-(3,3-dimethyl-1, 2-dioxopentyl)-2-pyrrolidinecarboxylate, in complex with FKBP12 was determined using heteronuclear NMR spectroscopy. The inhibitor exhibits a binding mode analogous to that observed for the macrocycle FK506, used widely as an immunosuppressant, with the prolyl ring replacing the pipecolyl moiety and the amide bond in a trans conformation. However, fewer favourable protein-ligand interactions are detected in the structure of the complex, suggesting weaker binding compared with the immunosuppressant drug. Indeed, a micromolar dissociation constant was estimated from the NMR ligand titration profile, in contrast to the previously published nanomolar inhibition activity. Although the inhibitor possesses a remarkable structural simplicity with respect to FK506, 15N relaxation studies show that it induces similar effects on the protein dynamics, stabilizing the conformation of solvent-exposed residues which are important for mediating the interaction of immunophilin/ligand complexes with molecular targets and potentially for the transmission of the neurotrophic action of FKBP12 inhibitors.     

Predicting structural effects in HIV-1 protease mutant complexes with flexible ligand docking and protein side-chain optimization

We present a computational approach for predicting structures of ligand-protein complexes and analyzing binding energy landscapes that combines Monte Carlo simulated annealing technique to determine the ligand bound conformation with the dead-end elimination algorithm for side-chain optimization of the protein active site residues. Flexible ligand docking and optimization of mobile protein side-chains have been performed to predict structural effects in the V32I/I47V/V82I HIV-1 protease mutant bound with the SB203386 ligand and in the V82A HIV-1 protease mutant bound with the A77003 ligand. The computational structure predictions are consistent with the crystal structures of these ligand-protein complexes. The emerging relationships between ligand docking and side-chain optimization of the active site residues are rationalized based on the analysis of the ligand-protein binding energy landscape. Proteins 33:295–310, 1998. © 1998 Wiley-Liss, Inc.     

Study on conformational homology of the HIV-1 gp120 protein V3 loop. Structural analysis of the HIV-RF and HIV-thailand viral strains

Based on NMR spectroscopy data, conformation of the HIV-RF gp120 protein V3 loop giving rise to the virus principal neutralizing determinant and also determinants of cell tropism and syncytium formation was calculated by computer modeling approaches. Elements of the HIV-RF V3 loop secondary structure and conformational states of its irregular stretches were determined. The calculated structure was compared with the conformation of the homologous stretch of the HIV-Thailand protein gp120 V3 loop, and structural elements preserved in the two viral strains were identified. Conservative elements of the HIV-1 V3 loop structure are considered to be promising targets for deriving chemically modified forms of this loop with the enhanced immunogenicity and cross-reactivity of neutralizing antibodies and also for creation of effective antiviral drugs on this base.     

Computer modeling of promising inhibitors of the HIV-1 subtype A replication as a framework for the rational anti-AIDS drug design

The model of the structural complex of cyclophilin B belonging to the immunophilins family with the HIV-1 subtype A V3 loop presenting the principal neutralizing determinant of the virus gp120 envelope protein as well as determinants of cell tropism and syncytium formation was generated by molecular docking methods. Based on the conformational and energy characteristics of the built complex, computer-aided design of the polypeptide able to block effectively the functionally crucial V3 segments was implemented. Analysis of the results obtained in this study and literature data suggests that the generated molecule represents a promising pharmacological substance, which may be used as the basis structure for realization of the protein engineering projects aimed to develop effective drugs for anti-AIDS therapy.     

Local structural properties of the V3 loop of Thailand HIV-1 isolate

The model of locally accurate conformation for the HIV-Thailand principal neutralizing determinant (PND) located within the V3 loop of the virus envelope protein gp120 was built in terms of NMR spectroscopy data. To this end, the NMR-based conformational analysis of synthetic molecule representing the peptide copy of the fragment under study was carried out using the published sequential d connectivity data and values of spin-spin coupling constants. As a result, (i) the local structure for the V3 loop from Thailand isolate was determined, (ii) the conformations of its irregular segments were analyzed, and the secondary structure elements identified, (iii) the ensemble of conformers matching the experimental and theoretical data was derived for the stretch forming the neutralizing epitope of the HIV-Thailand PND, (iv) to estimate the probability of realizing each of these conformers in solution, the results obtained were collated with the X-ray data for corresponding segments in synthetic molecules imitating the central region of the HIV-MN PND as well as for homologous segments 39-44 in Bence-Jonce REI protein (BJRP), 41-46 in immunoglobulin lambda (Ig lambda), and 50-55 in beta-chain of horse hemoglobin (HH), (v) to find the conserved structural motifs inside diverse HIV-1 isolates, the structure determined was compared with the one derived earlier for the HIV-MN PND from NMR spectroscopy data, (vi) on the basis of all data obtained, the 3D structure model describing the set of biologically relevant conformations, which may present different antigenic determinants to the immune system in various HIV-1 isolates, was proposed for the immunogenic crown of the V3 loop. The results obtained are discussed in conjunction with the data on the structure for the HIV-1 PND reported in literature.     

The third-dimensional structure of the complex between an Fv antibody fragment and an analogue of the main immunogenic region of the acetylcholine receptor: a combined two-dimensional NMR, homology, and molecular modeling approach

Binding of autoantibodies to the acetylcholine receptor (AChR) plays a major role in the autoimmune disease Myasthenia gravis (MG). In this paper, we propose a structure model of a putative immunocomplex that gives rise to the reduction of functional AChR molecules during the course of MG. The model complex consists of the [G(70), Nle(76)] decapeptide analogue of the main immunogenic region (MIR), representing the major antigenic epitope of AChR, and the single chain Fv fragment of monoclonal antibody 198, a potent MG autoantibody. The structure of the complexed decapeptide antigen [G(70), Nle(76)]MIR was determined using two-dimensional nmr, whereas the antibody structure was derived by means of homology modeling. The final complex was constructed using calculational docking and molecular dynamics. We termed this approach "directed modeling," since the known peptide structure directs the prestructured antibody binding site to its final conformation. The independently derived structures of the peptide antigen and antibody binding site already showed a high degree of surface complementarity after the initial docking calculation, during which the peptide was conformationally restrained. The docking routine was a soft algorithm, applying a combination of Monte Carlo simulation and energy minimization. The observed shape complementarity in the docking process suggested that the structure assessments already led to anti-idiotypic conformations of peptide antigen and antibody fragment. Refinement of the complex by dynamic simulation yielded improved surface adaptation by small rearrangements within antibody and antigen. The complex presented herein was analyzed in terms of antibody-antigen interactions, properties of contacting surfaces, and segmental mobility. The structural requirements for AChR complexation by autoantibodies were explored and compared with experimental data from alanine scans of the MIR peptides. The analysis revealed that the N-terminal loop of the peptide structure, which is indispensable for antibody recognition, aligns three hydrophobic groups in a favorable arrangement leading to the burial of 40% of the peptide surface in the binding cleft upon complexation. These data should be valuable in the rational design of an Fv mutant with much improved affinity for the MIR and AChR to be used in therapeutic approaches in MG.     

Solid-state NMR evidence for an antibody-dependent conformation of the V3 loop of HIV-1 gp120

Solid-state NMR measurements have been carried out on frozen solutions of the complex of a 24-residue peptide derived from the third variable (V3) loop of the HIV-1 envelope glycoprotein gp120 bound to the Fab fragment of an anti-gp120 antibody. The measurements place strong constraints on the conformation of the conserved central GPGR motif of the V3 loop in the antibody-bound state. In combination with earlier crystal structures of V3 peptide-antibody complexes and existing data on the cross-reactivity of the antibodies, the solid-state NMR measurements suggest that the Gly-Pro-Gly-Arg (GPGR) motif adopts an antibody-dependent conformation in the bound state and may be conformationally heterogeneous in unbound, full-length gp120. These measurements are the first application of solid-state NMR methods in a structural study of a peptide-protein complex.     

Modeling and docking of the three-dimensional structure of the human melanocortin 4 receptor

A three-dimensional structure of the human melanocortin 4 receptor (hMC4R) is constructed in this study using a computer-aided molecular modeling approach. Human melanocortin 4 receptor is a G Protein-Coupled Receptor (GPCR). We structurally aligned transmembrane helices with bovine rhodopsin transmembrane domains, simulated both intracellular and extracellular loop domains on homologous loop regions in other proteins of known 3D structure and modeled the C terminus on the corresponding part of bovine rhodopsin. Then tandem minimization and dynamics calculations were run to refine the crude structure. The simulative model was tested by docking with a triplet peptide (RFF) ligand. It was found that the ligand is located among transmembrane regions TM3, TM4, TM5, and TM6 of hMC4R. In consistence with mutational and biochemical data, binding site is mainly formed as a hydrophobic and negatively charged pocket. The model constructed here might provide a structural framework for making rational predictions in relevant fields.     

Phosphorylation-dependent structure of ATF4 peptides derived from a human ATF4 protein, a member of the family of transcription factors

ATF4 plays a crucial role in the cellular response to stress and the F-box protein β-TrCP, the receptor component of the SCF E3 ubiquitin ligase responsible for ATF4 degradation by the proteasome, binds to ATF4, and controls its stability. Association between the two proteins depends on ATF4 phosphorylation of serine residues 219 and 224 present in the context of DpSGXXXpS, which is similar but not identical to the DpSGXXpS motif found in most other substrates of β-TrCP. We used NMR spectroscopy to analyze the structure of the 23P-ATF4 peptide. The 3D structure of the ligand was determined on the basis of NOESY restraints that provide an hairpin loop structure. In contrast, no ordered structure was observed in the NMR experiments for the nonphosphorylated 23-ATF4 in solution. This structural study provides information, which could be used to study the β-TrCP receptor–ligand interaction in docking procedure. Docking studies showed that the binding epitope of the ligand, is represented by the DpSGIXXpSXE motif. 23P-ATF4 peptide fits the binding pocket of protein β-TrCP very well, considering that the DpSGIXXpSXE motif adopts an S-turning conformation contrary to the extended DpSGXXpS motif in the other known β-TrCP ligands.     

Three-dimensional model of the human urotensin-II receptor: docking of human urotensin-II and nonpeptide antagonists in the binding site and comparison with an antagonist pharmacophore model

Human urotensin-II (hU-II) is a cyclic peptide that plays a central role in cardiovascular homeostasis and is considered to be the most potent mammalian vasoconstrictor identified to date. It is a natural ligand of the human urotensin-II (hUT-II) receptor, a member of the family of rhodopsin-like G-protein-coupled receptors. To understand the molecular interactions of hU-II and certain antagonists with the hUT-II receptor, a model of the hUT-II receptor in an active conformation with all its connecting loops was constructed by homology modeling. The initial model was placed in a pre-equilibrated lipid bilayer and re-equilibrated by several procedures of energy minimization and molecular dynamics simulations. Docking studies were performed for hU-II and for a series of nonpeptide hUT-II receptor antagonists in the active site of the modeled receptor structure. Results of the hU-II docking study are in agreement with our previous work and with experimental data showing the contribution of the extracellular loops II and III to ligand recognition. The docking of hU-II nonpeptide antagonists allows identification of key molecular interactions and confirms a previously reported hU-II antagonist pharmacophore model. The results of the present studies will be used in structure-based drug design for developing novel antagonists for the hUT-II receptor.     

Alternative conformations of HIV-1 V3 loops mimic beta hairpins in chemokines,suggesting a mechanism for coreceptor selectivity

The V3 loop of the HIV-1 envelope glycoprotein gp120 is involved in binding to the CCR5 and CXCR4 coreceptors. The structure of an HIV-1(MN) V3 peptide bound to the Fv of the broadly neutralizing human monoclonal antibody 447-52D was solved by NMR and found to be a beta hairpin. This structure of V3(MN) was found to have conformation and sequence similarities to beta hairpins in CD8 and CCR5 ligands MIP-1alpha, MIP-1beta, and RANTES and differed from the beta hairpin of a V3(IIIB) peptide bound to the strain-specific murine anti-gp120(IIIB) antibody 0.5beta. In contrast to the structure of the bound V3(MN) peptide, the V3(IIIB) peptide resembles a beta hairpin in SDF-1, a CXCR4 ligand. These data suggest that the 447-52D-bound V3(MN) and the 0.5beta-bound V3(IIIB) structures represent alternative V3 conformations responsible for selective interactions with CCR5 and CXCR4, respectively.     

Molecular docking of thiamine reveals similarity in binding properties between the prion protein and other thiamine-binding proteins

Prion-induced diseases are a global health concern. The lack of effective therapy and 100 % mortality rates for such diseases have made the prion protein an important target for drug discovery. Previous NMR experimental work revealed that thiamine and its derivatives bind the prion protein in a pocket near the N-terminal loop of helix 1, and conserved intermolecular interactions were noted between thiamine and other thiamine-binding proteins. Furthermore, water-mediated interactions were observed in all of the X-ray crystallographic structures of thiamine-binding proteins, but were not observed in the thiamine–prion NMR study. To better understand the potential role of water in thiamine–prion binding, a docking study was employed using structural X-ray solvent. Before energy minimization, docked thiamine assumed a “V” shape similar to some of the known thiamine-dependent proteins. Following minimization with NMR-derived restraints, the “F” conformation was observed. Our findings confirmed that water is involved in ligand stabilization and phosphate group interaction. The resulting refined structure of thiamine bound to the prion protein allowed the 4-aminopyrimidine ring of thiamine to π-stack with Tyr150, and facilitated hydrogen bonding between Asp147 and the amino group of 4-aminopyrimidine. Investigation of the π-stacking interaction through mutation of the tyrosine residue further revealed its importance in ligand placement. The resulting refined structure is in good agreement with previous experimental restraints, and is consistent with the pharmacophore model of thiamine-binding proteins.     

NMR structure of the thromboxane A2 receptor ligand recognition pocket.

To overcome the difficulty of characterizing the structures of the extracellular loops (eLPs) of G protein-coupled receptors (GPCRs) other than rhodopsin, we have explored a strategy to generate a three-dimensional structural model for a GPCR, the thromboxane A(2) receptor. This three-dimensional structure was completed by the assembly of the NMR structures of the computation-guided constrained peptides that mimicked the extracellular loops and connected to the conserved seven transmembrane domains. The NMR structure-based model reveals the structural features of the eLPs, in which the second extracellular loop (eLP(2)) and the disulfide bond between the first extracellular loop (eLP(1)) and eLP(2) play a major role in forming the ligand recognition pocket. The eLP(2) conformation is dynamic and regulated by the oxidation and reduction of the disulfide bond, which affects ligand docking in the initial recognition. The reduced form of the thromboxane A(2) receptor experienced a decrease in ligand binding activity due to the rearrangement of the eLP(2) conformation. The ligand-bound receptor was, however, resistant to the reduction inactivation because the ligand covered the disulfide bond and stabilized the eLP(2) conformation. This molecular mechanism of ligand recognition is the first that may be applied to other prostanoid receptors and other GPCRs.