Molecular basis of the interaction for an essential subunit PA-PB1 in influenza virus RNA polymerase: insights from molecular dynamics simulation and free energy calculation

Inhibitory effect of thymoquinone from Nigella sativa against SARS-CoV-2 main protease. An in-silico study

Abstract Nigella sativa is known for the safety profile, containing a wealth of useful antiviral compounds. The main protease (Mpro, 3CLpro) of severe acute respiratory syndrome 2 (SARS-CoV-2) is being considered as one of the most attractive viral target, processing the polyproteins during viral pathogenesis and replication. In the current investigation we analyzed the potency of active component, thymoquinone (TQ) of Nigella sativa against SARS-CoV-2 Mpro. The structures of TQ and Mpro was retrieved from PubChem (CID10281) and Protein Data Bank (PDB ID 6MO3) respectively. The Mpro and TQ were docked and the complex was subjected to molecular dynamic (MD) simulations for a period 50ns. Protein folding effect was analyzed using radius of gyration (Rg) while stability and flexibility was measured, using root means square deviations (RMSD) and root means square fluctuation (RMSF) respectively. The simulation results shows that TQ is exhibiting good binding activity against SARS-CoV-2 Mpro, interacting many residues, present in the active site (His41, Cys145) and also the Glu166, facilitating the pocket shape. Further, experimental approaches are needed to validate the role of TQ against virus infection. The TQ is interfering with pocket maintaining residues as well as active site of virus Mpro which may be used as a potential inhibitor against SARS-CoV-2 for better management of COVID-19.

Elucidating the potential effects of point mutations on FGFR3 inhibitor resistance via combined molecular dynamics simulation and community network analysis

FGFR3 kinase mutations are associated with a variety of malignancies, but FGFR3 mutant inhibitors have rarely been studied. Furthermore, the mechanism of pan-FGFR inhibitors resistance caused by kinase domain mutations is still unclear. In this study, we try to explain the mechanism of drug resistance to FGFR3 mutation through global analysis and local analysis based on molecular dynamics simulation, binding free energy analysis, umbrella sampling and community network analysis. The results showed that FGFR3 mutations caused a decrease in the affinity between drugs and FGFR3 kinase, which was consistent with the reported experimental results. Possible mechanisms are that mutations affect drug-protein affinity by altering the environment of residues near the hinge region where the protein binds to the drug, or by affecting the A-loop and interfering with the allosteric communication networks. In conclusion, we systematically elucidated the underlying mechanism of pan-FGFR inhibitor resistance caused by FGFR3 mutation based on molecular dynamics simulation strategy, which provided theoretical guidance for the development of FGFR3 mutant kinase inhibitors.

In silico mutagenesis-based designing of oncogenic SHP2 peptide to inhibit cancer progression

Cancer is among the top causes of death, accounting for an estimated 9.6 million deaths in 2018, it appeared that approximately 500,000 people die from cancer in the United States alone annually. The SHP2 plays a major role in regulation of cell growth, proliferation, and differentiation, and functional upregulation of this enzyme is linked to oncogenesis and developmental disorders. SHP2 activity has been linked to several cancer types for which no drugs are currently available. In our study, we aimed to design peptide inhibitors against the SHP2 mutant. The crystal structure of the human Src SH2-PQpYEEIPI peptide mutant was downloaded from the protein databank. We generated several peptides from the native wild peptide using an in silico mutagenesis method, which showed that changes (P302W, Y304F, E306Q, and Q303A) might boost the peptide's affinity for binding to SHP2. Furthermore, the dynamical stability and binding affinities of the mutated peptide were confirmed using Molecular dynamics simulation and Molecular Mechanics with Generalized Born and Surface Area Solvation free energy calculations. The proposed substitution greatly enhanced the binding affinity at the residue level, according to a study that decomposed energy into its component residues. Our proposed peptide may prevent the spread of cancer by inhibiting SHP2, according to our detailed analyses of binding affinities.

Insight into non-nucleoside triazole-based systems as viral polymerases inhibitors

Viruses have been recognized as the etiological agents responsible for many pathological conditions ranging from asymptomatic infections to serious diseases, even leading to death. For this reason, many efforts have been made to identify selective viral targets with the aim of developing efficient therapeutic strategies, devoid of drug-resistance issues. Considering their crucial role in the viral life cycle, polymerases are very attractive targets. Among the classes of compounds explored as viral polymerases inhibitors, here we present an overview of non-nucleoside triazole-based compounds identified in the last fifteen years. Furthermore, the structure-activity relationships (SAR) of the different chemical entities are described in order to highlight the key chemical features required for the development of effective antiviral agents.

Molecular basis of ssDNA recognition by RBM45 protein of neurodegenerative disease from multiple molecular dynamics simulations and energy predictions

Neurodegenerative diseases (NDD) are a group of cognitive and behavioral disorders characterized by progressive loss of neuronal structure and function. As the population ages, the incidence is getting higher and higher, but there is currently no effective treatment. The details of RNA/DNA recognition by the RNA-binding protein RBM45 closely related to neurodegenerative diseases through its two tandem RNA-recognition domains at its N-terminus have important implications for structure-based drug discovery against degenerative diseases. To explore the key characteristics of interaction between ssDNA and RBM45, we performed multiple molecular dynamics (MD) simulations along with MM-PBSA energy prediction on the complexes of wild type (WT) and three mutant RBM45s (K100A, F124A/Y165A, and F29A/F70A/F124A/Y165A) with ssDNA, respectively. The findings suggest that these mutated residues of RBM45 modify the interaction of their surrounding residues with ssDNA, thereby affecting RBM45 protein binding to ssDNA. In contrast with WT RBM45 protein, variations in van der Waals and electrostatic interactions with ssDNA caused by these three RBM45 mutants are critical to affect binding between them. In addition, energy analysis showed that RBM45 is a specific ssDNA-binding protein. The results of our work provide valuable theoretical guidance for the design effective drugs of NDD.

Contemporary medicinal chemistry strategies for the discovery and optimization of influenza inhibitors targeting vRNP constituent proteins

Influenza is an acute respiratory infectious disease caused by the influenza virus, affecting people globally and causing significant social and economic losses. Due to the inevitable limitations of vaccines and approved drugs, there is an urgent need to discover new anti-influenza drugs with different mechanisms. The viral ribonucleoprotein complex (vRNP) plays an essential role in the life cycle of influenza viruses, representing an attractive target for drug design. In recent years, the functional area of constituent proteins in vRNP are widely used as targets for drug discovery, especially the PA endonuclease active site, the RNA-binding site of PB1, the cap-binding site of PB2 and the nuclear export signal of NP protein. Encouragingly, the PA inhibitor baloxavir has been marketed in Japan and the United States, and several drug candidates have also entered clinical trials, such as favipiravir. This article reviews the compositions and functions of the influenza virus vRNP and the research progress on vRNP inhibitors, and discusses the representative drug discovery and optimization strategies pursued.

Theoretical investigations on the effects of mutations in important residues of NS1B on its RNA-binding using molecular dynamics simulations

NS1B protein plays an important role in countering host antiviral defense and virulence of influenza virus B, considered as the promising target. The first experimental structure of the NS1B protein has recently been determined, was able to bind to double-stranded RNA (dsRNA). However, few studies attempt to investigate the RNA-binding mechanism of the NS1B. In this study, we provide our understanding of the structure-function relationship, dynamics and RNA-binding mechanism of the NS1B protein by performing molecular dynamics simulations combined and MM-GBSA calculations on the NS1B-dsRNA complex. 12 key residues are identified for RNA-binding by forming hydrogen bonds with the. Our results also demonstrate that mutations (R156A, K160A, R208A and K221A) can cause the local structure changes of NS1B CTD and the hydrogen bonds between NS1B CTD and RNA disappearance, which may be the main reasons for the decrease in RNA-binding affinity. These results mentioned will help us understanding the RNA-binding mechanism and could provide some medicinal chemistry insights chances for rational drug design targeting NS1B protein.

Lead Optimization of Influenza Virus RNA Polymerase Inhibitors Targeting PA-PB1 Interaction

Influenza viruses are responsible for contagious respiratory illnesses in humans and cause seasonal epidemics and occasional pandemics worldwide. Previously, we identified a quinolinone derivative PA-49, which inhibited the influenza virus RNA-dependent RNA polymerase (RdRp) by targeting PA-PB1 interaction. This paper reports the structure optimization of PA-49, which resulted in the identification of 3-((dibenzylamino)methyl)quinolinone derivatives with more potent anti-influenza virus activity. During the optimization, the hit compound 89, which was more active than PA-49, was identified. Further optimization and scaffold hopping of 89 led to the most potent compounds 100 and a 1,8-naphthyridinone derivative 118, respectively. We conclusively determined that compounds 100 and 118 suppressed the replication of influenza virus and exhibited anti-influenza virus activity against both influenza virus types A and B in the range of 50% effective concentration (EC₅₀) = 0.061-0.226 μM with low toxicity (50% cytotoxic concentration (CC₅₀) >10 μM).

Multiple Molecular Dynamics Simulations and Free-Energy Predictions Uncover the Susceptibility of Variants of HIV-1 Protease against Inhibitors Darunavir and KNI-1657

HIV-1 protease (PR) is considered to be the main targets of anti-AIDS drug design because of its role in the proteolytic processing of viral polyproteins. However, the emergence of drug-resistant HIV has become a major problem in the therapy of HIV-1-infected patients. Focused on the complexes of wild type (WT) PR and two mutant PRs (V32I/L33F/I54M/V82I and V32I/L33F/I54M/I84 V) with inhibitors Darunavir (DRV) and KNI-1657 (KNI), respectively, we have conducted research on the conformational dynamics and the resistance mechanism caused by residue mutations through multiple molecular dynamics (MD) simulations combined with an energy (MM-PBSA and solvated interaction energy (SIE)) prediction. The results indicate that mutated residues of PR alter the distance between flap regions and catalytic sites, the volume of the inner catalytic site, and the curling degree of the flap tips, thereby affecting DRV and KNI inhibitor binding to PR. These mutated residues reduced the binding affinity of the two mutant PRs to DRV, resulting in drug resistance, whereas the two mutant PRs increase the binding affinity with KNI, indicating they enhance the sensitivity to KNI. Compared with the WT PR, the changes in van der Waals interaction and electrostatic interaction in the two variant PRs play a vital part in the binding of PR with DRV and KNI. These results may supply valuable guidance for the design of anti-AIDS drugs targeting PR.

Emerging mutations in envelope protein of SARS-CoV-2 and their effect on thermodynamic properties

Structural proteins of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) are potential drug targets due to their role in the virus life cycle. The envelope (E) protein is one of the structural proteins; plays a critical role in virulency. However, the emergence of mutations oftenly leads to drug resistance and may also play a vital role in virus stabilization and evolution. In this study, we aimed to identify mutations in E proteins that affect the protein stability. About 0.3 million complete whole genome sequences were analyzed to screen mutations in E protein. All these mutations were subjected to stability prediction using the DynaMut server. The most common mutations that were detected at the C-terminal domain, Ser68Phe, Pro71Ser, and Leu73Phe, were examined through molecular dynamics (MD) simulations for a 100ns period. The sequence analysis shows the existence of 259 mutations in E protein. Interestingly, 16 of them were detected in the DFLV amino acid (aa) motif (aa72-aa75) that binds the host PALS1 protein. The results of root mean square deviation, fluctuations, radius of gyration, and free energy landscape show that Ser68Phe, Pro71Ser, and Leu73Phe are exhibiting a more stabilizing effect. However, a more comprehensive experimental study may be required to see the effect on virus pathogenicity. Potential antiviral drugs, and vaccines may be developed used after screening the genomic variations for better management of SARS-CoV-2 infections.

Inhibition of Influenza Virus Polymerase by Interfering with Its Protein-Protein Interactions

Influenza (flu) virus is a serious threat to global health with the potential to generate devastating pandemics. The availability of broad spectrum antiviral drugs is an unequaled weapon during pandemic events, especially when a vaccine is still not available. One of the most promising targets for the development of new antiflu therapeutics is the viral RNA-dependent RNA polymerase (RdRP). The assembly of the flu RdRP heterotrimeric complex from the individual polymerase acidic protein (PA), polymerase basic protein 1 (PB1), and polymerase basic protein 2 (PB2) subunits is a prerequisite for RdRP functions, such as mRNA synthesis and genome replication. In this Review, we report the known protein-protein interactions (PPIs) occurring by RdRP that could be disrupted by small molecules and analyze their benefits and drawbacks as drug targets. An overview of small molecules able to interfere with flu RdRP functions exploiting the PPI inhibition approach is described. In particular, an update on the most recent inhibitors targeting the well-consolidated RdRP PA-PB1 subunit heterodimerization is mainly reported, together with pioneer inhibitors targeting other virus-virus or virus-host interactions involving RdRP subunits. As demonstrated by the PA-PB1 interaction inhibitors discussed herein, the inhibition of flu RdRP functions by PPI disrupters clearly represents a valid means to identify compounds endowed with a broad spectrum of action and a reduced propensity to develop drug resistance, which are the main issues of antiviral drugs.

Triazolopyrimidine Nuclei: Privileged Scaffolds for Developing Antiviral Agents with a Proper Pharmacokinetic Profile

Viruses are a continuing threat to global health. The lack or limited therapeutic armamentarium against some viral infections and increasing drug resistance issues make the search for new antiviral agents urgent. In recent years, a growing literature highlighted the use of triazolopyrimidine (TZP) heterocycles in the development of antiviral agents, with numerous compounds that showed potent antiviral activities against different RNA and DNA viruses. TZP core represents a privileged scaffold for achieving biologically active molecules, thanks to: i) the synthetic feasibility that allows to variously functionalize TZPs in the different positions of the nucleus, ii) the ability of TZP core to establish multiple interactions with the molecular target, and iii) its favorable pharmacokinetic properties. In the present review, after mentioning selected examples of TZP-based compounds with varied biological activities, we will focus on those antivirals that appeared in the literature in the last 10 years. Approaches used for their identification, the hit-to-lead studies, and the emerged structure-activity relationship will be described. A mention of the synthetic methodologies to prepare TZP nuclei will also be given. In addition, their mechanism of action, the binding mode within the biological target, and pharmacokinetic properties will be analyzed, highlighting the strengths and weaknesses of compounds based on the TZP scaffold, which is increasingly used in medicinal chemistry.

Discovery of novel IDO1 inhibitors via structure-based virtual screening and biological assays

Indoleamine 2,3-dioxygenase 1 (IDO1) is a heme-containing enzyme that catalyzes the first and rate-limiting step in catabolism of tryptophan via the kynurenine pathway, which plays a pivotal role in the proliferation and differentiation of T cells. IDO1 has been proven to be an attractive target for many diseases, such as breast cancer, lung cancer, colon cancer, prostate cancer, etc. In this study, docking-based virtual screening and bioassays were conducted to identify novel inhibitors of IDO1. The cellular assay demonstrated that 24 compounds exhibited potent inhibitory activity against IDO1 at micromolar level, including 8 compounds with IC₅₀ values below 10 μM and the most potent one (compound 1) with IC₅₀ of 1.18 ± 0.04 μM. Further lead optimization based on similarity searching strategy led to the discovery of compound 28 as an excellent inhibitor with IC₅₀ of 0.27 ± 0.02 μM. Then, the structure-activity relationship of compounds 1, 2, 8 and 14 analogues is discussed. The interaction modes of two compounds against IDO1 were further explored through a Python Based Metal Center Parameter Builder (MCPB.py) molecular dynamics simulation, binding free energy calculation and electrostatic potential analysis. The novel IDO1 inhibitors of compound 1 and its analogues could be considered as promising scaffold for further development of IDO1 inhibitors.

Insights into RNA-dependent RNA Polymerase Inhibitors as Antiinfluenza Virus Agents

BACKGROUND

Influenza is a seasonal disease that affects millions of people every year and has a significant economic impact. Vaccines are the best strategy to fight this viral pathology, but they are not always available or administrable, prompting the search for antiviral drugs. RNA-dependent RNA polymerase (RdRp) recently emerged as a promising target because of its key role in viral replication and its high conservation among viral strains.

DISCUSSION

This review presents an overview of the most interesting RdRp inhibitors that have been discussed in the literature since 2000. Compounds already approved or in clinical trials and a selection of inhibitors endowed with different scaffolds are described, along with the main features responsible for their activity.

RESULTS

RdRp inhibitors are emerging as a new strategy to fight viral infections and the importance of this class of drugs has been confirmed by the FDA approval of baloxavir marboxil in 2018. Despite the complexity of the RdRp machine makes the identification of new compounds a challenging research topic, it is likely that in the coming years, this field will attract the interest of a number of academic and industrial scientists because of the potential strength of this therapeutic approach.

Revealing the Inhibition Mechanism of RNA-Dependent RNA Polymerase (RdRp) of SARS-CoV-2 by Remdesivir and Nucleotide Analogues: A Molecular Dynamics Simulation Study

Antiviral drug therapy against SARS-CoV-2 is not yet established and posing a serious global health issue. Remdesivir is the first antiviral compound approved by the US FDA for the SARS-CoV-2 treatment for emergency use, targeting RNA-dependent RNA polymerase (RdRp) enzyme. In this work, we have examined the action of remdesivir and other two ligands screened from the library of nucleotide analogues using docking and molecular dynamics (MD) simulation studies. The MD simulations have been performed for all the ligand-bound RdRp complexes for the 30 ns time scale. This is one of the earlier reports to perform the MD simulations studies using the SARS-CoV-2 RdRp crystal structure (PDB ID 7BTF). The MD trajectories were analyzed and Molecular Mechanics Poisson-Boltzmann Surface Area (MM-PBSA) calculations were performed to calculate the binding free energy. The binding energy data reveal that compound-17 (-59.6 kcal/mol) binds more strongly as compared to compound-8 (-46.3 kcal/mol) and remdesivir (-29.7 kcal/mol) with RdRp. The detailed analysis of trajectories shows that the remdesivir binds in the catalytic site and forms a hydrogen bond with the catalytic residues from 0 to 0.46 ns. Compound-8 binds in the catalytic site but does not form direct hydrogen bonds with catalytic residues. Compound-17 showed the formation of hydrogen bonds with catalytic residues throughout the simulation process. The MD simulation results such as hydrogen bonding, the center of mass distance analysis, snapshots at a different time interval, and binding energy suggest that compound-17 binds strongly with RdRp of SARS-CoV-2 and has the potential to develop as a new antiviral against COVID-19. Further, the frontier molecular orbital analysis and molecular electrostatic potential (MESP) iso-surface analysis using DFT calculations shed light on the superior binding of compound-17 with RdRp compared to remdesivir and compound-8. The computed as well as the experimentally reported pharmacokinetics and toxicity parameters of compound-17 is encouraging and therefore can be one of the potential candidates for the treatment of COVID-19.

Mechanistic analysis of A46V, H57Y, and D129N in pyrazinamidase associated with pyrazinamide resistance

Pyrazinamide (PZA) is a component of first-line drugs, active against latent Mycobacterium tuberculosis (MTB) isolates. The prodrug is activated into the active form, pyrazinoic acid (POA) via pncA gene-encoded pyrazinamidase (PZase). Mutations in pncA have been reported, most commonly responsible for PZA-resistance in more than 70% of the resistant cases. In our previous study, we detected many mutations in PZase among PZA-resistance MTB isolates including A46V, H71Y, and D129N. The current study was aimed to investigate the molecular mechanism of PZA-resistance behind mutants (MTs) A46V, H71Y, and D129N in comparison with the wild type (WT) through molecular dynamic (MD) simulation. MTB positive samples were subjected to PZA drug susceptibility testing (DST) against critical concentration (100ug/ml). The resistant samples were subjected to pncA sequencing. Thirty-six various mutations have been observed in the coding region of pncA of PZA-resistant isolates (GenBank accession No. MH461111) including A46V, H71Y, and D129N. The post-simulation analysis revealed a significant variation in MTs structural dynamics as compared to the WT. Root means square deviations (RMSD) and Root means square fluctuation (RMSF) has been found in variation between WT and MTs. Folding effect and pocket volume were altered in MTs when compared with WT. Geometric matching supports the effect of mutation A46V, H71Y, and D129N on PZase structure that may have an insight effect on PZase dynamics, making them vulnerable to convert pro-PZA into active form, POA. In conclusion, the current analyses will provide useful information behind PZA-resistance for better management of drug-resistant TB.

Synthesis and characterization of 1,2,4-triazolo[1,5-a]pyrimidine-2-carboxamide-based compounds targeting the PA-PB1 interface of influenza A virus polymerase

Influenza viruses (Flu) are responsible for seasonal epidemics causing high rates of morbidity, which can dramatically increase during severe pandemic outbreaks. Antiviral drugs are an indispensable weapon to treat infected people and reduce the impact on human health, nevertheless anti-Flu armamentarium still remains inadequate. In search for new anti-Flu drugs, our group has focused on viral RNA-dependent RNA polymerase (RdRP) developing disruptors of PA-PB1 subunits interface with the best compounds characterized by cycloheptathiophene-3-carboxamide and 1,2,4-triazolo[1,5-a]pyrimidine-2-carboxamide scaffolds. By merging these moieties, two very interesting hybrid compounds were recently identified, starting from which, in this paper, a series of analogues were designed and synthesized. In particular, a thorough exploration of the cycloheptathiophene-3-carboxamide moiety led to acquire important SAR insight and identify new active compounds showing both the ability to inhibit PA-PB1 interaction and viral replication in the micromolar range and at non-toxic concentrations. For few compounds, the ability to efficiently inhibit PA-PB1 subunits interaction did not translate into anti-Flu activity. Chemical/physical properties were investigated for a couple of compounds suggesting that the low solubility of compound 14, due to a strong crystal lattice, may have impaired its antiviral activity. Finally, computational studies performed on compound 23, in which the phenyl ring suitably replaced the cycloheptathiophene, suggested that, in addition to hydrophobic interactions, H-bonds enhanced its binding within the PA_C cavity.

Pyrazinamide resistance of novel mutations in pncA and their dynamic behavior

Pyrazinamide (PZA) is one of the essential anti-mycobacterium drugs, active against non-replicating Mycobacterium tuberculosis (MTB) isolates. PZA is converted into its active state, called pyrazinoic acid (POA), by action of pncA encoding pyrazinamidase (PZase). In the majority of PZA-resistance isolates, pncA harbored mutations in the coding region. In our recent report, we detected a number of novel variants in PZA-resistance (PZA^R) MTB isolates, whose resistance mechanisms were yet to be determined. Here we performed several analyses to unveil the PZA^R mechanism of R123P, T76P, G150A, and H71R mutants (MTs) through molecular dynamics (MD) simulations. In brief, culture positive MTB isolates were subjected to PZA susceptibility tests using the WHO recommended concentration of PZA (100 μg ml⁻¹). The PZA^R samples were screened for mutations in pncA along sensitive isolates through polymerase chain reactions and sequencing. A large number of variants (GeneBank accession no. MH461111), including R123P, T76P, G150A, and H71R, have been spotted in more than 70% of isolates. However, the mechanism of PZA^R for mutants (MTs) R123P, T76P, G150A, and H71R was unknown. For the MTs and native PZase structures (WT), thermodynamic properties were compared using molecular dynamics simulations for 100 ns. The MTs structural activity was compared to the WT. Folding effect and pocket volume variations have been detected when comparing between WT and MTs. Geometric matching further confirmed the effect of R123P, T76P, G150A, and H71R mutations on PZase dynamics, making them vulnerable for activating the pro-drug into POA. This study offers a better understanding for management of PZA^R TB. The results may be used as alternative diagnostic tools to infer PZA resistance at a structural dynamics level.

We performed several analyses to unveil the pyrazinamide-resistance mechanism of R123P, T76P, G150A, and H71R mutants through molecular dynamics simulations.

Bioinformatic and computational analysis for predominant mutations of the Nrf2/Keap1 complex in pediatric leukemia

The levels of reactive oxygen species (ROS) are tightly controlled and regulated by Nuclear Factor Erythroid-2-Like 2 (Nrf2) transcription factor, which is the main regulator of antioxidant responses and its suppressor protein Kelch-like ECH-associated protein 1 (Keap1). Our previous study has identified six novel changes in Nrf2/Keap1 pathway in pediatric ALL, which were described for the first time. These changes in the pathway are likely to alter the evolutionary process of amino acids and cause structural changes in the final products of genes. In this study, we aimed to compare the pathogenicity of eight determined mutations reported in our previous study by utilizing different programs with different algorithms and molecular dynamics simulation. Since it is too difficult to handle each existing mutation in a wet laboratory, in silico methods may give suggestion to choose the important mutations for further analysis and to establish the appropriate patient population and conduct wet laboratory studies. For this purpose, four different algorithms were used to evaluate the effects of single amino acid mutation. In addition, root-mean-square deviation, root-mean-square fluctuation and free-energy landscape analyses were performed to observe stability, flexibility and energetically favorable conformations, respectively, for each amino acid mutation. As a result, our study emphasizes the importance of Keap1 mutations in pediatric ALL Nrf2/Keap1 pathway, a total of eight mutations, two of which were shown for the first time in our study. Especially the mutations in the Keap1 Broad-Complex, Tramtrack and Bric-à-brac domain are worthy of attention.Communicated by Ramaswamy H. Sarma.

Gibbs Free Energy Calculation of Mutation in PncA and RpsA Associated With Pyrazinamide Resistance

A central approach for better understanding the forces involved in maintaining protein structures is to investigate the protein folding and thermodynamic properties. The effect of the folding process is often disturbed in mutated states. To explore the dynamic properties behind mutations, molecular dynamic (MD) simulations have been widely performed, especially in unveiling the mechanism of drug failure behind mutation. When comparing wild type (WT) and mutants (MTs), the structural changes along with solvation free energy (SFE), and Gibbs free energy (GFE) are calculated after the MD simulation, to measure the effect of mutations on protein structure. Pyrazinamide (PZA) is one of the first-line drugs, effective against latent Mycobacterium tuberculosis isolates, affecting the global TB control program 2030. Resistance to this drug emerges due to mutations in pncA and rpsA genes, encoding pyrazinamidase (PZase) and ribosomal protein S1 (RpsA) respectively. The question of how the GFE may be a measure of PZase and RpsA stabilities, has been addressed in the current review. The GFE and SFE of MTs have been compared with WT, which were already found to be PZA-resistant. WT structures attained a more stable state in comparison with MTs. The physiological effect of a mutation in PZase and RpsA may be due to the difference in energies. This difference between WT and MTs, depicted through GFE plots, might be useful in predicting the stability and PZA-resistance behind mutation. This study provides useful information for better management of drug resistance, to control the global TB problem.

Molecular Simulation Elaborating the Mechanism of 1β-Hydroxy Alantolactone Inhibiting Ubiquitin-Conjugating Enzyme UbcH5s

1β-hydroxy alantolactone, a sesquiterpene lactone, exhibits potent anti-inflammatory and anticancer activities. Recently, it has been found to target UbcH5s by covalently bonding with Cys85 specifically, but the exact molecular basis remains unclear. Here, we analyzed the structural specificity of the catalytic site of UbcH5s by comparing them with other E2 proteins. Molecular dynamics was performed to detect the structural stability of the catalytic site. Docking method was then used to predict conformations of ligand docked at the catalytic site of UbcH5s. The electrostatic surface and charge distribution of ligand and proteins were analyzed by quantitative calculation. Molecular dynamics was used to detect the stability of docking complexes of 1β-hydroxy alantolactone and UbcH5s, the covalently bonded intermediates and the products. The QM/MM methodology was used to calculate the free energy barrier of hydrogen transfer and formation of covalent bond between 15-position carbon of ligand and Cys85. Results revealed that the structure of the catalytic site is stable, and 1β-hydroxy alantolactone can dock at the catalytic site with correct conformation. Molecular dynamics further demonstrates that 1β-hydroxy alantolactone can steadily combine with UbcH5s. Intermediate and product of catalytic reaction are also certified to be stable. Besides, Asp112 and Asn114 function as anchors to fix ligand, ensuring it steadily docked at catalytic site to complete covalent reaction. More importantly, we have found that Cys85 of UbcH5c is more efficient to form a covalent bond with the ligand in comparison with UbcH5a and UbcH5b. Our results successfully explained the mechanism of 1β-hydroxy alantolactone covalently bonding with UbcH5s. Such molecular mechanism may provide a better insight into the molecular development or modification for ubiquitin-related drugs.

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.

Structural Dynamics Behind Clinical Mutants of PncA-Asp12Ala, Pro54Leu, and His57Pro of Mycobacterium tuberculosis Associated With Pyrazinamide Resistance

Pyrazinamide (PZA) is one of the main FDA approved drugs to be used as the first line of defense against Mycobacterium Tuberculosis (MTB). It is activated into pyrazinoic acid (POA) via MTB's pncA gene-encoded pyrazinamidase (PZase). Mutations are most commonly responsible for PZA-resistance in nearly 70% of the resistant samples. In the present work, MTB positive samples were chosen for PZA drug susceptibility testing (DST) against critical concentration (100 ug/ml) of PZA. The resistant samples were subjected to pncA sequencing. As a result, 36 various mutations have been observed in the PZA resistant samples, uploaded to the NCBI (GeneBank accession no. MH461111). Here we report the mechanism of PZA resistance behind the three mutants (MTs), Asp12Ala, Pro54Leu, and His57Pro in comparison with the wild type (WT) through molecular dynamics simulation to unveil how these mutations affect the overall conformational stability. The post-simulation analyses revealed notable deviations as compared to the WT structure. Molecular docking studies of PZA with MTs and WT, pocket volume inspection and overall shape complementarity analysis confirmed the deleterious nature of these mutations and gave an insight into the mechanism behind PZA-resistance. These analyses provide vital information regarding MTB drug resistance and could be extremely useful in therapy management and overcoming its global burden.

Structural dynamics behind variants in pyrazinamidase and pyrazinamide resistance

Pyrazinamide (PZA) is an important component of first-line anti-tuberculosis (anti-TB) drugs. The anti-TB agent is activated into an active form, pyrazinoic acid (POA), by Mycobacterium tuberculosis (MTB) pncA gene encoding pyrazinamidase (PZase). The major cause of PZA-resistance has been associated with mutations in the pncA gene. We have detected several novel mutations including V131F, Q141P, R154T, A170P, and V180F (GeneBank Accession No. MH461111) in the pncA gene of PZA-resistant isolates during PZA drug susceptibility testing followed by pncA gene sequencing. Here, we investigated molecular mechanism of PZA-resistance by comparing the results of experimental and molecular dynamics. The mutants (MTs) and wild type (WT) PZase structures in apo and complex with PZA were subjected to molecular dynamic simulations (MD) at the 40 ns. Multiple factors, including root mean square deviations (RMSD), binding pocket, total energy, dynamic cross correlation, and root mean square fluctuations (RMSF) of MTs and WT were compared. The MTs attained a high deviation and fluctuation compared to WT. Binding pocket volumes of the MTs, were found, lower than the WT, and the docking scores were high than WT while shape complementarity scores were lower than that of the WT. Residual motion in MTs are seemed to be dominant in anti-correlated motion. Mutations at locations, V131F, Q141P, R154T, A170P, and V180F, might be involved in the structural changes, possibly affecting the catalytic property of PZase to convert PZA into POA. Our study provides useful information that will enhance the understanding for better management of TB. AbbreviationsDSTdrug susceptibility testingΔelecelectrostatic energyLJLowenstein-Jensen mediumMGITmycobacterium growth indicator tubesMTsmutantsMDmolecular dynamic simulationsMTBMycobacterium tuberculosisNALC-NaOHN-acetyl-l-cysteine-sodium hydroxideNIHNational Institutes of HealthNPTamount of substance (N), pressure (P) temperature (T)NVTmoles (N), volume (V) temperature (T)PZasepyrazinamidaseΔpspolar solvation energyPTRLProvincial Tuberculosis Reference LaboratoryRMSDroot mean square deviationsRMSFroot mean square fluctuationsΔSASAsolvent accessible surface area energyTBtuberculosisGTotaltotal binding free energyΔvdWVan der Waals energyWTwild typeCommunicated by Ramaswamy H. Sarma.

Structural and free energy landscape of novel mutations in ribosomal protein S1 (rpsA) associated with pyrazinamide resistance

Resistance to key first-line drugs is a major hurdle to achieve the global end tuberculosis (TB) targets. A prodrug, pyrazinamide (PZA) is the only drug, effective in latent TB, recommended in drug resistance and susceptible Mycobacterium tuberculosis (MTB) isolates. The prodrug conversion into active form, pyrazinoic acid (POA), required the activity of pncA gene encoded pyrazinamidase (PZase). Although pncA mutations have been commonly associated with PZA resistance but a small number of resistance cases have been associated with mutationss in RpsA protein. Here in this study a total of 69 PZA resistance isolates have been sequenced for pncA mutations. However, samples that were found PZA resistant but pncA wild type (pncAWT), have been sequenced for rpsA and panD genes mutation. We repeated a drug susceptibility testing according to the WHO guidelines on 18 pncAWT MTB isolates. The rpsA and panD genes were sequenced. Out of total 69 PZA resistant isolates, 51 harbored 36 mutations in pncA gene (GeneBank Accession No. MH46111) while, fifteen different mutations including seven novel, were detected in the fourth S1 domain of RpsA known as C-terminal (MtRpsACTD) end. We did not detect any mutations in panD gene. Among the rpsA mutations, we investigated the molecular mechanism of resistance behind mutations, D342N, D343N, A344P, and I351F, present in the MtRpsACTD through molecular dynamic simulations (MD). WT showed a good drug binding affinity as compared to mutants (MTs), D342N, D343N, A344P, and I351F. Binding pocket volume, stability, and fluctuations have been altered whereas the total energy, protein folding, and geometric shape analysis further explored a significant variation between WT and MTs. In conclusion, mutations in MtRpsACTD might be involved to alter the RpsA activity, resulting in drug resistance. Such molecular mechanism behind resistance may provide a better insight into the resistance mechanism to achieve the global TB control targets.

Covalent vs. Non-Covalent Inhibition: Tackling Drug Resistance in EGFR - A Thorough Dynamic Perspective

A persistent challenge in the treatment of non-small cell lung cancer (NSCLC) with EGFR is the emergence of drug-resistant caused by somatic mutations. The EGFR L858R/T790 M double mutant (EGFR^DM ) was found to be the most alarming variant. Despite the development of a wide range of inhibitors, none of them could inhibit EGFR^DM effectively. Recently, 11h and 45a, have been found to be potent inhibitors against EGFR^DM through two distinctive mechanisms, non-covalent and covalent binding, respectively. However, the structural and dynamic implications of the two modes of inhibitions remain unexplored. Herein, two molecular dynamics simulation protocols, coupled with free-energy calculations, were applied to gain insight into the atomistic nature of each binding mode. The comparative analysis confirmed that there is a significant difference in the binding free energy between 11h and 45a (ΔΔG_bind =-21.17 kcal/mol). The main binding force that governs the binding of both inhibitors is vdW, with a higher contribution for 45a. Two residues ARG841 and THR854 were found to have curtailed role in the binding of 45a to EGFR^DM by stabilizing its flexible alcohol chain. The 45a binding to EGFR^DM induces structural rearrangement in the active site to allow easier accessibility of 45a to target residue CYS797. The findings of this work can substantially shed light on new strategies for developing novel classes of covalent and non-covalent inhibitors with increased specificity and potency.

Understanding the selectivity of inhibitors toward PI4KIIIα and PI4KIIIβ based molecular modeling

Molecular dynamics simulations and binding free energy calculations are combined to investigate the selectivity of inhibitors toward type III phosphatidylinositol 4 kinases.

The molecular mechanism of pH-regulating C3d-CR2 interactions: Insights from molecular dynamics simulation

The interactions of complement receptor 2 (CR2) and the degradation fragment C3d of complement component C3 mediate the innate and adaptive immune systems. Due to the importance of C3d-CR2 interaction in the design of vaccines, many studies have indicated the interactions are pH-dependent. Moreover, C3d-CR2 interactions at pH 5.0 are unknown. To investigate the molecular mechanism of pH-regulating C3d-CR2 interaction, molecular dynamics simulations for C3d-CR2 complex in different pH are performed. Our results revealed that the protonation of His9 in C3d at pH 6.0 slightly weakens C3d-CR2 association as reducing pH from 7.4 to 6.0, initiated from a key hydrogen bond formed between Gly270 and His9 in C3d at pH 6.0. When reducing pH from 6.0 to 5.0, the protonation of His33 in C3d weakens C3d-SCR1 association by changing the hydrogen-bond network of Asp36, Glu37, and Glu39 in C3d with Arg13 in CR2. In addition, the protonation of His90 significantly enhances C3d-SCR2 association. This is because the enhanced hydrogen-bond interactions of His90 with Glu63 and Ser69 of the linker change the conformations of the linker, Cys112-Asn116 and Pro87-Gly91 regions. This study uncovers the molecular mechanism of the mediation of pH on C3d-CR2 interaction, which is valuable for vaccine design.

Prediction of GluN2B-CT1290-1310/DAPK1 Interaction by Protein⁻Peptide Docking and Molecular Dynamics Simulation

The interaction of death-associated protein kinase 1 (DAPK1) with the 2B subunit (GluN2B) C-terminus of N-methyl-D-aspartate receptor (NMDAR) plays a critical role in the pathophysiology of depression and is considered a potential target for the structure-based discovery of new antidepressants. However, the 3D structures of C-terminus residues 1290⁻1310 of GluN2B (GluN2B-CT_1290-1310) remain elusive and the interaction between GluN2B-CT_1290-1310 and DAPK1 is unknown. In this study, the mechanism of interaction between DAPK1 and GluN2B-CT_1290-1310 was predicted by computational simulation methods including protein⁻peptide docking and molecular dynamics (MD) simulation. Based on the equilibrated MD trajectory, the total binding free energy between GluN2B-CT_1290-1310 and DAPK1 was computed by the mechanics generalized born surface area (MM/GBSA) approach. The simulation results showed that hydrophobic, van der Waals, and electrostatic interactions are responsible for the binding of GluN2B-CT_1290⁻1310/DAPK1. Moreover, through per-residue free energy decomposition and in silico alanine scanning analysis, hotspot residues between GluN2B-CT_1290-1310 and DAPK1 interface were identified. In conclusion, this work predicted the binding mode and quantitatively characterized the protein⁻peptide interface, which will aid in the discovery of novel drugs targeting the GluN2B-CT_1290-1310 and DAPK1 interface.

Synthesis and biological evaluation of a library of hybrid derivatives as inhibitors of influenza virus PA-PB1 interaction

The limited treatment options against influenza virus along with the growing public health concerns regarding the continuous emergence of drug-resistant viruses make essential the development of new anti-flu agents with novel mechanisms of action. One of the most attractive targets is the interaction between two subunits of the RNA-dependent RNA polymerase, PA and PB1. Herein we report the rational design of hybrid compounds starting from a 3-cyano-4,6-diphenylpyridine scaffold recently identified as disruptor of PA-PB1 interactions. Guided by the previously reported SAR data, a library of amino acid derivatives was synthesized. The biological evaluation led to the identification of new PA-PB1 inhibitors, that do not show appreciable toxicity. Molecular modeling shed further lights on the inhibition mechanism of these compounds.

Influenza A virus polymerase: an attractive target for next-generation anti-influenza therapeutics

The influenza RNA-dependent RNA polymerase (RdRP) is conserved among different types of influenza virus, playing an important part in transcription and replication. In this regard, influenza RdRP is an attractive target for novel anti-influenza drug discovery. Herein, we will introduce the structural and functional information of influenza polymerase; and an overview of inhibitors targeting the PA endonuclease and PB2 cap-binding site is provided, along with the approaches utilized for identification of these inhibitors. The protein-protein interactions (PPIs) of the three polymerase subunits: PA, PB1 and PB2, are described based on the published crystal structures, and inhibitors targeting the PA-PB1 interaction are introduced briefly.

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 Study on the Binding and Stabilization Mechanism of Antiprion Compounds to the "Hot Spot" Region of PrPC

Structural transitions in the prion protein from the cellular form, PrP^C, into the pathological isoform, PrP^Sc, are regarded as the main cause of the transmissible spongiform encephalopathies, also known as prion diseases. Hence, discovering and designing effective antiprion drugs that can inhibit PrP^C to PrP^Sc conversion is regarded as a promising way to cure prion disease. Among several strategies to inhibit PrP^C to PrP^Sc conversion, stabilizing the native PrP^C via specific binding is believed to be one of the valuable approaches and many antiprion compounds have been reported based on this strategy. However, the detailed mechanism to stabilize the native PrP^C is still unknown. As such, to unravel the stabilizing mechanism of these compounds to PrP^C is valuable for the further design and discovery of antiprion compounds. In this study, by molecular dynamics simulation method, we investigated the stabilizing mechanism of several antiprion compounds on PrP^C that were previously reported to have specific binding to the "hot spot" region of PrP^C. Our simulation results reveal that the stabilization mechanism of specific binding compounds can be summarized as (I) to stabilize both the flexible C-terminal of α2 and the hydrophobic core, such as BMD42-29 and GN8; (II) to stabilize the hydrophobic core, such as J1 and GJP49; (III) to stabilize the overall structure of PrP^C by high binding affinity, as NPR-056. In addition, as indicated by the H-bond analysis and decomposition analysis of binding free energy, the residues N159 and Q160 play an important role in the specific binding of the studied compounds and all these compounds interact with PrP^C in a similar way with the key interacting residues L130 in the β1 strand, P158, N159, Q160, etc. in the α1-β2 loop, and H187, T190, T191, etc. in the α2 C-terminus although the compounds have large structural difference. As a whole, our obtained results can provide some insights into the specific binding mechanism of main antiprion compounds to the "hot spot" region of PrP^C at the molecular level and also provide guidance for effective antiprion drug design in the future.

Host proteostasis modulates influenza evolution

Predicting and constraining RNA virus evolution require understanding the molecular factors that define the mutational landscape accessible to these pathogens. RNA viruses typically have high mutation rates, resulting in frequent production of protein variants with compromised biophysical properties. Their evolution is necessarily constrained by the consequent challenge to protein folding and function. We hypothesized that host proteostasis mechanisms may be significant determinants of the fitness of viral protein variants, serving as a critical force shaping viral evolution. Here, we test that hypothesis by propagating influenza in host cells displaying chemically-controlled, divergent proteostasis environments. We find that both the nature of selection on the influenza genome and the accessibility of specific mutational trajectories are significantly impacted by host proteostasis. These findings provide new insights into features of host–pathogen interactions that shape viral evolution, and into the potential design of host proteostasis-targeted antiviral therapeutics that are refractory to resistance.

Influenza viruses, commonly called flu, can evade our immune system and develop resistance to treatments by changing frequently. Specifically, mutations in their genome cause influenza proteins to change in ways that can help the virus evade our defences. However, these mutations come at a cost and can prevent the viral proteins from forming functional and stable three-dimensional shapes – a process known as protein folding – thereby hampering the virus’ ability to replicate.

In human cells, proteins called chaperones can help our other proteins fold properly. Influenza viruses do not have their own chaperones and, instead, hijack those of their host. Host chaperones are therefore crucial to the virus’ ability to replicate. However, until now, it was not known if host chaperones can influence how these viruses evolve.

Here, Phillips et al. used mammalian cells to study how host chaperones affect an evolving influenza population. First, cells were engineered to either have normal chaperone levels, elevated chaperone levels, or inactive chaperones. Next, the H3N2 influenza strain was grown in these different conditions for nearly 200 generations and sequenced to determine how the virus evolved in each distinctive host chaperone environment.

Phillips et al. discovered that host chaperones affect the rate at which mutations accumulate in the influenza population, and also the types of mutations in the influenza genome. For instance, when a chaperone called Hsp90 was inactivated, mutations became prevalent in the viral population more slowly than in cells with normal or elevated chaperone levels. Moreover, some specific mutations fared better in cells with high chaperone levels, whilst others worked better in cells with inactivated chaperones.

These results suggest that influenza evolution is affected by host chaperone levels in complex and important ways. Moreover, whether chaperones will promote or hinder the effects of any single mutation is difficult to predict ahead of time. This discovery is significant, as the chaperones available to influenza can vary in different tissues, organisms and infectious conditions, and may therefore influence the virus' ability to change and evolve in a context-specific manner.

The findings are likely to extend to other viruses such as HIV and Ebola, which also hijack host chaperones for the same purpose. More work is now needed to systematically quantify these effects so that we can better predict how specific chaperones will affect the ability of viruses to adapt, especially in pathologically relevant conditions like fever or viral host-switching. In the future, such insights could help shape the design of treatments to which viruses do not evolve resistance.

DNA binding strength increases the processivity and activity of a Y-Family DNA polymerase

DNA polymerase (pol) processivity, i.e., the bases a polymerase extends before falling off the DNA, and activity are important for copying difficult DNA sequences, including simple repeats. Y-family pols would be appealing for copying difficult DNA and incorporating non-natural dNTPs, due to their low fidelity and loose active site, but are limited by poor processivity and activity. In this study, the binding between Dbh and DNA was investigated to better understand how to rationally design enhanced processivity in a Y-family pol. Guided by structural simulation, a fused pol Sdbh with non-specific dsDNA binding protein Sso7d in the N-terminus was designed. This modification increased in vitro processivity 4-fold as compared to the wild-type Dbh. Additionally, bioinformatics was used to identify amino acid mutations that would increase stabilization of Dbh bound to DNA. The variant SdbhM76I further improved the processivity of Dbh by 10 fold. The variant SdbhKSKIP241–245RVRKS showed higher activity than Dbh on the incorporation of dCTP (correct) and dATP (incorrect) opposite the G (normal) or 8-oxoG(damaged) template base. These results demonstrate the capability to rationally design increases in pol processivity and catalytic efficiency through computational DNA binding predictions and the addition of non-specific DNA binding domains.

Molecular Simulation Studies on the Binding Selectivity of Type-I Inhibitors in the Complexes with ROS1 versus ALK

ROS1 and ALK are promising targets of anticancer drugs for non-small-cell lung cancer. Since they have 49% amide acid sequence homology in the kinases domain and 77% identity at the ATP binding area, some ALK inhibitors also showed some significant responses for ROS1 in the clinical trial, such as the type-I binding inhibitor crizotinib and PF-06463922. As a newly therapeutic target, the selective ROS1 inhibitor is relatively rare. Moreover, the molecular basis for the selectivity of ROS1 versus ALK still remains unclear. In order to disclose the binding preference toward ROS1 over ALK and to aid the design of selective ROS1 inhibitors, the specific interactions and difference of conformational changes in the dual and selective ROS1/ALK inhibitors systems were investigated by molecular dynamics (MD) simulation and principle component analysis (PCA) in our work. Afterward, binding free energies (MM/GBSA) and binding free energies decomposition analysis indicated that the dominating effect of Van der Waals interaction drives the specific binding process of the type-I inhibitor, and residues of the P-loop and the DFG motif would play an important role in selectivity. On the basis of the modeling results, the new designed compound 14c was verified as a selective ROS1 inhibitor versus ALK, and SMU-B was a dual ROS1/ALK inhibitor by the kinase inhibitory study. These results are expected to facilitate the discovery and rational design of novel and specific ROS1 inhibitors.

Polymerase Acidic Protein-Basic Protein 1 (PA-PB1) Protein-Protein Interaction as a Target for Next-Generation Anti-influenza Therapeutics

The limited therapeutic options against the influenza virus (flu) and increasing challenges in drug resistance make the search for next-generation agents imperative. In this context, heterotrimeric viral PA/PB1/PB2 RNA-dependent RNA polymerase is an attractive target for a challenging but strategic protein-protein interaction (PPI) inhibition approach. Since 2012, the inhibition of the polymerase PA-PB1 subunit interface has become an active field of research following the publication of PA-PB1 crystal structures. In this Perspective, we briefly discuss the validity of flu polymerase as a drug target and its inhibition through a PPI inhibition strategy, including a comprehensive analysis of available PA-PB1 structures. An overview of all of the reported PA-PB1 complex formation inhibitors is provided, and approaches used for identification of the inhibitors, the hit-to-lead studies, and the emerged structure-activity relationship are described. In addition to highlighting the strengths and weaknesses of all of the PA-PB1 heterodimerization inhibitors, we analyze their hypothesized binding modes and alignment with a pharmacophore model that we have developed.

Assessing the performance of MM/PBSA and MM/GBSA methods. 4. Accuracies of MM/PBSA and MM/GBSA methodologies evaluated by various simulation protocols using PDBbind data set

By using different evaluation strategies, we systemically evaluated the performance of Molecular Mechanics/Generalized Born Surface Area (MM/GBSA) and Molecular Mechanics/Poisson-Boltzmann Surface Area (MM/PBSA) methodologies based on more than 1800 protein-ligand crystal structures in the PDBbind database. The results can be summarized as follows: (1) for the one-protein-family/one-binding-ligand case which represents the unbiased protein-ligand complex sampling, both MM/GBSA and MM/PBSA methodologies achieve approximately equal accuracies at the interior dielectric constant of 4 (with rp = 0.408 ± 0.006 of MM/GBSA and rp = 0.388 ± 0.006 of MM/PBSA based on the minimized structures); while for the total dataset (1864 crystal structures), the overall best Pearson correlation coefficient (rp = 0.579 ± 0.002) based on MM/GBSA is better than that of MM/PBSA (rp = 0.491 ± 0.003), indicating that biased sampling may significantly affect the accuracy of the predicted result (some protein families contain too many instances and can bias the overall predicted accuracy). Therefore, family based classification is needed to evaluate the two methodologies; (2) the prediction accuracies of MM/GBSA and MM/PBSA for different protein families are quite different with rp ranging from 0 to 0.9, whereas the correlation and ranking scores (an averaged rp/rs over a list of protein folds and also representing the unbiased sampling) given by MM/PBSA (rp-score = 0.506 ± 0.050 and rs-score = 0.481 ± 0.052) are comparable to those given by MM/GBSA (rp-score = 0.516 ± 0.047 and rs-score = 0.463 ± 0.047) at the fold family level; (3) for the overall prediction accuracies, molecular dynamics (MD) simulation may not be quite necessary for MM/GBSA (rp-minimized = 0.579 ± 0.002 and rp-1ns = 0.564 ± 0.002), but is needed for MM/PBSA (rp-minimized = 0.412 ± 0.003 and rp-1ns = 0.491 ± 0.003). However, for the individual systems, whether to use MD simulation is depended. (4) both MM/GBSA and MM/PBSA may be unable to give successful predictions for the ligands with high formal charges, with the Pearson correlation coefficient ranging from 0.621 ± 0.003 (neutral ligands) to 0.125 ± 0.142 (ligands with a formal charge of 5). Therefore, it can be summarized that, although MM/GBSA and MM/PBSA perform similarly in the unbiased dataset, for the currently available crystal structures in the PDBbind database, compared with MM/GBSA, which may be used in multi-target comparisons, MM/PBSA is more sensitive to the investigated systems, and may be more suitable for individual-target-level binding free energy ranking. This study may provide useful guidance for the post-processing of docking based studies.

Molecular dynamics, monte carlo simulations, and langevin dynamics: a computational review

Macromolecular structures, such as neuraminidases, hemagglutinins, and monoclonal antibodies, are not rigid entities. Rather, they are characterised by their flexibility, which is the result of the interaction and collective motion of their constituent atoms. This conformational diversity has a significant impact on their physicochemical and biological properties. Among these are their structural stability, the transport of ions through the M2 channel, drug resistance, macromolecular docking, binding energy, and rational epitope design. To assess these properties and to calculate the associated thermodynamical observables, the conformational space must be efficiently sampled and the dynamic of the constituent atoms must be simulated. This paper presents algorithms and techniques that address the abovementioned issues. To this end, a computational review of molecular dynamics, Monte Carlo simulations, Langevin dynamics, and free energy calculation is presented. The exposition is made from first principles to promote a better understanding of the potentialities, limitations, applications, and interrelations of these computational methods.

Computational insights into the inhibition and destabilization of morin on the oligomer of full-length human islet amyloid polypeptide

In this work, we simulated the full-length human islet amyloid peptide (hIAPP) pentamer with and without morins to investigate the mechanism of inhibition and destabilization of this inhibitor on hIAPP oligomer, and identify its possible binding sites on hIAPP.

Importance of protein flexibility in ranking inhibitor affinities: modeling the binding mechanisms of piperidine carboxamides as Type I1/2 ALK inhibitors

Anaplastic lymphoma kinase (ALK) has gained increased attention as an attractive therapeutic target for the treatment of various cancers, especially non-small-cell lung cancer (NSCLC).

Molecular mechanism of the inhibition and remodeling of human islet amyloid polypeptide (hIAPP(1-37)) oligomer by resveratrol from molecular dynamics simulation

Natural polyphenols are one of the most actively investigated categories of amyloid inhibitors, and resveratrol has recently been reported to inhibit and remodel the human islet amyloid polypeptide (hIAPP) oligomers and fibrils. However, the exact mechanism of its action is still unknown, especially for the full-length hIAPP1-37. To this end, we performed all-atom molecular dynamics simulations for hIAPP1-37 pentamer with and without resveratrol. The obtained results show that the binding of resveratrol is able to cause remarkable conformational changes of hIAPP1-37 pentamer, in terms of secondary structures, order degree, and morphology. By clustering analysis, two possible binding sites of resveratrol on the hIAPP1-37 pentamer were found, located at the grooves of the top and bottom surfaces of β-sheet layer, respectively. After the binding free energy calculation and residue energy decomposition, it can be concluded that the bottom site is the more possible one, and that the nonpolar interactions act as the driving force for the binding of hIAPP1-37 to resveratrol. In addition, Arg11 is the most important residue for the binding of resveratrol. The full understanding of inhibitory mechanism of resveratrol on the hIAPP1-37 oligomer, and the identification of its binding sites on this protein are helpful for the future design and discovery of new amyloid inhibitors.