Understanding a point mutation signature D54K in the caspase activation recruitment domain of NOD1 capitulating concerted immunity via atomistic simulation.

Point mutation D54K in the human N-terminal caspase recruitment domain (CARD) of nucleotide-binding oligomerization domain -1 (NOD1) abrogates an imperative downstream interaction with receptor-interacting protein kinase (RIPK2) that entails combating bacterial infections and inflammatory dysfunction. Here, we addressed the molecular details concerning conformational changes and interaction patterns (monomeric-dimeric states) of D54K by signature-based molecular dynamics simulation. Initially, the sequence analysis prioritized D54K as a pathogenic mutation, among other variants, based on a sequence signature. Since the mutation is highly conserved, we derived the distant ortholog to predict the sequence and structural similarity between native and mutant. This analysis showed the utility of 33 communal core residues associated with structural-functional preservation and variations, concurrently served to infer the cryptic hotspots Cys39, Glu53, Asp54, Glu56, Ile57, Leu74, and Lys78 determining the inter helical fold forming homodimers for putative receptor interaction. Subsequently, the atomistic simulations with free energy (MM/PB(GB)SA) calculations predicted structural alteration that takes place in the N-terminal mutant CARD where coils changed to helices (45 α3- L4-α4-L6- α683) in contrast to native (45T2-L4-α4-L6-T483). Likewise, the C-terminal helices 93T1-α7105 connected to the loops distorted compared to native 93α6-L7105 may result in conformational misfolding that promotes functional regulation and activation. These structural perturbations of D54K possibly destabilize the flexible adaptation of critical homotypic NOD1CARD-CARDRIPK2 interactions (α4Asp42-Arg488α5 and α6Phe86-Lys471α4) is consistent with earlier experimental reports. Altogether, our findings unveil the conformational plasticity of mutation-dependent immunomodulatory response and may aid in functional validation exploring clinical investigation on CARD-regulated immunotherapies to prevent systemic infection and inflammation.Communicated by Ramaswamy H. Sarma.

Determination of potential combination of non-ß-lactam, ß-lactam, and ß-lactamase inhibitors/ß-lactam enhancer against class D oxacillinases producing Acinetobacter baumannii: Evidence from in-vitro, molecular docking and dynamics simulation.

Carbapenem-resistant Acinetobacter baumannii, a predominant nosocomial pathogen in hospitals of intensive care units, is associated with bacteremia and ventilator-associated pneumonia with a high-risk mortality rate. To increase the effectiveness of the ß-lactam (BL) antibiotics, the use of ß-lactamase inhibitors (BLI) acts as a booster when given in combination with BL antibiotics. To this aspect, we selected BL antibiotics of cefiderocol, cefepime, non-BL antibiotic eravacycline, BLI of durlobactam, avibactam, and a ß-lactam enhancer (BLE) of zidebactam. To prove our hypothesis, we determined the minimum inhibitory concentration (MIC) of various BL or non-BL/BLI or BLE combinations using broth microdilution method followed by in silico analysis of molecular docking, molecular dynamics (MD) simulation, and molecular mechanics Poisson-Boltzmann surface area (MM-PBSA) identifies the potential combination. In MIC testing, eravacycline, cefepime/zidebactam, cefiderocol/zidebactam, and eravacycline in combination with zidebactam or durlobactam were found to be effective against oxacillinases (OXAs) (OXA-23/24/58 like) expressing A. baumannii isolates. The docking results of the selected ligands toward OXA-23, OXA-24, and OXA-58 had an excellent binding score ranging from -5.8 to -9.3 kcal/mol. Further, the docked complexes were subjected and evaluated using gromacs for molecular dynamics simulation of 50 ns toward selected class D OXAs. The binding energies obtained from MM-PBSA shed light on the binding efficiencies of each non-BL, BL, and BLI/BLE, thereby helping us to propose the drug combinations. Based on the MD trajectories scoring acquired, we propose using eravacycline, cefepime/zidebactam, cefiderocol/zidebactam, and eravacycline in combination with durlobactam or zidebactam would be promising for treating OXA-23, OXA-24, and OXA-58 like expressing A. baumannii infections.

A new horizon in the phosphorylated sites of AGA: the structural impact of C163S mutation in aspartylglucosaminuria through molecular dynamics simulation.

Aspartylglucosaminuria (AGU) is a lysosomal storage disorder caused by insufficient aspartylglucosaminidase (AGA) activity leading to chronic neurodegeneration. We utilized the PhosphoSitePlus tool to identify the AGA protein's phosphorylation sites. The phosphorylation was induced on the specific residue of the three-dimensional AGA protein, and the structural changes upon phosphorylation were studied via molecular dynamics simulation. Furthermore, the structural behaviour of C163S mutation and C163S mutation with adjacent phosphorylation was investigated. We have examined the structural impact of phosphorylated forms and C163S mutation in AGA. Molecular dynamics simulations (200 ns) exposed patterns of deviation, fluctuation, and change in compactness of Y178 phosphorylated AGA protein (Y178-p), T215 phosphorylated AGA protein (T215-p), T324 phosphorylated AGA protein (T324-p), C163S mutant AGA protein (C163S), and C163S mutation with Y178 phosphorylated AGA protein (C163S-Y178-p). Y178-p, T215-p, and C163S mutation demonstrated an increase in intramolecular hydrogen bonds, leading to greater compactness of the AGA forms. Principle component analysis (PCA) and Gibbs free energy of the phosphorylated/C163S mutation structures exhibit transition in motion/orientation than Wild type (WT). T215-p may be more dominant among these than the other studied phosphorylated forms. It might contribute to hydrolyzing L-asparagine functioning as an asparaginase, thereby regulating neurotransmitter activity. This study revealed structural insights into the phosphorylation of Y178, T215, and T324 in AGA protein. Additionally, it exposed the structural changes of the C163S mutation and C163S-Y178-p of AGA protein. This research will shed light on a better understanding of AGA's phosphorylated mechanism.Communicated by Ramaswamy H. Sarma.

Structural localization of pathogenic mutations in the central nucleotide-binding domain (NBD) of nucleotide-binding oligomerization domain-2 (NOD2) protein and their inference in inflammatory disorders.

The human NBD domain which is centrally located in the NOD2 protein displays an essential role in oligomerization and initiates the immune response via CARD-RIPK2 interaction. The mutations associated with the NBD domain have been largely implicated in inflammatory disorders such as Blau syndrome and sarcoidosis. This study aims to determine the structural and phenotypic effect of a lethal mutation that occurs in the NBD domain which has an axiomatic impact on protein dysfunction. Initially, the most deleterious missense mutations were screened through various in silico analysis. Out of 33 variants, I-Mutant 3.0, SIFT, PolyPhen 2, Align GVGD, PHD SNP and SNP&GO have statistically identified 5 variants (R42W, D90E, E91K, G189D & W198L) as less stable, deleterious and damaging. Our predicted models have paved the way to understand the various structural properties such as physiochemical, secondary structural arrangements and stabilizing residues in folding associated with the native and mutant NBD domain especially of the functionally important regions. From the aforementioned results, R42W and G189D were found to be the more predominant among the mutants. Precisely, through molecular simulation, we have strongly justified the significant conformational disruption of R42W and G189D through the stabilization factors, folding and essential dynamics. Conclusively, these regions (α341-44, α13185-191 and ß6133-143ß7) seem to adopt such structures that are not conducive to wild-type-like functionality. Our prediction and validation of lethal mutations based on structural stability may be useful for conducting experimental studies in detail to uncover the protein deregulation leading to inflammatory disorders.

Delineating the folding perturbations and molecular mechanisms of Thr-Ala 642 mutation in Rab-GTPase activating protein Akt substrate of 160kDa and its impact on the aetiology of diabetes.

Diabetes Mellitus is a complex metabolic disorder with one of the highest prevalence rates in the world. The present study probes into the Thr-Ala 642 mutation of Akt substrate of 160 kDa (AS160) which has been implicated in diabetes by the dysregulation of glucose transported vesicle 4 (GLUT4) translocation. Our study provides a possible evidence on structural basis dysfunction of AS160 and how the association of phosphorylated AS160 with 14-3-3, a downstream binding partner regulating GLUT4 translocation got disrupted due to T642A mutation. We initially derived the disease-causing mutation (Thr642Ala) among others through in-silico based statistical analysis. Subsequently, we interpreted the perturbation induced in the structural arrangement and their impaired interaction in core regions. Due to mutation, the key interfacial interactions between AS160-14-3-3 were changing from Thr642-Asp756, Thr642-Asp757, and Thr642-Lys659 for phosphorylated form, to Ala642-Val681 for mutant. Further, for phosphorylated AS160 the hotspot residues observed were Glu-629, Gln-635, His-641, Lys-653 and Arg-842 which changed to Arg-637, His-641 for mutant. Eventually, the molecular dynamics analysis revealed that local region for phosphorylation site of AS160 is reducing the flexibility, whereas mutation is making the region more flexible. Principal component analysis and Free energy landscape analysis together reveals phosphorylated AS160 is occupying less phase space with more stable landscape when compared to mutant. Our study strongly confers the destabilizing effect of the point mutation at a conserved site providing novel insights right down to the residual level of the conformational dysregulation of AS160 implicating deadly disease.

Streptozocin; a GLUT2 binding drug, interacts with human serum albumin at loci h6DOM3-h7DOM3.

Streptozocin (STZ) is a broad range antibiotic, highly genotoxic, antineoplastic and hyperglycemic. HSA is the most abundant protein in physiology and it binds to almost all exogenic and endogenic ligands, including drugs. STZ-induced fluorescence quenching of HSA has been done at pH 7.4, pH 3.5 and at pH 7.4 with 4.5 M urea at temperatures 286 K, 291 K, and 306 K. Ksv found to be 103 M-1, binding constant 1.5X103M-1 and binding sites ~1. But, Ksv for HSA and glucopyranose interaction was found lesser than that of HSA-STZ binding. Binding of STZ/glucopyranose on HSA seems to result in complex formation as calculated Kq > 1010 M-1 s-1. The number of binding sites, binding constants, and binding energies were increased with temperature. The ΔG0, ΔH0, and ΔS0 for HSA-STZ interaction were found to be -17.7 × 103 J·mol-1; 2.34 × 105 J·mol-1 and 841 JK-1 mol-1 respectively at pH 7.4 and 291 K. The comparative bindings of N, F and I states of HSA with STZ and their molecular docking analyses indicate that IIIA-B junction (i.e., inter-helix h6DOM3-h7DOM3) is the probable binding site, a locus close to fatty acid binding site-5. These results could be useful for therapeutic and analytical exploitation of STZ, as albumin used as the vehicle for drug delivery.

Allura red rapidly induces amyloid-like fibril formation in hen egg white lysozyme at physiological pH.

Allura red (AR) is an artificial azo dye mostly used in food industries and has potential health risks. We examined the role of AR in amyloidogenesis using hen egg white lysozyme (HEWL) at pH 7.0. The amyloidogenic induction properties of AR in HEWL were identified by circular dichroism (CD), turbidity, intrinsic fluorescence, light scattering, transmission electron microscopy (TEM), and molecular dynamic simulation studies. Turbidity and light scattering measurements showed that HEWL becomes aggregated in the presence of 0.03-15.0 mM of AR at pH 7.0 but not at very low AR concentrations (0.01-0.28 mM). However, AR-induced aggregation is a kinetically rapid process, with no observable lag phase and saturation within 6 s. The kinetics results suggested that the HEWL aggregation induced by AR is very rapid. The CD results demonstrated that the total ß-sheet content of HEWL was increased in the AR treated samples. The TEM results are established that AR-induced aggregates had amyloid-like structures. Molecular dynamics simulations analysis showed that the bound AR-HEWL structures were highly favored compared to unbound structures. The mechanism of AR-induced amyloid fibril formation may involve electrostatic, hydrogen bonding, and hydrophobic interactions.

A molecular simulation analysis of vitamin D targets interleukin 13 (IL13) as an alternative to mometasone in asthma.

Asthma, a chronic lung disease characterized by obstruction of airway passage is characterized by inflammation and hyperresponsiveness with increase in the number of eosinophils. Interleukin-13, plays a significant role in causing inflammation during an asthmatic attack by bronchial constriction. Mometasone, a glucocorticoid has been used as the first line of administration for people affected with asthma for almost a decade. However, in several cases, people treated with mometasone have faced systemic and local side effects. To reduce these side effects, we hypothesized vitamin D that can be used as a substitute to mometasone. For this purpose, we employed the use of molecular docking and simulation studies for comparative study. The docking studies revealed the binding residues of interleukin-13 which are bound to the active site. Among all, we noticed three binding residue Leu83, His84 and Arg86 common for both mometasone and vitamin D. Also, the binding energies share a significant similarity between them. The docked complexes of mometasone and vitamin D with interleukin-13 were evaluated with molecular dynamics simulation. Consistently, the MD analysis uncovered the interesting note on conformational adaptation between the complexes as well as that vitamin D has the complementary binding efficiency to interleukin-13 as compared to mometasone. The substitution of vitamin D might provide a promising gateway to reduce the side effects caused by mometasone and also reduce the cost for treatment of asthma patients.

Expediting dynamics approach to understand the influence of 14-3-3ζ causing metastatic cancer through the interaction of YAP1 and ß-TRCP.

The 14-3-3ζ protein acts as a molecular switch in regulating the TGF-ß pathway, which alters from a tumor suppressor in the early stage of breast cancer to a promoter of metastasis in the late stage. This change is due to the binding of 14-3-3ζ with YAP1 and ß-TRCP in premalignant and cancer cells, respectively. Owing to this inappropriate role of 14-3-3ζ when involved in cancer and metastasis, we predicted that Gln15, Glu17, Tyr211, and Gln219 are hotspot residues of 14-3-3ζ during its interaction with YAP1 protein. Similarly, we identified Gln15, Tyr211, Leu216, and Leu220 as hotspot residues of 14-3-3ζ during its interaction with ß-TRCP protein. Targeting these residues of 14-3-3ζ can prevent cancer and metastasis caused by malfunctioning of the TGF-ß pathway. In this work, we also predicted that YAP1 is an intrinsically disordered protein (IDP), and such proteins bind with other proteins via either an induced fit or a conformational selection mechanism. Intuitively, we found that 14-3-3ζ has high affinity towards phosphorylated YAP1 at Ser127 rather than unphosphorylated YAP1, which is in close agreement with previously reported experimental works. Thus, we performed an analysis by molecular dynamics simulations to reveal the conformational changes in YAP1 after phosphorylation at the atomistic level. Our work clearly illustrates the effect of phosphorylation on YAP1 in terms of conformational changes and the regulation of its function.