Purification and characterization of the receptor‐binding domain of SARS‐CoV‐2 spike protein from Escherichia coli

SARS‐CoV‐2 is responsible for a disruptive worldwide viral pandemic, and renders a severe respiratory disease known as COVID‐19. Spike protein of SARS‐CoV‐2 mediates viral entry into host cells by binding ACE2 through the receptor‐binding domain (RBD). RBD is an important target for development of virus inhibitors, neutralizing antibodies, and vaccines. RBD expressed in mammalian cells suffers from low expression yield and high cost. E. coli is a popular host for protein expression, which has the advantage of easy scalability with low cost. However, RBD expressed by E. coli (RBD‐1) lacks the glycosylation, and its antigenic epitopes may not be sufficiently exposed. In the present study, RBD‐1 was expressed by E. coli and purified by a Ni Sepharose Fast Flow column. RBD‐1 was structurally characterized and compared with RBD expressed by the HEK293 cells (RBD‐2). The secondary structure and tertiary structure of RBD‐1 were largely maintained without glycosylation. In particular, the major β‐sheet content of RBD‐1 was almost unaltered. RBD‐1 could strongly bind ACE2 with a dissociation constant (K_D) of 2.98 × 10^–8 M. Thus, RBD‐1 was expected to apply in the vaccine development, screening drugs and virus test kit.     

Solution X‐ray scattering highlights discrepancies in Plasmodium multi‐aminoacyl‐tRNA synthetase complexes

tRip is a tRNA import protein specific to Plasmodium, the causative agent of malaria. In addition to its membrane localization and tRNA trafficking properties, tRip has the capacity to associate with three aminoacyl‐tRNA synthetases (aaRS), the glutamyl‐ (ERS), glutaminyl‐ (QRS), and methionyl‐ (MRS) tRNA synthetases. In eukaryotes, such multi‐aaRSs complexes (MSC) regulate the moonlighting activities of aaRSs. In Plasmodium, tRip and the three aaRSs all contain an N‐terminal GST‐like domain involved in the assembly of two independent complexes: the Q‐complex (tRip:ERS:QRS) and the M‐complex (tRip:ERS:MRS) with a 2:2:2 stoichiometry and in which the association of the GST‐like domains of tRip and ERS (tRip‐N:ERS‐N) is central. In this study, the crystal structure of the N‐terminal GST‐like domain of ERS was solved and made possible further investigation of the solution architecture of the Q‐ and M‐complexes by small‐angle x‐ray scattering (SAXS). This strategy relied on the engineering of a tRip‐N‐ERS‐N chimeric protein to study the structural scaffold of both Plasmodium MSCs and confirm the unique homodimerization pattern of tRip in solution. The biological impact of these structural arrangements is discussed.     

Cryo‐EM structure of the CENP‐A nucleosome in complex with phosphorylated CENP‐C

The CENP‐A nucleosome is a key structure for kinetochore assembly. Once the CENP‐A nucleosome is established in the centromere, additional proteins recognize the CENP‐A nucleosome to form a kinetochore. CENP‐C and CENP‐N are CENP‐A binding proteins. We previously demonstrated that vertebrate CENP‐C binding to the CENP‐A nucleosome is regulated by CDK1‐mediated CENP‐C phosphorylation. However, it is still unknown how the phosphorylation of CENP‐C regulates its binding to CENP‐A. It is also not completely understood how and whether CENP‐C and CENP‐N act together on the CENP‐A nucleosome. Here, using cryo‐electron microscopy (cryo‐EM) in combination with biochemical approaches, we reveal a stable CENP‐A nucleosome‐binding mode of CENP‐C through unique regions. The chicken CENP‐C structure bound to the CENP‐A nucleosome is stabilized by an intramolecular link through the phosphorylated CENP‐C residue. The stable CENP‐A‐CENP‐C complex excludes CENP‐N from the CENP‐A nucleosome. These findings provide mechanistic insights into the dynamic kinetochore assembly regulated by CDK1‐mediated CENP‐C phosphorylation.     

m6A‐mediated alternative splicing coupled with nonsense‐mediated mRNA decay regulates SAM synthetase homeostasis

Alternative splicing of pre‐mRNAs can regulate gene expression levels by coupling with nonsense‐mediated mRNA decay (NMD). In order to elucidate a repertoire of mRNAs regulated by alternative splicing coupled with NMD (AS‐NMD) in an organism, we performed long‐read RNA sequencing of poly(A)⁺ RNAs from an NMD‐deficient mutant strain of Caenorhabditis elegans, and obtained full‐length sequences for mRNA isoforms from 259 high‐confidence AS‐NMD genes. Among them are the S‐adenosyl‐L‐methionine (SAM) synthetase (sams) genes sams‐3 and sams‐4. SAM synthetase activity autoregulates sams gene expression through AS‐NMD in a negative feedback loop. We furthermore find that METT‐10, the orthologue of human U6 snRNA methyltransferase METTL16, is required for the splicing regulation in␣vivo, and specifically methylates the invariant AG dinucleotide at the distal 3′ splice site (3′SS) in␣vitro. Direct RNA sequencing coupled with machine learning confirms m⁶A modification of endogenous sams mRNAs. Overall, these results indicate that homeostasis of SAM synthetase in C. elegans is maintained by alternative splicing regulation through m⁶A modification at the 3′SS of the sams genes.     

Porin threading drives receptor disengagement and establishes active colicin transport through Escherichia coli OmpF

Bacteria deploy weapons to kill their neighbours during competition for resources and to aid survival within microbiomes. Colicins were the first such antibacterial system identified, yet how these bacteriocins cross the outer membrane (OM) of Escherichia coli is unknown. Here, by solving the structures of translocation intermediates via cryo‐EM and by imaging toxin import, we uncover the mechanism by which the Tol‐dependent nuclease colicin E9 (ColE9) crosses the bacterial OM. We show that threading of ColE9’s disordered N‐terminal domain through two pores of the trimeric porin OmpF causes the colicin to disengage from its primary receptor, BtuB, and reorganises the translocon either side of the membrane. Subsequent import of ColE9 through the lumen of a single OmpF subunit is driven by the proton‐motive force, which is delivered by the TolQ‐TolR‐TolA‐TolB assembly. Our study answers longstanding questions, such as why OmpF is a better translocator than OmpC, and reconciles the mechanisms by which both Tol‐ and Ton‐dependent bacteriocins cross the bacterial outer membrane.     

Protein yoga: Conformational versatility of the Hemolysin II C‐terminal domain detailed by NMR structures for multiple states

The C‐terminal domain of Bacillus cereus hemolysin II (HlyIIC), stabilizes the trans‐membrane‐pore formed by the HlyII toxin and may aid in target cell recognition. Initial efforts to determine the NMR structure of HlyIIC were hampered by cis/trans isomerization about the single proline at position 405 that leads to doubling of NMR resonances. We used the mutant P405M‐HlyIIC that eliminates the cis proline to determine the NMR structure of the domain, which revealed a novel fold. Here, we extend earlier studies to the NMR structure determination of the cis and trans states of WT‐HlyIIC that exist simultaneously in solution. The primary structural differences between the cis and trans states are in the loop that contains P405, and structurally adjacent loops. Thermodynamic linkage analysis shows that at 25 C the cis proline, which already has a large fraction of 20% in the unfolded protein, increases to 50% in the folded state due to coupling with the global stability of the domain. The P405M or P405A substitutions eliminate heterogeneity due to proline isomerization but lead to the formation of a new dimeric species. The NMR structure of the dimer shows that it is formed through domain‐swapping of strand β5, the last segment of secondary structure following P405. The presence of P405 in WT‐HlyIIC strongly disfavors the dimer compared to the P405M‐HlyIIC or P405A‐HlyIIC mutants. The WT proline may thus act as a “gatekeeper,” warding off aggregative misfolding.     

Characterization of an aminotransferase from Acanthamoeba polyphaga Mimivirus

Acanthamoeba polyphaga Mimivirus, a complex virus that infects amoeba, was first reported in 2003. It is now known that its DNA genome encodes for nearly 1,000 proteins including enzymes that are required for the biosynthesis of the unusual sugar 4‐amino‐4,6‐dideoxy‐d‐glucose, also known as d‐viosamine. As observed in some bacteria, the pathway for the production of this sugar initiates with a nucleotide‐linked sugar, which in the Mimivirus is thought to be UDP‐d‐glucose. The enzyme required for the installment of the amino group at the C‐4′ position of the pyranosyl moiety is encoded in the Mimivirus by the L136 gene. Here, we describe a structural and functional analysis of this pyridoxal 5′‐phosphate‐dependent enzyme, referred to as L136. For this analysis, three high‐resolution X‐ray structures were determined: the wildtype enzyme/pyridoxamine 5′‐phosphate/dTDP complex and the site‐directed mutant variant K185A in the presence of either UDP‐4‐amino‐4,6‐dideoxy‐d‐glucose or dTDP‐4‐amino‐4,6‐dideoxy‐d‐glucose. Additionally, the kinetic parameters of the enzyme utilizing either UDP‐d‐glucose or dTDP‐d‐glucose were measured and demonstrated that L136 is efficient with both substrates. This is in sharp contrast to the structurally related DesI from Streptomyces venezuelae, whose three‐dimensional architecture was previously reported by this laboratory. As determined in this investigation,DesI shows a profound preference in its catalytic efficiency for the dTDP‐linked sugar substrate. This difference can be explained in part by a hydrophobic patch in DesI that is missing in L136. Notably, the structure of L136 reported here represents the first three‐dimensional model for a virally encoded PLP‐dependent enzyme and thus provides new information on sugar aminotransferases in general.     

Cryo‐EM structures of perforin‐2 in isolation and assembled on a membrane suggest a mechanism for pore formation

Perforin‐2 (PFN2, MPEG1) is a key pore‐forming protein in mammalian innate immunity restricting intracellular bacteria proliferation. It forms a membrane‐bound pre‐pore complex that converts to a pore‐forming structure upon acidification; but its mechanism of conformational transition has been debated. Here we used cryo‐electron microscopy, tomography and subtomogram averaging to determine structures of PFN2 in pre‐pore and pore conformations in isolation and bound to liposomes. In isolation and upon acidification, the pre‐assembled complete pre‐pore rings convert to pores in both flat ring and twisted conformations. On membranes, in situ assembled PFN2 pre‐pores display various degrees of completeness; whereas PFN2 pores are mainly incomplete arc structures that follow the same subunit packing arrangements as found in isolation. Both assemblies on membranes use their P2 β‐hairpin for binding to the lipid membrane surface. Overall, these structural snapshots suggest a molecular mechanism for PFN2 pre‐pore to pore transition on a targeted membrane, potentially using the twisted pore as an intermediate or alternative state to the flat conformation, with the capacity to cause bilayer distortion during membrane insertion.     

Site‐directed mutagenesis of Serine‐72 reveals the location of the fructose 6‐phosphate regulatory site of the Agrobacterium tumefaciens ADP‐glucose pyrophosphorylase

The allosteric regulation of ADP–glucose pyrophosphorylase is critical for the biosynthesis of glycogen in bacteria and starch in plants. The enzyme from Agrobacterium tumefaciens is activated by fructose 6‐phosphate (Fru6P) and pyruvate (Pyr). The Pyr site has been recently found, but the site where Fru6P binds has remained unknown. We hypothesize that a sulfate ion previously found in the crystal structure reveals a part of the regulatory site mimicking the presence of the phosphoryl moiety of the activator Fru6P. Ser72 interacts with this sulfate ion and, if the hypothesis is correct, Ser72 would affect the interaction with Fru6P and activation of the enzyme. Here, we report structural, binding, and kinetic analysis of Ser72 mutants of the A. tumefaciens ADP‐glucose pyrophosphorylase. By X‐ray crystallography, we found that when Ser72 was replaced by Asp or Glu side chain carboxylates protruded into the sulfate‐binding pocket. They would present a strong steric and electrostatic hindrance to the phosphoryl moiety of Fru6P, while being remote from the Pyr site. In agreement, we found that Fru6P could not activate or bind to S72E or S72D mutants, whereas Pyr was still an effective activator. These mutants also blocked the binding of the inhibitor AMP. This could potentially have biotechnological importance in obtaining enzyme forms insensitive to inhibition. Other mutations in this position (Ala, Cys, and Trp) confirmed the importance of Ser72 in regulation. We propose that the ADP‐glucose pyrophosphorylase from A. tumefaciens have two distinct sites for Fru6P and Pyr working in tandem to regulate glycogen biosynthesis.     

Concordant geographic and genetic structure revealed by genotyping‐by‐sequencing in a New Zealand marine isopod

Population genetic structure in the marine environment can be influenced by life‐history traits such as developmental mode (biphasic, with distinct adult and larval morphology, and direct development, in which larvae resemble adults) or habitat specificity, as well as geography and selection. Developmental mode is thought to significantly influence dispersal, with direct developers expected to have much lower dispersal potential. However, this prediction can be complicated by the presence of geophysical barriers to dispersal. In this study, we use a panel of 8,020 SNPs to investigate population structure and biogeography over multiple spatial scales for a direct‐developing species, the New Zealand endemic marine isopod Isocladus armatus. Because our sampling range is intersected by two well‐known biogeographic barriers (the East Cape and the Cook Strait), our study provides an opportunity to understand how such barriers influence dispersal in direct developers. On a small spatial scale (20 km), gene flow between locations is extremely high, suggestive of an island model of migration. However, over larger spatial scales (600 km), populations exhibit a clear pattern of isolation‐by‐distance. Our results indicate that I. armatus exhibits significant migration across the hypothesized barriers and suggest that large‐scale ocean currents associated with these locations do not present a barrier to dispersal. Interestingly, we find evidence of a north‐south population genetic break occurring between Māhia and Wellington. While no known geophysical barrier is apparent in this area, it coincides with the location of a proposed border between bioregions. Analysis of loci under selection revealed that both isolation‐by‐distance and adaption may be contributing to the degree of population structure we have observed here. We conclude that developmental life history largely predicts dispersal in the intertidal isopod I. armatus. However, localized biogeographic processes can disrupt this expectation, and this may explain the potential meta‐population detected in the Auckland region.     

FastD: Fast detection of insecticide target‐site mutations and overexpressed detoxification genes in insect populations from RNA‐Seq data

Target‐site mutations and detoxification gene overexpression are two major mechanisms conferring insecticide resistance. Molecular assays applied to detect these resistance genetic markers are time‐consuming and with high false‐positive rates. RNA‐Seq data contains information on the variations within expressed genomic regions and expression of detoxification genes. However, there is no corresponding method to detect resistance markers at present. Here, we collected 66 reported resistance mutations of four insecticide targets (AChE, VGSC, RyR, and nAChR) from 82 insect species. Next, we obtained 403 sequences of the four target genes and 12,665 sequences of three kinds of detoxification genes including P450s, GSTs, and CCEs. Then, we developed a Perl program, FastD, to detect target‐site mutations and overexpressed detoxification genes from RNA‐Seq data and constructed a web server for FastD (http://www.insect-genome.com/fastd). The estimation of FastD on simulated RNA‐Seq data showed high sensitivity and specificity. We applied FastD to detect resistant markers in 15 populations of six insects, Plutella xylostella, Aphis gossypii, Anopheles arabiensis, Musca domestica, Leptinotarsa decemlineata and Apis mellifera. Results showed that 11 RyR mutations in P. xylostella, one nAChR mutation in A. gossypii, one VGSC mutation in A. arabiensis and five VGSC mutations in M. domestica were found to be with frequency difference >40% between resistant and susceptible populations including previously confirmed mutations G4946E in RyR, R81T in nAChR and L1014F in VGSC. And 49 detoxification genes were found to be overexpressed in resistant populations compared with susceptible populations including previously confirmed detoxification genes CYP6BG1, CYP6CY22, CYP6CY13, CYP6P3, CYP6M2, CYP6P4 and CYP4G16. The candidate target‐site mutations and detoxification genes were worth further validation. Resistance estimates according to confirmed markers were consistent with population phenotypes, confirming the reliability of this program in predicting population resistance at omics‐level.     

A heterotypic assembly mechanism regulates CHIP E3 ligase activity

CHIP (C‐terminus of Hsc70‐interacting protein) and its worm ortholog CHN‐1 are E3 ubiquitin ligases that link the chaperone system with the ubiquitin‐proteasome system (UPS). CHN‐1 can cooperate with UFD‐2, another E3 ligase, to accelerate ubiquitin chain formation; however, the basis for the high processivity of this E3s set has remained obscure. Here, we studied the molecular mechanism and function of the CHN‐1–UFD‐2 complex in Caenorhabditis elegans. Our data show that UFD‐2 binding promotes the cooperation between CHN‐1 and ubiquitin‐conjugating E2 enzymes by stabilizing the CHN‐1 U‐box dimer. However, HSP70/HSP‐1 chaperone outcompetes UFD‐2 for CHN‐1 binding, thereby promoting a shift to the autoinhibited CHN‐1 state by acting on a conserved residue in its U‐box domain. The interaction with UFD‐2 enables CHN‐1 to efficiently ubiquitylate and regulate S‐adenosylhomocysteinase (AHCY‐1), a key enzyme in the S‐adenosylmethionine (SAM) regeneration cycle, which is essential for SAM‐dependent methylation. Our results define the molecular mechanism underlying the synergistic cooperation of CHN‐1 and UFD‐2 in substrate ubiquitylation.     

Structure of cytosine transport protein CodB provides insight into nucleobase‐cation symporter 1 mechanism

CodB is a cytosine transporter from the Nucleobase‐Cation‐Symport‐1 (NCS1) transporter family, a member of the widespread LeuT superfamily. Previous experiments with the nosocomial pathogen Pseudomonas aeruginosa have shown CodB as also important for the uptake of 5‐fluorocytosine, which has been suggested as a novel drug to combat antimicrobial resistance by suppressing virulence. Here we solve the crystal structure of CodB from Proteus vulgaris, at 2.4 Å resolution in complex with cytosine. We show that CodB carries out the sodium‐dependent uptake of cytosine and can bind 5‐fluorocytosine. Comparison of the substrate‐bound structures of CodB and the hydantoin transporter Mhp1, the only other NCS1 family member for which the structure is known, highlight the importance of the hydrogen bonds that the substrates make with the main chain at the breakpoint in the discontinuous helix, TM6. In contrast to other LeuT superfamily members, neither CodB nor Mhp1 makes specific interactions with residues on TM1. Comparison of the structures provides insight into the intricate mechanisms of how these proteins transport substrates across the plasma membrane.     

Plasticity is a locally adapted trait with consequences for ecological dynamics in novel environments

Phenotypic plasticity is predicted to evolve in more variable environments, conferring an advantage on individual lifetime fitness. It is less clear what the potential consequences of that plasticity will have on ecological population dynamics. Here, we use an invertebrate model system to examine the effects of environmental variation (resource availability) on the evolution of phenotypic plasticity in two life history traits—age and size at maturation—in long‐running, experimental density‐dependent environments. Specifically, we then explore the feedback from evolution of life history plasticity to subsequent ecological dynamics in novel conditions. Plasticity in both traits initially declined in all microcosm environments, but then evolved increased plasticity for age‐at‐maturation, significantly so in more environmentally variable environments. We also demonstrate how plasticity affects ecological dynamics by creating founder populations of different plastic phenotypes into new microcosms that had either familiar or novel environments. Populations originating from periodically variable environments that had evolved greatest plasticity had lowest variability in population size when introduced to novel environments than those from constant or random environments. This suggests that while plasticity may be costly it can confer benefits by reducing the likelihood that offspring will experience low survival through competitive bottlenecks in variable environments. In this study, we demonstrate how plasticity evolves in response to environmental variation and can alter population dynamics—demonstrating an eco‐evolutionary feedback loop in a complex animal moderated by plasticity in growth.