Per-residue Energy Decomposition Research Articles

Bioinformatics-driven discovery of novel EGFR kinase inhibitors as anti-cancer therapeutics: In silico screening and in vitro evaluation.

Epidermal growth factor receptor EGFR inhibitors are widely used as first line therapy for the treatment of non-small-cell lung cancer (NSCLC) in patients harboring EGFR mutation. However, the acquisition of a second-site mutation (T790 M) limited the efficacy and developed resistance. Therefore, discovery and development of specific drug target for this mutation is of urgent needs. In our study we used the ChemDiv diversity database for receptor-based virtual screening to secure EGFR-TK inhibitors chemotherapeutics. We identified four compounds that bind to the ATP-binding region of the EGFR-TK using AutoDock 4.0 and AutoDock Vina1.1.2 and post-docking investigations. The ligand showed hydrophobic interactions to the hydrophobic region of the binding site and engaged in hydrogen bonding with Met793. The ligands also explored π-cation interactions between the π-system of the ligand-phenyl ring and the positive amino group of Lys745. Molecular mechanics Poisson-Boltzmann surface area MM/PBSA per-residue energy decomposition analyses revealed that Val726, Leu792, Met793, Gly796, Cys797, Leu798, and Thr844 contributed the most to the binding energy. Biological evaluation of the retrieved hit compounds showed suppressing activity against EGFR auto phosphorylation and selective apoptosis-induced effects toward lung cancer cells harboring the EGFR L858R/T790M double mutation. Our work anticipated into novel and specific EGFR-TKIs and identified new compounds with therapeutic potential against lung cancer.

From Biological Activity to Stereoselectivity: A Portrait of Molecular and Mechanistic Profiles of the Therapeutic Potential of G-1 and LNS8801 as GPER-1 Activator in the Treatment of Waldenström’s Macroglobulinemia

Therapeutic activation of G protein-coupled estrogen receptor 1 (GPER-1) shows promise for treating Waldenström's Macroglobulinemia (WM). Recently, the selective small-molecule agonist G-1 or its enantiomer LNS8801 were reported as potent GPER-1 pharmacological activators, however, the molecular events surrounding their chirality towards GPER-1 is still unexplored. This study aimed to explore the molecular events surrounding the chirality of the selective small-molecule agonists G-1 and LNS8801 towards GPER-1. Molecular docking and dynamics simulations revealed strong binding affinities and key amino acid residues involved in the activation process. The results revealed favourable binding affinities of -29.87kcal/mol, and -28.09kcal/mol for G-1 and LNS8801 towards GPER-1, respectively. Per-residue energy decomposition and time-dependent analysis proved that Arg253 and Arg254 are key amino acid residues in binding these activators towards GPER-1. The activators affected stability, flexibility, and structure of GPER-1. G-1 exhibited greater rigidity when bound to GPER-1 compared to LNS880. This study's findings illuminate chirality and the potential for optimizing the enantiomer ratio to enhance inhibitory effects. These results could also reveal the molecular specifics and binding mechanism of the G-1 and its enantiomer LNS8801 with GPER-1, laying the basis for the rational design of enhanced chiral molecules against WM.

Investigation of the Binding Interaction of Mfsd2a with NEDD4-2 via Molecular Dynamics Simulations.

Major facilitator superfamily domain-containing 2a (Mfsd2a) is a sodium-dependent lysophosphatidylcholine cotransporter that plays an important role in maintaining the integrity of the blood-brain barrier and neurological function. Abnormal degradation of Mfsd2a often leads to dysfunction of the blood-brain barrier, while upregulation of Mfsd2a can retrieve neurological damage. It has been reported that Mfsd2a can be specifically recognized and ubiquitinated by neural precursor cell-expressed developmentally downregulated gene 4 type 2 (NEDD4-2) ubiquitin ligase and finally degraded through the proteasome pathway. However, the structural basis for the specific binding of Mfsd2a to NEDD4-2 is unclear. In this work, we combined deep learning and molecular dynamics simulations to obtain a Mfsd2a structure with high quality and a stable Mfsd2a/NEDD4-2-WW3 interaction model. Moreover, molecular mechanics generalized Born surface area (MM-GBSA) methods coupled with per-residue energy decomposition studies were carried out to analyze the key residues that dominate the binding interaction. Based on these results, we designed three peptides containing the key residues by truncating the Mfsd2a sequences. One of them was found to significantly inhibit Mfsd2a ubiquitination, which was further validated in an oxygen-glucose deprivation (OGD) model in a human microvascular endothelial cell line. This work provides some new insights into the understanding of Mfsd2a and NEDD4-2 interaction and might promote further development of drugs targeting Mfsd2a ubiquitination.

In silico analysis, isolation, and cytotoxicity evaluation of the coumestans from Cullen corylifolium (L.) Medik

Cullen corylifolium is well known for diverse phytoconstituents that possess multifaceted pharmacology, and one such less explored class is coumestans, which have not been well explored for their anticancer activities. One of the popular cancer targets is the Epidermal Growth Factor Receptor, a tyrosine kinase involved in various cancers, especially breast and lung cancer hence, a crucial cancer target. This work is focussed on molecular docking and molecular simulation studies on coumestans against EGFR. The rigorous docking studies resulted in two coumestans (1 and 5) with binding energy less than Gefitinib and Erlotinib. Compounds 1 and 5 were subjected to molecular simulation, binding free energy calculation, per-residue energy decomposition, and in silico ADMET prediction. The best hit, compound 1 was evaluated for its cytotoxicity against MDA-MB-231 and A549 cells via in vitro assay. The ligand-protein complex exhibited good stability, binding free energies, better in silico pharmacokinetics, low toxicity, and good cytotoxicity.

Investigating the ibrutinib resistance mechanism of L528W mutation on Brutonʼs tyrosine kinase via molecular dynamics simulations

Drug resistance to Bruton's Tyrosine Kinase (BTK) inhibitors presents a challenge in treating B-cell malignancies, and the mechanism behind drug resistance remains unclear. In this study, we focused on the BTK L528W mutation and investigated the underlying mechanisms of resistance to ibrutinib (including prototype and its active metabolite from, PCI-45227) using a combination of bioinformatics analysis, and molecular dynamics (MD) simulations. Protein stability of wild type (WT) BTK and L528W mutant was predicted using DUET, PoPMuSiC, and I-Mutant2.0. We performed MD simulations of six systems, apo-WT, metabolite-WT, prototype-WT and their mutants, to analyze the significant conformational and BTK-inhibitor binding affinity changes induced by the L528W mutation. Results show that the L528W mutation reduces the conformational stability of BTK compared to the WT. Principal component analysis (PCA) based free energy landscape (FEL) analysis shows that the L528W mutant ensemble tends to form more conformation clusters and exhibit higher levels of local minima than the WT counterpart. The interaction analysis reveal that the L528W mutation disrupts the strong hydrogen bond between Cys481 and inhibitors and reduces the number of hydrogen bonds between inhibitors and BTK in the L528W mutant complex structures compared to the WT. Porcupine plot analysis in association with cross-correlation analysis show the high-intensity flexible motion exhibited by the P-loop region. MM/GBSA calculations show that the L528W mutation in metabolite-BTK and prototype-BTK complexes increases binding free energy compared to the WT, with a reduction in binding affinity confirmed by per-residue energy decomposition. Specifically, the binding free energy increases from −57.86 kcal/mol to −48.26 kcal/mol for the metabolite-BTK complex and from −62.04 kcal/mol to −50.55 kcal/mol for the prototype-BTK complex. Overall, our study finds that the L528W mutation reduces BTK stability, decreases binding affinity, and leads to drug resistance and potential disease recurrence.

Elucidating Residue-Level Determinants Affecting Dimerization of Ebola Virus Matrix Protein Using High-Throughput Site Saturation Mutagenesis and Biophysical Approaches.

The Ebola virus (EBOV) is a filamentous virus that acquires its lipid envelope from the plasma membrane of the host cell it infects. EBOV assembly and budding from the host cell plasma membrane are mediated by a peripheral protein, known as the matrix protein VP40. VP40 is a 326 amino acid protein with two domains that are loosely linked. The VP40 N-terminal domain (NTD) contains a hydrophobic α-helix, which mediates VP40 dimerization. The VP40 C-terminal domain has a cationic patch, which mediates interactions with anionic lipids and a hydrophobic region that mediates VP40 dimer-dimer interactions. The VP40 dimer is necessary for trafficking to the plasma membrane inner leaflet and interactions with anionic lipids to mediate the VP40 assembly and oligomerization. Despite significant structural information available on the VP40 dimer structure, little is known on how the VP40 dimer is stabilized and how residues outside the NTD hydrophobic portion of the α-helical dimer interface contribute to dimer stability. To better understand how VP40 dimer stability is maintained, we performed computational studies using per-residue energy decomposition and site saturation mutagenesis. These studies revealed a number of novel keystone residues for VP40 dimer stability just adjacent to the α-helical dimer interface as well as distant residues in the VP40 CTD that can stabilize the VP40 dimer form. Experimental studies with representative VP40 mutants in vitro and in cells were performed to test computational predictions that reveal residues that alter VP40 dimer stability. Taken together, these studies provide important biophysical insights into VP40 dimerization and may be useful in strategies to weaken or alter the VP40 dimer structure as a means of inhibiting the EBOV assembly.

Investigation into the Interaction Sites of the K84s and K102s Peptides with α-Synuclein for Understanding the Anti-Aggregation Mechanism: An In silico Study

Background: α-Synuclein has become the main therapeutic target in Parkinson's disease and related Synucleinopathies since the discovery of genetic associations between α-Synuclein and Parkinson's disease risk and the identification of aggregated α-Synuclein as the primary protein con-stituent of Lewy pathology two decades ago. The two new peptides K84s (FLVWGCLRGSAI-GECVVHGGPPSRH) and K102s (FLKRWARSTRWGTASCGGS) have recently been found to significantly reduce the oligomerization and aggregation of α-Synuclein. However, it is still unclear where these peptides interact with α-Synuclein at the moment. Objective: To examine the locations where K84s and K102s interact with α-Synuclein. Methods: In this investigation, the PEPFOLD3 server was used to generate the 3-D structures of the K84s and K102s peptides. Using the PatchDock web server, the two peptides were docked to the α-Synuclein molecule. After that, 50 ns of Molecular Dynamics (MD) simulations using the Am-berff99SBildn force field were performed on the two resulting docked complexes. The two com-plexes' structure, dynamics, energy profiles, and binding modes were identified through analysis of the respective MD simulation trajectories. By submitting the two complexes' lowest energy struc-ture to the PDBsum website, the interface residues in the two complexes were identified. The per residue energy decomposition (PRED) analysis using the MM-GBSA technique was used to calcu-late the contributions of each residue in the α-Synuclein of (α-Synuclein-K84s/K102s) complexes to the total binding free energy. Results: The binding of the two peptides with the α-Synuclein was demonstrated to have high bind-ing free energy. The binding free energies of the (α-Synuclein-K84s) and (α-Synuclein-K102s) complexes are -33.61 kcal/mol and -40.88 kcal/mol respectively. Using PDBsum server analysis, it was determined that in the (α-Synuclein-K84s) complex, the residues GLY 25, ALA 29, VAL 49, LEU 38, VAL 40, GLU 28, GLY 47, LYS 32, GLU 35, GLY 36, TYR 39, VAL 48 and VAL 26 (from α-Synuclein) and SER 23, LEU 7, ILE 12, HIS 25, PHE 1, HIS 18, CYS 6, ARG 24, PRO 21 and ARG 8 (from K84s peptide) were identified to be present at the interface. In the (α-Synuclein-K102s) complex, the residues VAL 40, GLY 36, GLU 35, TYR 39, LYS 45, LEU 38, LYS 43, VAL 37, THR 44, VAL 49, VAL 48, and GLU 46 (from α-Synuclein) and ARG 10, GLY 12, GLY 18, SER 15, THR 13, SER 19, TRP 11, ALA 14, CYS 16, ARG 7, ARG 4 and GLY 17 (from K102s peptide) were identified to be present at the interface. The PRED analysis revealed that the residues PHE 1, LEU 7, ILE 12, LEU 2, VAL 3, GLY 5, and PRO 21 of the K84s peptide and resi-dues VAL 48, ALA 29, VAL 40, TYR 39, VAL 49, VAL 26 and GLY 36 of α-Synuclein in the (α-Synuclein-K84s) complex are responsible for the intermolecular interaction. The residues ARG 4, ARG 10, TRP 11, ALA 14, SER 15, CYS 16 and SER 19 of the K102s peptide and residues GLU 46, LYS 45, VAL 49, GLU 35, VAL 48, TYR 39, and VAL 40 of α-Synuclein are responsible for the intermolecular interaction in the instance of the (α-Synuclein-K102s) complex. Additionally, it has been found that a sizable portion of the helical structure is preserved when α-Synuclein is in a complex form with the K84s and K102s peptides. Conclusion: Taken together the data implies that the two new peptides investigated here could be suitable candidates for future therapeutic development against α-Synuclein aggregation.

3D-QSAR, docking and molecular dynamics simulations of novel Pyrazolo-pyridazinone derivatives as covalent inhibitors of FGFR1: a scientific approach for possible anticancer agents

Developing highly potent covalent inhibitors of Fibroblast growth factor receptors 1 (FGFR1) has always been a challenging task. In the current study, various computational techniques, such as 3D-QSAR, covalent docking, fingerprinting analysis, MD simulation followed by MMGB/PBSA, and per-residue energy decomposition analysis were used to explore the binding mechanism of pyrazolo[3,4-d]pyridazinone derivatives to FGFR1. The high q2 and r2 values for the CoMFA and CoMSIA models, suggest that the constructed 3D-QSAR models could reliably predict the bioactivities of FGFR1 inhibitors. The structural requirements revealed by the model’s contour maps were strategically used to computationally create an in-house library of more than 100 new FGFR1 inhibitors using the R-group exploration technique implemented in the SparkTM software. The compounds from the in-house library were also mapped in the 3D-QSAR model that predicts comparable pIC50 values with the experimental values. A comparison between 3D-QSAR generated contours and molecular docking conformation of ligands was performed to reveal the fundamentals to design potent FGFR1 covalent inhibitors. The estimated binding free energies (MMGB/PBSA) for the selected compounds were in agreement with the experimental value ranking of their binding affinities towards FGFR1. Furthermore, per-residue energy decomposition analysis has identified Arg627 and Glu531 to contribute significantly in improved binding affinity of compound W16. During ADME analysis, the majority of in-house library compounds exhibited pharmacokinetic properties superior to those of experimentally produced compounds. These new compounds may help researchers better understand FGFR1 inhibition and lead to the creation of novel, potent FGFR1 inhibitors. Communicated by Ramaswamy H. Sarma

Molecular modeling of pyrrolo-pyrimidine based analogs as potential FGFR1 inhibitors: a scientific approach for therapeutic drugs

Fibroblast growth factor receptors 1 (FGFR1) is an emerging target for the development of anticancer drugs. Uncontrolled expression of FGFR1 is strongly associated with a number of different types of cancers. Apart from a few FGFR inhibitors, the FGFR family members have not been thoroughly studied to produce clinically effective anticancer drugs. The application of proper computational techniques may aid in understanding the mechanism of protein–ligand complex formation, which may provide a better notion for developing potent FGFR1 inhibitors. In this study, a variety of computational techniques, including 3D-QSAR, flexible docking and MD simulation followed by MMGB/PBSA, H-bonds and distance analysis, have been performed to systematically explore the binding mechanism of pyrrolo-pyrimidine derivatives against FGFR1. The 3D-QSAR model was generated to deduce the structural determinants of FGFR1 inhibition. The high q2 and r2 values for the CoMFA and CoMSIA models indicated that the created 3D-QSAR models could reliably predict the bioactivities of FGFR1 inhibitors. The computed binding free energies (MMGB/PBSA) for the selected compounds were consistent with the ranking of their experimental binding affinities against FGFR1. Furthermore, per-residue energy decomposition analysis revealed that the residues Lys514 in catalytic region, Asn568, Glu571 in solvent accessible portion and Asp641 in DFG motif exhibited a strong tendency to mediate ligand–protein interactions through the hydrogen bonding and Van Der Waals interactions. These findings may benefit researchers in gaining better knowledge of FGFR1 inhibition and may serve as a guideline for the development of novel and highly effective FGFR1 inhibitors.Communicated by Ramaswamy H. Sarma

Insights Into Resveratrol as an Inhibitor Against Aβ1-42 Peptide Aggregation: A Molecular Dynamics Simulation Study

Background: Resveratrol (RSV), a polyphenolic compound, is reported to have anti-aggregation properties against Amyloid-beta peptides. It is, therefore, significant to understand the mechanism of inhibition of Aβ1-42 peptide aggregation by the RSV at the molecular level. We have used Molecular docking along with Molecular dynamics (MD) simulation techniques to address the role of RSV in the inhibition of Aβ1-42 pep-tide aggregation. Objective: To understand the role of Resveratrol on the Aβ1-42 peptide aggregation. Method: In this computational study, we have docked the RSV to Aβ1-42 peptide using Molecular Docking software and then performed MD simulation for the Aβ1-42 peptide monomer Aβ1-42 peptide-RSV complex using the AMBER force field. From the analysis of MD trajectories, we obtained salient structural features and determined the Binding Free Energy(BFE) and Per-residue Energy Decomposition Analysis (PRED) using MM-PBSA/GBSA method. Results: The secondary structure and the conformational analysis obtained from MD trajectories show that the binding of RSV with the Aβ1-42 peptide monomer causes an increase in the helical content in the structure of the Aβ1-42 peptide. The BFE and PRED results show a high binding affinity (GBtotal=-11.07 kcal mol-1; PBtotal= -1.82 kcal mol-1) of RSV with Aβ1-42 peptide. Also, we found the RSV to interact with cru-cial residues (Asp 23 and Lys 28) of the Aβ1-42 peptide. These residues play a signif-icant role in facilitating the formation of toxic amyloid oligomers and amyloid fibrils. The salt bridge interaction between these residues D23–K28 was found to be destabi-lized in the Aβ1-42 peptide when it is complexed with RSV. Conclusion: In summary, it can be concluded that Resveratrol greatly aids the preven-tion of Aβ1-42 peptide aggregation. Therefore, it can be considered a possible drug candidate for therapeutic strategies for Alzheimer’s disease.

Bioinformatic analysis of highly consumed phytochemicals as P-gp binders to overcome drug-resistance.

P-glycoprotein (P-gp) is an adenosine triphosphate (ATP)-dependent membrane efflux pump for protecting cells against xenobiotic compounds. Unfortunately, overexpressed P-gp in neoplastic cells prevents cell entry of numerous chemotherapeutic agents leading to multidrug resistance (MDR). MDR cells may be re-sensitized to chemotherapeutic drugs via P-gp inhibition/modulation. Side effects of synthetic P-gp inhibitors encouraged the development of natural products. Molecular docking and density functional theory (DFT) calculations were used as fast and accurate computational methods to explore a structure binding relationship of some dietary phytochemicals inside distinctive P-gp binding sites (modulatory/inhibitory). For this purpose, top-scored docked conformations were subjected to per-residue energy decomposition analysis in the B3LYP level of theory with a 6-31g (d, p) basis set by Gaussian98 package. Consecutive application of computational techniques revealed binding modes/affinities of nutritive phytochemicals within dominant binding sites of P-gp. Blind docking scores for best-ranked compounds were superior to verapamil and rhodamine-123. Pairwise amino acid decomposition of superior docked conformations revealed Tyr303 as an important P-gp binding residue. DFT-based induced polarization analysis revealed major electrostatic fluctuations at the atomistic level and confirmed larger effects for amino acids with energy-favored binding interactions. Conformational analysis exhibited that auraptene and 7,4',7'',4'''-tetra-O-methylamentoflavone might not necessarily interact to P-gp binding sites through minimum energy conformations. Although there are still many hurdles to overcome, obtained results may propose a few nutritive phytochemicals as potential P-gp binding agents. Moreover; top-scored derivatives may have the chance to exhibit tumor chemo-sensitizing effects.

In silico analysis of SOD1 aggregation inhibition modes of tertiary amine pyrazolone and pyrano coumarin ferulate as ALS drug candidates.

Amyotrophic lateral sclerosis (ALS) is a fatal neurodegenerative disease, the familial form (fALS) of which is often cognate to mutations in the antioxidant enzyme Cu/Zn superoxide dismutase 1 (SOD1) leading to misfolding and aggregation. Two small molecules, a tertiary amine pyrazolone (TAP) and a pyrano coumarin ferulate (PCF) were suggested to be ALS drug candidates following experimental observation of their ability to inhibit SOD1 protein misfolding and aggregation. The present work aims at computational investigation of these experimentally proposed drug candidates to gain insight into their mechanism of SOD1 misfolding and aggregation inhibition. On the basis of molecular docking, molecular dynamics simulation, MM-PBSA and per-residue energy decomposition analysis, we examined the specific interactions of TAP and PCF with three probable binding sites of SOD1, namely, dimeric interface cavity, W32 and, UMP binding sites. Results suggest that the binding of TAP at W32 and at UMP sites are least probable due to absence of any favorable interaction. The binding of TAP to dimeric cavity is also unstable due to strong unfavorable interactions. In case of PCF, binding at the UMP site is least probable while binding at dimeric cavity is accompanied by unfavorable interactions. PCF, however, exhibits stable binding with the W32 binding site of SOD1 by stabilizing the solvent accessible hydrophobic residues, which otherwise would have acted as contact points for aggregation. Thus the results imply that compound PCF functions as an inhibitior of SOD1 misfolding/aggregation through direct interaction with the protein SOD1 at the W32 binding site. However, TAP is likely to act as an inhibitor through a different mechanism rather than direct interaction with the protein SOD1. These results apart from reinforcing previous experimental findings, shed light on the probable mechanism of action of the proposed drug candidates.

Atomistic insights into the binding of SARS-CoV-2 spike receptor binding domain with the human ACE2 receptor: The importance of residue 493

SARS-CoV-2 is a coronavirus that has created a global pandemic. The virus contains a spike protein which has been shown to bind to the ACE2 receptor on the surface of human cells. Vaccines have been developed that recognize elements of the SARS-CoV-2 spike protein and they have been successful in preventing infection. Recently, the Omicron variant of the SARS-CoV-2 virus was reported and quickly became a variant of concern due to its transmissibility. This variant contained an unusually large number (32) of point mutations, of which 15 of those mutations are in the receptor binding domain of the spike protein. While several computational and experimental investigations comparing the binding of the Omicron and wild type RBD to the human ACE2 receptor have been conducted, many of these report contradictory findings. In order to assess the differential binding ability, we conducted 2 μs of classical molecular dynamics (cMD) simulation to estimate the binding affinities and behaviors. Based upon MM-GBSA binding affinity, per-residue energy decomposition analysis, center of mass distance measurements, ensemble clustering, pairwise residue decomposition and hydrogen bonding analysis, our results suggest that a single point mutation is responsible for the enhanced binding of the Omicron mutant relative to the WT. While the 15-point mutations in the receptor binding domain contribute positively and negatively to the affinity of the spike protein for the human ACE2 receptor, it is the point mutation Q493R that confers enhanced binding while the Q493K mutation results in similar binding. The MM-GBSA binding estimations over a 2 μs trajectory, suggest that the wild type binds to ACE2 with a value of −29.69 kcal/mol while the Q493K and Q493R Omicron mutants bind with energy values of −26.67 and −34.56 kcal/mol, respectively. These values are significantly different, given the error estimates associated with the MM-GBSA method. In general, while some mutations increase binding, more mutations diminish binding, leading to an overall similar picture of binding for Q493K and enhanced binding for Q493R.

Identification of defactinib derivatives targeting focal adhesion kinase using ensemble docking, molecular dynamics simulations and binding free energy calculations

Focal adhesion kinase (FAK) belongs to the nonreceptor tyrosine kinases, which selectively phosphorylate tyrosine residues on substrate proteins. FAK is associated with bladder, esophageal, gastric, neck, breast, ovarian and lung cancers. Thus, FAK has been considered as a potential target for tumor treatment. Currently, there are six adenosine triphosphate (ATP)-competitive FAK inhibitors tested in clinical trials but no approved inhibitors targeting FAK. Defactinib (VS-6063) is a second-generation FAK inhibitor with an IC50 of 0.6 nM. The binding model of VS-6063 with FAK may provide a reference model for developing new antitumor FAK-targeting drugs. In this study, the VS-6063/FAK binding model was constructed using ensemble docking and molecular dynamics simulations. Furthermore, the molecular mechanics/generalized Born (GB) surface area (MM/GBSA) method was employed to estimate the binding free energy between VS-6063 and FAK. The key residues involved in VS-6063/FAK binding were also determined using per-residue energy decomposition analysis. Based on the binding model, VS-6063 could be separated into seven regions to enhance its binding affinity with FAK. Meanwhile, 60 novel defactinib-based compounds were designed and verified using ensemble docking. Overall, the present study improves our understanding of the binding mechanism of human FAK with VS-6063 and provides new insights into future drug designs targeting FAK. Communicated by Ramaswamy H. Sarma

Computational Insights into the Sequence-Activity Relationships of the NGF(1-14) Peptide by Molecular Dynamics Simulations.

The Nerve Growth Factor (NGF) belongs to the neurothrophins protein family involved in the survival of neurons in the nervous system. The interaction of NGF with its high-affinity receptor TrkA mediates different cellular pathways related to Alzheimer’s disease, pain, ocular dysfunction, and cancer. Therefore, targeting NGF-TrkA interaction represents a valuable strategy for the development of new therapeutic agents. In recent years, experimental studies have revealed that peptides belonging to the N-terminal domain of NGF are able to partly mimic the biological activity of the whole protein paving the way towards the development of small peptides that can selectively target specific signaling pathways. Hence, understanding the molecular basis of the interaction between the N-terminal segment of NGF and TrkA is fundamental for the rational design of new peptides mimicking the NGF N-terminal domain. In this study, molecular dynamics simulation, binding free energy calculations and per-residue energy decomposition analysis were combined in order to explore the molecular recognition pattern between the experimentally active NGF(1–14) peptide and TrkA. The results highlighted the importance of His4, Arg9 and Glu11 as crucial residues for the stabilization of NGF(1–14)-TrkA interaction, thus suggesting useful insights for the structure-based design of new therapeutic peptides able to modulate NGF-TrkA interaction.

Computational discovery of novel inhibitory candidates targeting versatile transcriptional repressor MBD2.

Genome methylation is a key epigenetic mechanism in various biological events such as development, cellular differentiation, cancer progression, aging, and iPSC reprogramming. Crosstalk between DNA methylation and gene expression is mediated by MBD2, known as the reader of DNA methylation and suggested as a drug target. Despite its magnitude of significance, a scarcely limited number of small molecules to be used as inhibitors have been detected so far. Therefore, we screened a comprehensive compound library to elicit MBD2 inhibitor candidates. Promising molecules were subjected to computational docking analysis by targeting the methylated DNA-binding domain of human MBD2. We could detect reasonable binding energies and docking residues, presumably located in druggable pockets. Docking results were also validated via MD simulation and per-residue energy decomposition calculation. Drug-likeness of these small molecules was assessed through ADMET prediction to foresee off-target side effects for future studies. All computational approaches notably highlighted two compounds named CID3100583 and 8,8-ethylenebistheophylline. These compounds have become prominent as novel candidates, possibly disrupting MBD2MBD-DNA interaction. Consequently, these compounds have been considered prospective inhibitors with the usage potential in a wide range of applications from cancer treatment to somatic cell reprogramming protocols.

Computational Design of Miniprotein Inhibitors Targeting SARS-CoV-2 Spike Protein.

The ongoing pandemic of COVID-19 caused by SARS-CoV-2 has become a global health problem. There is an urgent need to develop therapeutic drugs, effective therapies, and vaccines to prevent the spread of the virus. The virus first enters the host cell through the interaction between the receptor binding domain (RBD) of spike protein and the peptidase domain (PD) of the angiotensin-converting enzyme 2 (ACE2). Therefore, blocking the binding of RBD and ACE2 is a promising strategy to inhibit the invasion and infection of the virus in the host cell. In the study, we designed several miniprotein inhibitors against SARS-CoV-2 by single/double/triple-point mutant, based on the initial inhibitor LCB3. Molecular dynamics (MD) simulations and trajectory analysis were performed for an in-depth analysis of the structural stability, essential protein motions, and per-residue energy decomposition involved in the interaction of inhibitors with the RBD. The results showed that the inhibitors have adapted the protein RBD in the binding interface, thereby forming stable complexes. These inhibitors display low binding free energy in the MM/PBSA calculations, substantiating their strong interaction with RBD. Moreover, the binding affinity of the best miniprotein inhibitor, H6Y-M7L-L17F mutant, to RBD was ∼45 980 times (ΔG = RT ln Ki) higher than that of the initial inhibitor LCB3. Following H6Y-M7L-L17F mutant, the inhibitors with strong binding activity are successively H6Y-L17F, L17F, H6Y, and F30Y mutants. Our research proves that the miniprotein inhibitors can maintain their secondary structure and have a highly stable blocking (binding) effect on SARS-CoV-2. This study proposes novel miniprotein mutant inhibitors with enhanced binding to spike protein and provides potential guidance for the rational design of new SARS-CoV-2 spike protein inhibitors.

FastDRH: a webserver to predict and analyze protein-ligand complexes based on molecular docking and MM/PB(GB)SA computation.

Predicting the native or near-native binding pose of a small molecule within a protein binding pocket is an extremely important task in structure-based drug design, especially in the hit-to-lead and lead optimization phases. In this study, fastDRH, a free and open accessed web server, was developed to predict and analyze protein-ligand complex structures. In fastDRH server, AutoDock Vina and AutoDock-GPU docking engines, structure-truncated MM/PB(GB)SA free energy calculation procedures and multiple poses based per-residue energy decomposition analysis were well integrated into a user-friendly and multifunctional online platform. Benefit from the modular architecture, users can flexibly use one or more of three features, including molecular docking, docking pose rescoring and hotspot residue prediction, to obtain the key information clearly based on a result analysis panel supported by 3Dmol.js and Apache ECharts. In terms of protein-ligand binding mode prediction, the integrated structure-truncated MM/PB(GB)SA rescoring procedures exhibit a success rate of >80% in benchmark, which is much better than the AutoDock Vina (~70%). For hotspot residue identification, our multiple poses based per-residue energy decomposition analysis strategy is a more reliable solution than the one using only a single pose, and the performance of our solution has been experimentally validated in several drug discovery projects. To summarize, the fastDRH server is a useful tool for predicting the ligand binding mode and the hotspot residue of protein for ligand binding. The fastDRH server is accessible free of charge at http://cadd.zju.edu.cn/fastdrh/.

An Effective Chemo-Enzymatic method with An Evolved L-Threonine Aldolase for Preparing L-threo-4-Methylsulfonylphenylserine Ethyl Ester of High Optical Purity

L-threo-4-methylsulfonylphenylserine ethyl ester [(2S, 3R)-1d], a key building block of florfenicol which synthesized by esterification using ethanol and L-threo-4-methylsulfonylphenylserine [(2S, 3R)-1b] as substrates. L-threonine aldolase (LTA) is a promising biocatalyst for producing (2S, 3R)-1b through a one-step process taking 4-Methylsulphonyl benzaldehyde (1a) and glycine as substrates under the mild condition. However, the moderate Cβ-stereoselective blocked the industrial application of LTA. To address this issue, rational design by combination with prereaction state molecular dynamics (MD) simulation and per-residue energy decomposition algorithm are employed to engineer LTA for enhancing Cβ-stereoselective. As a result, a triple mutant N16A/E98S/Y314R (Mu3) is screened out to produce (2S, 3R)-1b with 93.7% de and 90.2% conversion under the 300 mM 1a substrate loading. Furthermore, esterification is applied to synthesize (2S, 3R)-1d with over 99% chiral purity and 99% product purity. The success of this study provides new insights for the rational design of LTA with improved Cβ-stereoselectivity and proves a chemo-enzymatic route for green synthesis of (2S, 3R)-1d.

Application of per-Residue Energy Decomposition to Design Peptide Inhibitors of PSD95 GK Domain.

Specific interaction between the postsynaptic density protein 95 (PSD95) and synapse-associated protein 90/postsynaptic density 95–associated protein (SAPAP) is crucial for excitatory synaptic development and plasticity. Designing inhibitors that target the guanylate kinase (GK) domain of PSD95, which is responsible for the interaction, is a promising manipulation tool for the investigation of the function of PSD95 GK and the etiology of its related psychiatric disorders. Herein, we designed new peptide inhibitors of PSD95 GK/SAPAP with higher binding affinity by using molecular dynamics simulations. First, the interactions between PSD95 GK and their reported phosphorylated and unphosphorylated peptides were explored by molecular dynamics simulations. Besides the hydrogen bonding interactions mediated by the phospho-serine (p-Ser) or corresponding phosphomimic residue Asp/Glu, the hydrophobic interactions from the other amino acids also contribute to the PSD95 GK/SAPAP interaction. As an unphosphorylated synthetic peptide with moderate binding affinity and relatively lower molecular weight, the QSF inhibitory peptide was selected for further modification. Based on per-residue energy decomposition results of the PSD95 GK/QSF complex, ten peptides were designed to enhance the binding interactions, especially the hydrophobic interactions. The top-ranked five peptides with lower binding energy were eventually synthesized. The binding affinities of the synthesized peptides were determined using fluorescence polarization (FP) assay. As expected, all peptides have higher binding affinity than the QSF peptide (K i = 5.64 ± 0.51 μM). Among them, F10W was the most potent inhibitor (K i = 0.75 ± 0.25 μM), suggesting that enhancement of the hydrophobic interactions is an important strategy for the design of new inhibitory peptides targeting PSD95 GK.