Molecular Forces Governing the Biological Function of Per-Arnt-Sim-B (PAS-B) Domains: A Comparative Computational Study

\ua9 2021 by the authors.Per-Arnt-Sim (PAS) domains are evolutionarily-conserved regions found in proteins in all living systems, involved in transcriptional regulation and the response to hypoxic and xenobiotic stress. Despite having low primary sequence similarity, they show an impressively high structural conservation. Nonetheless, understanding the underlying mechanisms that drive the biological function of the PAS domains remains elusive. In this work, we used molecular dynamics simulations and bioinformatics tools in order the investigate the molecular characteristics that govern the intrinsic dynamics of five PAS-B domains (human AhR receptor, NCOA1, HIF1α, and HIF2α transcription factors, and Drosophila Suzukii (D. Suzukii) juvenile hormone receptor JHR). First, we investigated the effects of different length of N and C terminal regions of the AhR PAS-B domain, showing that truncation of those segments directly affects structural stability and aggregation propensity of the domain. Secondly, using the recently annotated PAS-B located in the methoprene-tolerant protein/juvenile hormone receptor (JHR) from D. Suzukii, we have shown that the mutation of the highly conserved “gatekeeper” tyrosine to phenylalanine (Y322F) does not affect the stability of the domain. Finally, we investigated possible redox-regulation of the AhR PAS-B domain by focusing on the cysteinome residues within PAS-B domains. The cysteines in AhR PAS-B are directly regulating the dynamics of the small molecule ligand-gating loop (residues 305 to 326). In conclusion, we comprehensibly described several molecular features governing the behaviour of PAS-B domains in solution, which may lead to a better understanding of the forces driving their biological functions

Oxidation of SQSTM1/p62 mediates the link between redox state and protein homeostasis

Cellular homoeostatic pathways such as macroautophagy (hereinafter autophagy) are regulated by basic mechanisms that are conserved throughout the eukaryotic kingdom. However, it remains poorly understood how these mechanisms further evolved in higher organisms. Here we describe a modification in the autophagy pathway in vertebrates, which promotes its activity in response to oxidative stress. We have identified two oxidation-sensitive cysteine residues in a prototypic autophagy receptor SQSTM1/p62, which allow activation of pro-survival autophagy in stress conditions. The Drosophila p62 homologue, Ref(2)P, lacks these oxidation-sensitive cysteine residues and their introduction into the protein increases protein turnover and stress resistance of flies, whereas perturbation of p62 oxidation in humans may result in age-related pathology. We propose that the redox-sensitivity of p62 may have evolved in vertebrates as a mechanism that allows activation of autophagy in response to oxidative stress to maintain cellular homoeostasis and increase cell survival.Peer reviewe

Oxidation of SQSTM1/p62 mediates the link between redox state and protein homeostasis

Cellular homoeostatic pathways such as macroautophagy (hereinafter autophagy) are regulated by basic mechanisms that are conserved throughout the eukaryotic kingdom. However, it remains poorly understood how these mechanisms further evolved in higher organisms. Here we describe a modification in the autophagy pathway in vertebrates, which promotes its activity in response to oxidative stress. We have identified two oxidation-sensitive cysteine residues in a prototypic autophagy receptor SQSTM1/p62, which allow activation of pro-survival autophagy in stress conditions. The Drosophila p62 homologue, Ref(2)P, lacks these oxidation-sensitive cysteine residues and their introduction into the protein increases protein turnover and stress resistance of flies, whereas perturbation of p62 oxidation in humans may result in age-related pathology. We propose that the redox-sensitivity of p62 may have evolved in vertebrates as a mechanism that allows activation of autophagy in response to oxidative stress to maintain cellular homoeostasis and increase cell survival

Cosolvent Analysis Toolkit (CAT): a robust hotspot identification platform for cosolvent simulations of proteins to expand the druggable proteome

Cosolvent Molecular Dynamics (MD) simulations are increasingly popular techniques developed for prediction and characterisation of allosteric and cryptic binding sites, which can be rendered “druggable” by small molecule ligands. Despite their conceptual simplicity and effectiveness, the analysis of cosolvent MD trajectories relies on pocket volume data, which requires a high level of manual investigation and may introduce a bias. In this work, we present CAT (Cosolvent Analysis Toolkit): an open-source, freely accessible analytical tool, suitable for automated analysis of cosolvent MD trajectories. CAT is compatible with commonly used molecular graphics software packages such as UCSF Chimera and VMD. Using a novel hybrid empirical force field scoring function, CAT accurately ranks the dynamic interactions between the macromolecular target and cosolvent molecules. To benchmark, CAT was used for three validated protein targets with allosteric and orthosteric binding sites, using five chemically distinct cosolvent molecules. For all systems, CAT has accurately identified all known sites. CAT can thus assist in computational studies aiming at identification of protein “hotspots” in a wide range of systems. As an easy-to-use computational tool, we expect that CAT will contribute to an increase of the size of the potentially ‘druggable’ human proteome

Structure, dynamics, and ligand recognition of human-specific chrfam7a (Dupα7) nicotinic receptor linked to neuropsychiatric disorders

Cholinergic α7 nicotinic receptors encoded by the CHRNA7 gene are ligand-gated ion channels directly related to memory and immunomodulation. Exons 5-7 in CHRNA7 can be duplicated and fused to exons A-E of FAR7a, resulting in a hybrid gene known as CHRFAM7A, unique to humans. Its product, denoted herein as Dupα7, is a truncated subunit where the N-terminal 146 residues of the ligand binding domain of the α7 receptor have been replaced by 27 residues from FAM7. Dupα7 has a negative effect on the functioning of α7 receptors associated with neurological disorders, including Alzheimer\u27s diseases and schizophrenia. However, the stoichiometry for the α7 nicotinic receptor containing dupα7 monomers remains unknown. In this work, we developed computational models of all possible combinations of wild-type α7 and dupα7 pentamers and evaluated their stability via atomistic molecular dynamics and coarse-grain simulations. We assessed the effect of dupα7 subunits on the Ca2+ conductance using free energy calculations. We showed that receptors comprising of four or more dupα7 subunits are not stable enough to constitute a functional ion channel. We also showed that models with dupα7/ α7 interfaces are more stable and are less detrimental for the ion conductance in comparison to dupα7/ dupα7 interfaces. Based on these models, we used protein-protein docking to evaluate how such interfaces would interact with an antagonist, α-bungarotoxin, and amyloid Aβ42. Our findings show that the optimal stoichiometry of dupα7/ α7 functional pentamers should be no more than three dupα7 monomers, in favour of a dupα7/α7 interface in comparison to a homodimer dupα7/dupα7 interface. We also showed that receptors bearing dupα7 subunits are less sensitive to Aβ42 effects, which may shed light on the translational gap reported for strategies focused on nicotinic receptors in Alzheimer\u27s Disease research. </p