Model for the differential stabilities of rhinovirus and poliovirus to mild acidic pH, based on electrostatics calculations

The In Silico Prediction of Hotspot Residues that Contribute to the Structural Stability of Subunit Interfaces of a Picornavirus Capsid

The assembly of picornavirus capsids proceeds through the stepwise oligomerization of capsid protein subunits and depends on interactions between critical residues known as hotspots. Few studies have described the identification of hotspot residues at the protein subunit interfaces of the picornavirus capsid, some of which could represent novel drug targets. Using a combination of accessible web servers for hotspot prediction, we performed a comprehensive bioinformatic analysis of the hotspot residues at the intraprotomer, interprotomer and interpentamer interfaces of the Theiler’s murine encephalomyelitis virus (TMEV) capsid. Significantly, many of the predicted hotspot residues were found to be conserved in representative viruses from different genera, suggesting that the molecular determinants of capsid assembly are conserved across the family. The analysis presented here can be applied to any icosahedral structure and provides a platform for in vitro mutagenesis studies to further investigate the significance of these hotspots in critical stages of the virus life cycle with a view to identify potential targets for antiviral drug design.

The Critical Role of Codon Composition on the Translation Efficiency Robustness of the Hepatitis A Virus Capsid

Hepatoviruses show an intriguing deviated codon usage, suggesting an evolutionary signature. Abundant and rare codons in the cellular genome are scarce in the human hepatitis A virus (HAV) genome, while intermediately abundant host codons are abundant in the virus. Genotype–phenotype maps, or fitness landscapes, are a means of representing a genotype position in sequence space and uncovering how genotype relates to phenotype and fitness. Using genotype–phenotype maps of the translation efficiency, we have shown the critical role of the HAV capsid codon composition in regulating translation and determining its robustness. Adaptation to an environmental perturbation such as the artificial induction of cellular shutoff—not naturally occurring in HAV infection—involved movements in the sequence space and dramatic changes of the translation efficiency. Capsid rare codons, including abundant and rare codons of the cellular genome, slowed down the translation efficiency in conditions of no cellular shutoff. In contrast, rare capsid codons that are abundant in the cellular genome were efficiently translated in conditions of shutoff. Capsid regions very rich in slowly translated codons adapt to shutoff through sequence space movements from positions with highly robust translation to others with diminished translation robustness. These movements paralleled decreases of the capsid physical and biological robustness, and resulted in the diversification of capsid phenotypes. The deviated codon usage of extant hepatoviruses compared with that of their hosts may suggest the occurrence of a virus ancestor with an optimized codon usage with respect to an unknown ancient host.

Hepatitis A Virus Capsid Structure

Hepatitis A virus (HAV) has been enigmatic, evading detailed structural analysis for many years. Its recently determined high-resolution structure revealed an angular surface without the indentations often characteristic of receptor-binding sites. The viral protein 1 (VP1) β-barrel shows no sign of a pocket factor and the amino terminus of VP2 displays a "domain swap" across the pentamer interface, as in a subset of mammalian picornaviruses and insect picorna-like viruses. Structure-based phylogeny confirms this placement. These differences suggest an uncoating mechanism distinct from that of enteroviruses. An empty capsid structure reveals internal differences in VP0 and the VP1 amino terminus connected with particle maturation. An HAV/Fab complex structure, in which the antigen-binding fragment (Fab) appears to act as a receptor-mimic, clarifies some historical epitope mapping data, but some remain difficult to reconcile. We still have little idea of the structural features of enveloped HAV particles.

Molecular basis for the acid-initiated uncoating of human enterovirus D68

Enterovirus D68 (EV-D68) belongs to a group of enteroviruses that contain a single positive-sense RNA genome surrounded by an icosahedral capsid. Like common cold viruses, EV-D68 mainly causes respiratory infections and is acid-labile. The molecular mechanism by which the acid-sensitive EV-D68 virions uncoat and deliver their genome into a host cell is unknown. Using cryoelectron microscopy (cryo-EM), we have determined the structures of the full native virion and an uncoating intermediate [the A (altered) particle] of EV-D68 at 2.2- and 2.7-Å resolution, respectively. These structures showed that acid treatment of EV-D68 leads to particle expansion, externalization of the viral protein VP1 N termini from the capsid interior, and formation of pores around the icosahedral twofold axes through which the viral RNA can exit. Moreover, because of the low stability of EV-D68, cryo-EM analyses of a mixed population of particles at neutral pH and following acid treatment demonstrated the involvement of multiple structural intermediates during virus uncoating. Among these, a previously undescribed state, the expanded 1 ("E1") particle, shows a majority of internal regions (e.g., the VP1 N termini) to be ordered as in the full native virion. Thus, the E1 particle acts as an intermediate in the transition from full native virions to A particles. Together, the present work delineates the pathway of EV-D68 uncoating and provides the molecular basis for the acid lability of EV-D68 and of the related common cold viruses.

A generic approach to engineer antibody pH-switches using combinatorial histidine scanning libraries and yeast display

There is growing interest in the fast and robust engineering of protein pH-sensitivity that aims to reduce binding at acidic pH, compared to neutral pH. Here, we describe a novel strategy for the incorporation of pH-sensitive antigen binding functions into antibody variable domains using combinatorial histidine scanning libraries and yeast surface display. The strategy allows simultaneous screening for both, high affinity binding at pH 7.4 and pH-sensitivity, and excludes conventional negative selection steps. As proof of concept, we applied this strategy to incorporate pH-dependent antigen binding into the complementary-determining regions of adalimumab. After 3 consecutive rounds of separate heavy and light chain library screening, pH-sensitive variants could be isolated. Heavy and light chain mutations were combined, resulting in 3 full-length antibody variants that revealed sharp, reversible pH-dependent binding profiles. Dissociation rate constants at pH 6.0 increased 230- to 780-fold, while high affinity binding at pH 7.4 in the sub-nanomolar range was retained. Furthermore, binding to huFcRn and thermal stability were not affected by histidine substitutions. Overall, this study emphasizes a generalizable strategy for engineering pH-switch functions potentially applicable to a variety of antibodies and further proteins-based therapeutics.

Hepatitis A virus and the origins of picornaviruses

Hepatitis A virus (HAV) remains enigmatic, despite 1.4 million cases worldwide annually. It differs radically from other picornaviruses, existing in an enveloped form and being unusually stable, both genetically and physically, but has proved difficult to study. Here we report high-resolution X-ray structures for the mature virus and the empty particle. The structures of the two particles are indistinguishable, apart from some disorder on the inside of the empty particle. The full virus contains the small viral protein VP4, whereas the empty particle harbours only the uncleaved precursor, VP0. The smooth particle surface is devoid of depressions that might correspond to receptor-binding sites. Peptide scanning data extend the previously reported VP3 antigenic site, while structure-based predictions suggest further epitopes. HAV contains no pocket factor and can withstand remarkably high temperature and low pH, and empty particles are even more robust than full particles. The virus probably uncoats via a novel mechanism, being assembled differently to other picornaviruses. It utilizes a VP2 'domain swap' characteristic of insect picorna-like viruses, and structure-based phylogenetic analysis places HAV between typical picornaviruses and the insect viruses. The enigmatic properties of HAV may reflect its position as a link between 'modern' picornaviruses and the more 'primitive' precursor insect viruses; for instance, HAV retains the ability to move from cell-to-cell by transcytosis.

A histidine residue of the influenza virus hemagglutinin controls the pH dependence of the conformational change mediating membrane fusion

The conformational change of the influenza virus hemagglutinin (HA) protein mediating the fusion between the virus envelope and the endosomal membrane was hypothesized to be induced by protonation of specific histidine residues since their pKas match the pHs of late endosomes (pK(a) of ∼ 6.0). However, such critical key histidine residues remain to be identified. We investigated the highly conserved His184 at the HA1-HA1 interface and His110 at the HA1-HA2 interface of highly pathogenic H5N1 HA as potential pH sensors. By replacing both histidines with different amino acids and analyzing the effect of these mutations on conformational change and fusion, we found that His184, but not His110, plays an essential role in the pH dependence of the conformational change of HA. Computational modeling of the protonated His184 revealed that His184 is central in a conserved interaction network possibly regulating the pH dependence of conformational change via its pKa. As the propensity of histidine to get protonated largely depends on its local environment, mutation of residues in the vicinity of histidine may affect its pK(a). The HA of highly pathogenic H5N1 viruses carries a Glu-to-Arg mutation at position 216 close to His184. By mutation of residue 216 in the highly pathogenic as well as the low pathogenic H5 HA, we observed a significant influence on the pH dependence of conformational change and fusion. These results are in support of a pK(a)-modulating effect of neighboring residues.

IMPORTANCE

The main pathogenic determinant of influenza viruses, the hemagglutinin (HA) protein, triggers a key step of the infection process: the fusion of the virus envelope with the endosomal membrane releasing the viral genome. Whereas essential aspects of the fusion-inducing mechanism of HA at low pH are well understood, the molecular trigger of the pH-dependent conformational change inducing fusion has been unclear. We provide evidence that His184 regulates the pH dependence of the HA conformational change via its pK(a). Mutations of neighboring residues which may affect the pK(a) of His184 could play an important role in virus adaptation to a specific host. We suggest that mutation of neighboring residue 216, which is present in all highly pathogenic phenotypes of H5N1 influenza virus strains, contributed to the adaptation of these viruses to the human host via its effect on the pKa of His184.

Redox specificity of 2-hydroxyacid-coupled NAD(+)/NADH dehydrogenases: a study exploiting "reactive" arginine as a reporter of protein electrostatics

With "reactive" arginine as a kinetic reporter, 2-hydroxyacid dehydrogenases are assessed in basis of their specialization as NAD(+)-reducing or NADH-oxidizing enzymes. Specifically, M4 and H4 lactate dehydrogenases (LDHs) and cytoplasmic and mitochondrial malate dehydrogenases (MDHs) are compared to assess if their coenzyme specificity may involve electrostatics of cationic or neutral nicotinamide structure as the basis. The enzymes from diverse eukaryote and prokaryote sources thus are assessed in "reactivity" of functionally-critical arginine as a function of salt concentration and pH. Electrostatic calculations were performed on "reactive" arginines and found good correspondence with experiment. The reductive and oxidative LDHs and MDHs are assessed in their count over ionizable residues and in placement details of the residues in their structures as proteins. The variants found to be high or low in ΔpKa of "reactive" arginine are found to be also strong or weak cations that preferentially oxidize NADH (neutral nicotinamide structure) or reduce NAD(+) (cationic nicotinamide structure). The ionized groups of protein structure may thus be important to redox specificity of the enzyme on basis of electrostatic preference for the oxidized (cationic nicotinamide) or reduced (neutral nicotinamide) coenzyme. Detailed comparisons of isozymes establish that the residues contributing in their redox specificity are scrambled in structure of the reductive enzyme.

The Role of histidine residues in low-pH-mediated viral membrane fusion

A central event in the invasion of a host cell by an enveloped virus is the fusion of viral and cell membranes. For many viruses, membrane fusion is driven by specific viral surface proteins that undergo large-scale conformational rearrangements, triggered by exposure to low pH in the endosome upon internalization. Here, we present evidence suggesting that in both class I (helical hairpin proteins) and class II (beta-structure-rich proteins) pH-dependent fusion proteins the protonation of specific histidine residues triggers fusion via an analogous molecular mechanism. These histidines are located in the vicinity of positively charged residues in the prefusion conformation, and they subsequently form salt bridges with negatively charged residues in the postfusion conformation. The molecular surfaces involved in the corresponding structural rearrangements leading to fusion are highly conserved and thus might provide a suitable common target for the design of antivirals, which could be active against a diverse range of pathogenic viruses.

Crystal structure of Swine vesicular disease virus and implications for host adaptation

Swine vesicular disease virus (SVDV) is an Enterovirus of the family Picornaviridae that causes symptoms indistinguishable from those of foot-and-mouth disease virus. Phylogenetic studies suggest that it is a recently evolved genetic sublineage of the important human pathogen coxsackievirus B5 (CBV5), and in agreement with this, it has been shown to utilize the coxsackie and adenovirus receptor (CAR) for cell entry. The 3.0-A crystal structure of strain UK/27/72 SVDV (highly virulent) reveals the expected similarity in core structure to those of other picornaviruses, showing most similarity to the closest available structure to CBV5, that of coxsackievirus B3 (CBV3). Features that help to cement together and rigidify the protein subunits are extended in this virus, perhaps explaining its extreme tolerance of environmental factors. Using the large number of capsid sequences available for both SVDV and CBV5, we have mapped the amino acid substitutions that may have occurred during the supposed adaptation of SVDV to a new host onto the structure of SVDV and a model of the SVDV/CAR complex generated by reference to the cryo-electron microscopy-visualized complex of CBV3 and CAR. The changes fall into three clusters as follows: one lines the fivefold pore, a second maps to the CAR-binding site and partially overlaps the site for decay accelerating factor (DAF) to bind to echovirus 7 (ECHO7), and the third lies close to the fivefold axis, where the low-density lipoprotein receptor binds to the minor group of rhinoviruses. Later changes in SVDV (post-1971) map to the first two clusters and may, by optimizing recognition of a pig CAR and/or DAF homologue, have improved the adaptation of the virus to pigs.

Molecular modelling in structural biology

Molecular modelling is a powerful methodology for analysing the three dimensional structure of biological macromolecules. There are many ways in which molecular modelling methods have been used to address problems in structural biology. It is not widely appreciated that modelling methods are often an integral component of structure determination by NMR spectroscopy and X-ray crystallography. In this review we consider some of the numerous ways in which modelling can be used to interpret and rationalise experimental data and in constructing hypotheses that can be tested by experiment. Genome sequencing projects are producing a vast wealth of data describing the protein coding regions of the genome under study. However, only a minority of the protein sequences thus identified will have a clear sequence homology to a known protein. In such cases valuable three-dimensional models of the protein coding sequence can be constructed by homology modelling methods. Threading methods, which used specialised schemes to relate protein sequences to a library of known structures, have been shown to be able to identify the likely protein fold even in cases where there is no clear sequence homology. The number of protein sequences that cannot be assigned to a structural class by homology or threading methods, simply because they belong to a previously unidentified protein folding class, will decrease in the future as collaborative efforts in systematic structure determination begin to develop. For this reason, modelling methods are likely to become increasingly useful in the near future. The role of the blind prediction contests, such as the Critical Assessment of techniques for protein Structure Prediction (CASP), will be briefly discussed. Methods for modelling protein-ligand and protein-protein complexes are also described and examples of their applications given.

Role of an electrostatic network of residues in the enzymatic action of the Rhizomucor miehei lipase family

We have used continuum electrostatic methods to investigate the role of electrostatic interactions in the structure, function, and pH-dependent stability of the fungal Rhizomucor miehei lipase (RmL) family. We identify a functionally important electrostatic network which includes residues S144, D203, H257, Y260, H143, Y28, R80, and D91 (residue numbering is from RmL). This network consists of residues belonging to the catalytic triad (S144, D203, H257), residues located in proximity to the active site (Y260), residues stabilizing the geometry of the active site (Y28, H143), and residues located in the lid (D91) or close to the first hinge (R80). The lid and the first hinge are associated with the interfacial activation of lipases, where an alpha-helical lid opens up by rotating around two hinge regions. All network residues are well conserved in a set of 12 lipase homologues, and 6 of the network residues are located in sequence motifs. We observe that the effects of modeled mutations R86L, D91N, and H257F on the pH-dependent electrostatic free energies differ significantly in the closed and open conformations of RmL. Mutation R86L is especially interesting since it stabilizes the closed conformation but destabilizes the open one. Site-site electrostatic interaction energies reveal that interactions between R86 and D61, D113, and E117 stabilize the open conformation.

Charge and solvation effects in anion recognition centers: an inquiry exploiting reactive arginines

Following a long-standing suggestion of Riordan et al. [Riordan, J. F., McElvany, K. D., and Borders, C. L., Jr. (1977) Science 195, 884-885], we sought to exploit chemically activated arginines as probes to characterize the microenvironmental effects in enzymes that mediate the recognition of anionic substrates. A micellar simulation study establishes that octylguanidine (OGn) becomes chemically activated upon incorporation into both cetyltrimethylammonium bromide (CTAB) and Triton X-100 micelles and that the activations correlate with the pKa diminutions induced in its guanidinium group by the effects of electrostatic or nonelectrostatic nature as reflected in the results of pH and salt titration experiments. Next, a protein modification study establishes that the modifiable arginines in a number of enzymes also have diminished pKa's, again due to effects of electrostatic or nonelectrostatic nature as reflected in the results of pH and salt titration experiments. Warwicker's finite difference Poisson--Boltzmann algorithm [Warwicker, J. (1992) J. Mol. Biol. 223, 247-257] is applied to several of the enzymes with available crystal structure coordinates, and indeed, their chemically activated arginines are found to be in an electrostatic microenvironment that can diminish their pKa's, with the magnitudes of these diminutions matching closely the diminutions measured experimentally. Finally, the chemically activated arginines are examined with respect to their atomic atmosphere and are thus found to occur in a local microenvironment that would facilitate their roles as anion anchors. Thus, electrostatic and solvation effects are found to be critical determinants of the arginine role as an anion anchor.

Titration calculations of foot-and-mouth disease virus capsids and their stabilities as a function of pH

Foot-and-mouth disease virus (FMDV), a non-enveloped picornavirus, is sensitive to acidic conditions. At pH values below 7 the icosahedral virus capsid, formed from 60 copies of a protomer containing four polypeptides (VP1 to 4), dissociates into 12 pentamers, releasing the viral RNA. Evidence suggests that this acid lability may assist FMDV cell entry via an endosomal pathway. Calculations of titration curves and pH-stability profiles are presented for three different strains of FMDV, O1BFS, A10(61) and A22 Iraq, and compared with experimental data for complete virions and empty capsids (which lack RNA). The finite difference Poisson-Boltzmann method was used for the calculation of electrostatic free energies with the solvent treated as a dielectric continuum. The inter-pentamer interface in the virus is formed by two protomers related by 2-fold icosahedral symmetry. As a simple model for inter-pentamer interactions, a dimer and two separate protomers were compared. The association free energy was computed by integrating the difference between the titration curves of the two species. The calculations reproduced the observed decrease in capsid stability at acidic pH but not the difference in pH sensitivities of the two type A viruses. It is shown that only residues within 15 A of the interface play a significant role in determining acid lability. For the experimentally studied pH range (5 to 7.6), histidine residues were found to dominate the pH-dependence of the stability. Two histidine residues in VP3, H142 and H145, are shown to have the greatest effect by virtue of their interactions with many polar residues across the inter-pentamer interface; the interaction of H142 with an alpha-helix in the opposite pentamer contributes only a small proportion of the destabilization energy.

A review of protein engineering for the food industry

In this review I briefly describe the technique of protein engineering and indicate how the present state of knowledge allows proteins to be mutated to increase or decrease stability. I discuss experiments on both model proteins and those of relevance to the food industry and show how hydrophobic forces are a major driving force for folding as well as having a major role in thermostability. I also indicate the large contribution that hydrogen bonding, electrostatic interactions and, in a less well predicted way, disulfide bridges make to thermostability.

Characterization of an acid-resistant mutant of foot-and-mouth disease virus

A foot-and-mouth disease virus mutant which is stable at pH 6.4 has been isolated from a virus of serotype A. In contrast to the parent (P) virus, which gave a mixture of large and small plaques in BHK21 cells and in a bovine kidney cell line, the acid-resistant (AR) virus gave small plaques which did not increase markedly in size after 24 hr. The infectivity titer of the acid-resistant virus was about 100-fold lower in suckling mice than in BHK21 cells, whether the inoculation was made intraperitoneally or intracerebrally, whereas the parent virus gave similar titers in both systems. Furthermore, in mice the AR virus reached its end point two to three times more slowly. The diameter of the AR virus was almost 20% less than that of the P virus and it had a more distinct topography, but the two viruses cosedimented in sucrose gradients. However, the buoyant density in CsCl of the AR virus was slightly lower (1.42 compared with 1.43 g/cc) in coruns. The RNAs and capsid proteins of the two viruses gave similar profiles in sucrose gradients and by SDS-PAGE, respectively. However, isoelectric focusing of the capsid proteins revealed considerable differences between the two viruses. Whereas the P virus gave four protein bands, corresponding to VP1-VP4, the AR virus gave one band for VP4, two for VP3, two for VP2, and four for VP1. Sequence analysis of the genes coding for the capsid protein regions of the two viruses showed four changes (one silent), resulting in an Ala-3-->Ser substitution in VP1 and Glu-131-->Lys and Asp-133-->Ser substitutions in VP2.

Viral RNA modulates the acid sensitivity of foot-and-mouth disease virus capsids

Foot-and-mouth disease virus (FMDV) manifests an extreme sensitivity to acid, which is thought to be important for entry of the RNA genome into the cell. We have compared the low-pH-induced disassembly in vitro of virions and natural empty capsids of three subtypes of serotype A FMDV by enzyme-linked immunosorbent assay and sucrose gradient sedimentation analysis. For all three subtypes (A22 Iraq 24/64, A10(61), and A24 Cruzeiro), the empty capsid was more stable by 0.5 pH unit on average than the corresponding virion. Unexpectedly, in the natural empty capsids used in this study, the precursor capsid protein VP0 was found largely to be cleaved into VP2 and VP4. For picornaviruses the processing of VP0 is closely associated with encapsidation of viral RNA, which is considered likely to play a catalytic role in the cleavage. Investigation of the cleavage of VP0 in natural empty capsids failed to implicate the viral RNA. However, it remains possible that these particles arise from abortive attempts to encapsidate RNA. Empty capsids expressed from a vaccinia virus recombinant showed essentially the same acid lability as natural empty capsids, despite differing considerably in the extent of VP0 processing, with the synthetic particles containing almost exclusively uncleaved VP0. These results indicate that it is the viral RNA that modulates acid lability in FMDV. In all cases the capsids dissociate at low pH directly into pentameric subunits. Comparison of the three viruses indicates that FMDV A22 Iraq is about 0.5 pH unit more sensitive to low pH than types A10(61) and A24 Cruzeiro. Sequence analysis of the three subtypes identified several differences at the interface between pentamers and highlighted a His-alpha-helix dipole interaction which spans the pentamer interface and appears likely to influence the acid lability of the virus.