d-lactate modulates M2 tumor-associated macrophages and remodels immunosuppressive tumor microenvironment for hepatocellular carcinoma.

The polarization of tumor-associated macrophages (TAMs) from M2 to M1 phenotype demonstrates great potential for remodeling the immunosuppressive tumor microenvironment (TME) of hepatocellular carcinoma (HCC). d-lactate (DL; a gut microbiome metabolite) acts as an endogenous immunomodulatory agent that enhances Kupffer cells for clearance of pathogens. In this study, the potential of DL for transformation of M2 TAMs to M1 was confirmed, and the mechanisms underlying such polarization were mainly due to the modulation of phosphatidylinositol 3-kinase/protein kinase B pathway. A poly(lactide-co-glycolide) nanoparticle (NP) was used to load DL, and the DL-loaded NP was modified with HCC membrane and M2 macrophage-binding peptide (M2pep), forming a nanoformulation (DL@NP-M-M2pep). DL@NP-M-M2pep transformed M2 TAMs to M1 and remodeled the immunosuppressive TME in HCC mice, promoting the efficacy of anti-CD47 antibody for long-term animal survival. These findings reveal a potential TAM modulatory function of DL and provide a combinatorial strategy for HCC immunotherapy.

How does multiple substrate binding lead to substrate inhibition of CYP2D6 metabolizing dextromethorphan? A theoretical study.

CYP2D6 is one of the most important metalloenzymes involved in the biodegradation of many drug molecules in the human body. It has been found that multiple substrate binding can lead to substrate inhibition of CYP2D6 metabolizing dextromethorphan (DM), but the corresponding theoretical mechanism is rarely reported. Therefore, we chose DM as the probe and performed molecular dynamics simulations and quantum mechanical calculations on CYP2D6-DM systems to investigate the mechanism of how the multiple substrate binding leads to the substrate inhibition of CYP2D6 metabolizing substrates. According to our results, three gate residues (Arg221, Val374, and Phe483) for the catalytic pocket are determined. We also found that the multiple substrate binding can lead to substrate inhibition by reducing the stability of CYP2D6 binding DM and increasing the reactive activation energy of the rate-determining step. Our findings would help to understand the substrate inhibition of CYP2D6 metabolizing the DM and enrich the knowledge of the drug-drug interactions for the cytochrome P450 superfamily.

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.

The regioselectivity of the interaction between dextromethorphan and CYP2D6.

CYP2D6 is an important enzyme of the cytochrome P450 superfamily, and catalyzes nearly 25% of the drugs sold in the market. For decades, the interactions and metabolism between CYP2D6 and substrates have been a hot topic. However, the key factors of the catalytic regioselectivity for CYP2D6 still remain controversial. Here, we construct four systems to explore the interaction between dextromethorphan (DM) and CYP2D6. A new binding mode of CYP2D6 is defined, and two key residues (residue Asp301 and residue Glu216) are discovered working simultaneously to stabilize the DM at the reactive site by forming water bridge hydrogen bonds when CYP2D6 binds DM. Our results also indicate that the substrate concentration could mediate the binding mode between the substrate and CYP2D6 by decreasing the volume of the catalytic pocket, which is not conducive to the O-demethylation of DM but benefits the N-demethylation of DM. These results could shed light on the process of CYP2D6 binding to the substrate, and help to better understand the regioselectivity of CYP2D6 catalyzing the substrates.

Predicting a Kind of Unusual Multiple-States Dimerization-Modes Transformation in Protein PD-L1 System by Computational Investigation and a Generalized Rate Theory.

The new cancer immunotherapy has been carried out with an almost messianic zeal, but its molecular basis remains unclear due to the complexity of programmed death ligand 1 (PD-L1) dimerization. In this study, a new and integral multiple dimerization-modes transformation process of PD-L1s (with a new PD-L1 dimerization mode and a new transformation path discovered) and the corresponding mechanism are predicted using theoretical and computational methods. The results of the state analysis show that 5 stable binding states exist in system. A generalized inter-state transformation rate (GITR) theory is also proposed in such multiple-states self-assembly system to explore the kinetic characteristics of inter-state transformation. A "drug insertion" path was identified as the dominant path of the PD-L1 dimerization-modes transformation. Above results can provide supports for both the relative drug design and other multiple-states self-assembly system from the theoretical chemistry perspective.

Exploring the Distinct Binding and Activation Mechanisms for Different CagA Oncoproteins and SHP2 by Molecular Dynamics Simulations.

CagA is a major virulence factor of Helicobacter pylori. H. pylori CagA is geographically subclassified into East Asian CagA and Western CagA, which are characterized by the presence of a EPIYA-D or EPIYA-C segment. The East Asian CagA is more closely associated with gastric cancer than the Western CagA. In this study, molecular dynamic (MD) simulations were performed to investigate the binding details of SHP2 and EPIYA segments, and to explore the allosteric regulation mechanism of SHP2. Our results show that the EPIYA-D has a stronger binding affinity to the N-SH2 domain of SHP2 than EPIYA-C. In addition, a single EPIYA-D binding to N-SH2 domain of SHP2 can cause a deflection of the key helix B, and the deflected helix B could squeeze the N-SH2 and PTP domains to break the autoinhibition pocket of SHP2. However, a single EPIYA-C binding to the N-SH2 domain of SHP2 cannot break the autoinhibition of SHP2 because the secondary structure of the key helix B is destroyed. However, the tandem EPIYA-C not only increases its binding affinity to SHP2, but also does not significantly break the secondary structure of the key helix B. Our study can help us better understand the mechanism of gastric cancer caused by Helicobacter pylori infection.

In Silico Analysis Revealed a Unique Binding but Ineffective Mode of Amantadine to Influenza Virus B M2 Channel.

The M2 proton channel of influenza A (AM2) and B (BM2) have a highly conserved function motif, considered as the effective target. As yet, there is no effective drug against BM2. Research showed that AM2 channel blocker, amantadine (AMT), was able to bind to BM2 channel, but AMT lacked inhibition against BM2. Nevertheless, the study of the binding but ineffective mode of AMT to BM2 is challenging. To resolve the challenge and obtain more information for drug design of inhibitors targeting BM2, multiple molecular dynamics simulations were performed. We discovered AMT mainly adopted up binding mode in BM2, involved in a transition flipping from down mode to up mode. Furthermore, we discovered a new key factor to explain ineffective inhibition of AMT to BM2 because of the unmatched spatial geometry between AMT and BM2. Our work could enrich structural feature information on BM2 and provide a new perspective for rational drug design of anti-influenza B.

Exploring the Allosteric Mechanism of Src Homology-2 Domain-Containing Protein Tyrosine Phosphatase 2 (SHP2) by Molecular Dynamics Simulations.

The Src homology-2 (SH2) domain-containing protein tyrosine phosphatase 2 (SHP2, encoded by PTPN11) is a critical allosteric phosphatase for many signaling pathways. Programmed cell death 1 (PD-1) could be phosphorylated at its immunoreceptor tyrosine-based inhibitory motif (ITIM) and immunoreceptor tyrosine-based switch motif (ITSM) and can bind to SHP2 to initiate T cell inactivation. Although the interaction of SHP2-PD-1 plays an important role in the immune process, the complex structure and the allosteric regulation mechanism remain unknown. In this study, molecular dynamics (MD) simulations were performed to study the binding details of SHP2 and PD-1, and explore the allosteric regulation mechanism of SHP2. The results show that ITIM has a preference to bind to the N-SH2 domain and ITSM has almost the same binding affinity to the N-SH2 and C-SH2 domain. Only when ITIM binds to the N-SH2 domain and ITSM binds to the C-SH2 domain can the full activation of SHP2 be obtained. The binding of ITIM and ITSM could change the motion mode of SHP2 and switch it to the activated state.

In Silico Study of Membrane Lipid Composition Regulating Conformation and Hydration of Influenza Virus B M2 Channel.

The proton conduction of transmembrane influenza virus B M2 (BM2) proton channel is possibly mediated by the membrane environment, but the detailed molecular mechanism is challenging to determine. In this work, how membrane lipid composition regulates the conformation and hydration of BM2 channel is elucidated in silico. The appearance of several important hydrogen-bond networks has been discovered, as the addition of negatively charged lipid palmitoyloleoyl phosphatidylglycerol (POPG) and cholesterol reduces membrane fluidity and augments membrane rigidity. A more rigid membrane environment is beneficial to expand the channel, allow more water to enter the channel, promote channel hydration, and then even affect the proton conduction facilitated by the hydrated channel. Thus, membrane environment could be identified as an important influence factor of conformation and hydration of BM2. These findings can provide a unique perspective for understanding the mechanism of membrane lipid composition regulating conformation and hydration of BM2 and have important significance to the further study of anti-influenza virus B 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.

Infection by human immunodeficiency virus type 1 (HIV-1) not only destroys the immune system bringing about acquired immune deficiency syndrome (AIDS), but also induces serious neurological diseases including behavioral abnormalities, motor dysfunction, toxoplasmosis, and HIV-1 associated dementia. The emergence of HIV-1 multidrug-resistant mutants has become a major problem in the therapy of patients with HIV-1 infection. Focusing on the wild type (WT) and G48T/L89M mutated forms of HIV-1 protease (HIV-1 PR) in complex with amprenavir (APV), indinavir (IDV), ritonavir (RTV), and nelfinavir (NFV), we have investigated the conformational dynamics and the resistance mechanism due to the G48T/L89M mutations by conducting a series of molecular dynamics (MD) simulations and free energy (MM-PBSA and solvated interaction energy (SIE)) analyses. The simulation results indicate that alterations in the side-chains of G48T/L89M mutated residues cause the inner active site to increase in volume and induce more curling of the flap tips, which provide the main contributions to weaker binding of inhibitors to the HIV-1 PR. The results of energy analysis reveal that the decrease in van der Waals interactions of inhibitors with the mutated PR relative to the wild-type (WT) PR mostly drives the drug resistance of mutations toward these four inhibitors. The energy decomposition analysis further indicates that the drug resistance of mutations can be mainly attributed to the change in van der Waals and electrostatic energy of some key residues (around Ala28/Ala28' and Ile50/Ile50'). Our work can give significant guidance to design a new generation of anti-AIDS inhibitors targeting PR in the therapy of patients with HIV-1 infection.

Exploring the allosteric mechanism of protein tyrosine phosphatase 1B by molecular dynamics simulations.

Communicated by Ramaswamy H. Sarma.

What are the effects of the serine triad on proton conduction of an influenza B M2 channel? An investigation by molecular dynamics simulations.

The tetrameric influenza B M2 channel (BM2), an acid activated proton channel, is important in the influenza virus B lifecycle. A conserved HxxxW motif is responsible for proton conduction and channel gating. In this study, to explore the effects of the serine triad (S9, S12 and S16) on proton conduction, we performed classical molecular dynamics (CMD) simulations and adaptive steered molecular dynamics (ASMD) simulations at different protonation states of the H19 tetrad. The results of the pore radius and the C-terminal tilt angle show that the electrostatic repulsion induced by protonated H19 is the key driving force for opening the BM2 channel. The open states could be stabilized by the hydrogen bonds between S16 and protonated H19. The solvent accessible surface area and water density indicate that the polar hydrophilic environment provided by the serine triad facilitates the formation of a water wire, and then exhibits favourable effects on proton conduction. The mutant research verifies and supports these views. Our work clarifies the effects of the serine triad on proton conduction in the BM2 channel, which would help us deeply understand the proton conduction mechanism in BM2 and provides a new perspective for antiviral drug design against BM2.

Molecular dynamics simulations study of influence of Tyr422Ala mutation on transcriptional enhancer activation domain 4 (TEAD4) and transcription co-activators complexes.

Transcriptional enhancer activation domain (TEAD) proteins are the downstream transcriptional factor of the Hippo pathway. The transcription co-activators Yes-associated protein (YAP) and its paralog transcription co-activators with PDZ-binding motif (TAZ), binding to TEAD to promote transcription of genes in cell proliferation and anti-apoptosis, are key effectors of the Hippo pathway. TEAD4, one member of TEAD proteins, is specifically required in embryo implantation. The recently reported crystal structure of TEAD4-TAZ complex (PDB Code 5GN0) in mouse reveals that the interactions between the two helices of YAP/TAZ and TEAD4 are highly conserved. Point mutation of the residue Tyr422 of TEAD4 protein would disrupt the relevant hydrogen bond and even abolish the interaction. However, detailed information affected by the mutation at the atom level are still unrevealed. Molecular dynamics (MD) simulations and the molecular mechanics/Generalized-Born surface area (MM/GBSA) free energy calculations were used to explore the effects of mutation Tyr422Ala on the structural flexibility and conformational dynamics. The non-polar interactions play an indispensable role in the binding process of TEAD4 and YAP/TAZ. The helices α1 and α2 of YAP/TAZ provide a primary function to anchor YAP/TAZ well bound to TEAD4. The mutation Tyr422Ala disrupts the hydrogen-bonding network but do not obviously influence the secondary structure stability of TEAD4. The binding conformation of YAP/TAZ distorted by decreased non-polar interaction and the lost hydrogen bonds would lead to reduced interaction activity. The present study would provide important insights into the structure-function relationships of TEAD protein and give a new explanation for the affinity of YAP/TAZ with TEAD.

Molecular dynamics investigation on the Asciminib resistance mechanism of I502L and V468F mutations in BCR-ABL.

Asciminib, a highly selective non-ATP competitive inhibitor of BCR-ABL, has demonstrated to be a promising drug for patients with chronic myeloid leukemia. It is a pity that two resistant mutations (I502L and V468F) have been found during the clinical trial, which is a challenge for the curative effect of Asciminib. In this study, molecular dynamics simulations and molecular mechanics generalized Born surface area (MM-GB/SA) calculations were performed to investigate the molecular mechanism of Asciminib resistance induced by the two mutants. The obtained results indicate that the mutations have adversely influence on the binding of Asciminib to BCR-ABL, as the nonpolar contributions decline in the two mutants. In addition, I502L mutation causes α-helix I' (αI') to shift away from the helical bundle composed of αE, αF, and αH, making the distance between αI' and Asciminib increased. For V468F mutant, the side chain of Phe468 occupies the bottom of the myristoyl pocket (MP), which drives Asciminib to shift toward the outside of MP. Our results provide the molecular insights of Asciminib resistance mechanism in BCR-ABL mutants, which may help the design of novel inhibitors.

MD Simulation Investigation on the Binding Process of Smoke-Derived Germination Stimulants to Its Receptor.

Karrikins (KARs) are a class of smoke-derived seed germination stimulants with great significance in both agriculture and plant biology. By means of direct binding to the receptor protein KAI2, the compounds can initiate the KAR signal transduction pathway, hence triggering germination of the dormant seeds in the soil. In the research, several molecular dynamics (MD) simulation techniques were properly integrated to investigate the binding process of KAR1 to KAI2 and reveal the details of the whole binding event. The calculated binding free energy, -7.00 kcal/mol, is in good agreement with the experimental measurement, -6.83 kcal/mol. The obtained PMF profile indicates the existence of three intermediate states in the binding process. The analysis of the simulation trajectories demonstrates that, in the intermediate structures, KAR1 is stabilized by some hydrophobic residues (Phe26, Phe134, Leu142, Trp153, Phe157, Leu160, Phe194), along with several bridging water molecules, and meanwhile, the significant shifting occurs in the local conformation of the protein as the ligand's binding. A series of the residues (Gln141-Phe157) on the so-called "cap domain" are proposed to be responsible for capturing the ligand at the initial stage of the binding. Besides, the changes of the ligand's poses are also quantitatively characterized by the proper choice of the coordinate system. Our work will contribute to the more penetrating understanding of the ligand binding process and the receptor affinity difference between several members in the KAR family and help design new, more effective germination stimulants.

Insight on mutation-induced resistance to anaplastic lymphoma kinase inhibitor ceritinib from molecular dynamics simulations.

Ceritinib, an advanced anaplastic lymphoma kinase (ALK) next-generation inhibitor, has been proved excellent antitumor activity in the treatment of ALK-associated cancers. However, the accumulation of acquired resistance mutations compromise the therapeutic efficacy of ceritinib. Despite abundant mutagenesis data, the structural determinants for reduced ceritinib binding in mutants remains elusive. Focusing on the G1123S and F1174C mutations, we applied molecular dynamics (MD) simulations to study possible reasons for drug resistance caused by these mutations. The MD simulations predict that the studied mutations allosterically impact the configurations of the ATP-binding pocket. An important hydrophobic cluster is identified that connects P-loop and the αC-helix, which has effects on stabilizing the conformation of ATP-binding pocket. It is suggested, in this study, that the G1123S and F1174C mutations can induce the conformational change of P-loop thereby causing the reduced ceritinib affinity and causing drug resistance.

Exploring the interactional details between aldose reductase (AKR1B1) and 3-Mercapto-5H-1,2,4-triazino[5,6-b]indole-5-acetic acid through molecular dynamics simulations.

Aldose reductase (AKR1B1) has been considered as a significant target for designing drugs to counteract the development of diabetic complications. In the present study, molecular dynamics (MD) simulations and molecular mechanics generalized Born surface area (MM-GB/SA) calculations were performed to make sure which tautomer is the preferred one among three tautomeric forms (Mtia1, Mtia2, and Mtia3) of 3-Mercapto-5H-1,2,4-triazino[5,6-b]indole-5-acetic acid (Mtia) for binding to AKR1B1. The overall structural features and the results of calculated binding free energies indicate that Mtia1 and Mtia2 have more superiority than Mtia3 in terms of binding to AKR1B1. Furtherly, the local active site conformational characteristics and non-covalent interaction analysis were identified. The results indicate that the combination of Mtia2 and AKR1B1 is more stable than that of Mtia1. Furthermore, two extra hydrogen bonds between AKR1B1 and Mtia2 are found with respect to Mtia1. In addition, Mtia2 makes slightly stronger electrostatic interaction with the positively charged nicotinamide group of NADP+ than Mtia1. Based on the results above, Mtia2 is the preferred tautomeric form among the three tautomers. Our study can provide an insight into the details of the interaction between AKR1B1 and Mtia at the atomic level, and will be helpful for the further design of AKR1B1 inhibitors.

Recognition mechanism of Wilms' tumour suppressor protein and DNA triplets: insights from molecular dynamics simulation and free energy analysis.

The Wilms' tumour suppressor protein (WT1) plays a multifaceted role in human cancer processes. Mutations on its DNA recognition domain could lead to Denys-Drash syndrome, and alternate splicing results in insertion of the tripeptide Lys-Thr-Ser (KTS) between the third and fourth zinc fingers (ZFs), leading to changes in the DNA-binding function. However, detailed recognition mechanisms of the WT1-DNA complex have not been explored. To clarify the mutational effects upon WT1 towards DNA binding at the atomic level, molecular dynamics simulations and the molecular mechanics/Poisson Boltzmann surface area (MM/PBSA) method were employed. The simulation results indicate that mutations in ZF domains (E427Q and Q369H) may weaken the binding affinity, and the statistical analyses of the hydrogen bonds and hydrophobic interactions show that eight residues (Lys351, Arg366, Arg375, Arg376, Lys399, Arg403, Arg424 and Arg430) have a significant influence on recognition and binding to DNA. Insertion of the tripeptide KTS could form an immobilized hydrogen-bonding network with Arg403, affecting the flexibility and angle of the linker between ZF3 and ZF4, thus influencing the recognition between the protein and the DNA triplet at its 5' terminus. These results represent the first step towards a thorough characterization of the WT1 recognition mechanisms, providing a better understanding of the structure-function relationship of WT1 and its mutants.

The influence of residue in the position of 116 on the inhibitory potency of TH588 for MTH1.

As one of the first-in-class inhibitor, TH588 was found to be efficient in the suppression of MutT homolog1 (MTH1). A recent work shows that the inhibitory potency of TH588 against human MTH1 (hsMTH1) is approximately 20-fold over that of mouse MTH1 (mmMTH1) and identifies residue in position 116 in MTH1 has an important contribution to TH588 affinity. But the effect of residue Leu or Met in position 116 on the binding affinity remains unclear. In this study, molecular dynamics (MD) simulations and free energy calculations were used to elucidate the mechanism about the effect of residue 116 to the different inhibitory potency of TH588 against MTH1. The binding free energy of TH588 in M116 complexes predicated by the Molecular Mechanics/Generalized Born Surface Area (MM/GBSA) is much lower than that in L116 complexes, which is consistent with the experiment results. The analysis of the individual energy terms suggests that the non-polar interactions are important for distinguishing the binding of TH588. The MD results show that the Leu116 disrupts the interactions between Asn33 and TH588, thus induces the conformational changes of Asn33 as well as TH588. The altered interactions between TH588 and mmMTH1 change the flexibility of TH588, which could induce the remarkable conformational fluctuation of mmMTH1. The conformations of the two loops covering the binding pocket have obvious influence on the opening or closure of the active site. The more open binding site may explain the lower inhibitor potency of TH588 against mmMTH1 than hsMTH1. Our results provide mechanistic insight into the effect of different residue Leu or Met in position 116 on the binding affinity of TH588 for MTH1, which is expected to contribute to the further rational design of more potent inhibitors.

Exploring the inhibition mechanism on HIF-2 by inhibitor PT2399 and 0X3 using molecular dynamics simulations.

Targeting transcription factors HIF-2 is currently considered to be the most direct way for the therapy of clear cell renal cell carcinoma. The preclinical inhibitor PT2399 and artificial inhibitor 0X3 have been identified as promising on-target inhibitors to inhibit the heterodimerization of HIF-2. However, the inhibition mechanism of PT2399 and 0X3 on HIF-2 remains unclear. To this end, molecular dynamics (MD) simulations and molecular docking were applied to investigate the effects of 2 inhibitors on structural motifs and heterodimerization of HIF-2. Our simulation results reveal that the binding of inhibitors disrupts the crucial hydrogen bond and hydrophobic interactions of interdomain of HIF-2 heterodimer due to the local conformational changes of binding interface, confirming the hypothesis that the perturbation of few residues is sufficient to disrupt the heterodimerization of HIF-2. In addition, it can be found that PT2399 with dominant substituents (cyano, fluorine, sulfuryl, and hydroxyl) is more preferred than 0X3 as HIF-2 inhibitor and these substituents play a crucial role in involving more hydrogen bond interactions with residues of interface and then cause the larger structural change of protein. This study may provide a deeper atomic-level insight into the effect of on-target inhibitors on HIF-2 heterodimer, which is expected to contribute to further rational design of effective clear cell renal cell carcinoma drugs.