Deglycosylated RBD produced in Pichia pastoris as a low-cost sera COVID-19 diagnosis tool and a vaccine candidate.

During the COVID-19 outbreak, numerous tools including protein-based vaccines have been developed. The methylotrophic yeast Pichia pastoris (synonymous to Komagataella phaffii) is an eukaryotic cost-effective and scalable system for recombinant protein production, with the advantages of an efficient secretion system and the protein folding assistance of the secretory pathway of eukaryotic cells. In a previous work, we compared the expression of SARS-CoV-2 Spike Receptor Binding Domain in P. pastoris with that in human cells. Although the size and glycosylation pattern was different between them, their protein structural and conformational features were indistinguishable. Nevertheless, since high mannose glycan extensions in proteins expressed by yeast may be the cause of a nonspecific immune recognition, we deglycosylated RBD in native conditions. This resulted in a highly pure, homogenous, properly folded and monomeric stable protein. This was confirmed by circular dichroism and tryptophan fluorescence spectra and by SEC-HPLC, which were similar to those of RBD proteins produced in yeast or human cells. Deglycosylated RBD was obtained at high yields in a single step, and it was efficient in distinguishing between SARS-CoV-2-negative and positive sera from patients. Moreover, when the deglycosylated variant was used as an immunogen, it elicited a humoral immune response ten times greater than the glycosylated form, producing antibodies with enhanced neutralizing power and eliciting a more robust cellular response. The proposed approach may be used to produce at a low cost, many antigens that require glycosylation to fold and express, but do not require glycans for recognition purposes.

Transcription factor NAC1 activates expression of peptidase-encoding AtCEPs in roots to limit root hair growth.

Plant genomes encode a unique group of papain-type Cysteine EndoPeptidases (CysEPs) containing a KDEL endoplasmic reticulum (ER) retention signal (KDEL-CysEPs or CEPs). CEPs process the cell-wall scaffolding EXTENSIN (EXT) proteins that regulate de novo cell-wall formation and cell expansion. Since CEPs cleave EXTs and EXT-related proteins, acting as cell-wall-weakening agents, they may play a role in cell elongation. The Arabidopsis (Arabidopsis thaliana) genome encodes 3 CEPs (AtCPE1-AtCEP3). Here, we report that the genes encoding these 3 Arabidopsis CEPs are highly expressed in root-hair (RH) cell files. Single mutants have no evident abnormal RH phenotype, but atcep1-3 atcep3-2 and atcep1-3 atcep2-2 double mutants have longer RHs than wild-type (Wt) plants, suggesting that expression of AtCEPs in root trichoblasts restrains polar elongation of the RH. We provide evidence that the transcription factor NAC1 (petunia NAM and Arabidopsis ATAF1, ATAF2, and CUC2) activates AtCEPs expression in roots to limit RH growth. Chromatin immunoprecipitation indicates that NAC1 binds to the promoter of AtCEP1, AtCEP2, and, to a lower extent, AtCEP3 and may directly regulate their expression. Inducible NAC1 overexpression increases AtCEP1 and AtCEP2 transcript levels in roots and leads to reduced RH growth while the loss of function nac1-2 mutation reduces AtCEP1-AtCEP3 gene expression and enhances RH growth. Likewise, expression of a dominant chimeric NAC1-SRDX repressor construct leads to increased RH length. Finally, we show that RH cell walls in the atcep1-3 atcep3-2 double mutant have reduced levels of EXT deposition, suggesting that the defects in RH elongation are linked to alterations in EXT processing and accumulation. Our results support the involvement of AtCEPs in controlling RH polar growth through EXT processing and insolubilization at the cell wall.

From in-silico screening to in-vitro evaluation: Enhancing the detection of Microcystins with engineered PP1 mutant variants.

Cyanotoxins produced during harmful algal blooms (CyanoHABs) have become a worldwide issue of concern. Microcystins (MC) are the most ubiquitous group of cyanotoxins and have known carcinogenic and hepatotoxic effects. The protein phosphatase inhibition assays (PPIAs), based on the inhibition of Protein Phosphatase 1/2A (PP1/PP2A) by MC, are one of the most cost-effective options for detecting MC. In this work, we aimed to design in-silico and evaluate in-vitro mutant variants of the PP1 protein, in order to enhance their capabilities as a MC biosensor. To this end, we performed an in-silico active site-saturated mutagenesis screening, followed by stability and docking affinity calculation with the MCLR cyanotoxin. Candidates with improved both affinity and stability were further tested in a fully flexible active-site docking. The best-scored mutations (19) were individually analysed regarding their locations and interactions. Four of them (p.D197F; p.Q249Y; p.S129W; p.D220Q) were selected for in-vitro expression and evaluation. Mutant p.D197F, exhibited a significant increment in inhibition by MCLR with respect to the WT, while showing a non-significant difference in stability nor activity. This successful PP1 inhibition enhancement suggests the potential of the p.D197F variant for practical MC detection applications.

Class III Peroxidases PRX01, PRX44, and PRX73 Control Root Hair Growth in Arabidopsis thaliana.

Root hair cells are important sensors of soil conditions. They grow towards and absorb water-soluble nutrients. This fast and oscillatory growth is mediated by continuous remodeling of the cell wall. Root hair cell walls contain polysaccharides and hydroxyproline-rich glycoproteins, including extensins (EXTs). Class-III peroxidases (PRXs) are secreted into the apoplastic space and are thought to trigger either cell wall loosening or polymerization of cell wall components, such as Tyr-mediated assembly of EXT networks (EXT-PRXs). The precise role of these EXT-PRXs is unknown. Using genetic, biochemical, and modeling approaches, we identified and characterized three root-hair-specific putative EXT-PRXs, PRX01, PRX44, and PRX73. prx01,44,73 triple mutation and PRX44 and PRX73 overexpression had opposite effects on root hair growth, peroxidase activity, and ROS production, with a clear impact on cell wall thickness. We use an EXT fluorescent reporter with contrasting levels of cell wall insolubilization in prx01,44,73 and PRX44-overexpressing background plants. In this study, we propose that PRX01, PRX44, and PRX73 control EXT-mediated cell wall properties during polar expansion of root hair cells.

An update on genetic variants of the NKX2-5.

NKX2-5 is a homeodomain transcription factor that plays a crucial role in heart development. It is the first gene where a single genetic variant (GV) was found to be associated with congenital heart diseases in humans. In this study, we carried out a comprehensive survey of NKX2-5 GVs to build a unified, curated, and updated compilation of all available GVs. We retrieved a total of 1,380 unique GVs. From these, 970 had information on their frequency in the general population and 143 have been linked to pathogenic phenotypes in humans. In vitro effect was ascertained for 38 GVs. The homeodomain had the biggest cluster of pathogenic variants in the protein: 49 GVs in 60 residues, 23 in its third α-helix, where 11 missense variants may affect protein-DNA interaction or the hydrophobic core. We also pinpointed the likely location of pathogenic GVs in four linear motifs. These analyses allowed us to assign a putative explanation for the effect of 90 GVs. This study pointed to reliable pathogenicity for GVs in helix 3 of the homeodomain and may broaden the scope of functional and structural studies that can be done to better understand the effect of GVs in NKX2-5 function.

Developing synthetic biology in Argentina: the Latin American TECNOx community as an alternative way for growth of the field.

Synthetic biology emerged in the USA and Europe twenty years ago and quickly developed innovative research and technology as a result of continued funding. Synthetic biology is also growing in many developing countries of Africa, Asia and Latin America, where it could have a large economic impact by helping its use of genetic biodiversity in order to boost existing industries. Starting in 2011, Argentine synthetic biology developed along an idiosyncratic path. In 2011-2012, the main focus was not exclusively research but also on community building through teaching and participation in iGEM, following the template of the early "MIT school" of synthetic biology. In 2013-2015, activities diversified and included society-centered projects, social science studies on synthetic biology and bioart. Standard research outputs such as articles and industrial applications helped consolidate several academic working groups. Since 2016, the lack of a critical mass of researchers and a funding crisis were partially compensated by establishing links with Latin American synthetic biologists and with other socially oriented open technology collectives. The TECNOx community is a central node in this growing research and technology network. The first four annual TECNOx meetings brought together synthetic biologists with other open science and engineering platforms and explored the relationship of Latin American technologies with entrepreneurship, open hardware, ethics and human rights. In sum, the socioeconomic context encouraged Latin American synthetic biology to develop in a meandering and diversifying manner. This revealed alternative ways for growth of the field that may be relevant to other developing countries.

Thermal Stability of Globins: Implications of Flexibility and Heme Coordination Studied by Molecular Dynamics Simulations.

Proteins are sensitive to temperature, and abrupt changes in the normal temperature conditions can have a profound impact on both structure and function, leading to protein unfolding. However, the adaptation of certain organisms to extreme conditions raises questions about the structural features that permit the structure and function of proteins to be preserved under these adverse conditions. To gain insight into the molecular basis of protein thermostability in the globin family, we have examined three representative examples: human neuroglobin, horse heart myoglobin, and Drosophila hemoglobin, which differ in their melting temperatures and coordination states of the heme iron in the absence of external ligands. In order to elucidate the possible mechanisms that govern the thermostability of these proteins, microsecond-scale classical molecular dynamics simulations were performed at different temperatures. Structural fluctuations and essential dynamics were analyzed, indicating that the flexibility of the CD region, which includes the two short C and D helixes and the connecting CD loop, is directly related to the thermostability. We observed that a larger inherent flexibility of the protein produces higher thermostability, probably concentrating the thermal fluctuations observed at high temperature in flexible regions, preventing unfolding. Globally, the results of this work improve our understanding of thermostability in the globin family.

CYP21A2 mutation update: Comprehensive analysis of databases and published genetic variants.

Congenital adrenal hyperplasia (CAH) is a group of autosomal recessive disorders of adrenal steroidogenesis. Disorders in steroid 21-hydroxylation account for over 95% of patients with CAH. Clinically, the 21-hydroxylase deficiency has been classified in a broad spectrum of clinical forms, ranging from severe or classical, to mild late onset or non-classical. Known allelic variants in the disease causing CYP21A2 gene are spread among different sources. Until recently, most variants reported have been identified in the clinical setting, which presumably bias described variants to pathogenic ones, as those found in the CYPAlleles database. Nevertheless, a large number of variants are being described in massive genome projects, many of which are found in dbSNP, but lack functional implications and/or their phenotypic effect. In this work, we gathered a total of 1,340 GVs in the CYP21A2 gene, from which 899 variants were unique and 230 have an effect on human health, and compiled all this information in an integrated database. We also connected CYP21A2 sequence information to phenotypic effects for all available mutations, including double mutants in cis. Data compiled in the present work could help physicians in the genetic counseling of families affected with 21-hydroxylase deficiency.

Engineered chimeras reveal the structural basis of hexacoordination in globins: a case study of neuroglobin and myoglobin.

BACKGROUND: Myoglobin (Mb) and neuroglobin (Ngb) are representative members of pentacoordinated and bis-histidyl, hexacoordinated globins. In spite of their low sequence identity, they show surprisingly similar three-dimensional folds. The ability of Ngb to form a hexacoordinated bis-histidyl complex with the distal HisE7 has a strong impact on ligand affinity. The factors governing such different behaviors have not been completely understood yet, even though they are extremely relevant to establish structure-function relationships within the globin superfamily. METHODS: In this work we generated chimeric proteins by swapping a previously identified regulatory segment - the CD region - and evaluated comparatively the structural and functional properties of the resulting proteins by molecular dynamics simulations, and spectroscopic and kinetic investigations. RESULTS: Our results show that chimeric proteins display heme coordination properties displaced towards those expected for the corresponding CD region. In particular, in the absence of exogenous ligands, chimeric Mb is found as a partially hexacoordinated bis-histidyl species, whereas chimeric Ngb shows a lower equilibrium constant for forming a hexacoordinated bis-histidyl species. CONCLUSIONS: While these results confirm the regulatory role of the CD region for bis-histidyl hexacoordination, they also suggest that additional sources contribute to fine tune the equilibrium. General significance Globins constitute a ubiquitous group of heme proteins widely found in all kingdoms of life. These findings raise challenging questions regarding the structure-function relationships in these proteins, as bis-histidyl hexacoordination emerges as a novel regulatory mechanism of the physiological function of globins.

Circadian period integrates network information through activation of the BMP signaling pathway.

Living organisms use biological clocks to maintain their internal temporal order and anticipate daily environmental changes. In Drosophila, circadian regulation of locomotor behavior is controlled by ∼150 neurons; among them, neurons expressing the PIGMENT DISPERSING FACTOR (PDF) set the period of locomotor behavior under free-running conditions. To date, it remains unclear how individual circadian clusters integrate their activity to assemble a distinctive behavioral output. Here we show that the BONE MORPHOGENETIC PROTEIN (BMP) signaling pathway plays a crucial role in setting the circadian period in PDF neurons in the adult brain. Acute deregulation of BMP signaling causes period lengthening through regulation of dClock transcription, providing evidence for a novel function of this pathway in the adult brain. We propose that coherence in the circadian network arises from integration in PDF neurons of both the pace of the cell-autonomous molecular clock and information derived from circadian-relevant neurons through release of BMP ligands.

Molecular basis of thermal stability in truncated (2/2) hemoglobins.

BACKGROUND: Understanding the molecular mechanism through which proteins are functional at extreme high and low temperatures is one of the key issues in structural biology. To investigate this phenomenon, we have focused on two instructive truncated hemoglobins from Thermobifida fusca (Tf-trHbO) and Mycobacterium tuberculosis (Mt-trHbO); although the two proteins are structurally nearly identical, only the former is stable at high temperatures. METHODS: We used molecular dynamics simulations at different temperatures as well as thermal melting profile measurements of both wild type proteins and two mutants designed to interchange the amino acid residue, either Pro or Gly, at E3 position. RESULTS: The results show that the presence of a Pro at the E3 position is able to increase (by 8°) or decrease (by 4°) the melting temperature of Mt-trHbO and Tf-trHbO, respectively. We observed that the ProE3 alters the structure of the CD loop, making it more flexible. CONCLUSIONS: This gain in flexibility allows the protein to concentrate its fluctuations in this single loop and avoid unfolding. The alternate conformations of the CD loop also favor the formation of more salt-bridge interactions, together augmenting the protein's thermostability. GENERAL SIGNIFICANCE: These results indicate a clear structural and dynamical role of a key residue for thermal stability in truncated hemoglobins.

Protein dynamics and ligand migration interplay as studied by computer simulation.

Since proteins are dynamic systems in living organisms, the employment of methodologies contemplating this crucial characteristic results fundamental to allow revealing several aspects of their function. In this work, we present results obtained using classical mechanical atomistic simulation tools applied to understand the connection between protein dynamics and ligand migration. Firstly, we will present a review of the different sampling schemes used in the last years to obtain both ligand migration pathways and the thermodynamic information associated with the process. Secondly, we will focus on representative examples in which the schemes previously presented are employed, concerning the following: i) ligand migration, tunnels, and cavities in myoglobin and neuroglobin; ii) ligand migration in truncated hemoglobin members; iii) NO escape and conformational changes in nitrophorins; iv) ligand selectivity in catalase and hydrogenase; and v) larger ligand migration: the P450 and haloalkane dehalogenase cases. This article is part of a Special Issue entitled: Protein Dynamics: Experimental and Computational Approaches.

Using protein design algorithms to understand the molecular basis of disease caused by protein-DNA interactions: the Pax6 example.

Quite often a single or a combination of protein mutations is linked to specific diseases. However, distinguishing from sequence information which mutations have real effects in the protein's function is not trivial. Protein design tools are commonly used to explain mutations that affect protein stability, or protein-protein interaction, but not for mutations that could affect protein-DNA binding. Here, we used the protein design algorithm FoldX to model all known missense mutations in the paired box domain of Pax6, a highly conserved transcription factor involved in eye development and in several diseases such as aniridia. The validity of FoldX to deal with protein-DNA interactions was demonstrated by showing that high levels of accuracy can be achieved for mutations affecting these interactions. Also we showed that protein-design algorithms can accurately reproduce experimental DNA-binding logos. We conclude that 88% of the Pax6 mutations can be linked to changes in intrinsic stability (77%) and/or to its capabilities to bind DNA (30%). Our study emphasizes the importance of structure-based analysis to understand the molecular basis of diseases and shows that protein-DNA interactions can be analyzed to the same level of accuracy as protein stability, or protein-protein interactions.

Transcription factor specificity limits the number of DNA-binding motifs.

We study the limits imposed by transcription factor specificity on the maximum number of binding motifs that can coexist in a gene regulatory network, using the SwissRegulon Fantom5 collection of 684 human transcription factor binding sites as a model. We describe transcription factor specificity using regular expressions and find that most human transcription factor binding site motifs are separated in sequence space by one to three motif-discriminating positions. We apply theorems based on the pigeonhole principle to calculate the maximum number of transcription factors that can coexist given this degree of specificity, which is in the order of ten thousand and would fully utilize the space of DNA subsequences. Taking into account an expanded DNA alphabet with modified bases can further raise this limit by several orders of magnitude, at a lower level of sequence space usage. Our results may guide the design of transcription factors at both the molecular and system scale.

Apoplastic class III peroxidases PRX62 and PRX69 promote Arabidopsis root hair growth at low temperature.

Root Hairs (RHs) growth is influenced by endogenous and by external environmental signals that coordinately regulate its final cell size. We have recently determined that RH growth was unexpectedly boosted when Arabidopsis thaliana seedlings are cultivated at low temperatures. It was proposed that RH growth plasticity in response to low temperature was linked to a reduced nutrient availability in the media. Here, we explore the molecular basis of this RH growth response by using a Genome Wide Association Study (GWAS) approach using Arabidopsis thaliana natural accessions. We identify the poorly characterized PEROXIDASE 62 (PRX62) and a related protein PRX69 as key proteins under moderate low temperature stress. Strikingly, a cell wall protein extensin (EXT) reporter reveals the effect of peroxidase activity on EXT cell wall association at 10 °C in the RH apical zone. Collectively, our results indicate that PRX62, and to a lesser extent PRX69, are key apoplastic PRXs that modulate ROS-homeostasis and cell wall EXT-insolubilization linked to RH elongation at low temperature.

Structure and dynamics of Antarctic fish neuroglobin assessed by computer simulations.

Neuroglobin (Ngb) is a heme protein, highly conserved along evolution, predominantly found in the nervous system. It is upregulated by hypoxia and ischemia and may have a neuroprotective role under hypoxic stress. Although many other roles have been proposed, the physiological function is still unclear. Antarctic icefishes lack hemoglobin and some species also lack myoglobin, but all have Ngb and thus may help the elucidation of Ngb function. We present the first theoretically derived structure of fish Ngb and describe its behavior using molecular dynamics simulations. Specifically, we sequenced and analyzed Ngbs from a colorless-blooded Antarctic icefish species Chaenocephalus aceratus and a related red-blooded species (Dissostichus mawsoni). Both fish Ngbs are 6-coordinated but have some peculiarities that differentiate them from mammalian counterparts: they have extensions in the N and C termini that can interact with the EF loop, and a gap in the alignment that changes the CD-region structure/dynamics that has been found to play a key role in human neuroglobin. Our results suggest that a single mutation between both fish Ngbs is responsible for significant difference in the behavior of the proteins. The functional role of these characteristics is discussed.

Crystal structure of I-DmoI in complex with its target DNA provides new insights into meganuclease engineering.

Homing endonucleases, also known as meganucleases, are sequence-specific enzymes with large DNA recognition sites. These enzymes can be used to induce efficient homologous gene targeting in cells and plants, opening perspectives for genome engineering with applications in a wide series of fields, ranging from biotechnology to gene therapy. Here, we report the crystal structures at 2.0 and 2.1 A resolution of the I-DmoI meganuclease in complex with its substrate DNA before and after cleavage, providing snapshots of the catalytic process. Our study suggests that I-DmoI requires only 2 cations instead of 3 for DNA cleavage. The structure sheds light onto the basis of DNA binding, indicating key residues responsible for nonpalindromic target DNA recognition. In silico and in vivo analysis of the I-DmoI DNA cleavage specificity suggests that despite the relatively few protein-base contacts, I-DmoI is highly specific when compared with other meganucleases. Our data open the door toward the generation of custom endonucleases for targeted genome engineering using the monomeric I-DmoI scaffold.

Asymmetric bifunctional protein nanoparticles through redesign of self-assembly.

Engineering oligomeric protein self-assembly is an attractive approach to fabricate nanostructures with well-defined geometries, stoichiometry and functions. The homodecamer Brucella Lumazine Synthase (BLS) is a highly stable and immunogenic protein nanoparticle (PNP). Here, we engineered the BLS protein scaffold to display two functions in spatially opposite regions of its structure yielding a Janus-like nanoparticle. An in silico analysis of the BLS head-to-head dimer of homopentamers shows major inter-pentameric interactions located in the equatorial interface. Based on this analysis, two BLS protomer variants were designed to interrupt pentamer self-dimerization and promote heteropentameric dimers. This strategy enabled us to generate a decameric particle with two distinct sides formed by two independent pentamers. The versatility of this new self-assembly nanofabrication strategy is illustrated with two example applications. First, a bifunctional BLS bearing Alexa Fluor 488 fluorophores on one side and sialic acid binding domains on the other side was used for labelling murine and human cells and analyzed by flow cytometry and confocal microscopy. Second, multichromophoric FRET nanoparticles were fabricated and characterized at the single molecule level, showing discrete energy transfer events. The engineered BLS variants constitute a general platform for displaying two functions in a controlled manner within the same PNP with potential applications in various areas such as biomedicine, biotechnology and nanotechnology.

Exploring the molecular basis of heme coordination in human neuroglobin.

Neuroglobin (Ngb), a recently discovered ancient heme protein, presents the typical globin fold and is around 20% identical to myoglobin (Mb). In contrast with Mb, however, its heme is hexacoordinated (6c). It is expressed in the nervous system and has been the subject of numerous investigations in the last years, but its function is still unclear. The proposed roles include oxygen transport, reactive oxygen species (ROS) detoxification, hypoxia protection, and redox state sensing. All proposed functions require distal histidine dissociation from the heme to yield a reactive iron. With the aim of understanding the 6c to 5c transition, we have performed molecular dynamics simulations for ferrous Ngb in the 6c, 5c, and oxy states. We also computed free energy profiles associated with the transition employing an advanced sampling technique. Finally, we studied the effect of the redox state of CysCD7 and CysD5, which are known to form a disulfide bridge. Our results show that protein oxidation promotes a stabilization of the pentacoordinated species, thus favoring the protein to adopt the more reactive state and supporting the existence of a molecular mechanism whereby O2 would be released under hypoxic conditions, thereby suggesting an O(2) storage function for Ngb. Taken together, our results provide structural information not available experimentally which may shed light on the protein proposed functions, particularly as a redox sensor.