Biochimica et biophysica acta. Proteins and proteomics最新文献

Molecular dynamics (MD) simulations are a powerful approach to studying the structure and dynamics of proteins related to health and disease. Advances in the MD field allow modeling proteins with high accuracy. However, modeling metal ions and their interactions with proteins is still challenging. NPL4 is a zinc-binding protein and works as a cofactor for p97 to regulate protein homeostasis. NPL4 is of biomedical importance and has been proposed as the target of disulfiram, a drug recently repurposed for cancer treatment. Experimental studies proposed that the disulfiram metabolites, bis-(diethyldithiocarbamate)‑copper and cupric ions, induce NPL4 misfolding and aggregation. However, the molecular details of their interactions with NPL4 and consequent structural effects are still elusive. Here, biomolecular simulations can help to shed light on the related structural details. To apply MD simulations to NPL4 and its interaction with copper the first important step is identifying a suitable force field to describe the protein in its zinc-bound states. We examined different sets of non-bonded parameters because we want to study the misfolding mechanism and cannot rule out that the zinc may detach from the protein during the process and copper replaces it. We investigated the force-field ability to model the coordination geometry of the metal ions by comparing the results from MD simulations with optimized geometries from quantum mechanics (QM) calculations using model systems of NPL4. Furthermore, we investigated the performance of a force field including bonded parameters to treat copper ions in NPL4 that we obtained based on QM calculations.

Cancer pathologies are associated with the unfolding and aggregation of most recurring mutations in the DNA Binding Domain (DBD) of p53 that coordinate the destabilization of protein. Substitution at the 175th codon with arginine to histidine (R175H, a mutation of large to small side-chain amino acid) destabilizes the DBD by 3 kcal/mol and triggers breasts, lung cancer, etc. Stabilizing the p53 mutant by small molecules offers an attractive drug-targeted anti-cancer therapy. The thiosemicarbazone (TSC) molecules NPC and DPT are known to act as zinc-metallochaperones to reactivate p53R175H. Here, a combination of LESMD simulations for 10 TSC conformations with a p53R175H receptor, single ligand-protein conformation MD, and ensemble docking with multiple p53R175H conformations observed during simulations is suggested to identify the potential binding site of the target protein in light of their importance for the direct TSC – p53R175H binding. NPC binds mutant R175H in the loop region L2-L3, forming pivotal hydrogen bonds with HIS175, pi‑sulfur bonds with TYR163, and pi-alkyl linkages with ARG174 and PRO190, all of which are contiguous to the zinc-binding native site on p53DBD. DPT, on the other hand, was primarily targeting alternative binding sites such as the loop-helix L1/H2 region and the S8 strand. The similar structural characteristics of TSC-bound p53R175H complexes with wild-type p53DBD are thought to be attributable to involved interactions that favour binding free energy contributions of TSC ligands. Our findings may be useful in the identification of novel pockets with druggable properties.

Clustered regularly interspaced short palindromic repeats (CRISPR) in many prokaryotes functions as an adaptive immune system against mobile genetic elements. A heterologous ribonucleoprotein silencing complex composed of CRISPR-associated (Cas) proteins and a CRISPR RNA (crRNA) neutralizes the incoming mobile genetic elements. The type I and III silencing complexes commonly include a protein-helical backbone of several copies of identical subunits, for example, Cas7 in the type I silencing complex.

In this study, we structurally characterized type I-B Cas7 (Csh2 from Thermobaculum terrenum; TterCsh2). The revealed crystal structure of TterCsh2 shows a typical glove-like architecture of Cas7, which consists of a palm, a thumb, and a finger domain. Csh2 proteins have 5 conserved sequence motifs that are arranged to form a presumable crRNA-binding site in the TterCsh2 structure. This crRNA binding site of TterCsh2 is structurally and potentially comparable to those observed in helix-forming Cas7 structures in other sub-types. Analysis of the reported Cas7 structures and their sequences suggests that Cas7s can be divided into at least two sub-classes. These data will broaden our understanding on the Cascade complex of CRISPR/Cas systems.

The short peptidoglycan recognition protein (PGRP-S) of the innate immune system recognizes the invading microbes through binding to their cell wall molecules. In order to understand the mode of binding of PGRP-S to bacterial cell wall molecules, the structure of the complex of camel PGRP-S (CPGRP-S) with hexanoic acid has been determined at 2.07 Å resolution. Previously, we had reported the structures of CPGRP-S in the native unbound state as well as in the complexed forms with the components of various bacterial cell wall molecules such as peptidoglycan (PGN), lipopolysaccharide (LPS), lipoteichoic acid (LTA), mycolic acid (MA) and other fatty acids. These structures revealed that CPGRP-S formed two homodimers which were designated as A-B and CD dimers. It also showed that the fatty acids bind to CPGRP-S in the binding site at the A-B dimer while the non-fatty acids were shown to bind at the interfaces of both A-B and CD dimers. The present structure of the complex of CPGRP-S with hexanoic acid (HA) showed that HA binds to CPGRP-S at the interface of CD dimer. HA was located in the same groove at the CD interface which was occupied by non-fatty acids such as PGN, LPS and LTA and interacts with residues from both C and D molecules. HA is firmly held in the groove with several hydrogen bonds and a number of van der Waals contacts. This is the first structure which reports the binding of a fatty acid in the cleft at the interface of CD dimer.

Metabolic stability of proteins plays a vital role in various dedicated cellular processes. Traditional methods of measuring the metabolic stability are time-consuming and expensive. Therefore, we developed a more efficient computational approach to understand the protein dynamic action mechanisms in biological process networks. In this study, we collected 341 short-lived proteins and 824 non-short-lived proteins from U2OS; 342 short-lived proteins and 821 non-short-lived proteins from HEK293T; 424 short-lived proteins and 1153 non-short-lived proteins from HCT116; and 384 short-lived proteins and 992 non-short-lived proteins from RPE1. The proteins were encoded by GO and KEGG enrichment scores based on the genes and their neighbors in STRING, resulting in 20,681 GO term features and 297 KEGG pathway features. We also incorporated the protein interaction information from STRING into the features and obtained 19,247 node features. Boruta and mRMR methods were used for feature filtering, and IFS method was used to obtain the best feature subsets and create the models with the highest performance. The present study identified 42 features that did not appear in previous studies and classified them into eight groups according to their functional annotation. By reviewing the literature, we found that the following three functional groups were critical in determining the stability of proteins: synaptic transmission, post-translational modifications, and cell fate determination. These findings may serve as a valuable reference for developing drugs that target protein stability.

Hesperidin and narirutin are the major citrus flavanones. Several studies have associated these compounds with pancreatic β-cell survival through their capacity to reduce oxidative stress, inflammation, and inhibit apoptosis. However, the molecular mechanisms of action of flavanones in pancreatic β-cells under high-glycemic stress is still largely unknown. Therefore, this study aimed to decipher molecular mechanisms of flavanone metabolites in pancreatic β-cells treated with high glucose concentration using untargeted shotgun proteomics. We identified 569 proteins differentially expressed in cells exposed to hesperetin 7-glucuronide (H7G) and 265 in cells exposed to 3-(4′-hydroxyphenyl) propanoic acid (PA). Comparison of global proteomic profiles suggest that these metabolites could counteract changes in protein expression induced by high glucose stress. The bioinformatic analyses suggested that H7G and PA modulated the expression of proteins involved in cell adhesion, cell signaling, metabolism, inflammation, and protein processing in endoplasmic reticulum (ER) pathways. Taken together, this study suggests that H7G and PA can modulate the expression of proteins that may prevent dysfunction of pancreatic β-cells under stress induced by high glucose.

Inhibition of highly ordered cross-β-sheet-rich aggregates of misfolded amyloid proteins using rationally designed sequence-based short peptides is a promising therapeutic strategy for the treatment of neurodegenerative diseases. Here, we have explored the anti-amyloidogenic potency of a rationally designed hexapeptide (Tyr-Pro-Gln-Ile-Pro-Asn) on in vitro hen egg white lysozyme (HEWL) amyloid fibril formation at acidic pH and physiological pH using computational docking as well as various biophysical techniques such as fluorescence spectroscopy, UV–vis spectroscopy, FTIR spectroscopy, confocal microscopy and TEM. The peptide was designed based on the aggregation-prone region (APR) of HEWL and thus referred to as SqP1 (Sequence-based Peptide 1). SqP1 showed over 70% inhibition of HEWL amyloid formation at pH 2.2 and approximately 50% inhibition at pH 7.5. We propose that SqP1 binds to the APR of HEWL and interacts strongly with the Trp62/Trp63, ultimately stabilizing monomeric HEWL at both the pH conditions and preventing conformation changes in the structure of HEWL, leading to the formation of amyloidogenic fibrillar structures. A sequence-based peptide inhibitor of HEWL amyloid formation was not reported previously, making this a critical study that will further emphasize the importance of short synthetic peptides as amyloid inhibitors.

Natural and artificial nucleases have extensive applications in biotechnology and biomedicine. The exploration of protein with potential DNA cleavage activity also inspires the design of artificial nuclease and helps to understand the physiological process of DNA damage. In this study, we engineered four human cytochrome c (Cyt c) mutants (N52S, N52A, I81N, and I81D Cyt c), which showed enhanced DNA cleavage activity and degradation in comparison with WT Cyt c, especially under acidic conditions. The mechanism assays revealed that the superoxide (O₂^•−) plays an important role in the nuclease reaction. The kinetic assays showed that the peroxidase activity of the I81D Cyt c mutant enhanced up to 9-fold at pH 5. This study suggests that the mutations of Ile81 and Asn52 in Ω-loop C/D are critical for the nuclease activity of Cyt c, which may have physiological significance in DNA damage and potential applications in biomedicine.

UDP-glucuronosyltransferase 2B15 (UGT2B15) is a crucial phase II drug-metabolizing enzyme, which glucuronidates various compounds, including clinical drugs and hormones. Mutants might affect glucuronidation, leading to a disruption of drug metabolism in vivo and decrease of therapeutic effect. Here, we mainly analyzed two representative mutants, H401P and L446S, on UGT2B15 activity using glucuronidation assays, molecular dynamic (MD) simulation and X-ray diffraction methods. The enzyme activity of L446S obviously increased six-fold than the wild type, although the enzyme activities of P191L, T374A, and H401P were lost apparently. Furthermore, we used MD simulations to calculate the energy change in the catalytic process of H401P and L446S, and the results indicated the free binding energies of H401P mutant to oxazepam and UDPGA were −30.98 ± 1.00 kcal/mol and −36.42 ± 1.04 kcal/mol, respectively, increased obviously compared to wild type, suggesting the mutation on position 401 had a crucial effect on the catalysis. Moreover, the three-dimensional structure of UGT2B15 C-terminal domain L446S was determined through protein crystallography and X-ray diffraction technology and the results suggested that one more hydrogen bonding between S446 and K410 was formed in the S446 crystal structure, compared to the wild type. Isothermal titration calorimetry assay further revealed the K_d values of C-terminal domain of UGT2B15 harbored L446S towards the cofactor UDPGA was similar to the value of wild type. Above all, our results pointed out that H401P and L446S affected the enzyme activity by different mechanism. Our work provided a helpful mechanism for variance explained in the UGTs catalyzation process.

LexA, a well-characterized transcriptional repressor of SOS genes in heterotrophic bacteria, has been shown to regulate diverse genes in cyanobacteria. An earlier study showed that LexA overexpression in a cyanobacterium, Anabaena sp. PCC7120 reduces its tolerance to Cd stress. This was later shown to be due to modulation of photosynthetic redox poising by LexA under Cd stress. However, due to the global regulatory nature of LexA and the prior prediction of AnLexA-box in a few heavy metal-responsive genes, we speculated that LexA has a broad role in Cd tolerance, with regulation over a variety of Cd stress-responsive genes in addition to photosynthetic genes. Thus, to further expand the knowledge on the regulatory role of LexA in Cd stress tolerance, a cytosolic proteome profiling of Anabaena constitutively overexpressing LexA upon Cd stress was performed. The proteomic study revealed 25 differentially accumulated proteins (DAPs) in response to the combined effect of LexA overexpression and Cd stress, and the other 11 DAPs exclusively in response to either LexA overexpression or Cd stress. The 36 identified proteins were related with a variety of functions, including photosynthesis, C-metabolism, antioxidants, protein turnover, post-transcriptional modifications, and a few unknown and hypothetical proteins. The regulation of LexA on corresponding genes, and six previously reported Cd efflux transporters, was further validated by the presence of AnLexA-boxes, transcript, and/or promoter analyses. In a nutshell, this study identifies the regulation of Anabaena LexA on several Cd stress-responsive genes of various functions, hence expanding the regulatory role of LexA under Cd stress.