Search

Single-point mutations in the transmembrane (TM) region of receptor tyrosine kinases (RTKs) can lead to abnormal ligand-independent activation. We use a combination of computational modeling, NMR spectroscopy and cell experiments to analyze in detail the mechanism of how TM domains contribute to the activation of wild-type (WT) PDGFRA and its oncogenic V536E mutant. Using a computational framework, we scan all positions in PDGFRA TM helix for identification of potential functional mutations for the WT and the mutant and reveal the relationship between the receptor activity and TM dimerization via different interfaces. This strategy also allows us design a novel activating mutation in the WT (I537D) and a compensatory mutation in the V536E background eliminating its constitutive activity (S541G). We show both computationally and experimentally that single-point mutations in the TM region reshape the TM dimer ensemble and delineate the structural and dynamic determinants of spontaneous activation of PDGFRA via its TM domain. Our atomistic picture of the coupling between TM dimerization and PDGFRA activation corroborates the data obtained for other RTKs and provides a foundation for developing novel modulators of the pathological activity of PDGFRA.

The epidermal growth factor receptor (EGFR) family is an important class of receptor tyrosine kinases, mediating a variety of cellular responses in normal biological processes and in pathological states of multicellular organisms. Different modes of dimerization of the human EGFR transmembrane domain (TMD) in different membrane mimetics recently prompted us to propose a novel signal transduction mechanism based on protein-lipid interaction. However, the experimental evidence for it was originally obtained with slightly different TMD fragments used in the two different mimetics, compromising the validity of the comparison. To eliminate ambiguity, we determined the nuclear magnetic resonance (NMR) structure of the bicelle-incorporated dimer of the EGFR TMD fragment identical to the one previously used in micelles. The NMR results augmented by molecular dynamics simulations confirm the mutual influence of the TMD and lipid environment, as is required for the proposed lipid-mediated activation mechanism. They also reveal the possible functional relevance of a subtle interplay between the concurrent processes in the lipid and protein during signal transduction.

PR homology domain-containing member 12 (PRDM12) belongs to a family of conserved transcription factors implicated in cell fate decisions. Here we show that PRDM12 is a key regulator of sensory neuronal specification in Xenopus. Modeling of human PRDM12 mutations that cause hereditary sensory and autonomic neuropathy (HSAN) revealed remarkable conservation of the mutated residues in evolution. Expression of wild-type human PRDM12 in Xenopus induced the expression of sensory neuronal markers, which was reduced using various human PRDM12 mutants. In Drosophila, we identified Hamlet as the functional PRDM12 homolog that controls nociceptive behavior in sensory neurons. Furthermore, expression analysis of human patient fibroblasts with PRDM12 mutations uncovered possible downstream target genes. Knockdown of several of these target genes including thyrotropin-releasing hormone degrading enzyme (TRHDE) in Drosophila sensory neurons resulted in altered cellular morphology and impaired nociception. These data show that PRDM12 and its functional fly homolog Hamlet are evolutionary conserved master regulators of sensory neuronal specification and play a critical role in pain perception. Our data also uncover novel pathways in multiple species that regulate evolutionary conserved nociception.

Recently, the ability to interact with messenger RNA (mRNA) has been reported for a number of known RNA-binding proteins, but surprisingly also for different proteins without recognizable RNA binding domains including several transcription factors and metabolic enzymes. Moreover, direct binding to cognate mRNAs has been detected for multiple proteins, thus creating a strong impetus to search for functional significance and basic physico-chemical principles behind such interactions. Here, we derive interaction preferences between amino acids and RNA bases by analyzing binding interfaces in the known 3D structures of protein–RNA complexes. By applying this tool to human proteome, we reveal statistically significant matching between the composition of mRNA sequences and base-binding preferences of protein sequences they code for. For example, purine density profiles of mRNA sequences mirror guanine affinity profiles of cognate protein sequences with quantitative accuracy (median Pearson correlation coefficient R = −0.80 across the entire human proteome). Notably, statistically significant anti-matching is seen only in the case of adenine. Our results provide strong evidence for the stereo-chemical foundation of the genetic code and suggest that mRNAs and cognate proteins may in general be directly complementary to each other and associate, especially if unstructured.

Mutations in transmembrane (TM) domains of receptor tyrosine kinases are shown to cause a number of inherited diseases and cancer development. Here, we use a combined molecular modeling approach to understand molecular mechanism of effect of G380R and A391E mutations on dimerization of TM domains of human fibroblast growth factor receptor 3 (FGFR3). According to results of Monte Carlo conformational search in the implicit membrane and further molecular dynamics simulations, TM dimer of this receptor is able to form a number of various conformations, which differ significantly by the free energy of association in a full-atom model bilayer. The aforementioned mutations affect dimerization efficiency of TM segments and lead to repopulation of conformational ensemble for the dimer. Particularly, both mutations do not change the dimerization free energy of the predominant (putative "non-active") symmetric conformation of TM dimer, while affect dimerization efficiency of its asymmetric ("intermediate") and alternative symmetric (putative "active") models. Results of our simulations provide novel atomistic prospective of the role of G380 and A391E mutations in dimerization of TM domains of FGFR3 and their consecutive contributions to the activation pathway of the receptor.

A central intermediate in purine catabolism, the inosine nucleobase hypoxanthine is also one of the most abundant modified nucleobases in RNA and plays key roles in the regulation of gene expression and determination of cell fate. It is known that hypoxanthine acts as guanine when interacting with other nucleobases and base pairs most favorably with cytosine. However, its preferences when it comes to interactions with amino acids remain unknown. Here we present for the first time the absolute binding free energies and the associated interaction modes between hypoxanthine and all standard, non-glycyl/non-prolyl amino acid side chain analogs as derived from molecular dynamics simulations and umbrella sampling in high- and low-dielectric environments. We illustrate the biological relevance of the derived affinities by providing a quantitative explanation for the specificity of hypoxanthine-guanine phosphoribosyltransferase, a key enzyme in the purine salvage pathway. Our results demonstrate that in its affinities for protein side chains, hypoxanthine closely matches guanine, much more so than its precursor adenine.

Association of transmembrane (TM) helices taking place in the cell membrane has an important contribution to the biological function of bitopic proteins, among which receptor tyrosine kinases represent a typical example and a potent target for medical applications. Since this process depends on a complex interplay of different factors (primary structures of TM domains and juxtamembrane regions, composition and phase of the local membrane environment, etc.), it is still far from being fully understood. Here, we present a computational modeling framework, which we have applied to systematically analyze dimerization of 18 TM helical homo- and heterodimers of different bitopic proteins, including the family of epidermal growth factor receptors (ErbBs). For this purpose, we have developed a novel surface-based modeling approach, which not only is able to predict particular conformations of TM dimers in good agreement with experiment but also provides screening of their conformational heterogeneity. Using all-atom molecular dynamics simulations of several of the predicted dimers in different model membranes, we have elucidated a putative role of the environment in selection of particular conformations. Simulation results clearly show that each particular bilayer preferentially stabilizes one of possible dimer conformations, and that the energy gain depends on the interplay between structural properties of the protein and the membrane. Moreover, the character of protein-driven perturbations of the bilayer is reflected in the contribution of a particular membrane to the free energy gain. We have found that the approximated dimerization strength for ErbBs family can be related to their oncogenic ability.

A potential connection between physico-chemical properties of mRNAs and cognate proteins, with implications concerning both the origin of the genetic code and mRNA–protein interactions, is unexplored. We compare pyrimidine content of naturally occurring mRNA coding sequences with the propensity of cognate protein sequences to interact with pyrimidines. The latter is captured by polar requirement, a measure of solubility of amino acids in aqueous solutions of pyridines, heterocycles closely related to pyrimidines. We find that the higher the pyrimidine content of an mRNA, the stronger the average propensity of its cognate protein’s amino acids to interact with pyridines. Moreover, window-averaged pyrimidine profiles of individual mRNAs strongly mirror polar-requirement profiles of cognate protein sequences. For example, 4953 human proteins exhibit a correlation between the two with |R| > 0.8. In other words, pyrimidine-rich mRNA regions quantitatively correspond to regions in cognate proteins containing residues soluble in pyrimidine mimetics and vice versa. Finally, by studying randomized genetic code variants we show that the universal genetic code is highly optimized to preserve these correlations. Overall, our findings redefine the stereo-chemical hypothesis concerning code’s origin and provide evidence of direct complementary interactions between mRNAs and cognate proteins before development of ribosomal decoding, but also presently, especially if both are unstructured.

We performed microsecond long coarse-grained molecular dynamics simulations to elucidate the lateral structure and domain dynamics of a phospha- tidylethanolamine (PE)/phosphatidylglycerol (PG) mixed bilayer (7/3), mimicking the inner membrane of gram-negative bacteria. Specifically, we address the effect of surface bound antimicrobial peptides (AMPs) on the lateral organization of the membrane.Wefindthat, intheabsenceofthe peptides,theminor PGfractiononly forms small clusters, but that these clusters grow in size upon binding of the cationic AMPs. The presence of AMPs systematically affects the dynamics and induces long-range order in the structure of PG domains, stabilizing the separation between the two lipid fractions. Our results help in understanding the initial stages of destabilization of cytoplasmic bacterial membranes below the critical peptide concentration necessary for disruption, and provide a possible explanation for the multimodal character of AMP activity. SECTION Surfactants, Membranes

Cotranscriptional ubiquitination of histone H2B is key to gene regulation. The yeast E3 ubiquitin ligase Bre1 (human RNF20/40) pairs with the E2 ubiquitin conjugating enzyme Rad6 to monoubiquitinate H2B at Lys123. How this single lysine residue on the nucleosome core particle (NCP) is targeted by the Rad6-Bre1 machinery is unknown. Using chemical cross-linking and mass spectrometry, we identified the functional interfaces of Rad6, Bre1, and NCPs in a defined in vitro system. The Bre1 RING domain cross-links exclusively with distinct regions of histone H2B and H2A, indicating a spatial alignment of Bre1 with the NCP acidic patch. By docking onto the NCP surface in this distinct orientation, Bre1 positions the Rad6 active site directly over H2B Lys123. The Spt-Ada-Gcn5 acetyltransferase (SAGA) H2B deubiquitinase module competes with Bre1 for binding to the NCP acidic patch, indicating regulatory control. Our study reveals a mechanism that ensures site-specific NCP ubiquitination and fine-tuning of opposing enzymatic activities.

The cell membrane is "stuffed" with proteins, whose transmembrane (TM) helical domains spontaneously associate to form functionally active complexes. For a number of membrane receptors, a modulation of TM domains' oligomerization has been shown to contribute to the development of severe pathological states, thus calling for detailed studies of the atomistic aspects of the process. Despite considerable progress achieved so far, several crucial questions still remain: How do the helices recognize each other in the membrane? What is the driving force of their association? Here, we assess the dimerization free energy of TM helices along with a careful consideration of the interplay between the structure and dynamics of protein and lipids using atomistic molecular dynamics simulations in the hydrated lipid bilayer for three different model systems - TM fragments of glycophorin A, polyalanine and polyleucine peptides. We observe that the membrane driven association of TM helices exhibits a prominent entropic character, which depends on the peptide sequence. Thus, a single TM peptide of a given composition induces strong and characteristic perturbations in the hydrophobic core of the bilayer, which may facilitate the initial "communication" between TM helices even at the distances of 20-30 Å. Upon tight helix-helix association, the immobilized lipids accommodate near the peripheral surfaces of the dimer, thus disturbing the packing of the surrounding. The dimerization free energy of the modeled peptides corresponds to the strength of their interactions with lipids inside the membrane being the lowest for glycophorin A and similarly higher for both homopolymers. We propose that the ability to accommodate lipid tails determines the dimerization strength of TM peptides and that the lipid matrix directly governs their association.

Early-stage evolutionary development of the universal genetic code remains a fundamental, open problem. One of the possible scenarios suggests that the code evolved in response to direct interactions between peptides and RNA oligonucleotides in the primordial environment. Recently, we have revealed a strong matching between base-binding preferences of modern protein sequences and the composition of their cognate mRNA coding sequences. These results point directly at the physicochemical foundation behind the code’s origin, but also support the possibility of direct complementary interactions between proteins and their cognate mRNAs, especially if the two are unstructured. Here, we analyze molecular-surface mapping of knowledge-based amino-acid/nucleobase interaction preferences for a set of complete, high-resolution protein structures and show that the connection between the two biopolymers could remain relevant even for structured, folded proteins. Specifically, protein surface loops are strongly enriched in residues with a high binding propensity for guanine and cytosine, while adenine- and uracil-preferring residues are uniformly distributed throughout protein structures. Moreover, compositional complementarity of cognate protein and mRNA sequences remains strong even after weighting protein sequence profiles by residue solvent exposure. Our results support the possibility that protein/mRNA sequence complementarity may also translate to cognate interactions between structured biopolymers.

To gain success in the evolutionary "arms race," venomous animals such as scorpions produce diverse neurotoxins selected to hit targets in the nervous system of prey. Scorpion α-toxins affect insect and/or mammalian voltage-gated sodium channels (Na(v)s) and thereby modify the excitability of muscle and nerve cells. Although more than 100 α-toxins are known and a number of them have been studied into detail, the molecular mechanism of their interaction with Na(v)s is still poorly understood. Here, we employ extensive molecular dynamics simulations and spatial mapping of hydrophobic/hydrophilic properties distributed over the molecular surface of α-toxins. It is revealed that despite the small size and relatively rigid structure, these toxins possess modular organization from structural, functional, and evolutionary perspectives. The more conserved and rigid "core module" is supplemented with the "specificity module" (SM) that is comparatively flexible and variable and determines the taxon (mammal versus insect) specificity of α-toxin activity. We further show that SMs in mammal toxins are more flexible and hydrophilic than in insect toxins. Concomitant sequence-based analysis of the extracellular loops of Na(v)s suggests that α-toxins recognize the channels using both modules. We propose that the core module binds to the voltage-sensing domain IV, whereas the more versatile SM interacts with the pore domain in repeat I of Na(v)s. These findings corroborate and expand the hypothesis on different functional epitopes of toxins that has been reported previously. In effect, we propose that the modular structure in toxins evolved to match the domain architecture of Na(v)s. ispartof: Journal of Biological Chemistry vol:288 issue:26 pages:19014-27 ispartof: location:United States status: published

Bacterial cell wall is targeted by many antibiotics. Among them are lantibiotics, which realize their function via interaction with plasma membrane lipid-II molecule — a chemically conserved part of the cell wall synthesis pathway. To investigate structural and dynamic properties of this molecule, we have performed a series of nearly microsecond-long molecular dynamics simulations of lipid-II and some of its analogs in zwitterionic single component and charged mixed simulated phospholipid bilayers (the reference and the mimic of the bacterial plasma membrane, respectively). Extensive analysis revealed that lipid-II forms a unique “amphiphilic pattern” exclusively on the surface of the simulated bacterial membrane (and not in the reference one). We hypothesize that many lantibiotics exploit the conserved features of lipid-II along with characteristic modulation of the bacterial membrane as the “landing site”. This putative recognition mechanism opens new opportunities for studies on lantibiotics action and design of novel armament against resistant bacterial strains.

Activating mutations in the platelet-derived growth factor (PDGF) receptor alpha (PDGFRA) have been described in patients with gastrointestinal stromal tumors or myeloid malignancies associated with hypereosinophilia. These patients respond well to imatinib mesylate, raising the question as to whether patients with a PDGF receptor mutation in other tumor types should receive a tyrosine kinase inhibitor treatment. We characterized 10 novel somatic point mutations in PDGFRA that have been reported in isolated cases of glioblastoma, melanoma, acute myeloid leukemia, peripheral nerve sheath tumors and neuroendocrine carcinoma. The PDGFRA transmembrane domain mutation V536E stimulated Ba/F3 cell growth and signaling via ERK and STAT5 in the absence of ligand. This mutant, identified in glioblastoma, was strongly inhibited by imatinib. Modeling suggested that the mutation modulates the packing of the transmembrane domain helices in the receptor dimer. By contrast, two mutations in highly conserved residues affected the receptor traffic to the cell surface or kinase activity, thereby preventing the response to PDGF. The other mutations had no significant impact on the receptor activity. This functional analysis matched the predictions of SIFT and PolyPhen for only five mutations and these algorithms do not discriminate gain from loss of function. Finally, an E996K variant that had been identified in a melanoma cell line was not expressed in these cells. Altogether, several newly identified PDGFRA mutations do not activate the receptor and may therefore represent passenger mutations. Our results also underline the importance of characterizing novel kinase alterations in cancer patients.

Influenza virus hemagglutinin is a homotrimeric spike glycoprotein crucial for virions' attachment, membrane fusion, and assembly reactions. X-ray crystallography data are available for hemagglutinin ectodomains of various types/subtypes but not for anchoring segments. To get structural information for the linker and transmembrane regions of hemagglutinin, influenza A (H1–H16 subtypes except H8 and H15) and B viruses were digested with bromelain or subtilisin Carlsberg, either within virions or in non-ionic detergent micelles. Proteolytical fragments were analyzed by sodium dodecyl sulfate–polyacrylamide gel electrophoresis and matrix-assisted laser desorption/ionization time-of-flight mass spectrometry. Within virions, hemagglutinins of most influenza A/Group-1 and type B virus strains were more susceptible to digestion with bromelain and/or subtilisin compared to A/Group-2 hemagglutinins. The cleavage sites were always located in the hemagglutinin linker sequence. In detergent, 1) bromelain cleaved hemagglutinin of every influenza A subtype in the linker region; 2) subtilisin cleaved Group-2 hemagglutinins in the linker region; 3) subtilisin cleaved Group-1 hemagglutinins in the transmembrane region; 4) both enzymes cleaved influenza B virus hemagglutinin in the transmembrane region. We propose that the A/Group-2 hemagglutinin linker and/or transmembrane regions are more tightly associated within trimers than type A/Group-1 and particularly type B ones. This hypothesis is underpinned by spatial trimeric structure modeling performed for transmembrane regions of both Group-1 and Group-2 hemagglutinin representatives. Differential S-acylation of the hemagglutinin C-terminal anchoring segment with palmitate/stearate residues possibly contributes to fine tuning of transmembrane trimer packing and stabilization since decreased stearate amount correlated with deeper digestion of influenza B and some A/Group-1 hemagglutinins.

In silico structural analyses of sets of α-helical antimicrobial peptides (AMPs) are performed. Differences between hemolytic and non-hemolytic AMPs are revealed in organization of their N-terminal region. A parameter related to hydrophobicity of the N-terminal part is proposed as a measure of the peptide propensity to exhibit hemolytic and other unwanted cytotoxic activities. Based on the information acquired, a rational approach for selective removal of these properties in AMPs is suggested. A proof of concept is gained through engineering specific mutations that resulted in elimination of the hemolytic activity of AMPs (latarcins) while leaving the beneficial antimicrobial effect intact.

Being related by the genetic code, messenger RNAs (mRNAs) and cognate proteins are polymers with mutually interdependent compositions, which further implies the possibility of a direct connection between their general physicochemical properties. How efficiently do different characteristics of mRNA coding regions reflect the features of cognate proteins and is it possible for the cell to obtain information about proteins from their mRNAs without first reading them on the ribosome? We address these issues in a theoretical proteome-wide analysis and show that average protein hydrophobicity, calculated from either sequences or 3D structures, can be encoded in an analogue fashion by many different mRNA sequence properties with the only constraint being that pyrimidine and purine bases be clearly distinguishable on average [1]. Moreover, average characteristics of mRNA sequences allow for a reasonable discrimination between human proteins with different cellular localization and, in particular, cytosolic and membrane proteins, even in the absence of topogenic signal-based mechanisms. We discuss our findings in the context of protein and mRNA localization and propose that this cellular process may be partly determined by basic physicochemical rationales and interdependencies between the two biomolecules.[1] Polyansky AA, Hlevnjak M, Zagrovic B. (2013) Nat. Commun. 4, 2784.