Optogenetic control of YAP cellular localisation and function.

YAP, an effector of the Hippo signalling pathway, promotes organ growth and regeneration. Prolonged YAP activation results in uncontrolled proliferation and cancer. Therefore, exogenous regulation of YAP activity has potential translational applications. We present a versatile optogenetic construct (optoYAP) for manipulating YAP localisation, and consequently its activity and function. We attach a LOV2 domain that photocages a nuclear localisation signal (NLS) to the N-terminus of YAP. In 488 nm light, the LOV2 domain unfolds, exposing the NLS, which shuttles optoYAP into the nucleus. Nuclear import of optoYAP is reversible and tuneable by light intensity. In cell culture, activated optoYAP promotes YAP target gene expression and cell proliferation. Similarly, optofYap can be used in zebrafish embryos to modulate target genes. We demonstrate that optoYAP can override a cell's response to substrate stiffness to generate anchorage-independent growth. OptoYAP is functional in both cell culture and in vivo, providing a powerful tool to address basic research questions and therapeutic applications in regeneration and disease.

The Hippo pathway effector Wwtr1 regulates cardiac wall maturation in zebrafish.

Cardiac trabeculation is a highly regulated process that starts with the delamination of compact layer cardiomyocytes. The Hippo signaling pathway has been implicated in cardiac development but many questions remain. We have investigated the role of Wwtr1, a nuclear effector of the Hippo pathway, in zebrafish and find that its loss leads to reduced cardiac trabeculation. However, in mosaic animals, wwtr1-/- cardiomyocytes contribute more frequently than wwtr1+/- cardiomyocytes to the trabecular layer of wild-type hearts. To investigate this paradox, we examined the myocardial wall at early stages and found that compact layer cardiomyocytes in wwtr1-/- hearts exhibit disorganized cortical actin structure and abnormal cell-cell junctions. Accordingly, wild-type cardiomyocytes in mosaic mutant hearts contribute less frequently to the trabecular layer than when present in mosaic wild-type hearts, indicating that wwtr1-/- hearts are not able to support trabeculation. We also found that Nrg/Erbb2 signaling, which is required for trabeculation, could promote Wwtr1 nuclear export in cardiomyocytes. Altogether, these data suggest that Wwtr1 establishes the compact wall architecture necessary for trabeculation, and that Nrg/Erbb2 signaling negatively regulates its nuclear localization and therefore its activity.

An autism-linked missense mutation in SHANK3 reveals the modularity of Shank3 function.

Genome sequencing has revealed an increasing number of genetic variations that are associated with neuropsychiatric disorders. Frequently, studies limit their focus to likely gene-disrupting mutations because they are relatively easy to interpret. Missense variants, instead, have often been undervalued. However, some missense variants can be informative for developing a more profound understanding of disease pathogenesis and ultimately targeted therapies. Here we present an example of this by studying a missense variant in a well-known autism spectrum disorder (ASD) causing gene SHANK3. We analyzed Shank3's in vivo phosphorylation profile and identified S685 as one phosphorylation site where one ASD-linked variant has been reported. Detailed analysis of this variant revealed a novel function of Shank3 in recruiting Abelson interactor 1 (ABI1) and the WAVE complex to the post-synaptic density (PSD), which is critical for synapse and dendritic spine development. This function was found to be independent of Shank3's other functions such as binding to GKAP and Homer. Introduction of this human ASD mutation into mice resulted in a small subset of phenotypes seen previously in constitutive Shank3 knockout mice, including increased allogrooming, increased social dominance, and reduced pup USV. Together, these findings demonstrate the modularity of Shank3 function in vivo. This modularity further indicates that there is more than one independent pathogenic pathway downstream of Shank3 and correcting a single downstream pathway is unlikely to be sufficient for clear clinical improvement. In addition, this study illustrates the value of deep biological analysis of select missense mutations in elucidating the pathogenesis of neuropsychiatric phenotypes.

Induction of interferon-stimulated genes and cellular stress pathways by morpholinos in zebrafish.

The phenotypes caused by morpholino-mediated interference of gene function in zebrafish are often not observed in the corresponding mutant(s). We took advantage of the availability of a relatively large collection of transcriptomic datasets to identify common signatures that characterize morpholino-injected animals (morphants). In addition to the previously reported activation of tp53 expression, we observed increased expression of the interferon-stimulated genes (ISGs), isg15 and isg20, the cell death pathway gene casp8, and other cellular stress response genes including phlda3, mdm2 and gadd45aa. Studies involving segmentation stage embryos were more likely to show upregulation of these genes. We also found that the expression of these genes could be upregulated by increasing doses of an egfl7 morpholino, or even high doses of the standard control morpholino. Thus, these data show that morpholinos can induce the expression of ISGs in zebrafish embryos and further our understanding of morpholino effects.

A new parameter-rich structure-aware mechanistic model for amino acid substitution during evolution.

Improvements in the description of amino acid substitution are required to develop better pseudo-energy-based protein structure-aware models for use in phylogenetic studies. These models are used to characterize the probabilities of amino acid substitution and enable better simulation of protein sequences over a phylogeny. A better characterization of amino acid substitution probabilities in turn enables numerous downstream applications, like detecting positive selection, ancestral sequence reconstruction, and evolutionarily-motivated protein engineering. Many existing Markov models for amino acid substitution in molecular evolution disregard molecular structure and describe the amino acid substitution process over longer evolutionary periods poorly. Here, we present a new model upgraded with a site-specific parameterization of pseudo-energy terms in a coarse-grained force field, which describes local heterogeneity in physical constraints on amino acid substitution better than a previous pseudo-energy-based model with minimum cost in runtime. The importance of each weight term parameterization in characterizing underlying features of the site, including contact number, solvent accessibility, and secondary structural elements was evaluated, returning both expected and biologically reasonable relationships between model parameters. This results in the acceptance of proposed amino acid substitutions that more closely resemble those observed site-specific frequencies in gene family alignments. The modular site-specific pseudo-energy function is made available for download through the following website: https://liberles.cst.temple.edu/Software/CASS/index.html.

Sequence, structure, and cooperativity in folding of elementary protein structural motifs.

Residue-level unfolding of two helix-turn-helix proteins--one naturally occurring and one de novo designed--is reconstructed from multiple sets of site-specific (13)C isotopically edited infrared (IR) and circular dichroism (CD) data using Ising-like statistical-mechanical models. Several model variants are parameterized to test the importance of sequence-specific interactions (approximated by Miyazawa-Jernigan statistical potentials), local structural flexibility (derived from the ensemble of NMR structures), interhelical hydrogen bonds, and native contacts separated by intervening disordered regions (through the Wako-Saitô-Muñoz-Eaton scheme, which disallows such configurations). The models are optimized by directly simulating experimental observables: CD ellipticity at 222 nm for model proteins and their fragments and (13)C-amide I' bands for multiple isotopologues of each protein. We find that data can be quantitatively reproduced by the model that allows two interacting segments flanking a disordered loop (double sequence approximation) and incorporates flexibility in the native contact maps, but neither sequence-specific interactions nor hydrogen bonds are required. The near-identical free energy profiles as a function of the global order parameter are consistent with expected similar folding kinetics for nearly identical structures. However, the predicted folding mechanism for the two motifs is different, reflecting the order of local stability. We introduce free energy profiles for "experimental" reaction coordinates--namely, the degree of local folding as sensed by site-specific (13)C-edited IR, which highlight folding heterogeneity and contrast its overall, average description with the detailed, local picture.

DNA-damage induced cell death in yap1;wwtr1 mutant epidermal basal cells.

In a previous study, it was reported that Yap1 and Wwtr1 in zebrafish regulates the morphogenesis of the posterior body and epidermal fin fold (Kimelman et al., 2017). We report here that DNA damage induces apoptosis of epidermal basal cells (EBCs) in zebrafish yap1-/-;wwtr1-/- embryos. Specifically, these mutant EBCs exhibit active Caspase-3, Caspase-8, and γH2AX, consistent with DNA damage serving as a stimulus of the extrinsic apoptotic pathway in epidermal cells. Live imaging of zebrafish epidermal cells reveals a steady growth of basal cell size in the developing embryo, but this growth is inhibited in mutant basal cells followed by apoptosis, leading to the hypothesis that factors underscoring cell size play a role in this DNA damage-induced apoptosis phenotype. We tested two of these factors using cell stretching and substrate stiffness assays, and found that HaCaT cells cultured on stiff substrates exhibit more numerous γH2AX foci compared to ones cultured on soft substrates. Thus, our experiments suggest that substrate rigidity may modulate genomic stress in epidermal cells, and that Yap1 and Wwtr1 promotes their survival.

Effect of Mutations on the Global and Site-Specific Stability and Folding of an Elementary Protein Structural Motif.

Understanding the folding mechanism of proteins requires detailed knowledge of the roles of individual amino acid residues in stabilization of specific elements and local segments of the native structure. Recently, we have utilized the combination of circular dichroism (CD) and site-specific 13C isotopically edited infrared spectroscopy (IR) coupled with the Ising-like model for protein folding to map the thermal unfolding at the residue level of a de novo designed helix-turn-helix motif αtα. Here we use the same methodology to study how the sequence of local thermal unfolding is affected by selected mutations introduced into the most and least stable parts of the motif. Seven different mutants of αtα are screened to find substitutions with the most pronounced effects on the overall stability. Subsequently, thermal unfolding of two mutated αtα sequences is studied with site-specific resolution, using four distinct 13C isotopologues of each. The data are analyzed with the Ising-like model, which builds on a previous parametrization for the original αtα sequence and tests different ways of incorporating the amino acid substitution. We show that for both more and less stable mutants only the adjustment of all interaction parameters of the model can yield a satisfactory fit to the experimental data. The stabilizing and destabilizing mutations result, respectively, in a similar increase and decrease of the stability of all probed local segments, irrespective of their position with respect to the mutation site. Consequently, the relative order of their unfolding remains essentially unchanged. These results underline the importance of the interconnectivity of the stabilizing interaction network and cooperativity of the protein structure, which is evident even in a small motif with apparently noncooperative, heterogeneous unfolding. Overall, our findings are consistent with the native structure being the dominant factor in determining the folding mechanism, regardless of the details of its overall or local thermodynamic stabilization.

Loss of pyruvate kinase M2 limits growth and triggers innate immune signaling in endothelial cells.

Despite their inherent proximity to circulating oxygen and nutrients, endothelial cells (ECs) oxidize only a minor fraction of glucose in mitochondria, a metabolic specialization that is poorly understood. Here we show that the glycolytic enzyme pyruvate kinase M2 (PKM2) limits glucose oxidation, and maintains the growth and epigenetic state of ECs. We find that loss of PKM2 alters mitochondrial substrate utilization and impairs EC proliferation and migration in vivo. Mechanistically, we show that the NF-κB transcription factor RELB is responsive to PKM2 loss, limiting EC growth through the regulation of P53. Furthermore, S-adenosylmethionine synthesis is impaired in the absence of PKM2, resulting in DNA hypomethylation, de-repression of endogenous retroviral elements (ERVs) and activation of antiviral innate immune signalling. This work reveals the metabolic and functional consequences of glucose oxidation in the endothelium, highlights the importance of PKM2 for endothelial growth and links metabolic dysfunction with autoimmune activation in ECs.

Enhancing Structure Prediction and Design of Soluble and Membrane Proteins with Explicit Solvent-Protein Interactions.

Solvent molecules interact intimately with proteins and can profoundly regulate their structure and function. However, accurately and efficiently modeling protein solvation effects at the molecular level has been challenging. Here, we present a method that improves the atomic-level modeling of soluble and membrane protein structures and binding by efficiently predicting de novo protein-solvent molecule interactions. The method predicted with unprecedented accuracy buried water molecule positions, solvated protein conformations, and challenging mutational effects on protein binding. When applied to homology modeling, solvent-bound membrane protein structures, pockets, and cavities were recapitulated with near-atomic precision even from distant homologs. Blindly refined atomic-level structures of evolutionary distant G protein-coupled receptors imply strikingly different functional roles of buried solvent between receptor classes. The method should prove useful for refining low-resolution protein structures, accurately modeling drug-binding sites in structurally uncharacterized receptors, and designing solvent-mediated protein catalysis, recognition, ligand binding, and membrane protein signaling.