Thermodynamic coupling of the tandem RRM domains of hnRNP A1 underlie its pleiotropic RNA binding functions.

The functional properties of RNA binding proteins (RBPs) require allosteric regulation through interdomain communication. Despite the importance of allostery to biological regulation, only a few studies have been conducted to describe the biophysical nature by which interdomain communication manifests in RBPs. Here, we show for hnRNP A1 that interdomain communication is vital for the unique stability of its amino-terminal domain, which consists of two RNA recognition motifs (RRMs). These RRMs exhibit drastically different stability under pressure. RRM2 unfolds as an individual domain but remains stable when appended to RRM1. Variants that disrupt interdomain communication between the tandem RRMs show a significant decrease in stability. Carrying these mutations over to the full-length protein for in vivo experiments revealed that the mutations affected the ability of the disordered carboxyl-terminal domain to engage in protein-protein interactions and influenced the protein's RNA binding capacity. Collectively, this work reveals that thermodynamic coupling between the tandem RRMs of hnRNP A1 accounts for its allosteric regulatory functions.

APIS: an updated parentage assignment software managing triploids induced from diploid parents.

In aquaculture, sterile triploids are commonly used for production as sterility gives them potential gains in growth, yields, and quality. However, they cannot be reproduced, and DNA parentage assignment to their diploid or tetraploid parents is required to estimate breeding values for triploid phenotypes. No publicly available software has the ability to assign triploids to their parents. Here, we updated the R package APIS to support triploids induced from diploid parents. First, we created new exclusion and likelihood tables that account for the double allelic contribution of the dam and the recombination that can occur during female meiosis. As the effective recombination rate of each marker with the centromere is usually unknown, we set it at 0.5 and found that this value maximizes the assignment rate even for markers with high or low recombination rates. The number of markers needed for a high true assignment rate did not strongly depend on the proportion of missing parental genotypes. The assignment power was however affected by the quality of the markers (minor allele frequency, call rate). Altogether, 96-192 SNPs were required to have a high parentage assignment rate in a real rainbow trout dataset of 1,232 triploid progenies from 288 parents. The likelihood approach was more efficient than exclusion when the power of the marker set was limiting. When more markers were used, exclusion was more advantageous, with sensitivity reaching unity, very low false discovery rate (<0.01), and excellent specificity (0.96-0.99). Thus, APIS provides an efficient solution to assign triploids to their diploid parents.

Elucidation of the Mechanisms of Inter-domain Coupling in the Monomeric State of Enzyme I by High-pressure NMR.

The catalytic cycle of Enzyme I (EI), a phosphotransferase enzyme responsible for converting phosphoenolpyruvate (PEP) into pyruvate, is characterized by a series of local and global conformational rearrangements. This multistep process includes a monomer-to-dimer transition, followed by an open-to-closed rearrangement of the dimeric complex upon PEP binding. In the present study, we investigate the thermodynamics of EI dimerization using a range of high-pressure solution NMR techniques complemented by SAXS experiments. 1H-15N TROSY and 1H-13C methyl TROSY NMR spectra combined with 15N relaxation measurements revealed that a native-like engineered variant of full-length EI fully dissociates into stable monomeric state above 1.5 kbar. Conformational ensembles of EI monomeric state were generated via a recently developed protocol combining coarse-grained molecular simulations with experimental backbone residual dipolar coupling measurements. Analysis of the structural ensembles provided detailed insights into the molecular mechanisms driving formation of the catalytically competent dimeric state, and reveals that each step of EI catalytical cycle is associated with a significant reduction in either inter- or intra-domain conformational entropy. Altogether, this study completes a large body work conducted by our group on EI and establishes a comprehensive structural and dynamical description of the catalytic cycle of this prototypical multidomain, oligomeric enzyme.

The intrinsically disordered cytoplasmic tail of a dendrite branching receptor uses two distinct mechanisms to regulate the actin cytoskeleton.

Dendrite morphogenesis is essential for neural circuit formation, yet the molecular mechanisms underlying complex dendrite branching remain elusive. Previous studies on the highly branched Caenorhabditis elegans PVD sensory neuron identified a membrane co-receptor complex that links extracellular signals to intracellular actin remodeling machinery, promoting high-order dendrite branching. In this complex, the claudin-like transmembrane protein HPO-30 recruits the WAVE regulatory complex (WRC) to dendrite branching sites, stimulating the Arp2/3 complex to polymerize actin. We report here our biochemical and structural analysis of this interaction, revealing that the intracellular domain (ICD) of HPO-30 is intrinsically disordered and employs two distinct mechanisms to regulate the actin cytoskeleton. First, HPO-30 ICD binding to the WRC requires dimerization and involves the entire ICD sequence, rather than a short linear peptide motif. This interaction enhances WRC activation by the GTPase Rac1. Second, HPO-30 ICD directly binds to the sides and barbed end of actin filaments. Binding to the barbed end requires ICD dimerization and inhibits both actin polymerization and depolymerization, resembling the actin capping protein CapZ. These dual functions provide an intriguing model of how membrane proteins can integrate distinct mechanisms to fine-tune local actin dynamics.

Thermodynamic Coupling of the tandem RRM domains of hnRNP A1 underlie its Pleiotropic RNA Binding Functions.

The functional properties of RNA-binding proteins (RBPs) require allosteric regulation through inter-domain communication. Despite the foundational importance of allostery to biological regulation, almost no studies have been conducted to describe the biophysical nature by which inter-domain communication manifests in RBPs. Here, we show through high-pressure studies with hnRNP A1 that inter-domain communication is vital for the unique stability of its N- terminal domain containing a tandem of RNA Recognition Motifs (RRMs). Despite high sequence similarity and nearly identical tertiary structures, the two RRMs exhibit drastically different stability under pressure. RRM2 unfolds completely under high-pressure as an individual domain, but when appended to RRM1, it remains stable. Variants in which inter-domain communication is disrupted between the tandem RRMs show a large decrease in stability under pressure. Carrying these mutations over to the full-length protein for in vivo experiments revealed that the mutations affected the ability of the disordered C-terminus to engage in protein-protein interactions and more importantly, they also influenced the RNA binding capacity. Collectively, this work reveals that thermodynamic coupling between the tandem RRMs of hnRNP A1 accounts for its allosteric regulatory functions.

Evidence of Disrupted-in Schizophrenia 1 (DISC1) as an arsenic binding protein and implications regarding its role as a translational activator.

Disrupted-in-schizophrenia-1 (DISC1) is a scaffold protein that plays a pivotal role in orchestrating signaling pathways involved in neurodevelopment, neural migration, and synaptogenesis. Among those, it has recently been reported that the role DISC1 in the Akt/mTOR pathway can shift from a global translational repressor to a translational activator in response to oxidative stress induced by arsenic. In this study we are providing evidence that DISC1 can directly bind arsenic via a C-terminal cysteine motif (C-X-C-X-C). A series of fluorescence-based binding assays were conducted with a truncated C-terminal domain construct of DISC1 and a of series of single, double, and triple cysteine mutants. We found that arsenous acid, a trivalent arsenic derivative, specifically binds to the C-terminal cysteine motif of DISC1 with low micromolar affinity. All three cysteines of the motif are required for high-affinity binding. Electron microscopy experiments combined with in silico structural predictions revealed that that the C-terminal of DISC1 forms an elongated tetrameric complex. The cysteine motif is consistently predicted to be located within a loop, fully exposed to solvent, providing a simple molecular framework to explain the high-affinity of DISC1 toward arsenous acid. This study sheds light on a novel functional facet of DISC1 as an arsenic binding protein and highlights its potential role as both a sensor and translational modulator within the Akt/mTOR pathway.

Evidence of DISC1 as an arsenic binding protein and implications regarding its role as a translational activator.

Disrupted-in-schizophrenia-1 (DISC1) is a scaffolding protein that plays a pivotal role in orchestrating signaling pathways involved in neurodevelopment, neural migration, and synaptogenesis. Among those, it has recently been reported that the role of DISC1 in the Akt/mTOR pathway can shift from a global translational repressor to a translational activator in response to oxidative stress induced by arsenic. In this study we provide evidence that DISC1 can directly bind arsenic via a C-terminal cysteine motif (C-X-C-X-C). A series of fluorescence-based binding assays were conducted with a truncated C-terminal domain construct of DISC1 and a series of single, double, and triple cysteine mutants. We found that arsenous acid, a trivalent arsenic derivative, specifically binds to the C-terminal cysteine motif of DISC1 with low micromolar affinity. All three cysteines of the motif are required for high-affinity binding. Electron microscopy experiments combined with in silico structural predictions reveal that the C-terminal of DISC1 forms an elongated tetrameric complex. The cysteine motif is consistently predicted to be located within a loop, fully exposed to solvent, providing a simple molecular framework to explain the high-affinity of DISC1 toward arsenous acid. This study sheds light on a novel functional facet of DISC1 as an arsenic binding protein and highlights its potential role as both a sensor and translational modulator within Akt/mTOR pathway.

Intricate coupling between the transactivation and basic-leucine zipper domains governs phosphorylation of transcription factor ATF4 by casein kinase 2.

Most transcription factors possess at least one long intrinsically disordered transactivation domain that binds to a variety of coactivators and corepressors and plays a key role in modulating the transcriptional activity. Despite the crucial importance of these domains, the structural and functional basis of transactivation remains poorly understood. Here, we focused on activating transcription factor 4 (ATF4)/cAMP response element-binding protein-2, an essential transcription factor for cellular stress adaptation. Bioinformatic sequence analysis of the ATF4 transactivation domain sequence revealed that the first 125 amino acids have noticeably less propensity for structural disorder than the rest of the domain. Using solution nuclear magnetic resonance spectroscopy complemented by a range of biophysical methods, we found that the isolated transactivation domain is predominantly yet not fully disordered in solution. We also observed that a short motif at the N-terminus of the transactivation domain has a high helical propensity. Importantly, we found that the N-terminal region of the transactivation domain is involved in transient long-range interactions with the basic-leucine zipper domain involved in DNA binding. Finally, in vitro phosphorylation assays with the casein kinase 2 show that the presence of the basic-leucine zipper domain is required for phosphorylation of the transactivation domain. This study uncovers the intricate coupling existing between the transactivation and basic-leucine zipper domains of ATF4, highlighting its potential regulatory significance.

MbnC Is Not Required for the Formation of the N-Terminal Oxazolone in the Methanobactin from Methylosinus trichosporium OB3b.

Methanobactins (MBs) are ribosomally synthesized and posttranslationally modified peptides (RiPPs) produced by methanotrophs for copper uptake. The posttranslational modification that defines MBs is the formation of two heterocyclic groups with associated thioamines from X-Cys dipeptide sequences. Both heterocyclic groups in the MB from Methylosinus trichosporium OB3b (MB-OB3b) are oxazolone groups. The precursor gene for MB-OB3b is mbnA, which is part of a gene cluster that contains both annotated and unannotated genes. One of those unannotated genes, mbnC, is found in all MB operons and, in conjunction with mbnB, is reported to be involved in the formation of both heterocyclic groups in all MBs. To determine the function of mbnC, a deletion mutation was constructed in M. trichosporium OB3b, and the MB produced from the ΔmbnC mutant was purified and structurally characterized by UV-visible absorption spectroscopy, mass spectrometry, and solution nuclear magnetic resonance (NMR) spectroscopy. MB-OB3b from the ΔmbnC mutant was missing the C-terminal Met and was also found to contain a Pro and a Cys in place of the pyrrolidinyl-oxazolone-thioamide group. These results demonstrate MbnC is required for the formation of the C-terminal pyrrolidinyl-oxazolone-thioamide group from the Pro-Cys dipeptide, but not for the formation of the N-terminal 3-methylbutanol-oxazolone-thioamide group from the N-terminal dipeptide Leu-Cys. IMPORTANCE A number of environmental and medical applications have been proposed for MBs, including bioremediation of toxic metals and nanoparticle formation, as well as the treatment of copper- and iron-related diseases. However, before MBs can be modified and optimized for any specific application, the biosynthetic pathway for MB production must be defined. The discovery that mbnC is involved in the formation of the C-terminal oxazolone group with associated thioamide but not for the formation of the N-terminal oxazolone group with associated thioamide in M. trichosporium OB3b suggests the enzymes responsible for posttranslational modification(s) of the two oxazolone groups are not identical.

High-Pressure NMR Experiments for Detecting Protein Low-Lying Conformational States.

High-pressure is a well-known perturbation method that can be used to destabilize globular proteins and dissociate protein complexes in a reversible manner. Hydrostatic pressure drives thermodynamical equilibria toward the state(s) with the lower molar volume. Increasing pressure offers, therefore, the opportunities to finely tune the stability of globular proteins and the oligomerization equilibria of protein complexes. High-pressure NMR experiments allow a detailed characterization of the factors governing the stability of globular proteins, their folding mechanisms, and oligomerization mechanisms by combining the fine stability tuning ability of pressure perturbation and the site resolution offered by solution NMR spectroscopy. Here we present a protocol to probe the local folding stability of a protein via a set of 2D 1H-15N experiments recorded from 1 bar to 2.5 kbar. The steps required for the acquisition and analysis of such experiments are illustrated with data acquired on the RRM2 domain of hnRNPA1.

Structure elucidation of the elusive Enzyme I monomer reveals the molecular mechanisms linking oligomerization and enzymatic activity.

Enzyme I (EI) is a phosphotransferase enzyme responsible for converting phosphoenolpyruvate (PEP) into pyruvate. This reaction initiates a five-step phosphorylation cascade in the bacterial phosphotransferase (PTS) transduction pathway. Under physiological conditions, EI exists in an equilibrium between a functional dimer and an inactive monomer. The monomer-dimer equilibrium is a crucial factor regulating EI activity and the phosphorylation state of the overall PTS. Experimental studies of EI's monomeric state have yet been hampered by the dimer's high thermodynamic stability, which prevents its characterization by standard structural techniques. In this study, we modified the dimerization domain of EI (EIC) by mutating three amino acids involved in the formation of intersubunit salt bridges. The engineered variant forms an active dimer in solution that can bind and hydrolyze PEP. Using hydrostatic pressure as an additional perturbation, we were then able to study the complete dissociation of the variant from 1 bar to 2.5 kbar in the absence and the presence of EI natural ligands. Backbone residual dipolar couplings collected under high-pressure conditions allowed us to determine the conformational ensemble of the isolated EIC monomeric state in solution. Our calculations reveal that three catalytic loops near the dimerization interface become unstructured upon monomerization, preventing the monomeric enzyme from binding its natural substrate. This study provides an atomic-level characterization of EI's monomeric state and highlights the role of the catalytic loops as allosteric connectors controlling both the activity and oligomerization of the enzyme.

Thermodynamic stability of hnRNP A1 low complexity domain revealed by high-pressure NMR.

We have investigated the pressure- and temperature-induced conformational changes associated with the low complexity domain of hnRNP A1, an RNA-binding protein able to phase separate in response to cellular stress. Solution NMR spectra of the hnRNP A1 low-complexity domain fused with protein-G B1 domain were collected from 1 to 2500 bar and from 268 to 290 K. While the GB1 domain shows the typical pressure-induced and cold temperature-induced unfolding expected for small globular domains, the low-complexity domain of hnRNP A1 exhibits unusual pressure and temperature dependences. We observed that the low-complexity domain is pressure sensitive, undergoing a major conformational transition within the prescribed pressure range. Remarkably, this transition has the inverse temperature dependence of a typical folding-unfolding transition. Our results suggest the presence of a low-lying extended and fully solvated state(s) of the low-complexity domain that may play a role in phase separation. This study highlights the exquisite sensitivity of solution NMR spectroscopy to observe subtle conformational changes and illustrates how pressure perturbation can be used to determine the properties of metastable conformational ensembles.

Disorder Mediated Oligomerization of DISC1 Proteins Revealed by Coarse-Grained Molecular Dynamics Simulations.

Disrupted-in-schizophrenia-1 (DISC1) is a scaffold protein of significant importance for neuro-development and a prominent candidate protein in the etiology of mental disorders. In this work, we investigate the role of conformational heterogeneity and local structural disorder in the oligomerization pathway of the full-length DISC1 and of two truncation variants. Through extensive coarse-grained molecular dynamics simulations with a predictive energy landscape-based model, we shed light on the interplay of local and global disorder which lead to different oligomerization pathways. We found that both global conformational heterogeneity and local structural disorder play an important role in shaping distinct oligomerization trends of DISC1. This study also sheds light on the differences in oligomerization pathways of the full-length protein compared to the truncated variants produced by a chromosomal translocation associated with schizophrenia. We report that oligomerization of full-length DISC1 sequence works in a nonadditive manner with respect to truncated fragments that do not mirror the conformational landscape or binding affinities of the full-length unit.

Local unfolding of the HSP27 monomer regulates chaperone activity.

The small heat-shock protein HSP27 is a redox-sensitive molecular chaperone that is expressed throughout the human body. Here, we describe redox-induced changes to the structure, dynamics, and function of HSP27 and its conserved α-crystallin domain (ACD). While HSP27 assembles into oligomers, we show that the monomers formed upon reduction are highly active chaperones in vitro, but are susceptible to self-aggregation. By using relaxation dispersion and high-pressure nuclear magnetic resonance (NMR) spectroscopy, we observe that the pair of ß-strands that mediate dimerisation partially unfold in the monomer. We note that numerous HSP27 mutations associated with inherited neuropathies cluster to this dynamic region. High levels of sequence conservation in ACDs from mammalian sHSPs suggest that the exposed, disordered interface present in free monomers or oligomeric subunits may be a general, functional feature of sHSPs.

Exploring Protein Conformational Landscapes Using High-Pressure NMR.

Protein conformational landscapes define their functional properties as well as their proteostasis. Hence, detailed mapping of these landscapes is necessary to understand and modulate protein conformation. The combination of high pressure and NMR provides a particularly powerful approach to characterizing protein conformational transitions. First, pressure, because its effects on protein structure arise from elimination of solvent excluded void volume, represents a more subtle perturbation than chemical denaturants, favoring the population of intermediates. Second, the residue-specific and multifaceted nature of NMR observables informs on many local structural properties of proteins, aiding in the characterization of intermediate and excited states.

Lessons from pressure denaturation of proteins.

Although it is now relatively well understood how sequence defines and impacts global protein stability in specific structural contexts, the question of how sequence modulates the configurational landscape of proteins remains to be defined. Protein configurational equilibria are generally characterized by using various chemical denaturants or by changing temperature or pH. Another thermodynamic parameter which is less often used in such studies is high hydrostatic pressure. This review discusses the basis for pressure effects on protein structure and stability, and describes how the unique mechanisms of pressure-induced unfolding can provide unique insights into protein conformational landscapes.

Pushing the Limits of Structure-Based Models: Prediction of Nonglobular Protein Folding and Fibrils Formation with Go-Model Simulations.

The development of computational efficient models is essential to obtain a detailed characterization of the mechanisms underlying the folding of proteins and the formation of amyloid fibrils. Structure-based computational models (Go-model) with Cα or all-atom resolutions have been able to successfully delineate the mechanisms of folding of several globular proteins and offer an interesting alternative to computationally intensive simulations with explicit solvent description. Here, we explore the limits of Go-model predictions by analyzing the folding of the nonglobular repeat domain proteins Notch Ankyrin and p16INK4 and the formation of human islet amyloid polypeptide (hIAPP) fibrils. Folding trajectories of the repeat domain proteins revealed that an all-atom resolution is required to capture the folding pathways and cooperativity reported in experimental studies. The all-atom Go-model was also successful in predicting the free-energy landscape of hIAPP fibrillation, suggesting a "dock and lock" mechanism of fibril elongation. We finally explored how mutations can affect the co-assembly of hIAPP fibrils by simulating a heterogeneous system composed of wild-type and mutated hIAPP peptides. Overall, this study shows that all-atom Go-model-based simulations have the potential of discerning the effects of mutations and post-translational modifications in protein folding and association and may help in resolving the dichotomy between experimental and theoretical studies on protein folding and amyloid fibrillation.

Prediction of nearest neighbor effects on backbone torsion angles and NMR scalar coupling constants in disordered proteins.

Using fine-tuned hydrogen bonding criteria, a library of coiled peptide fragments has been generated from a large set of high-resolution protein X-ray structures. This library is shown to be an improved representation of ϕ/ψ torsion angles seen in intrinsically disordered proteins (IDPs). The ϕ/ψ torsion angle distribution of the library, on average, provides good agreement with experimentally observed chemical shifts and 3 JHN-Hα coupling constants for a set of five disordered proteins. Inspection of the coil library confirms that nearest-neighbor effects significantly impact the ϕ/ψ distribution of residues in the coil state. Importantly, 3 JHN-Hα coupling constants derived from the nearest-neighbor modulated backbone ϕ distribution in the coil library show improved agreement to experimental values, thereby providing a better way to predict 3 JHN-Hα coupling constants for IDPs, and for identifying locations that deviate from fully random behavior.

ScanFold: an approach for genome-wide discovery of local RNA structural elements-applications to Zika virus and HIV.

In addition to encoding RNA primary structures, genomes also encode RNA secondary and tertiary structures that play roles in gene regulation and, in the case of RNA viruses, genome replication. Methods for the identification of functional RNA structures in genomes typically rely on scanning analysis windows, where multiple partially-overlapping windows are used to predict RNA structures and folding metrics to deduce regions likely to form functional structure. Separate structural models are produced for each window, where the step size can greatly affect the returned model. This makes deducing unique local structures challenging, as the same nucleotides in each window can be alternatively base paired. We are presenting here a new approach where all base pairs from analysis windows are considered and weighted by favorable folding. This results in unique base pairing throughout the genome and the generation of local regions/structures that can be ranked by their propensity to form unusually thermodynamically stable folds. We applied this approach to the Zika virus (ZIKV) and HIV-1 genomes. ZIKV is linked to a variety of neurological ailments including microcephaly and Guillain-Barré syndrome and its (+)-sense RNA genome encodes two, previously described, functionally essential structured RNA regions. HIV, the cause of AIDS, contains multiple functional RNA motifs in its genome, which have been extensively studied. Our approach is able to successfully identify and model the structures of known functional motifs in both viruses, while also finding additional regions likely to form functional structures. All data have been archived at the RNAStructuromeDB (www.structurome.bb.iastate.edu), a repository of RNA folding data for humans and their pathogens.

Monitoring protein folding through high pressure NMR spectroscopy.

High-pressure is a well-known perturbation method used to destabilize globular proteins. It is perfectly reversible, which is essential for a proper thermodynamic characterization of a protein equilibrium. In contrast to other perturbation methods such as heat or chemical denaturant that destabilize protein structures uniformly, pressure exerts local effects on regions or domains of a protein containing internal cavities. When combined with NMR spectroscopy, hydrostatic pressure offers the possibility to monitor at a residue level the structural transitions occurring upon unfolding and to determine the kinetic properties of the process. High-pressure NMR experiments can now be routinely performed, owing to the recent development of commercially available high-pressure sample cells. This review summarizes recent advances and some future directions of high-pressure NMR techniques for the characterization at atomic resolution of the energy landscape of protein folding.