Mutagenesis of DsbAss is Crucial for the Signal Recognition Particle Mechanism in Escherichia coli: Insights from Molecular Dynamics Simulations

Structural insights into the mechanism of resistance to bicalutamide by the clinical mutations in androgen receptor in chemo-treatment resistant prostate cancer

Despite advanced diagnosis and detection technologies, prostate cancer (PCa) is the most prevalent neoplasms in males. Dysregulation of the androgen receptor (AR) is centrally involved in the tumorigenesis of PCa cells. Acquisition of drug resistance due to modifications in AR leads to therapeutic failure and relapse in PCa. An overhaul of comprehensive catalogues of cancer-causing mutations and their juxta positioning on 3D protein can help in guiding the exploration of small drug molecules. Among several well-studied PCa-specific mutations, T877A, T877S and H874Y are the most common substitutions in the ligand-binding domain (LBD) of the AR. In this study, we combined structure as well as dynamics-based in silico approaches to infer the mechanistic effect of amino acid substitutions on the structural stability of LBD. Molecular dynamics simulations allowed us to unveil a possible drug resistance mechanism that acts through structural alteration and changes in the molecular motions of LBD. Our findings suggest that the resistance to bicalutamide is partially due to increased flexibility in the H12 helix, which disturbs the compactness, thereby reducing the affinity for bicalutamide. In conclusion, the current study helps in understanding the structural changes caused by mutations and could assist in the drug development process.Communicated by Ramaswamy H. Sarma.

Structural and functional implications of SLC13A3 and SLC9A6 mutations: an in silico approach to understanding intellectual disability

BACKGROUND

Intellectual disability (ID) is a condition that varies widely in both its clinical presentation and its genetic underpinnings. It significantly impacts patients' learning capacities and lowers their IQ below 70. The solute carrier (SLC) family is the most abundant class of transmembrane transporters and is responsible for the translocation of various substances across cell membranes, including nutrients, ions, metabolites, and medicines. The SLC13A3 gene encodes a plasma membrane-localized Na+/dicarboxylate cotransporter 3 (NaDC3) primarily expressed in the kidney, astrocytes, and the choroid plexus. In addition to three Na + ions, it brings four to six carbon dicarboxylates into the cytosol. Recently, it was discovered that patients with acute reversible leukoencephalopathy and a-ketoglutarate accumulation (ARLIAK) carry pathogenic mutations in the SLC13A3 gene, and the X-linked neurodevelopmental condition Christianson Syndrome is caused by mutations in the SLC9A6 gene, which encodes the recycling endosomal alkali cation/proton exchanger NHE6, also called sodium-hydrogen exchanger-6. As a result, there are severe impairments in the patient's mental capacity, physical skills, and adaptive behavior.

METHODS AND RESULTS

Two Pakistani families (A and B) with autosomal recessive and X-linked intellectual disorders were clinically evaluated, and two novel disease-causing variants in the SLC13A3 gene (NM 022829.5) and the SLC9A6 gene (NM 001042537.2) were identified using whole exome sequencing. Family-A segregated a novel homozygous missense variant (c.1478 C > T; p. Pro493Leu) in the exon-11 of the SLC13A3 gene. At the same time, family-B segregated a novel missense variant (c.1342G > A; p.Gly448Arg) in the exon-10 of the SLC9A6 gene. By integrating computational approaches, our findings provided insights into the molecular mechanisms underlying the development of ID in individuals with SLC13A3 and SLC9A6 mutations.

CONCLUSION

We have utilized in-silico tools in the current study to examine the deleterious effects of the identified variants, which carry the potential to understand the genotype-phenotype relationships in neurodevelopmental disorders.

Autosomal recessive variants c.953A>C and c.97-1G>C in NSUN2 causing intellectual disability: a molecular dynamics simulation study of loss-of-function mechanisms

Introduction

Intellectual disability (ID) is a clinically and genetically heterogeneous disorder. It drastically affects the learning capabilities of patients and eventually reduces their IQ level below 70.

Methods

The current genetic study ascertained two consanguineous Pakistani families suffering from autosomal recessive intellectual developmental disorder-5 (MRT5). We have used exome sequencing followed by Sanger sequencing to identify the disease-causing variants.

Results and discussion

Genetic analysis using whole exome sequencing in these families identified two novel mutations in the NSUN2 (NM_017755.5). Family-A segregated a novel missense variant c.953A>C; p.Tyr318Ser in exon-9 of the NSUN2. The variant substituted an amino acid Tyr318, highly conserved among different animal species and located in the functional domain of NSUN2 known as "SAM-dependent methyltransferase RsmB/NOP2-type". Whereas in family B, we identified a novel splice site variant c.97-1G>C that affects the splice acceptor site of NSUN2. The identified splice variant (c.97-1G>C) was predicted to result in the skipping of exon-2, which would lead to a frameshift followed by a premature stop codon (p. His86Profs^*16). Furthermore, it could result in the termination of translation and synthesis of dysfunctional protein, most likely leading to nonsense-mediated decay. The dynamic consequences of NSUN2 missense variant was further explored together with wildtype through molecular dynamic simulations, which uncovered the disruption of NSUN2 function due to a gain in structural flexibility. The present molecular genetic study further extends the mutational spectrum of NSUN2 to be involved in ID and its genetic heterogeneity in the Pakistani population.

Wound Healing Potential and In Silico Appraisal of Convolvulus arvensis L. Methanolic Extract

Convolvulus arvensis L. is rich in phenolic compounds and traditionally used to treat wounds, skin ulcer, and inflammation. The current study is aimed at scientifically potentiating its traditional wound healing use. The methanolic extract of C. arvensis stem (CaME) was analyzed for HPLC and GC-MS analyses. The binding modes of active compounds were investigated against protein targets glycogen synthase kinase-3β (GSK-3β), transforming growth factor-beta (TGF-β), c-myc, and β-catenin by molecular docking followed by molecular dynamic simulations which revealed some conserved mode of binding as reported in crystal structures. The antioxidant potential of CaME was evaluated by in vitro methods such as 2,2-diphenyl-1-picrylhydrazyl (DPPH) radical scavenging, hydrogen peroxide scavenging, and ferric reducing power assays. Ointment formulations of 10 and 20% CaME were applied topically and evaluated for wound healing potency against the excisional wound on the skin of Wistar rats. Gentamycin (0.1%) served as standard therapy. The healing process was observed for 20 days in the form of wound size and epithelialization followed by histopathological evaluation of the wound area. Chemical characterization showed the presence of 7-hexadecenoic acid, 2-hexadecylicosan-1-ol, quercetin, gallic acid, ferulic acid, and other compounds. The plant extract exhibited significant in vitro antioxidant activity. The animals treated with 10% ointment showed moderate healing, whereas the treatment with 20% CaME revealed healing potential comparable to the standard 0.1% gentamycin as coevidenced from histopathological evaluation of skin. The study corroborates promising potential of C. arvensis on the healing of wounds, which possibly will be attributed to its antioxidant activity, fatty acids, quercetin, and gallic and caffeic acids, and binding potential of its phytoconstituents (phenolic acids) with wound healing targets.

Combining virtual screening and in vitro evaluation for the discovery of potential CYP11B2 inhibitors

Aim: To search for highly bioactive hits for CYP11B2 inhibitors by virtual screening and in vitro evaluation. Materials & methods: Virtual screening of potential CYP11B2 inhibitors was performed by molecular docking and molecular dynamics simulation. Compound activity was determined by in vitro evaluation using MTT and ELISA assays. Results & conclusion: Based on the results of molecular docking and molecular dynamics simulation, nine lead hits were selected for in vitro biochemical testing. All hits in in vitro experiments had lower inhibitory effects on cell proliferation and certain inhibitory effects on aldosterone secretion. These hits may be excellent candidates for CYP11B2 inhibitors.

Discovery of Rift Valley fever virus natural pan-inhibitors by targeting its multiple key proteins through computational approaches

The Rift Valley fever virus (RVFV) is a zoonotic arbovirus and pathogenic to both humans and animals. Currently, no proven effective RVFV drugs or licensed vaccine are available for human or animal use. Hence, there is an urgent need to develop effective treatment options to control this viral infection. RVFV glycoprotein N (GN), glycoprotein C (GC), and nucleocapsid (N) proteins are attractive antiviral drug targets due to their critical roles in RVFV replication. In present study, an integrated docking-based virtual screening of more than 6000 phytochemicals with known antiviral activities against these conserved RVFV proteins was conducted. The top five hit compounds, calyxin C, calyxin D, calyxin J, gericudranins A, and blepharocalyxin C displayed optimal binding against all three target proteins. Moreover, multiple parameters from the molecular dynamics (MD) simulations and MM/GBSA analysis confirmed the stability of protein-ligand complexes and revealed that these compounds may act as potential pan-inhibitors of RVFV replication. Our computational analyses may contribute toward the development of promising effective drugs against RVFV infection.

Structural Elucidation of Rift Valley Fever Virus L Protein towards the Discovery of Its Potential Inhibitors

Rift valley fever virus (RVFV) is the causative agent of a viral zoonosis that causes a significant clinical burden in domestic and wild ruminants. Major outbreaks of the virus occur in livestock, and contaminated animal products or arthropod vectors can transmit the virus to humans. The viral RNA-dependent RNA polymerase (RdRp; L protein) of the RVFV is responsible for viral replication and is thus an appealing drug target because no effective and specific vaccine against this virus is available. The current study reported the structural elucidation of the RVFV-L protein by in-depth homology modeling since no crystal structure is available yet. The inhibitory binding modes of known potent L protein inhibitors were analyzed. Based on the results, further molecular docking-based virtual screening of Selleckchem Nucleoside Analogue Library (156 compounds) was performed to find potential new inhibitors against the RVFV L protein. ADME (Absorption, Distribution, Metabolism, and Excretion) and toxicity analysis of these compounds was also performed. Besides, the binding mechanism and stability of identified compounds were confirmed by a 50 ns molecular dynamic (MD) simulation followed by MM/PBSA binding free energy calculations. Homology modeling determined a stable multi-domain structure of L protein. An analysis of known L protein inhibitors, including Monensin, Mycophenolic acid, and Ribavirin, provide insights into the binding mechanism and reveals key residues of the L protein binding pocket. The screening results revealed that the top three compounds, A-317491, Khasianine, and VER155008, exhibited a high affinity at the L protein binding pocket. ADME analysis revealed good pharmacodynamics and pharmacokinetic profiles of these compounds. Furthermore, MD simulation and binding free energy analysis endorsed the binding stability of potential compounds with L protein. In a nutshell, the present study determined potential compounds that may aid in the rational design of novel inhibitors of the RVFV L protein as anti-RVFV drugs.

The discovery of Zika virus NS2B-NS3 inhibitors with antiviral activity via an integrated virtual screening approach

With expanding recent outbreaks and a lack of treatment options, the Zika virus (ZIKV) poses a severe health concern. The availability of ZIKV NS2B-NS3 co-crystallized structures paved the way for rational drug discovery. A computer-aided structure-based approach was used to screen a diverse library of compounds against ZIKV NS2B-NS3 protease. The top hits were selected based on various binding free energy calculations followed by per-residue decomposition analysis. The selected hits were then evaluated for their biological potential with ZIKV protease inhibition assay and antiviral activity. Among 26 selected compounds, 8 compounds showed promising activity against ZIKV protease with a percentage inhibition of greater than 25 and 3 compounds displayed ∼50% at 10 µM, which indicates an enrichment rate of approximately 36% (threshold IC₅₀ < 10 µM) in the ZIKV-NS2B-NS3 protease inhibition assay. Of these, only one compound (23) produced whole-cell anti-ZIKV activity, and the binding mode of 23 was extensively analyzed through long-run molecular dynamics simulations. The current study provides a promising starting point for the further development of novel compounds against ZIKV.

Pharmacoinformatics and molecular dynamics simulation studies reveal potential covalent and FDA-approved inhibitors of SARS-CoV-2 main protease 3CLpro

The SARS-CoV-2 was confirmed to cause the global pandemic of coronavirus disease 2019 (COVID-19). The 3-chymotrypsin-like protease (3CLpro), an essential enzyme for viral replication, is a valid target to combat SARS-CoV and MERS-CoV. In this work, we present a structure-based study to identify potential covalent inhibitors containing a variety of chemical warheads. The targeted Asinex Focused Covalent (AFCL) library was screened based on different reaction types and potential covalent inhibitors were identified. In addition, we screened FDA-approved protease inhibitors to find candidates to be repurposed against SARS-CoV-2 3CLpro. A number of compounds with significant covalent docking scores were identified. These compounds were able to establish a covalent bond (C-S) with the reactive thiol group of Cys145 and to form favorable interactions with residues lining the substrate-binding site. Moreover, paritaprevir and simeprevir from FDA-approved protease inhibitors were identified as potential inhibitors of SARS-CoV-2 3CLpro. The mechanism and dynamic stability of binding between the identified compounds and SARS-CoV-2 3CLpro were characterized by molecular dynamics (MD) simulations. The identified compounds are potential inhibitors worthy of further development as COVID-19 drugs. Importantly, the identified FDA-approved anti-hepatitis-C virus (HCV) drugs paritaprevir and simeprevir could be ready for clinical trials to treat infected patients and help curb COVID-19. Communicated by Ramaswamy H. Sarma.

Enhanced Thermostability and Enzymatic Activity of Cel6A Variants from Thermobifida fusca by Empirical Domain Engineering (Short Title: Domain Engineering of Cel6A)

Cellulases are a set of lignocellulolytic enzymes, capable of producing eco-friendly low-cost renewable bioethanol. However, low stability and hydrolytic activity limit their wide-scale applicability at the industrial scale. In this work, we report the domain engineering of endoglucanase (Cel6A) of Thermobifida fusca to improve their catalytic activity and thermal stability. Later, enzymatic activity and thermostability of the most efficient variant named as Cel6A.CBC was analyzed by molecular dynamics simulations. This variant demonstrated profound activity against soluble and insoluble cellulosic substrates like filter paper, alkali-treated bagasse, regenerated amorphous cellulose (RAC), and bacterial microcrystalline cellulose. The variant Cel6A.CBC showed the highest catalysis of carboxymethyl cellulose (CMC) and other related insoluble substrates at a pH of 6.0 and a temperature of 60 °C. Furthermore, a sound rationale was observed between experimental findings and molecular modeling of Cel6A.CBC which revealed thermostability of Cel6A.CBC at 26.85, 60.85, and 74.85 °C as well as structural flexibility at 126.85 °C. Therefore, a thermostable derivative of Cel6A engineered in the present work has enhanced biological performance and can be a useful construct for the mass production of bioethanol from plant biomass.

In Vitro Antidiabetic, Anti-Obesity and Antioxidant Analysis of Ocimum basilicum Aerial Biomass and in Silico Molecular Docking Simulations with Alpha-Amylase and Lipase Enzymes

The present study explored phytochemicals, porcine pancreatic α-amylase (PPA) and lipase (PPL) inhibitory activities and antioxidant potential of polar and nonpolar extracts of the leaves and flowers of Ocimum basilicum and the in-silico mode of interaction between these enzymes and the major chemical constituents of the herb. The hexane extract (HE) and hydro-ethanolic extract (EE) obtained sequentially were used to estimate PPA and PPL inhibitory and antioxidant activities, total phenolic content (TPC) and total flavonoid content (TFC). Chemical constituents of the essential oils and HE were determined by GC-MS (Gas Chromatography-Mass Spectrometry). For PPA inhibition, IC₅₀ (µg/mL) of the extracts were 0.27-0.37, which were close to 0.24 of acarbose, while for PPL inhibition, IC₅₀ (µg/mL) of the extracts were 278.40-399.65, and that of Orlistat 145.72. The flowers EE was most potent antioxidant followed by leaves EE. The leaves EE had highest TPC and TFC followed of flowers EE. The essential oil of flowers had higher estragole (55%) than linalool (37%), while the essential oil of the leaves had higher linalool (42%) than estragole (38%). The HE of the flowers contained higher estragole (42%) than linalool (23%), while of the HE of the leaves too had higher estragole (65%) than linalool (18%). The in-silico molecular docking study showed linalool and estragole to have considerable PPA and PPL binding potential, which were further investigated through molecular dynamics simulations and binding free energy calculations. The PPA and PPL inhibitory activities of O. basilicum extracts and their notable antioxidant potential propose the herb as a multi-target complimentary medicine for diabetes, obesity and oxidative stress.

Deleterious Variants in WNT10A, EDAR, and EDA Causing Isolated and Syndromic Tooth Agenesis: A Structural Perspective from Molecular Dynamics Simulations

The dental abnormalities are the typical features of many ectodermal dysplasias along with congenital malformations of nails, skin, hair, and sweat glands. However, several reports of non-syndromic/isolated tooth agenesis have also been found in the literature. The characteristic features of hypohidrotic ectodermal dysplasia (HED) comprise of hypodontia/oligodontia, along with hypohidrosis/anhidrosis, and hypotrichosis. Pathogenic variants in EDA, EDAR, EDARADD, and TRAF6, cause the phenotypic expression of HED. Genetic alterations in EDA and WNT10A cause particularly non-syndromic/isolated oligodontia. In the current project, we recruited 57 patients of 17 genetic pedigrees (A-Q) from different geographic regions of the world, including Pakistan, Egypt, Saudi Arabia, and Syria. The molecular investigation of different syndromic and non-syndromic dental conditions, including hypodontia, oligodontia, generalized odontodysplasia, and dental crowding was carried out by using exome and Sanger sequencing. We have identified a novel missense variant (c.311G>A; p.Arg104His) in WNT10A in three oligodontia patients of family A, two novel sequence variants (c.207delinsTT, p.Gly70Trpfs*25 and c.1300T>G; p.Try434Gly) in EDAR in three patients of family B and four patients of family C, respectively. To better understand the structural and functional consequences of missense variants in WNT10A and EDAR on the stability of the proteins, we have performed extensive molecular dynamic (MD) simulations. We have also identified three previously reported pathogenic variants (c.1076T>C; p.Met359Thr), (c.1133C>T; p.Thr378Met) and (c.594_595insC; Gly201Argfs*39) in EDA in family D (four patients), E (two patients) and F (one patient), correspondingly. Presently, our data explain the genetic cause of 18 syndromic and non-syndromic tooth agenesis patients in six autosomal recessive and X-linked pedigrees (A-F), which expand the mutational spectrum of these unique clinical manifestations.