Raghavan Varadarajan

Design of an HA2-based Escherichia coli expressed influenza immunogen that protects mice from pathogenic challenge

Influenza HA is the primary target of neutralizing antibodies during infection, and its sequence undergoes genetic drift and shift in response to immune pressure. The receptor binding HA1 subunit of HA shows much higher sequence variability relative to the metastable, fusion-active HA2 subunit, presumably because neutralizing antibodies are primarily targeted against the former in natural infection. We have designed an HA2-based immunogen using a protein minimization approach that incorporates designed mutations to destabilize the low pH conformation of HA2. The resulting construct (HA6) was expressed in Escherichia coli and refolded from inclusion bodies. Biophysical studies and mutational analysis of the protein indicate that it is folded into the desired neutral pH conformation competent to bind the broadly neutralizing HA2 directed monoclonal 12D1, not the low pH conformation observed in previous studies. HA6 was highly immunogenic in mice and the mice were protected against lethal challenge by the homologous A/HK/68 mouse-adapted virus. An HA6-like construct from another H3 strain (A/Phil/2/82) also protected mice against A/HK/68 challenge. Regions included in HA6 are highly conserved within a subtype and are fairly well conserved within a clade. Targeting the highly conserved HA2 subunit with a bacterially produced immunogen is a vaccine strategy that may aid in pandemic preparedness.

Residue depth: A novel parameter for the analysis of protein structure and stability

BACKGROUND Accessible surface area is a parameter that is widely used in analyses of protein structure and stability. Accessible surface area does not, however, distinguish between atoms just below the protein surface and those in the core of the protein. In order to differentiate between such buried residues we describe a computational procedure for calculating the depth of a residue from the protein surface. RESULTS Residue depth correlates significantly better than accessibility with effects of mutations on protein stability and on protein-protein interactions. The deepest residues in the native state invariably undergo hydrogen exchange by global unfolding of the protein and are often significantly protected in the corresponding molten-globule states. CONCLUSIONS Depth is often a more useful gage of residue burial than accessibility. This is probably related to the fact that the protein interior and surrounding solvent differ significantly in polarity and packing density. Hence, the strengths of van der Waals and electrostatic interactions between residues in a protein might be expected to depend on the distance of the residue(s) from the protein surface.

Elucidation of determinants of protein stability through genome sequence analysis

Sequences of putative soluble proteins from complete genomes of eight thermophiles and 12 mesophiles were analyzed to gain insight into determinants of protein thermostability. The predator algorithm was used to assign secondary structures to each protein sequence. Based on simple statistical tests, a set of stabilizing factors was identified. These include reduced protein size, increases in number of residues involved in hydrogen bonding, β‐strand content and helix stabilization through ion pairs. There are also significant increases in the relative amounts of charged and hydrophobic β‐branched amino acids and decreases in uncharged polar amino acids in proteins from thermophiles relative to mesophilic organisms. Factors such as the relative proportion of residues in loops, proline and glycine content and helix capping do not appear to be important.

Influenza hemagglutinin stem-fragment immunogen elicits broadly neutralizing antibodies and confers heterologous protection.

Significance Hemagglutinin (HA), the major influenza virus envelope glycoprotein, is the principal target of neutralizing antibodies. Wide diversity and variation of HA entails annual vaccination, as current vaccines typically fail to elicit/boost cross-reactive, broadly neutralizing antibodies (bnAbs). Although several bnAbs bind at the conserved stem of HA making it an attractive universal vaccine candidate, the metastable conformation of this domain imposes challenges in designing a stable, independently folding HA stem immunogen. We rationally designed a stem-fragment immunogen, mimicking the native HA stem that binds conformation-specific bnAbs with high affinity. The immunogen elicited bnAbs and conferred robust protection against lethal, heterologous virus challenge in vivo. Additionally, soluble bacterial expression of such a thermotolerant, disulfide-free immunogen allows for rapid scale-up during pandemic outbreak. Influenza hemagglutinin (HA) is the primary target of the humoral response during infection/vaccination. Current influenza vaccines typically fail to elicit/boost broadly neutralizing antibodies (bnAbs), thereby limiting their efficacy. Although several bnAbs bind to the conserved stem domain of HA, focusing the immune response to this conserved stem in the presence of the immunodominant, variable head domain of HA is challenging. We report the design of a thermotolerant, disulfide-free, and trimeric HA stem-fragment immunogen which mimics the native, prefusion conformation of HA and binds conformation specific bnAbs with high affinity. The immunogen elicited bnAbs that neutralized highly divergent group 1 (H1 and H5 subtypes) and 2 (H3 subtype) influenza virus strains in vitro. Stem immunogens designed from unmatched, highly drifted influenza strains conferred robust protection against a lethal heterologous A/Puerto Rico/8/34 virus challenge in vivo. Soluble, bacterial expression of such designed immunogens allows for rapid scale-up during pandemic outbreaks.

Binding of nucleobases with single-walled carbon nanotubes: Theory and experiment

We report the binding energy of various nucleobases (guanine (G), adenine (A), thymine (T) and cytosine (C)) with (5,5) single-walled carbon nanotube (SWNT) calculated using first-principle Hartre–Fock method (HF) together with classical force field. The binding energy without including the solvation effects of water decreases in the order G>A>T>C. The inclusion of solvation energy changes the order of binding preference to be G>T>A>C. Using isothermal titration (micro) calorimetry experiments, we also show the relative binding affinity to be T>A>C, in agreement with our calculations.

Design of Escherichia coli-Expressed Stalk Domain Immunogens of H1N1 Hemagglutinin That Protect Mice from Lethal Challenge

ABSTRACT The hemagglutinin protein (HA) on the surface of influenza virus is essential for viral entry into the host cells. The HA1 subunit of HA is also the primary target for neutralizing antibodies. The HA2 subunit is less exposed on the virion surface and more conserved than HA1. We have previously designed an HA2-based immunogen derived from the sequence of the H3N2 A/HK/68 virus. In the present study, we report the design of an HA2-based immunogen from the H1N1 subtype (PR/8/34). This immunogen (H1HA0HA6) and its circular permutant (H1HA6) were well folded and provided complete protection against homologous viral challenge. Antisera of immunized mice showed cross-reactivity with HA proteins of different strains and subtypes. Although no neutralization was observable in a conventional neutralization assay, sera of immunized guinea pigs competed with a broadly neutralizing antibody, CR6261, for binding to recombinant Viet/04 HA protein, suggesting that CR6261-like antibodies were elicited by the immunogens. Stem domain immunogens from a seasonal H1N1 strain (A/NC/20/99) and a recent pandemic strain (A/Cal/07/09) provided cross-protection against A/PR/8/34 viral challenge. HA2-containing stem domain immunogens therefore have the potential to provide subtype-specific protection.

Depth: a web server to compute depth, cavity sizes, detect potential small-molecule ligand-binding cavities and predict the pKa of ionizable residues in proteins

Residue depth accurately measures burial and parameterizes local protein environment. Depth is the distance of any atom/residue to the closest bulk water. We consider the non-bulk waters to occupy cavities, whose volumes are determined using a Voronoi procedure. Our estimation of cavity sizes is statistically superior to estimates made by CASTp and VOIDOO, and on par with McVol over a data set of 40 cavities. Our calculated cavity volumes correlated best with the experimentally determined destabilization of 34 mutants from five proteins. Some of the cavities identified are capable of binding small molecule ligands. In this study, we have enhanced our depth-based predictions of binding sites by including evolutionary information. We have demonstrated that on a database (LigASite) of ∼200 proteins, we perform on par with ConCavity and better than MetaPocket 2.0. Our predictions, while less sensitive, are more specific and precise. Finally, we use depth (and other features) to predict pKas of GLU, ASP, LYS and HIS residues. Our results produce an average error of just <1 pH unit over 60 predictions. Our simple empirical method is statistically on par with two and superior to three other methods while inferior to only one. The DEPTH server (http://mspc.bii.a-star.edu.sg/depth/) is an ideal tool for rapid yet accurate structural analyses of protein structures.

Characterization of gp120 and its single-chain derivatives, gp120-CD4D12 and gp120-M9: implications for targeting the CD4i epitope in human immunodeficiency virus vaccine design.

ABSTRACT Single-chain derivatives of JRFL gp120 linked to the first two domains of human CD4 (gp120-CD4D12) or to the CD4 miniprotein analog CD4M9 (gp120-M9), have been constructed. Biacore studies revealed that gp120-CD4D12 and gp120-M9 bound to antibody 17b with dissociation constants of 0.8 and 25 nM, respectively, at pH 7.0, while gp120 alone did not bind. The binding of gp120-CD4D12 to 17b is not affected by the addition of excess soluble CD4D12, while the binding of gp120-M9 is enhanced. This finding indicates that the M9 component of the single chain interacts relatively weakly with gp120 and can be displaced by soluble CD4D12. Immunogenicity studies of gp120, gp120-CD4D12, and gp120-M9 were carried out with guinea pigs. All three molecules were highly immunogenic. The resulting antisera were examined for neutralizing activities against various human immunodeficiency virus type 1 isolates. Broadly neutralizing activity was observed only with sera generated against gp120-CD4D12. These antisera were depleted of anti-CD4D12 antibodies by being passed over a column containing immobilized CD4D12. The depleted sera showed a loss of broadly neutralizing activity. Sera that were affinity purified over a column containing immobilized gp120-M9 also lacked such neutralizing activity. This finding suggests that the broadly neutralizing response observed is exclusively due to anti-CD4 antibodies. Competition experiments showed that only antisera generated against gp120-CD4D12 competed with the CD4i antibody 17b and that this activity was not affected by depletion of anti-CD4 antibodies. The data indicate that although antibodies targeting the CD4i epitope were generated by the gp120-CD4D12 immunogen, these antibodies were nonneutralizing.

DEPTH: a web server to compute depth and predict small-molecule binding cavities in proteins

Depth measures the extent of atom/residue burial within a protein. It correlates with properties such as protein stability, hydrogen exchange rate, protein–protein interaction hot spots, post-translational modification sites and sequence variability. Our server, DEPTH, accurately computes depth and solvent-accessible surface area (SASA) values. We show that depth can be used to predict small molecule ligand binding cavities in proteins. Often, some of the residues lining a ligand binding cavity are both deep and solvent exposed. Using the depth-SASA pair values for a residue, its likelihood to form part of a small molecule binding cavity is estimated. The parameters of the method were calibrated over a training set of 900 high-resolution X-ray crystal structures of single-domain proteins bound to small molecules (molecular weight <1.5 KDa). The prediction accuracy of DEPTH is comparable to that of other geometry-based prediction methods including LIGSITE, SURFNET and Pocket-Finder (all with Matthew’s correlation coefficient of ∼0.4) over a testing set of 225 single and multi-chain protein structures. Users have the option of tuning several parameters to detect cavities of different sizes, for example, geometrically flat binding sites. The input to the server is a protein 3D structure in PDB format. The users have the option of tuning the values of four parameters associated with the computation of residue depth and the prediction of binding cavities. The computed depths, SASA and binding cavity predictions are displayed in 2D plots and mapped onto 3D representations of the protein structure using Jmol. Links are provided to download the outputs. Our server is useful for all structural analysis based on residue depth and SASA, such as guiding site-directed mutagenesis experiments and small molecule docking exercises, in the context of protein functional annotation and drug discovery.

Thermodynamic Effects of Proline Introduction on Protein Stability

The amino acid Pro is more rigid than other naturally occurring amino acids and, in proteins, lacks an amide hydrogen. To understand the structural and thermodynamic effects of Pro substitutions, it was introduced at 13 different positions in four different proteins, leucine–isoleucine–valine binding protein, maltose binding protein, ribose binding protein, and thioredoxin. Three of the maltose binding protein mutants were characterized by X‐ray crystallography to confirm that no structural changes had occurred upon mutation. In the remaining cases, fluorescence and CD spectroscopy were used to show the absence of structural change. Stabilities of wild type and mutant proteins were characterized by chemical denaturation at neutral pH and by differential scanning calorimetry as a function of pH. The mutants did not show enhanced stability with respect to chemical denaturation at room temperature. However, 6 of the 13 single mutants showed a small but significant increase in the free energy of thermal unfolding in the range of 0.3–2.4 kcal/mol, 2 mutants showed no change, and 5 were destabilized. In five of the six cases, the stabilization was because of reduced entropy of unfolding. However, the magnitude of the reduction in entropy of unfolding was typically several fold larger than the theoretical estimate of −4 cal K−1 mol−1 derived from the relative areas in the Ramachandran map accessible to Pro and Ala residues, respectively. Two double mutants were constructed. In both cases, the effects of the single mutations on the free energy of thermal unfolding were nonadditive. Proteins 2007.