Solvent constraints for biopolymer folding and evolution in extraterrestrial environments.

We propose that spontaneous folding and molecular evolution of biopolymers are two universal aspects that must concur for life to happen. These aspects are fundamentally related to the chemical composition of biopolymers and crucially depend on the solvent in which they are embedded. We show that molecular information theory and energy landscape theory allow us to explore the limits that solvents impose on biopolymer existence. We consider 54 solvents, including water, alcohols, hydrocarbons, halogenated solvents, aromatic solvents, and low molecular weight substances made up of elements abundant in the universe, which may potentially take part in alternative biochemistries. We find that along with water, there are many solvents for which the liquid regime is compatible with biopolymer folding and evolution. We present a ranking of the solvents in terms of biopolymer compatibility. Many of these solvents have been found in molecular clouds or may be expected to occur in extrasolar planets.

Associating protein sequence positions with the modulation of quantitative phenotypes.

Although protein sequences encode the information for folding and function, understanding their link is not an easy task. Unluckily, the prediction of how specific amino acids contribute to these features is still considerably impaired. Here, we developed a simple algorithm that finds positions in a protein sequence with potential to modulate the studied quantitative phenotypes. From a few hundred protein sequences, we perform multiple sequence alignments, obtain the per-position pairwise differences for both the sequence and the observed phenotypes, and calculate the correlation between these last two quantities. We tested our methodology with four cases: archaeal Adenylate Kinases and the organisms optimal growth temperatures, microbial rhodopsins and their maximal absorption wavelengths, mammalian myoglobins and their muscular concentration, and inhibition of HIV protease clinical isolates by two different molecules. We found from 3 to 10 positions tightly associated with those phenotypes, depending on the studied case. We showed that these correlations appear using individual positions but an improvement is achieved when the most correlated positions are jointly analyzed. Noteworthy, we performed phenotype predictions using a simple linear model that links per-position divergences and differences in the observed phenotypes. Predictions are comparable to the state-of-art methodologies which, in most of the cases, are far more complex. All of the calculations are obtained at a very low information cost since the only input needed is a multiple sequence alignment of protein sequences with their associated quantitative phenotypes. The diversity of the explored systems makes our work a valuable tool to find sequence determinants of biological activity modulation and to predict various functional features for uncharacterized members of a protein family.

The evolving copiotrophic/oligotrophic dichotomy: From Winogradsky to physiology and genomics.

Nearly 100 years ago, Winogradsky published a classic communication in which he described two groups of microbes, zymogenic and autochthonous. When organic matter penetrates the soil, zymogenic microbes quickly multiply and degrade it, then giving way to the slow combustion of autochthonous microbes. Although the text was originally written in French, it is often cited by English-speaking authors. We undertook a complete translation of the 1924 publication, which we provide as Supporting information. Here, we introduce the translation and describe how the zymogenic/autochthonous dichotomy shaped research questions in the study of microbial diversity and physiology. We also identify in the literature three additional and closely related dichotomies, which we propose to call exclusive copiotrophs/oligotrophs, coexisting copiotrophs/oligotrophs and fast-growing/slow-growing microbes. While Winogradsky focussed on a successional view of microbial populations over time, the current discussion is focussed on the differences in the specific growth rate of microbes as a function of the concentration of a given limiting substrate. In the future, it will be relevant to keep in mind both nutrient-focussed and time-focussed microbial dichotomies and to design experiments with both isolated laboratory cultures and multi-species communities in the spirit of Winogradsky's direct method.

Ecology theory disentangles microbial dichotomies.

Microbes are often discussed in terms of dichotomies such as copiotrophic/oligotrophic and fast/slow-growing microbes, defined using the characterisation of microbial growth in isolated cultures. The dichotomies are usually qualitative and/or study-specific, sometimes precluding clear-cut results interpretation. We can unravel microbial dichotomies as life history strategies by combining ecology theory with Monod curves, a laboratory mathematical tool of bacterial physiology that relates the specific growth rate of a microbe with the concentration of a limiting nutrient. Fitting of Monod curves provides quantities that directly correspond to key parameters in ecological theories addressing species coexistence and diversity, such as r/K selection theory, resource competition and community structure theory and the CSR triangle of life strategies. The resulting model allows us to reconcile the copiotrophic/oligotrophic and fast/slow-growing dichotomies as different subsamples of a life history strategy triangle that also includes r/K strategists. We also used the number of known carbon sources together with community structure theory to partially explain the diversity of heterotrophic microbes observed in metagenomics experiments. In sum, we propose a theoretical framework for the study of natural microbial communities that unifies several existing proposals. Its application would require the integration of metagenomics, metametabolomics, Monod curves and carbon source data.

r/K selection of GC content in prokaryotes.

The guanine/cytosine (GC) content of prokaryotic genomes is species-specific, taking values from 16% to 77%. This diversity of selection for GC content remains contentious. We analyse the correlations between GC content and a range of phenotypic and genotypic data in thousands of prokaryotes. GC content integrates well with these traits into r/K selection theory when phenotypic plasticity is considered. High GC-content prokaryotes are r-strategists with cheaper descendants thanks to a lower average amino acid metabolic cost, colonize unstable environments thanks to flagella and a bacillus form and are generalists in terms of resource opportunism and their defence mechanisms. Low GC content prokaryotes are K-strategists specialized for stable environments that maintain homeostasis via a high-cost outer cell membrane and endospore formation as a response to nutrient deprivation, and attain a higher nutrient-to-biomass yield. The lower proteome cost of high GC content prokaryotes is driven by the association between GC-rich codons and cheaper amino acids in the genetic code, while the correlation between GC content and genome size may be partly due to functional diversity driven by r/K selection. In all, molecular diversity in the GC content of prokaryotes may be a consequence of ecological r/K selection.

NGOME-Lite: Proteome-wide prediction of spontaneous protein deamidation highlights differences between taxa.

Asparagines in proteins deamidate spontaneously, which changes the chemical structure of a protein and often affects its function. Current prediction algorithms for asparagine deamidation require a structure as an input or are too slow to be applied at a proteomic scale. We present NGOME-Lite, a new version of our sequence-based predictor for spontaneous asparagine deamidation that is faster by over two orders of magnitude at a similar degree of accuracy. The algorithm takes into account intrinsic sequence propensities and slowing down of deamidation by local structure. NGOME-Lite can run in a proteomic analysis mode that provides the half-time of the intact form of each protein, predicted by taking into account sequence propensities and structural protection or sequence propensities only, and a structure protection factor. The detailed analysis mode also provides graphical output for all Asn residues in the query sequence. We applied NGOME-Lite to over 257,000 sequences in 38 proteomes and found that different taxa differ in their predicted deamidation dynamics. Spontaneous protein deamidation is faster in Eukarya than in Bacteria because of a higher degree of structural protection in the latter. Predicted protein deamidation half-lifes correlate with protein turnover in human, mouse, rat, C. elegans and budding yeast but not in two plants and two bacteria. NGOME-Lite is implemented in a docker container available at https://ngome.proteinphysiologylab.org.

Deamidation drives molecular aging of the SARS-CoV-2 spike protein receptor-binding motif.

The spike protein is the main protein component of the SARS-CoV-2 virion surface. The spike receptor-binding motif mediates recognition of the human angiotensin-converting enzyme 2 receptor, a critical step in infection, and is the preferential target for spike-neutralizing antibodies. Posttranslational modifications of the spike receptor-binding motif have been shown to modulate viral infectivity and host immune response, but these modifications are still being explored. Here we studied asparagine deamidation of the spike protein, a spontaneous event that leads to the appearance of aspartic and isoaspartic residues, which affect both the protein backbone and its charge. We used computational prediction and biochemical experiments to identify five deamidation hotspots in the SARS-CoV-2 spike protein. Asparagine residues 481 and 501 in the receptor-binding motif deamidate with a half-life of 16.5 and 123 days at 37 °C, respectively. Deamidation is significantly slowed at 4 °C, indicating a strong dependence of spike protein molecular aging on environmental conditions. Deamidation of the spike receptor-binding motif decreases the equilibrium constant for binding to the human angiotensin-converting enzyme 2 receptor more than 3.5-fold, yet its high conservation pattern suggests some positive effect on viral fitness. We propose a model for deamidation of the full SARS-CoV-2 virion illustrating how deamidation of the spike receptor-binding motif could lead to the accumulation on the virion surface of a nonnegligible chemically diverse spike population in a timescale of days. Our findings provide a potential mechanism for molecular aging of the spike protein with significant consequences for understanding virus infectivity and vaccine development.

Topology Dictates Evolution of Regulatory Cysteines in a Family of Viral Oncoproteins.

Redox regulation in biology is largely operated by cysteine chemistry in response to a variety of cell environmental and intracellular stimuli. The high chemical reactivity of cysteines determines their conservation in functional roles, but their presence can also result in harmful oxidation limiting their general use by proteins. Papillomaviruses constitute a unique system for studying protein sequence evolution since there are hundreds of anciently evolved stable genomes. E7, the viral transforming factor, is a dimeric, cysteine-rich oncoprotein that shows both conserved structural and variable regulatory cysteines constituting an excellent model for uncovering the mechanism that drives the acquisition of redox-sensitive groups. By analyzing over 300 E7 sequences, we found that although noncanonical cysteines show no obvious sequence conservation pattern, they are nonrandomly distributed based on topological constrains. Regulatory residues are strictly excluded from six positions stabilizing the hydrophobic core while they are enriched in key positions located at the dimerization interface or around the Zn+2 ion. Oxidation of regulatory cysteines is linked to dimer dissociation, acting as a reversible redox-sensing mechanism that triggers a conformational switch. Based on comparative sequence analysis, molecular dynamics simulations and biophysical analysis, we propose a model in which the occurrence of cysteine-rich positions is dictated by topological constrains, providing an explanation to why a degenerate pattern of cysteines can be achieved in a family of homologs. Thus, topological principles should enable the possibility to identify hidden regulatory cysteines that are not accurately detected using sequence based methodology.

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.

Hidden Structural Codes in Protein Intrinsic Disorder.

Intrinsic disorder is a major structural category in biology, accounting for more than 30% of coding regions across the domains of life, yet consists of conformational ensembles in equilibrium, a major challenge in protein chemistry. Anciently evolved papillomavirus genomes constitute an unparalleled case for sequence to structure-function correlation in cases in which there are no folded structures. E7, the major transforming oncoprotein of human papillomaviruses, is a paradigmatic example among the intrinsically disordered proteins. Analysis of a large number of sequences of the same viral protein allowed for the identification of a handful of residues with absolute conservation, scattered along the sequence of its N-terminal intrinsically disordered domain, which intriguingly are mostly leucine residues. Mutation of these led to a pronounced increase in both α-helix and ß-sheet structural content, reflected by drastic effects on equilibrium propensities and oligomerization kinetics, and uncovers the existence of local structural elements that oppose canonical folding. These folding relays suggest the existence of yet undefined hidden structural codes behind intrinsic disorder in this model protein. Thus, evolution pinpoints conformational hot spots that could have not been identified by direct experimental methods for analyzing or perturbing the equilibrium of an intrinsically disordered protein ensemble.

The E. coli thioredoxin folding mechanism: the key role of the C-terminal helix.

In this work, the unfolding mechanism of oxidized Escherichia coli thioredoxin (EcTRX) was investigated experimentally and computationally. We characterized seven point mutants distributed along the C-terminal α-helix (CTH) and the preceding loop. The mutations destabilized the protein against global unfolding while leaving the native structure unchanged. Global analysis of the unfolding kinetics of all variants revealed a linear unfolding route with a high-energy on-pathway intermediate state flanked by two transition state ensembles TSE1 and TSE2. The experiments show that CTH is mainly unfolded in TSE1 and the intermediate and becomes structured in TSE2. Structure-based molecular dynamics are in agreement with these experiments and provide protein-wide structural information on transient states. In our model, EcTRX folding starts with structure formation in the ß-sheet, while the protein helices coalesce later. As a whole, our results indicate that the CTH is a critical module in the folding process, restraining a heterogeneous intermediate ensemble into a biologically active native state and providing the native protein with thermodynamic and kinetic stability.

The eukaryotic linear motif resource ELM: 10 years and counting.

The eukaryotic linear motif (ELM http://elm.eu.org) resource is a hub for collecting, classifying and curating information about short linear motifs (SLiMs). For >10 years, this resource has provided the scientific community with a freely accessible guide to the biology and function of linear motifs. The current version of ELM contains ∼200 different motif classes with over 2400 experimentally validated instances manually curated from >2000 scientific publications. Furthermore, detailed information about motif-mediated interactions has been annotated and made available in standard exchange formats. Where appropriate, links are provided to resources such as switches.elm.eu.org and KEGG pathways.

Conformational Heterogeneity Determined by Folding and Oligomer Assembly Routes of the Interferon Response Inhibitor NS1 Protein, Unique to Human Respiratory Syncytial Virus.

The nonstructural NS1 protein is an essential virulence factor of the human respiratory syncytial virus, with a predominant role in the inhibition of the host antiviral innate immune response. This inhibition is mediated by multiple protein-protein interactions and involves the formation of large oligomeric complexes. There is neither a structure nor sequence or functional homologues of this protein, which points to a distinctive mechanism for blocking the interferon response among viruses. The NS1 native monomer follows a simple unfolding kinetics via a nativelike transition state ensemble, with a half-life of 45 min, in agreement with a highly stable core structure at equilibrium. Refolding is a complex process that involves several slowly interconverting species compatible with proline isomerization. However, an ultrafast folding event with a half-life of 0.2 ms is indicative of a highly folding compatible species within the unfolded state ensemble. On the other hand, the oligomeric assembly route from the native monomer, which does not involve unfolding, shows a monodisperse and irreversible end-point species triggered by a mild temperature change, with half-lives of 160 and 26 min at 37 and 47 °C, respectively, and at a low protein concentration (10 µM). A large secondary structure change into ß-sheet structure and the formation of a dimeric nucleus precede polymerization by the sequential addition of monomers at the surprisingly low rate of one monomer every 34 s. The polymerization phase is followed by the binding to thioflavin-T indicative of amyloid-like, albeit soluble, repetitive ß-sheet quaternary structure. The overall process is reversible only up until ~8 min, a time window in which most of the secondary structure change takes place. NS1's multiple binding activities must be accommodated in a few binding interfaces at most, something to be considered remarkable given its small size (15 kDa). Thus, conformational heterogeneity, and in particular oligomer formation, may provide a means of expand its binding repertoire. These equilibria will be determined by variables such as macromolecular crowding, protein-protein interactions, expression levels, turnover, or specific subcellular localization. The irreversible and quasi-spontaneous nature of the oligomer assembly, together with the fact that NS1 is the most abundant viral protein in infected cells, makes its accumulation highly conceivable under conditions compatible with the cellular milieu. The implications of NS1 oligomers in the viral life cycle and the inhibition of host innate immune response remain to be determined.

Amino Acid metabolism conflicts with protein diversity.

The 20 protein-coding amino acids are found in proteomes with different relative abundances. The most abundant amino acid, leucine, is nearly an order of magnitude more prevalent than the least abundant amino acid, cysteine. Amino acid metabolic costs differ similarly, constraining their incorporation into proteins. On the other hand, a diverse set of protein sequences is necessary to build functional proteomes. Here, we present a simple model for a cost-diversity trade-off postulating that natural proteomes minimize amino acid metabolic flux while maximizing sequence entropy. The model explains the relative abundances of amino acids across a diverse set of proteomes. We found that the data are remarkably well explained when the cost function accounts for amino acid chemical decay. More than 100 organisms reach comparable solutions to the trade-off by different combinations of proteome cost and sequence diversity. Quantifying the interplay between proteome size and entropy shows that proteomes can get optimally large and diverse.

Cysteine-rich positions outside the structural zinc motif of human papillomavirus E7 provide conformational modulation and suggest functional redox roles.

The E7 protein from high-risk human papillomavirus is essential for cell transformation in cervical, oropharyngeal, and other HPV-related cancers, mainly through the inactivation of the retinoblastoma (Rb) tumor suppressor. Its high cysteine content (~7%) and the observation that HPV-transformed cells are under oxidative stress prompted us to investigate the redox properties of the HPV16 E7 protein under biologically compatible oxidative conditions. The seven cysteines in HPV16 E7 remain reduced in conditions resembling the basal reduced state of a cell. However, under oxidative stress, a stable disulfide bridge forms between cysteines 59 and 68. Residue 59 has a protective effect on the other cysteines, and its mutation leads to an overall increase in the oxidation propensity of E7, including cysteine 24 central to the Rb binding motif. Gluthationylation of Cys 24 abolishes Rb binding, which is reversibly recovered upon reduction. Cysteines 59 and 68 are located 18.6 Å apart, and the formation of the disulfide bridge leads to a large structural rearrangement while retaining strong Zn association. These conformational and covalent changes are fully reversible upon restoration of the reductive environment. In addition, this is the first evidence of an interaction between the N-terminal intrinsically disordered and the C-terminal globular domains, known to be highly and separately conserved among human papillomaviruses. The significant conservation of such noncanonical cysteines in HPV E7 proteins leads us to propose a functional redox activity. Such an activity adds to the previously discovered chaperone activity of E7 and supports the picture of a moonlighting pathological role of this paradigmatic viral oncoprotein.

Folding of a cyclin box: linking multitarget binding to marginal stability, oligomerization, and aggregation of the retinoblastoma tumor suppressor AB pocket domain.

The retinoblastoma tumor suppressor (Rb) controls the proliferation, differentiation, and survival of cells in most eukaryotes with a role in the fate of stem cells. Its inactivation by mutation or oncogenic viruses is required for cellular transformation and eventually carcinogenesis. The high conservation of the Rb cyclin fold prompted us to investigate the link between conformational stability and ligand binding properties of the RbAB pocket domain. RbAB unfolding presents a three-state transition involving cooperative secondary and tertiary structure changes and a partially folded intermediate that can oligomerize. The first transition corresponds to unfolding of the metastable B subdomain containing the binding site for the LXCXE motif present in cellular and viral targets, and the second transition corresponds to the stable A subdomain. The low thermodynamic stability of RbAB translates into a propensity to rapidly oligomerize and aggregate at 37 °C (T50 = 28 min) that is suppressed by human papillomavirus E7 and E2F peptide ligands, suggesting that Rb is likely stabilized in vivo through binding to target proteins. We propose that marginal stability and associated oligomerization may be conserved for function as a "hub" protein, allowing the formation of multiprotein complexes, which could constitute a robust mechanism to retain its cell cycle regulatory role throughout evolution. Decreased stability and oligomerization are shared with the p53 tumor suppressor, suggesting a link between folding and function in these two essential cell regulators that are inactivated in most cancers and operate within multitarget signaling pathways.

Minute time scale prolyl isomerization governs antibody recognition of an intrinsically disordered immunodominant epitope.

Conformational rearrangements in antibody·antigen recognition are essential events where kinetic discrimination of isomers expands the universe of combinations. We investigated the interaction mechanism of a monoclonal antibody, M1, raised against E7 from human papillomavirus, a prototypic viral oncoprotein and a model intrinsically disordered protein. The mapped 12-amino acid immunodominant epitope lies within a "hinge" region between the N-terminal intrinsically disordered and the C-terminal globular domains. Kinetic experiments show that despite being within an intrinsically disordered region, the hinge E7 epitope has at least two populations separated by a high energy barrier. Nuclear magnetic resonance traced the origin of this barrier to a very slow (t(1/2)∼4 min) trans-cis prolyl isomerization event involving changes in secondary structure. The less populated (10%) cis isomer is the binding-competent species, thus requiring the 90% of molecules in the trans configuration to isomerize before binding. The association rate for the cis isomer approaches 6 × 10(7) M(-1) s(-1), a ceiling for antigen-antibody interactions. Mutagenesis experiments showed that Pro-41 in E7Ep was required for both binding and isomerization. After a slow postbinding unimolecular rearrangement, a consolidated complex with K(D) = 1.2 × 10(-7) M is reached. Our results suggest that presentation of this viral epitope by the antigen-presenting cells would have to be "locked" in the cis conformation, in opposition to the most populated trans isomer, in order to select the specific antibody clone that goes through affinity and kinetic maturation.

TDP-43 proteinopathy in ALS is triggered by loss of ASRGL1 and associated with HML-2 expression.

TAR DNA-binding protein 43 (TDP-43) proteinopathy in brain cells is the hallmark of amyotrophic lateral sclerosis (ALS) but its cause remains elusive. Asparaginase-like-1 protein (ASRGL1) cleaves isoaspartates, which alter protein folding and susceptibility to proteolysis. ASRGL1 gene harbors a copy of the human endogenous retrovirus HML-2, whose overexpression contributes to ALS pathogenesis. Here we show that ASRGL1 expression was diminished in ALS brain samples by RNA sequencing, immunohistochemistry, and western blotting. TDP-43 and ASRGL1 colocalized in neurons but, in the absence of ASRGL1, TDP-43 aggregated in the cytoplasm. TDP-43 was found to be prone to isoaspartate formation and a substrate for ASRGL1. ASRGL1 silencing triggered accumulation of misfolded, fragmented, phosphorylated and mislocalized TDP-43 in cultured neurons and motor cortex of female mice. Overexpression of ASRGL1 restored neuronal viability. Overexpression of HML-2 led to ASRGL1 silencing. Loss of ASRGL1 leading to TDP-43 aggregation may be a critical mechanism in ALS pathophysiology.

Experimental snapshots of a protein-DNA binding landscape.

Protein recognition of DNA sites is a primary event for gene function. Its ultimate mechanistic understanding requires an integrated structural, dynamic, kinetic, and thermodynamic dissection that is currently limited considering the hundreds of structures of protein-DNA complexes available. We describe a protein-DNA-binding pathway in which an initial, diffuse, transition state ensemble with some nonnative contacts is followed by formation of extensive nonnative interactions that drive the system into a kinetic trap. Finally, nonnative contacts are slowly rearranged into native-like interactions with the DNA backbone. Dissimilar protein-DNA interfaces that populate along the DNA-binding route are explained by a temporary degeneracy of protein-DNA interactions, centered on "dual-role" residues. The nonnative species slow down the reaction allowing for extended functionality.

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.