Computational analysis of binding of P1 variants to trypsin

The binding of P1 variants of bovine pancreatic trypsin inhibitor (BPTI) to trypsin has been investigated by means of molecular dynamics simulations. The specific interaction formed between the amino acid at the primary binding (P1) position of the binding loop of BPTI and the specificity pocket of trypsin was estimated by use of the linear interaction energy (LIE) method. Calculations for 13 of the naturally occurring amino acids at the P1 position were carried out, and the results obtained were found to correlate well with the experimental binding free energies. The LIE calculations rank the majority of the 13 variants correctly according to the experimental association energies and the mean error between calculated and experimental binding free energies is only 0.38 kcal/mole, excluding the Glu and Asp variants, which are associated with some uncertainties regarding protonation and the possible presence of counter-ions. The three-dimensional structures of the complex with three of the P1 variants (Asn, Tyr, and Ser) included in this study have not at present been solved by any experimental techniques and, therefore, were modeled on the basis of experimental data from P1 variants of similar size. Average structures were calculated from the MD simulations, from which specific interactions explaining the broad variation in association energies were identified. The present study also shows that explicit treatment of the complex water-mediated hydrogen bonding network at the protein-protein interface is of crucial importance for obtaining reliable binding free energies. The successful reproduction of relative binding energies shows that this type of methodology can be very useful as an aid in rational design and redesign of biologically active macromolecules.     

Study of the insulin dimerization: binding free energy calculations and per-residue free energy decomposition

A calculation of the binding free energy for the dimerization of insulin has been performed using the molecular mechanics-generalized Born surface area approach. The calculated absolute binding free energy is -11.9 kcal/mol, in approximate agreement with the experimental value of -7.2 kcal/mol. The results show that the dimerization is mainly due to nonpolar interactions. The role of the hydrogen bonds between the 2 monomers appears to give the direction of the interactions. A per-atom decomposition of the binding free energy has been performed to identify the residues contributing most to the self association free energy. Residues B24-B26 are found to make the largest favorable contributions to the dimerization. Other residues situated at the interface between the 2 monomers were found to make favorable but smaller contributions to the dimerization: Tyr B16, Val B12, and Pro B28, and to an even lesser extent, Gly B23. The energy decomposition on a per-residue basis is in agreement with experimental alanine scanning data. The results obtained from a single trajectory (i.e., the dimer trajectory is also used for the monomer analysis) and 2 trajectories (i.e., separate trajectories are used for the monomer and dimer) are similar.     

Identification of potential inhibitors of Fasciola gigantica thioredoxin1: computational screening,molecular dynamics simulation,and binding free energy studies

Fasciola gigantica is the causative organism of fascioliasis and is responsible for major economic losses in livestock production globally. F. gigantica thioredoxin1 (FgTrx1) is an important redox-active enzyme involved in maintaining the redox homeostasis in the cell. To identify a potential anti-fasciolid compound, we conducted a structure-based virtual screening of natural compounds from the ZINC database (n = 1,67,740) against the FgTrx1 structure. The ligands were docked against FgTrx1 and 309 ligands were found to have better docking score. These compounds were evaluated for Lipinski and ADMET prediction, and 30 compounds were found to fit well for re-docking studies. After refinement by molecular docking and drug-likeness analysis, three potential inhibitors (ZINC15970091, ZINC9312362, and ZINC9312661) were identified. These three ligands were further subjected to molecular dynamics simulation (MDS) to compare the dynamics and stability of the protein structure after binding of the ligands. The binding free energy analyses were calculated to determine the intermolecular interactions. The results suggested that the two compounds had a binding free energy of –82.237, and –109.52 kJ.mol^?1 for compounds with IDs ZINC9312362 and ZINC9312661, respectively. These predicted compounds displayed considerable pharmacological and structural properties to be drug candidates. We concluded that these two compounds could be potential drug candidates to fight against F. gigantica parasites.     

Binding modes of Bcl-2 homology 3 (BH3) peptides with anti-apoptotic protein A1 and redesign of peptide inhibitors: a computational study

The interaction between protein and peptide ligand is a challenging problem in molecular biology and drug design. The binding of the Bcl-2 homology 3 (BH3) peptide to the anti-apoptotic protein A1 was revealed as a critical step in the regulation of apoptosis. These BH3 peptides hold high structural similarity, but are diverse in their regulation abilities. Based on molecular simulations and MM-P(G)BSA methods, this work presented a detailed analysis on binding mechanism of the BH3 peptides derived from PUMA and BMF. Residue-level energy decomposition showed that the core regions of BH3 peptides maintain in stable helical conformations and the four conserved hydrophobic residues together with an invariant aspartic acid contribute the major driving force for binding, whereas their two terminal segments exhibit obvious flexibility and various binding modes. Such kind of behavior was suggested as the reason for binding diversity and selectivity of BH3 peptides. As a further step, several BH3-mimetic peptides have been redesigned by computational mutation. Those new peptides showed not only stronger affinities when binding to protein A1, as well transferable binding patterns at some specific positions. A long-range coupling effect was disclosed for BH3 peptides, side-chain orientation and binding contribution of terminal residues were even affected by mutations at large sequence interval. Overall, this work reports that the binding modes of BH3 peptides are primarily dependent on its two terminal segments. The computational methods applied herein are also demonstrated to be of great assistance in the rational design of peptide inhibitors.     

Free energy calculations on binding and catalysis by alpha-lytic protease: the role of substrate size in the P1 pocket

We present free energy calculations using molecular dynamics on different substrates of alpha-lytic protease in the gas phase, in solution, while forming a noncovalent Michaelis complex with the enzyme, and in a tetrahedral structure representing a transition state/intermediate for acylation by the enzyme. Various P1 substrates were studied, with P1 = Gly, Ala, Val, and Leu. In qualitative agreement with experiment, the enzyme was calculated to bind and catalyze most effectively substrates with P1 = Ala over those with P1 = Gly, Val or Leu. Also, the calculated relative solvation free energies of Gly----Ala and Ala----Val were in qualitative agreement with experimental values in corresponding model systems. However, the level of quantitative agreement with experiment achieved in our earlier study of relative binding and catalysis of native subtilisin and an Asn-155----Ala mutant was not achieved. We surmise that this is due to the greater difficulty in quantitatively simulating effects that are predominantly van der Waals and hydrophobic compared to those that are hydrogen bonding/electrostatic.     

Flavonoids as potent allosteric inhibitors of protein tyrosine phosphatase 1B: molecular dynamics simulation and free energy calculation

Protein tyrosine phosphatase 1B (PTP1B) is a member of the PTP superfamily which is considered to be a negative regulator of insulin receptor (IR) signaling pathway. PTP1B is a promising drug target for the treatment of type 2 diabetes, obesity, and cancer. The existence of allosteric site in PTP1B has turned the researcher’s attention to an alternate strategy for inhibition of this enzyme. Herein, the molecular interactions between the allosteric site of PTP1B with three non-competitive flavonoids, (MOR), (MOK), and (DPO) have been investigated. Three ligands were docked into allosteric site of the enzyme. The resulting protein–ligand complexes were used for molecular dynamics studies. Principal component and free-energy landscape (FEL) as well as cluster analyses were used to investigate the conformational and dynamical properties of the protein and identify representative enzyme substrates bounded to the inhibitors. Per residue energy decomposition analysis attributed dissimilar affinities of three inhibitors to the several hydrogen bonds and non-bonded interactions. In conclusion, our results exhibited an inhibitory pattern of the ligands against PTP1B.     

The binding mechanism of a novel nicotinamide isostere inhibiting with TNKSs: a molecular dynamic simulation and binding free energy calculation

Tankyrases (TNKSs), a member of human poly (ADP-ribose) polymerase (PARP) protein superfamily, plays a key role in regulation of cell proliferation. Among the representative proteins of the PARPs family, it is found that the inhibitors have high selectivity for Tankyrase1 (TNKS1). The specific binding modes are investigated between the TNKS1 protein and nicotinamide isostere (ISX) which functions as an inhibitor of TNKS1. The stabilities of ISX-TNKS1 and AVA939-TNKS1 complexes are estimated by molecular dynamics (MD) simulations and free energy calculations; a good agreement with experimental results is reached. On the basis of the calculated results of MD simulations, we found that the inhibitors influence the conformational flexibility of TNKS1 and the XAV939 binding drive the peptide Ile1228-Gly1229-Gly1230 to form a helical structure while the ISX binding drive the peptide to form a turn structure. Moreover, the formed important hydrogen bonds of Tyr1203 residue with XVA939 and WAT1551 with ISX enhance stabilities of the complexes, and the electrostatic interactions in XAV939-TNKS1 and van der Waals interactions in ISX-TNKS1 system are main driving forces for affinity. According to the results of the decomposition of binding free energy, it is obvious that the residues Try1224 and Lys1220 make the most favorable contributions to the binding in, respectively, ISX and XAV939 complexes. Taken together, the obtained results are useful for studying the binding mechanisms of TNKSs and inhibitors and for designing potent inhibitors.     

Investigation on the binding mechanism of loratinib with the c-ros oncogene 1 (ROS1) receptor tyrosine kinase via molecular dynamics simulation and binding free energy calculations

The c-ros oncogene 1 (ROS1) has proven to be an important cancer target for the treatment of various human cancers. The anaplastic lymphoma kinase inhibitor crizotinib has been granted approval for the treatment of patients with ROS1 positive metastatic non-small-cell lung cancer by the Food and Drug Administration on 2016. However, serious resistance due to the secondary mutation of glycine 2032 to arginine (G2032R) was developed in clinical studies. Loratinib (PF-06463922), a macrocyclic analog of crizotinib, showed significantly improved inhibitory activity against wild–type (WT) ROS1 and ROS1^G2032R mutant. To provide insights into the inhibition mechanism, molecular dynamics simulations and free energy calculations were carried out for the complexes of loratinib with WT and G2032R mutated ROS1. The apo-ROS1^WT and apo-ROS1^G2032R systems showed similar RMSF distributions, while ROS1^G2032R-loratinib showed significantly higher than that of WT ROS1-loratinib, which revealed that the binding of loratinib to ROS1^G2032R significantly interfered the ?uctuation of protein. Calculations of binding free energies indicate that G2032R mutation significantly reduces the binding affinity of loratinib for ROS1, which arose mostly from the increase of conformation entropy and the decrease of solvation energy. Furthermore, detailed per-residue binding free energies highlighted the increased and decreased contributions of some residues in the G2032R mutated systems. The present study revealed the detailed inhibitory mechanism of loratinib as potent WT and G2032R mutated ROS1 inhibitor, which was expected to provide a basis for rational drug design.     

中华按蚊气味结合蛋白AsinOBP2的基因克隆、表达谱及与人体气味物质的结合特性分析

[目的]克隆中华按蚊Anopheles sinensis气味结合蛋白2(odorant binding protein 2,OBP2)基因AsinOBP2,分析该基因的表达及其重组蛋白与人体气味物质的结合能力.[方法]采用RT-PCR和RACE技术克隆AsinOBP2的全长cDNA序列,通过qPCR分析AsinOBP2...     

Calculating proton uptake/release and binding free energy taking into account ionization and conformation changes induced by protein-inhibitor association: application to plasmepsin, cathepsin D and endothiapepsin-pepstatin complexes

The protein-inhibitor binding energies of enzymes are often pH dependent, and binding induces either proton uptake or proton release. The proton uptake/release and the binding energy for three complexes with available experimental data were numerically studied: pepstatin-cathepsin D, pepstatin-plasmepsin II and pepstatin-endothiapepsin. Very good agreement with the experimental data was achieved when conformational changes were taken into account. The role of the desolvation energy and the conformational changes was revealed by modeling the complex, the separated molecules in the complex conformation and the free molecules. It was shown that the conformational changes induced by the complex formation are as important for the proton transfer as the loss of solvation energy caused by the burial of interface residues. The residues responsible for the proton transfer were identified and their contribution to the proton uptake/release calculated. These residues were found to be scattered along the whole protein rather than being localized only at the active site. In the case of cathepsin D, these residues were found to be highly conserved among the cathepsin D sequences of other species. It was shown that conformation and ionization changes induced by the complex formation are critical for the correct calculation of the binding energy. Taking into account the electrostatics and the van der Waals (vdW) energies within the Boltzmann distribution of energies and allowing ionization and conformation changes to occur makes the calculated binding energy more realistic and closer to the experimental value. The interplay between electrostatic and vdW forces makes the pH dependence of the binding energy smoother, because the vdW force acts in reaction to the changes of the electrostatic energy. It was found that a small fraction of the ionizable groups remain uncharged in both the free and complexed molecules. The sequence and structural position of these groups aligns well within the three proteases, suggesting that these may have specific role.     

Computational design of small molecular modulators of protein–protein interactions with a novel thermodynamic cycle: Allosteric inhibitors of HIV‐1 integrase

Targeting protein–protein interactions for therapeutic development involves designing small molecules to either disrupt or enhance a known PPI. For this purpose, it is necessary to compute reliably the effect of chemical modifications of small molecules on the protein–protein association free energy. Here we present results obtained using a novel thermodynamic free energy cycle, for the rational design of allosteric inhibitors of HIV‐1 integrase (ALLINI) that act specifically in the early stage of the infection cycle. The new compounds can serve as molecular probes to dissect the multifunctional mechanisms of ALLINIs to inform the discovery of new allosteric inhibitors. The free energy protocol developed here can be more broadly applied to study quantitatively the effects of small molecules on modulating the strengths of protein–protein interactions.     

斜纹夜蛾普通气味结合蛋白GOBP1基因的表达定位分析

气味调控斜纹夜蛾Spodoptera litura（鳞翅目,夜蛾科）的交配和产卵行为,而嗅觉气味结合蛋白（OBP）是昆虫与外界环境进行化学信息交流中的一种重要蛋白。本研究基于已报道的斜纹夜蛾普通气味结合蛋白GOBP1基因cDNA序列（GenBank登录号为EF159978）设计引物,通过RT-PCR法分析了GOBP1 mRNA的组织特异性表达情况,并对GOBP1基因条带进行了克隆、测序和Blast比对。此外,再通过RNA原位杂交的方法进一步分析了GOBP1 mRNA在触角中的表达定位。结果表明:GOBP1只在斜纹夜蛾的触角组织中表达,而且GOBP1转录本较多分布于靠近嗅觉感受器即触角边缘部位。这些结果进一步说明GOBP1是斜纹夜蛾嗅觉过程中的重要蛋白,同时为深入研究GOBP1与其他蛋白的相互空间定位关系、GOBP1的功能等奠定了基础。     

Thermodynamic characterization of two homologous protein complexes: Associations of the semaphorin receptor plexin‐B1 RhoGTPase binding domain with Rnd1 and active Rac1

Plexin receptors function in response to semaphorin guidance cues in a variety of developmental processes involving cell motility. Interactions with Rho, as well as Ras family small GTPases are critical events in the cell signaling mechanism. We have recently determined the structure of a cytoplasmic domain (RBD) of plexin‐B1 and mapped its binding interface with several Rho‐GTPases, Rac1, Rnd1, and RhoD. All three GTPases associate with a similar region of this plexin domain, but show different functional behavior in cells. To understand whether thermodynamic properties of the GTPase–RBD interaction contribute to such different behavior, we have examined the interaction at different temperatures, buffer, and pH conditions. Although the binding affinity of both Rnd1 and Rac1 with the plexin‐B1 RBD is similar, the detailed thermodynamic properties of the interactions are considerably different. These data suggest that on Rac1 binding to the plexin‐B1 RBD, the proteins become more rigid in the complex. By contrast, Rnd1 binding is consistent with unchanged or slightly increased flexibility in one or both proteins. Both GTPases show an appreciable reduction in affinity for the dimeric plexin‐B1 RBD indicating that GTPase binding is not cooperative with dimer formation, but that a partial steric hindrance destabilizes the dimer. However, a reduced affinity binding mode to a disulphide stabilized model for the dimeric RBD is also possible. Consistent with cellular studies, the interaction thermodynamics imply that further levels of regulation involving additional binding partners and/or regions outside of the RhoGTPase binding domain are required for receptor activation.     

Hierarchy of simulation models in predicting molecular recognition mechanisms from the binding energy landscapes: structural analysis of the peptide complexes with SH2 domains.

Computer simulations using the simplified energy function and simulated tempering dynamics have accurately determined the native structure of the pYVPML, SVLpYTAVQPNE, and SPGEpYVNIEF peptides in the complexes with SH2 domains. Structural and equilibrium aspects of the peptide binding with SH2 domains have been studied by generating temperature-dependent binding free energy landscapes. Once some native peptide-SH2 domain contacts are constrained, the underlying binding free energy profile has the funnel-like shape that leads to a rapid and consistent acquisition of the native structure. The dominant native topology of the peptide-SH2 domain complexes represents an extended peptide conformation with strong specific interactions in the phosphotyrosine pocket and hydrophobic interactions of the peptide residues C-terminal to the pTyr group. The topological features of the peptide-protein interface are primarily determined by the thermodynamically stable phosphotyrosyl group. A diversity of structurally different binding orientations has been observed for the amino-terminal residues to the phosphotyrosine. The dominant native topology for the peptide residues carboxy-terminal to the phosphotyrosine is tolerant to flexibility in this region of the peptide-SH2 domain interface observed in equilibrium simulations. The energy landscape analysis has revealed a broad, entropically favorable topology of the native binding mode for the bound peptides, which is robust to structural perturbations. This could provide an additional positive mechanism underlying tolerance of the SH2 domains to hydrophobic conservative substitutions in the peptide specificity region.     

Molecular dynamics simulation,binding free energy calculation and unbinding pathway analysis on selectivity difference between FKBP51 and FKBP52: Insight into the molecular mechanism of isoform selectivity

As co‐chaperones of the 90‐kDa heat shock protein(HSP90), FK506 binding protein 51 (FKBP51) and FK506 binding protein 52 (FKBP52) modulate the maturation of steroid hormone receptor through their specific FK1 domains (FKBP12‐like domain 1). The inhibitors targeting FK1 domains are potential therapies for endocrine‐related physiological disorders. However, the structural conservation of the FK1 domains between FKBP51 and FKBP52 make it difficult to obtain satisfactory selectivity in FK506‐based drug design. Fortunately, a series of iFit ligands synthesized by Hausch et al exhibited excellent selectivity for FKBP51, providing new opportunity for design selective inhibitors. We performed molecular dynamics simulation, binding free energy calculation and unbinding pathway analysis to reveal selective mechanism for the inhibitor iFit4 binding with FKBP51 and FKBP52. The conformational stability evaluation of the “Phe67‐in” and “Phe67‐out” states implies that FKBP51 and FKBP52 have different preferences for “Phe67‐in” and “Phe67‐out” states, which we suggest as the determinant factor for the selectivity for FKBP51. The binding free energy calculations demonstrate that nonpolar interaction is favorable for the inhibitors binding, while the polar interaction and entropy contribution are adverse for the inhibitors binding. According to the results from binding free energy decomposition, the electrostatic difference of residue 85 causes the most significant thermodynamics effects on the binding of iFit4 to FKBP51 and FKBP52. Furthermore, the importance of substructure units on iFit4 were further evaluated by unbinding pathway analysis and residue‐residue contact analysis between iFit4 and the proteins. The results will provide new clues for the design of selective inhibitors for FKBP51.     

Computational analysis of molecular basis of 1:1 interactions of NRG-1beta wild-type and variants with ErbB3 and ErbB4

The neuregulin/ErbB system is a growth factor/receptor cascade that has been proven to be essential in the development of the heart and the sympathetic nervous system. However, the basis of the specificity of ligand-receptor recognition remains to be elucidated. In this study, the structures of NRG-1beta/ErbB3 and NRG-1beta/ErbB4 complexes were modeled based on the available structures of the homologous proteins. The binding free energies of NRG-1beta to ErbB3 and ErbB4 were calculated using the molecular mechanics Poisson-Boltzmann surface area (MM-PBSA) computational method. In addition, computational alanine-scanning mutagenesis was performed in the binding site of NRG-1beta and the difference in the binding free energies between NRG-1beta mutants and the receptors was calculated. The results specify the contribution of each residue at the interaction interfaces to the binding affinity of NRG-1beta with ErbB3 and ErbB4, identifying several important interaction residue pairs that are in agreement with previously acquired experimental data. This indicates that the presented structural models of NRG-1beta/ErbB3 and NRG-1beta/ErbB4 complexes are reliable and could be used to guide future studies, such as performing desirable mutations on NRG-1beta to increase the binding affinity and selectivity to the receptor and discovering new therapeutic agents for the treatment of heart failure.