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.

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.

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.

Drug repurposing studies targeting SARS-CoV-2: an ensemble docking approach on drug target 3C-like protease (3CLpro).

The COVID-19 pandemic has been responsible for several deaths worldwide. The causative agent behind this disease is the Severe Acute Respiratory Syndrome - novel Coronavirus 2 (SARS-CoV-2). SARS-CoV-2 belongs to the category of RNA viruses. The main protease, responsible for the cleavage of the viral polyprotein is considered as one of the hot targets for treating COVID-19. Earlier reports suggest the use of HIV anti-viral drugs for targeting the main protease of SARS-CoV, which caused SARS in the year 2002-2003. Hence, drug repurposing approach may prove to be useful in targeting the main protease of SARS-CoV-2. The high-resolution crystal structure of the main protease of SARS-CoV-2 (PDB ID: 6LU7) was used as the target. The Food and Drug Administration approved and SWEETLEAD database of drug molecules were screened. The apo form of the main protease was simulated for a cumulative of 150 ns and 10 µs open-source simulation data was used, to obtain conformations for ensemble docking. The representative structures for docking were selected using RMSD-based clustering and Markov State Modeling analysis. This ensemble docking approach for the main protease helped in exploring the conformational variation in the drug-binding site of the main protease leading to the efficient binding of more relevant drug molecules. The drugs obtained as top hits from the ensemble docking possessed anti-bacterial and anti-viral properties. This in silico ensemble docking approach would support the identification of potential candidates for repurposing against COVID-19.Communicated by Ramaswamy H. Sarma.

Investigating DNA Binding and Conformational Variation in Temperature Sensitive p53 Cancer Mutants Using QM-MM Simulations.

The tp53 gene is found to be mutated in 50% of all the cancers. The p53 protein, a product of tp53 gene, is a multi-domain protein. It consists of a core DNA binding domain (DBD) which is responsible for its binding and transcription of downstream target genes. The mutations in p53 protein are responsible for creating cancerous conditions and are found to be occurring at a high frequency in the DBD region of p53. Some of these mutations are also known to be temperature sensitive (ts) in nature. They are known to exhibit partial or strong binding with DNA in the temperature range (298-306 K). Whereas, at 310 K and above they show complete loss in binding. We have analyzed the changes in binding and conformational behavior at 300 K and 310 K for three of the ts-mutants viz., V143A, R249S and R175H. QM-MM simulations have been performed on the wild type and the above mentioned ts-mutants for 30 ns each. The optimal estimate of free energy of binding for a particular number of interface hydrogen bonds was calculated using the maximum likelihood method as described by Chodera et. al (2007). This parameter has been observed to be able to mimic the binding affinity of the p53 ts-mutants at 300 K and 310 K. Thus the correlation between MM-GBSA free energy of binding and hydrogen bonds formed by the interface residues between p53 and DNA has revealed the temperature dependent nature of these mutants. The role of main chain dihedrals was obtained by performing dihedral principal component analysis (PCA). This analysis, suggests that the conformational variations in the main chain dihedrals (ϕ and ψ) of the p53 ts-mutants may have caused reduction in the overall stability of the protein. The solvent exposure of the side chains of the interface residues were found to hamper the binding of the p53 to the DNA. Solvent Accessible Surface Area (SASA) also proved to be a crucial property in distinguishing the conformers obtained at 300 K and 310 K for the three ts-mutants from the wild type at 300 K.

QM-MM simulations on p53-DNA complex: a study of hot spot and rescue mutants.

p53 is a transcription factor involved in the expression of a number of downstream genes in response to genotoxic stress. It is activated through post translation modifications in normal as well as cancerous cells. However, due to mutations occurring in p53 in cancer cells it is not able to perform its function of DNA binding which leads to cell proliferation. It is found to be mutated in 50% of the cancers. These mutations occur at a high frequency in the DNA binding region of the p53. Among the known seven hot spot cancer mutations G245S, R249S, and R273C have been studied here using quantum mechanics and molecular mechanics (QM-MM) simulations. These mutations along with their experimentally proven rescue mutations have also been included in the present work. A comparative study of these cancer mutations along with wild type and their rescue mutations has been performed. A computational measure based on the free energy changes occurring in the binding of the p53 to the DNA has been presented. A correlation between the DNA binding property and important interaction between p53 and DNA has been observed for all the mutants. The keys residues which contribute to the binding of p53 to DNA by forming crucial hydrogen bonds have also been discussed in detail. A 30 ns simulation study was analyzed to observe the local structural changes and DNA binding property of p53 in case of wild type, cancer and rescue mutants.