Genome-wide RNA structure changes during human neurogenesis modulate gene regulatory networks.

The distribution, dynamics, and function of RNA structures in human development are under-explored. Here, we systematically assayed RNA structural dynamics and their relationship with gene expression, translation, and decay during human neurogenesis. We observed that the human ESC transcriptome is globally more structurally accessible than differentiated cells and undergoes extensive RNA structure changes, particularly in the 3' UTR. Additionally, RNA structure changes during differentiation are associated with translation and decay. We observed that RBP and miRNA binding is associated with RNA structural changes during early neuronal differentiation, and splicing is associated during later neuronal differentiation. Furthermore, our analysis suggests that RBPs are major factors in structure remodeling and co-regulate additional RBPs and miRNAs through structure. We demonstrated an example of this by showing that PUM2-induced structure changes on LIN28A enable miR-30 binding. This study deepens our understanding of the widespread and complex role of RNA-based gene regulation during human development.

Flavivirus genome recoding by codon optimisation confers genetically stable in vivo attenuation in both mice and mosquitoes.

Virus genome recoding is an attenuation method that confers genetically stable attenuation by rewriting a virus genome with numerous silent mutations. Prior flavivirus genome recoding attempts utilised codon deoptimisation approaches. However, these codon deoptimisation approaches act in a species dependent manner and were unable to confer flavivirus attenuation in mosquito cells or in mosquito animal models. To overcome these limitations, we performed flavivirus genome recoding using the contrary approach of codon optimisation. The genomes of flaviviruses such as dengue virus type 2 (DENV2) and Zika virus (ZIKV) contain functional RNA elements that regulate viral replication. We hypothesised that flavivirus genome recoding by codon optimisation would introduce silent mutations that disrupt these RNA elements, leading to decreased replication efficiency and attenuation. We chose DENV2 and ZIKV as representative flaviviruses and recoded them by codon optimising their genomes for human expression. Our study confirms that this recoding approach of codon optimisation does translate into reduced replication efficiency in mammalian, human, and mosquito cells as well as in vivo attenuation in both mice and mosquitoes. In silico modelling and RNA SHAPE analysis confirmed that DENV2 recoding resulted in the extensive disruption of genomic structural elements. Serial passaging of recoded DENV2 resulted in the emergence of rescue or adaptation mutations, but no reversion mutations. These rescue mutations were unable to rescue the delayed replication kinetics and in vivo attenuation of recoded DENV2, demonstrating that recoding confers genetically stable attenuation. Therefore, our recoding approach is a reliable attenuation method with potential applications for developing flavivirus vaccines.

Phosphorylation of the compartmentalized PKA substrate TAF15 regulates RNA-protein interactions.

Spatiotemporal-controlled second messengers alter molecular interactions of central signaling nodes for ensuring physiological signal transmission. One prototypical second messenger molecule which modulates kinase signal transmission is the cyclic-adenosine monophosphate (cAMP). The main proteinogenic cellular effectors of cAMP are compartmentalized protein kinase A (PKA) complexes. Their cell-type specific compositions precisely coordinate substrate phosphorylation and proper signal propagation which is indispensable for numerous cell-type specific functions. Here we present evidence that TAF15, which is implicated in the etiology of amyotrophic lateral sclerosis, represents a novel nuclear PKA substrate. In cross-linking and immunoprecipitation experiments (iCLIP) we showed that TAF15 phosphorylation alters the binding to target transcripts related to mRNA maturation, splicing and protein-binding related functions. TAF15 appears to be one of multiple PKA substrates that undergo RNA-binding dynamics upon phosphorylation. We observed that the activation of the cAMP-PKA signaling axis caused a change in the composition of a collection of RNA species that interact with TAF15. This observation appears to be a broader principle in the regulation of molecular interactions, as we identified a significant enrichment of RNA-binding proteins within endogenous PKA complexes. We assume that phosphorylation of RNA-binding domains adds another layer of regulation to binary protein-RNAs interactions with consequences to RNA features including binding specificities, localization, abundance and composition.

A Reversible Association between Smc Coiled Coils Is Regulated by Lysine Acetylation and Is Required for Cohesin Association with the DNA.

Cohesin is a ring-shaped protein complex that is capable of embracing DNA. Most of the ring circumference is comprised of the anti-parallel intramolecular coiled coils of the Smc1 and Smc3 proteins, which connect globular head and hinge domains. Smc coiled coil arms contain multiple acetylated and ubiquitylated lysines. To investigate the role of these modifications, we substituted lysines for arginines to mimic the unmodified state and uncovered genetic interaction between the Smc arms. Using scanning force microscopy, we show that wild-type Smc arms associate with each other when the complex is not on DNA. Deacetylation of the Smc1/Smc3 dimers promotes arms' dissociation. Smc arginine mutants display loose packing of the Smc arms and, although they dimerize at the hinges, fail to connect the heads and associate with the DNA. Our findings highlight the importance of a "collapsed ring," or "rod," conformation of cohesin for its loading on the chromosomes.

Antibody affinity versus dengue morphology influences neutralization.

Different strains within a dengue serotype (DENV1-4) can have smooth, or "bumpy" surface morphologies with different antigenic characteristics at average body temperature (37°C). We determined the neutralizing properties of a serotype cross-reactive human monoclonal antibody (HMAb) 1C19 for strains with differing morphologies within the DENV1 and DENV2 serotypes. We mapped the 1C19 epitope to E protein domain II by hydrogen deuterium exchange mass spectrometry, cryoEM and molecular dynamics simulations, revealing that this epitope is likely partially hidden on the virus surface. We showed the antibody has high affinity for binding to recombinant DENV1 E proteins compared to those of DENV2, consistent with its strong neutralizing activities for all DENV1 strains tested regardless of their morphologies. This finding suggests that the antibody could out-compete E-to-E interaction for binding to its epitope. In contrast, for DENV2, HMAb 1C19 can only neutralize when the epitope becomes exposed on the bumpy-surfaced particle. Although HMAb 1C19 is not a suitable therapeutic candidate, this study with HMAb 1C19 shows the importance of choosing a high-affinity antibody that could neutralize diverse dengue virus morphologies for therapeutic purposes.

Mutation-oriented profiling of autoinhibitory kinase conformations predicts RAF inhibitor efficacies.

Kinase-targeted therapies have the potential to improve the survival of patients with cancer. However, the cancer-specific spectrum of kinase alterations exhibits distinct functional properties and requires mutation-oriented drug treatments. Besides post-translational modifications and diverse intermolecular interactions of kinases, it is the distinct disease mutation which reshapes full-length kinase conformations, affecting their activity. Oncokinase mutation profiles differ between cancer types, as it was shown for BRAF in melanoma and non-small-cell lung cancers. Here, we present the target-oriented application of a kinase conformation (KinCon) reporter platform for live-cell measurements of autoinhibitory kinase activity states. The bioluminescence-based KinCon biosensor allows the tracking of conformation dynamics of full-length kinases in intact cells and real time. We show that the most frequent BRAF cancer mutations affect kinase conformations and thus the engagement and efficacy of V600E-specific BRAF inhibitors (BRAFi). We illustrate that the patient mutation harboring KinCon reporters display differences in the effectiveness of the three clinically approved BRAFi vemurafenib, encorafenib, and dabrafenib and the preclinical paradox breaker PLX8394. We confirmed KinCon-based drug efficacy predictions for BRAF mutations other than V600E in proliferation assays using patient-derived lung cancer cell lines and by analyzing downstream kinase signaling. The systematic implementation of such conformation reporters will allow to accelerate the decision process for the mutation-oriented RAF-kinase cancer therapy. Moreover, we illustrate that the presented kinase reporter concept can be extended to other kinases which harbor patient mutations. Overall, KinCon profiling provides additional mechanistic insights into full-length kinase functions by reporting protein-protein interaction (PPI)-dependent, mutation-specific, and drug-driven changes of kinase activity conformations.

Dengue Virus Capsid Protein Facilitates Genome Compaction and Packaging.

Dengue virus (DENV) is a single-stranded (+)-sense RNA virus that infects humans and mosquitoes, posing a significant health risk in tropical and subtropical regions. Mature virions are composed of an icosahedral shell of envelope (E) and membrane (M) proteins circumscribing a lipid bilayer, which in turn contains a complex of the approximately 11 kb genomic RNA with capsid (C) proteins. Whereas the structure of the envelope is clearly defined, the structure of the packaged genome in complex with C proteins remains elusive. Here, we investigated the interactions of C proteins with viral RNA, in solution and inside mature virions, via footprinting and cross-linking experiments. We demonstrated that C protein interaction with DENV genomes saturates at an RNA:C protein ratio below 1:250. Moreover, we also showed that the length of the RNA genome interaction sites varies, in a multimodal distribution, consistent with the C protein binding to each RNA site mostly in singlets or pairs (and, in some instances, higher numbers). We showed that interaction sites are preferentially sites with low base pairing, as previously measured by 2'-acetylation analyzed by primer extension (SHAPE) reactivity indicating structuredness. We found a clear association pattern emerged: RNA-C protein binding sites are strongly associated with long-range RNA-RNA interaction sites, particularly inside virions. This, in turn, explains the need for C protein in viral genome packaging: the protein has a chief role in coordinating these key interactions, promoting proper packaging of viral RNA. Such sites are, thus, highly consequential for viral assembly, and, as such, may be targeted in future drug development strategies against these and related viruses.

Computational modelling of flavivirus dynamics: The ins and outs.

Enveloped viruses such as the flaviviruses represent a significant burden to human health around the world, with hundreds of millions of people each year affected by dengue alone. In an effort to improve our understanding of the molecular basis for the infective mechanisms of these viruses, extensive computational modelling approaches have been applied to elucidate their conformational dynamics. Multiscale protocols have been developed to simulate flavivirus envelopes in close accordance with biophysical data, in particular derived from cryo-electron microscopy, enabling high-resolution refinement of their structures and elucidation of the conformational changes associated with adaptation both to host environments and to immunological factors such as antibodies. Likewise, integrative modelling efforts combining data from biophysical experiments and from genome sequencing with chemical modification are providing unparalleled insights into the architecture of the previously unresolved nucleocapsid complex. Collectively, this work provides the basis for the future rational design of new antiviral therapeutics and vaccine development strategies targeting enveloped viruses.

Thrombin-derived C-terminal fragments aggregate and scavenge bacteria and their proinflammatory products.

Thrombin-derived C-terminal peptides (TCPs), including a major 11-kDa fragment (TCP96), are produced through cleavage by human neutrophil elastase and aggregate lipopolysaccharide (LPS) and the Gram-negative bacterium Escherichia coli However, the physiological roles of TCP96 in controlling bacterial infections and reducing LPS-induced inflammation are unclear. Here, using various biophysical methods, in silico molecular modeling, microbiological and cellular assays, and animal models, we examined the structural features and functional roles of recombinant TCP96 (rTCP96) in the aggregation of multiple bacteria and the Toll-like receptor (TLR) agonists they produce. We found that rTCP96 aggregates both Gram-negative and Gram-positive bacteria, including Staphylococcus aureus and Pseudomonas aeruginosa, and their cell-wall components LPS, lipid A, and lipoteichoic acid (LTA). The Gram-negative bacteria E. coli and P. aeruginosa were particularly sensitive to aggregation-induced bacterial permeabilization and killing. As a proof of concept, we show that rTCP96 reduces LPS-induced NF-κB activation in human monocytes, as well as in mouse models of LPS-induced subcutaneous inflammation. Moreover, in a mouse model of subcutaneous inoculation with P. aeruginosa, rTCP96 reduced bacterial levels. Together, these results link TCP-mediated aggregation of endotoxins and bacteria in vitro to attenuation of inflammation and bacterial levels in vivo.

3D reconstruction and flexibility of the hybrid engine Acetobacterium woodii F-ATP synthase.

The Na+-translocating F1FO ATP synthase from Acetobacterium woodii (AwF-ATP synthase) with a subunit stoichiometry of α3:ß3:γ:δ:ε:a:b2:(c2/3)9:c1 represents an evolutionary path between ATP-synthases and vacuolar ATPases, by containing a heteromeric rotor c-ring, composed of subunits c1, c2 and c3, and an extra loop (γ195-211) within the rotary γ subunit. Here, the recombinant AwF-ATP synthase was subjected to negative stain electron microscopy and single particle analysis. The reference free 2D class averages revealed high flexibility of the enzyme, wherein the F1 and FO domains distinctively bended to adopt multiple conformations. Moreover, both the F1 and FO domains tilted relative to each other to a maximum extent of 28° and 30°, respectively. The first 3D reconstruction of the AwF-ATP synthase was determined which accommodates well the modelled structure of the AwF-ATP synthase as well as the γ195-211-loop. Molecular simulations of the enzyme underlined the bending features and flexibility observed in the electron micrographs, and enabled assessment of the dynamics of the extra γ195-211-loop.

The Molecular Basis for Purine Binding Selectivity in the Bacterial ATP Synthase ϵ Subunit.

The ϵ subunit of ATP synthases has been proposed to regulate ATP hydrolysis in bacteria. Prevailing evidence supports the notion that when the ATP concentration falls below a certain threshold, the ϵ subunit changes its conformation from a non-inhibitory down-state to an extended up-state that then inhibits enzymatic ATP hydrolysis by binding to the catalytic domain. It has been demonstrated that the ϵ subunit from Bacillus PS3 is selective for ATP over other nucleotides, including GTP. In this study, the purine triphosphate selectivity is rationalized by using results from MD simulations and free energy calculations for the R103A/R115A mutant of the ϵ subunit from Bacillus PS3, which binds ATP more strongly than the wild-type protein. Our results are in good agreement with experimental data, and the elucidated molecular basis for selectivity could help to guide the design of novel GTP sensors.

Alpha-Carbonic Acid Revisited: Carbonic Acid Monomethyl Ester as a Solid and its Conformational Isomerism in the Gas Phase.

In this work, earlier studies reporting α-H2 CO3 are revised. The cryo-technique pioneered by Hage, Hallbrucker, and Mayer (HHM) is adapted to supposedly prepare carbonic acid from KHCO3 . In methanolic solution, methylation of the salt is found, which upon acidification transforms to the monomethyl ester of carbonic acid (CAME, HO-CO-OCH3 ). Infrared spectroscopy data both of the solid at 210 K and of the evaporated molecules trapped and isolated in argon matrix at 10 K are presented. The interpretation of the observed bands on the basis of carbonic acid [as suggested originally by HHM in their publications from 1993-1997 and taken over by Winkel et al., J. Am. Chem. Soc. 2007 and Bernard et al., Angew. Chem. Int. Ed. 2011] is inferior compared with the interpretation on the basis of CAME. The assignment relies on isotope substitution experiments, including deuteration of the OH- and CH3 - groups as well as 12 C and 13 C isotope exchange and on variation of the solvents in both preparation steps. The interpretation of the single molecule spectroscopy experiments is aided by a comprehensive calculation of high-level ab initio frequencies for gas-phase molecules and clusters in the harmonic approximation. This analysis provides evidence for the existence of not only single CAME molecules but also CAME dimers and water complexes in the argon matrix. Furthermore, different conformational CAME isomers are identified, where conformational isomerism is triggered in experiments through UV irradiation. In contrast to earlier studies, this analysis allows explanation of almost every single band of the complex spectra in the range between 4000 and 600 cm-1 .

How Ligand Binding Affects the Dynamical Transition Temperature in Proteins.

The biochemical functions of proteins are activated at the protein glass transition temperature, which has been proposed to be dependent upon protein-water interactions. However, at the molecular level it is unclear how ligand binding to well-defined binding sites can influence this transition temperature. We thus report molecular dynamics (MD) simulations of the ϵ subunit from thermophilic Bacillus PS3 in the ATP-free and ligand-bound states over a range of temperatures from 20 to 300 K, to study the influence of ligand association upon the transition temperature. We also measure the protein mean square displacement (MSD) in each state, which is well established as a means to quantify this dynamical temperature dependence. We find that the transition temperature is largely unaffected by ligand association, but the MSD beyond the transition temperature increases more rapidly in the ATP-free state. Our data suggests that ligands can effectively "shield" a binding site from solvent, and hence stabilize protein domains with increasing temperature.

Extending the Martini Coarse-Grained Force Field to N-Glycans.

Glycans play a vital role in a large number of cellular processes. Their complex and flexible nature hampers structure-function studies using experimental techniques. Molecular dynamics (MD) simulations can help in understanding dynamic aspects of glycans if the force field parameters used can reproduce key experimentally observed properties. Here, we present optimized coarse-grained (CG) Martini force field parameters for N-glycans, calibrated against experimentally derived binding affinities for lectins. The CG bonded parameters were obtained from atomistic (ATM) simulations for different glycan topologies including high mannose and complex glycans with various branching patterns. In the CG model, additional elastic networks are shown to improve maintenance of the overall conformational distribution. Solvation free energies and octanol-water partition coefficients were also calculated for various N-glycan disaccharide combinations. When using standard Martini nonbonded parameters, we observed that glycans spontaneously aggregated in the solution and required down-scaling of their interactions for reproduction of ATM model radial distribution functions. We also optimized the nonbonded interactions for glycans interacting with seven lectin candidates and show that a relatively modest scaling down of the glycan-protein interactions can reproduce free energies obtained from experimental studies. These parameters should be of use in studying the role of glycans in various glycoproteins and carbohydrate binding proteins as well as their complexes, while benefiting from the efficiency of CG sampling.

Aggregation of thrombin-derived C-terminal fragments as a previously undisclosed host defense mechanism.

Effective control of endotoxins and bacteria is crucial for normal wound healing. During injury, the key enzyme thrombin is formed, leading to generation of fibrin. Here, we show that human neutrophil elastase cleaves thrombin, generating 11-kDa thrombin-derived C-terminal peptides (TCPs), which bind to and form amorphous amyloid-like aggregates with both bacterial lipopolysaccharide (LPS) and gram-negative bacteria. In silico molecular modeling using atomic resolution and coarse-grained simulations corroborates our experimental observations, altogether indicating increased aggregation through LPS-mediated intermolecular contacts between clusters of TCP molecules. Upon bacterial aggregation, recombinantly produced TCPs induce permeabilization of Escherichia coli and phagocytic uptake. TCPs of about 11 kDa are present in acute wound fluids as well as in fibrin sloughs from patients with infected wounds. We noted aggregation and colocalization of LPS with TCPs in such fibrin material, which indicates the presence of TCP-LPS aggregates under physiological conditions. Apart from identifying a function of proteolyzed thrombin and its fragments, our findings provide an interesting link between the coagulation system, innate immunity, LPS scavenging, and protein aggregation/amyloid formation.

Structure and subunit arrangement of Mycobacterial F1FO ATP synthase and novel features of the unique mycobacterial subunit δ.

In contrast to other prokaryotes, the Mycobacterial F1FO ATP synthase (α3:ß3:γ:δ:ε:a:b:b':c9) is essential for growth. The mycobacterial enzyme is also unique as a result of its 111 amino acids extended δ subunit, whose gene is fused to the peripheral stalk subunit b. Recently, the crystallographic structures of the mycobacterial α3:ß3:γ:ε-domain and c subunit ring were resolved. Here, we report the first purification protocol of the intact M. smegmatis F1FO ATP synthase including the F1-domain, the entire membrane-embedded FO sector, and the stator subunits b' and the fused b-δ. This enzyme purification enabled the determination of the first projected 2D- and 3D structure of the intact M. smegmatis F1FO ATP synthase by electron microscopy (EM) and single particle analysis. Expression and purification of the fused mycobacterial b-δ24-446 construct, excluding the membrane-embedded N-terminal amino acids, provided insight into its secondary structure. By combining these data with homology and ab-initio modeling techniques, a model of the mycobacterial peripheral stalk subunits b-δ and b' was generated. Superposition of the 3D M. smegmatis F-ATP synthase EM-structure, the α3:ß3:γ:ε and c-ring, and the derived structural models of the peripheral stalk enabled a clear assignment of all F-ATP synthase subunits, in particular with respect to the unique mycobacterial peripheral stalk subunit b' and the elongated δ fused with subunit b. The arrangement of δ relative to the N-termini of the catalytic α3ß3-headpiece and its potential as a drug target are discussed.

Multiscale modeling of innate immune receptors: Endotoxin recognition and regulation by host defense peptides.

The innate immune system provides a first line of defense against foreign microorganisms, and is typified by the Toll-like receptor (TLR) family. TLR4 is of particular interest, since over-stimulation of its pathway by excess lipopolysaccharide (LPS) molecules from the outer membranes of Gram-negative bacteria can result in sepsis, which causes millions of deaths each year. In this review, we outline our use of molecular simulation approaches to gain a better understanding of the determinants of LPS recognition, towards the search for novel immunotherapeutics. We first describe how atomic-resolution simulations have enabled us to elucidate the regulatory conformational changes in TLR4 associated with different LPS analogues, and hence a means to rationalize experimental structure-activity data. Furthermore, multiscale modelling strategies have provided a detailed description of the thermodynamics and intermediate structures associated with the entire TLR4 relay - which consists of a number of transient receptor/coreceptor complexes - allowing us trace the pathway of LPS transfer from bacterial membranes to the terminal receptor complex at the plasma membrane surface. Finally, we describe our efforts to leverage these computational models, in order to elucidate previously undisclosed anti-inflammatory mechanisms of endogenous host-defense peptides found in wounds. Collectively, this work represents a promising avenue for the development of novel anti-septic treatments, inspired by nature's innate defense strategies.

Energetic Fingerprinting of Ligand Binding to Paralogous Proteins: The Case of the Apoptotic Pathway.

Networks of biological molecules are key to cellular function, governing processes ranging from signal cascade propagation to metabolic pathway regulation. Genetic duplication processes give rise to sets of regulatory proteins that have evolved from a common ancestor, leading to interactomes whose dysregulation is often associated with disease. A better understanding of the determinants of specificity at interfaces shared by functionally related proteins is crucial to the rational design of novel pharmacotherapeutic agents. To this end, a comprehensive data set of drug and drug-like binders was assembled for the Bcl-xL and Bcl-2 antiapoptotic proteins-archetypal examples of regulatory systems governed by evolutionarily conserved protein-protein interactions. These were first used to derive a two-dimensional quantitative structure-activity relationship (2D QSAR) model, predicting ligand specificity for these homologous proteins. The strengths and weaknesses of high-throughput 2D QSAR were then compared and contrasted to those of theoretically rigorous thermodynamic integration calculations performed on 14 complexes of Bcl-xL-specific, Bcl-2-specific, and potent dual binders bound to the Bcl-xL and Bcl-2 proteins. We demonstrate that free energy calculations provide an added layer of essential information, which traditional QSAR cannot capture. Moreover, we show that protein energetic responses to different ligands, expressed as per-residue energy values, can be used to fingerprint the protein-ligand interaction, extending the framework of four-dimensional molecular dynamics/quantitative structure-activity relationships (4D-MD/QSAR) toward the facilitation of future drug design strategies.

Hydrocarbons Are Essential for Optimal Cell Size, Division, and Growth of Cyanobacteria.

Cyanobacteria are intricately organized, incorporating an array of internal thylakoid membranes, the site of photosynthesis, into cells no larger than other bacteria. They also synthesize C15-C19 alkanes and alkenes, which results in substantial production of hydrocarbons in the environment. All sequenced cyanobacteria encode hydrocarbon biosynthesis pathways, suggesting an important, undefined physiological role for these compounds. Here, we demonstrate that hydrocarbon-deficient mutants of Synechococcus sp. PCC 7002 and Synechocystis sp. PCC 6803 exhibit significant phenotypic differences from wild type, including enlarged cell size, reduced growth, and increased division defects. Photosynthetic rates were similar between strains, although a minor reduction in energy transfer between the soluble light harvesting phycobilisome complex and membrane-bound photosystems was observed. Hydrocarbons were shown to accumulate in thylakoid and cytoplasmic membranes. Modeling of membranes suggests these compounds aggregate in the center of the lipid bilayer, potentially promoting membrane flexibility and facilitating curvature. In vivo measurements confirmed that Synechococcus sp. PCC 7002 mutants lacking hydrocarbons exhibit reduced thylakoid membrane curvature compared to wild type. We propose that hydrocarbons may have a role in inducing the flexibility in membranes required for optimal cell division, size, and growth, and efficient association of soluble and membrane bound proteins. The recent identification of C15-C17 alkanes and alkenes in microalgal species suggests hydrocarbons may serve a similar function in a broad range of photosynthetic organisms.

Systematic analysis of protein identity between Zika virus and other arthropod-borne viruses.

OBJECTIVE: To analyse the proportions of protein identity between Zika virus and dengue, Japanese encephalitis, yellow fever, West Nile and chikungunya viruses as well as polymorphism between different Zika virus strains. METHODS: We used published protein sequences for the Zika virus and obtained protein sequences for the other viruses from the National Center for Biotechnology Information (NCBI) protein database or the NCBI virus variation resource. We used BLASTP to find regions of identity between viruses. We quantified the identity between the Zika virus and each of the other viruses, as well as within-Zika virus polymorphism for all amino acid k-mers across the proteome, with k ranging from 6 to 100. We assessed accessibility of protein fragments by calculating the solvent accessible surface area for the envelope and nonstructural-1 (NS1) proteins. FINDINGS: In total, we identified 294 Zika virus protein fragments with both low proportion of identity with other viruses and low levels of polymorphisms among Zika virus strains. The list includes protein fragments from all Zika virus proteins, except NS3. NS4A has the highest number (190 k-mers) of protein fragments on the list. CONCLUSION: We provide a candidate list of protein fragments that could be used when developing a sensitive and specific serological test to detect previous Zika virus infections.