Middle East respiratory syndrome coronavirus (MERS-CoV) originated in bats and spread to humans via zoonotic transmission from camels. We analyzed the evolution of the spike (S) gene in betacoronaviruses (betaCoVs) isolated from different mammals, in bat coronavirus populations, as well as in MERS-CoV strains from the current outbreak. Results indicated several positively selected sites located in the region comprising the two heptad repeats (HR1 and HR2) and their linker. Two sites (R652 and V1060) were positively selected in the betaCoVs phylogeny and correspond to mutations associated with expanded host range in other coronaviruses. During the most recent evolution of MERS-CoV, adaptive mutations in the HR1 (Q/R/H1020) arose in camels or in a previous host and spread to humans. We determined that different residues at position 1020 establish distinct inter- and intra-helical interactions and affect the stability of the six-helix bundle formed by the HRs. A similar effect on stability was observed for a nearby mutation (T1015N) that increases MERS-CoV infection efficiency in vitro. Data herein indicate that the heptad repeat region was a major target of adaptive evolution in MERS-CoV-related viruses; these results are relevant for the design of fusion inhibitor peptides with antiviral function.

Density functional theory has been used to study diiron dithiolates [HFe2(xdt)(PR3)n(CO)5-nX] (n = 0, 2, 4; R = H, Me, Et; X = CH3S(-), PMe3, NHC = 1,3-dimethylimidazol-2-ylidene; xdt = adt, pdt; adt = azadithiolate; pdt = propanedithiolate). These species are related to the [FeFe]-hydrogenases catalyzing the 2H(+) + 2e(-) ↔ H2 reaction. Our study is focused on the reduction step following protonation of the Fe2(SR)2 core. Fe(H)s detected in solution are terminal (t-H) and bridging (μ-H) hydrides. Although unstable versus μ-Hs, synthetic t-Hs feature milder reduction potentials than μ-Hs. Accordingly, attempts were previously made to hinder the isomerization of t-H to μ-H. Herein, we present another strategy: in place of preventing isomerization, μ-H could be made a stronger oxidant than t-H (E°μ-H > E°t-H). The nature and number of PR3 unusually affect ΔE°t-H-μ-H: 4PEt3 models feature a μ-H with a milder E° than t-H, whereas the 4PMe3 analogues behave oppositely. The correlation ΔE°t-H-μ-H ↔ stereoelectronic features arises from the steric strain induced by bulky Et groups in 4PEt3 derivatives. One-electron reduction alleviates intramolecular repulsions only in μ-H species, which is reflected in the loss of bridging coordination. Conversely, in t-H, the strain is retained because a bridging CO holds together the Fe2 core. That implies that E°μ-H > E°t-H in 4-PEt3 species but not in 4PMe3 analogues. Also determinant to observe E°μ-H > E°t-H is the presence of a Fe apical σ-donor because its replacement with a CO yields E°μ-H < E°t-H even in 4PEt3 species. Variants with neutral NHC and PMe3 in place of CH3S(-) still feature E°μ-H > E°t-H. Replacing pdt with (Hadt)(+) lowers E° but yields E°μ-H < E°t-H, indicating that μ-H activation can occur to the detriment of the overpotential increase. In conclusion, our results indicate that the electron richness of the Fe2 core influences ΔE°t-H-μ-H, provided that (i) the R size of PR3 must be greater than that of Me and (ii) an electron donor must be bound to Fe apically.

Saccharomyces cerevisiae cells respond to hypotonic stress (HTS) by a cytosolic calcium rise, either generated by an influx of calcium from extracellular medium, when calcium is available, or by a release from intracellular stores in scarcity of extracellular calcium. Calcium release from intracellular compartments is peculiarly inhibited by external calcium in a calcineurin-independent and Cch1-, but not Mid1-, driven manner. HTS-induced calcium release is also negatively regulated by the ER protein Cls2 and involves a poorly characterized protein, FLC2/YAL053W gene product, previously proposed to be required for FAD transport in the ER, albeit, due to its molecular features, it was also previously classified as an ion transporter. A computational analysis revealed that this gene and its three homologs in S. cerevisiae, together with previously identified Schizosaccharomyces pombe pkd2 and Neurospora crassa calcium-related spray protein, belong to a fungal branch of TRP-like ion transporters related to human mucolipin and polycystin 2 calcium transporters. Moreover, disruption of FLC2 gene confers severe sensitivity to Calcofluor white and hyper-activation of the cell wall integrity MAPK cascade, suggesting a role in cell wall maintenance as previously suggested for the fission yeast homolog. Perturbation in cytosolic resting calcium concentration and hyper-activation of calcineurin in exponentially growing cells suggest a role for this transporter in calcium homeostasis in yeast.

The antigenic repertoire presented by MHC molecules is generated by the antigen processing and presentation (APP) pathway. We analyzed the evolutionary history of 45 genes involved in APP at the inter- and intra-species level. Results showed that 11 genes evolved adaptively in mammals. Several positively selected sites involve positions of fundamental importance to the protein function (e.g. the TAP1 peptide-binding domains, the sugar binding interface of langerin, and the CD1D trafficking signal region). In CYBB, all selected sites cluster in two loops protruding into the endosomal lumen; analysis of missense mutations responsible for chronic granulomatous disease (CGD) showed the action of different selective forces on the very same gene region, as most CGD substitutions involve aminoacid positions that are conserved in all mammals. As for ERAP2, different computational methods indicated that positive selection has driven the recurrent appearance of protein-destabilizing variants during mammalian evolution. Application of a population-genetics phylogenetics approach showed that purifying selection represented a major force acting on some APP components (e.g. immunoproteasome subunits and chaperones) and allowed identification of positive selection events in the human lineage.We also investigated the evolutionary history of APP genes in human populations by developing a new approach that uses several different tests to identify the selection target, and that integrates low-coverage whole-genome sequencing data with Sanger sequencing. This analysis revealed that 9 APP genes underwent local adaptation in human populations. Most positive selection targets are located within noncoding regions with regulatory function in myeloid cells or act as expression quantitative trait loci. Conversely, balancing selection targeted nonsynonymous variants in TAP1 and CD207 (langerin). Finally, we suggest that selected variants in PSMB10 and CD207 contribute to human phenotypes. Thus, we used evolutionary information to generate experimentally-testable hypotheses and to provide a list of sites to prioritize in follow-up analyses.

No abstract