Computational anti-AIDS drug design based on the analysis of the specific interactions between immunophilins and the HIV-1 gp120 V3 loop. Application to the FK506-binding protein

The object of the study was to model the structural complex of the FK506-binding protein (FKBP) with the CRK peptide imitating the central region of the HIV-1 V3 loop, as well as to define the FKBP stretch giving rise to the binding site for V3 the synthetic copy of which, on the assumption of preserving the spatial peptide structure in the free state, can be considered as a promising applicant for the role of antiviral drug. To this end, the following successive steps were carried out: (i) the NMR-based conformational analysis of CRK was put into practice, and, in the light of the results derived, the best energy CRK structure meeting the requirements of the input NMR data was identified; (ii) molecular docking of the CRK structure with the X-ray FKBP conformation was implemented, and energy refining the simulated structural complex was realized; (iii) the matrix of distances between amino acids of the ligand and receptor was computed to specify the FKBP stretch keeping in touch with CRK followed by analyzing the types of interactions stabilizing the over-molecular ensemble; (iv) 3D structure of this stretch in the unbound status referred to as the FKBP peptide was predicted, and its collation with the X-ray conformation of the identical FKBP site was performed; (v) the potential energy function and its constituents were studied for the structural complex generated by molecular docking of the CRK molecule with the FKBP peptide; and (vi) from all evidence, the virtual FKBP-derived peptide was submitted to be utilized as a prospective structural framework in the anti-HIV-1 drug design. Summing up the results obtained, the following principal conclusion was drawn: a high affinity of the V3 loop peptide to the FKBP is based on the principle of "mirror similarity" that implies the near resemblance of 3D structures for the two individual fragments of the receptor and ligand, which, most likely, accounts for recognizing the immunophilin by V3 and determines the specificity of their efficacious interactions arising from the experimental observations.     

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.     

3D Structure Model of the Principal Neutralizing Epitope of Minnesota HIV-1 Isolate

Abstract

A hierarchical procedure, using a “bottom-up” strategy and combining (i) a probabilistic approach for estimating all possible starting structures, (ii) restrained molecular mechanics algorithms for preliminary selection of all energetically preferred conformers, as well as (iii) quantum chemical computations for refining their geometry, was used to study the structural properties of the HIV-MN neutralizing epitope in terms of NMR spectroscopy data. As a result, only one of initial structures matching the experimental and theoretical data was found to be well-ground for implementing the function of immunoreactive conformation of the virus immunogenic crown. The geometric parameters of this structure in water solution were shown to correspond to a double β-turn conformation similar to that revealed in crystal for synthetic molecules imitating the central region of the HIV-MN V3 loop. The following conclusion was drawn from the comparative analysis of simulated structure with the one computed previously: the HIV-MN immunogenic tip has some inherent conformational flexibility that manifests at the alterations of hexapeptide environment and leads to the structural transitions changing the local conformation of the stretch of interest but retaining its spatial main chain fold. As a matter of record, the high resolution 3D structure model for the HIV-MN principal neutralization site was constructed, and its geometric parameters were compared with the corresponding characteristics of conformers derived earlier for describing the conformational features of immunogenic tip of gp120 from Thailand HIV-1 isolate.

The results are discussed in the light of literature data on HIV-1 neutralizing epitope structure.     

3D structure model of the principal neutralizing epitope of Minnesota HIV-1 isolate

A hierarchical procedure, using a "bottom-up" strategy and combining (i). a probabilistic approach for estimating all possible starting structures, (ii). restrained molecular mechanics algorithms for preliminary selection of all energetically preferred conformers, as well as (iii). quantum chemical computations for refining their geometry, was used to study the structural properties of the HIV-MN neutralizing epitope in terms of NMR spectroscopy data. As a result, only one of initial structures matching the experimental and theoretical data was found to be well-ground for implementing the function of immunoreactive conformation of the virus immunogenic crown. The geometric parameters of this structure in water solution were shown to correspond to a double beta-turn conformation similar to that revealed in crystal for synthetic molecules imitating the central region of the HIV-MN V3 loop. The following conclusion was drawn from the comparative analysis of simulated structure with the one computed previously: the HIV-MN immunogenic tip has some inherent conformational flexibility that manifests at the alterations of hexapeptide environment and leads to the structural transitions changing the local conformation of the stretch of interest but retaining its spatial main chain fold. As a matter of record, the high resolution 3D structure model for the HIV-MN principal neutralization site was constructed, and its geometric parameters were compared with the corresponding characteristics of conformers derived earlier for describing the conformational features of immunogenic tip of gp120 from Thailand HIV-1 isolate.     

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.     

Recognition of the gluconeogenic enzyme,Pck1, via the Gid4 E3 ligase: An in silico perspective

Gluconeogenesis, the reverse process of glycolysis, is a favorable mechanism at conditions of glucose deprivation. Pck1 is a rate‐limiting gluconeogenic enzyme, where its deficiency or mutation contributes to serious clinical situations as neonatal hypoglycemia and liver failure. A recent report confirms that Pck1 is a target for proteasomal degradation through its proline residue at the penultimate position, recognized by Gid4 E3 ligase, but with a lack of informative structural details. In this study, we delineate the localized sequence motif, degron, that specifically interact with Gid4 ligase and unravel the binding mode of Pck1 to the Gid4 ligase by using molecular docking and molecular dynamics. The peptide/protein docking HPEPDOCK web server along with molecular dynamic simulations are applied to demonstrate the binding mode and interactions of a Pck1 wild type (SPSK) and mutant (K4V) with the recently solved structure of Gid4 ligase. Results unveil a distinct binding mode of the mutated peptide compared with the wild type despite having comparable binding affinities to Gid4. Moreover, the four‐residue peptide is found insufficient for Gid4 binding, while the seven‐residue peptide suffices for binding to Gid4. The amino acids S134, K135, and N137 in the loop L1 (between β1 and β2) of the Gid4 are essential for the stabilization of the seven‐residue peptide in the binding site of the ligase. The presence of Val⁴ instead of Lys⁴ smashes the H‐bonds that are formed between Lys⁴ and Gid4 in the wild type peptide, making the peptide prone to bind with the other side of the binding pocket (L4 loop of Gid4). The dynamics of Gid4 L3 loop is affected dramatically once K4V mutant Pck1 peptide is introduced. This opens the door to explore the mutation effects on the binding mode and smooth the path to target protein degradation by design competitive and non‐competitive inhibitors.     

Discovery of novel anti-HIV-1 agents based on a broadly neutralizing antibody against the envelope gp120 V3 loop: a computational study

A computer-aided search for novel anti-HIV-1 agents able to mimic the pharmacophore properties of broadly neutralizing antibody (bNAb) 3074 was carried out based on the analysis of X-ray complexes of this antibody Fab with the MN, UR29, and VI191 peptides from the V3 loop of the HIV envelope protein gp120. Using these empirical data, peptidomimetic candidates of bNAb 3074 were identified by a public, web-oriented virtual screening platform (pepMMsMIMIC) and models of these candidates bound to the above V3 peptides were generated by molecular docking. The docking calculations identified four molecules exhibiting a high affinity to all of the V3 peptides. These molecules were selected as the most probable peptidomimetics of bNAb 3074. Finally, the stability of the complexes of these molecules with the MN, UR29, and VI191 V3 peptides was estimated by molecular dynamics and free energy simulations. Specific binding to the V3 loop was accomplished primarily by π–π interactions between the aromatic rings of the peptidomimetics and the conserved Phe-20 and/or Tyr-21 of the V3 immunogenic crown. In a mechanism similar to that of bNAb 3074, these compounds were found to block the tip of the V3 loop forming its invariant structural motif that contains residues critical for cell tropism. Based on these findings, the compounds selected are considered as promising basic structures for the rational design of novel, potent, and broad-spectrum anti-HIV-1 therapeutics.     

Structural Analysis of the HIV-1 gp120 V3 Loop: Application to the HIV-Haiti Isolates

Abstract

The model describing the structure and conformational preferences of the HIV-Haiti V3 loop in the geometric spaces of Cartesian coordinates and dihedral angles was generated in terms of NMR spectroscopy data published in literature. To this end, the following successive steps were put into effect: (i) the NMR-based 3D structure for the HIV-Haiti V3 loop in water was built by computer modeling methods; (ii) the conformations of its irregular segments were analyzed and the secondary structure elements identified; and (iii) to reveal a common structural motifs in the HIV-Haiti V3 loop regardless of its environment variability, the simulated structure was collated with the one deciphered previously for the HIV-Haiti V3 loop in a water/trifluoroethanol (TFE) mixed solvent.

As a result, the HIV-Haiti V3 loop was found to offer the highly variable fragment of gp120 sensitive to its environment whose changes trigger the large-scale structural rearrangements, bringing in substantial altering the secondary and tertiary structures of this functionally important site of the virus envelope. In spite of this fact, over half of amino acid residues that reside, for the most part, in the functionally important regions of the gp120 protein and may present promising targets for AIDS drug researches, were shown to preserve their conformational states in the structures under review. In particular, the register of these amino acids holds Asn-25 that is critical for the virus binding with primary cell receptor CD4 as well as Arg-3 that is critical for utilization of CCR5 co-receptor and heparan sulfate proteoglycans. The conservative structural motif embracing one of the potential sites of the gp120 N-linked glycosylation was detected, which seems to be a promising target for the HIV-1 drug design.

The implications are discussed in conjunction with the literature data on the biological activity of the individual amino acids for the HIV-1 gp120 V3 loop.     

192 Computer-aided search for novel anti-HIV-1 agents presenting peptidomimetics of broadly neutralizing antibodies against the virus envelope GP120 V3 loop

Computer-aided search for novel anti-HIV-1 agents that are able to imitate the pharmacophore properties of the antigen-binding site of a broadly neutralizing mAb 3074 against the envelope gp120 V3 loop was carried out followed by evaluation of their potential inhibitory activity by molecular modeling. In doing so, the following problems were solved: (1) the mAb 3074 amino acid residues responsible for specific binding to the HIV-1 V3 loop were identified from the X-ray structures of this antibody Fab in complexes with the MN, UR29, and VI191 V3 peptides (Jiang et al., 2010); (2) using these data, 2039 possible mAb-3074 peptidomimetics were found by pepMMsMIMIC presenting a public, web-oriented virtual screening platform (Floris et al., 2011); (3) the complexes of these compounds with the above V3 peptides were built by molecular docking and, based on their analysis, the four molecules exhibiting a high affinity to V3 in the in silico studies were selected as the most probable peptidomimetics of mAb 3074 (Figure 1); and (4) stability of the complexes of these molecules with the MN, UR29, and VI191 V3 peptides was estimated by molecular dynamics and free energy simulations. As a result, a key role in specific binding of the selected compounds to the V3 loop was shown to belong to π-π interactions between their aromatic rings and the conserved Phe20 and/or Tyr21 of the V3 immunogenic crown. Similarly to mAb 3074, these compounds were found to block the tip of the V3 loop forming its invariant structural motif, which contains residues critical for cell tropism (Andrianov et al., 2011; Andrianov et al., 2012). In addition, the complexes of interest do not undergo significant changes within the molecular dynamics calculations, exhibiting the low values of free energy of their formation. In this context, the compounds given in Figure 1 are considered as the promising basic structures for the design of novel, potent, and broad anti-HIV-1 drugs.     

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.     

Modeling the intact form of the alpha 1-proteinase inhibitor

The structure of the intact form of the serpin alpha 1-proteinase inhibitor has been modeled based on the assumption that the central strand s4A of the six-stranded beta-sheet A of the cleaved inhibitor is not incorporated into the sheet of intact alpha 1-proteinase inhibitor. This strand was removed from its position in the center of the sheet by suitable rotations about the backbone dihedrals of Lys343 using molecular graphics. The resulting structure was then annealed using molecular dynamics (MD) while applying progressive distance restraints to the reactive peptide bond (Met358-Ser359) for 50 ps. During this time, the disrupted beta-sheet reformed to create a five-stranded beta-sheet with strands 3 and 5 in a parallel arrangement. This change and accompanying structural rearrangements are largely confirmed by the X-ray structure of plakalbumin, whose structure reflects the overall structure of intact serpins. The successful modeling experiment demonstrates the utility of MD for making gross structural predictions based on related structures. The binding loop of the intact form is modeled to allow docking with serine proteinases, in particular thrombin, which most highly constrains the possible conformations of the binding loop.     

Determination of structurally conservative amino acids of the HIV-1 protein gp120 V3 loop as promising targets for drug design by protein engineering approaches

Based on the published NMR spectroscopy data, three-dimensional structures of the HIV-1 gp120 protein V3 loop were obtained by computer modeling in the viral strains HIV-Haiti and HIV-MN. In both cases, the secondary structure elements and conformations of irregular stretches were determined for the fragment representing the principal antigenic determinant of the virus, as well as determinants of the cellular tropism and syncytium formation. Notwithstanding the high variability of the amino acid sequence of gp120 protein, more than 50% of the V3 loop residues retained their conformations in the different HIV-1 virions. The combined analysis of the findings and the literature data on the biological activity of the individual residues of the HIV-1 V3 loop resulted in identification of its structurally conservative amino acids, which seem to be promising targets for antiviral drug design by protein engineering approaches.     

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.     

Flexible protein-flexible ligand docking with disrupted velocity simulated annealing

By docking flexible balanol to a rigid model of protein kinase A (PKA), we found that a new simulated annealing protocol termed disrupted velocity simulated annealing (DIVE-SA) outperformed the replica-exchange method and the traditional simulated annealing method in identifying the correct docking pose. In this protocol, the atomic velocities were reassigned periodically to encourage the system to sample a large conformational space. We also found that scaling potential energy surface to reduce structural transition barriers could further facilitate docking. The DIVE-SA method was then evaluated on its ability to perform flexible ligand-flexible protein docking of three ligands (balanol, a balanol analog, and ATP) to PKA. To reduce computational time and to avoid possible unphysical structural changes resulting from the use of nonoptimal force fields, a soft restrain was applied to keep the root-mean-square-deviation (RMSD) between instantaneous protein structures and a chosen reference structure small. Because the restrain was applied to the overall RMSD rather than to individual atoms, a protein could still experience relatively large conformational changes during docking. To examine the impact of applying such a restrain on docking, we constructed two semi-flexible protein models by choosing two different crystal structures as reference. Both the balanol analog and ATP were able to dock to either one of these semi-flexible protein models. On the other hand, balanol could only dock well to one of them. Further analysis indicated that the restrain on the glycine-rich loop was too strong, preventing it to adjust its structure to accommodate balanol in the binding pocket of PKA. Removing the restrain on the glycine-rich loop resulted in much better docking poses. This finding demonstrates the important role that the flexibility of the glycine-rich loop play in accepting different ligands and should profitably not be restrained in molecular docking so that more diverse ligands can be studied.     

Systematic search for putative new domain families in Mycoplasma gallisepticum genome

Background

Stem cell factor (SCF) receptor c-Kit is recognized as a key signaling molecule, which transduces signals for the proliferation, differentiation and survival of stem cells. Binding of SCF to its receptor triggers transactivation, leading to the recruitment of kinases and phosphatases to the docking platforms of c-Kit catalytic domain. Tyrosine phosphatase-1 (Shp-1) deactivates/attenuates 'Kit' kinase activity. Whereas, Asp816Val mutation in the Kit activation loop transforms kinase domain to a constitutively activated state (switch off-to-on state), in a ligand-independent manner. This phenomenon completely abrogates negative regulation of Shp-1. To predict the possible molecular basis of interaction between c-Kit and Shp-1, we have performed an in silico protein-protein docking study between crystal structure of activated c-Kit (phosphorylated c-Kit) and full length crystal structure of Shp-2, a close structural counterpart of Shp-1.

Findings

Study revealed a stretch of conserved amino acids (Lys818 to Ser821) in the Kit activation domain, which makes decisive H-bonds with N-sh2 and phosphotyrosine binding pocket residues of the phosphatase. These H-bonds may impose an inhibitory steric hindrance to the catalytic domain of c-Kit, there by blocking further interaction of the activation loop molecules with incoming kinases. We have also predicted a phosphotyrosine binding pocket in SH2 domains of Shp-1, which is found to be predominantly closer to a catalytic groove like structure in c-Kit kinase domain.

Conclusions

This study predicts that crucial hydrogen bonding between N-sh2 domain of Shp-1 and Kit activation loop can modulate the negative regulation of c-Kit kinase by Shp-1. Thus, this finding is expected to play a significant role in designing suitable gain-of-function c-Kit mutants for inducing conditional proliferation of hematopoietic stem cells.     

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.     

Identification of the potential regions of Epap-1 that interacts with V3 loop of HIV-1 gp120

Early pregnancy associated protein-1 (Epap-1), a 90 kDa glycoprotein present in first trimester placental tissue, inhibits HIV-1 entry through interaction with HIV-1 gp120 at V3 and C5 regions. In the present study, we have identified the specific 32 mer region of Epap-1 that can interact with V3 loop. This was achieved by docking between Epap-1 molecular model and gp120 and studying the interaction of peptides with gp120 in vitro. Out of four peptides analyzed, two peptides (P-2 and P-3) showed significant interaction with V3 domain (N = 8; N = 7) of gp120. In the studies conducted using soluble gp120 and virus, peptide P-2 has shown conserved interaction at V3 loop regions recognized by 257D and F425 antibodies and higher anti-viral activity. Also, P-2 inhibited cell fusion mediated dye transfer between gp120 expressing HL2/3 and CD4 expressing Sup T1 cells suggesting its inhibition of viral entry, which is further confirmed by its action on HIV infection mediated by Tat activated beta gal expression in TZM-bl cells. Further optimization of P-2 peptide showed that the anti-viral activity and gp120 interaction residues lie in the N-terminal region of the peptide. These results together suggest that P-2 inhibits viral entry through specific interaction at V3 loop region.     

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.     

Identification of Contact Sites between Ankyrin and Band 3 in the Human Erythrocyte Membrane

The red cell membrane is stabilized by a spectrin/actin-based cortical cytoskeleton connected to the phospholipid bilayer via multiple protein bridges. By virtue of its interaction with ankyrin and adducin, the anion transporter, band 3 (AE1), contributes prominently to these bridges. In a previous study, we demonstrated that an exposed loop comprising residues 175-185 of the cytoplasmic domain of band 3 (cdB3) constitutes a critical docking site for ankyrin on band 3. In this paper, we demonstrate that an adjacent loop, comprising residues 63-73 of cdB3, is also essential for ankyrin binding. Data that support this hypothesis include the following. (1) Deletion or mutation of residues within the latter loop abrogates ankyrin binding without affecting cdB3 structure or its other functions. (2) Association of cdB3 with ankyrin is inhibited by competition with the loop peptide. (3) Resealing of the loop peptide into erythrocyte ghosts alters membrane morphology and stability. To characterize cdB3-ankyrin interaction further, we identified their interfacial contact sites using molecular docking software and the crystal structures of D(3)D(4)-ankyrin and cdB3. The best fit for the interaction reveals multiple salt bridges and hydrophobic contacts between the two proteins. The most important ion pair interactions are (i) cdB3 K69-ankyrin E645, (ii) cdB3 E72-ankyrin K611, and (iii) cdB3 D183-ankyrin N601 and Q634. Mutation of these four residues on ankyrin yielded an ankyrin with a native CD spectrum but little or no affinity for cdB3. These data define the docking interface between cdB3 and ankyrin in greater detail.     

Discovery of potential sclerostin inhibitors from plants with loop2 region of sclerostin inhibition by interacting with residues outside Pro-Asn-Ala-Ile-Gly motif

Abstract

Sclerostin, an antagonist of the Wnt/β-catenin signaling pathway, was discovered as a potential therapeutic target for stimulating bone formation in osteoporosis. In this study, molecular docking was employed to predict the binding of 29 herbal compounds, which were reported as bone formation stimulators, to the loop2 region of sclerostin. Then, the 50 ns molecular dynamics (MD) simulation of the complexes between sclerostin and the top 10 hits obtained from molecular docking were carried out. Root mean square deviations (RMSDs) analysis of MD trajectories pointed out that all ligands-complexes remain stable throughout the duration of MD simulations. In addition, the molecular mechanics/generalized born surface area (MM/GBSA) binding free energy and energy decomposition analyses were determined. The results here suggested that baicalin is the most promising inhibitor of sclerostin. Interestingly, baicalin binds to sclerostin via the hydrophobic interaction with the amino acid residues on loop2 region but outside the Pro-Asn-Ala-Ile-Gly (PNAIG) motif, particularly the Arg-Gly-Lys-Trp-Trp-Arg (RGKWWR) motif. This finding could be a novel strategy for developing new sclerostin inhibitors in the future.

Communicated by Ramaswamy H. Sarma