In silico identification of the potential drug resistance sites over 2009 influenza A (H1N1) virus neuraminidase

Pharmacophore-Assisted Covalent Docking Identifies a Potential Covalent Inhibitor for Drug-Resistant Genotype 3 Variants of Hepatitis C Viral NS3/4A Serine Protease

The emergence of drug-resistance-inducing mutations in Hepatitis C virus (HCV) coupled with genotypic heterogeneity has made targeting NS3/4A serine protease difficult. In this work, we investigated the mutagenic variations in the binding pocket of Genotype 3 (G3) HCV NS3/4A and evaluated ligands for efficacious inhibition. We report mutations at 14 positions within the ligand-binding residues of HCV NS3/4A, including H57R and S139P within the catalytic triad. We then modelled each mutational variant for pharmacophore-based virtual screening (PBVS) followed by covalent docking towards identifying a potential covalent inhibitor, i.e., cpd-217. The binding stability of cpd-217 was then supported by molecular dynamic simulation followed by MM/GBSA binding free energy calculation. The free energy decomposition analysis indicated that the resistant mutants alter the HCV NS3/4A-ligand interaction, resulting in unbalanced energy distribution within the binding site, leading to drug resistance. Cpd-217 was identified as interacting with all NS3/4A G3 variants with significant covalent docking scores. In conclusion, cpd-217 emerges as a potential inhibitor of HCV NS3/4A G3 variants that warrants further in vitro and in vivo studies. This study provides a theoretical foundation for drug design and development targeting HCV G3 NS3/4A.

Discovery of Novel HPK1 Inhibitors Through Structure-Based Virtual Screening

Hematopoietic progenitor kinase (HPK1) is a negative regulator of T-cell receptor and B-cell signaling, which has been recognized as a novel antitumor target for immunotherapy. In this work, Glide docking-based virtual screening and kinase inhibition assay were performed to identify novel HPK1 inhibitors. The kinase inhibition assay results demonstrated five compounds with IC50 values below 20 μM, and the most potent one (compound M074-2865) had an IC50 value of 2.93 ± 0.09 μM. Molecular dynamics (MD) simulations were performed to delve into the interaction of sunitinib and the identified compound M074-2865 with the kinase domain of HPK1. The five compounds identified in this work could be considered promising hit compounds for further development of HPK1 inhibitors for immunotherapy.

Revealing the binding and drug resistance mechanism of amprenavir, indinavir, ritonavir, and nelfinavir complexed with HIV-1 protease due to double mutations G48T/L89M by molecular dynamics simulations and free energy analyses

MD simulations, MM-PBSA, and SIE analyses were used to investigate the drug resistance mechanisms of two mutations G48T and L89M in HIV-1 protease toward four inhibitors.

Insight into the selective mechanism of phosphoinositide 3-kinase γ with benzothiazole and thiazolopiperidine γ-specific inhibitors by in silico approaches

The phosphoinositide 3-kinase γ (PI3Kγ) has been verified to be a potential drug target for the treatments of various human physical disorders. Although received lots of attention, the development of PI3Kγ-selective inhibitors is still a challenging subject because of its unique protein structural features. Aiming to uncover the interaction mechanism between the selective inhibitors and PI3Kγ, a series of benzothiazole and thiazolopiperidine PI3Kγ isoform-selective inhibitors were studied with an integrated in silico strategy by combining molecular docking, molecular dynamic simulations, binding free energy calculations, and decomposition analysis. Firstly, three molecular docking models, including rigid receptor docking, induced fit docking (IFD), and quantum mechanical-polarized ligand docking, were respectively, built, and the IFD preliminarily predicted the docking poses of all studied inhibitors and roughly analyzed the binding mechanism. Secondly, four binding complexes with representative inhibitors were selected to perform molecular dynamic simulations and free energy calculations. The predicted binding energies were consistent with the experimental bioactivities and different binding patterns between potent and weak inhibitors were uncovered. Finally, through the Molecular Mechanics/Generalized Born Surface Area binding free energy decomposition, residue-inhibitor interactions spectra were obtained and several key residues contributing to favorable binding were highlighted, which provides valuable information for rational PI3Kγ inhibitor design and modification.

Drug-resistance mechanisms of three mutations in anaplastic lymphoma kinase against two inhibitors based on MM/PBSA combined with interaction entropy

As a promising drug target in the treatment of lung cancer, anaplastic lymphoma kinase (ALK) and its mutations have been studied widely. This work explored the origin of the resistance mechanism of the ALK mutants again two inhibitors.

Understanding the structural and energetic basis of PD-1 and monoclonal antibodies bound to PD-L1: A molecular modeling perspective

BACKGROUND

The inhibitors blocking the interaction between programmed cell death protein 1(PD-1) and programmed death-ligand 1(PD-L1) can activate the immune response of T cell and eliminate cancer cells. The crystallographic studies have provided structural insights of the interactive interfaces between PD-L1 and its protein ligands. However, the hotspot residues on PD-L1 as well as structural and energetic basis for different protein ligands still need to be further investigated.

METHODS

Molecular modeling methods including molecular dynamics simulation, per-residue free energy decomposition, virtual alanine scanning mutagenesis and residue-residue contact analysis were used to qualitatively and quantitatively analyze the interactions between PD-L1 and different protein ligands.

RESULTS

The results of virtual alanine scanning mutagenesis suggest that Y56, Q66, M115, D122, Y123, R125 are the hotspot residues on PD-L1. The residue-residue contact analysis further shows that PD-1 interacts with PD-L1 mainly by F and G strands while monoclonal antibodies like avelumab and BMS-936559 mainly interact with PD-L1 by CDR2 and CDR3 loops of the heavy chain.

CONCLUSIONS

A structurally similar β-hairpin peptide with 13 or 14 residues was extracted from each protein ligand and these β-hairpin peptides were found tightly binding to the putative hotspot residues on PD-L1.

GENERAL SIGNIFICANCE

This study recognizes the hotspot residues on PD-L1 and uncovers the common structural and energetic basis of different protein ligands binding to PD-L1. These results will be valuable for the design of small molecule or peptide inhibitors targeting on PD-L1.

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.

Theoretical studies on the selective mechanisms of GSK3β and CDK2 by molecular dynamics simulations and free energy calculations

Glycogen synthase kinase 3 (GSK3) is a serine/threonine protein kinase which is widely involved in cell signaling and controls a broad number of cellular functions. GSK3 contains α and β isoforms, and GSK3β has received more attention and becomes an attractive drug target for the treatment of several diseases. The binding pocket of cyclin-dependent kinase 2 (CDK2) shares high sequence identity to that of GSK3β, and therefore, the design of highly selective inhibitors toward GSK3β remains a big challenge. In this study, a computational strategy, which combines molecular docking, molecular dynamics simulations, free energy calculations, and umbrella sampling simulations, was employed to explore the binding mechanisms of two selective inhibitors to GSK3β and CDK2. The simulation results highlighted the key residues critical for GSK3β selectivity. It was observed that although GSK3β and CDK2 share the conserved ATP-binding pockets, some different residues have significant contributions to protein selectivity. This study provides valuable information for understanding the GSK3β-selective binding mechanisms and the rational design of selective GSK3β inhibitors.

Role of R292K mutation in influenza H7N9 neuraminidase toward oseltamivir susceptibility: MD and MM/PB(GB)SA study

The H7N9 avian influenza virus is a novel re-assortment from at least four different strains of virus. Neuraminidase, which is a glycoprotein on the surface membrane, has been the target for drug treatment. However, some H7N9 strains that have been isolated from patient after drug treatment have a R292K mutation in neuraminidase. This substitution was found to facilitate drug resistance using protein- and virus- assays, in particular it gave a high resistance to the most commonly used drug, oseltamivir. The aim of this research is to understand the source of oseltamivir resistance using MD simulations and the MM/PB(GB)SA binding free energy approaches. Both methods can predict the reduced susceptibility of oseltamivir in good agreement to the IC ₅₀ binding energy, although MM/GBSA underestimates this prediction compared to the MM/PBSA calculation. Electrostatic interaction is the main contribution for oseltamivir binding in terms of both interaction and solvation. We found that the source of the drug resistance is a decrease in the binding interaction combined with the reduction of the dehydration penalty. The smaller K292 mutated residue has a larger binding pocket cavity compared to the wild-type resulting in the loss of drug carboxylate-K292 hydrogen bonding and an increased accessibility for water molecules around the K292 mutated residue. In addition, oseltamivir does not bind well to the R292K mutant complex as shown by the high degree of fluctuation in ligand RMSD during the simulation and the change in angular distribution of bulky side chain groups.

Binding mechanism of CDK5 with roscovitine derivatives based on molecular dynamics simulations and MM/PBSA methods

Roscovitine derivatives are potent inhibitors of cyclin-dependent kinase 5 (CDK5), but they exhibit different activities, which has not been understood clearly up to now. On the other hand, the task of drug design is difficult because of the fuzzy binding mechanism. In this context, the methods of molecular docking, molecular dynamics (MD) simulation, and binding free energy analysis are applied to investigate and reveal the detailed binding mechanism of four roscovitine derivatives with CDK5. The electrostatic and van der Waals interactions of the four inhibitors with CDK5 are analyzed and discussed. The calculated binding free energies in terms of MM-PBSA method are consistent with experimental ranking of inhibitor effectiveness for the four inhibitors. The hydrogen bonds of the inhibitors with Cys83 and Lys33 can stabilize the inhibitors in binding sites. The van der Waals interactions, especially the pivotal contacts with Ile10 and Leu133 have larger contributions to the binding free energy and play critical roles in distinguishing the variant bioactivity of four inhibitors. In terms of binding mechanism of the four inhibitors with CDK5 and energy contribution of fragments of each inhibitor, two new CDK5 inhibitors are designed and have stronger inhibitory potency.

Exploring the mechanism how Marburg virus VP35 recognizes and binds dsRNA by molecular dynamics simulations and free energy calculations

Filoviruses often cause terrible infectious disease which has not been successfully dealt with pharmacologically. All filoviruses encode a unique protein termed VP35 which can mask doubled-stranded RNA to deactivate interferon. The interface of VP35-dsRNA would be a feasible target for structure-based antiviral agent design. To explore the essence of VP35-dsRNA interaction, molecular dynamics simulation combined with MM-GBSA calculations were performed on Marburg virus VP35-dsRNA complex and several mutational complexes. The energetic analysis indicates that nonpolar interactions provide the main driving force for the binding process. Although the intermolecular electrostatic interactions play important roles in VP35-dsRNA interaction, the whole polar interactions are unfavorable for binding which result in a low binding affinity. Compared with wild type VP35, the studied mutants F228A, R271A, and K298A have obviously reduced binding free energies with dsRNA reflecting in the reduction of polar or nonpolar interactions. The results also indicate that the loss of binding affinity for one dsRNA strand would abolish the total binding affinity. Three important residues Arg271, Arg294, and Lys298 which makes the largest contribution for binding in VP35 lose their binding affinity significantly in mutants. The uncovering of VP35-dsRNA recognition mechanism will provide some insights for development of antiviral drug.

Identification of potential quinoxalinone-based aldose reductase inhibitors by 3D-QSAR, molecular docking and molecular dynamics

The 3D-QSAR model of aldose reductase (ARIs) inhibitors is built to gain insights into the key structural factors affecting the inhibitory activity. Based on the model, six new potential ARIs were designed.

Mixed Inhibition of cPEPCK by Genistein, Using an Extended Binding Site Located Adjacent to Its Catalytic Cleft

Cytosolic phosphoenolpyruvate carboxykinase (cPEPCK) is a critical enzyme involved in gluconeogenesis, glyceroneogenesis and cataplerosis. cPEPCK converts oxaloacetic acid (OAA) into phosphoenol pyruvate (PEP) in the presence of GTP. cPEPCK is known to be associated with type 2 diabetes. Genistein is an isoflavone compound that shows anti-diabetic and anti-obesitic properties. Experimental studies have shown a decrease in the blood glucose level in the presence of genistein by lowering the functional activity of cPEPCK, an enzyme of gluconeogenesis. Using computational techniques such as molecular modeling, molecular docking, molecular dynamics simulation and binding free energy calculations, we identified cPEPCK as a direct target of genistein. We studied the molecular interactions of genistein with three possible conformations of cPEPCK-unbound cPEPCK (u_cPEPCK), GTP bound cPEPCK (GTP_cPEPCK) and GDP bound cPEPCK (GDP_cPEPCK). Binding of genistein was also compared with an already known cPEPCK inhibitor. We analyzed the interactions of genistein with cPEPCK enzyme and compared them with its natural substrate (OAA), product (PEP) and known inhibitor (3-MPA). Our results demonstrate that genistein uses the mechanism of mixed inhibition to block the functional activity of cPEPCK and thus can serve as a potential anti-diabetic and anti-obesity drug candidate. We also identified an extended binding site in the catalytic cleft of cPEPCK which is used by 3-MPA to inhibit cPEPCK non-competitively. We demonstrate that extended binding site of cPEPCK can further be exploited for designing new drugs against cPEPCK.

Exploring resistance mechanisms of HCV NS3/4A protease mutations to MK5172: insight from molecular dynamics simulations and free energy calculations

Mutations at a number of key positions (Ala156, Asp168 and Arg155) of the HCV NS3/4A protease can induce medium to high resistance to MK5172.

Molecular basis of R294K mutation effects of H7N9 neuraminidases with drugs and cyclic peptides: an in silico and experimental study

An overview of Shanghai N9/cyclic peptide I complex structure.

Exploring the prominent performance of CX-4945 derivatives as protein kinase CK2 inhibitors by a combined computational study

Protein kinase CK2, also known as casein kinase II, is related to various cellular events and is a potential target for numerous cancers. In this study, we attempted to gain more insight into the inhibition process of CK2 by a series of CX-4945 derivatives through an integrated computational study that combines molecular docking, molecular dynamics (MD) simulations, and binding free energy calculations. Based on the binding poses predicted by molecular docking, the MD simulations were performed to explore the dynamic binding processes for ten selected inhibitors. Then, both Molecular Mechanics/Poisson Boltzmann Surface Area (MM/PBSA) and Molecular Mechanics/Generalized Born Surface Area (MM/GBSA) techniques were employed to predict the binding affinities of the studied systems. The predicted binding energies of the selected inhibitors correlate well with their experimental activities (r(2) = 0.78). The van der Waals term is the most favorable component for the total energies. The free energy decomposition on a per residue basis reveals that the residue K68 is essential for the electrostatic interactions between CK2 and the studied inhibitors and numerous residues, including L45, V53, V66, F113, M163 and I174, play critical roles in forming van der Waals interactions with the inhibitors. Finally, a number of new derivatives were designed and the binding affinity and the predicted binding free energies of each designed molecule were obtained on the basis of molecular docking and MM/PBSA. It is expected that our research will benefit the future rational design of novel and potent inhibitors of CK2.

Understanding the recognition mechanisms of Zα domain of human editing enzyme ADAR1 (hZα(ADAR1)) and various Z-DNAs from molecular dynamics simulation

The Z-DNA-binding domain of human double-stranded RNA adenosine deaminase I (hZαADAR1) can specifically recognize the left-handed Z-DNA which preferentially occurs at alternating purine-pyrimidine repeats, especially the CG-repeats. The interactions of hZαADAR1 and Z-DNAs in different sequence contexts can affect many important biological functions including gene regulation and chromatin remodeling. Therefore it is of great necessity to fully understand their recognition mechanisms. However, most existing studies are aimed at the standard CG-repeat Z-DNA rather than the non-CG-repeats, and whether the molecular basis of hZαADAR1 binding to various Z-DNAs are identical or not is still unclear on the atomic level. Here, based on the recently determined crystal structures of three representative non-CG-repeat Z-DNAs (d(CACGTG)2, d(CGTACG)2 and d(CGGCCG)2) in complex with hZαADAR1, 40 ns molecular dynamics simulation together with binding free energy calculation were performed for each system. For comparison, the standard CG-repeat Z-DNA (d(CGCGCG)2) complexed with hZαADAR1 was also simulated. The consistent results demonstrate that nonpolar interaction is the driving force during the protein-DNA binding process, and that polar interaction mainly from helix α3 also provides important contributions. Five common hot-spot residues were identified, namely Lys169, Lys170, Asn173, Arg174 and Tyr177. Hydrogen bond analysis coupled with surface charge distribution further reveal the interfacial information between hZαADAR1 and Z-DNA in detail. All of the analysis illustrate that four complexes share the common key features and the similar binding modes irrespective of Z-DNA sequences, suggesting that Z-DNA recognition by hZαADAR1 is conformation-specific rather than sequence-specific. Additionally, by analyzing the conformational changes of hZαADAR1, we found that the binding of Z-DNA could effectively stabilize hZαADAR1 protein. Our study can provide some valuable information for better understanding the binding mechanism between hZαADAR1 or even other Z-DNA-binding protein and Z-DNA.

Theoretical study on the interaction of pyrrolopyrimidine derivatives as LIMK2 inhibitors: insight into structure-based inhibitor design

LIM kinases (LIMKs), downstream of Rho-associated protein kinases (ROCKs) and p21-activated protein kinases (PAKs), are shown to be promising targets for the treatment of cancers. In this study, the inhibition mechanism of 41 pyrrolopyrimidine derivatives as LIMK2 inhibitors was explored through a series of theoretical approaches. First, a model of LIMK2 was generated through molecular homology modeling, and the studied inhibitors were docked into the binding active site of LIMK2 by the docking protocol, taking into consideration the flexibility of the protein. The binding poses predicted by molecular docking for 17 selected inhibitors with different bioactivities complexed with LIMK2 underwent molecular dynamics (MD) simulations, and the binding free energies for the complexes were predicted by using the molecular mechanics/generalized born surface area (MM/GBSA) method. The predicted binding free energies correlated well with the experimental bioactivities (r(2) = 0.63 or 0.62). Next, the free energy decomposition analysis was utilized to highlight the following key structural features related to biological activity: (1) the important H-bond between Ile408 and pyrrolopyrimidine, (2) the H-bonds between the inhibitors and Asp469 and Gly471 which maintain the stability of the DFG-out conformation, and (3) the hydrophobic interactions between the inhibitors and several key residues (Leu337, Phe342, Ala345, Val358, Lys360, Leu389, Ile408, Leu458 and Leu472). Finally, a variety of LIMK2 inhibitors with a pyrrolopyrimidine scaffold were designed, some of which showed improved potency according to the predictions. Our studies suggest that the use of molecular docking with MD simulations and free energy calculations could be a powerful tool for understanding the binding mechanism of LIMK2 inhibitors and for the design of more potent LIMK2 inhibitors.

Insights from QM/MM Modeling the 3D Structure of the 2009 H1N1 Influenza A Virus Neuraminidase and Its Binding Interactions with Antiviral Drugs

Control of drug release through the inhibition of neuraminidase-1 has been identified as a potential target for the treatment of H1N1 influenza; however, the drug binding mode of neuraminidase is not yet completely understood. In this work, we propose a model for a neuraminidase-1 complex based on four known X-ray structures of drug/neuraminidase-1 complexes. Specifically, H1N1 neuraminidase-1 complexed with 4 drugs (zanamivir, laninamivir, laninamivir octanoate and oseltamivir) was first investigated using a combined quantum mechanical and molecular mechanical (QM/MM) approach. Based on these structures, a model for the H1N1 neuraminidase-1 complex was proposed and simulated using the same computational protocol. Implications to drug/H1N1 neuraminidase-1 binding modes are discussed. From our simulations, the predicted binding free energies of the four drugs are in good agreement with the experimental results, with the correlation coefficient being 0.84.

Molecular modeling and residue interaction network studies on the mechanism of binding and resistance of the HCV NS5B polymerase mutants to VX-222 and ANA598

Hepatitis C virus (HCV) NS5B protein is an RNA-dependent RNA polymerase (RdRp) with essential functions in viral genome replication and represents a promising therapeutic target to develop direct-acting antivirals (DAAs). Multiple nonnucleoside inhibitors (NNIs) binding sites have been identified within the polymerase. VX-222 and ANA598 are two NNIs targeting thumb II site and palm I site of HCV NS5B polymerase, respectively. These two molecules have been shown to be very effective in phase II clinical trials. However, the emergence of resistant HCV replicon variants (L419M, M423T, I482L mutants to VX-222 and M414T, M414L, G554D mutants to ANA598) has significantly decreased their efficacy. To elucidate the molecular mechanism about how these mutations influenced the drug binding mode and decreased drug efficacy, we studied the binding modes of VX-222 and ANA598 to wild-type and mutant polymerase by molecular modeling approach. Molecular dynamics (MD) simulations results combined with binding free energy calculations indicated that the mutations significantly altered the binding free energy and the interaction for the drugs to polymerase. The further per-residue binding free energy decomposition analysis revealed that the mutations decreased the interactions with several key residues, such as L419, M423, L474, S476, I482, L497, for VX-222 and L384, N411, M414, Y415, Q446, S556, G557 for ANA598. These were the major origins for the resistance to these two drugs. In addition, by analyzing the residue interaction network (RIN) of the complexes between the drugs with wild-type and the mutant polymerase, we found that the mutation residues in the networks involved in the drug resistance possessed a relatively lower size of topology centralities. The shift of betweenness and closeness values of binding site residues in the mutant polymerase is relevant to the mechanism of drug resistance of VX-222 and ANA598. These results can provide an atomic-level understanding about the mechanisms of drug resistance conferred by the studied mutations and will be helpful to design more potent inhibitors which could effectively overcome drug resistance of antivirus agents.

Assessing the performance of MM/PBSA and MM/GBSA methods. 5. Improved docking performance using high solute dielectric constant MM/GBSA and MM/PBSA rescoring

We have thoroughly investigated the performance of MM/GBSA and MM/PBSA methodologies on virtual screening based on various protocols for kinase targets.

Theoretical studies on beta and delta isoform-specific binding mechanisms of phosphoinositide 3-kinase inhibitors

Phosphoinositide 3-kinase (PI3K) is known to be closely related to tumorigenesis and cell proliferation, and controls a variety of cellular processes, including proliferation, growth, apoptosis, migration, metabolism, etc. The PI3K family comprises eight catalytic isoforms, which are subdivided into three classes. Recently, the discovery of inhibitors that block a single isoform of PI3K has continued to attract special attention because they may have higher selectivity for certain tumors and less toxicity for healthy cells. The PI3Kβ and PI3Kδ share fewer studies than α/γ, and therefore, in this work, the combination of molecular dynamics simulations and free energy calculations was employed to explore the binding of three isoform-specific PI3K inhibitors (COM8, IC87114, and GDC-0941) to PI3Kβ or PI3Kδ. The isoform specificities of the studied inhibitors derived from the predicted binding free energies are in good agreement with the experimental data. In addition, the key residues critical for PI3Kβ or PI3Kδ selectivity were highlighted by decomposing the binding free energies into the contributions from individual residues. It was observed that although PI3Kβ and PI3Kδ share the conserved ATP-binding pockets, individual residues do behave differently, particularly the residues critical for PI3Kβ or PI3Kδ selectivity. It can be concluded that the inhibitor specificity between PI3Kβ and PI3Kδ is determined by the additive contributions from multiple residues, not just a single one. This study provides valuable information for understanding the isoform-specific binding mechanisms of PI3K inhibitors, and should be useful for the rational design of novel and selective PI3K inhibitors.

Insight into crizotinib resistance mechanisms caused by three mutations in ALK tyrosine kinase using free energy calculation approaches

As a safe and efficacious drug, crizotinib was approved by the U.S. Food and Drug Administration (FDA) in 2011 for the treatment of advanced fusion-type nonsmall-cell lung cancer. Although high response ratio was detected from the patients treated with crizotinib, the cancer has eventually conferred resistance to crizotinib. Several drug resistance mutations have been found in the anaplastic lymphoma kinase (ALK) tyrosine kinase domain as the target for crizotinib, but the drug resistance mechanisms remain unclear. Therefore, in this study, the adaptive biasing force (ABF) method and two-end-state free energy calculation approaches were employed to elucidate the resistance mechanisms of crizotinib induced by the mutations L1152R, G1202R, and S1206Y. The ABF simulation results suggest that the reaction coordinates for the unbinding processes of crizotinib from the binding pockets of the mutated ALKs is different from that of the wild type ALK. The potentials of mean force for the crizotinib unbinding and the binding free energies predicted by the two-end-state free energy calculations are consistent with the experimental data. Our results indicate that the three mutations weaken the binding affinity of crizotinib obviously and lead to drug resistance. The free energy decomposition analysis illustrates the importance of the loss of two important H-bonds in the L1152R and S1206Y mutants on drug resistance. The entropy analysis shows that the entropy term plays a critical role in the substantial change of the conformational entropies of G1202R and L1152R. Our results reveal the mechanisms of drug resistance and provide vital clues for the development of new inhibitors to combat drug resistance.

Molecular distribution of amino acid substitutions on neuraminidase from the 2009 (H1N1) human influenza pandemic virus

The pandemic influenza AH1N1 (2009) caused an outbreak of human infection that spread to the world. Neuraminidase (NA) is an antigenic surface glycoprotein, which is essential to the influenza infection process, and is the target of anti-flu drugs oseltamivir and zanamivir. Currently, NA inhibitors are the pillar pharmacological strategy against seasonal and global influenza. Although mutations observed after NA-inhibitor treatment are characterized by changes in conserved amino acids of the enzyme catalytic site, it is possible that specific amino acid substitutions (AASs) distant from the active site such as H274Y, could confer oseltamivir or zanamivir resistance. To better understand the molecular distribution pattern of NA AASs, we analyzed NA AASs from all available reported pandemic AH1N1 NA sequences, including those reported from America, Africa, Asia, Europe, Oceania, and specifically from Mexico. The molecular distributions of the AASs were obtained at the secondary structure domain level for both the active and catalytic sites, and compared between geographic regions. Our results showed that NA AASs from America, Asia, Europe, Oceania and Mexico followed similar molecular distribution patterns. The compiled data of this study showed that highly conserved amino acids from the NA active site and catalytic site are indeed being affected by mutations. The reported NA AASs follow a similar molecular distribution pattern worldwide. Although most AASs are distributed distantly from the active site, this study shows the emergence of mutations affecting the previously conserved active and catalytic site. A significant number of unique AASs were reported simultaneously on different continents.

Computational study on the drug resistance mechanism against HCV NS3/4A protease inhibitors vaniprevir and MK-5172 by the combination use of molecular dynamics simulation, residue interaction network, and substrate envelope analysis

Hepatitis C virus (HCV) NS3/4A protease is an important and attractive target for anti-HCV drug development and discovery. Vaniprevir (phase III clinical trials) and MK-5172 (phase II clinical trials) are two potent antiviral compounds that target NS3/4A protease. However, the emergence of resistance to these two inhibitors reduced the effectiveness of vaniprevir and MK-5172 against viral replication. Among the drug resistance mutations, three single-site mutations at residues Arg155, Ala156, and Asp168 in NS3/4A protease are especially important due to their resistance to nearly all inhibitors in clinical development. A detailed understanding of drug resistance mechanism to vaniprevir and MK-5172 is therefore very crucial for the design of novel potent agents targeting viral variants. In this work, molecular dynamics (MD) simulation, binding free energy calculation, free energy decomposition, residue interaction network (RIN), and substrate envelope analysis were used to study the detailed drug resistance mechanism of the three mutants R155K, A156T, and D168A to vaniprevir and MK-5172. MD simulation was used to investigate the binding mode for these two inhibitors to wild-type and resistant mutants of HCV NS3/4A protease. Binding free energy calculation and free energy decomposition analysis reveal that drug resistance mutations reduced the interactions between the active site residues and substituent in the P2 to P4 linker of vaniprevir and MK-5172. Furthermore, RIN and substrate envelope analysis indicate that the studied mutations of the residues are located outside the substrate (4B5A) binding site and selectively decrease the affinity of inhibitors but not the activity of the enzyme and consequently help NS3/4A protease escape from the effect of the inhibitors without influencing the affinity of substrate binding. These findings can provide useful information for understanding the drug resistance mechanism against vaniprevir and MK-5172. The results can also provide some potential clues for further design of novel inhibitors that are less susceptible to drug resistance.

Structural basis of the interactions between CXCR4 and CXCL12/SDF-1 revealed by theoretical approaches

The G protein-coupled chemokine receptor CXCR4 is implicated in a variety of physiological responses that also share several downstream effectors involved in multiple pathological processes. The interaction between CXCR4 and its natural ligand CXCL12/stromal-derived factor-1 (SDF-1) plays important roles in cancer metastasis, HIV-1 infection, and inflammatory diseases. Therefore, investigating the CXCR4-CXCL12 interaction is critical for understanding the molecular mechanisms of the modulation of chemokine-receptor functions and designing new pharmaceutical agents to target the CXCR4-CXCL12 pathway. Based on known experimental data, the interaction between CXCR4 and CXCL12 was predicted by an integrated protocol, which combines protein-protein docking, molecular dynamics (MD) simulations, Molecular Mechanics/Generalized Born Surface Area (MM/GBSA) binding free energy calculations, and MM/GBSA binding free energy decomposition analysis. The predicted CXCR4-CXCL12 binding pattern is in good agreement with the experimental data. Analysis of the binding structure reveals an obvious electrostatic complementarity between CXCR4 and CXCL12. Moreover, significant conformational rearrangements were observed during the 50 ns MD simulations. In particular, the basic Lys1 at the CXCL12 N-terminus, an essential residue in receptor activation, forms a strong polar interaction with the Glu32 in the CXCR4 extracellular region. It facilitates the significant movement of TM5 and TM6 in the conformational transition, which is coupled to the association with the intracellular signal transduction pathways via heterotrimer G-protein. Based on the dynamic and energetic analyses, a two-site binding model was proposed. We believe that our study provides useful information for understanding the mechanisms of CXCR4 ligand binding and structure-based drug design of CXCR4.

Combining molecular docking and molecular dynamics to predict the binding modes of flavonoid derivatives with the neuraminidase of the 2009 H1N1 influenza A virus

Control of flavonoid derivatives inhibitors release through the inhibition of neuraminidase has been identified as a potential target for the treatment of H1N1 influenza disease. We have employed molecular dynamics simulation techniques to optimize the 2009 H1N1 influenza neuraminidase X-ray crystal structure. Molecular docking of the compounds revealed the possible binding mode. Our molecular dynamics simulations combined with the solvated interaction energies technique was applied to predict the docking models of the inhibitors in the binding pocket of the H1N1 influenza neuraminidase. In the simulations, the correlation of the predicted and experimental binding free energies of all 20 flavonoid derivatives inhibitors is satisfactory, as indicated by R² = 0.75.