Molecular Mechanism of Small-Molecule Inhibitors in Blocking the PD-1/PD-L1 Pathway through PD-L1 Dimerization

Programmed cell death-1 (PD-1), which is a molecule involved in the inhibitory signal in the immune system and is important due to blocking of the interactions between PD-1 and programmed cell death ligand-1 (PD-L1), has emerged as a promising immunotherapy for treating cancer. In this work, molecular dynamics simulations were performed on complex systems consisting of the PD-L1 dimer with (S)-BMS-200, (R)-BMS-200 and (MOD)-BMS-200 (i.e., S, R and MOD systems) to systematically evaluate the inhibitory mechanism of BMS-200-related small-molecule inhibitors in detail. Among them, (MOD)-BMS-200 was modified from the original (S)-BMS-200 by replacing the hydroxyl group with a carbonyl to remove its chirality. Binding free energy analysis indicates that BMS-200-related inhibitors can promote the dimerization of PD-L1. Meanwhile, no significant differences were observed between the S and MOD systems, though the R system exhibited a slightly higher energy. Residue energy decomposition, nonbonded interaction, and contact number analyses show that the inhibitors mainly bind with the C, F and G regions of the PD-L1 dimer, while nonpolar interactions of key residues Ile54, Tyr56, Met115, Ala121 and Tyr123 on both PD-L1 monomers are the dominant binding-related stability factors. Furthermore, compared with (S)-BMS-200, (R)-BMS-200 is more likely to form hydrogen bonds with charged residues. Finally, free energy landscape and protein–protein interaction analyses show that the key residues of the PD-L1 dimer undergo remarkable conformational changes induced by (S)-BMS-200, which boosts its intimate interactions. This systematic investigation provides a comprehensive molecular insight into the ligand recognition process, which will benefit the design of new small-molecule inhibitors targeting PD-L1 for use in anticancer therapy.     

聚合物阻垢剂阻垢机理的分子动力学研究

用分子模拟方法构建了碳酸钙六方体晶胞结构,优化了丙烯酸与丙烯酸甲酯共聚物的构型,模拟了单体不同配比时各共聚物与碳酸钙晶体之间的相互作用,并计算了其作用能量的变化。计算结果显示:所有阻垢剂分子都逐渐接近方解石晶体,活性基团占据晶体的晶格生长点或嵌入晶体内部,同时伴随小幅度的分子结构变形,抑制晶体成核而有效阻止垢的形成,或使晶体内部形成空洞,导致晶格畸变而使垢松软。同时随着丙烯酸单体的比例增大,阻垢作用更明显。     

Evaluation of the Molecular Landscape in PD-L1 Positive Metastatic NSCLC: Data from Campania,Italy

Background: Immune-checkpoint inhibitors (ICIs) have increased and improved the treatment options for patients with non-oncogene-addicted advanced stage non-small cell lung cancer (NSCLC). However, the role of ICIs in oncogene-addicted advanced stage NSCLC patients is still debated. In this study, in an attempt to fill in the informational gap on the effect of ICIs on other driver mutations, we set out to provide a molecular landscape of clinically relevant oncogenic drivers in programmed death-ligand 1 (PD-L1) positive NSCLC patients. Methods: We retrospectively reviewed data on 167 advanced stage NSCLC PD-L1 positive patients (≥1%) who were referred to our clinic for molecular evaluation of five driver oncogenes, namely, EGFR, KRAS, BRAF, ALK and ROS1. Results: Interestingly, n = 93 (55.7%) patients showed at least one genomic alteration within the tested genes. Furthermore, analyzing a subset of patients with PD-L1 tumor proportion score (TPS) ≥ 50% and concomitant gene alterations (n = 8), we found that n = 3 (37.5%) of these patients feature clinical benefit with ICIs administration, despite the presence of a concomitant KRAS gene alteration. Conclusions: In this study, we provide a molecular landscape of clinically relevant biomarkers in NSCLC PD-L1 positive patients, along with data evidencing the clinical benefit of ICIs in patient NSCLC PD-L1 positive alterations.     

Synthesis of Dihydrobenzofuro[3,2-b]chromenes as Potential 3CLpro Inhibitors of SARS-CoV-2: A Molecular Docking and Molecular Dynamics Study

The recent emergence of pandemic of coronavirus (COVID-19) caused by SARS-CoV-2 has raised significant global health concerns. More importantly, there is no specific therapeutics currently available to combat against this deadly infection. The enzyme 3-chymotrypsin-like cysteine protease (3CLpro) is known to be essential for viral life cycle as it controls the coronavirus replication. 3CLpro could be a potential drug target as established before in the case of severe acute respiratory syndrome coronavirus (SARS-CoV) and Middle East respiratory syndrome coronavirus (MERS-CoV). In the current study, we wanted to explore the potential of fused flavonoids as 3CLpro inhibitors. Fused flavonoids (5a,10a-dihydro-11H-benzofuro[3,2-b]chromene) are unexplored for their potential bioactivities due to their low natural occurrences. Their synthetic congeners are also rare due to unavailability of general synthetic methodology. Here we designed a simple strategy to synthesize 5a,10a-dihydro-11H-benzofuro[3,2-b]chromene skeleton and it's four novel derivatives. Our structural bioinformatics study clearly shows excellent potential of the synthesized compounds in comparison to experimentally validated inhibitor N3. Moreover, in-silico ADMET study displays excellent druggability and extremely low level of toxicity of the synthesized molecules. Further, for better understanding, the molecular dynamic approach was implemented to study the change in dynamicity after the compounds bind to the protein. A detailed investigation through clustering analysis and distance calculation gave us sound comprehensive data about their molecular interaction. In summary, we anticipate that the currently synthesized molecules could not only be a potential set of inhibitors against 3CLpro but also the insights acquired from the current study would be instrumental in further developing novel natural flavonoid based anti-COVID therapeutic spectrums.     

A Combination of 3D-QSAR, Molecular Docking and Molecular Dynamics Simulation Studies of Benzimidazole-Quinolinone Derivatives as iNOS Inhibitors

Inducible Nitric Oxide Synthase (iNOS) has been involved in a variety of diseases, and thus it is interesting to discover and optimize new iNOS inhibitors. In previous studies, a series of benzimidazole-quinolinone derivatives with high inhibitory activity against human iNOS were discovered. In this work, three-dimensional quantitative structure-activity relationships (3D-QSAR), molecular docking and molecular dynamics (MD) simulation approaches were applied to investigate the functionalities of active molecular interaction between these active ligands and iNOS. A QSAR model with R² of 0.9356, Q² of 0.8373 and Pearson-R value of 0.9406 was constructed, which presents a good predictive ability in both internal and external validation. Furthermore, a combined analysis incorporating the obtained model and the MD results indicates: (1) compounds with the proper-size hydrophobic substituents at position 3 in ring-C (R₃ substituent), hydrophilic substituents near the X₆ of ring-D and hydrophilic or H-bond acceptor groups at position 2 in ring-B show enhanced biological activities; (2) Met368, Trp366, Gly365, Tyr367, Phe363, Pro344, Gln257, Val346, Asn364, Met349, Thr370, Glu371 and Tyr485 are key amino acids in the active pocket, and activities of iNOS inhibitors are consistent with their capability to alter the position of these important residues, especially Glu371 and Thr370. The results provide a set of useful guidelines for the rational design of novel iNOS inhibitors.     

Combining Molecular Docking and Molecular Dynamics to Predict the Binding Modes of Flavonoid Derivatives with the Neuraminidase of the 2009 H1N1 Influenza A Virus

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

Molecular Docking and Molecular Dynamics Simulation Studies of Triterpenes from Vernonia patula with the Cannabinoid Type 1 Receptor

A molecular docking approach was employed to evaluate the binding affinity of six triterpenes, namely epifriedelanol, friedelin, α-amyrin, α-amyrin acetate, β-amyrin acetate, and bauerenyl acetate, towards the cannabinoid type 1 receptor (CB1). Molecular docking studies showed that friedelin, α-amyrin, and epifriedelanol had the strongest binding affinity towards CB1. Molecular dynamics simulation studies revealed that friedelin and α-amyrin engaged in stable non-bonding interactions by binding to a pocket close to the active site on the surface of the CB1 target protein. The studied triterpenes showed a good capacity to penetrate the blood–brain barrier. These results help to provide some evidence to justify, at least in part, the previously reported antinociceptive and sedative properties of Vernonia patula.     

Focus on PD-1/PD-L1 as a Therapeutic Target in Ovarian Cancer

Ovarian cancer is considered one of the most aggressive and deadliest gynecological malignancies worldwide. Unfortunately, the therapeutic methods that are considered the gold standard at this moment are associated with frequent recurrences. Survival in ovarian cancer is associated with the presence of a high number of intra tumor infiltrating lymphocytes (TILs). Therefore, immunomodulation is considered to have an important role in cancer treatment, and immune checkpoint inhibitors may be useful for restoring T cell-mediated antitumor immunity. However, the data presented in the literature until now are not sufficient to allow for the identification and selection of patients who really respond to immunotherapy among those with ovarian cancer. Although there are some studies with favorable results, more prospective trials are needed in this sense. This review focuses on the current and future perspectives of PD-1/L1 blockade in ovarian cancer and analyzes the most important immune checkpoint inhibitors used, with the aim of achieving optimal clinical outcomes. Future studies and trials are needed to maximize the efficacy of immune checkpoint blockade therapy in ovarian cancer, as well as in all cancers, in general.     

气体在水中的分子动力学模拟

采用分子动力学（ＭＤ）模拟的方法在常温及工业应用背景条件下对ＣＨ４、ＮＨ３、ＣＯ２、Ｏ２这些气体在水中的结构及扩散情形进行了研究。ＭＤ模拟可以为这些涉及到气体在水中的工业应用情形的机理提供分子水平的解释,同时ＭＤ模拟还可为一些不易实验测定扩散性质的体系提供工程初步设计和过程开发所需的数据。     

Identification of the FDA-Approved Drug Pyrvinium as a Small-Molecule Inhibitor of the PD-1/PD-L1 Interaction

Immune checkpoint blockade involving inhibition of the PD-1/PD-L1 interaction has provided unprecedented clinical benefits in treating a variety of tumors. To date, a total of six antibodies that bind to either PD-1 or PD-L1 protein and in turn inhibit the PD-1/PD-L1 interaction have received clinical approvals. Despite being highly effective, these expensive large biotherapeutics possess several inherent pharmacokinetic limitations that can be successfully overcome through the use of low-molecular-weight inhibitors. One such promising approach involves small-molecule induced dimerization and sequestration of PD-L1, leading to effective PD-1/PD-L1 inhibition. Herein, we present the discovery of such potential bioactive PD-L1 dimerizers through a structure- and ligand-based screening of a focused library of approved and investigational drugs worldwide. Pyrvinium, an FDA-approved anthelmintic drug, showed the highest activity in our study with IC₅₀ value of ∼29.66 μM. It is noteworthy that Pyrvinium, being an approved drug, may prove especially suitable as a good starting point for further medicinal chemistry efforts, leading to design and development of even more potent structural analogs as selective PD-1/PD-L1 inhibitors. Furthermore, the adopted integrated virtual screening protocol may prove useful in screening other larger databases of lead- and drug-like molecules for hit identification in the domain of small-molecule PD-1/PD-L1 inhibitors.     

Interactions of Borneol with DPPC Phospholipid Membranes: A Molecular Dynamics Simulation Study

Borneol, known as a “guide” drug in traditional Chinese medicine, is widely used as a natural penetration enhancer in modern clinical applications. Despite a large number of experimental studies on borneol’s penetration enhancing effect, the molecular basis of its action on bio-membranes is still unclear. We carried out a series of coarse-grained molecular dynamics simulations with the borneol concentration ranging from 3.31% to 54.59% (v/v, lipid-free basis) to study the interactions of borneol with aDPPC(1,2-dipalmitoylsn-glycero-3-phosphatidylcholine) bilayer membrane, and the temperature effects were also considered. At concentrations below 21.89%, borneol’s presence only caused DPPC bilayer thinning and an increase in fluidity; A rise in temperature could promote the diffusing progress of borneol. When the concentration was 21.89% or above, inverted micelle-like structures were formed within the bilayer interior, which led to increased bilayer thickness, and an optimum temperature was found for the interaction of borneol with the DPPC bilayer membrane. These findings revealed that the choice of optimal concentration and temperature is critical for a given application in which borneol is used as a penetration enhancer. Our results not only clarify some molecular basis for borneol’s penetration enhancing effects, but also provide some guidance for the development and applications of new preparations containing borneol.     

KEAP1 Cancer Mutants: A Large-Scale Molecular Dynamics Study of Protein Stability

We have performed 280 μs of unbiased molecular dynamics (MD) simulations to investigate the effects of 12 different cancer mutations on Kelch-like ECH-associated protein 1 (KEAP1) (G333C, G350S, G364C, G379D, R413L, R415G, A427V, G430C, R470C, R470H, R470S and G476R), one of the frequently mutated proteins in lung cancer. The aim was to provide structural insight into the effects of these mutants, including a new class of ANCHOR (additionally NRF2-complexed hypomorph) mutant variants. Our work provides additional insight into the structural dynamics of mutants that could not be analyzed experimentally, painting a more complete picture of their mutagenic effects. Notably, blade-wise analysis of the Kelch domain points to stability as a possible target of cancer in KEAP1. Interestingly, structural analysis of the R470C ANCHOR mutant, the most prevalent missense mutation in KEAP1, revealed no significant change in structural stability or NRF2 binding site dynamics, possibly indicating an covalent modification as this mutant’s mode of action.     

Molecular Dynamics Simulation of Tryptophan Hydroxylase-1: Binding Modes and Free Energy Analysis to Phenylalanine Derivative Inhibitors

Serotonin is a neurotransmitter that modulates many central and peripheral functions. Tryptophan hydroxylase-1 (TPH1) is a key enzyme of serotonin synthesis. In the current study, the interaction mechanism of phenylalanine derivative TPH1 inhibitors was investigated using molecular dynamics (MD) simulations, free energy calculations, free energy decomposition analysis and computational alanine scanning. The predicted binding free energies of these complexes are consistent with the experimental data. The analysis of the individual energy terms indicates that although the van der Waals and electrostatics interaction contributions are important in distinguishing the binding affinities of these inhibitors, the electrostatic contribution plays a more crucial role in that. Moreover, it is observed that different configurations of the naphthalene substituent could form different binding patterns with protein, yet lead to similar inhibitory potency. The combination of different molecular modeling techniques is an efficient way to interpret the interaction mechanism of inhibitors and our work could provide valuable information for the TPH1 inhibitor design in the future.     

Reconstitution of Monocyte Subsets and PD-L1 Expression but Not T Cell PD-1 Expression in Obstructive Sleep Apnea Patients upon PAP Therapy

Obstructive sleep apnea (OSA) is characterized by nocturnal breathing intermissions resulting in oxidative stress and eventually, a low-grade systemic inflammation. The study aimed to investigate the impact of positive airway pressure (PAP) therapy on the inflammatory milieu as measured by monocyte and T cell phenotypic alterations. Participants were assessed for their OSA severity before PAP therapy and about six months later, including patient-reported outcome and therapy usage by telemetry readout. The distributions of the CD14/CD16-characterized monocyte subsets as well as the CD4/CD8-characterized effector T cell subsets with regard to their PD-1 and PD-L1 expression were analyzed by flow cytometry from blood samples. Data of 25 patients revealed a significant reconstitution of the monocyte subset distribution and a decrease in PD-L1 expression on pan-monocytes and CD8⁺ T cells without an association to initial AHI and overweight. The PD-1 expression was still increased on T cell subsets, especially on CD4⁺ TH17/22 cells. We conclude that PAP therapy might have a rapid effect on the monocyte phenotype and overall PD-L1 expression levels. However, T cell immune alterations especially on TH17/22 cells persist longer, indicating an ongoing disturbance of the adaptive immune system.     

NPT系综分子动力学模拟水的密度和扩散系数

采用常温常压(NPT)系综分子动力学方法模拟了298.15K、0.1013MPa条件下水的密度和自扩散系数,模拟结果与实验值一致,相对误差分别为1.74%和3.83%。     

DNTF对NC塑化特性的分子动力学模拟及实验研究

为研究DNTF对CMDB推进剂力学性能的作用机理,通过分子动力学模拟方法建立了NC纯物质和NC/DNTF共混物分子模型,研究了DNTF对NC塑化过程中微观结构的影响,采用拉伸试验研究了DNTF对CMDB推进剂力学性能的影响。结果表明,DNTF与NC分子形成氢键,使NC内分子的氢键减弱,NC分子的刚性降低;加入DNTF后,DNTF-CMDB推进剂的抗拉强度降低、延伸率增加;20℃时延伸率由8.69%增加到33.6%,50℃时延伸率由14.86%增加到45.6%。分子动力学模拟计算结果与拉伸试验结果一致。     

氨的自扩散系数的分子动力学模拟研究

采用分子动力学模拟方法研究了氨在较宽温度和压力范围的分子自扩散系数。从常温到高温,自扩散系数的模拟值与实验值吻合得很好,这表明可以采用分子动力学模拟来代替实验,获得高温高压条件下实验难以测量的自扩散系数。     

甲烷的自扩散系数的分子动力学模拟研究

运用分子动力学模拟方法研究了甲烷在较宽温度和压力范围的分子自扩散系数。结果表明,甲烷自扩散系数的模拟值与文献实验值吻合得较好。因此,可以采用分子动力学模拟来代替实验,获得高温高压条件下实验难以测量的甲烷自扩散系数。     

Interaction of Classical Platinum Agents with the Monomeric and Dimeric Atox1 Proteins: A Molecular Dynamics Simulation Study

We carried out molecular dynamics simulations and free energy calculations for a series of binary and ternary models of the cisplatin, transplatin and oxaliplatin agents binding to a monomeric Atox1 protein and a dimeric Atox1 protein to investigate their interaction mechanisms. All three platinum agents could respectively combine with the monomeric Atox1 protein and the dimeric Atox1 protein to form a stable binary and ternary complex due to the covalent interaction of the platinum center with the Atox1 protein. The results suggested that the extra interaction from the oxaliplatin ligand–Atox1 protein interface increases its affinity only for the OxaliPt + Atox1 model. The binding of the oxaliplatin agent to the Atox1 protein might cause larger deformation of the protein than those of the cisplatin and transplatin agents due to the larger size of the oxaliplatin ligand. However, the extra interactions to facilitate the stabilities of the ternary CisPt + 2Atox1 and OxaliPt + 2Atox1 models come from the α1 helices and α2-β4 loops of the Atox1 protein–Atox1 protein interface due to the cis conformation of the platinum agents. The combinations of two Atox1 proteins in an asymmetric way in the three ternary models were analyzed. These investigations might provide detailed information for understanding the interaction mechanism of the platinum agents binding to the Atox1 protein in the cytoplasm.