RESUMO
Several mammalian genes have originated from the domestication of retrotransposons, selfish mobile elements related to retroviruses. Some of the proteins encoded by these genes have maintained virus-like features; including self-processing, capsid structure formation, and the generation of different isoforms through -1 programmed ribosomal frameshifting. Using quantitative approaches in molecular evolution and biophysical analyses, we studied 28 retrotransposon-derived genes, with a focus on the evolution of virus-like features. By analyzing the rate of synonymous substitutions, we show that the -1 programmed ribosomal frameshifting mechanism in three of these genes (PEG10, PNMA3, and PNMA5) is conserved across mammals and originates alternative proteins. These genes were targets of positive selection in primates, and one of the positively selected sites affects a B-cell epitope on the spike domain of the PNMA5 capsid, a finding reminiscent of observations in infectious viruses. More generally, we found that retrotransposon-derived proteins vary in their intrinsically disordered region content and this is directly associated with their evolutionary rates. Most positively selected sites in these proteins are located in intrinsically disordered regions and some of them impact protein posttranslational modifications, such as autocleavage and phosphorylation. Detailed analyses of the biophysical properties of intrinsically disordered regions showed that positive selection preferentially targeted regions with lower conformational entropy. Furthermore, positive selection introduces variation in binary sequence patterns across orthologues, as well as in chain compaction. Our results shed light on the evolutionary trajectories of a unique class of mammalian genes and suggest a novel approach to study how intrinsically disordered region biophysical characteristics are affected by evolution.
Assuntos
Evolução Molecular , Retroelementos , Animais , Seleção Genética , Mamíferos/genética , Mamíferos/virologia , Proteínas Intrinsicamente Desordenadas/genética , Mudança da Fase de Leitura do Gene Ribossômico , HumanosRESUMO
Laccases (EC 1.10.3.2) are multicopper oxidases with the capability to oxidize diverse phenolic and non-phenolic substrates. While the molecular mechanism of their activity towards phenolic substrates is well-established, their reactivity towards non-phenolic substrates, such as polycyclic aromatic hydrocarbons (PAHs), remains unclear. To elucidate the oxidation mechanism of PAHs, particularly the activation mechanism of the sp2 aromatic C-H bond, we conducted a density functional theory investigation on the oxidation of two PAHs (anthracene and benzo[a]pyrene) using an extensive model of the T1 copper catalytic site of the fungal laccase from Trametes versicolor.
Assuntos
Antracenos , Benzo(a)pireno , Cobre , Lacase , Oxirredução , Lacase/metabolismo , Lacase/química , Antracenos/química , Antracenos/metabolismo , Cobre/química , Cobre/metabolismo , Benzo(a)pireno/metabolismo , Benzo(a)pireno/química , Teoria da Densidade Funcional , Modelos Moleculares , Polyporaceae/enzimologia , Domínio Catalítico , Polyporales/enzimologia , Polyporales/metabolismo , Trametes/enzimologiaRESUMO
The novel dinuclear complex related to the [FeFe]-hydrogenases active site, [Fe2(µ-pdt)(κ2-dmpe)2(CO)2] (1), is highly reactive toward chlorinated compounds CHxCl4-x (x = 1, 2) affording selectively terminal or bridging chloro diiron isomers through a C-Cl bond activation. DFT calculations suggest a cooperative mechanism involving a formal concerted regioselective chloronium transfer depending on the unrotated or rotated conformation of two isomers of 1.
Assuntos
Hidrogenase , Proteínas Ferro-Enxofre , Hidrogenase/química , Proteínas Ferro-Enxofre/química , Isomerismo , Domínio Catalítico , Teoria da Densidade FuncionalRESUMO
The salt [K(18-crown-6)]2[Ru(CN)2(CO)3] ([K(18-crown-6)]2[1]) was generated by the reaction of Ru(C2H4)(CO)4 with [K(18-crown-6)]CN. An initial thermal reaction gives [Ru(CN)(CO)4]-, which, upon ultraviolet (UV) irradiation, reacts with a second equiv of CN-. Protonation of [1]2- gave [HRu(CN)2(CO)3]- ([H1]-), which was isolated as a single isomer with mutually trans cyanide ligands. The complex cis,cis,cis-[Ru(pdt)(CN)2(CO)2]2- ([2]2-) was prepared by the UV-induced reaction of [1]2- with propanedithiol (pdtH2). The corresponding iron complex cis,cis,cis-[Fe(pdt)(CN)2(CO)2]2- ([3]2-) was prepared similarly. The pdt complexes [2]2- and [3]2- were treated with Fe(benzylideneacetone)(CO)3 to give, respectively, [RuFe (µ-pdt)(CN)2(CO)4]2- ([5]2-) and [Fe2(µ-pdt)(CN)2(CO)4]2- ([4]2-). The pathway from [3]2- to Fe2 complex [4]2- implicates intermetallic migration of CN-. In contrast, the formation of [5]2- leaves the Ru(CN)2(CO) center intact, as confirmed by X-ray crystallography. The structure of [5]2- features a "rotated" square-pyramidal Fe(CO)2(µ-CO) site. NMR measurements indicate that the octahedral Ru site is stereochemically rigid, whereas the Fe site dynamically undergoes turnstile rotation. 57Fe Mössbauer spectral parameters are very similar for rotated [5]2- and unrotated Fe2 complex [4]2-, indicating the insensitivity of that technique to both the geometry and the oxidation state of the Fe site. According to cyclic voltammetry, [5]2- oxidizes at E1/2 â¼ -0.8 V vs Fc+/0. Electron paramagnetic resonance (EPR) measurements show that 1e- oxidation of [5]2- gives an S = 1/2 rhombic species, consistent with the formulation Ru(II)Fe(I), related to the Hox state of the [FeFe] hydrogenases. Density functional theory (DFT) studies reproduce the structure, 1H NMR shifts, and infrared (IR) spectra observed for [5]2-. Related homometallic complexes with both cyanides on a single metal are predicted to not adopt rotated structures. These data suggest that [5]2- is best described as Ru(II)Fe(0). This conclusion raises the possibility that for some reduced states of the [FeFe]-hydrogenases, the [2Fe]H site may be better described as Fe(II)Fe(0) than Fe(I)Fe(I).
RESUMO
Flavodoxins are enzymes that contain the redox-active flavin mononucleotide (FMN) cofactor and play a crucial role in numerous biological processes, including energy conversion and electron transfer. Since the redox characteristics of flavodoxins are significantly impacted by the molecular environment of the FMN cofactor, the evaluation of the interplay between the redox properties of the flavin cofactor and its molecular surroundings in flavoproteins is a critical area of investigation for both fundamental research and technological advancements, as the electrochemical tuning of flavoproteins is necessary for optimal interaction with redox acceptor or donor molecules. In order to facilitate the rational design of biomolecular devices, it is imperative to have access to computational tools that can accurately predict the redox potential of both natural and artificial flavoproteins. In this study, we have investigated the feasibility of using non-equilibrium thermodynamic integration protocols to reliably predict the redox potential of flavodoxins. Using as a test set the wild-type flavodoxin from Clostridium Beijerinckii and eight experimentally characterized single-point mutants, we have computed their redox potential. Our results show that 75% (6 out of 8) of the calculated reaction free energies are within 1 kcal/mol of the experimental values, and none exceed an error of 2 kcal/mol, confirming that non-equilibrium thermodynamic integration is a trustworthy tool for the quantitative estimation of the redox potential of this biologically and technologically significant class of enzymes.
Assuntos
Clostridium beijerinckii , Flavodoxina , Termodinâmica , Flavoproteínas , Transporte de ElétronsRESUMO
In endemic regions (West Africa and the Congo Basin), the genetic diversity of monkeypox virus (MPXV) is geographically structured into two major clades (Clades I and II) that differ in virulence and host associations. Clade IIb is closely related to the B.1 lineage, which is dominating a worldwide outbreak initiated in 2022. Lineage B.1 has however accumulated mutations of unknown significance that most likely result from apolipoprotein B mRNA editing catalytic polypeptide-like 3 (APOBEC3) editing. We applied a population genetics-phylogenetics approach to investigate the evolution of MPXV during historical viral spread in Africa and to infer the distribution of fitness effects. We observed a high preponderance of codons evolving under strong purifying selection, particularly in viral genes involved in morphogenesis and replication or transcription. However, signals of positive selection were also detected and were enriched in genes involved in immunomodulation and/or virulence. In particular, several genes showing evidence of positive selection were found to hijack different steps of the cellular pathway that senses cytosolic DNA. Also, a few selected sites in genes that are not directly involved in immunomodulation are suggestive of antibody escape or other immune-mediated pressures. Because orthopoxvirus host range is primarily determined by the interaction with the host immune system, we suggest that the positive selection signals represent signatures of host adaptation and contribute to the different virulence of Clade I and II MPXVs. We also used the calculated selection coefficients to infer the effects of mutations that define the predominant human MPXV1 (hMPXV1) lineage B.1, as well as the changes that have been accumulating during the worldwide outbreak. Results indicated that a proportion of deleterious mutations were purged from the predominant outbreak lineage, whose spread was not driven by the presence of beneficial changes. Polymorphic mutations with a predicted beneficial effect on fitness are few and have a low frequency. It remains to be determined whether they have any significance for ongoing virus evolution.
RESUMO
Nitrogen gas is a highly inert molecule and its activation under mild conditions represents a crucial goal in current research. In a recent study, the discovery of low-valence Ca(I) compounds capable of coordinating and reducing N2 was reported [B. Rösch, T. X. Gentner, J. Langer, C. Färber, J. Eyselein, L. Zhao, C. Ding, G. Frenking and S. Harder, Science, 2021, 371, 1125]. The study of low-valence alkaline earth complexes represents a new horizon in inorganic chemistry and demonstrates examples of spectacular reactivity. For example, complexes of the [BDI]2Mg2 type are selective reducing reagents in both organic and inorganic synthesis reactions. To date, however, no activity of Mg(I) complexes in the activation of the nitrogen molecule has been reported. By computational studies, in the present work, we investigated the analogies and differences of low-valence Ca(I) and Mg(I) complexes in the coordination, activation and protonation of N2. We have shown that the possibility of alkaline earth metals to employ atomic orbitals of the d type is reflected in the differences in the N2 binding energy and its coordination mode (end-on vs. side-on), as well as in the spin state of the resulting adduct (singlet vs. triplet). These divergences are finally observed in the subsequent protonation reaction, which turned out to be prohibitive in the presence of Mg.
RESUMO
Three hexacarbonyl diiron dithiolate complexes [Fe2 (CO)6 (µ-(SCH2 )2 X)] with different substituted bridgeheads (X=CH2 , CEt2 , CBn2 (Bn=CH2 C6 H5 )), have been studied under the same experimental conditions by cyclic voltammetry in dichloromethane [NBu4 ][PF6 ] 0.2â M. DFT calculations were performed to rationalize the mechanism of reduction of these compounds. The three complexes undergo a two-electron transfer whose the mechanism depends on the bulkiness of the dithiolate bridge, which involves a different timing of the structural changes (Fe-S bond cleavage, inversion of conformation and CO bridging) vs redox steps. The introduction of a bulky group in the dithiolate linker has obviously an effect on normally ordered (as for propanedithiolate (pdt)) or inverted (pdtEt2 , pdtBn2 ) reduction potentials. EtâBn replacement is not theoretically predicted to alter the geometry and energy of the most stable mono-reduced and bi-reduced forms but such a replacement alters the kinetics of the electron transfer vs the structural changes.
Assuntos
Hidrogenase , Proteínas Ferro-Enxofre , Hidrogenase/química , Proteínas Ferro-Enxofre/química , Biomimética , Transporte de Elétrons , OxirreduçãoRESUMO
Molecular modeling techniques have become indispensable in many fields of molecular sciences in which the details related to mechanisms and reactivity need to be studied at an atomistic level. This review article provides a collection of computational modeling works on a topic of enormous interest and urgent relevance: the properties of metalloenzymes involved in the degradation and valorization of natural biopolymers and synthetic plastics on the basis of both circular biofuel production and bioremediation strategies. In particular, we will focus on lytic polysaccharide monooxygenase, laccases, and various heme peroxidases involved in the processing of polysaccharides, lignins, rubbers, and some synthetic polymers. Special attention will be dedicated to the interaction between these enzymes and their substrate studied at different levels of theory, starting from classical molecular docking and molecular dynamics techniques up to techniques based on quantum chemistry.
Assuntos
Plásticos , Polissacarídeos , Plásticos/metabolismo , Simulação de Acoplamento Molecular , Oxirredução , Polissacarídeos/metabolismo , Lignina/metabolismo , Estresse Oxidativo , Biopolímeros/metabolismoRESUMO
De novo metalloprotein design is a remarkable approach to shape protein scaffolds toward specific functions. Here, we report the design and characterization of Due Rameâ 1 (DR1), a de novo designed protein housing a di-copper site and mimicking the Typeâ 3 (T3) copper-containing polyphenol oxidases (PPOs). To achieve this goal, we hierarchically designed the first and the second di-metal coordination spheres to engineer the di-copper site into a simple four-helix bundle scaffold. Spectroscopic, thermodynamic, and functional characterization revealed that DR1 recapitulates the T3 copper site, supporting different copper redox states, and being active in the O2 -dependent oxidation of catechols to o-quinones. Careful design of the residues lining the substrate access site endows DR1 with substrate recognition, as revealed by Hammet analysis and computational studies on substituted catechols. This study represents a premier example in the construction of a functional T3 copper site into a designed four-helix bundle protein.
Assuntos
Cobre , Metaloproteínas , Cobre/química , Catecóis/química , Metaloproteínas/química , OxirreduçãoRESUMO
This work originated from the need to functionalize surfactant-coated inorganic nanoparticles for biomedical applications, a process that is limited by excess unbound surfactant. These limitations are connected to the bioconjugation of targeting molecules that are often in equilibrium between the free aliquot in solution and that which binds the surface of the nanoparticles. The excess in solution can play a role in the biocompatability in vitro and in vivo of the final nanoparticles stock. For this purpose, we tested the ability of common surfactants - monothiolated polyethylene glycol and amphiphilic polymers - to colloidally stabilize nanoparticles as excess surfactant is removed and compared them to newly appearing multidentate surfactants endowed with high avidity for inorganic nanoparticles. Our results showed that monothiolated polyethylene glycol or amphiphilic polymers have an insufficient affinity to the nanoparticles and as the excess surfactant is removed the colloidal stability is lost, while multidentate high-avidity surfactants excel in the same regard, possibly allowing improvement in an array of nanoparticle applications, especially in those stated.
Assuntos
Nanopartículas Metálicas , Surfactantes Pulmonares , Tensoativos , Ouro , Polietilenoglicóis , PolímerosRESUMO
Autosomal dominant sleep-related hypermotor epilepsy (ADSHE) is the familial form of a focal epilepsy characterized by hyperkinetic focal seizures, mainly arising during non-rapid eye movements (NREM) sleep. Mutations associated with ADSHE account for a small proportion of the genetically determined cases, suggesting the existence of other disease-causing genes. Here, we reported the results obtained by performing trio-based whole-exome sequencing (WES) in an Italian family showing ADSHE and investigated the structural impact of putative variants by in silico modeling analysis. We identified a p.(Trp276Gly) variant in MOXD1 gene encoding the monooxigenase DBH like 1 protein, cosegregating with the disease and annotated as VUS under the ACMG recommendations. Structural bioinformatic analysis predicted a high destabilizing effect of this variant, due to the loss of important hydrophilic bonds and an expansion of cavity volume in the protein hydrophobic core. Although our data support a functional effect of the p.(Trp276Gly) variant, we highlight the need to identify additional families carrying MOXD1 mutations or functional analyses in suitable models to clarify its role in ADSHE pathogenesis. Moreover, we discuss the importance of VUS reporting due to the low rate of pathogenic variant identification by NGS in epilepsy and for future reinterpretation studies.
Assuntos
Epilepsia , Exoma , Humanos , Sequenciamento do Exoma/métodos , Mutação , LinhagemRESUMO
In the structure from motion, the viewing graph is a graph where the vertices correspond to cameras (or images) and the edges represent the fundamental matrices. We provide a new formulation and an algorithm for determining whether a viewing graph is solvable, i.e., uniquely determines a set of projective cameras. The known theoretical conditions either do not fully characterize the solvability of all viewing graphs, or are extremely difficult to compute because they involve solving a system of polynomial equations with a large number of unknowns. The main result of this paper is a method to reduce the number of unknowns by exploiting cycle consistency. We advance the understanding of solvability by (i) finishing the classification of all minimal graphs up to 9 nodes, (ii) extending the practical verification of solvability to minimal graphs with up to 90 nodes, (iii) finally answering an open research question by showing that finite solvability is not equivalent to solvability, and (iv) formally drawing the connection with the calibrated case (i.e., parallel rigidity). Finally, we present an experiment on real data that shows that unsolvable graphs may appear in practice.
RESUMO
Determining the redox potentials of protein cofactors and how they are influenced by their molecular neighborhoods is essential for basic research and many biotechnological applications, from biosensors and biocatalysis to bioremediation and bioelectronics. The laborious determination of redox potential with current experimental technologies pushes forward the need for computational approaches that can reliably predict it. Although current computational approaches based on quantum and molecular mechanics are accurate, their large computational costs hinder their usage. In this work, we explored the possibility of using more efficient QSPR models based on machine learning (ML) for the prediction of protein redox potential, as an alternative to classical approaches. As a proof of concept, we focused on flavoproteins, one of the most important families of enzymes directly involved in redox processes. To train and test different ML models, we retrieved a dataset of flavoproteins with a known midpoint redox potential (Em) and 3D structure. The features of interest, accounting for both short- and long-range effects of the protein matrix on the flavin cofactor, have been automatically extracted from each protein PDB file. Our best ML model (XGB) has a performance error below 1 kcal/mol (â¼36 mV), comparing favorably to more sophisticated computational approaches. We also provided indications on the features that mostly affect the Em value, and when possible, we rationalized them on the basis of previous studies.
Assuntos
Flavinas , Flavoproteínas , Flavinas/química , Flavinas/metabolismo , Flavoproteínas/química , Aprendizado de Máquina , OxirreduçãoRESUMO
The behaviour of triazolylidene ligands coordinated at a {Fe2(CO)5(µ-dithiolate)} core related to the active site of [FeFe]-hydrogenases have been considered to determine whether such carbenes may act as redox electron-reservoirs, with innocent or non-innocent properties. A novel complex featuring a mesoionic carbene (MIC) [Fe2(CO)5(Pmpt)(µ-pdt)] (1; Pmpt = 1-phenyl-3-methyl-4-phenyl-1,2,3-triazol-5-ylidene; pdt = propanedithiolate) was synthesized and characterized by IR, 1H, 13C{1H} NMR spectroscopies, elemental analyses, X-ray diffraction ,and cyclic voltammetry. Comparison with the spectroscopic characteristics of its analogue [Fe2(CO)5(Pmbt)(µ-pdt)] (2; Pmbt = 1-phenyl-3-methyl-4-butyl-1,2,3-triazol-5-ylidene) showed the effect of the replacement of a n-butyl by a phenyl group in the 1,2,3-triazole heterocycle. A DFT study was performed to rationalize the electronic behaviour of 1, 2 upon the transfer of two electrons and showed that such carbenes do not behave as redox ligands. With highly perfluorinated carbenes, electronic communication between the di-iron site and the triazole cycle is still limited, suggesting low redox properties of MIC ligands used in this study. Finally, although the catalytic performances of 2 towards proton reduction are weak, the protonation process after a two-electron reduction of 2 was examined by DFT and revealed that the protonation process is favoured by S-protonation but the stabilized diprotonated intermediate featuring a {Fe-Hâ¯H-S} interaction does not facilitate the release of H2 and may explain low efficiency towards HER (Hydrogen Evolution Reaction).
Assuntos
Hidrogenase , Proteínas Ferro-Enxofre , Hidrogenase/química , Ferro/química , Proteínas Ferro-Enxofre/química , Ligantes , Prótons , TriazóisRESUMO
Despite the high levels of interest in the synthesis of bio-inspired [FeFe]-hydrogenase complexes, H2 oxidation, which is one specific aspect of hydrogenase enzymatic activity, is not observed for most reported complexes. To attempt H-H bond cleavage, two disubstituted diiron dithiolate complexes in the form of [Fe2(µ-pdt)L2(CO)4] (L: PMe3, dmpe) have been used to play the non-biomimetic role of a Lewis base, with frustrated Lewis pairs (FLPs) formed in the presence of B(C6F5)3 Lewis acid. These unprecedented FLPs, based on the bimetallic Lewis base partner, allow the heterolytic splitting of the H2 molecule, forming a protonated diiron cation and hydrido-borate anion. The substitution, symmetrical or asymmetrical, of two phosphine ligands at the diiron dithiolate core induces a strong difference in the H2 bond cleavage abilities, with the FLP based on the first complex being more efficient than the second. DFT investigations examined the different mechanistic pathways involving each accessible isomer and rationalized the experimental findings. One of the main DFT results highlights that the iron site acting as a Lewis base for the asymmetrical complex is the {Fe(CO)3} subunit, which is less electron-rich than the {FeL(CO)2} site of the symmetrical complex, diminishing the reactivity towards H2. Calculations relating to the different mechanistic pathways revealed the presence of a terminal hydride intermediate at the apical site of a rotated {Fe(CO)3} site, which is experimentally observed, and a semi-bridging hydride intermediate from H2 activation at the Fe-Fe site; these are responsible for a favourable back-reaction, reducing the conversion yield observed in the case of the asymmetrical complex. The use of two equivalents of Lewis acid allows for more complete and faster H2 bond cleavage due to the encapsulation of the hydrido-borate species by a second borane, favouring the reactivity of each FLP, in agreement with DFT calculations.
RESUMO
Four endemic coronaviruses infect humans and cause mild symptoms. Because previous analyses were based on a limited number of sequences and did not control for effects that affect molecular dating, we re-assessed the timing of endemic coronavirus emergence. After controlling for recombination, selective pressure, and molecular clock model, we obtained similar tMRCA (time to the most recent common ancestor) estimates for the four coronaviruses, ranging from 72 (HCoV-229E) to 54 (HCoV-NL63) years ago. The split times of HCoV-229E and HCoV-OC43 from camel alphacoronavirus and bovine coronavirus were dated ~268 and ~99 years ago. The split times of HCoV-HKU1 and HCoV-NL63 could not be calculated, as their zoonoticic sources are unknown. To compare the timing of coronavirus emergence to that of another respiratory virus, we recorded the occurrence of influenza pandemics since 1500. Although there is no clear relationship between pandemic occurrence and human population size, the frequency of influenza pandemics seems to intensify starting around 1700, which corresponds with the initial phase of exponential increase of human population and to the emergence of HCoV-229E. The frequency of flu pandemics in the 19th century also suggests that the concurrence of HCoV-OC43 emergence and the Russian flu pandemic may be due to chance.
Assuntos
Coronavirus Humano 229E , Infecções por Coronavirus , Coronavirus Humano NL63 , Coronavirus Humano OC43 , Influenza Humana , Animais , Bovinos , Coronavirus Humano 229E/genética , Infecções por Coronavirus/epidemiologia , Coronavirus Humano OC43/genética , Humanos , Fatores de TempoRESUMO
Amyotrophic lateral sclerosis (ALS) is a fatal neurodegenerative disease characterized by the loss of the upper and lower motor neurons (MNs). About 10% of patients have a family history (familial, fALS); however, most patients seem to develop the sporadic form of the disease (sALS). SOD1 (Cu/Zn superoxide dismutase-1) is the first studied gene among the ones related to ALS. Mutant SOD1 can adopt multiple misfolded conformation, lose the correct coordination of metal binding, decrease structural stability, and form aggregates. For all these reasons, it is complicated to characterize the conformational alterations of the ALS-associated mutant SOD1, and how they relate to toxicity. In this work, we performed a multilayered study on fibroblasts derived from two ALS patients, namely SOD1L145F and SOD1S135N, carrying the p.L145F and the p.S135N missense variants, respectively. The patients showed diverse symptoms and disease progression in accordance with our bioinformatic analysis, which predicted the different effects of the two mutations in terms of protein structure. Interestingly, both mutations had an effect on the fibroblast energy metabolisms. However, while the SOD1L145F fibroblasts still relied more on oxidative phosphorylation, the SOD1S135N fibroblasts showed a metabolic shift toward glycolysis. Our study suggests that SOD1 mutations might lead to alterations in the energy metabolism.
RESUMO
Coronaviruses (CoVs) have complex genomes that encode a fixed array of structural and nonstructural components, as well as a variety of accessory proteins that differ even among closely related viruses. Accessory proteins often play a role in the suppression of immune responses and may represent virulence factors. Despite their relevance for CoV phenotypic variability, information on accessory proteins is fragmentary. We applied a systematic approach based on homology detection to create a comprehensive catalogue of accessory proteins encoded by CoVs. Our analyses grouped accessory proteins into 379 orthogroups and 12 super-groups. No orthogroup was shared by the four CoV genera and very few were present in all or most viruses in the same genus, reflecting the dynamic evolution of CoV genomes. We observed differences in the distribution of accessory proteins in CoV genera. Alphacoronaviruses harboured the largest diversity of accessory open reading frames (ORFs), deltacoronaviruses the smallest. However, the average number of accessory proteins per genome was highest in betacoronaviruses. Analysis of the evolutionary history of some orthogroups indicated that the different CoV genera adopted similar evolutionary strategies. Thus, alphacoronaviruses and betacoronaviruses acquired phosphodiesterases and spike-like accessory proteins independently, whereas horizontal gene transfer from reoviruses endowed betacoronaviruses and deltacoronaviruses with fusion-associated small transmembrane (FAST) proteins. Finally, analysis of accessory ORFs in annotated CoV genomes indicated ambiguity in their naming. This complicates cross-communication among researchers and hinders automated searches of large data sets (e.g., PubMed, GenBank). We suggest that orthogroup membership is used together with a naming system to provide information on protein function.
Assuntos
Coronavirus , Sequência de Aminoácidos , Coronavirus/química , Coronavirus/genética , Evolução Molecular , Genoma Viral/genética , Fases de Leitura Aberta/genéticaRESUMO
Four coronaviruses (HCoV-OC43, HCoV-HKU1, HCoV-NL63, and HCoV-229E) are endemic in human populations. All these viruses are seasonal and generate short-term immunity. Like the highly pathogenic coronaviruses, the endemic coronaviruses have zoonotic origins. Thus, understanding the evolutionary dynamics of these human viruses might provide insight into the future trajectories of SARS-CoV-2 evolution. Because the zoonotic sources of HCoV-OC43 and HCoV-229E are known, we applied a population genetics-phylogenetic approach to investigate which selective events accompanied the divergence of these viruses from the animal ones. Results indicated that positive selection drove the evolution of some accessory proteins, as well as of the membrane proteins. However, the spike proteins of both viruses and the hemagglutinin-esterase (HE) of HCoV-OC43 represented the major selection targets. Specifically, for both viruses, most positively selected sites map to the receptor-binding domains (RBDs) and are polymorphic. Molecular dating for the HCoV-229E spike protein indicated that RBD Classes I, II, III, and IV emerged 3-9 years apart. However, since the appearance of Class V (with much higher binding affinity), around 25 years ago, limited genetic diversity accumulated in the RBD. These different time intervals are not fully consistent with the hypothesis that HCoV-229E spike evolution was driven by antigenic drift. An alternative, not mutually exclusive possibility is that strains with higher affinity for the cellular receptor have out-competed strains with lower affinity. The evolution of the HCoV-OC43 spike protein was also suggested to undergo antigenic drift. However, we also found abundant signals of positive selection in HE. Whereas such signals might result from antigenic drift, as well, previous data showing co-evolution of the spike protein with HE suggest that optimization for human cell infection also drove the evolution of this virus. These data provide insight into the possible trajectories of SARS-CoV-2 evolution, especially in case the virus should become endemic.