T54R mutation destabilizes the dimer of superoxide dismutase 1T54R by inducing steric clashes at the dimer interface

Identification and characterization of nuclear localization signals in the circumsporozoite protein of Plasmodium falciparum

Secretory proteins of Plasmodium exhibit differential spatial and functional activity within the host cell nucleus. However, the nuclear localization signals (NLSs) for these proteins remain largely uncharacterized. In this study, we have identified and characterized two NLSs in the circumsporozoite protein of Plasmodium falciparum (Pf-CSP). Both NLSs in the Pf-CSP contain clusters of lysine and arginine residues essential for specific interactions with the conserved tryptophan and asparagine residues of importin-α, facilitating nuclear translocation of Pf-CSP. While the two NLSs of Pf-CSP function independently and are both crucial for nuclear localization, a single NLS of Pf-CSP leads to weak nuclear localization. These findings shed light on the mechanism of nuclear penetrability of secretory proteins of Plasmodium proteins.

Unveiling local and global conformational changes and allosteric communications in SOD1 systems using molecular dynamics simulation and network analyses

BACKGROUND

Amyotrophic lateral sclerosis (ALS) is a serious neurodegenerative disorder affecting nerve cells in the brain and spinal cord that is caused by mutations in the superoxide dismutase 1 (SOD1) enzyme. ALS-related mutations cause misfolding, dimerisation instability, and increased formation of aggregates. The underlying allosteric mechanisms, however, remain obscure as far as details of their fundamental atomistic structure are concerned. Hence, this gap in knowledge limits the development of novel SOD1 inhibitors and the understanding of how disease-associated mutations in distal sites affect enzyme activity.

METHODS

We combined microsecond-scale based unbiased molecular dynamics (MD) simulation with network analysis to elucidate the local and global conformational changes and allosteric communications in SOD1 Apo (unmetallated form), Holo, Apo_CallA (mutant and unmetallated form), and Holo_CallA (mutant form) systems. To identify hotspot residues involved in SOD1 signalling and allosteric communications, we performed network centrality, community network, and path analyses.

RESULTS

Structural analyses showed that unmetallated SOD1 systems and cysteine mutations displayed large structural variations in the catalytic sites, affecting structural stability. Inter- and intra H-bond analyses identified several important residues crucial for maintaining interfacial stability, structural stability, and enzyme catalysis. Dynamic motion analysis demonstrated more balanced atomic displacement and highly correlated motions in the Holo system. The rationale for structural disparity observed in the disulfide bond formation and R143 configuration in Apo and Holo systems were elucidated using distance and dihedral probability distribution analyses.

CONCLUSION

Our study highlights the efficiency of combining extensive MD simulations with network analyses to unravel the features of protein allostery.

Characterizing proteomic and transcriptomic features of missense variants in amyotrophic lateral sclerosis genes

Within recent years, there has been a growing number of genes associated with amyotrophic lateral sclerosis (ALS), resulting in an increasing number of novel variants, particularly missense variants, many of which are of unknown clinical significance. Here, we leverage the sequencing efforts of the ALS Knowledge Portal (3864 individuals with ALS and 7839 controls) and Project MinE ALS Sequencing Consortium (4366 individuals with ALS and 1832 controls) to perform proteomic and transcriptomic characterization of missense variants in 24 ALS-associated genes. The two sequencing datasets were interrogated for missense variants in the 24 genes, and variants were annotated with gnomAD minor allele frequencies, ClinVar pathogenicity classifications, protein sequence features including Uniprot functional site annotations, and PhosphoSitePlus post-translational modification site annotations, structural features from AlphaFold predicted monomeric 3D structures, and transcriptomic expression levels from Genotype-Tissue Expression. We then applied missense variant enrichment and gene-burden testing following binning of variation based on the selected proteomic and transcriptomic features to identify those most relevant to pathogenicity in ALS-associated genes. Using predicted human protein structures from AlphaFold, we determined that missense variants carried by individuals with ALS were significantly enriched in β-sheets and α-helices, as well as in core, buried or moderately buried regions. At the same time, we identified that hydrophobic amino acid residues, compositionally biased protein regions and regions of interest are predominantly enriched in missense variants carried by individuals with ALS. Assessment of expression level based on transcriptomics also revealed enrichment of variants of high and medium expression across all tissues and within the brain. We further explored enriched features of interest using burden analyses and identified individual genes were indeed driving certain enrichment signals. A case study is presented for SOD1 to demonstrate proof-of-concept of how enriched features may aid in defining variant pathogenicity. Our results present proteomic and transcriptomic features that are important indicators of missense variant pathogenicity in ALS and are distinct from features associated with neurodevelopmental disorders.

Role of distal sites in enzyme engineering

The limitations associated with natural enzyme catalysis have triggered the rise of the field of protein engineering. Traditional rational design was based on the analysis of protein structural information and catalytic mechanisms to identify key active sites or ligand binding sites to reshape the substrate pocket. The role and significance of functional sites in the active center have been studied extensively. With a deeper understanding of the structure-catalysis relationship map, the entire protein molecule can be filled with residues that play a substantial role in its structure and function. However, the catalytic mechanism underlying distal mutations remains unclear. The aim of this review was to highlight the criticality of the distal site in enzyme engineering based on the following three aspects: What can distal mutations exert on function from mutability landscape? How do distal sites influence enzyme function? How to predict and design distal mutations? This review provides insights into the catalytic mechanism of enzymes from the global interaction network, knowledge from sequence-structure-dynamics-function relationships, and strategies for distal mutation-based protein engineering.

Computational insight into in silico analysis and molecular dynamics simulation of the dimer interface residues of ALS-linked hSOD1 forms in apo/holo states: a combined experimental and bioinformatic perspective

The aggregation of misfolded SOD1 proteins in neurodegenerative illnesses is a key pathological hallmark in amyotrophic lateral sclerosis (ALS). SOD1 is stabilized and enzymatically activated after binding to Cu/Zn and forming intramolecular disulfide. SOD1 aggregation/oligomerization is triggered by the dissociation of Cu and/or Zn ions. Therefore, we compared the possible effects of ALS-associated point mutations of the holo/apo forms of WT/I149T/V148G SOD1 variants located at the dimer interface to determine structural characterization using spectroscopic methods, computational approaches as well as molecular dynamics (MD) simulations. Predictive results of computational analysis of single-nucleotide polymorphisms (SNPs) suggested that mutant SOD1 has a deleterious effect on activity and structure destabilization. MD data analysis indicated that changes in flexibility, stability, hydrophobicity of the protein as well as increased intramolecular interactions of apo-SOD1 were more than holo-SOD1. Furthermore, a decrease in enzymatic activity in apo-SOD1 was observed compared to holo-SOD1. Comparative intrinsic and ANS fluorescence results of holo/apo-WT-hSOD1 and mutants indicated structural alterations in the local environment of tryptophan residue and hydrophobic patches, respectively. Experimental and MD data supported that substitution effect and metal deficiency of mutants (apo forms) in the dimer interface may promote the tendency to protein mis-folding and aggregation, consequently disrupting the dimer-monomer equilibrium and increased propensity to dissociation dimer into SOD-monomer ultimately leading to loss of stability and function. Overall, data analysis of apo/holo SOD1 forms on protein structure and function using computational and experimental studies will contribute to a better understanding of ALS pathogenicity.

Amyotrophic lateral sclerosis disease-related mutations disrupt the dimerization of superoxide dismutase 1 - A comparative molecular dynamics simulation study

More than 150 genes are involved in amyotrophic lateral sclerosis (ALS), with superoxide dismutase 1 (SOD1) being one of the most studied. Mutations in SOD1 gene, which encodes the enzyme SOD1 is the second most prevalent and studied cause of familial ALS. SOD1 is a ubiquitous, homodimeric metalloenzyme that forms a critical component of the cellular defense against reactive oxygen species. Several mutations in the SOD1 enzyme cause misfolding, dimerization instability, and increased aggregate formation in ALS. However, there is a lack of information on the dimerization of SOD1 monomers and the mechanistic underpinnings on how the pathogenic mutations disrupt the dimerization mechanism. Here, we presented microsecond-scale molecular dynamics (MD) simulations to unravel how interface-based mutations compromise SOD1 dimerization and provide mechanistic understanding into the corresponding process using WT and three interface-based mutant systems (A4V, T54R, and I113T). Structural stability analysis showed that the mutant systems displayed disparate variations in the catalytic sites which may directly alter the stability and activity of the SOD1 enzyme. Based on the dynamic network analysis and principal component analysis, it has been identified that the mutations weakened the correlated motions along the dimer interface and altered the protein conformational behavior, thus weakening the stability of dimer formation. Moreover, the simulation results identified crucial residues such as G51, D52, G114, I151, and Q153 in establishing the dimerization interaction network, which were weakened or absent in the presence of interfacial mutants. Surface potential analysis on mutant systems also displayed changes in the dimerization potential, thus showing the unfavorable dimer formation. Furthermore, network analysis identified the hotspot residues necessary for SOD1 signal transduction which were surprisingly found in the catalytic sites rather than the anticipated dimerization interface.

HYPK coordinates degradation of polyneddylated proteins by autophagy

Selective degradation of protein aggregates by macroautophagy/autophagy is an essential homeostatic process of safeguarding cells from the effects of proteotoxicity. Among the ubiquitin-like proteins, NEDD8 conjugation to misfolded proteins is prominent in stress-induced protein aggregates, albeit the function of neddylation in autophagy is unclear. Here, we report that polyneddylation functions as a post-translational modification for autophagic degradation of proteotoxic-stress induced protein aggregates. We also show that HYPK functions as an autophagy receptor in the polyneddylation-dependent aggrephagy. The scaffolding function of HYPK is facilitated by its C-terminal ubiquitin-associated domain and N-terminal tyrosine-type LC3-interacting region which bind to NEDD8 and LC3 respectively. Both NEDD8 and HYPK are positive modulators of basal and proteotoxicity-induced autophagy, leading to protection of cells from protein aggregates, such as aggregates of mutant HTT exon 1. Thus, we propose an indispensable and additive role of neddylation and HYPK in clearance of protein aggregates by autophagy, resulting in cytoprotective effect during proteotoxic stress.Abbreviations: ATG5, autophagy related 5; ATG12, autophagy related 12; ATG14, autophagy related 14; BECN1, beclin 1; CBL, casitas B-lineage lymphoma; CBLB, Cbl proto-oncogene B; GABARAP, GABA type A receptor-associated protein; GABARAPL1, GABA type A receptor associated protein like 1; GABARAPL2, GABA type A receptor associated protein like 2; GFP, green fluorescent protein; HTT, huntingtin; HTT97Q exon 1, huntingtin 97-glutamine exon 1; HUWE1, HECT, UBA and WWE domain containing E3 ubiquitin protein ligase 1; HYPK, huntingtin interacting protein K; IgG, immunoglobulin G; IMR-32, Institute for Medical Research-32; KD, knockdown; K_d, dissociation constant; LAMP1, lysosomal associated membrane protein 1; LIR, LC3 interacting region; MAP1LC3/LC3, microtubule associated protein 1 light chain 3; MAP1LC3A/LC3A, microtubule associated protein 1 light chain 3 alpha; MAP1LC3B/LC3B, microtubule associated protein 1 light chain 3 beta; MARK1, microtubule affinity regulating kinase 1; MARK2, microtubule affinity regulating kinase 2; MARK3, microtubule affinity regulating kinase 3; MARK4, microtubule affinity regulating kinase 4; MCF7, Michigan Cancer Foundation-7; MTOR, mechanistic target of rapamycin kinase; NAE1, NEDD8 activating enzyme E1 subunit 1; NBR1, NBR1 autophagy cargo receptor; NEDD8, NEDD8 ubiquitin like modifier; Ni-NTA, nickel-nitrilotriacetic acid; NUB1, negative regulator of ubiquitin like proteins 1; PIK3C3, phosphatidylinositol 3-kinase catalytic subunit type 3; PolyQ, poly-glutamine; PSMD8, proteasome 26S subunit, non-ATPase 8; RAD23A, RAD23 homolog A, nucleotide excision repair protein; RAD23B, RAD23 homolog B, nucleotide excision repair protein; RFP, red fluorescent protein; RPS27A, ribosomal protein S27a; RSC1A1, regulator of solute carriers 1; SNCA, synuclein alpha; SIK1, salt inducible kinase 1; siRNA, small interfering ribonucleic acid; SOD1, superoxide dismutase 1; SPR, surface plasmon resonance; SQSTM1, sequestosome 1; SUMO1, small ubiquitin like modifier 1; TAX1BP1, Tax1 binding protein 1; TDRD3, tudor domain containing 3; TNRC6C, trinucleotide repeat containing adaptor 6C; TOLLIP, toll interacting protein; TUBA, tubulin alpha; TUBB, tubulin beta class I; UBA, ubiquitin-associated; UBA1, ubiquitin like modifier activating enzyme 1; UBA5, ubiquitin like modifier activating enzyme 5; UBAC1, UBA domain containing 1; UBAC2, UBA domain containing 2; UBAP1, ubiquitin associated protein 1; UBAP2, ubiquitin associated protein 2; UBASH3B, ubiquitin associated and SH3 domain containing B; UBD/FAT10, ubiquitin D; UBE2K, ubiquitin conjugating enzyme E2 K; UBLs, ubiquitin-like proteins; UBL7, ubiquitin like 7; UBQLN1, ubiquilin 1; UBQLN2, ubiquilin 2; UBQLN3, ubiquilin 3; UBQLN4, ubiquilin 4; UBXN1, UBX domain protein 1; ULK1, unc-51 like autophagy activating kinase 1; URM1, ubiquitin related modifier 1; USP5, ubiquitin specific peptidase 5; USP13, ubiquitin specific peptidase 13; VPS13D, vacuolar protein sorting 13 homolog D.

Phenothiazine as novel human superoxide dismutase modulators: discovery, optimization, and biological evaluation

Superoxide dismutases (SODs) are regarded as important antioxidants for protecting cells against damage arising from oxidative stress. Much research is focused on finding new chemicals with an ability to boost human SOD activity. In the research described herein a structure-based approach was used to identify new human Cu-Zn superoxide dismutase (SOD1) modulators based on previously reported plasmodium falciparum iron SOD inhibitors using induced fit docking and molecular dynamic (MD) protocols. The compound with the highest docking binding energy was selected for further structure simplification followed by structural similarity and MD in order to find a new activator/inhibitor scaffold of the SOD1 enzyme. According to the docking survey of the mentioned series, 1,4-bis(3-(1,4,8-trichloro-10Hphenothiazin-10-yl) propyl) piperazine (DS88) was the top scoring compound interacting with the SOD1 active site channel. Following structure simplification and similarity search, the most promising scaffold which is closely related to the phenothiazine antipsychotic class, was identified. Compared with the normal blood SOD1 activity, the percent of O₂ production increased with trifluoperazine, while it decreased with the chlorpromazine. The molecular dynamic investigation shows that trifluoperazine exerts its SOD1 activating effect by stabilizing electrostatic loop while chlorpromazine employs SOD1 inhibition activity through repositioning of the electrostatic loop and increasing its distance from the catalytic metal site which diminished substrate specificity and catalytic activity of the SOD1 enzyme. The results identified the preferred region, orientation, and types of interaction for each activator or inhibitor compound.

Differential Stabilities of Mefloquine-Bound Human and Plasmodium falciparum Acyl-CoA-Binding Proteins

Toxic effects of pharmacological drugs restrict their robust application against human diseases. Although used as a drug in the combinatorial therapy to treat malaria, the use of mefloquine is not highly recommended because of its adverse effects in humans. Mefloquine inhibits the binding of acyl-CoAs to acyl-CoA-binding proteins of Plasmodium falciparum (PfACBPs) and human (hACBP). In this study, we have used molecular dynamics simulation and other computational approaches to investigate the differences of stabilities of mefloquine-PfACBP749 and mefloquine-hACBP complexes. The stability of mefloquine in the binding cavity of PfACBP749 is less than its stability in the binding pocket of hACBP. Although the essential tyrosine residues (tyrosine-30 and tyrosine-33 of PfACBP749 and tyrosine-29 and tyrosine-32 of hACBP) mediate the initial binding of mefloquine to the proteins by π-stacking interactions, additional temporally longer interactions between mefloquine and aspartate-22 and methionine-25 of hACBP result in stronger binding of mefloquine to hACBP. The higher fluctuation of mefloquine-binding residues of PfACBP749 contributes to the instability of mefloquine in the binding cavity of the protein. On the contrary, in the mefloquine-bound state, the stability of hACBP protein is less than the stability of PfACBP749. The helix-to-coil transition of the N-terminal hydrophobic region of hACBP has a destabilizing effect upon the protein's structure. This causes the induction of aggregation properties in the hACBP in the mefloquine-bound state. Taken together, we describe the mechanistic features that affect the differential dynamic stabilities of mefloquine-bound PfACBP749 and hACBP proteins.

RNA-protein interactions: Central players in coordination of regulatory networks

Changes in the abundance of protein and RNA molecules can impair the formation of complexes in the cell leading to toxicity and death. Here we exploit the information contained in protein, RNA and DNA interaction networks to provide a comprehensive view of the regulation layers controlling the concentration-dependent formation of assemblies in the cell. We present the emerging concept that RNAs can act as scaffolds to promote the formation ribonucleoprotein complexes and coordinate the post-transcriptional layer of gene regulation. We describe the structural and interaction network properties that characterize the ability of protein and RNA molecules to interact and phase separate in liquid-like compartments. Finally, we show that presence of structurally disordered regions in proteins correlate with the propensity to undergo liquid-to-solid phase transitions and cause human diseases. Also see the video abstract here https://youtu.be/kfpqibsNfS0.

The metastable states of proteins

The intriguing process of protein folding comprises discrete steps that stabilize the protein molecules in different conformations. The metastable state of protein is represented by specific conformational characteristics, which place the protein in a local free energy minimum state of the energy landscape. The native-to-metastable structural transitions are governed by transient or long-lived thermodynamic and kinetic fluctuations of the intrinsic interactions of the protein molecules. Depiction of the structural and functional properties of metastable proteins is not only required to understand the complexity of folding patterns but also to comprehend the mechanisms of anomalous aggregation of different proteins. In this article, we review the properties of metastable proteins in context of their stability and capability of undergoing atypical aggregation in physiological conditions.