Markov State Modeling Analysis Captures Changes in the Temperature-Sensitive N-Terminal and ß-Turn Regions of the p53 DNA-Binding Domain.

The transcription factor p53 is one of the most widely studied cancer proteins. Its temperature-sensitive nature suggests reduction in functionality at physiological temperatures. Temperature-induced conformational variations and their impact on its functional ability still remain unexplored. A total of 20.8 µs molecular dynamics simulations of wildtype p53 in the apo and the DNA-bound states have been performed at 300 K and 310 K. Further, Markov State Modeling (MSM) analyses were performed, considering Cα-Cα distances as reaction coordinates. Filtering of these distances based on correlation with the time-independent components (tICs) resulted in 16 and 32 distances for apo and DNA-bound systems, respectively. Individual MSM analyses using these filtered distances were performed for both p53 systems. These Cα-Cα residue pairs belonged to the N-terminal, S6/7 ß-turn, loop L2, loop L3, and hydrophobic core residues. At physiological temperatures, apo-p53 exhibits exposure of its hydrophobic core, where the temperature-sensitive hotspot residues were also located. This exposure was the result of the S6/7 ß-turn and N-terminal moving apart. In the DNA-bound p53 system, loop L1 attains an open conformation at physiological temperatures, which weakens the DNA binding. It is already known that p53 mutants that lack DNA binding also tend to show similar conformational variations. The S6/7 ß-turn along with the already known functionally important loop L2 may pose as regions to be targeted to overcome the loss in binding of temperature-sensitive wildtype p53. Rescue strategies directed toward these temperature-sensitive regions may be useful to recuperate its strong binding at physiological temperatures.

An insight into the inhibitory mechanism of phytochemicals and FDA-approved drugs on the ACE2-Spike complex of SARS-CoV-2 using computational methods.

The S-glycoprotein (Spike) of the SARS-CoV-2 forms a complex with the human transmembrane protein angiotensin-converting enzyme 2 (ACE2) during infection. It forms the first line of contact with the human cell. The FDA-approved drugs and phytochemicals from Indian medicinal plants were explored. Molecular docking and simulations of these molecules targeting the ACE2-Spike complex were performed. Rutin DAB10 and Swertiapuniside were obtained as the top-scored drugs as per the docking protocol. The MD simulations of ligand-free, Rutin DAB10-bound, and Swertiapuniside-bound ACE2-Spike complex revealed abrogation of the hydrogen bonding network between the two proteins. The principal component and dynamic cross-correlation analysis pointed out conformational changes in both the proteins unique to the ligand-bound systems. The interface residues, His34, and Lys353 from ACE2 and Arg403, and Tyr495 from the Spike protein formed significant strong interactions with the ligand molecules, inferring the inhibition of ACE2-Spike complex. Few novel interactions specific to Rutin-DAB10 and Swertiapuniside were also identified. The conformational flexibility of the drug-binding pocket was captured using the RMSD-based clustering of the ligand-free simulations. Ensemble docking was performed wherein the FDA-approved database and phytochemical dataset were docked on each of the cluster representatives of the ACE2-Spike. The phytochemicals identified belonged to Withania somnifera, Swertia chirayita, Tinospora cordifolia and Rutin DAB10, fulvestrant, elbasvir from FDA. Supplementary Information: The online version contains supplementary material available at 10.1007/s11696-021-01680-1.

TANGO: A high through-put conformation generation and semiempirical method-based optimization tool for ligand molecules.

Lead optimization is one of the crucial steps in the drug discovery pipeline. After identifying the lead molecule and obtaining its 2D geometry, understanding the best conformation it would attain in 3D still remains one of the most challenging steps in drug discovery. There have been multiple methods and algorithms that are directed toward achieving best conformation for the lead molecules. TANGO focuses on conformation generation and its optimization using semiempirical energy calculations. The conformation generation is based on torsion angle rotation of the exocyclic bonds. The energy calculations are performed using MOPAC. The unique feature of this tool lies in the implementation of Message Passing Interface (MPI) for conformation generation and semiempirical-based optimization. A well-defined architecture handling the input and output generation has been used. The master and slave approach to handle operations involved in torsion angle rotation and energy calculations has helped in load balancing the process of conformation generation. The benchmarking results suggest that TANGO scales significantly well across eight nodes with each node utilizing 16 cores. This tool may prove to very useful in high throughput generation of semiempirically optimized small molecule conformations. The use of semiempirical methods for optimization generates a conformational ensemble thereby helping to obtain stable and alternate stable conformers for a given ligand molecule. © 2018 Wiley Periodicals, Inc.

Acetylcholinesterase and Aß Aggregation Inhibition by Heterometallic Ruthenium(II)-Platinum(II) Polypyridyl Complexes.

Two heteronuclear ruthenium(II)-platinum(II) complexes [Ru(bpy)2(BPIMBp)PtCl2]2+ (3) and [Ru(phen)2(BPIMBp)PtCl2]2+ (4), where bpy = 2,2'-bipyridine, phen = 1,10-phenanthroline, and BPIMBp = 1,4'-bis[(2-pyridin-2-yl)-1H-imidazol-1-ylmethyl]-1,1'-biphenyl, have been designed and synthesized from their mononuclear precursors [Ru(bpy)2(BPIMBp)]2+ (1) and [Ru(phen)2(BPIMBp)]2+ (2) as multitarget molecules for Alzheimer's disease (AD). The inclusion of the cis-PtCl2 moiety facilitates the covalent interaction of Ru(II) polypyridyl complexes with amyloid ß (Aß) peptide. These multifunctional complexes act as inhibitors of acetylcholinesterase (AChE), Aß aggregation, and Cu-induced oxidative stress and protect neuronal cells against Aß-toxicity. The study highlights the design of metal based anti-Alzheimer's disease (AD) systems.

Auto-regulation of Slug mediates its activity during epithelial to mesenchymal transition.

Slug, a five C2H2 zinc finger (ZF) motif transcription factor mediates cell migration in development, adult tissue repair and regeneration, as well as during tumor metastases through epithelial to mesenchymal transition. At the molecular level, this involves interactions with E-box (CACC/GGTG) consensus elements within target gene promoters to achieve transcriptional repression. However, precise elucidation of events involved in this DNA recognition and binding of specific promoters to regulate target genes have not been achieved. In the present study, we show that besides transcriptional repression, Slug can also directly activate its own expression by preferential binding to specific E-box elements in the distal binding region of its promoter. Our findings suggest that while the first ZF does not contribute to the transcription-associated functions of Slug, all the remaining four ZFs are involved in regulating the expression of target genes with ZF3 and ZF4 being more crucial than ZF2 or ZF5. We also report that recognition and binding preferences of ZFs are defined through intrinsic differences in the E-box core base pairs and/or flanking sequences, with the S2 E-box element being most critical during autoregulation. However, specific target E-box recognition and binding are also defined by the cellular context, which implies that in silico and/or biochemical DNA binding preferences may not necessarily be able to accurately predict in situ events. Our studies thus constitute a novel understanding of transcriptional regulation.

Traversing the folding pathway of proteins using temperature-aided cascade molecular dynamics with conformation-dependent charges.

Protein folding is a multi-micro second time scale event and involves many conformational transitions. Crucial conformational transitions responsible for biological functions of biomolecules are difficult to capture using current state-of-the-art molecular dynamics (MD) simulations. Protein folding, being a stochastic process, witnesses these transitions as rare events. Many new methodologies have been proposed for observing these rare events. In this work, a temperature-aided cascade MD is proposed as a technique for studying the conformational transitions. Folding studies for Engrailed homeodomain and Immunoglobulin domain B of protein A have been carried out. Using this methodology, the unfolded structures with RMSD of 20 Å were folded to a structure with RMSD of 2 Å. Three sets of cascade MD runs were carried out using implicit solvation, explicit solvation, and charge updation scheme. In the charge updation scheme, charges based on the conformation obtained are calculated and are updated in the topology file. In all the simulations, the structure of 2 Å was reached within a few nanoseconds using these methods. Umbrella sampling has been performed using snapshots from the temperature-aided cascade MD simulation trajectory to build an entire conformational transition pathway. The advantage of the method is that the possible pathways for a particular reaction can be explored within a short duration of simulation time and the disadvantage is that the knowledge of the start and end state is required. The charge updation scheme adds the polarization effects in the force fields. This improves the electrostatic interaction among the atoms, which may help the protein to fold faster.

Probing the wild-type HRas activation mechanism using steered molecular dynamics, understanding the energy barrier and role of water in the activation.

Ras is one of the most common oncogenes in human cancers. It belongs to a family of GTPases that functions as binary conformational switches by timely switching of their conformations from GDP to GTP and vice versa. It attains the final active state structure via an intermediate GTP-bound state. The transition between these states is a millisecond-time-scale event. This makes studying this mechanism beyond the scope of classical molecular dynamics. In the present study, we describe the activation pathway of the HRas protein complex along the distance-based reaction coordinate using steered molecular dynamics. Approximately ~720 ns of MD simulations using CMD and SMD was performed. We demonstrated the change in orientation and arrangement of the two switch regions and the role of various hydrogen bonds during the activation process. The weighted histogram analysis method was also performed, and the potential of mean force was calculated between the inactive and active via the intermediate state (state 1) of HRas. The study indicates that water seems to play a crucial role in the activation process and to transfer the HRas protein from its intermediate state to the fully active state. The implications of our study hereby suggest that the HRas activation mechanism is a multistep process. It starts from the inactive state to an intermediate state 1 followed by trapping of water molecules and flipping of the Thr35 residue to form a fully active state (state 2). This state 2 also comprises Gly60, Thr35, GTP, Mg(2+) and water-forming stable interactions.

Structural insights into the activation mechanism of phosphoinositide 3-kinase alpha.

Phosphoinositide 3-kinases (PI3Ks) are lipid kinases known to regulate important cellular functions by phosphorylating the inositol ring of inositol-phospholipids (PtdIns) at 3' position. The PI3Kα is a heterodimer and the activation of the catalytic subunit (p110α) is regulated by its regulatory subunit (p85α). The current work deals with studying the activation mechanism of the PI3Kα using multi micro-second molecular dynamic simulations. Structural changes involved in activation mechanism is studied by gradually releasing the inhibitory effects of different domains of regulatory subunit namely, n-terminal SH2 (nSH2) and inter SH2 (iSH2). The observation shows that even in the presence of n-terminal and inter SH2 domain (niSH2) of regulatory subunit, the catalytic domain has some intrinsic activation activity and the presence of c-terminal SH2 (cSH2) domain may be required for complete inhibition. The release of nSH2 domain leads to loss of interactions between iSH2 domain (regulatory subunit) and C2 and kinase domain (catalytic subunit). The study shows that early events in the activation mechanism involve the movement of the ABD domain of the catalytic subunit along with the linker region between ABD and RBD region which may lead to movement of ABD closer to the CLobe of the kinase domain. This movement is essentially as it triggers the rearrangement of CLobe especially the catalytic loop and activation loop which bring catalytic important residues closer to ATP and PIP2(phosphatidylinositol-4,5-bisphosphate). Water mediated interaction analysis reveal that water may be playing an important role in the transfer of phosphate from ATP to PIP2. The study shows that initial signal for release of inhibitory effect of the regulatory subunit might be propagated through the linker region between ABD and RBD through allosteric effect to different regions of the protein. These understanding of early events during the activation mechanism may help in the design of better therapeutic targeting PI3K.

MolToxPred: small molecule toxicity prediction using machine learning approach.

Different types of chemicals and products may exhibit various health risks when administered into the human body. For toxicity reasons, the number of new drugs entering the market through the conventional drug development process has been reduced over the years. However, with the advent of big data and artificial intelligence, machine learning techniques have emerged as a potential solution for predicting toxicity and ensuring efficient drug development and chemical safety. An ML model for toxicity prediction can reduce experimental costs and time while addressing ethical concerns by drastically reducing the need for animals and clinical trials. Herein, MolToxPred, an ML-based tool, has been developed using a stacked model approach to predict the potential toxicity of small molecules and metabolites. The stacked model consists of random forest, multi-layer perceptron, and LightGBM as base classifiers and Logistic Regression as the meta classifier. For training and validation purposes, a comprehensive set of toxic and non-toxic molecules is curated. Different structural and physicochemical-based features in the form of molecular descriptors and fingerprints were employed. MolToxPred utilizes a comprehensive feature selection process and optimizes its hyperparameters through Bayesian optimization with stratified 5-fold cross-validation. In the evaluation phase, MolToxPred achieved an AUROC of 87.76% on the test set and 88.84% on an external validation set. The McNemar test was used as the post-hoc test to determine if the stacked models' performance was significantly different compared to the base learners. The developed stacked model outperformed its base classifiers and an existing tool in the literature, reaffirming its better performance. The hypothesis is that the incorporation of a diverse set of data, the subsequent feature selection, and a stacked ensemble approach give MolToxPred the edge over other methods. In addition to this, an attempt has been made to identify structural alerts responsible for endpoints of the Tox21 data to determine the association of a molecule with a plausible downstream pathway of action. MolToxPred may be helpful for drug discovery and regulatory pipelines in pharmaceutical and other industries for in silico toxicity prediction of small molecule candidates.

Studying early structural changes in SOS1 mediated KRAS activation mechanism.

KRAS activation is known to be modulated by a guanine nucleotide exchange factor (GEF), namely, Son of Sevenless1 (SOS1). SOS1 facilitates the exchange of GDP to GTP thereby leading to activation of KRAS. The binding of GDP/GTP to KRAS at the REM/allosteric site of SOS1 regulates the activation of KRAS at CDC25/catalytic site by facilitating its exchange. Different aspects of the allosteric activation of KRAS through SOS1 are still being explored. To understand the SOS1 mediated activation of KRAS, molecular dynamics simulations for a total of nine SOS1 complexes (KRAS-SOS1-KRAS) were performed. These nine systems comprised different combinations of KRAS-bound nucleotides (GTP/GDP) at REM and CDC25 sites of SOS1. Various conformational and thermodynamic parameters were analyzed for these simulation systems. MMPBSA free energy analysis revealed that binding at CDC25 site of SOS1 was significantly low for GDP-bound KRAS as compared to that of GTP-bound KRAS. It was observed that presence of either GDP/GTP bound KRAS at the REM site of SOS1 affected the activation related changes in the KRAS present at CDC25 site. The conformational changes at the catalytic site of SOS1 resulting from GDP/GTP-bound KRAS at the allosteric changes may hint at KRAS activation through different pathways (slow/fast/rare). The allosteric effect on activation of KRAS at CDC25 site may be due to conformations adopted by switch-I, switch-II, beta2 regions of KRAS at REM site. The effect of structural rearrangements occurring at allosteric KRAS may have led to increased interactions between SOS1 and KRAS at both the sites. The SOS1 residues involved in these important interactions with KRAS at the REM site were R694, S732 and K735. Whereas the ones interacting with KRAS at CDC25 site were S807, W809 and K814. This may suggest the crucial role of these residues in guiding the allosteric activation of KRAS at CDC25 site. The conformational shifts observed in the switch-I, switch-II and alpha3 regions of KRAS at CDC25 site may be attributed to be a part of allosteric activation. The binding affinities, interacting residues and conformational dynamics may provide an insight into development of inhibitors targeting the SOS1 mediated KRAS activation.

The multi-dimensional impact of captopril modification on human serum albumin.

Captopril is a thiol drug, widely used for the management of hypertension and cardiovascular diseases. Reactive thiols are found to covalently modify the cysteines of plasma proteins and affect their structure and function. Human serum albumin (HSA) is prone to undergo modification by various low molecular weight compounds, including drugs. Cysteine34 (Cys34) in HSA has a free thiol group with antioxidant properties, considered to be the most redox-sensitive amino acid in plasma. Through mass-spectrometric analysis, we demonstrate for the first time that captopril forms a disulfide adduct at Cys34 residue and increases the protease susceptibility of HSA to trypsin. As evidenced by our biophysical and electron microscopy studies, HSA undergoes structural alteration, aggregation and morphological changes when treated with different captopril concentrations. Molecular dynamics studies further revealed the regions of secondary structural changes in HSA due to disulfide adduct formation by captopril at Cys34. It also elucidated the residues involved in the noncovalent interactions with captopril. It is envisaged that structural change in HSA may influence the efficacy of drug delivery as well as its own biological function. These findings may thus provide significant insights into the field of pharmacology intriguing further investigation into the effects of long-term captopril treatment.

Evaluating therapeutic potential of AYUSH-64 constituents against omicron variant of SARS-CoV-2 using ensemble docking and simulations.

The COVID-19 pandemic in the later phase showed the presence of the B.1.1.529 variant of the SARS-CoV-2 designated as Omicron. AYUSH-64 a poly herbal drug developed by Central Council for Research in Ayurvedic Sciences (CCRAS) has been recommended by Ministry of Ayush in asymptomatic, mild to moderate COVID-19 patients. One of the earlier, in-silico study has shown the binding of the constituents of AYUSH-64 to the main protease (Mpro) of the SARS-CoV-2. This study enlisted four phytochemicals of AYUSH-64, which were found to have significant binding with the Mpro. In continuation to the same, the current study proposes to understand the binding of these four phytochemicals to main protease (Mpro) and receptor binding domain (RBD) of spike protein of the Omicron variant. An enhanced molecular docking methodology, namely, ensemble docking has been used to find the most efficiently binding phytochemical. Using molecular dynamics (MD) simulations and clustering approach it was observed that the Mpro and RBD Spike of Omicron variant of SARS-CoV-2 in complex with human ACE2 tends to attain 4 and 8 conformational respectively. Based on the docking studies, the best binding phytochemical of the AYUSH-64, akummicine N-oxide was selected for MD simulations. MD simulations of akummicine N-oxide bound to omicron variant of Mpro and RBD Spike-ACE complex was performed. The conformational, interaction and binding energy analysis suggested that the akummicine N-oxide binds well with Mpro and RBD Spike-ACE2 complex. The interaction between RBD Spike and ACE2 was observed to weaken in the presence of akummicine N-oxide. Hence, it can be inferred that, these phytochemicals from AYUSH-64 formulation may have the potential to act against the Omicron variant of SARS-CoV-2.

In-silico evaluation of Oroxylum indicum vent compounds in the plausible treatment and prevention of nasopharyngeal cancer.

BACKGROUND: Shyonaka (Oroxylum indicum Vent) is widely used in Ayurveda and in ethnomedical practice for the treatment of inflammation, pain, diarrhea, non-healing ulcers, and cancer. Owing to the high prevalence of Epstein-Barr virus (EBV) infection in Nasopharyngeal carcinoma (NPC) patients, simultaneous targeting of proteins involved in both EBV replication and NPC proliferation might help to manage the disease effectively. OBJECTIVES: This study is designed to identify potential dual targeting inhibitors from Oroxylum indicum having the potential to inhibit both EBV and NPC. This study also attempted quantitative analysis of Shyonaka Bark Decoction (SBD) to confirm the presence of Baicalein and Chrysin which are predominant marker compounds of Shyonaka. METHODOLOGY: The HPLC analysis of stem bark and root bark of Oroxylum indicum was done to estimate the presence of marker compounds Baicalein and Chrysalin. The in-silico analysis included ADMET analysis followed by molecular docking of known compounds from Oroxylum indicum (retrieved from IMPPAT database) onto the target proteins of EBV (BHRF1, NEC1, dUTPase, Uracil DNA glycosylase) and NPC (COX-2, EGFR, and MDM2) using DOCK6 tool. Further validations were done using the molecular dynamics simulations of top screened molecules onto the selected target proteins using AMBER20 package and their corresponding MMGBSA binding free-energy values were calculated. RESULTS: The molecular docking revealed that the key molecules from the plant, scutellarein 7-rutinoside (S7R), scutellarin (SCU) and 6-hydroxyluteolin, Baicalein and 5,7-Dihydroxy-2-phenyl-6-[3,4,5-trihydroxy-6-(hydroxymethyl)oxan-2-yl]oxychromen-4-one (57D) are effectively intervening with the target proteins of EBV, one of the key causative factors of NPC and the NPC specific targets which have the potential to reduce tumor size and other consequences of NPC. The molecular dynamics simulations of S7R, Baicalein and 57D, Baicalein with MDM-2 protein and dUTPase protein, respectively, showed stable interactions between them which were further assessed by the binding energy calculations. CONCLUSION: Overall, the in-silico evaluation of these phytochemicals with target proteins indicates their potential to inhibit both EBV and NPC which needs further in-vitro and in-vivo validations.

Probing the structure of human tRNA3Lys in the presence of ligands using docking, MD simulations and MSM analysis.

The tRNA3Lys, which acts as a primer for human immunodeficiency virus type 1 (HIV-1) reverse transcription, undergoes structural changes required for the formation of a primer-template complex. Small molecules have been targeted against tRNA3Lys to inhibit the primer-template complex formation. The present study aims to understand the kinetics of the conformational landscape spanned by tRNA3Lys in apo form using molecular dynamics simulations and Markov state modeling. The study is taken further to investigate the effect of small molecules like 1,4T and 1,5T on structural conformations and kinetics of tRNA3Lys, and comparative analysis is presented. Markov state modeling of tRNA3Lys apo resulted in three metastable states where the conformations have shown the non-canonical structures of the anticodon loop. Based on analyses of ligand-tRNA3Lys interactions, crucial ion and water mediated H-bonds and free energy calculations, it was observed that the 1,4-triazole more strongly binds to the tRNA3Lys compared to 1,5-triazole. However, the MSM analysis suggest that the 1,5-triazole binding to tRNA3Lys has brought rigidity not only in the binding pocket (TΨC arm, D-TΨC loop) but also in the whole structure of tRNA3Lys. This may affect the easy opening of primer tRNA3Lys required for HIV-1 reverse transcription.

Design of LNA Analogues Using a Combined Density Functional Theory and Molecular Dynamics Approach for RNA Therapeutics.

Antisense therapeutics treat a wide spectrum of diseases, many of which cannot be addressed with the current drug technologies. In the quest to design better antisense oligonucleotide drugs, we propose five novel LNA analogues (A1-A5) for modifying antisense oligonucleotides and establishing each with the five standard nucleic acids: adenine (A), guanine (G), cytosine (C), thymine (T), and uracil (U). Monomer nucleotides of these modifications were considered for a detailed Density Functional Theory (DFT)-based quantum chemical analysis to determine their molecular-level structural and electronic properties. A detailed MD simulation study was done on a 14-mer ASO (5'-CTTAGCACTGGCCT-3') containing these modifications targeting PTEN mRNA. Results from both molecular- and oligomer-level analysis clearly depicted LNA-level stability of the modifications, the ASO/RNA duplexes maintaining stable Watson-Crick base pairing preferring RNA-mimicking A-form duplexes. Notably, monomer MO isosurfaces for both purines and pyrimidines were majorly distributed on the nucleobase region in modifications A1 and A2 and in the bridging unit in modifications A3, A4, and A5, suggesting that A3/RNA, A4/RNA, and A5/RNA duplexes interact more with the RNase H and solvent environment. Accordingly, solvation of A3/RNA, A4/RNA, and A5/RNA duplexes was higher compared to that of LNA/RNA, A1/RNA, and A2/RNA duplexes. This study has resulted in a successful archetype for creating advantageous nucleic acid modifications tailored for particular needs, fulfilling a useful purpose of designing novel antisense modifications, which may overcome the drawbacks and improve the pharmacokinetics of existing LNA antisense modifications.

Novel alkyl-substituted 4-methoxy benzaldehyde thiosemicarbazones: Multi-target directed ligands for the treatment of Alzheimer's disease.

Alzheimer's disease (AD) is a devastating neurodegenerative disorder affecting mental ability and interrupts neurocognitive functions. Treating multifactorial conditions of AD with a single-target-directed drug is highly difficult. Thus, a multi-target-directed ligand (MTDL) development strategy has been developed as a promising approach for the treatment of AD. Herein, we have synthesized two novel thiosemicarbazones as MTDLs and reported their bioactivities against diverse neuropathological events involved in AD. In vitro studies revealed that both compounds exhibited promising anticholinesterase activity (AChE, IC50 = 15.98 µM, MZET and IC50 = 30.23 µM, MZMT), well supported by a detailed computational study. Both analogs have shown good thermodynamic behaviour and stability through interactions with characteristic amino acid residues throughout simulation of 100 ns against acetylcholinesterase enzyme. In an electrophysiology assay, these analogs have shown a characteristic inhibitory response against the GluN1-1a + GluN2B subunit of N-methyl-D-aspartate receptors. Pre-treatment of BV-2 microglial cells with MZET effectively decreased nitrite production compared to nitrite produced by lipopolysaccharide-treated cells alone. Further, the effect of MZMT and MZET on autophagy regulation was determined using stably transfected SH-SY5Y neuroblastoma cells. MZET significantly enhanced the autophagy flux in neuroblastoma cells. A significant decrease in copper-catalysed oxidation of amyloid-ß in presence of synthesized thiosemicarbazones was also observed. Collectively, our findings indicated that these analogs have potential as effective anti-AD candidates and can be used as a prototype to develop more safer multi-targeted anti-AD drugs.

Structural insight into the binding interactions of NTPs and nucleotide analogues to RNA dependent RNA polymerase of SARS-CoV-2.

RNA dependent RNA polymerase (RdRP) from positive-stranded RNA viruses has always been a hot target for designing of new drugs. Major class of drugs that are targeted against RdRP are nucleotide analogues. Extensive docking and molecular dynamics study describing the binding of natural nucleotides (NTPs) and its analogues leading to significant structural variation in the RdRP has been presented here. RdRP simulations in its apo, NTP-bound, and analogue-bound form have been performed. Nucleotide analogues included in this study were, favipiravir, galidesivir, lamivudine, ribavirin, remdesivir and sofosbuvir. The conformational flexibility of the RdRP molecule has been explored using principal component (PCA) and Markov state modeling (MSM) analysis. PCA inferred the presence of correlated motions among the conserved motifs of RdRP. Inter-domain distances between the finger and thumb subdomain flanking the nascent RNA template entry site sampled open and closed conformations. The ligand and template binding motifs F and G showed negatively correlated motions. K551, R553, and R555, a part of motif F appear to form strong interactions with the ligand molecules. R836, a primer binding residue was observed to strongly bind to the analogues. MSM analysis helped to extract statistically distinct conformations explored by the RdRP. Ensemble docking of the ligands on the Markov states also suggested the involvement of the above residues in ligand interactions. Markov states obtained clearly demarcated the open/closed conformations of the template entry site. These observations on residues from the conserved motifs involved in binding to the ligands may provide an insight into designing new inhibitors.Communicated by Ramaswamy H. Sarma.

Understanding the binding affinities between SFRP1CRD, SFRP1Netrin, Wnt5B and frizzled receptors 2, 3 and 7 using MD simulations.

cWnt-signalling plays a crucial role in stem cell maintenance and tissue homeostasis. Secreted frizzled-related proteins(SFRP), Wnt inhibitors consist of the N-terminal cysteine rich domain(CRD) and the C-terminal netrin(NTR) domain. SFRP1 binds to the Wnt ligands and frizzled receptors(FZ) either through its SFRP1CRD or through its SFRP1Netrin domains; however, very little is known on these binding affinities. Here, we attempted to understand the interactions and binding affinities of SFRP1-Wnt5B, SFRP1-FZ(2, 3 & 7) and Wnt5B-FZ(2, 3 & 7) that are mainly expressed in murine hair follicle stem cells. SFRP1CRD, SFRP1Netrin, Wnt5B and FZ(2, 3 & 7) structures were built using homology modelling, followed by their molecular dynamics simulations. SFRP1CRD showed lower fluctuation when in complex with FZ2, FZ3 and FZ7 and Wnt5B as compared to SFRP1Netrin using RMSF and RMSD. However, free energy showed SFRP1Netrin was energetically more stable than SFRP1CRD. SFRP1Netrin formed more number of interactions with FZ as compared to SFRP1CRD. Importantly, SFRP1Netrin favoured binding to the FZ receptors(FZ3 > FZ7 > FZ2) as compared to Wnt5B ligand. Conversely, the SFRP1CRD showed more affinity towards the Wnt5B ligand as compared to FZ receptors. Wnt5B showed the best binding affinity with FZ3 followed by SFRP1CRD and SFRP1Netrin. Therefore, SFRP1Netrin can bind to the FZ3 with higher binding affinity and may inhibit non-canonical Wnt-signalling pathway. Our study provides the comprehensive information on the binding affinities among the Wnt5B, SFRP1CRD/Netrin and FZ(2, 3 & 7). Thus, this information might also help in designing novel strategies to inhibit aberrant Wnt-signalling.Communicated by Ramaswamy H. Sarma.

Molecular dynamics of hERG channel: insights into understanding the binding of small molecules for detuning cardiotoxicity.

Evaluation of cardiotoxicity potential of new chemical entities (NCEs) has lately become one of the stringent filters in the drug discovery and development process. Cardiotoxicity is caused mainly by the inhibition of human ether-a-go-go related gene (hERG) channel protein. Inhibition of the hERG channel leads to a life-threatening condition known as cardiac arrhythmia. Knowledge of the structural behaviour of the hERG would aid greatly in the design of new drug molecules that do not interact with the protein and add to the safety index. In this study, a computational model for the active-state of hERG was developed. This model was equilibrated by performing the molecular dynamics simulations for 100 ns followed by clustering and selection of a representative structure based on the largest populated cluster. To study the changes in the protein structure on inhibition, three inhibitory ligands, namely, dofetilide, cisapride and terfenadine were docked, followed by molecular dynamics simulations of 200 ns for the apo and each ligand-bound structure. It was observed that docking and simulation studies of the hERG model exhibited noticeable conformational changes in the protein upon ligand-binding. A significant change in the kink of the S6-transmembrane helix was observed. Inter-chain distances between the crucial residues Y652 and F656 (present below the ion-selectivity filter), their side-chain orientation and hydrogen bonding indicated a probable collapse of the pore. These changes may infer the initiation in transition of hERG from an open to an inactive state. Hence, these findings would help in designing compounds devoid of hERG inhibition with reduced cardiotoxicity.Communicated by Ramaswamy H. Sarma.

Destabilization potential of beta sheet breaker peptides on Abeta fibril structure: an insight from molecular dynamics simulation study.

Alzheimer's disease is characterized by amyloid-ß aggregation. Currently, all the approved medications are to treat the symptoms but there is no clinically approved treatment for the cure or to prevent the progression of Alzheimer's disease (AD). Earlier reports suggest the use of small molecules and peptides to target and destabilize the amyloid fibril. The use of Beta Sheet Breaker (BSB) peptides seems to be a promising and attractive therapeutic approach as it can strongly bind and destabilize the preformed amyloid fibril. There are experimental studies describing the destabilization role of various BSB peptides, but the exact mechanism remains elusive. In the current work, an attempt is made to study the destabilization mechanism of different BSB peptides on preformed amyloid protofibril using molecular docking and simulations. Molecular docking of eight different BSB peptides of varying length (5-mer to 10-mer) on the Abeta protofibril was done. Docking was followed by multiple sets of molecular simulations for the Abeta protofibril-BSB peptide complex for each of the top ranked poses of the eight BSB peptides. As a control, multiple sets of simulations for the Abeta protofibril (APO) were also carried out. An increase in the RMSD, decrease in the number of interchain hydrogen bonds, destabilization of important salt bridge interactions (D23-K28), and destabilization of interchain hydrophobic interactions suggested the destabilization of Abeta protofibril by BSB peptides. The MM-GBSA free energy of binding for each of the BSB peptides was calculated to measure the binding affinity of BSB peptides to Abeta protofibril. Further residue wise contribution of free energy of binding was also calculated. The study showed that 7-mer peptides tend to bind strongly to Abeta protofibril as compared to other BSB peptides. The KKLVFFA peptide showed better destabilization potential as compared to the other BSB peptides. The details about the destabilization mechanism of BSB peptides will help in the design of other peptides for the therapeutic intervention for AD.