Benchmarked molecular docking integrated molecular dynamics stability analysis for prediction of SARS-CoV-2 papain-like protease inhibition by olive secoiridoids

Objectives

We performed a virtual screening of olive secoiridoids of the OliveNet^TM library to predict SARS-CoV-2 PLpro inhibition. Benchmarked molecular docking protocol that evaluated the performance of two docking programs was applied to execute virtual screening. Molecular dynamics stability analysis of the top-ranked olive secoiridoid docked to PLpro was also carried out.

Methods

Benchmarking virtual screening used two freely available docking programs, AutoDock Vina 1.1.2. and AutoDock 4.2.1. for molecular docking of olive secoiridoids to a single PLpro structure. Screening also included benchmark structures of known active and decoy molecules from the DEKOIS 2.0 library. Based on the predicted binding energies, the docking programs ranked the screened molecules. We applied the usual performance evaluation metrices to evaluate the docking programs using the predicted ranks. Molecular dynamics of the top-ranked olive secoiridoid bound to PLpro and computation of MM-GBSA energy using three iterations during the last 50 ps of the analysis of the dynamics in Desmond supported the stability prediction.

Results and discussions

Predictiveness curves suggested that AutoDock Vina has a better predictive ability than AutoDock, although there was a moderate correlation between the active molecules rankings (Kendall's correlation of rank (τ) = 0.581). Interestingly, two same molecules, Demethyloleuropein aglycone, and Oleuroside enriched the top 1 % ranked olive secoiridoids predicted by both programs. Demethyloleuropein aglycone bound to PLpro obtained by docking in AutoDock Vina when analyzed for stability by molecular dynamics simulation for 50 ns displayed an RMSD, RMSF<2 Å, and MM-GBSA energy of -94.54 ± 6.05 kcal/mol indicating good stability. Molecular dynamics also revealed the interactions of Demethyloleuropein aglycone with binding sites 2 and 3 of PLpro, suggesting a potent inhibition. In addition, for 98 % of the simulation time, two phenolic hydroxy groups of Demethyloleuropein aglycone maintained two hydrogen bonds with Asp302 of PLpro, specifying the significance of the groups in receptor binding.

Conclusion

AutoDock Vina retrieved the active molecules accurately and predicted Demethyloleuropein aglycone as the best inhibitor of PLpro. The Arabian diet consisting of olive products rich in secoiridoids benefits from the PLpro inhibition property and reduces the risk of viral infection.

PP2A is activated by cytochrome c upon formation of a diffuse encounter complex with SET/TAF-Iβ

•

Cytochrome c activates PP2A by disrupting SET/TAF-Iβ-mediated PP2A inhibition.

•

The “earmuff” domains of SET/TAF-Iβ are sufficient to bind PP2A or cytochrome c.

•

SET/TAF-Iβ ΔC exhibits high flexibility and intrinsic dynamics.

•

Cytochrome c and SET/TAF-Iβ ΔC form a polyconformational encounter complex.

Intrinsic protein flexibility is of overwhelming relevance for intermolecular recognition and adaptability of highly dynamic ensemble of complexes, and the phenomenon is essential for the understanding of numerous biological processes. These conformational ensembles—encounter complexes—lack a unique organization, which prevents the determination of well-defined high resolution structures. This is the case for complexes involving the oncoprotein SET/template-activating factor-Iβ (SET/TAF-Iβ), a histone chaperone whose functions and interactions are significantly affected by its intrinsic structural plasticity. Besides its role in chromatin remodeling, SET/TAF-Iβ is an inhibitor of protein phosphatase 2A (PP2A), which is a key phosphatase counteracting transcription and signaling events controlling the activity of DNA damage response (DDR) mediators. During DDR, SET/TAF-Iβ is sequestered by cytochrome c (Cc) upon migration of the hemeprotein from mitochondria to the cell nucleus. Here, we report that the nuclear SET/TAF-Iβ:Cc polyconformational ensemble is able to activate PP2A. In particular, the N-end folded, globular region of SET/TAF-Iβ (a.k.a. SET/TAF-Iβ ΔC)—which exhibits an unexpected, intrinsically highly dynamic behavior—is sufficient to be recognized by Cc in a diffuse encounter manner. Cc-mediated blocking of PP2A inhibition is deciphered using an integrated structural and computational approach, combining small-angle X-ray scattering, electron paramagnetic resonance, nuclear magnetic resonance, calorimetry and molecular dynamics simulations.

In silico epitope-based vaccine design against influenza a neuraminidase protein: Computational analysis established on B- and T-cell epitope predictions

Objective

Influenza A virus belongs to the most studied virus and its mutant initiates epidemic and pandemics outbreaks. Inoculation is the significant foundation to diminish the risk of infection. To prevent an incidence of influenza from the transmission, various practical approaches require more advancement and progress. More efforts and research must take in front to enhance vaccine efficacy.

Methods

The present research emphasizes the development and expansion of a universal vaccine for the influenza virus. Research focuses on vaccine design with high efficacy. In this study, numerous computational approaches were used, covering a wide range of elements and ideas in bioinformatics methodology. Various B and T-cell epitopic peptides derived from the Neuraminidase protein N1 are recognized by these approaches. With the implementation of numerous obtained databases and bioinformatics tools, the different immune framework methods of the conserved sequences of N1 neuraminidase were analyzed. NCBI databases were employed to retrieve amino acid sequences. The antigenic nature of the neuraminidase sequence was achieved by the VaxiJen server and Kolaskar and Tongaonkar method. After screening of various B and T cell epitopes, one efficient peptide each from B cell epitope and T cell epitopes was assessed for their antigenic determinant vaccine efficacy. Identical two B cell epitopes were recognized from the N1 protein when analyzed using B-cell epitope prediction servers. The detailed examination of amino acid sequences for interpretation of B and T cell epitopes was achieved with the help of the ABCPred and Immune Epitope Database.

Results

Computational immunology via immunoinformatic study exhibited RPNDKTG as having its high conservancy efficiency and demonstrated as a good antigenic, accessible surface hydrophilic B-cell epitope. Among T cell epitope analysis, YVNISNTNF was selected for being a conserved epitope. T cell epitope was also analyzed for its allergenicity and cytotoxicity evaluation. YVNISNTNF epitope was found to be a non-allergen and not toxic for cells as well. This T-cell epitope with maximum world populace coverages was scrutinized for its association with the HLA-DRB1*0401 molecule. Results from docking simulation analyses showed YVNISNTNF having lower binding energy, the radius of gyration (Rg), RMSD values, and RMSE values which make the protein structure more stable and increase its ability to become an epitopic peptide for influenza virus vaccination.

Conclusions

We propose that this epitope analysis may be successfully used as a measurement tool for the robustness of an antigen-antibody reaction between mutant strains in the annual design of the influenza vaccine.

The fatty acid site is coupled to functional motifs in the SARS-CoV-2 spike protein and modulates spike allosteric behaviour

The SARS-CoV-2 spike protein is the first contact point between the SARS-CoV-2 virus and host cells and mediates membrane fusion. Recently, a fatty acid binding site was identified in the spike (Toelzer et al. Science 2020). The presence of linoleic acid at this site modulates binding of the spike to the human ACE2 receptor, stabilizing a locked conformation of the protein. Here, dynamical-nonequilibrium molecular dynamics simulations reveal that this fatty acid site is coupled to functionally relevant regions of the spike, some of them far from the fatty acid binding pocket. Removal of a ligand from the fatty acid binding site significantly affects the dynamics of distant, functionally important regions of the spike, including the receptor-binding motif, furin cleavage site and fusion-peptide-adjacent regions. Simulations of the D614G mutant show differences in behaviour between these clinical variants of the spike: the D614G mutant shows a significantly different conformational response for some structural motifs relevant for binding and fusion. The simulations identify structural networks through which changes at the fatty acid binding site are transmitted within the protein. These communication networks significantly involve positions that are prone to mutation, indicating that observed genetic variation in the spike may alter its response to linoleate binding and associated allosteric communication.

Neuroprotective derivatives of tacrine that target NMDA receptor and acetyl cholinesterase - Design, synthesis and biological evaluation

The complex and multifactorial nature of neuropsychiatric diseases demands multi-target drugs that can intervene with various sub-pathologies underlying disease progression. Targeting the impairments in cholinergic and glutamatergic neurotransmissions with small molecules has been suggested as one of the potential disease-modifying approaches for Alzheimer’s disease (AD). Tacrine, a potent inhibitor of acetylcholinesterase (AChE) is the first FDA approved drug for the treatment of AD. Tacrine is also a low affinity antagonist of N-methyl-D-aspartate receptor (NMDAR). However, tacrine was withdrawn from its clinical use later due to its hepatotoxicity. With an aim to develop novel high affinity multi-target directed ligands (MTDLs) against AChE and NMDAR, with reduced hepatotoxicity, we performed in silico structure-based modifications on tacrine, chemical synthesis of the derivatives and in vitro validation of their activities. Nineteen such derivatives showed inhibition with IC₅₀ values in the range of 18.53 ± 2.09 – 184.09 ± 19.23 nM against AChE and 0.27 ± 0.05 – 38.84 ± 9.64 μM against NMDAR. Some of the selected compounds also protected rat primary cortical neurons from glutamate induced excitotoxicity. Two of the tacrine derived MTDLs, 201 and 208 exhibited in vivo efficacy in rats by protecting against behavioral impairment induced by administration of the excitotoxic agent, monosodium glutamate. Additionally, several of these synthesized compounds also exhibited promising inhibitory activitiy against butyrylcholinesterase. MTDL-201 was also devoid of hepatotoxicity in vivo. Given the therapeutic potential of MTDLs in disease-modifying therapy, our studies revealed several promising MTDLs among which 201 appears to be a potential candidate for immediate preclinical evaluations.

A two-stage computational approach to predict novel ligands for a chemosensory receptor

Olfactory receptor (OR) 1A2 is the member of largest superfamily of G protein-coupled receptors (GPCRs). OR1A2 is an ectopically expressed receptor with only 13 known ligands, implicated in reducing hepatocellular carcinoma progression, with enormous therapeutic potential. We have developed a two-stage screening approach to identify novel putative ligands of OR1A2. We first used a pharmacophore model based on atomic property field (APF) to virtually screen a library of 5942 human metabolites. We then carried out structure-based virtual screening (SBVS) for predicting the potential agonists, based on a 3D homology model of OR1A2. This model was developed using a biophysical approach for template selection, based on multiple parameters including hydrophobicity correspondence, applied to the complete set of available GPCR structures to pick the most appropriate template. Finally, the membrane-embedded 3D model was refined by molecular dynamics (MD) simulations in both the apo and holo forms. The refined model in the apo form was selected for SBVS. Four novel small molecules were identified as strong binders to this olfactory receptor on the basis of computed binding energies.

Screening of inhibitors against SARS-CoV-2 spike protein and their capability to block the viral entry mechanism: A viroinformatics study

SARS-CoV-2, previously named 2019 novel coronavirus (2019-nCoV), has been associated with the global pandemic of acute respiratory distress syndrome. First reported in December 2019 in the Wuhan province of China, this new RNA virus has several folds higher transmission among humans than its other family member (SARS-CoV and MERS-CoV). The SARS-CoV-2 spike receptor-binding domain (RBD) is the region mediating the binding of the virus to host cells via Angiotensin-converting enzyme 2 (ACE2), a critical step of viral. Here in this study, we have utilized in silico approach for the virtual screening of antiviral library extracted from the Asinex database against the Receptor binding domain (RBD) of the S1 subunit of the SARS-CoV-2 spike glycoprotein. Further, the molecules were ranked based on their binding affinity against RBD, and the top 15 molecules were selected. The affinity of these selected molecules to interrupt the ACE2-Spike interaction was also studied. It was found that the chosen molecules were demonstrating excellent binding affinity against spike protein, and these molecules were also very effectively interrupting the ACE2-RBD interaction. Furthermore, molecular dynamics (MD) simulation studies were utilized to investigate the top 3 selected molecules' stability in the ACE2-RBD complexes. To the best of our knowledge, this is the first study where molecules' inhibitory potential against the Receptor binding domain (RBD) of the S1 subunit of the SARS-CoV-2 spike glycoprotein and their inhibitory potential against the ACE2-Spike has been studied. We believe that these compounds can be further tested as a potential therapeutic option against COVID-19.

Nature's therapy for COVID-19: Targeting the vital non-structural proteins (NSP) from SARS-CoV-2 with phytochemicals from Indian medicinal plants

BACKGROUND

Containing COVID-19 is still a global challenge. It has affected the "normal" world by targeting its economy and health sector. The effect is shifting of focus of research from life threatening diseases like cancer. Thus, we need to develop a medical solution at the earliest. The purpose of this present work was to understand the efficacy of 22 rationally screened phytochemicals from Indian medicinal plants obtained from our previous work, following drug-likeness properties, against 6 non-structural-proteins (NSP) from SARS-CoV-2.

METHODS

100 ns molecular dynamics simulations were performed, and relative binding free energies were computed by MM/PBSA. Further, principal component analysis, dynamic cross correlation and hydrogen bond occupancy were analyzed to characterize protein-ligand interactions. Biological pathway enrichment analysis was also carried out to elucidate the therapeutic targets of the phytochemicals in comparison to SARS-CoV-2.

RESULTS

The potential binding modes and favourable molecular interaction profile of 9 phytochemicals, majorly from Withania somnifera with lowest free binding energies, against the SARS-CoV-2 NSP targets were identified. It was understood that phytochemicals and 2 repurposed drugs with steroidal moieties in their chemical structures formed stable interactions with the NSPs. Additionally, human target pathway analysis for SARS-CoV-2 and phytochemicals showed that cytokine mediated pathway and phosphorylation pathways were with the most significant p-value.

CONCLUSIONS

To summarize this work, we suggest a global approach of targeting multiple proteins of SARS-CoV-2 with phytochemicals as a natural alternative therapy for COVID-19. We also suggest that these phytochemicals need to be tested experimentally to confirm their efficacy.

Determining the unbinding events and conserved motions associated with the pyrazinamide release due to resistance mutations of Mycobacterium tuberculosis pyrazinamidase

•

WT and 82 Mtb PZase-PZA complexes were established using molecular modelling.

•

AMBER force field parameters for the PZase Fe²⁺ cofactor were established.

•

Statistically guided network analysis described late PZA unbinding events.

•

MBS and nearby residue mutations lead to premature PZA release from PZase.

•

Cofactor destabilisation was found in only some PZase mutants.

Pyrazinamide (PZA) is the only first-line antitubercular drug active against latent Mycobacterium tuberculosis (Mtb). It is activated to pyrazinoic acid by the pncA-encoded pyrazinamidase enzyme (PZase). Despite the emergence of PZA drug resistance, the underlying mechanisms of resistance remain unclear. This study investigated part of these mechanisms by modelling a PZA-bound wild type and 82 mutant PZase structures before applying molecular dynamics (MD) with an accurate Fe²⁺ cofactor coordination geometry. After observing nanosecond-scale PZA unbinding from several PZase mutants, an algorithm was developed to systematically detect ligand release via centre of mass distances (COM) and ligand average speed calculations, before applying the statistically guided network analysis (SGNA) method to investigate conserved protein motions associated with ligand unbinding. Ligand and cofactor perspectives were also investigated. A conserved pair of lid-destabilising motions was found. These consisted of (1) antiparallel lid and side flap motions; (2) the contractions of a flanking region within the same flap and residue 74 towards the core. Mutations affecting the hinge residues (H51 and H71), nearby residues or L19 were found to destabilise the lid. Additionally, other metal binding site (MBS) mutations delocalised the Fe²⁺ cofactor, also facilitating lid opening. In the early stages of unbinding, a wider variety of PZA poses were observed, suggesting multiple exit pathways. These findings provide insights into the late events preceding PZA unbinding, which we found to occur in some resistant PZase mutants. Further, the algorithm developed here to identify unbinding events coupled with SGNA can be applicable to other similar problems.

Oxidation State Dependent Conformational Changes of HMGB1 Regulate the Formation of the CXCL12/HMGB1 Heterocomplex

High-mobility Group Box 1 (HMGB1) is an abundant protein present in all mammalian cells and involved in several processes. During inflammation or tissue damage, HMGB1 is released in the extracellular space and, depending on its redox state, can form a heterocomplex with CXCL12. The heterocomplex acts exclusively via the chemokine receptor CXCR4 enhancing leukocyte recruitment.

Here, we used multi-microsecond molecular dynamics (MD) simulations to elucidate the effect of the disulfide bond on the structure and dynamics of HMGB1.

The results of the MD simulations show that the presence or lack of the disulfide bond between Cys23 and Cys45 modulates the conformational space explored by HMGB1, making the reduced protein more suitable to form a complex with CXCL12.

Computational Simulations Identified Two Candidate Inhibitors of Cdk5/p25 to Abrogate Tau-associated Neurological Disorders

Deregulation of Cdk5 is a hallmark in neurodegenerative diseases and its complex with p25 forms Cdk5/p25, thereby causes severe neuropathological insults. Cdk5/p25 abnormally phosphorylates tau protein, and induces tau-associated neurofibrillary tangles in neurological disorders. Therefore, the pharmacological inhibition of Cdk5/p25 alleviates tau-associated neurological disorders. Herein, computational simulations probed two candidate inhibitors of Cdk5/p25. Structure-based pharmacophore investigated the essential complementary chemical features of ATP-binding site of Cdk5 in complex with roscovitine. Resultant pharmacophore harbored polar interactions with Cys83 and Asp86 residues and non-polar interactions with Ile10, Phe80, and Lys133 residues of Cdk5. The chemical space of selected pharmacophore was comprised of two hydrogen bond donors, one hydrogen bond acceptor, and three hydrophobic features. Decoy test validation of pharmacophore obtained highest Guner-Henry score (0.88) and enrichment factor score (7.23). The screening of natural product drug-like databases by validated pharmacophore retrieved 1126 compounds as candidate inhibitors of Cdk5/p25. The docking of candidate inhibitors filtered 10 molecules with docking score >80.00 and established polar and non-polar interactions with the ATP-binding site residues of Cdk5/p25. Finally, molecular dynamics simulation and binding free energy analyses identified two candidate inhibitors of Cdk5/p25. During 30 ns simulation, the candidate inhibitors established <3.0 Å root mean square deviation and stable hydrogen bond interactions with the ATP-binding site residues of Cdk5/p25. The final candidate inhibitors obtained lowest binding free energies of -122.18 kJ/mol and - 117.26 kJ/mol with Cdk5/p25. Overall, we recommend two natural product candidate inhibitors to target the pharmacological inhibition of Cdk5/p25 in tau-associated neurological disorders.

Asymmetric Spontaneous Intercalation of Lutein into a Phospholipid Bilayer, a Computational Study

Lutein, a hydroxylated carotenoid, is a pigment synthesised by plants and bacteria. Animals are unable to synthesise lutein, nevertheless, it is present in animal tissues, where its only source is dietary intake. Both in plants and animals, carotenoids are associated mainly with membranes where they carry out important physiological functions. Their trafficking to and insertion into membranes are not well recognised due to experimental difficulties. In this paper, a computational approach is used to elucidate details of the dynamics and energetics of lutein intercalation from the water to the phospholipid bilayer phase. The dynamics is studied using molecular dynamics simulation, and the energetics using umbrella sampling. Lutein spontaneous insertion into the bilayer and translocation across it proceed via formation of hydrogen bonds between its hydroxyl groups and water and/or phospholipid oxygen atoms as well as desolvation of its polyene chain. As lutein molecule is asymmetric, its bilayer intercalation is also asymmetric. The course of events and timescale of the intercalation are different from those of helical peptides. The time of full lutein intercalation ranges from 20 to 100 ns and its final orientation is predominately vertical. Nevertheless, some lutein molecules are in the final horizontal position and some aggregate in the water phase and remain there for the whole simulation time. The highest energy barrier for the intercalation process is ~2.2 kcal/mol and the energy gain is ~18 kcal/mol. The results obtained for lutein can be applied to other xanthophylls and molecules of a similar structure.

•

Lutein as an amphiphilic molecule readily intercalates into a lipid bilayer.

•

MD simulations with dense sampling revealed details of the intercalation process.

•

Time of full lutein membrane intercalation ranges from 20 to 100 ns.

•

Lutein intercalation from the ε ring end is less probable than from the β ring end.

•

Horizontal position of lutein in the bilayer is less probable than vertical.

An Energetically Favorable Ligand Entrance Gate of a Multidrug Transporter Revealed by Partial Nudged Elastic Band Simulations

P-glycoprotein (P-gp) is a multidrug transporter, which harnesses the chemical energy of ATP to power the efflux of diverse chemotherapeutics out of cells and thus contributes to the development of multidrug resistance (MDR) in cancer. It has been proved that the ligand-binding pocket of P-gp is located at the transmembrane domains (TMDs). However, the access of ligands into the binding pocket remains to be elucidated, which definitely hinder the development of P-gp inhibitors. Herein, the access pathways of a well-known substrate rhodamine-123 and a cyclopeptide inhibitor QZ-Leu were characterized by time-independent partial nudged elastic band (PNEB) simulations. The decreasing free energies along the PNEB-optimized access pathway indicated that TM4/6 cleft may be an energetically favorable entrance gate for ligand entry into the binding pocket of P-gp. The results can be reconciled with a range of experimental studies, further corroborating the reliability of the gate revealed by computational simulations. Our atomic level description of the ligand access pathway provides valuable insights into the gating mechanism for drug uptake and transport by P-gp and other multidrug transporters.

•

P-gp contributes to the development of multidrug resistance in cancer.

•

The entrance of drugs into P-gp binding pocket has yet to be elucidated.

•

An energetically favorable entrance gate was revealed by PNEB simulations.

•

The computational results were reconciled with the experimental data.

•

The atomic simulations provide insights into the gating mechanism of P-gp.

Interacting mechanism of ID3 HLH domain towards E2A/E12 transcription factor - An Insight through molecular dynamics and docking approach

Inhibitor of DNA binding protein 3 (ID3) has long been characterized as an oncogene that implicates its functional role through its Helix-Loop-Helix (HLH) domain upon protein-protein interaction. An insight into the dimerization brought by this domain helps in identifying the key residues that favor the mechanism behind it. Molecular dynamics (MD) simulations were performed for the HLH proteins ID3 and Transcription factor E2-alpha (E2A/E12) and their ensemble complex (ID3-E2A/E12) to gather information about the HLH domain region and its role in the interaction process. Further evaluation of the results by Principal Component Analysis (PCA) and Free Energy Landscape (FEL) helped in revealing residues of E2A/E12: Lys570, Ala595, Val598, and Ile599 and ID3: Glu53, Gln63, and Gln66 buried in their HLH motifs imparting key roles in dimerization process. Furthermore the T-pad analysis results helped in identifying the key fluctuations and conformational transitions using the intrinsic properties of the residues present in the domain region of the proteins thus specifying their crucial role towards molecular recognition. The study provides an insight into the interacting mechanism of the ID3-E2A/E12 complex and maps the structural transitions arising in the essential conformational space indicating the key structural changes within the helical regions of the motif. It thereby describes how the internal dynamics of the proteins might regulate their intrinsic structural features and its subsequent functionality.

Low-mass molecular dynamics simulation for configurational sampling enhancement: More evidence and theoretical explanation

It has been reported recently that classical, isothermal-isobaric molecular dynamics (NTP MD) simulations at a time step of 1.00 fs of the standard-mass time (Δt=1.00 fssmt) and a temperature of ≤340 K using uniformly reduced atomic masses by tenfold offers better configurational sampling than standard-mass NTP MD simulations at the same time step. However, it has long been reported that atomic masses can also be increased to improve configurational sampling because higher atomic masses permit the use of a longer time step. It is worth investigating whether standard-mass NTP MD simulations at Δt=2.00 or 3.16 fssmt can offer better or comparable configurational sampling than low-mass NTP MD simulations at Δt=1.00 fssmt. This article reports folding simulations of two β-hairpins showing that the configurational sampling efficiency of NTP MD simulations using atomic masses uniformly reduced by tenfold at Δt=1.00 fssmt is statistically equivalent to and better than those using standard masses at Δt=3.16 and 2.00 fssmt, respectively. The results confirm that, relative to those using standard masses at routine Δt=2.00 fssmt, the low-mass NTP MD simulations at Δt=1.00 fssmt are a simple and generic technique to enhance configurational sampling at temperatures of ≤340 K.