In silico and structure-based evaluation of deleterious mutations identified in human Chk1, Chk2, and Wee1 protein kinase.

Checkpoint kinases Chk1, Chk2, Wee1 are playing a key role in DNA damage response and genomic integrity. Cancer-associated mutations identified in human Chk1, Chk2, and Wee1 were retrieved to understand the function associated with the mutation and also alterations in the folding pattern. Therefore, an attempt has been made to identify deleterious effect of variants using in silico and structure-based approach. Variants of uncertain significance for Chk1, Chk2, and Wee1 were retrieved from different databases and four prediction servers were employed to predict pathogenicity of mutations. Further, Interpro, I-Mutant 3.0, Consurf, TM-align, and have (y)our protein explained were used for comprehensive study of the deleterious effects of variants. The sequences of Chk1, Chk2, and Wee1 were analyzed using Clustal Omega, and the three-dimensional structures of the proteins were aligned using TM-align. The molecular dynamics simulations were performed to explore the differences in folding pattern between Chk1, Chk2, Wee1 wild-type, and mutant protein and also to evaluate the structural integrity. Thirty-six variants in Chk1, 250 Variants in Chk2, and 29 in Wee1 were categorized as pathogenic using in silico prediction tools. Furthermore, 25 mutations in Chk1, 189 in Chk2, and 14 in Wee1 were highly conserved, possessing deleterious effect and also influencing the protein structure and function. These identified mutations may provide underlying genetic intricacies to serve as potential targets for therapeutic inventions and clinical management.

Phytocompounds from Phyllanthus acidus (L.) Skeels in the management of Monkeypox Virus infections.

Monkeypox is a communicable disease similar to smallpox, primarily occurring in African countries. However, recently it has spread to countries outside Africa and may arise as the next threat after COVID-pandemic. The causative organism, i.e. Monkeypox Virus (MPV) spreads from one individual to another primarily through inhalation of respiratory droplets or through contact with skin lesions of infected individuals. No known drugs are available specifically for MPV. Due to its similarity with smallpox, treatment of monkeypox is being attempted through the administration of the smallpox vaccine. Therefore, we evaluated the efficacy of the plant Phyllanthus acidus against MPV since it is traditionally used in the treatment of chickenpox and smallpox. Through functional annotation, PASS prediction and Network pharmacology analysis, the effectiveness of these chosen P. acidus-derived phytocompounds against MPV was confirmed. Target prediction of the phytocompounds identified in GC-MS analysis of the plant extract showed them to be associated with 76 human proteins. The compounds also show good binding affinity with selected viral proteins: DNA polymerase (DNApol), Putative Virulence Factor (vPVF) and Cytokine Binding Protein. Prediction of Activity Spectra for Substances (PASS) and functional annotation of the target proteins further support their antiviral nature through interaction with these proteins. The compounds were found to modulate pathways related to symptoms of viral infection and this may help in maintaining homeostasis. Our study demonstrates antiviral activity as well as the therapeutic potential of the plant against MPV infection.Communicated by Ramaswamy H. Sarma.

Evaluation and analysis of novel germline variants in ethanol metabolism pathway genes predisposition to liver disease.

The pathogenetic events of liver disease are seemingly determined by factors linked to ethanol metabolism. The variations in genes encoding enzymes of the ethanol metabolic pathway can influence exposure to alcohol and thus may act as risk factors for the development of liver disease. The present study aimed to understand the genetic aspect of germline variations in ethanol metabolic pathway genes in two major categories of liver disease i.e. ALD and NAFLD. Targeted Re-sequencing was performed in the two disease categories along with healthy control followed by an assessment and evaluation of the variants in a case vs control manner. The pathogenicity prediction was evaluated using SIFT, PolyPhen, PROVEN, LRT, CADD, FATHMM, EIGEN, REVEL and VarSome, while MD simulation of a novel significant variant was performed using the GROMACS 5.1.4 package. The annotation of targeted re-sequencing results revealed 2172 variants in different locations of the genes. Upon recurrent assessment predominantly focusing on exonic missense variants from these genes of the alcohol metabolism pathway, the ALDH1L2 [c.337C > G, p.Pro113Ala, (rs199841702)] variant was found highly significant with comprehensive results. The amino acid substitution tool that predicted protein stability due to a point mutation showed a decrease in stability. The genotyping distribution of the identified novel variant in the population revealed that heterozygosity is significantly distributed in ALD patients. However, the predominant association between the inherited variant and the cause of developing disease needs further robust study.

Evaluation of Pathogenicity and Structural Alterations for the Mutations Identified in the Conserved Region of the C-Terminal Kinase Domain of Human-Ribosomal S6 Kinase 1.

Human-ribosomal s6 kinase 1 (h-RSK1) is an effector kinase of the Ras/MAPK signaling pathway, which is involved in the regulation of the cell cycle, proliferation, and survival. RSKs comprise two functionally distinct kinase domains at the N-terminal (NTKD) and C-terminal (CTKD) separated by a linker region. The mutations in RSK1 may have the potential to provide an extra benefit to the cancer cell to proliferate, migrate, and survive. The present study focuses on evaluating the structural basis for the missense mutations identified at the C-terminal kinase domain of human-RSK1. A total of 139 mutations reported on RSK1 were retrieved from cBioPortal, where 62 were located at the CTKD region. Furthermore, 10 missense mutations Arg434Pro, Thr701Met, Ala704Thr, Arg725Trp, Arg726Gln, His533Asn, Pro613Leu, Ser720Cys, Arg725Gln, and Ser732Phe were predicted to be deleterious using in silico tools. To our observation, these mutations are located in the evolutionarily conserved region of RSK1 and shown to alter the inter- and intramolecular interactions and also the conformational stability of RSK1-CTKD. The molecular dynamics (MD) simulation study further revealed that the five mutations Arg434Pro, Thr701Met, Ala704Thr, Arg725Trp, and Arg726Gln showed maximum structural alterations in RSK1-CTKD. Thus, based on the in silico and MD simulation analysis, it can be concluded that the reported mutations may serve as potential candidates for further functional studies.

Deciphering the role of aquaporin 1 in the adaptation of the stinging catfish Heteropneustes fossilis to environmental hypertonicity by molecular dynamics simulation studies.

A thorough investigation of the water permeability of H. fossilis aquaporin 1 (hfAQP1) in a hypertonic environment can provide a useful insight into the understanding of the underlying molecular mechanism of its high tolerance to salinity. Here, we constructed a 3 D homology model of hfAQP1 by taking Bos taurus AQP1, AQP0, and human AQP2 as templates using I-TASSER. The model obtained has similar structural organizations with mammalian AQP1s in all aspects. We investigated the water permeability of the modeled hfAQP1 in a 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC) membrane under neutral and 100 mM hypersalinity by subjecting each system to a 100 ns molecular dynamics simulation. Our results show that hypersalinity hinders water permeation across the membrane through the hfAQP1 channel. A change in the intermolecular distance between key residues of the ar/R selectivity filter along with charge redistribution resulted in the accommodation of only 2-6 water molecules inside the channel at once under hypersaline conditions. We investigated the mRNA expression pattern of hfaqp1 in osmoregulatory organs of H. fossilis in response to 100 mM hypertonicity by using qPCR analysis. The transcript was downregulated in kidney and GI tract, but upregulated in the Gills. Thus, the catfish survive in a hypertonic environment by reducing the transport of water in its cellular systems and downregulating the expression of the hfaqp1 gene. The results observed in our study can shed more light on the functionality of AQP1 in catfishes under salinity stress and aid in future researches on solving more gating mechanisms involved in its regulation.Communicated by Ramaswamy H. Sarma.

Computational approaches for evaluation of isobavachin as potential inhibitor against t877a and w741l mutations in prostate cancer.

Prostate cancer is the World's second most common cancer, with the fifth-highest male mortality rate. Point mutations such as T877A and W741L are frequently seen in advanced prostate cancer patients, conferring drug-resistance and hence driving cancer growth. Such occurrence of drug resistance in prostate cancer necessitates designing of suitable ligands to ensure better interactions with the receptors which can block the progression of the disease. The present study focus on the modification of plant-derived flavonoids that might act as inhibitors against such point mutations namely, T877A and W741L. In T877A mutation threonine is substituted by alanine at the 877 codon and W741L mutation, tryptophan is substituted by lysine at the 741 codon in prostate cancer. The study revolved on the aspect of the evaluation of Isobavachin and its derivatives as a potential agent to tackle such point mutations by using the in silico approach. A total of 98 molecular dockings were performed to find the ligand-receptor complexes with the lowest binding energy employing Autodock Software to conduct the blind and site-specific docking. Additionally, ligands were screened for Drug-likeness and toxicity using several tools yielding eight possible drug candidates. Based on the results of Molecular Docking, Drug-likeness, and ADMET testing, ten structures, including six complexes and three receptors were subjected to molecular dynamics simulation of 100 ns covering RMSD, RMSF, Rg, and MM/PBSA. Based on the simulation results, Isobavachin, IsoMod4, and IsoMod7 were concluded to be stable and exhibited potential properties for developing a novel drug to combat prostate cancer and its associated drug-resistance.Communicated by Ramaswamy H. Sarma.

In Silico and Structure-Based Assessment of Similar Variants Discovered in Tandem Repeats of BRCT Domains of BRCA1 and BARD1 To Characterize the Folding Pattern.

BRCA1 and BARD1 are important proteins in the homologous DNA damage repair pathways. Different genetic variants identified in these proteins have been clinically correlated with the occurrence of hereditary breast and ovarian cancer (HBOC). Variants of unknown significance (VUS) reported in the BRCT domains of BRCA1 and BARD1 substantiate the importance of BRCT domain-containing proteins for genomic integrity. To classify the pathogenicity of variants, in silico, structural and molecular dynamics (MD)-based approaches were explored. Different variants reported in the BRCT region were retrieved from cBioPortal, LOVD3, BRCA Exchange, and COSMIC databases to evaluate the pathogenicity. Multiple sequence alignment and superimposition of the structures of BRCA1 BRCT and BARD1 BRCT domains were performed to compare alterations in folding patterns. From 11 in silico predictions servers, variants reported to be pathogenic by 70% of the servers were considered for structural analysis. To our observations, four residue pairs of both the proteins were reported, harboring 11 variants, H1686Y, W1718L, P1749L, P1749S, and W1837L variants for BRCA1 BRCT and H606D, H606N, W635L, P657L, P657S, and W762F for BARD1 BRCT. MD simulations of the BRCT repeat regions of these variants and wild-type proteins were performed to evaluate the differences of folding patterns. Root mean square deviation (RMSD), R g, solvent-accessible surface area (SASA), and root mean square fluctuation (RMSF) of variants showed slight differences in the folding patterns from the wild-type proteins. Furthermore, principal components analysis of H1686Y, P1749S, and W1718L variants of BRCA1 showed less flexibility than the wild type, whereas that of H606D, W635L, and W762F of BARD1 showed more flexibility than the wild type. Normal mode analysis of the energy minima from the simulation trajectories revealed that most of the variants do not show much differences in the flexibility compared to the wild-type proteins, except for the discrete regions in the BRCT repeats, most prominently in the 1798-1801 amino acid region of BRCA1 and at the residue 744 in BARD1.

Multimodal approach to characterize the tetrameric form of human PML-RBCC domain and ATO-mediated conformational changes.

RING-B box-coiled coil (RBCC) domain of promyelocytic leukemia (PML) comprises a zinc finger motif that is targeted by arsenic trioxide (ATO) to treat acute promyelocytic leukemia (APL) pathogenesis. Preliminary evidence suggests that the PML-RBCC has different functional characteristics, but no structural details have been reported despite its importance in differential expression and cell-cycle regulation. Therefore, the recombinant h-PML-RBCC protein was purified to its homogeneity, and characterized for oligomeric behaviour which indicated that RBCC domain exists as a tetramer in solution. Furthermore, nano-DSF and circular-dichroism demonstrated that the tetrameric form preserves its native conformation along with thermal stability (Tm = 83.2 °C). In-silico-based PML-RBCC structure was used to perform the molecular dynamics simulation for 300 ns in the presence of zinc atoms, which demonstrated the differential dynamic of PML-RBCC tetrameric chains. MMPBSA analysis also indicated the role of hydrophobic interactions that favours stable tetrameric structure of PML-RBCC. ATO-induced secondary and tertiary structure changes were observed in PML-RBCC using circular dichroism and fluorescence spectroscopy. Dynamic light scattering and transmission electron microscopy revealed ATO-induced higher-order oligomerization and aggregation of PML-RBCC. The unique oligomeric nature of the h-PML-RBCC protein and its interactions with ATO will help to understand the mechanism of APL pathogenesis and drug resistance.

Understanding the binding mechanism for potential inhibition of SARS-CoV-2 Mpro and exploring the modes of ACE2 inhibition by hydroxychloroquine.

As per the World Health Organization report, around 226 844 344 confirmed positive cases and 4 666 334 deaths are reported till September 17, 2021 due to the recent viral outbreak. A novel coronavirus (severe acute respiratory syndrome coronavirus 2 [SARS-CoV-2]) is responsible for the associated coronavirus disease (COVID-19), which causes serious or even fatal respiratory tract infection and yet no approved therapeutics or effective treatment is currently available to combat the outbreak. Due to the emergency, the drug repurposing approach is being explored for COVID-19. In this study, we attempt to understand the potential mechanism and also the effect of the approved antiviral drugs against the SARS-CoV-2 main protease (Mpro). To understand the mechanism of inhibition of the malaria drug hydroxychloroquine (HCQ) against SARS-CoV-2, we performed molecular interaction studies. The studies revealed that HCQ docked at the active site of the Human ACE2 receptor as a possible way of inhibition. Our in silico analysis revealed that the three drugs Lopinavir, Ritonavir, and Remdesivir showed interaction with the active site residues of Mpro. During molecular dynamics simulation, based on the binding free energy contributions, Lopinavir showed better results than Ritonavir and Remdesivir.

Identification of potential plant-based inhibitor against viral proteases of SARS-CoV-2 through molecular docking, MM-PBSA binding energy calculations and molecular dynamics simulation.

The Coronavirus disease 2019 (COVID-19), caused by the novel coronavirus, SARS-CoV-2, has recently emerged as a pandemic. Here, an attempt has been made through in-silico high throughput screening to explore the antiviral compounds from traditionally used plants for antiviral treatments in India namely, Tea, Neem and Turmeric, as potential inhibitors of two widely studied viral proteases, main protease (Mpro) and papain-like protease (PLpro) of the SARS-CoV-2. Molecular docking study using BIOVIA Discovery Studio 2018 revealed, (-)-epicatechin-3-O-gallate (ECG), a tea polyphenol has a binding affinity toward both the selected receptors, with the lowest CDocker energy - 46.22 kcal mol-1 for SARS-CoV-2 Mpro and CDocker energy - 44.72 kcal mol-1 for SARS-CoV-2 PLpro, respectively. The SARS-CoV-2 Mpro complexed with (-)-epicatechin-3-O-gallate, which had shown the best binding affinity was subjected to molecular dynamics simulations to validate its binding affinity, during which, the root-mean-square-deviation values of SARS-CoV-2 Mpro-Co-crystal ligand (N3) and SARS-CoV-2 Mpro- (-)-epicatechin-3-O-gallate systems were found to be more stable than SARS-CoV-2 Mpro system. Further, (-)-epicatechin-3-O-gallate was subjected to QSAR analysis which predicted IC50 of 0.3281 nM against SARS-CoV-2 Mpro. Overall, (-)-epicatechin-3-O-gallate showed a potential binding affinity with SARS-CoV-2 Mpro and could be proposed as a potential natural compound for COVID-19 treatment.

Potential anti-viral activity of approved repurposed drug against main protease of SARS-CoV-2: an in silico based approach.

The Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) which was first reported in Wuhan province of China, has become a deadly pandemic causing alarmingly high morbidity and mortality. In the absence of new targeted drugs and vaccines against SARS-CoV-2 at present, the choices for effective treatments are limited. Therefore, considering the exigency of the situation, we focused on identifying the available approved drugs as potential inhibitor against the promising Coronavirus drug target, the Main Protease, using computer-aided methods. We created a library of U. S. Food and Drug Administration approved anti-microbial drugs and virtually screened it against the available crystal structures of Main Protease of the virus. The study revealed that Viomycin showed the highest -CDocker energy after docking at the active site of SARS-CoV-2 Main Protease. It is noteworthy that Viomycin showed higher -CDocker energy as compared to the drugs currently under clinical trial for SARS-CoV-2 treatment viz. Ritonavir and Lopinavir. Additionally, Viomycin formed higher number of H-bonds with SARS-CoV-2 Main Protease than its co-crystallised inhibitor compound N3. Molecular dynamics simulation further showed that Viomycin embedded deeply inside the binding pocket and formed robust binding with SARS-CoV-2 Main Protease. Therefore, we propose that Viomycin may act as a potential inhibitor of the Main Protease of SARS-CoV-2. Further optimisations with the drug may support the much-needed rapid response to mitigate the pandemic.Communicated by Ramaswamy H. Sarma.

Composition and in silico structural analysis of fibroin from liquid silk of non-mulberry silkworm Antheraea assamensis.

Silk is spun from the liquid precursor known as liquid silk secreted from the posterior part and stored in the silk gland lumen with occurrence of many momentary events. The liquid silk in the silk gland is transformed to the spun silk fibre. In this study the elucidation of the protein components of liquid silk from the posterior part of the silk gland (PSG) of saturniid silkworm Antheraea assamensis along with its structural characterization has been reported. The 3D model of the N-terminal amorphous portion with some repeat crystalline motifs (19-255) of core protein fibroin has also been constructed. 1D and 2D electrophoresis revealed the homo-dimeric structure of the silk protein. Secondary structure analysis by Circular dichroism, FTIR spectroscopy showed α helical structural component as predominant conformation in the liquid silk. The crystalline structure investigated through X ray diffraction (XRD) analysis also revealed the presence of less ordered amorphous α helical conformation in the liquid silk. The 3D structural model proposed of the residues from 19 to 255 has revealed structural stability throughout the molecular dynamics simulation process. This study will provide the detailed structural information and in silico analysis of the core protein present in the liquid silk of PSG.

Molecular dynamics approach to probe the antigenicity of PagN - an outer membrane protein of Salmonella Typhi.

PagN is a highly immunogenic 27-kDa outer membrane adhesin present in Salmonella Typhi. It plays a major role in the pathogenesis of typhoid fever and has emerged as a strong vaccine candidate. In this report, we predict the three-dimensional structure of PagN and describe the conformational dynamics associated with its four extracellular loops based on two 100-ns molecular dynamics simulations at 300 and 310 K. The formation and deformation of the secondary structures on these loops were also investigated during the simulations which revealed loops L1 and L2 to be highly flexible, whereas the relative flexibility of loops L3 and L4 was minimal. Essential dynamics and principal component analysis deciphered more realistic dynamic behaviours of the loops, particularly at 310 K. Moreover, our epitope predictions suggest that the antigenic peptides for B-cell recognition are located within the loops L1 and L2, while those for T-cell recognition are located within the loops L3 and L4. The binding specificities of the antigenic peptides towards specific human MHC-I and MHC-II HLA alleles closely resembled the stability of the loops L3 and L4 inferred from the simulations. Finally, we identified potential antigenic peptides in the flexible (L1 and L2) as well as stable (L3 and L4) regions of PagN for both B- and T-cell recognitions, which can help in developing effective sub-unit vaccines.

Understanding the molecular basis of stability in Kunitz (STI) family of inhibitors in terms of a conserved core tryptophan residue: A theoretical investigation.

ß-trefoil is one of the superfolds among proteins. Important classes of proteins like Interleukins (ILs), FibroblastGrowth Factors (FGFs), Kunitz (STI) family of inhibitors etc. belong to this fold. Kunitz (STI) family of inhibitors of proteins possess a highly conserved and structurally important Trytophan 91 (W91) residue, which stitches the top layer of the barrel with the lid. In this article we have investigated the molecular insights of the involvement of this W91 residue in the stability and folding pathway of Kunitz (STI) family. Winged bean Chymotrypsin inhibitor (WCI), a member of Kunitz (STI) family was chosen as a model system for carrying out the work. Molecular dynamics (MD) simulations were run with a set of total six proteins, including wild type WCI (WT) & five mutants namely W91F, W91M, W91A, W91H and W91I. Among all of them the coordinates of four proteins were taken from their crystal structures deposited in the Protein Data Bank (PDB), where as the coordinates for the rest two was generated using in-silico modelling. Our results suggest that truly this W91 residue plays a determining role in stability and folding pathway of Kunitz (STI) family. The mutants are less stable and more susceptible to quicker unfolding at higher temperatures compared to the wild type WCI. These effects are most pronounced for the smallest mutants namely W91H and W91A, indicating more is the cavity created by mutation at W91 position more the proteins becomes unstable.

Molecular dynamics approach to probe the allosteric inhibition of PTP1B by chlorogenic and cichoric acid.

Protein tyrosine phosphatase 1B (PTP1B), a major negative regulator of the insulin and leptin signaling pathway, is a potential target for therapeutic intervention against diabetes and obesity. The recent discovery of an allosteric site in PTP1B has created an alternate strategy in the development of PTP1B targeted therapy. The current study investigates the molecular interactions between the allosteric site of PTP1B with two caffeoyl derivatives, chlorogenic acid (CGA) and cichoric acid (CHA), using computational strategies. Molecular docking analysis with CGA and CHA at the allosteric site of PTP1B were performed and the resulting protein-ligand complexes used for molecular dynamics simulation studies for a time scale of 10 ns. Results show stable binding of CGA and CHA at the allosteric site of PTP1B. The flexibility of the WPD loop was observed to be constrained by CGA and CHA in the open (inactive), providing molecular mechanism of allosteric inhibition. The allosteric inhibition of CGA and CHA of PTP1B was shown to be favorable due to no restriction by the α-7 helix in the binding of CGA and CHA at the allosteric binding site. In conclusion, our results exhibit an inhibitory pattern of CGA and CHA against PTP1B through potent binding at the allosteric site.

Tackling multiple antibiotic resistance in enteropathogenic Escherichia coli (EPEC) clinical isolates: a diarylheptanoid from Alpinia officinarum shows promising antibacterial and immunomodulatory activity against EPEC and its lipopolysaccharide-induced inflammation.

Antibiotic treatment for infectious diseases commonly leads to host inflammatory responses. Molecules with bifunctional antibacterial and anti-inflammatory properties could provide a solution for such clinical manifestations. Here we report such bifunctional activity for a diarylheptanoid (5-hydroxy-7-(4''-hydroxy-3-methoxyphenyl)-1-phenyl-3-heptanone) isolated from Alpinia officinarum, a medicinal plant belonging to the Zingiberaceae family, against enteropathogenic Escherichia coli (EPEC). The diarylheptanoid showed inhibitory and bactericidal activity against EPEC clinical isolates and efficiently suppressed EPEC lipopolysaccharide-induced inflammation in human peripheral blood mononuclear cells. In silico docking analysis revealed that the diarylheptanoid could interact with subunit A of E. coli DNA gyrase. Such molecules with bifunctional activity may be potential therapeutics for infectious diseases.

A modeled structure for amidase-03 from Bacillus anthracis.

Homology models of amidase-03 from Bacillus anthracis were constructed using Modeller (9v2). Modeller constructs protein models using an automated approach for comparative protein structure modeling by the satisfaction of spatial restraints. A template structure of Listeria monocytogenes bacteriophage PSA endolysin PlyPSA (PDB ID: 1XOV) was selected from protein databank (PDB) using BLASTp with BLOSUM62 sequence alignment scoring matrix. We generated five models using the Modeller default routine in which initial coordinates are randomized and evaluated by pseudo-energy parameters. The protein models were validated using PROCHECK and energy minimized using the steepest descent method in GROMACS 3.2 (flexible SPC water model in cubic box of size 1 Å instead of rigid SPC model). We used G43a1 force field in GROMACS for energy calculations and the generated structure was subsequently analyzed using the VMD software for stereo-chemistry, atomic clash and misfolding. A detailed analysis of the amidase-03 model structure from Bacillus anthracis will provide insight to the molecular design of suitable inhibitors as drug candidates.