MITF is a novel transcriptional regulator of the calcium sensor STIM1: Significance in physiological melanogenesis.

Stromal Interaction Molecule1 (STIM1) is an endoplasmic reticulum membrane-localized calcium (Ca2+) sensor that plays a critical role in the store-operated Ca2+ entry (SOCE) pathway. STIM1 regulates a variety of physiological processes and contributes to a plethora of pathophysiological conditions. Several disease states and enhanced biological phenomena are associated with increased STIM1 levels and activity. However, molecular mechanisms driving STIM1 expression remain largely unappreciated. We recently reported that STIM1 expression augments during pigmentation. Nonetheless, the molecular choreography regulating STIM1 expression in melanocytes is completely unexplored. Here, we characterized the molecular events that regulate STIM1 expression during pigmentation. We demonstrate that physiological melanogenic stimuli α-melanocyte stimulating hormone (αMSH) increases STIM1 mRNA and protein levels. Further, αMSH stimulates STIM1 promoter-driven luciferase activity, thereby suggesting transcriptional upregulation of STIM1. We show that downstream of αMSH, microphthalmia-associated transcription factor (MITF) drives STIM1 expression. By performing knockdown and overexpression studies, we corroborated that MITF regulates STIM1 expression and SOCE. Next, we conducted extensive bioinformatics analysis and identified MITF-binding sites on the STIM1 promoter. We validated significance of the MITF-binding sites in controlling STIM1 expression by performing ChIP and luciferase assays with truncated STIM1 promoters. Moreover, we confirmed MITF's role in regulating STIM1 expression and SOCE in primary human melanocytes. Importantly, analysis of publicly available datasets substantiates a positive correlation between STIM1 and MITF expression in sun-exposed tanned human skin, thereby highlighting physiological relevance of this regulation. Taken together, we have identified a novel physiologically relevant molecular pathway that transcriptionally enhances STIM1 expression.

A rigorous framework for detecting SARS-CoV-2 spike protein mutational ensemble from genomic and structural features.

The recent release of SARS-CoV-2 genomic data from several countries has provided clues into the potential antigenic drift of the coronavirus population. In particular, the genomic instability observed in the spike protein necessitates immediate action and further exploration in the context of viral-host interactions. By temporally tracking 527,988 SARS-CoV-2 genomes, we identified invariant and hypervariable regions within the spike protein. We evaluated combination of mutations from SARS-CoV-2 lineages and found that maximum number of lineage-defining mutations were present in the N-terminal domain (NTD). Based on mutant 3D-structural models of known Variants of Concern (VOCs), we found that structural properties such as accessibility, secondary structural type, and intra-protein interactions at local mutation sites are greatly altered. Further, we observed significant differences between intra-protein networks of wild-type and Delta mutant, with the latter showing dense intra-protein contacts. Extensive molecular dynamics simulations of D614G mutant spike structure with hACE2 further revealed dynamic features with 47.7% of mutations mapping on flexible regions of spike protein. Thus, we propose that significant changes within spike protein structure have occurred that may impact SARS-CoV-2 pathogenesis, and repositioning of vaccine candidates is required to contain the spread of COVID-19 pathogen.

Integrated genomic view of SARS-CoV-2 in India.

Background: India first detected SARS-CoV-2, causal agent of COVID-19 in late January 2020, imported from Wuhan, China. From March 2020 onwards, the importation of cases from countries in the rest of the world followed by seeding of local transmission triggered further outbreaks in India. Methods: We used ARTIC protocol-based tiling amplicon sequencing of SARS-CoV-2 (n=104) from different states of India using a combination of MinION and MinIT sequencing from Oxford Nanopore Technology to understand how introduction and local transmission occurred. Results: The analyses revealed multiple introductions of SARS-CoV-2 genomes, including the A2a cluster from Europe and the USA, A3 cluster from Middle East and A4 cluster (haplotype redefined) from Southeast Asia (Indonesia, Thailand and Malaysia) and Central Asia (Kyrgyzstan). The local transmission and persistence of genomes A4, A2a and A3 was also observed in the studied locations. The most prevalent genomes with patterns of variance (confined in a cluster) remain unclassified, and are here proposed as A4-clade based on its divergence within the A cluster. Conclusions: The viral haplotypes may link their persistence to geo-climatic conditions and host response. Multipronged strategies including molecular surveillance based on real-time viral genomic data is of paramount importance for a timely management of the pandemic.

Human LC3 and GABARAP subfamily members achieve functional specificity via specific structural modulations.

Autophagy is a conserved adaptive cellular pathway essential to maintain a variety of physiological functions. Core components of this machinery are the six human Atg8 orthologs that initiate formation of appropriate protein complexes. While these proteins are routinely used as indicators of autophagic flux, it is presently not possible to discern their individual biological functions due to our inability to predict specific binding partners. In our attempts towards determining downstream effector functions, we developed a computational pipeline to define structural determinants of human Atg8 family members that dictate functional diversity. We found a clear evolutionary separation between human LC3 and GABARAP subfamilies and also defined a novel sequence motif responsible for their specificity. By analyzing known protein structures, we observed that functional modules or microclusters reveal a pattern of intramolecular network, including distinct hydrogen bonding of key residues (F52/Y49; a subset of HP2) that may directly modulate their interaction preferences. Multiple molecular dynamics simulations were performed to characterize how these proteins interact with a common protein binding partner, PLEKHM1. Our analysis showed remarkable differences in binding modes via intrinsic protein dynamics, with PLEKHM1-bound GABARAP complexes showing less fluctuations and higher number of contacts. We further mapped 373 genomic variations and demonstrated that distinct cancer-related mutations are likely to lead to significant structural changes. Our findings present a quantitative framework to establish factors underlying exquisite specificity of human Atg8 proteins, and thus facilitate the design of precise modulators.Abbreviations: Atg: autophagy-related; ECs: evolutionary constraints; GABARAP: GABA type A receptor-associated protein; HsAtg8: human Atg8; HP: hydrophobic pocket; KBTBD6: kelch repeat and BTB domain containing 6; LIR: LC3-interacting region; MAP1LC3/LC3: microtubule associated protein 1 light chain 3; MD: molecular dynamics; HIV-1 Nef: human immunodeficiency virus type 1 negative regulatory factor; PLEKHM1: pleckstrin homology and RUN domain containing M1; RMSD: root mean square deviation; SQSTM1/p62: sequestosome 1; WDFY3/ALFY: WD repeat and FYVE domain containing 3.

Revealing Conformational Transitions in G-Protein-Coupled Receptor Rhodopsin upon Phosphorylation.

G-protein-coupled receptors (GPCRs) have evolved as highly specialized cellular machinery that can dictate biological outcomes in response to diverse stimuli. Specifically, they induce multiple pathway responses upon structural perturbations induced at local protein sites. GPCRs utilize a concurrent strategy involving a central transmembrane topology and biochemical modifications for precise functional implementation. However, the specific role of the latter is not known due to the lack of precise probing techniques that can characterize receptor dynamics upon biochemical modifications. Phosphorylation is known to be one of the critical biochemical modifications in GPCRs that aids in receptor desensitization via arrestin binding. Here, we carry out all-atom molecular dynamics simulations of rhodopsin in a membrane environment to study its conformational dynamics induced upon phosphorylation. Interestingly, our comparative analysis of non-phosphorylated and phosphorylated rhodopsin structure demonstrated enhanced receptor stability upon phosphorylation at the C-terminal region that leads to the opening of the extracellular part of the transmembrane helices. In addition, monitoring the distinct number of phosphorylation states showed that having fewer phosphorylated residues does not bring about appropriate conformational changes in the extracellular region. Since phosphorylation results in receptor desensitization and recycling of the ligand, our findings provide significant insights into the conformational dynamics of the mechanism of ligand exit from the receptor.

A novel mutation alters the stability of PapA2 resulting in the complete abrogation of sulfolipids in clinical mycobacterial strains.

The analysis of whole genomes has revealed specific geographical distribution of Mycobacterium tuberculosis (Mtb) strains across the globe suggestive of unique niche dependent adaptive mechanisms. We provide an important correlation of a genome-based mutation to a molecular phenotype across two predominant clinical Mtb lineages of the Indian subcontinent. We have identified a distinct lineage specific mutation-G247C, translating into an alanine-proline conversion in the papA2 gene of Indo-oceanic lineage 1 (L1) Mtb strains, and restoration of cell wall sulfolipids by simple genetic complementation of papA2 from lineage 3 (L3) or from H37Rv (lineage 4-L4) attributed the loss of this glycolipid to this specific mutation in Indo-Oceanic L1 Mtb. The investigation of structure of Mtb PapA2 revealed a distinct nonribosomal peptide synthetase (NRPS) C domain conformation with an unconventional presence of a zinc binding motif. Surprisingly, the A83P mutation did not map to either the catalytic center in the N-terminal subdomain or any of the substrate-binding region of the protein. On the contrary, the inherent ability of mutant PapA2 to form insoluble aggregates and molecular simulations with the wild-type/mutant (Wt/mut) PapA2 purports an important role for the surface associated 83rd residue in protein conformation. This study demonstrates the importance of a critical structural residue in the papA2 protein of Mtb and helps establish a link between observed genomic alteration and its molecular consequence in the successful human pathogen Mtb. Significance We demonstrate the effect of a unique SNP in PapA2 gene of Indo-oceanic Mycobacterium tuberculosis (Mtb) strains leading to the loss of sulfolipid from these strains. By X-ray crystallographic analysis and molecular dynamics (MD) simulations, we show the importance of this residue in the global PapA2 structure. The presence of a Zn atom has not been reported before for this class of proteins. Here, we provide an important link between genomic alteration and its molecular consequence in Mtb highlighting one of the many adaptive mechanisms that have contributed to its success as a human pathogen. A high degree of identity with PapA1, 3, or 4 would help in interpreting the structure of these PapA proteins and other acyl transferases of other biological systems.

A conserved degron containing an amphipathic helix regulates the cholesterol-mediated turnover of human squalene monooxygenase, a rate-limiting enzyme in cholesterol synthesis.

Cholesterol biosynthesis in the endoplasmic reticulum (ER) is tightly controlled by multiple mechanisms to regulate cellular cholesterol levels. Squalene monooxygenase (SM) is the second rate-limiting enzyme in cholesterol biosynthesis and is regulated both transcriptionally and post-translationally. SM undergoes cholesterol-dependent proteasomal degradation when cholesterol is in excess. The first 100 amino acids of SM (designated SM N100) are necessary for this degradative process and represent the shortest cholesterol-regulated degron identified to date. However, the fundamental intrinsic characteristics of this degron remain unknown. In this study, we performed a series of deletions, point mutations, and domain swaps to identify a 12-residue region (residues Gln-62-Leu-73), required for SM cholesterol-mediated turnover. Molecular dynamics and circular dichroism revealed an amphipathic helix within this 12-residue region. Moreover, 70% of the variation in cholesterol regulation was dependent on the hydrophobicity of this region. Of note, the earliest known Doa10 yeast degron, Deg1, also contains an amphipathic helix and exhibits 42% amino acid similarity with SM N100. Mutating SM residues Phe-35/Ser-37/Leu-65/Ile-69 into alanine, based on the key residues in Deg1, blunted SM cholesterol-mediated turnover. Taken together, our results support a model whereby the amphipathic helix in SM N100 attaches reversibly to the ER membrane depending on cholesterol levels; with excess, the helix is ejected and unravels, exposing a hydrophobic patch, which then serves as a degradation signal. Our findings shed new light on the regulation of a key cholesterol synthesis enzyme, highlighting the conservation of critical degron features from yeast to humans.

Lipidated proteins: Spotlight on protein-membrane binding interfaces.

The covalent attachment of lipids to proteins is a fundamental property of all living cells. These lipidated or lipid-modified proteins are directly targeted to the membranes and display diverse functional roles that are critical to cell function such as membrane signaling and trafficking. All lipidated proteins have been classified by the type of chemical moieties that are attached to them mainly palmitoylation, myristoylation, prenylation, cholesterylation, and addition of the Glycosylphosphatidyl inositol (GPI) anchor. Although the distinct hydrophobic lipid moiety largely dictates the functional compartmentalization, it also facilitates membrane trafficking and triggers a wide range of signaling pathways in cellular growth. In this review, we will focus on mechanistic insights underlying their membrane attachment, with a view to understand the regions that contribute to protein-membrane interface specificity. We also present a computational case study to report how different membrane lipids modulate insertion of the lipidated LC3 protein. Finally, in this review, we demonstrate the potential of regulating lipid modifications that are essential for cell survival and we discuss the concepts regarding their role in disease and therapeutics.

Structural signatures of DRD4 mutants revealed using molecular dynamics simulations: Implications for drug targeting.

Human Dopamine Receptor D4 (DRD4) orchestrates several neurological functions and represents a target for many psychological disorders. Here, we examined two rare variants in DRD4; V194G and R237L, which elicit functional alterations leading to disruption of ligand binding and G protein coupling, respectively. Using atomistic molecular dynamics (MD) simulations, we provide in-depth analysis to reveal structural signatures of wild and mutant complexes with their bound agonist and antagonist ligands. We constructed intra-protein network graphs to discriminate the global conformational changes induced by mutations. The simulations also allowed us to elucidate the local side-chain dynamical variations in ligand-bound mutant receptors. The data suggest that the mutation in transmembrane V (V194G) drastically disrupts the organization of ligand binding site and causes disorder in the native helical arrangement. Interestingly, the R237L mutation leads to significant rewiring of side-chain contacts in the intracellular loop 3 (site of mutation) and also affects the distant transmembrane topology. Additionally, these mutations lead to compact ICL3 region compared to the wild type, indicating that the receptor would be inaccessible for G protein coupling. Our findings thus reveal unreported structural determinants of the mutated DRD4 receptor and provide a robust framework for design of effective novel drugs.

Structure and dynamics of DRD4 bound to an agonist and an antagonist using in silico approaches.

Human dopamine receptor D4 (DRD4), a member of G-protein coupled receptor (GPCR) family, plays a central role in cell signaling and trafficking. Dysfunctional activity of DRD4 can lead to several psychiatric conditions and, therefore, represents target for many neurological disorders. However, lack of atomic structure impairs our understanding of the mechanism regulating its activity. Here, we report the modeled structure of DRD4 alone and in complex with dopamine and spiperone, its natural agonist and antagonist, respectively. To assess the conformational dynamics induced upon ligand binding, all-atom explicit solvent molecular dynamics simulations in membrane environment were performed. Comprehensive analyses of simulations reveal that agonist binding triggers a series of conformational changes in the transmembrane region, including rearrangement of residues, characteristic of transmission and tyrosine toggle molecular switches. Further, the trajectories indicate that a loop region in the intracellular region--ICL3, is significantly dynamic in nature, mainly due to the side-chain movements of conserved proline residues involved in SH3 binding domains. Interestingly, in dopamine-bound receptor simulation, ICL3 represents an open conformation ideal for G protein binding. The structural and dynamical information presented here suggest a mode of activation of DRD4, upon ligand binding. Our study will help in further understanding of receptor activation, as acquiring structural information is crucial for the design of highly selective DRD4 ligands.

Human Dopamine Receptors Interaction Network (DRIN): a systems biology perspective on topology, stability and functionality of the network.

Dopamine receptors (DR) are one of the major neurotransmitter receptors present in human brain. Malfunctioning of these receptors is well established to trigger many neurological and psychiatric disorders. Taking into consideration that proteins function collectively in a network for most of the biological processes, the present study is aimed to depict the interactions between all dopamine receptors following a systems biology approach. To capture comprehensive interactions of candidate proteins associated with human dopamine receptors, we performed a protein-protein interaction network (PPIN) analysis of all five receptors and their protein partners by mapping them into human interactome and constructed a human Dopamine Receptors Interaction Network (DRIN). We explored the topology of dopamine receptors as molecular network, revealing their characteristics and the role of central network elements. More to the point, a sub-network analysis was done to determine major functional clusters in human DRIN that govern key neurological pathways. Besides, interacting proteins in a pathway were characterized and prioritized based on their affinity for utmost drug molecules. The vulnerability of different networks to the dysfunction of diverse combination of components was estimated under random and direct attack scenarios. To the best of our knowledge, the current study is unique to put all five dopamine receptors together in a common interaction network and to understand the functionality of interacting proteins collectively. Our study pinpointed distinctive topological and functional properties of human dopamine receptors that have helped in identifying potential therapeutic drug targets in the dopamine interaction network.

Pharmacophore modeling and virtual screening studies to design potential COMT inhibitors as new leads.

Catechol-O-methyltransferase (COMT) catalyzes the methylation of catecholamines, including neurotransmitters like dopamine, epinephrine and norepinephrine, leading to their degradation. COMT has been a subject of study for its implications in numerous neurological disorders like Parkinson's disease (PD), schizophrenia, and depression. The COMT gene is associated with many allelic variants, the Val108Met polymorphism being the most clinically significant. Availability of crystal structure of both 108V and 108M forms of human soluble-COMT (S-COMT) facilitated us to use structure-based virtual screening approach to obtain new hits by screening a library of CNS permeable compounds from ZINC database. In this study, E-pharmacophore was also used to generate pharmacophore models based on a series of known COMT inhibitors. A five-point pharmacophore model consisting of one hydrogen-bond acceptor (A), two hydrogen bond donors (D), and two aromatic rings (R) was generated for both the polymorphic forms of COMT. These models were then used for filtering ZINC-CNS permeable library to obtain new hits. Physicochemical properties were also calculated for all the hits obtained from both the approaches for favorable ADME properties. These identified hits maybe of interest for further structural optimization and biological evaluation assays.

Inhibitors of catechol-O-methyltransferase in the treatment of neurological disorders.

Catechol-O-methyltransferase (COMT) is the enzyme which catalyzes the transfer of a methyl group from S-adenosylmethionine to catechols and catecholamines, like the neurotransmitters dopamine, epinephrine and norepinephrine. COMT has implications in many neurological and psychiatric disorders like schizophrenia, Parkinson's disease (PD), bipolar disorders, etc. and therefore, it serves as an important drug target. Since its characterization in 1957, many inhibitors were designed where the first generation inhibitors were found to be highly toxic, short acting and had poor bioavailability. Currently, two of the second generation inhibitors, tolcapone and entacapone have been used for treatment of PD but are associated with various dopaminergic and gastro-intestinal side-effects. There have been several approaches for the design of novel COMT inhibitors with a good and safe therapeutic profile. The focus of this article is to review the current knowledge on COMT and the role of COMT inhibitors in the treatment of neurological disorders. The inhibitors have been classified into six different classes based on the structural framework. A historical overview of the discovery and development of COMT inhibitors is presented with a special emphasis on new generation of inhibitors till date.

Molecular modeling studies of Fatty acyl-CoA synthetase (FadD13) from Mycobacterium tuberculosis--a potential target for the development of antitubercular drugs.

Tuberculosis (TB) is a global health problem and the situation has become more precarious due to the advent of HIV infections and continuous rise in the number of multi-drug resistant strains of Mycobacterium tuberculosis (M. tb). Biochemical studies on Fatty Acyl-CoA Synthetases (FadD13), one of the gene products of mymA operon, have provided insights into the involvement of this protein in the activation of fatty acids. Due to non-availability of the crystal structure of FadD13, we have employed in silico approaches to resolve and characterize the structure of this important protein of M. tb. A three dimensional model of M. tb FadD13 was predicted by a de novo structure prediction server that integrates fragment assembly with SimFold energy function. With the aid of molecular mechanics and dynamics methods, the final model was obtained and assessed subsequently for global and local accuracy by various assessment programs. With this model, a flexible docking study with the substrates was performed. Results of ligand interactions with key amino acids in the binding site are also summarized. The molecular model for the M. tb FadD13 obtained sheds light on the topographical features of the binding pocket of the protein and provides atomic insight into the possible modes of substrate recognition. The three-dimensional model of FadD13 presented here would be helpful in guiding both enzymatic studies as well as design of specific inhibitors.