COVID-19: The Effect of Host Genetic Variations on Host-Virus Interactions.

Spurred into action by the COVID-19 pandemic, the global scientific community has, in a short of period of time, made astonishing progress in understanding and combating COVID-19. Given the known human protein machinery for (a) SARS-CoV-2 entry, (b) the host innate immune response, and (c) virus-host interactions (protein-protein and RNA-protein), the potential effects of human genetic variation in this machinery, which may contribute to clinical differences in SARS-CoV-2 pathogenesis and help determine individual risk for COVID-19 infection, are explored. The Genome Aggregation Database (gnomAD) was used to show that several rare germline exome variants of proteins in these pathways occur in the human population, suggesting that carriers of these rare variants (especially for proteins of innate immunity pathways) are at risk for severe symptoms (like the severe symptoms in patients who are known to be rare variant carriers), whereas carriers of other variants could have a protective advantage against infection. The occurrence of genetic variation is thus expected to motivate the experimental probing of natural variants to understand the mechanistic differences in SARS-CoV-2 pathogenesis from one individual to another.

A simple and sensitive SYBR Gold-based assay to quantify DNA-protein interactions.

A simple, accessible, and inexpensive assay to quantify the strength of DNA-protein interactions was developed. The assay relies on capturing DNA-protein complexes using an affinity resin that binds tagged, recombinant proteins. Sequential washes with filtration spin cups and centrifugation remove non-specific interactions in a gentle, uniform manner and a final elution isolates specific DNA-protein complexes. SYBR Gold nucleic acid stain is added to the eluted product and the fluorescence intensity accurately quantifies the amount of captured DNA, ultimately illustrating the relative strength of the DNA-protein interaction. The major utility of the assay resides in the versatility and quantitative nature of the SYBR Gold:nucleic acid interaction, eliminating the need for customized or labeled oligos and permitting relatively inexpensive quantification of binding capacity. The assay also employs DNA-protein complex capture by the very common purification tag, 6xHis, but other tags could likely be utilized. Further, SYBR Gold fluorescence is compatible with a wide variety of instruments, including UV transilluminators, a staple to any molecular biology laboratory. This assay was used to compare the binding capacities of different auxin response factor (ARF) transcription factors to various dsDNA targets, including the classical AuxRE motif and several divergent sequences. Results from dose-response assays suggest that different ARF proteins might show distinct comparative affinities for AuxRE variants, emphasizing that specific ARF-AuxRE binding strengths likely contribute to the complex and fine-tuned cellular auxin response.

A Comprehensive Analysis of Anion-Quadrupole Interactions in Protein Structures.

The edgewise interactions of anions with phenylalanine (Phe) aromatic rings in proteins, known as anion-quadrupole interactions, have been well studied. However, the anion-quadrupole interactions of the tyrosine (Tyr) and tryptophan (Trp) rings have been less well studied, probably because these have been considered weaker than interactions of anions hydrogen bonded to Trp/Tyr side chains. Distinguishing such hydrogen bonding interactions, we comprehensively surveyed the edgewise interactions of certain anions (aspartate, glutamate, and phosphate) with Trp, Tyr, and Phe rings in high-resolution, nonredundant protein single chains and interfaces (protein-protein, DNA/RNA-protein, and membrane-protein). Trp/Tyr anion-quadrupole interactions are common, with Trp showing the highest propensity and average interaction energy for this type of interaction. The energy of an anion-quadrupole interaction (-15.0 to 0.0 kcal/mol, based on quantum mechanical calculations) depends not only on the interaction geometry but also on the ring atom. The phosphate anions at DNA/RNA-protein interfaces interact with aromatic residues with energies comparable to that of aspartate/glutamate anion-quadrupole interactions. At DNA-protein interfaces, the frequency of aromatic ring participation in anion-quadrupole interactions is comparable to that of positive charge participation in salt bridges, suggesting an underappreciated role for anion-quadrupole interactions at DNA-protein (or membrane-protein) interfaces. Although less frequent than salt bridges in single-chain proteins, we observed highly conserved anion-quadrupole interactions in the structures of remote homologues, and evolutionary covariance-based residue contact score predictions suggest that conserved anion-quadrupole interacting pairs, like salt bridges, contribute to polypeptide folding, stability, and recognition.

Characteristics of a PHD Finger Subtype.

Although the plant homeodomain (PHD) finger superfamily is known for its site-specific readouts of histone tails, the origins of the mechanistic differences in histone H3 readout by different PHD subtypes remain less clear. We show that sequences containing the xCDxCDx motif in the PHD treble clef (xCDxCDx-PHD) constitute a distinct subtype, based on the following observations: (i) the amino acid composition of the binding site is strikingly different from other subtypes due to position-specific enrichment of negatively charged and bulky nonpolar residues, (ii) the binding site positions are mutually and positively correlated, and this correlation is absent in other subtypes, and (iii) there are only small structural deviations, despite low sequence similarity. The xCDxCDx-PHD constitutes ∼20% of the PHD family, and the double PHD fingers (DPFs) are 10% of the total number of xCDxCDx-PHDs. This subtype originated early in the evolution of eukaryotes but has diversified within the metazoan lineage. Despite sequence diversification, the positions of the enriched nonpolar residues, in particular, show very small structural deviations, suggesting critical contributions of nonpolar residues in the binding mechanism of this subtype. Using mutagenesis, we probed the contributions of the binding-site positions enriched in nonpolar residues in four xCDxCDx-PHD proteins and found that they contribute to the tight packing of the H3 residues. This effect may potentially be exploited, as we observed affinity enhancement upon substituting a bulky nonpolar residue at the same binding site in another histone reader. Overall, we present a detailed characterization of PHD subtypes.

Histone Peptide Recognition by KDM5B-PHD1: A Case Study.

A detailed understanding of the energetic contributions to histone peptide recognition would be valuable for a better understanding of chromatin anchoring mechanisms and histone diagnostic design. Here, we probed the energetic contributions to recognize the same unmodified histone H3 by three different plant homeodomain (PHD) H3K4me0 readers: hKDM5B-PHD1 (first PHD finger of hKDM5B), hBAZ2A-PHD, and hAIRE-PHD1. The energetic contributions of residues differ significantly from one complex to the next. For example, H3K4A substitution completely aborts the formation of the hAIRE-histone peptide complex, while it has only a small destabilizing effect on binding of the other readers, even though H3K4 methylation disrupts all three complexes. Packing density suggests that methylation of more tightly packed Lys/Arg residues can disrupt binding, even if the energetic contribution is small. The binding behavior of hKDM5B-PHD1 and hBAZ2A-PHD is similar, and like PHD H3R2 readers, both possess a pair of Asp residues in the treble clef for interaction with H3R2. PHD subtype sequences, especially the tandem PHD-PHD fingers, show enrichment in the treble clef Asp residues, suggesting that it is a subtype-specific property. These Asp residues make significant energetic contributions to the formation of the hKDM5B-histone peptide complex, suggesting that there are interactions in addition to those reported in the recent NMR structure. However, the presence of the treble clef Asp in PHD sequences may not always be sufficient for histone peptide binding. This study showcases reader-histone peptide interactions in the context of residue conservation, energetic contributions, interfacial packing, and sequence-based reader subtype predictability.

Isolation of a novel swine influenza virus from Oklahoma in 2011 which is distantly related to human influenza C viruses.

Of the Orthomyxoviridae family of viruses, only influenza A viruses are thought to exist as multiple subtypes and has non-human maintenance hosts. In April 2011, nasal swabs were collected for virus isolation from pigs exhibiting influenza-like illness. Subsequent electron microscopic, biochemical, and genetic studies identified an orthomyxovirus with seven RNA segments exhibiting approximately 50% overall amino acid identity to human influenza C virus. Based on its genetic organizational similarities to influenza C viruses this virus has been provisionally designated C/Oklahoma/1334/2011 (C/OK). Phylogenetic analysis of the predicted viral proteins found that the divergence between C/OK and human influenza C viruses was similar to that observed between influenza A and B viruses. No cross reactivity was observed between C/OK and human influenza C viruses using hemagglutination inhibition (HI) assays. Additionally, screening of pig and human serum samples found that 9.5% and 1.3%, respectively, of individuals had measurable HI antibody titers to C/OK virus. C/OK virus was able to infect both ferrets and pigs and transmit to naive animals by direct contact. Cell culture studies showed that C/OK virus displayed a broader cellular tropism than a human influenza C virus. The observed difference in cellular tropism was further supported by structural analysis showing that hemagglutinin esterase (HE) proteins between two viruses have conserved enzymatic but divergent receptor-binding sites. These results suggest that C/OK virus represents a new subtype of influenza C viruses that currently circulates in pigs that has not been recognized previously. The presence of multiple subtypes of co-circulating influenza C viruses raises the possibility of reassortment and antigenic shift as mechanisms of influenza C virus evolution.

Genomic and evolutionary characterization of a novel influenza-C-like virus from swine.

We recently described the isolation of a novel influenza virus from swine exhibiting respiratory disease in the United States that is distantly related to human influenza C virus. Based on genetic, biochemical and morphological analysis, the virus was provisionally classified as C/swine/Oklahoma/1334/2011 (C/OK). To further understand the genetics and evolution of this novel pathogen, we performed a comprehensive analysis of its sequence and phylogeny. The results demonstrated that C/OK and human influenza C viruses share a conserved array of predicted functional domains in the viral RNA genome replication and viral entry machinery but vary at key functional sites. Furthermore, our evolutionary analysis showed that homologous genes of C/OK and human influenza C viruses diverged from each other an estimated several hundred to several thousand years ago. Taken together, the findings described in this study support and extend our previous observations that C/OK is a genetically and evolutionarily distinct influenza virus in the family Orthomyxoviridae.

Assessing the functional roles of coevolving PHD finger residues.

Although in silico folding based on coevolving residue constraints in the deep-learning era has transformed protein structure prediction, the contributions of coevolving residues to protein folding, stability, and other functions in physical contexts remain to be clarified and experimentally validated. Herein, the PHD finger module, a well-known histone reader with distinct subtypes containing subtype-specific coevolving residues, was used as a model to experimentally assess the contributions of coevolving residues and to clarify their specific roles. The results of the assessment, including proteolysis and thermal unfolding of wildtype and mutant proteins, suggested that coevolving residues have varying contributions, despite their large in silico constraints. Residue positions with large constraints were found to contribute to stability in one subtype but not others. Computational sequence design and generative model-based energy estimates of individual structures were also implemented to complement the experimental assessment. Sequence design and energy estimates distinguish coevolving residues that contribute to folding from those that do not. The results of proteolytic analysis of mutations at positions contributing to folding were consistent with those suggested by sequence design and energy estimation. Thus, we report a comprehensive assessment of the contributions of coevolving residues, as well as a strategy based on a combination of approaches that should enable detailed understanding of the residue contributions in other large protein families.

AQcalc: A web server that identifies weak molecular interactions in protein structures.

Weak molecular interactions play an important role in protein structure and function. Computational tools that identify weak molecular interactions are, therefore, valuable for the study of proteins. Here, we present AQcalc, a web server (https://aqcalcbiocomputing.com/) that can be used to identify anion-quadrupole (AQ) interactions, which are weak interactions involving aromatic residue (Trp, Tyr, and Phe) ring edges and anions (Asp, Glu, and phosphate ion) both within proteins and at their interfaces (protein-protein, protein-nucleic acids, and protein-lipid bilayer). AQcalc identifies AQ interactions as well as clusters involving AQ, cation-π, and salt bridges, among others. Utilizing AQcalc we analyzed weak interactions in protein models, even in the absence of experimental structures, to understand the contributions of weak interactions to deleterious structural changes, including those associated with oncogenic and germline disease variants. We identified several deleterious variants with disrupted AQ interactions (comparable in frequency to cation-π disruptions). Amyloid fibrils utilize AQ to bury anions at frequencies that far exceed those observed for globular proteins. AQ interactions were detected three and five times more frequently than the hydrogen-bonded AQ (HBAQ) in fibril structures and protein-lipid bilayer interfaces, respectively. By contrast, AQ and HBAQ interactions were detected with similar frequencies in globular proteins. Collectively, these findings suggest AQcalc will be effective in facilitating fine structural analysis. As other web utilities designed to identify protein residue interaction networks do not report AQ interactions, wide use of AQcalc will enrich our understanding of residue interaction networks and facilitate hypothesis testing by identifying and experimentally characterizing these comparably weak but important interactions.

A single polymorphism in HIV-1 subtype C SP1 is sufficient to confer natural resistance to the maturation inhibitor bevirimat.

3-O-(3',3'-Dimethylsuccinyl) betulinic acid (DSB), also known as PA-457, bevirimat (BVM), or MPC-4326, is a novel HIV-1 maturation inhibitor. Unlike protease inhibitors, BVM blocks the cleavage of the Gag capsid precursor (CA-SP1) to mature capsid (CA) protein, resulting in the release of immature, noninfectious viral particles. Despite the novel mechanism of action and initial progress made in small-scale clinical trials, further development of bevirimat has encountered unexpected challenges, because patients whose viruses contain genetic polymorphisms in the Gag SP1 (positions 6 to 8) protein do not generally respond well to BVM treatment. To better define the role of amino acid residues in the HIV-1 Gag SP1 protein that are involved in natural polymorphisms to confer resistance to the HIV-1 maturation inhibitor BVM, a series of Gag SP1 chimeras involving BVM-sensitive (subtype B) and BVM-resistant (subtype C) viruses was generated and characterized for sensitivity to BVM. We show that SP1 residue 7 of the Gag protein is a primary determinant of SP1 polymorphism-associated drug resistance to BVM.

Structure and site-specific recognition of histone H3 by the PHD finger of human autoimmune regulator.

Human autoimmune regulator (AIRE) functions to control thymic expression of tissue-specific antigens via sequence-specific histone H3 recognition by its plant homeodomain (PHD) finger. Mutations in the AIRE PHD finger have been linked to autoimmune polyendocrinopathy-candidiasis-ectodermal dystrophy (APECED). Here we report the three-dimensional solution structure of the first PHD finger of human AIRE bound to a histone H3 peptide. The structure reveals a detailed network of interactions between the protein and the amino-terminal residues of histone H3, and particularly key electrostatic interactions of a conserved aspartic acid 297 in AIRE with the unmodified lysine 4 of histone H3 (H3K4). NMR binding study with H3 peptides carrying known posttranslational modifications flanking H3K4 confirms that transcriptional regulation by AIRE through its interactions with histone H3 is confined to the first N-terminal eight residues in H3. Our study offers a molecular explanation for the APECED mutations and helps define a subclass of the PHD finger family proteins that recognize histone H3 in a sequence-specific manner.

Systematic analysis of the effect of multiple templates on the accuracy of comparative models of protein structure.

BACKGROUND: Although multiple templates are frequently used in comparative modeling, the effect of inclusion of additional template(s) on model accuracy (when compared to that of corresponding single-template based models) is not clear. To address this, we systematically analyze two-template models, the simplest case of multiple-template modeling. For an existing target-template pair (single-template modeling), a two-template based model of the target sequence is constructed by including an additional template without changing the original alignment to measure the effect of the second template on model accuracy. RESULTS: Even though in a large number of cases a two-template model showed higher accuracy than the corresponding one-template model, over the entire dataset only a marginal improvement was observed on average, as there were many cases where no change or the reverse change was observed. The increase in accuracy due to the structural complementarity of the templates increases at higher alignment accuracies. The combination of templates showing the highest potential for improvement is that where both templates share similar and low (less than 30%) sequence identity with the target, as well as low sequence identity with each other. The structural similarity between the templates also helps in identifying template combinations having a higher chance of resulting in an improved model. CONCLUSION: Inclusion of additional template(s) does not necessarily improve model quality, but there are distinct combinations of the two templates, which can be selected a priori, that tend to show improvement in model quality over the single template model. The benefit derived from the structural complementarity is dependent on the accuracy of the modeling alignment. The study helps to explain the observation that a careful selection of templates together with an accurate target:template alignment are necessary to the benefit from using multiple templates in comparative modeling and provides guidelines to maximize the benefit from using multiple templates. This enables formulation of simple template selection rules to rank targets of a protein family in the context of structural genomics.

FRETting about the affinity of bimolecular protein-protein interactions.

Fluorescence resonance energy transfer (FRET) is a powerful tool to study macromolecular interactions such as protein-protein interactions (PPIs). Fluorescent protein (FP) fusions enable FRET-based PPI analysis of signaling pathways and molecular structure in living cells. Despite FRET's importance in PPI studies, FRET has seen limited use in quantifying the affinities of PPIs in living cells. Here, we have explored the relationship between FRET efficiency and PPI affinity over a wide range when expressed from a single plasmid system in Escherichia coli. Using live-cell microscopy and a set of 20 pairs of small interacting proteins, belonging to different structural folds and interaction affinities, we demonstrate that FRET efficiency can reliably measure the dissociation constant (KD ) over a range of mM to nM. A 10-fold increase in the interaction affinity results in 0.05 unit increase in FRET efficiency, providing sufficient resolution to quantify large affinity differences (> 10-fold) using live-cell FRET. This approach provides a rapid and simple strategy for assessment of PPI affinities over a wide range and will have utility for high-throughput analysis of protein interactions.

Accuracy of structure-derived properties in simple comparative models of protein structures.

The accuracy of comparative models of proteins is addressed here. A set of 12,732 single-template models of sequences of known high-resolution structures was built by an automated procedure. Accuracy of several structure-derived properties, such as surface area, residue accessibility, presence of pockets, electrostatic potential and others, was determined as a function of template:target sequence identity by comparing models with their corresponding experimental structures. As expected, the average accuracy of structure-derived properties always increases with higher template:target sequence identity, but the exact shape of this relationship can differ from one property to another. A comparison of structure-derived properties measured from NMR and X-ray structures of the same protein shows that for most properties, the NMR/X-ray difference is of the same order as the error in models based on approximately 40% template:target sequence identity. The exact sequence identity at which properties reach that accuracy varies between 25 and 50%, depending on the property being analyzed. A general characteristic of simple comparative models is that their surface has increased area as a consequence of being more rugged than that of experimental structures. This suggests that including solvent effects during model building or refinement could significantly improve the accuracy of surface properties in comparative models.

Systematic analysis of added-value in simple comparative models of protein structure.

Added-value is the additional information that a model carries with respect to the template structure used for model building. Thousands of single-template models, corresponding to proteins of known structure, were analyzed. The accuracy of structure-derived properties, such as residue accessibility, surface area, electrostatic potential, and others, was determined as a function of template:target sequence identity by comparing the models with their corresponding experimental structures. Added-value was determined by comparing the accuracy in models with that from templates. Geometry-dependent properties such as neighborhood of buried residues and accessible surface area showed low added-value. Properties that also depend on the protein sequence, such as presence of polar areas and electrostatic potential, showed high added-value. In general added-value increases when template:target sequence identity decreases, but it is also affected by alignment errors. This study justifies the use of models instead of the use of templates to estimate structure-derived properties of a target protein.

Phenethylisothiocyanate alters site- and promoter-specific histone tail modifications in cancer cells.

Site-specific histone modifications are important epigenetic regulators of gene expression. As deregulation of genes often results in complex disorders, corrective modulation of site-specific histone marks could be a powerful therapeutic or disease-preventive strategy. However, such modulation by dietary compounds and the resulting impact on disease risk remain relatively unexplored. Here we examined phenethylisothiocyanate (PEITC), a common dietary compound derived from cruciferous vegetables with known chemopreventive properties under experimental conditions, as a possible modulator of histone modifications in human colon cancer cells. The present study reports novel, dynamic, site-specific chemical changes to histone H3 in a gene-promoter-specific manner, associated with PEITC exposure in human colon tumor-derived SW480 epithelial cells. In addition, PEITC attenuated cell proliferation in a concentration- and time-dependent manner, likely mediated by caspase-dependent apoptotic signalling. The effects of PEITC on histone modifications and gene expression changes were achieved at low, non-cytotoxic concentrations, in contrast to the higher concentrations necessary to halt cancer cell proliferation. Increased understanding of specific epigenetic alterations by dietary compounds may provide improved chemopreventive strategies for reducing the healthcare burden of cancer and other human diseases.

Structural, antigenic, and evolutionary characterizations of the envelope protein of newly emerging Duck Tembusu Virus.

Since the first reported cases of ducks infected with a previously unknown flavivirus in eastern China in April 2010, the virus, provisionally designated Duck Tembusu Virus (DTMUV), has spread widely in domestic ducks in China and caused significant economic losses to poultry industry. In this study, we examined in detail structural, antigenic, and evolutionary properties of envelope (E) proteins of six DTMUV isolates spanning 2010-2012, each being isolated from individual farms with different geographical locations where disease outbreaks were documented. Structural analysis showed that E proteins of DTMUV and its closely related flavivirus (Japanese Encephalitis Virus) shared a conserved array of predicted functional domains and motifs. Among the six DTMUV strains, mutations were observed only at thirteen amino acid positions across three separate domains of the E protein. Interestingly, these genetic polymorphisms resulted in no detectable change in viral neutralization properties as demonstrated in a serum neutralization assay. Furthermore, phylogenetic analysis of the nucleotide sequences of the E proteins showed that viruses evolved into two distinct genotypes, termed as DTMUV.I and DTMUV.II, with II emerging as the dominant genotype. New findings described here shall give insights into the antigenicity and evolution of this new pathogen and provide guidance for further functional studies of the E protein for which no effective vaccine has yet been developed.

"η6"-Type anion-π in biomolecular recognition.

Theoretical and compelling experimental evidence indicates that the interaction between an anion and an aromatic π system when the anion is directly above the ring face ("η(6)"-type anion-π), can be attractive. This may play an important role in the formation and recognition of biomolecular structures. We examined high-resolution structures of proteins and nucleic acids for the presence of "η(6)"-type anion-π. Though less frequent than its counterpart cation-π, "η(6)"-type anion-π is observed unambiguously, occurring in protein/nucleic acid loops and often involving conserved/coevolving sites in proteins, suggesting it plays an important role in macromolecular folding and function.

Migration of the swine influenza virus δ-cluster hemagglutinin N-linked glycosylation site from N142 to N144 results in loss of antibody cross-reactivity.

Routine antigenic characterization of swine influenza virus isolates in a high-throughput serum neutralization (HTSN) assay found that approximately 20% of isolates were not neutralized by a panel of reference antisera. Genetic analysis revealed that nearly all of the neutralization-resistant isolates possessed a seasonal human-lineage hemagglutinin (HA; δ cluster). Subsequent sequencing analysis of full-length HA identified a conserved N144 residue present only in neutralization-resistant strains. N144 lies in a predicted N-linked glycosylation consensus sequence, i.e., N-X-S/T (where X is any amino acid except for proline). Interestingly, neutralization-sensitive viruses all had predicted N-linked glycosylation sites at N137 or N142, with threonine (T) occupying position 144 of HA. Consistent with the HTSN assay, hemagglutination inhibition (HI) and serum neutralization (SN) assays demonstrated that migration of the potential N-linked glycosylation site from N137 or N142 to N144 resulted in a >8-fold decrease in titers. These results were further confirmed in a reverse genetics system where syngeneic viruses varying only by predicted N-glycosylation sites at either N142 or N144 exhibited distinct antigenic characteristics like those observed in field isolates. Molecular modeling of the hemagglutinin protein containing N142 or N144 in complex with a neutralizing antibody suggested that N144-induced potential glycosylation may sterically hinder access of antibodies to the hemagglutinin head domain, allowing viruses to escape neutralization. Since N-linked glycosylation at these sites has been implicated in genetic and antigenic evolution of human influenza A viruses, we conclude that the relocation of the hemagglutinin N-linked glycosylation site from N142 to N144 renders swine influenza virus δ-cluster viruses resistant to antibody-mediated neutralization.

Systematic assessment of accuracy of comparative model of proteins belonging to different structural fold classes.

In the absence of experimental structures, comparative modeling continues to be the chosen method for retrieving structural information on target proteins. However, models lack the accuracy of experimental structures. Alignment error and structural divergence (between target and template) influence model accuracy the most. Here, we examine the potential additional impact of backbone geometry, as our previous studies have suggested that the structural class (all-α, αß, all-ß) of a protein may influence the accuracy of its model. In the twilight zone (sequence identity ≤ 30%) and at a similar level of target-template divergence, the accuracy of protein models does indeed follow the trend all-α > αß > all-ß. This is mainly because the alignment accuracy follows the same trend (all-α > αß > all-ß), with backbone geometry playing only a minor role. Differences in the diversity of sequences belonging to different structural classes leads to the observed accuracy differences, thus enabling the accuracy of alignments/models to be estimated a priori in a class-dependent manner. This study provides a systematic description of and quantifies the structural class-dependent effect in comparative modeling. The study also suggests that datasets for large-scale sequence/structure analyses should have equal representations of different structural classes to avoid class-dependent bias.